Papers for Wednesday, Aug 18 2021

The existence of primordial black holes (PBHs), which may form from the collapse of matter overdensities shortly after the Big Bang, is still under debate. Among the potential signatures of PBHs are gravitational waves (GWs) emitted from binary black hole (BBH) mergers at redshifts $z\gtrsim 30$, where the formation of astrophysical black holes is unlikely. Future ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe equal-mass BBH mergers with total mass of $\mathcal{O}(10-100)~M_{\odot}$ at such distances. In this work, we investigate whether the redshift measurement of a single BBH source can be precise enough to establish its primordial origin. We simulate BBHs of different masses, mass ratios and orbital orientations. We show that for BBHs with total masses between $20~M_{\odot}$ and $40~M_{\odot}$ merging at $z \geq 40$ one can infer $z>30$ at up to 97\% credibility, with a network of one Einstein Telescope, one 40-km Cosmic Explorer in the US and one 20-km Cosmic Explorer in Australia. A smaller network made of one Einstein Telescope and one 40-km Cosmic Explorer in the US measures $z>30$ at larger than 90\% credibility for roughly half of the sources than the larger network. We then assess the dependence of this result on the Bayesian redshift priors used for the analysis, specifically on the relative abundance of the BBH mergers originated from the first stars, and the primordial BBH mergers.

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly-localized events have increased the overall number of detections, of which well-localized, "gold-plated" events make a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers.

We present observations of SN 2021csp, a unique supernova (SN) which displays evidence for interaction with H- and He- poor circumstellar material (CSM) at early times. Using high-cadence spectroscopy taken over the first week after explosion, we show that the spectra of SN 2021csp are dominated by C III lines with a velocity of 1800 km s$^{-1}$. We associate this emission with CSM lost by the progenitor prior to explosion. Subsequently, the SN displays narrow He lines before metamorphosing into a broad-lined Type Ic SN. We model the bolometric light curve of SN 2021csp, and show that it is consistent with the energetic ($4\times10^{51}$ erg) explosion of a stripped star, producing 0.4 M$_\odot$ of 56Ni within a $\sim$1 M$_\odot$ shell of CSM extending out to 400 R$_\odot$.

Radio interferometric experiments aim to constrain the reionization model parameters by measuring the 21-cm signal statistics, primarily the power spectrum. However the Epoch of Reionization (EoR) 21-cm signal is highly non-Gaussian, and this non-Gaussianity encodes important information about this era. The bispectrum is the lowest order statistic able to capture this inherent non-Gaussianity. Here we are the first to demonstrate that bispectra for large and intermediate length scales and for all unique $k$-triangle shapes provide tighter constraints on the EoR parameters compared to the power spectrum or the bispectra for a limited number of shapes of $k$-triangles. We use the Bayesian inference technique to constrain EoR parameters. We have also developed an Artificial Neural Network (ANN) based emulator for the EoR 21-cm power spectrum and bispectrum which we use to remarkably speed up our parameter inference pipeline. Here we have considered the sample variance and the system noise uncertainties corresponding to $1000$ hrs of SKA-Low observations for estimating errors in the signal statistics. We find that using all unique $k$-triangle bispectra improves the constraints on parameters by a factor of $2-4$ (depending on the stage of reionization) over the constraints that are obtained using power spectrum alone.

Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results.

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.

Relatively little information is available about the Universe at ultra-low radio frequencies, i.e. below 50 MHz (ULF), although the ULF spectral window contains a wealth of unique diagnostics for studying galactic and extragalactic phenomena. Sub-arcsecond resolution imaging at these frequencies is extremely difficult, due to the long baselines (>1000 km) required and large ionospheric perturbations. We have conducted a pilot project to investigate the ULF performance and potential of the International LOFAR Telescope (ILT), a trans-European interferometric array with baselines up to ~2000 km and observing frequencies down to 10 MHz. We have successfully produced images with sub-arcsecond resolution for 6 radio sources at frequencies down to 30 MHz. This is more than an order of magnitude better resolution than pre-ILT observations at similar frequencies. The six targets that we have imaged (3C 196, 3C 225, 3C 273, 3C 295, 3C 298 and 3C 380) are bright radio sources with compact structures. By comparing our data of 3C 196 and 3C 273 with observations at higher frequencies, we investigate their spatially resolved radio spectral properties. Our success shows that at frequencies down to 30 MHz, sub-arcsecond imaging with the ILT is possible. Further analysis is needed to determine the feasibility of observations of fainter sources or sources with less compact emission.

The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity.

Since its discovery in 1963, 3C273 has become one of the most widely studied quasars with investigations spanning the electromagnetic spectrum. While much has been discovered about this historically notable source, its low-frequency emission is far less well understood. Observations in the MHz regime have traditionally lacked the resolution required to explore small-scale structures that are key to understanding the processes that result in the observed emission. In this paper we use the first sub-arcsecond images of 3C273 at MHz frequencies to investigate the morphology of the compact jet structures and the processes that result in the observed spectrum. Using the full complement of LOFAR's international stations, we produce $0.31 \times 0.21$ arcsec images of 3C273 at 150 MHz to determine the jet's kinetic power, place constraints on the bulk speed and inclination angle of the jets, and look for evidence of the elusive counter-jet at 150 MHz. Using ancillary data at GHz frequencies, we fit free-free absorption (FFA) and synchrotron self-absorption (SSA) models to determine their validity in explaining the observed spectra. The images presented display for the first time that robust, high-fidelity imaging of low-declination complex sources is now possible with the LOFAR international baselines. We show that the main small-scale structures of 3C273 match those seen at higher frequencies and that absorption is present in the observed emission. We determine the kinetic power of the jet to be in the range of $3.5 \times 10^{43}$ - $1.5 \times 10^{44}$ erg s$^{-1}$ which agrees with estimates made using higher frequency observations. We derive lower limits for the bulk speed and Lorentz factor of $\beta \gtrsim 0.55$ and $\Gamma \geq 1.2$ respectively. The counter-jet remains undetected at $150$ MHz, placing a limit on the peak brightness of $S_\mathrm{cj\_150} < 40$ mJy beam$^{-1}$.

Active galactic nuclei (AGNs) show episodic activity, evident in galaxies that exhibit restarted radio jets. These restarted jets can interact with their environment, leaving signatures on the radio spectral energy distribution. Tracing these signatures requires resolved spectral index measurements over a broad frequency range including low frequencies. We present such a study for the radio galaxy 3C 293. Using the International LOFAR telescope (ILT) we probed spatial scales as fine as ~0.2" at 144 MHz, and to constrain the spectrum we combined these data with Multi-Element Radio Linked Interferometer Network (MERLIN) and Very Large Array (VLA) archival data. In the inner lobes (~2 kpc), we detect the presence of a spectral turnover that peaks at ~225 MHz and is most likely caused by free-free absorption from the rich surrounding medium. We confirm that these inner lobes are part of a jet-dominated young radio source (spectral age $\lesssim$0.17 Myr), which is strongly interacting with the rich interstellar medium (ISM) of the host galaxy. The outer lobes (~100 kpc) have a spectral index of $\alpha$~0.6-0.8 from 144-4850 MHz with a remarkably uniform spatial distribution and only mild spectral curvature ($\Delta\alpha\lesssim$ 0.2). We propose that intermittent fuelling and jet flow disruptions are powering the mechanisms that keep the spectral index in the outer lobes from steepening and maintain the spatial uniformity of the spectral index. Overall, it appears that 3C 293 has gone through multiple (two to three) epochs of activity. This study adds 3C 293 to the new sub-group of restarted galaxies with short interruption time periods. This is the first time a spatially resolved study simultaneously studies a young source as well as the older outer lobes at such low frequencies. This illustrates the potential of the ILT to expand such studies to a larger sample of radio galaxies.

Radio sources with steep spectra are preferentially associated with the most distant galaxies, the $\alpha-z$ relation, but the reason for this relation is an open question. The spatial distribution of spectra in high-z radio sources can be used to study this relation, and low-frequency observations are particularly important in understanding the particle acceleration and injection mechanisms. However, the small angular sizes of high-z sources together with the inherently low resolution of low-frequency radio telescopes until now has prevented high angular resolution low-frequency observations of distant objects. Here we present subarcsecond observations of a $z = 2.4$ radio galaxy at frequencies between $121$ MHz and $166$ MHz. We measure the spatial distribution of spectra, and discuss the implications for models of the $\alpha-z$ relation. We targeted 4C 43.15 with the High Band Antennas (HBAs) of the \textit{International LOFAR Telescope} (ILT) with a range of baselines up to $1300\ \mathrm{km}$. At the central frequency of $143$ MHz we achieve an angular resolution of $\sim 0.3''$. By complementing our data with archival \textit{Very Large Array} (VLA) data we study the spectral index distribution across 4C 43.15 between $55\ \mathrm{MHz}$ and $8.4\ \mathrm{GHz}$ at resolutions of $0.4''$ and $0.9''$. With a magnetic field strength of $B = 5.2$ nT and fitted injection indices of $\alpha^\mathrm{north}_\mathrm{inj} = -0.8$ and $\alpha^\mathrm{south}_\mathrm{inj} = -0.6$, fitting a Tribble spectral ageing model results in a spectral age of $\tau_\mathrm{spec} = 1.1 \pm 0.1$ Myr. We conclude that our data on 4C 43.15 indicates that inverse Compton losses could become comparable to or exceed synchrotron losses at higher redshifts and that inverse Compton losses could be a viable explanation for the $\alpha-z$ relation.

We study for the first time the low-frequency ($\sim$150 MHz) radio brightness distribution of Arp~299 at subarcsecond resolution, tracing in both compact and extended emission regions the local spectral energy distribution (SED) in order to characterize the dominant emission and absorption processes. We analysed the spatially resolved emission of Arp 299 revealed by 150 MHz international baseline Low-Frequency Array (LOFAR) and 1.4, 5.0, and 8.4 GHz Very Large Array (VLA) observations. We present the first subarcsecond (0.4"$\sim$100~pc) image of the whole Arp~299 system at 150~MHz. The high surface brightness sensitivity of our LOFAR observations ($\sim$100 $\mu$Jy/beam) allowed us to detect all of the nuclear components detected at higher frequencies, as well as the extended steep-spectrum emission surrounding the nuclei. We obtained spatially resolved, two-point spectral index maps for the whole galaxy: the compact nuclei show relatively flat spectra, while the extended, diffuse component shows a steep spectrum. We fitted the radio SED of the nuclear regions using two different models: a continuous free-free medium model and a clumpy model. The continuous model can explain the SED of the nuclei assuming a population of relativistic electrons subjected to synchrotron, bremsstrahlung, and ionization losses. The clumpy model fits assuming relativistic electrons with negligible energy losses, and thermal fractions that are more typical of star-forming galaxies than those required for the continuous model. Our results confirm the usefulness of combining spatially resolved radio imaging at both MHz and GHz frequencies to characterize in detail the radio emission properties of LIRGs from the central 100 pc out to the kiloparsec galaxy-wide scales.

We present Low-Frequency Array (LOFAR) telescope observations of the radio-loud gravitational lens systems MG 0751+2716 and CLASS B1600+434. These observations produce images at 300 milliarcseconds (mas) resolution at 150 MHz. In the case of MG 0751+2716, lens modelling is used to derive a size estimate of around 2 kpc for the low-frequency source, which is consistent with a previous 27.4 GHz study in the radio continuum with Karl G. Jansky Very Large Array (VLA). This consistency implies that the low-frequency radio source is cospatial with the core-jet structure that forms the radio structure at higher frequencies, and no significant lobe emission or further components associated with star formation are detected within the magnified region of the lens. CLASS B1600+434 is a two-image lens where one of the images passes through the edge-on spiral lensing galaxy, and the low radio frequency allows us to derive limits on propagation effects, namely scattering, in the lensing galaxy. The observed flux density ratio of the two lensed images is 1.19 +/- 0.04 at an observed frequency of 150 MHz. The widths of the two images give an upper limit of 0.035 kpc m^-20/3 on the integrated scattering column through the galaxy at a distance approximately 1 kpc above its plane, under the assumption that image A is not affected by scattering. This is relatively small compared to limits derived through very long baseline interferometry (VLBI) studies of differential scattering in lens systems. These observations demonstrate that LOFAR is an excellent instrument for studying gravitational lenses. We also report on the inability to calibrate three further lens observations: two from early observations that have less well determined station calibration, and a third observation impacted by phase transfer problems.

3C295 is a bright, compact steep spectrum source with a well-studied integrated radio spectral energy distribution (SED) from 132 MHz to 15 GHz. However, spatially resolved spectral studies have been limited due to a lack of high resolution images at low radio frequencies. These frequencies are crucial for measuring absorption processes, and anchoring the overall spectral modelling of the radio SED. In this paper, we use International LOFAR (LOw-Frequency ARray) Telescope (ILT) observations of 3C295 to study its spatially resolved spectral properties with sub-arcsecond resolution at 132 MHz. Combining our new 132 MHz observation with archival data at 1.6 GHz, 4.8 GHz, and 15 GHz, we are able to carry out a resolved radio spectral analysis. The spectral properties of the hotspots provides evidence for low frequency flattening. In contrast, the spectral shape across the lobes is consistent with a JP spectral ageing model. Using the integrated spectral information for each component, we then fit low-frequency absorption models to the hotspots, finding that both free-free absorption and synchrotron self-absorption models provide a better fit to the data than a standard power law. Although we can say there is low-frequency absorption present in the hotspots of 3C295, future observations with the Low Band Antenna of the ILT at 55 MHz may allow us to distinguish the type of absorption.

The functional renormalisation group is employed to study the non-linear regime of late-time cosmic structure formation. This framework naturally allows for non-perturbative approximation schemes, usually guided by underlying symmetries or a truncation of the theory space. An extended symmetry that is related to Galilean invariance is studied and corresponding Ward identities are derived. These are used to obtain (formally) closed renormalisation group flow equations for two-point correlation functions in the limit of large wave numbers (small scales). The flow equations are analytically solved in an approximation that is connected to the 'sweeping effect' previously described in the context of fluid turbulence.

Water is a molecule that is tightly related to many facets of star and planet formation. Water's abundance and distribution, especially the location of it's snowline has thus been the subject of much study. While water is seen to be abundant in the inner region of proto-planetary disks in infrared spectroscopy, detections of water in the disk in the sub-millimeter are rare, with only one detection towards AS 205. Here we put the multitude of non-detections and the single detection into context of recent physico-chemical models. We find that the 321.2257 GHz (10(2,9) - 9(3,6)) line detection towards AS 205 is inconsistent with a normal inner disk temperature structure and that the observed line must be masing. Furthermore, the emitting area derived from the line width, together with published analyses on water in disks around T-Tauri stars implies that the water snowline in the disk surface is at the same location as the snowline in the mid-plane. We propose that this is caused by vertical mixing continuously sequestering water from the warm surface layers into the cold disk midplane.

We present an investigation of the quiescent and transient X-ray binaries (XRBs) of the Galactic Center (GC). We extended our Chandra analysis of the non-thermal X-ray sources, located in the central parsec, from Hailey et al. (2018), using an additional 4.6 Msec of ACIS-S data obtained in 2012-2018. The individual Chandra spectra of the 12 sources fit to an absorbed power-law model with a mean photon index $\Gamma$~2 and show no Fe emission lines. Long-term variability was detected from nine of them, confirming that a majority are quiescent XRBs. Frequent X-ray monitoring of the GC revealed that the 12 non-thermal X-ray sources, as well as four X-ray transients have shown at most a single outburst over the last two decades. They are distinct from the six known neutron star LMXBs in the GC, which have all undergone multiple outbursts with <~ 5 year recurrence time on average. Based on the outburst history data of the broader population of X-ray transients, we conclude that the 16 sources represent a population of ~250-650 tightly-bound BH-LMXBs with ~4-12 hour orbital periods, consistent with the stellar/binary dynamics modelling in the vicinity of Sgr A*. The distribution of the 16 BH-LMXB candidates is disk-like (at 87% CL) and aligned with the nuclear star cluster. Our results have implications for XRB formation and the rate of gravitational wave events in other galactic nuclei.

The recent discovery of electromagnetic signals in coincidence with gravitational waves from neutron-star mergers has solidified the importance of multimessenger campaigns for studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as the LIGO/Virgo gravitational wave detectors and the IceCube neutrino observatory, record many candidate signals that fall short of the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by these candidates may, when combined with time-coincident gamma-ray observations, lead to significant detections. In this contribution, I describe our use of sub-threshold binary neutron star merger candidates identified in Advanced LIGO's first observing run (O1) to search for transient events in very-high-energy gamma rays using archival observations from the VERITAS imaging atmospheric Cherenkov telescope array. I describe the promise of this technique for future joint sub-threshold searches.

The development of new techniques for characterizing atmospheric optical turbulence (OT) has become an active topic of research again in recent years. In order to facilitate these studies, we reconsidered known theoretical results and obtained some new practically useful conclusions. We introduce a dimensionless Fresnel filter, which allows us to approximate a polychromatic weighting function (WF) by a monochromatic one with a typical precision of several percent. A so-called dimensionless WF can be easily scaled for a receiving aperture of any size. For the case of a circular aperture and monochromatic radiation, an analytical expression for the WF was found. The WFs for a square aperture and for a circular aperture match with relative difference less than 0.01 if the circular aperture diameter is 1.15 times larger than the square aperture side. A linear digital filter can be applied to the scintillation signal from an image detector. As an example of digital filtering, we considered the power law filter $\propto f^{5/3}$ with the WF being constant in a wide range of altitudes. We discuss the main limitations of this approach for measuring OT integral: finite pixel size, aliasing, and finite image detector size.

A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $\delta \rho_{max}$=30\%(50\%) with respect to underdense regions, and $\delta \rho_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers.

Solar system bodies with surface and sub-surface volatiles will show observational evidence of activity when they reach a temperature where those volatiles change from solid to gas and are released. This is most frequently seen in comets, where activity is driven by the sublimation of water, carbon dioxide, or carbon monoxide ices. However, some bodies (notably the asteroid (3200) Phaethon) show initiation of activity at very small heliocentric distances, long after they have reached the sublimation temperatures of these ices. We investigate whether the sodium present in the mineral matrix could act as the volatile element responsible for this activity. We conduct theoretical modeling which indicates that sodium has the potential to sublimate in the conditions that Phaethon experiences, depending on the mineral phase it is held in. To test this, we then exposed samples of the carbonaceous chondrite Allende to varying heating events similar to what would be experienced by low perihelion asteroids. We measured the change in sodium present in each sample, and find that the highest temperature samples show a significant loss of sodium from specific mineral phases over a single heating event, comparable to a day on the surface of Phaethon. Under specific thermal histories possible for Phaethon, this outgassing could be sufficient to explain this object's observed activity. This effect would also be expected to be observed for other low-perihelia asteroids as well, and may act as a critical step in the process of disrupting small low-albedo asteroids.

In an EHT study of a Jy-level target, Safarzadeh et al. (2019) show how astrometric monitoring could constrain massive black hole binaries with the wide separations that make them long-lived against gravitational wave losses, and with the small mass ratios expected from merged satellite galaxies. With this ngVLA study, we show how such frontier topics could be explored for the more numerous mJy-level targets, such as NGC\,4472. We also discuss how ngVLA astrometric monitoring could test the upper limits from pulsar timing arrays on gravitational waves from NGC\,4472.

Cool clouds are expected to be destroyed and incorporated into hot supernova-driven galactic winds. The mass-loading of a wind by the cool medium modifies the bulk velocity, temperature, density, entropy, and abundance profiles of the hot phase relative to an un-mass-loaded outflow. We provide general equations and limits for this physics that can be used to infer the rate of cool gas entrainment from X-ray observations, accounting for non-spherical expansion. In general, mass-loading flattens the density and temperature profiles, decreases the velocity and increases the entropy if the Mach number is above a critical value. We first apply this model to a recent high-resolution galactic outflow simulation where the mass-loading can be directly inferred. We show that the temperature, entropy, and composition profiles are well-matched, providing evidence that this physics sets the bulk hot gas profiles. We then model the diffuse X-ray emission from the local starburst M82. The non-spherical (more cylindrical) outflow geometry is directly taken from the observed X-ray surface brightness profile. These models imply a total mass-loading rate that is about equal to that injected in the starburst, $\simeq 10$ M$_\odot$ yr$^{-1}$, and they predict an asymptotic hot wind velocity of $\sim 1000\,{\rm km \ s^{-1}}$ that is $\sim1.5-2$ times smaller than previous predictions. We also show how the observed entropy profile can be used to constrain the outflow velocity, making predictions for future missions like XRISM. We argue that the observed X-ray limb-brightening may be explained by mass-loading at the outflow's edges.

We present a suite of 16 high-resolution hydrodynamic simulations of an isolated dwarf galaxy (gaseous and stellar disk plus a stellar bulge) within an initially cuspy dark matter (DM) halo, including self-interactions between the DM particles (SIDM); as well as stochastic star formation and subsequent supernova feedback (SNF), implemented using the stellar feedback model SMUGGLE. SIDM momentum transfer cross section and star formation threshold are varied between simulations. The DM halo forms a constant density core of similar size and shape for several combinations of those two parameters. Haloes with cores which are formed due to SIDM (adiabatic cusp-core transformation) have velocity dispersion profiles which are closer to isothermal than those of haloes with cores which are formed due to SNF in simulations with bursty star formation (impulsive cusp-core transformation). Impulsive SNF can generate steep stellar age gradients and increase random motion in the gas at the centre of the galaxy. Simulated galaxies in haloes with cores which were formed adiabatically are spatially more extended, with stellar metallicity gradients that are shallower (at late times) than those of galaxies in other simulations. Such observable properties of the gas and the stars, which indicate either an adiabatic or an impulsive evolution of the gravitational potential, may be used to determine whether observed cores in DM haloes are formed through self-interactions between the DM particles or in response to impulsive SNF.

In this work we analyze the orbital evolution and the dynamical stability in the vicinity of the small Saturnian moons Aegaeon, Methone, Anthe and Pallene. We numerically resolve the exact equations of motions to investigate the orbital motion of thousands of test particles within and near to the domain of the 7/6, 14/15, 10/11 mean motion resonances of Aegaeon, Methone and Anthe with Mimas, respectively. We show that, for massless small moons, the orbits of particles initially restricted to the resonant domains remain stable for at least $10^4$ yr. We also conduct numerical simulations considering Aegaeon, Methone, Anthe and Pallene as massive bodies. The results show that most particles undergo significant perturbations in their orbital motions, ultimately destabilizing in timescales of a few hundreds of years or even less through collisions with the four small moons. In addition, we also simulate the orbital evolution of test particles initially distributed in form of arc around Aegaeon, Methone and Anthe. We show that the initial arcs are dynamically eroded on timescales of hundreds of years, allowing us to constraint the timescales for which gravitational forces operate to remove particles from the observed arcs.

Models of solar-like oscillators yield acoustic modes at different frequencies than would be seen in actual stars possessing identical interior structure, due to modelling error near the surface. This asteroseismic "surface term" must be corrected when mode frequencies are used to infer stellar structure. Subgiants exhibit oscillations of mixed acoustic ($p$-mode) and gravity ($g$-mode) character, which defy description by the traditional $p$-mode asymptotic relation. Since nonparametric diagnostics of the surface term rely on this description, they cannot be applied to subgiants directly. In Paper I, we generalised such nonparametric methods to mixed modes, and showed that traditional surface-term corrections only account for mixed-mode coupling to, at best, first order in a perturbative expansion. Here, we apply those results, modelling subgiants using asteroseismic data. We demonstrate that, for grid-based inference of subgiant properties using individual mode frequencies, neglecting higher-order effects of mode coupling in the surface term results in significant systematic differences in the inferred stellar masses, and measurable systematics in other fundamental properties. While these systematics are smaller than those resulting from other choices of model construction, they persist for both parametric and nonparametric formulations of the surface term. This suggests that mode coupling should be fully accounted for when correcting for the surface term in seismic modelling with mixed modes, irrespective of the choice of correction used. The inferred properties of subgiants, in particular masses and ages, also depend on the choice of surface-term correction, in a different manner from both main-sequence and red giant stars.

The vector vortex coronagraph (VVC) performance in the laboratory and in ground-based observatories has earned it a spot on the NASA mission concepts HabEx and LUVOIR. The VVC induces a phase ramp through the manipulation of the polarization state. Left- and right-circular polarizations get imprinted a phase ramp of opposite signs, which prevents model-based focal plane wavefront sensing and control strategies in natural light. We thus have to work with a polarization state than ensures circularly polarized light at the VVC mask. However, achieving this polarization state can be non trivial if there are optics that add phase retardance of any kind between the circular polarizer and the focal plane mask. Here we present the method currently used at the Caltech high contrast spectroscopy testbed (HCST) to achieve the proper circular polarization state for a VVC, which only uses the deformable mirror and appropriate rotation of the circular polarizer and analyzer optics. At HCST we achieve raw contrast levels of \tentoe~for broadband light with a VVC.

It is shown that, from the two independent approaches of McCrea-Milne and of Zeldovich, one can fully recover the set equations corresponding to relativistic equations of the expanding universe of Friedmann-Lemaitre-Robertson-Walker geometry. Although similar, the Newtonian and relativistic set of equations, have principal difference in the content and hence define two flows, local and global ones, thus naturally exposing the Hubble tension at presence of the cosmological constant \Lambda. From that, we obtain "absolute" constraints on the lower and upper values for the local Hubble parameter, \sqrt{\Lambda c^2/3} \simeq 56.2 and \sqrt{\Lambda c^2} \simeq 97.3$ (km/sec Mpc^{-1}), respectively. The link to the so-called "maximum force/tension" issue in cosmological models is revealed.

Quasi-periodic oscillations (QPOs) are an important key to understand the dynamic behavior of astrophysical objects during transient events like gamma-ray bursts, solar flares, and magnetar flares. Searches for QPOs often use the periodogram of the time series and perform spectral density estimation using a Whittle likelihood function. However, the Whittle likelihood is only valid if the time series is stationary since the frequency bins are otherwise not statistically independent. We show that if time series are non-stationary, the significance of QPOs can be highly overestimated and estimates of the central frequencies and QPO widths can be overconstrained. The effect occurs if the QPO is only present for a fraction of the time series and the noise level is varying throughout the time series. This can occur for example if background noise from before or after the transient is included in the time series or if the low-frequency noise profile varies strongly over the time series. We confirm the presence of this bias in previously reported results from solar flare data and show that significance can be highly overstated. Finally, we provide some suggestions that help identify if an analysis is affected by this bias.

We present a census of neutral gas in the Milky Way disk and halo down to limiting column densities of $N$(HI)$\sim10^{14}$ cm$^{-2}$ using measurements of HI Lyman-series absorption from the Far Ultraviolet Spectroscopic Explorer (FUSE). Our results are drawn from an analysis of 25 AGN sightlines spread evenly across the sky with Galactic latitude |b|$\gtrsim 20^{\circ}$. By simultaneously fitting multi-component Voigt profiles to 11 Lyman-series absorption transitions covered by FUSE (Ly$\beta$-Ly$\mu$) plus HST measurements of Ly$\alpha$, we derive the kinematics and column densities of a sample of 152 HI absorption components. While saturation prevents accurate measurements of many components with column densities 17$\lesssim$log$N$(HI)$\lesssim$19, we derive robust measurements at log$N$(HI)$\lesssim$17 and log$N$(HI)$\gtrsim$19. We derive the first ultraviolet HI column density distribution function (CDDF) of the Milky Way, both globally and for low-velocity (ISM), intermediate-velocity clouds (IVCs), and high-velocity clouds (HVCs). We find that IVCs and HVCs show statistically indistinguishable CDDF slopes, with $\beta_{\rm IVC}=$ $-1.01_{-0.14}^{+0.15}$ and $\beta_{\rm HVC}=$ $-1.05_{-0.06}^{+0.07}$. Overall, the CDDF of the Galactic disk and halo appears shallower than that found by comparable extragalactic surveys, suggesting a relative abundance of high-column density gas in the Galactic halo. We derive the sky covering fractions as a function of HI column density, finding an enhancement of IVC gas in the northern hemisphere compared to the south. We also find evidence for an excess of inflowing HI over outflowing HI, with $-$0.88$\pm$0.40 M$_\odot$ yr$^{-1}$ of HVC inflow versus 0.20$\pm$0.10 M$_\odot$ yr$^{-1}$ of HVC outflow, confirming an excess of inflowing HVCs seen in UV metal lines.

The lack of adequate X-ray observing capability is seriously impeding the progress in understanding the hot phase of circumgalactic medium (CGM), which is predicted to extend to the virial radius of a galaxy or beyond, and thus in acquiring key boundary conditions for studying galaxy evolution. To this end, the Hot Universe Baryon Surveyor (\textit{HUBS}) is proposed. \textit{HUBS} is designed to probe hot CGM by detecting its emission or absorption lines with a non-dispersive X-ray spectrometer of high resolution and high throughput. The spectrometer consists of a $60\times60$ array of microcalorimeters, with each detector providing an energy resolution of $2~\mathrm{eV}$, and is placed in the focal plane of an X-ray telescope of $1^{\circ}$ field-of-view. With such a design, the spectrometer is also expected to enable studies of intra-group medium (IGrM) and the outer region of intra-cluster medium (ICM). To assess the scientific potential of \textit{HUBS}, we created mock observations of galaxies, groups, and clusters at different redshifts with the \tng simulation. Focusing exclusively on emission studies in this work, we took into account the effects of light cone, Galactic foreground emission, and background AGN contribution in the mock observations. From the observations, we made mock X-ray images and spectra, analyzed them to derive the properties of the emitting gas in each case, and compared the results with the input parameters from the simulation. The results show that \textit{HUBS} is well suited for studying hot CGM at low redshifts. The redshift range is significantly extended for measuring IGrM and ICM.

Some surviving companions of type Ia supernovae (SNe Ia) from the white-dwarf + main-sequence (WD + MS) channel may evolve to hot subdwarfs. In this paper, we preformed stellar evolution calculations for surviving companions of close WD + MS systems in the spin-up/spin-down model and the canonical non-rotating model to map out the initial parameter spaces in the orbital period - secondary mass plane in which the surviving companions can evolve to hot subdwarfs. Based on these results, we carried out a series of binary population synthesis calculation to obtain the Galactic birth rate of the hot subdwarfs from the WD + MS channel, which is $2.3-6\times10^{\rm -4}\,{\rm yr}^{\rm -1}$ for the spin-up/spin-down model and $0.7-3\times10^{\rm -4}\,{\rm yr}^{\rm -1}$ for the canonical non-rotating model. We also show the distributions of some integral properties of the hot subdwarfs, e.g. the mass and the space velocity, for different models. In addition, comparing our results with the observations of the intermediate helium-rich (iHe-rich) hot subdwarfs, the hot subdwarfs from the WD + MS channel may explain some observational features of the iHe-rich hot subdwarfs, especially for those from the spin-up/spin-down model. Although we expect that the SN Ia channel can only contribute a small fraction of the iHe-rich hot subdwarf population, some of these may help to explain cases with unusual kinematics.

In this paper, we extend previous works on the relation between mass and the inner slope in dark matter density profiles. We calculate that relation in the mass range going from dwarf galaxies to cluster of galaxies. This was done thanks to a modeling of energy transfer via SN and AGN feedback, as well as via dynamical friction of baryon clumps. We show that, in the mass range above galaxy masses (Groups and clusters), the inner slope-mass relation changes its trend. It flattens (towards less cuspy profile) around masses corresponding to groups of galaxies and steepens again for large galaxy cluster masses. The flattening is produced by the AGN outflows (AGN feedback). The one-$ \sigma$ scatter on $\alpha$ is approximately constant in all the mass range ($\Delta\alpha\simeq 0.3$). This is the first paper extending the inner density profile slope-mass relationship to clusters of galaxies, accounting for the role of baryons. The result can be used to obtain a complete density profile, also taking baryons into account. Such kind of density profile was previously only available for galaxies.

Galactic black hole candidate MAXI J1910-057/Swift J1910.2-0546 was simultaneously discovered by MAXI/GSC and Swift/BAT satellites during its first outburst in 2012. We study detailed spectral and temporal properties of the source in a broad energy range using archival data from Swift/XRT, MAXI/GSC and Swift/BAT instruments. Low frequency quasi periodic oscillations are observed on several instances of observations. The combined 1-50 keV spectra are analyzed using the transonic flow solution based Two Component Advective Flow (TCAF) model. Based on the variation of soft and hard X-ray fluxes, their hardness ratios and the variation of the spectral model parameters, we found that the source has evolved through six spectral states, which is unusual in the case of transient black holes. We interpret this unusual spectral state evolution to be a result of two overlapping outbursts where the leftover matter from a primary outburst is released from the pile-up radius due to a sudden rise of viscosity causing a secondary outburst. From the spectral analysis with TCAF model, we estimate the mass of the black hole to be $9.98^{+3.54}_{-3.09}$ $M_\odot$ , and the source distance is estimated to be $4.6-14.6$~kpc from transition luminosity considerations.

We present far-infrared (FIR) spectroscopy of supernova remnants (SNRs) based on the archival data of the Infrared Space Observatory ($ISO$) taken with the Long Wavelength Spectrometer (LWS). Our sample includes previously unpublished profiles of line and continuum spectra for 20 SNRs in the Galaxy and Magellanic Clouds. In several SNRs including G21.5-0.9, G29.7-0.3, the Crab Nebula, and G320.4-1.2, we find evidence for broad [O I], [O III], [N II], and [C II] lines with velocity dispersions up to a few 10$^3$ km s$^{-1}$, indicating that they are associated with high-velocity SN ejecta. Our detection of Doppler-broadened atomic emission lines and a bright FIR continuum hints at the presence of newly formed dust in SN ejecta. For G320.4-1.2, we present the first estimate of an ejecta-dust mass of 0.1 - 0.2 M$_\odot$, which spatially coincides with the broad line emission, by applying a blackbody model fit with components of the SNR and background emission. Our sample includes raster maps of 63, 145 $\mu$m [O I] and 158 $\mu$m [C II] lines toward SNRs Kes 79, CTB 109, and IC 443. Based on these line intensities, we suggest interacting shock types in these SNRs. Finally, we compare our LWS spectra of our sample SNRs with the spectra of several HII regions, and discuss their FIR line intensity ratios and continuum properties. Follow-up observations with modern instruments (e.g. $JWST$ and $SOFIA$) with higher spatial and spectral resolution are encouraged for an extensive study of the SN ejecta and the SN dust.

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) started median-resolution spectroscopic (MRS, R$\sim$7500) survey since October 2018. The main scientific goals of MRS, including binary stars, pulsators, and other variable stars are launched with a time-domain spectroscopic survey. However, the systematic errors, including the bias induced from wavelength calibration and the systematic difference between different spectrographs have to be carefully considered during radial velocity measurement. In this work, we provide a technique to correct the systematics in the wavelength calibration based on the relative radial velocity measurements from LAMOST MRS spectra. We show that, for the stars with multi-epoch spectra, the systematic bias which is induced from the exposures of different nights can be well corrected for LAMOST MRS in each spectrograph. And the precision of radial velocity zero-point of multi-epoch time-domain observations reaches below 0.5 km/s . As a by-product, we also give the constant star candidates, which can be the secondary radial-velocity standard star candidates of LAMOST MRS time-domain surveys.

A scenario for the cosmological evolution of self-interacting Bose-Einstein condensed (SIBEC) dark matter (DM) as the final product of a transition from an initial cold DM (CDM)-like phase is considered, motivated by suggestions in the literature that a cold DM gas might have undergone a Bose-Einstein condensate phase transition. The phenomenological model employed for the cold-SIBEC transition introduces three additional parameters to those already present in $\Lambda$CDM; the strength of the DM self-interaction in the SIBEC phase, the time of the transition, and the rate of the transition. Constraints on these extra parameters are obtained from large-scale observables, using the cosmic microwave background (CMB), baryonic acoustic oscillations (BAO) and growth factor measurements, and type Ia supernovae (SNIa) distances. The standard cosmological parameters are found to be unchanged from $\Lambda$CDM, and upper bounds on the SIBEC-DM self-interaction for the various transition times and rates are obtained. If, however, SIBEC-DM is responsible for the tendency of low-mass halos to be cored rather than cuspy, then cold-SIBEC transition times around matter-radiation equality and earlier are ruled out.

The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution must be calibrated using high resolution numerical simulations, We exploit our state of the art graphics processing unit accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a detailed and statistically robust calibration of the VOS model. We rely on the largest set of high resolution simulations carried out to date, with a wide range of cosmological expansion rates, and explore the impact of key numerical parameters, including the dynamic range (comparing box sizes from $1024^3$ to $4096^3$), the lattice spacing, and the choice of numerical estimators for the string velocity. We explore the sensitivity of the VOS model parameters to these numerical parameters, with a particular emphasis on the observationally crucial loop chopping efficiency, and also identify key differences between the equation of state and conjugate momentum estimators for the string velocities, showing that the latter one is more reliable for fast expansion rates (while in Minkowski space the opposite has been previously shown). Finally, we briefly illustrate how our results impact observational constraints on cosmic strings.

In this work, we extend our previous study of the bulk viscosity of hot and dense $npe$ matter induced by the Urca process in the neutrino trapped regime to $npe\mu$ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as $\zeta\sim T^{-2}$ at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range $\gtrsim$ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

One of the new astrophysical phase spaces is exploring sky in the optical (320 nm - 650 nm) range within millisecond to nanosecond timescales known as ultra-fast astronomy (UFA). For this purpose, we developed our own customized readout system for silicon photomultiplier (SiPM) to scan the sky as fast as possible. SiPMs are capable of single-photon detection in the visible light range. The developed readout system for these detectors consists of 16 channels of 14-bit data logging. Each channel includes a 50-dB gain pre-amplifier, signal shaping circuits, analog front end, analogue to digital converter and Xilinx UltraScale+ MPSoC board for data-logging. The results of our readout system show that, scans can be done with 16 ns time frame and a power consumption of 250~mW per channel.

The number of muons in extensive air showers predicted using LHC-tuned hadronic interaction models, such as EPOS-LHC and QGSJetII-04, is smaller than observed in showers recorded by leading cosmic ray experiments. In this paper, we present a new method to derive muon rescaling factors by analyzing reconstructions of simulated showers. The z-variable used (difference of initially simulated and reconstructed total signal in detectors) is connected to the muon signal and is roughly independent of the zenith angle but depends on the mass of primary cosmic ray. The performance of the method is tested using Monte Carlo shower simulations for the hybrid detector of the Pierre Auger Observatory. Having an individual z-value from each simulated hybrid event, the corresponding signal at 1000 m from the shower axis, and using a parametrization of the muon fraction in simulated showers, we can calculate the multiplicative rescaling parameters of the muon signals in the ground detector even for an individual event. We can also study its dependence as a function of zenith angle and the mass of primary cosmic ray. This gives a possibility not only to test/calibrate the hadronic interaction models, but also to derive the $\beta$-exponent, describing an increase of the number of muons as a function of primary energy and mass of the cosmic ray. Detailed simulations show dependence of the $\beta$-exponent on hadronic interaction properties, thus the determination of this parameter is important for understanding the muon deficit problem.

We analyze properties of non-thermal radio emission from the Central Molecular Zone (CMZ) and individual molecular clouds, and argue that the observed features can be interpreted in the framework of our recent theory of self-modulation of cosmic rays (CRs) penetrating dense molecular regions. For clouds with gas column densities of $\sim10^{23}$ cm$^{-2}$, the theory predicts depletion of sub-GeV CR electrons, occurring due to self-modulation of CR protons and leading to harder synchrotron spectra in the sub-GHz range. The predicted imprints of electron depletion in the synchrotron spectra agree well with the spectral hardening seen in available radio observations of the CMZ. A similar, but even stronger effect on the synchrotron emission is predicted for individual (denser) CMZ clouds, such as the Sgr B2. However, the emission at frequencies above $\sim$ GHz, where observational data are available, is completely dominated by the thermal component, and therefore new observations at lower frequencies are needed to verify the predictions.

Radiative transfer calculations in strong (few $\times 10^{12}$ G) magnetic fields, observed in X-ray pulsars, require accurate resonant differential scattering cross sections. Such cross sections exist, but they are quite cumbersome. Here we compare the classical (non-relativistic) with the quantum-mechanical (relativistic) resonant differential scattering cross sections and offer a prescription for the use of the much simpler classical expressions with impressively accurate results. We have expanded the quantum-mechanical differential cross sections and kept terms up to first order in $\epsilon \equiv E/m_ec^2$ and $B \equiv {\cal B}/{\cal B}_{cr}$, where $E$ is the photon energy and ${\cal B}_{cr}$ is the critical magnetic field, and recovered the classical differential cross sections plus terms that are due to spin flip, which is a pure quantum-mechanical phenomenon. Adding by hand the spin-flip terms to the polarization-dependent classical differential cross sections, we find that they are in excellent agreement with the quantum mechanical ones for all energies near resonance and all angles. We have plotted both of them and the agreement is impressive. We give a prescription for the use of the classical differential cross sections that guarantees very accurate results.

Accretion onto magnetic neutron stars results in X-ray spectra that often exhibit a cyclotron resonance scattering feature (CRSF) and, sometimes, higher harmonics of it. Two places are suspect for the formation of a CRSF: the surface of the neutron star and the radiative shock in the accretion column. Here we explore the first possibility: reflection at the neutron-star surface of the continuum produced at the radiative shock. It has been proposed that for high-luminosity sources, as the luminosity increases, the height of the radiative shock increases, thus a larger polar area is illuminated, and as a consequence the energy of the CRSF decreases because the dipole magnetic field decreases by a factor of two from the pole to the equator. We used a Monte Carlo code to compute the reflected spectrum from the atmosphere of a magnetic neutron star, when the incident spectrum is a power-law one. We restricted ourselves to cyclotron energies $\ll m_ec^2$ and used polarization-dependent scattering cross sections, allowing for polarization mode change. As expected, a prominent CRSF is produced in the reflected spectra if the incident photons are in a pencil beam, which hits the neutron-star surface at a point with a well-defined magnetic field strength. However, the incident beam from the radiative shock has a finite width and thus various magnetic field strengths are sampled. As a result of overlap, the reflected spectra have a CRSF, which is close to that produced at the magnetic pole, independent of the height of the radiative shock. Reflection at the surface of a magnetic neutron star cannot explain the observed decrease in the CRSF energy with luminosity in the high-luminosity X-ray pulsar V 0332+53.

KCDC, the 'KASCADE Cosmic-ray Data Centre', is a web-based interface where initially the scientific data from the completed air-shower experiment KASCADE-Grande was made available for the astroparticle community as well as for the interested public. Over the past 7 years, we have continuously extended the data shop with various releases and increased both the number of detector components from the KASCADE-Grande experiment and the data sets and corresponding simulations. With the latest releases we added a new and independent data shop for a specific KASCADE-Grande event selection and by that created the technology for integrating further data shops and data of other experiments, like the data of the air-shower experiment MAKET-ANI in Armenia. In addition, we made available educational examples how to use the data, more than 100 cosmic ray energy spectra from various experiments, and recently attached a public server with access to Jupyter notebooks. In this paper we present a brief history of KCDC, the main features of the recent release as well as will discuss future development plans.

We investigate the origin of intergalactic light (IGL) in close groups of galaxies. IGL is hypothesized to be the byproduct of interaction and merger within compact groups. Comparing the X-ray point source population in our sample of compact groups that have intergalactic light with compact groups without IGL, we find marginal evidence for a small increase in ultra-luminous X-ray sources (ULXs). There is also a significant bias towards lower luminosity high mass X-ray binaries (HMXRBs). We interpret this as an indication that groups with visible IGL represent a later evolutionary phase than other compact groups. They have galaxies characterized by quenching of star formation (lower star formation rate (SFR) inferred from lower HMXRB luminosity) after stellar material has been removed from the galaxies into the intergalactic medium, which is the source of the IGL. We conclude that the presence of an increased fraction of ULXs is due to past interaction and mergers within groups that have IGL.

Context. The Vela complex is a region of the sky that gathers several stellar and interstellar structures in a few hundred square degrees. Aims. Gaia data now allows us to obtain a 3D view of the Vela interstellar structures through the dust extinction. Methods. We used the FEDReD (Field Extinction-Distance Relation Deconvolver) algorithm on near-infrared 2MASS data, cross-matched with the Gaia DR2 catalogue, to obtain a 3D cube of extinction density. We applied the FellWalker algorithm on this cube to locate clumps and dense structures. Results. We analysed 18 million stars on $450~\mathrm{deg}^2$ to obtain the extinction density of the Vela complex from 0.5 to 8~kpc at $\ell\in[250\degr,280\degr]$ and $b\in[$-10$\degr,5\degr]$. This cube reveals the complete morphology of known structures and relations between them. In particular, we show that the Vela Molecular Ridge is more likely composed by three substructures instead of four, as suggested by the 2D densities. These substructures form the shell of a large cavity. This cavity is visually aligned with the Vela supernova remnant but located at a larger distance. We provide a catalogue of location, distance, size and total dust content of ISM clumps that we extracted from the extinction density cube.

We present spectroscopic follow-up observations of 68 red, faint candidates from our multi-epoch, multi-wavelength, previously published survey of NGC 2264. Using near-infrared spectra from VLT/KMOS, we measure spectral types and extinction for 32 young low-mass sources. We confirm 13 as brown dwarfs in NGC 2264, with spectral types between M6 and M8, corresponding to masses between 0.02 and 0.08$M_{\odot}$. These are the first spectroscopically confirmed brown dwarfs in this benchmark cluster. 19 more objects are found to be young M-type stars of NGC 2264 with masses of 0.08 to 0.3$\,M_{\odot}$. 7 of the confirmed brown dwarfs as well as 15 of the M-stars have IR excess caused by a disc. Comparing with isochrones, the typical age of the confirmed brown dwarfs is $<$0.5 to 5Myr. More than half of the newly identified brown dwarfs and very low mass stars have ages $<$0.5Myr, significantly younger than the bulk of the known cluster population. Based on the success rate of our spectroscopic follow-up, we estimate that NGC 2264 hosts 200-600 brown dwarfs in total (in the given mass range). This would correspond to a star-to-brown dwarf ratio between 2.5:1 and 7.5:1. We determine the slope of the substellar mass function as $\alpha = 0.43^{+0.41}_{-0.56}$, these values are consistent with those measured for other young clusters. This points to a uniform substellar mass function across all star forming environments.

Context. Radar scattering from meteor trails depends on several poorly constrained quantities, such as electron line density, q, initial trail radius, r0, and ambipolar diffusion coefficient, D. Aims. The goal is to apply a numerical model of full wave backscatter to triple frequency echo measurements to validate theory and constrain estimates of electron radial distribution, initial trail radius, and the ambipolar diffusion coefficient. Methods. A selection of 50 transversely polarized and 50 parallel polarized echoes with complete trajectory information were identified from simultaneous tri-frequency echoes recorded by the Canadian Meteor Orbit Radar (CMOR). The amplitude-time profile of each echo was fit to our model using three different choices for the radial electron distribution assuming a Gaussian, parabolicexponential, and 1-by-r2 electron line density model. The observations were manually fit by varying, q, r0, and D per model until all three synthetic echo-amplitude profiles at each frequency matched observation. Results. The Gaussian radial electron distribution was the most successful at fitting echo power profiles, followed by the 1-by-r2. We were unable to fit any echoes using a profile where electron density varied from the trail axis as an exponential-parabolic distribution. While fewer than 5% of all examined echoes had self-consistent fits, the estimates of r0 and D as a function of height obtained were broadly similar to earlier studies, though with considerable scatter. Most meteor echoes are found to not be described well by the idealized full wave scattering model.

Imaging atmospheric Cherenkov telescopes, such as the Very Energetic Radiation Imaging Tele-scope Array System (VERITAS), are uniquely suited to resolve the detailed morphology ofextended regions of gamma-ray emission. However, standard VERITAS data analysis techniques have insufficient sensitivity to gamma-ray sources spanning the VERITAS field of view (3.5{\deg}),due to difficulties with background estimation. For analysis of such spatially extended sources with 0.5{\deg} to greater than 2{\deg} radius, we developed the Matched Runs Method. This method derives background estimations for observations of extended sources using matched separate observations of known point sources taken under similar observing conditions. Our technique has been validated by application to archival VERITAS data. Here we present a summary of the Matched RunsMethod and multiple validation studies on different gamma-ray sources using VERITAS data.

The C/O-ratio as traced with C$_2$H emission in protoplanetary disks is fundamental for constraining the formation mechanisms of exoplanets and our understanding of volatile depletion in disks, but current C$_2$H observations show an apparent bimodal distribution which is not well understood, indicating that the C/O distribution is not described by a simple radial dependence. The transport of icy pebbles has been suggested to alter the local elemental abundances in protoplanetary disks, through settling, drift and trapping in pressure bumps resulting in a depletion of volatiles in the surface and an increase of the elemental C/O. We combine all disks with spatially resolved ALMA C$_2$H observations with high-resolution continuum images and constraints on the CO snowline to determine if the C$_2$H emission is indeed related to the location of the icy pebbles. We report a possible correlation between the presence of a significant CO-icy dust reservoir and high C$_2$H emission, which is only found in disks with dust rings outside the CO snowline. In contrast, compact dust disks (without pressure bumps) and warm transition disks (with their dust ring inside the CO snowline) are not detected in C$_2$H, suggesting that such disks may never have contained a significant CO ice reservoir. This correlation provides evidence for the regulation of the C/O profile by the complex interplay of CO snowline and pressure bump locations in the disk. These results demonstrate the importance of including dust transport in chemical disk models, for a proper interpretation of exoplanet atmospheric compositions, and a better understanding of volatile depletion in disks, in particular the use of CO isotopologues to determine gas surface densities.

Dust is observed in the polar regions of nearby AGN and it is known to contribute substantially to their mid-IR emission and to the obscuration of their UV to optical emission. We selected a sample of 1275 BLAGN in the XMM-XXL field, with optical to infrared photometric data. These AGN are seen along their polar direction and we expect a maximal impact of dust located around the poles when it is present. We used X-CIGALE, which introduces a dust component to account for obscuration along the polar directions, modeled as a foreground screen, and an extinction curve that is chosen as it steepens significantly at short wavelengths or is much grayer. By comparing the results of different fits, we are able to define subsamples of sources with positive statistical evidence in favor of or against polar obscuration and described using the gray or steep extinction curve. We find a similar fraction of sources with positive evidence for and against polar dust. Applying statistical corrections, we estimate that half of our sample could contain polar dust and among them, 60% exhibit a steep extinction curve and 40% a flat extinction curve; although these latter percentages are found to depend on the adopted extinction curves. The obscuration in the V-band is not found to correlate with the X-ray column density, while A_V/N_H ratios span a large range of values and higher dust temperatures are found with the flat, rather than with the steep extinction curve. Ignoring this polar dust component in the fit of the spectral energy distribution of these composite systems leads to an overestimation of the stellar contribution. A single fit with a polar dust component described with an SMC extinction curve efficiently overcomes this issue but it fails at identifying all the AGN with polar dust obscuration.

Sub-millimeter emission lines produced by the interstellar medium (ISM) are strong tracers of star formation and are some of the main targets of line intensity mapping (LIM) surveys. In this work we present an empirical multi-line emission model that simultaneously covers the mean, scatter, and correlations of [CII], CO J=1-0 to J=5-4, and [CI] lines in redshift range $1\leq z\leq9$. We assume the galaxy ISM line emission luminosity versus halo mass relations can be described by double power laws with redshift-dependent log normal scatter. The model parameters are then derived by fitting to the state of the art semi-analytic simulation results that have successfully reproduced multiple sub-millimeter line observations at $0\leq z\lesssim6$. We cross check the line emission statistics predicted by the semi-analytic simulation and our empirical model, finding that at $z\geq1$ our model reproduces the simulated line intensities with fractional error less than about 10%. The fractional difference is less than 25% for the power spectra. Grounded on physically-motivated and self-consistent galaxy simulations, this computationally efficient model will be helpful in forecasting ISM emission line statistics for upcoming LIM surveys.

During a tracking mode observation, every telescope with an alt-azimuthal mount shows a rotation in the field of view (FoV) due to the diurnal motion of the Earth. The angular extension of the rotation depends mainly on the time length of the observation, but also on the telescope's latitude and pointing, because it is determined by the evolution of the parallactic angle of the target, which is a function of those two parameters. In many cases, the rotation of the FoV can be exploited to assess some optomechanical properties of the telescope, e.g. the alignment of the optical elements or the motors' precision during the tracking. As a consequence, it could happen that a proper simulation of the FoV rotation is crucial to program an observation aiming at calibrating the whole system. We present a tool to simulate the apparent rotation of the FoV, calculating the actual "star coverage" exploitable for scientific goals. Given the FoV and the pointing direction, the software calculates the angular extension of the rotation, considering only the stars observable by the telescope below the magnitude limit. This tool will be adopted to schedule the pointing calibration runs of the innovative ASTRI-Horn Cherenkov telescope, developed by INAF for gamma-ray ground-based astronomy, but with the potentiality to produce sky images as an ancillary output, using the so-called Variance method. By exploiting the FoV rotation with the Variance method, the critical assessment of the camera axis can be successfully performed.

Magnetic helicity (MH) is a measure of twist and shear of magnetic field. MH is injected in the active region (AR) corona through photospheric footpoint motions causing twisted and sheared magnetic fields. From the conservation property of the helicity, it was conjectured that an already twisted flux rope (FR) with continuous injection of MH inevitably erupts to remove the excess accumulated coronal helicity. Therefore, understanding the nature and evolution of the photospheric helicity flux transfer is crucial to reveal the intensity of the flare/CME activity. Using the time-sequence vector-magnetograms of \textit{Helioseismic Magnetic Imager}, we study the evolution of MH injection in emerging AR 12257. The photospheric flux motions in this AR inject positive helicity in the first 2.5 days followed by negative helicity later. This successive injection of opposite helicity is consistent with the sign of mean force-free twist parameter ($\alpha_{av}$), orientation of magnetic-tongues. Also, the extrapolated AR magnetic structure exhibits transformation of global-shear without a twisted FR in the core of the AR. No CMEs are launched from this AR but C-class flaring activity is observed predominantly in the second half of the evolution period. The ARs with sign reversal of the MH injection are not favorable to twisted FR formation with excess coronal helicity and therefore are important to identify CME-less ARs readily. A possible scenario in these ARs is that when one sign of helicity flux is replaced by opposite sign, the magnetic field of different connectivity with opposite shear undergoes reconnection at different scales giving rise to both intermittent flares and enhanced coronal heating.

Supermassive black holes (SMBHs) are ubiquitously found at the centers of most massive galaxies. Measuring SMBH mass is important for understanding the origin and evolution of SMBHs. However, traditional methods require spectroscopic data which is expensive to gather. We present an algorithm that weighs SMBHs using quasar light time series, circumventing the need for expensive spectra. We train, validate, and test neural networks that directly learn from the Sloan Digital Sky Survey (SDSS) Stripe 82 light curves for a sample of $38,939$ spectroscopically confirmed quasars to map out the nonlinear encoding between SMBH mass and multi-color optical light curves. We find a 1$\sigma$ scatter of 0.37 dex between the predicted SMBH mass and the fiducial virial mass estimate based on SDSS single-epoch spectra, which is comparable to the systematic uncertainty in the virial mass estimate. Our results have direct implications for more efficient applications with future observations from the Vera C. Rubin Observatory. Our code, \textsf{AGNet}, is publicly available at {\color{red} \url{https://github.com/snehjp2/AGNet}}.

We present the results of a high-cadence spectroscopic and imaging monitoring campaign of the active galactic nucleus (AGN) of NGC 4395. High signal-to-noise-ratio spectra were obtained at the Gemini-N 8 m telescope using the GMOS integral field spectrograph (IFS) on 2019 March 7, and at the Keck-I 10 m telescope using the Low-Resolution Imaging Spectrometer (LRIS) with slitmasks on 2019 March 3 and April 2. Photometric data were obtained with a number of 1 m-class telescopes during the same nights. The narrow-line region (NLR) is spatially resolved; therefore, its variable contributions to the slit spectra make the standard procedure of relative flux calibration impractical. We demonstrate that spatially-resolved data from the IFS can be effectively used to correct the slit-mask spectral light curves. While we obtained no reliable lag owing to the lack of strong variability pattern in the light curves, we constrain the broad line time lag to be less than 3 hr, consistent with the photometric lag of $\sim80$ min reported by Woo et al. (2019). By exploiting the high-quality spectra, we measure the second moment of the broad component of the H$\alpha$ emission line to be $586\pm19$ km s$^{-1}$, superseding the lower value reported by Woo et al. (2019). Combining the revised line dispersion and the photometric time lag, we update the black hole mass as $(1.7\pm 0.3)\times10^4$ M$_{\odot}$.

We present fully three-dimensional equations to describe the rotations of a body made of a deformable mantle and a fluid core. The model in its essence is similar to that used by INPOP (Integration Plan\'{e}taire de l'Observatoire de Paris), e.g. Viswanathan et al. (2019), and by JPL (Jet Propulsion Laboratory), e.g. Folkner et al. (2014), to represent the Moon. The intended advantages of our model are: straightforward use of any linear-viscoelastic model for the rheology of the mantle; easy numerical implementation in time-domain (no time lags are necessary); all parameters, including those related to the "permanent deformation", have a physical interpretation. The paper also contains: 1) A physical model to explain the usual lack of hydrostaticity of the mantle (permanent deformation). 2) Formulas for free librations of bodies in and out-of spin-orbit resonance that are valid for any linear viscoelastic rheology of the mantle. 3) Formulas for the offset between the mantle and the idealized rigid-body motion (Peale's Cassini states). 4) Applications to the librations of Moon, Earth, and Mercury that are used for model validation.

Muons from extensive air showers appear as rings in images taken with imaging atmospheric Cherenkov telescopes, such as VERITAS. These muon-ring images are used for the calibration of the VERITAS telescopes, however the calibration accuracy can be improved with a more efficient muon-identification algorithm. Convolutional neural networks (CNNs) are used in many state-of-the-art image-recognition systems and are ideal for muon image identification, once trained on a suitable dataset with labels for muon images. However, by training a CNN on a dataset labelled by existing algorithms, the performance of the CNN would be limited by the suboptimal muon-identification efficiency of the original algorithms. Muon Hunters 2 is a citizen science project that asks users to label grids of VERITAS telescope images, stating which images contain muon rings. Each image is labelled 10 times by independent volunteers, and the votes are aggregated and used to assign a `muon' or `non-muon' label to the corresponding image. An analysis was performed using an expert-labelled dataset in order to determine the optimal vote percentage cut-offs for assigning labels to each image for CNN training. This was optimised so as to identify as many muon images as possible while avoiding false positives. The performance of this model greatly improves on existing muon identification algorithms, identifying approximately 30 times the number of muon images identified by the current algorithm implemented in VEGAS (VERITAS Gamma-ray Analysis Suite), and roughly 2.5 times the number identified by the Hough transform method, along with significantly outperforming a CNN trained on VEGAS-labelled data.

The Universe's initial conditions, in particular baryon and cold dark matter (CDM) isocurvature perturbations, are poorly constrained on sub-Mpc scales. In this paper, we develop a new formalism to compute the effect of small-scale baryon perturbations on the mean free-electron abundance, thus on cosmic microwave background (CMB) anisotropies. Our framework can accommodate perturbations with arbitrary time and scale dependence. We apply this formalism to four different combinations of baryon and CDM isocurvature modes, and use Planck CMB-anisotropy data to probe their initial amplitude. We find that Planck data is consistent with no small-scale isocurvature perturbations, and that this additional ingredient does not help alleviate the Hubble tension. We set upper bounds to the dimensionless initial power spectrum $\Delta_{\mathcal{I}}^2(k)$ of these isocurvature modes at comoving wavenumbers $1~\textrm{Mpc}^{-1} \le k \le 10^3$ Mpc$^{-1}$, for several parameterizations. For a scale-invariant power spectrum, our 95% confidence-level limits on $\Delta_{\mathcal{I}}^2$ are 0.023 for pure baryon isocurvature, 0.099 for pure CDM isocurvature, 0.026 for compensated baryon-CDM perturbations, and 0.009 for joint baryon-CDM isocurvature perturbations. Using a Fisher analysis generalized to non-analytic parameter dependence, we forecast that a CMB Stage-4 experiment would be able to probe small-scale isocurvature perturbations with initial power 3 to 10 times smaller than Planck limits. The formalism introduced in this work is very general and can be used more widely to probe any physical processes or initial conditions sourcing small-scale baryon perturbations.

We numerically construct a symmetric wormhole solution in pure Einstein gravity supported by a massive $3$-form field with a potential that contains a quartic self-interaction term. The wormhole spacetimes have only a single throat and they are everywhere regular and asymptotically flat. Furthermore, their mass and throat circumference increase almost linearly as the coefficient of the quartic self-interaction term $\Lambda$ increases. The amount of violation of the null energy condition (NEC) is proportional to the magnitude of $3$-form, thus the NEC is less violated as $\Lambda$ increases, since the magnitude of $3$-form decreases with $\Lambda$. In addition, we investigate the geodesics of particles moving around the wormhole. The unstable photon orbit is located at the throat. We also find that the wormhole can cast a shadow whose apparent size is smaller than that cast by the Schwarzschild black hole, but reduces to it when $\Lambda$ acquires a large value. The behavior of the innermost stable circular orbit around this wormhole is also discussed. The results of this paper hint toward the possibility that the 3-form wormholes could be potential black hole mimickers, as long as $\Lambda$ is sufficiently large, precisely when NEC is weakly violated.

The acceleration of charged particles by interplanetary shocks can drain a non-negligible fraction of the upstream ram pressure. For a sample of shocks observed in-situ at 1 AU by the ACE and Wind spacecraft, time-series of the non-Maxwellian components (supra-thermal and higher-energy) of the ion and electron energy spectra were acquired for each event. These were averaged for one hour before and after the time of the shock passage to determine their partial pressure. Using the MHD Rankine-Hugoniot jump conditions, we find that the fraction of the total upstream energy flux transferred to non-Maxwellian downstream particles is typically about 2-16%. Notably, our sample shows that neither the fast magnetosonic Mach number nor the angle between the shock normal and average upstream magnetic field are correlated with non-Maxwellian particle pressure.

We consider a modified gravity framework for inflation by adding to the Einstein-Hilbert action a direct $f(\phi)T$ term, where $\phi$ is identified as the inflaton and $T$ is the trace of the energy-momentum tensor. The framework goes to Einstein gravity naturally when inflaton decays out. We investigate inflation dynamics in this $f(\phi)T$ gravity (not to be confused with torsion-scalar coupled theories) on a general basis, and then apply it to three well-motivated inflationary models. We find that the predictions for the spectral tilt and the tensor-to-scalar ratio are sensitive to this new $f(\phi)T$ term. This $f(\phi)T$ gravity brings both chaotic and natural inflation into better agreement with data. For Starobinsky inflation, the coupling constant $\alpha$ in $[-0.0026,0.0031]$ for $N=60$ is in Planck-allowed $2\sigma$ region.

We consider the possibility of the lightest sterile neutrino dark matter which has dipole interaction with heavier sterile neutrinos. The lifetime can be long enough to be a dark matter candidate without violating other constraints and the correct amount of relic abundance can be produced in the early Universe. We find that a sterile neutrino with the mass of around MeV and the dimension-five non-renormalisable dipole interaction suppressed by $\Lambda_5 \gtrsim 10^{15}$ GeV can be a good candidate of dark matter, while heavier sterile neutrinos with masses of the order of GeV can explain the active neutrino oscillations.

We consider the quantum electrodynamic corrections to the two-stream instability. We find these corrections vanish at first order unless a guiding magnetic field $\mathbf{B}_0$ is considered. With respect to the classical version of the instability, quantum electrodynamic effects reduce the most unstable wave vector and its growth rate by a factor $\sqrt{1+\xi}$, with $\xi = \frac{\alpha}{9\pi} (B_0/B_{cr})^2$, where $\alpha$ is the fine-structure constant and $B_{cr}$ the Schwinger critical magnetic field. Although derived for a cold system, these results are valid for the kinetic case. The results are valid in the range $\xi \ll 1$ and, actually, up to linear corrections in $\xi$.

The universe has evolved through several phases as its various constituents dominated its energy content. Candidate dark matter particles may have undergone freeze-out during any such phase. While the standard freeze-out scenarios have been explored during the radiation-dominated era, and more recently during scalar field decay, this work extends the study of dark matter freeze-out to a potential early period during which the universe is matter-dominated and its evolution adiabatic. Decoupling during an adiabatic matter dominated era changes the freeze-out dynamics, since the Hubble rate is parametrically different for matter and radiation domination. Furthermore, for successful Big Bang Nucleosynthesis the state dominating the early universe energy density must decay, this dilutes (or repopulates) the dark matter. As a result, the masses and couplings required to match the observed dark matter relic density can differ signifcantly from radiation dominated freeze-out.