-
Intergalactic medium rotation measure of primordial magnetic fields
Authors:
Salome Mtchedlidze,
Paola Domínguez-Fernández,
Xiaolong Du,
Ettore Carretti,
Franco Vazza,
Shane Patrick O'Sullivan,
Axel Brandenburg,
Tina Kahniashvili
Abstract:
The Faraday rotation effect, quantified by the Rotation Measure (RM), is a powerful probe of the large-scale magnetization of the Universe - tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic Medium (IGM; referred to as $\mathrm{RM}_{\text{IGM}}$). The redshift dependence of the latter has extensively been explored with observations. It has also been…
▽ More
The Faraday rotation effect, quantified by the Rotation Measure (RM), is a powerful probe of the large-scale magnetization of the Universe - tracing magnetic fields not only on galaxy and galaxy cluster scales but also in the intergalactic Medium (IGM; referred to as $\mathrm{RM}_{\text{IGM}}$). The redshift dependence of the latter has extensively been explored with observations. It has also been shown that this relation can help to distinguish between different large-scale magnetization scenarios. We study the evolution of this $\mathrm{RM}_{\text{IGM}}$ for different primordial magnetogenesis scenarios to search for the imprints of primordial magnetic fields (PMFs; magnetic fields originating in the early Universe) on the redshift-dependence of $\mathrm{RM}_{\text{IGM}}$. We use cosmological magnetohydrodynamic (MHD) simulations for evolving PMFs during large-scale structure formation, coupled to the light cone analysis to produce a realistic statistical sample of mock $\mathrm{RM}_{\text{IGM}}$ images. We study the predicted behavior for the cosmic evolution of $\mathrm{RM}_{\text{IGM}}$ for different correlation lengths of PMFs, and provide fitting functions for their dependence on redshifts. We compare these mock RM trends with the recent analysis of the the LOw-Frequency ARray (LOFAR) RM Grid and find that large-scale-correlated PMFs should have (comoving) strengths $\lesssim 0.75$ nanoGauss, if originated during inflation with the scale invariant spectrum and (comoving) correlation length $\sim 19$ cMpc/h or $ \lesssim 30$ nanoGauss if they originated during phase-transition epochs with the comoving correlation length $\sim 1$ cMpc/h. Our findings agree with previous observations and confirm the results of semi-analytical studies, showing that upper limits on the PMF strength decrease as their coherence scales increase.
△ Less
Submitted 23 June, 2024;
originally announced June 2024.
-
JWST Lensed quasar dark matter survey II: Strongest gravitational lensing limit on the dark matter free streaming length to date
Authors:
Ryan E. Keeley,
Anna M. Nierenberg,
Daniel Gilman,
Charles Gannon,
Simon Birrer,
Tommaso Treu,
Andrew J. Benson,
Xiaolong Du,
K. N. Abazajian,
T. Anguita,
V. N. Bennert,
S. G. Djorgovski,
K. K. Gupta,
S. F. Hoenig,
A. Kusenko,
C. Lemon,
M. Malkan,
V. Motta,
L. A. Moustakas,
M. S. H. Oh,
D. Sluse,
D. Stern,
R. H. Wechsler
Abstract:
This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to pertu…
▽ More
This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$α$ forest and Milky Way satellite galaxies.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Tidal evolution of cored and cuspy dark matter halos
Authors:
Xiaolong Du,
Andrew Benson,
Zhichao Carton Zeng,
Tommaso Treu,
Annika H. G. Peter,
Charlie Mace,
Fangzhou Jiang,
Shengqi Yang,
Charles Gannon,
Daniel Gilman,
Anna. M. Nierenberg,
Ethan O. Nadler
Abstract:
The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first cons…
▽ More
The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first consider idealized N-body simulations of subhalos within their host halo, using a generalized mass density profile that describes their properties in a variety of dark matter models at infall, including popular warm, cold, and self-interacting ones. Using these simulations we construct tidal "tracks" for the evolution of subhalos based on their conditions at infall. Second, we use the results of these simulations to build semi-analytic models for tidal effects, including stripping and heating and implement them within the code GALACTICUS. Our semi-analytic models can accurately predict the tidal evolution of both cored and cuspy subhalos, including the bound mass and density profiles, providing a powerful and efficient tool for studying the post-infall properties of subhalos in different dark matter models.
△ Less
Submitted 15 July, 2024; v1 submitted 14 March, 2024;
originally announced March 2024.
-
Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei
Authors:
DRHBc Mass Table Collaboration,
Peng Guo,
Xiaojie Cao,
Kangmin Chen,
Zhihui Chen,
Myung-Ki Cheoun,
Yong-Beom Choi,
Pak Chung Lam,
Wenmin Deng,
Jianmin Dong,
Pengxiang Du,
Xiaokai Du,
Kangda Duan,
Xiaohua Fan,
Wei Gao,
Lisheng Geng,
Eunja Ha,
Xiao-Tao He,
Jinniu Hu,
Jingke Huang,
Kun Huang,
Yanan Huang,
Zidan Huang,
Kim Da Hyung,
Hoi Yat Chan
, et al. (58 additional authors not shown)
Abstract:
The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne…
▽ More
The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $α$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed.
△ Less
Submitted 10 June, 2024; v1 submitted 5 February, 2024;
originally announced February 2024.
-
Convergence Tests of Self-Interacting Dark Matter Simulations
Authors:
Charlie Mace,
Zhichao Carton Zeng,
Annika H. G. Peter,
Xiaolong Du,
Shengqi Yang,
Andrew Benson,
Mark Vogelsberger
Abstract:
Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this process can suffer from numerical errors when simulation parameters are selected incorrectly. Optimal choices for simulation parameters are well studied for cold dar…
▽ More
Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this process can suffer from numerical errors when simulation parameters are selected incorrectly. Optimal choices for simulation parameters are well studied for cold dark matter (CDM), but are not deeply understood when self-interactions are included. In order to perform reliable N-body simulations and model core-collapse accurately we must understand the potential numerical errors, how to diagnose them, and what parameter selections must be made to reduce them. We use the \texttt{Arepo} N-body code to perform convergence tests of core-collapsing SIDM halos across a range of halo concentrations and SIDM cross-sections, and quantify potential numerical issues related to mass resolution, timestep size, and gravitational softening length. Our tests discover that halos with fewer than $10^5$ simulation particles, a resolution typically not met by subhalos in N-body simulations, suffer from significant discreteness noise that leads to variation and extreme outliers in the collapse rate. At our lowest resolution of $N=10^4$ particles, this collapse time variation can reach as high as 20\%. At this low resolution we also find a bias in collapse times and a small number of extreme outliers. Additionally, we find that simulations which run far beyond the age of the Universe, which have been used to calibrate SIDM gravothermal fluid models in previous work, have a sensitivity to the timestep size that is not present in shorter simulations or simulations using only CDM. Our work shows that choices of simulation parameters that yield converged results for some halo masses and SIDM models do not necessarily yield convergence for others.
△ Less
Submitted 12 April, 2024; v1 submitted 2 February, 2024;
originally announced February 2024.
-
Till the core collapses: the evolution and properties of self-interacting dark matter subhalos
Authors:
Zhichao Carton Zeng,
Annika H. G. Peter,
Xiaolong Du,
Shengqi Yang,
Andrew Benson,
Francis-Yan Cyr-Racine,
Fangzhou Jiang,
Charlie Mace,
R. Benton Metcalf
Abstract:
One of the hottest questions in the cosmology of self-interacting dark matter (SIDM) is whether scatterings can induce detectable core-collapse in halos by the present day. Because gravitational tides can accelerate core-collapse, the most promising targets to observe core-collapse are satellite galaxies and subhalo systems. However, simulating small subhalos is computationally intensive, especial…
▽ More
One of the hottest questions in the cosmology of self-interacting dark matter (SIDM) is whether scatterings can induce detectable core-collapse in halos by the present day. Because gravitational tides can accelerate core-collapse, the most promising targets to observe core-collapse are satellite galaxies and subhalo systems. However, simulating small subhalos is computationally intensive, especially when subhalos start to core-collapse. In this work, we present a hierarchical framework for simulating a population of SIDM subhalos, which reduces the computation time to linear order in the total number of subhalos. With this method, we simulate substructure lensing systems with multiple velocity-dependent SIDM models, and show how subhalo evolution depends on the SIDM model, subhalo mass and orbits. We find that an SIDM cross section of $\gtrsim 200$ cm$^2$/g at velocity scales relevant for subhalos' internal heat transfer is needed for a significant fraction of subhalos to core-collapse in a typical lens system at redshift $z=0.5$, and that core-collapse has unique observable features in lensing. We show quantitatively that core-collapse in subhalos is typically accelerated compared to field halos, except when the SIDM cross section is non-negligible ($\gtrsim \mathcal{O}(1)$ cm$^2$/g) at subhalos' orbital velocities, in which case evaporation by the host can delay core-collapse. This suggests that substructure lensing can be used to probe velocity-dependent SIDM models, especially if line-of-sight structures (field halos) can be distinguished from lens-plane subhalos. Intriguingly, we find that core-collapse in subhalos can explain the recently reported ultra-steep density profiles of substructures found by lensing with the \emph{Hubble Space Telescope}
△ Less
Submitted 4 November, 2023; v1 submitted 15 October, 2023;
originally announced October 2023.
-
JWST lensed quasar dark matter survey I: Description and First Results
Authors:
A. M. Nierenberg,
R. E. Keeley,
D. Sluse,
D. Gilman,
S. Birrer,
T. Treu,
K. N. Abazajian,
T. Anguita,
A. J. Benson,
V. N. Bennert,
S. G. Djorgovski,
X. Du,
C. D. Fassnacht,
S. F. Hoenig,
A. Kusenko,
C. Lemon,
M. Malkan,
V. Motta,
L. A. Moustakas,
D. Stern,
R. H. Wechsler
Abstract:
The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light…
▽ More
The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$μ$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
A comprehensive model for the formation and evolution of the faintest Milky Way dwarf satellites
Authors:
Niusha Ahvazi,
Andrew Benson,
Laura V. Sales,
Ethan O. Nadler,
Sachi Weerasooriya,
Xiaolong Du,
Mia Sauda Bovill
Abstract:
In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the…
▽ More
In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size$-$mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy-halo connection. We predict a total of $300 ^{+75} _{-99}$ satellites with an absolute $V$-band magnitude (M$_{V}$) less than $0$ within $300$ kpc from our Milky Way-analogs. The fraction of subhalos that host a galaxy at least this bright drops to $50\%$ by a halo peak mass of $\sim 8.9 \times 10^{7}$ M$_{\odot}$, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population.
△ Less
Submitted 18 March, 2024; v1 submitted 25 August, 2023;
originally announced August 2023.
-
Change of rotation measure during eclipse of a black widow PSR J2051$-$0827
Authors:
S. Q. Wang,
J. B. Wang,
D. Z. Li,
J. M. Yao,
R. N. Manchester,
G. Hobbs,
N. Wang,
S. Dai,
H. Xu,
R. Luo,
Y. Feng,
W. Y. Wang,
D. Li,
Y. W. Yu,
Z. X. Du,
C. H. Niu,
S. B. Zhang,
C. M. Zhang
Abstract:
Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051$-$0827 using observations made with the Five-hundred-meter Ape…
▽ More
Black widows are millisecond pulsars ablating their companions. The material blown from the companion blocks the radio emission, resulting in radio eclipses. The properties of the eclipse medium are poorly understood. Here, we present direct evidence of the existence of magnetic fields in the eclipse medium of the black widow PSR J2051$-$0827 using observations made with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We detect a regular decrease in rotation measure (RM) in the egress of eclipse, changing from $60\,\rm rad\,m^{-2}$ to $-28.7\,\rm rad\,m^{-2}$. The RM gradually changes back to normal when the line-of-sight moves away from the eclipse. The estimated line-of-sight magnetic field strength in the eclipse medium is $\sim 0.1$ G. The RM reversal could be caused by a change of the magnetic field strength along the line of sight due to binary orbital motion. The RM reversal phenomenon has also been observed in some repeating fast radio bursts (FRBs), and the study of spider pulsars may provide additional information about the origin of FRBs.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
The MOSDEF-LRIS Survey: Detection of Inflowing Gas Towards Three Star-forming Galaxies at z ~ 2
Authors:
Andrew Weldon,
Naveen A. Reddy,
Michael W. Topping,
Alice E. Shapley,
Xinnan Du,
Sedona H. Price,
Ryan L. Sanders,
Alison L. Coil,
Bahram Mobasher,
Mariska Kriek,
Brian Siana,
Saeed Rezaee
Abstract:
We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left<z\right>\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated,…
▽ More
We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left<z\right>\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
A quantitative comparison between velocity dependent SIDM cross sections constrained by the gravothermal and isothermal models
Authors:
Shengqi Yang,
Fangzhou Jiang,
Andrew Benson,
Yi-Ming Zhong,
Charlie Mace,
Xiaolong Du,
Zhichao Carton Zeng,
Annika H. G. Peter,
Moritz S. Fischer
Abstract:
One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these proc…
▽ More
One necessary step for probing the nature of self-interacting dark matter (SIDM) particles with astrophysical observations is to pin down any possible velocity dependence in the SIDM cross section. Major challenges for achieving this goal include eliminating, or mitigating, the impact of the baryonic components and tidal effects within the dark matter halos of interest -- the effects of these processes can be highly degenerate with those of dark matter self-interactions at small scales. In this work we select 9 isolated galaxies and brightest cluster galaxies (BCGs) with baryonic components small enough such that the baryonic gravitational potentials do not significantly influence the halo gravothermal evolution processes. We then constrain the parameters of Rutherford and Moller scattering cross section models with the measured rotation curves and stellar kinematics through the gravothermal fluid formalism and isothermal method. Cross sections constrained by the two methods are consistent at $1σ$ confidence level, but the isothermal method prefers cross sections greater than the gravothermal approach constraints by a factor of $\sim3$.
△ Less
Submitted 26 April, 2024; v1 submitted 8 May, 2023;
originally announced May 2023.
-
Gravitational Bose-Einstein Condensation of Vector/Hidden Photon Dark Matter
Authors:
Jiajun Chen,
Xiaolong Du,
Mingzhen Zhou,
Andrew Benson,
David J. E. Marsh
Abstract:
We study the gravitational Bose-Einstein condensation of a massive vector field in the kinetic regime and the non-relativistic limit using non-linear dynamical numerical methods. Gravitational condensation leads to the spontaneous formation of solitons. We measure the condensation time and growth rate, and compare to analytical models in analogy to the scalar case. We find that the condensation ti…
▽ More
We study the gravitational Bose-Einstein condensation of a massive vector field in the kinetic regime and the non-relativistic limit using non-linear dynamical numerical methods. Gravitational condensation leads to the spontaneous formation of solitons. We measure the condensation time and growth rate, and compare to analytical models in analogy to the scalar case. We find that the condensation time of the vector field depends on the correlation between its different components. For fully correlated configurations, the condensation time is the same as that for a scalar field. On the other hand, uncorrelated or partially correlated configurations condense slower than the scalar case. As the vector soliton grows, it can acquire a net spin angular momentum even if the total spin angular momentum of the initial conditions is zero.
△ Less
Submitted 25 May, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
Axion Star Explosions: A New Source for Axion Indirect Detection
Authors:
Miguel Escudero,
Charis Kaur Pooni,
Malcolm Fairbairn,
Diego Blas,
Xiaolong Du,
David J. E. Marsh
Abstract:
If dark matter is composed of axions, then axion stars form in the cores of dark matter halos. These stars are unstable above a critical mass, decaying to radio photons that heat the intergalactic medium, offering a new channel for axion indirect detection. We recently provided the first accurate calculation of the axion decay rate due to axion star mergers. In this work we show how existing data…
▽ More
If dark matter is composed of axions, then axion stars form in the cores of dark matter halos. These stars are unstable above a critical mass, decaying to radio photons that heat the intergalactic medium, offering a new channel for axion indirect detection. We recently provided the first accurate calculation of the axion decay rate due to axion star mergers. In this work we show how existing data concerning the CMB optical depth leads to strong constraints on the axion photon coupling in the mass range $10^{-14}\,{\rm eV}\lesssim m_a\lesssim 10^{-8}\,{\rm eV}$. Axion star decays lead to efficient reionization of the intergalactic medium during the dark ages. By comparing this non-standard reionization with Planck legacy measurements of the Thompson optical width, we show that couplings in the range $10^{-14}\,{\rm GeV}^{-1} \lesssim g_{aγγ} \lesssim 10^{-10}\,{\rm GeV}^{-1}$ are excluded for our benchmark model of axion star abundance. Future measurements of the 21cm emission of neutral hydrogen at high redshift could improve this limit by an order of magnitude or more, providing complementary indirect constraints on axion dark matter in parameter space also targeted by direct detection haloscopes.
△ Less
Submitted 8 February, 2024; v1 submitted 20 February, 2023;
originally announced February 2023.
-
Soliton Merger Rates and Enhanced Axion Dark Matter Decay
Authors:
Xiaolong Du,
David J. E. Marsh,
Miguel Escudero,
Andrew Benson,
Diego Blas,
Charis Kaur Pooni,
Malcolm Fairbairn
Abstract:
Solitons are observed to form in simulations of dark matter (DM) halos consisting of bosonic fields. We use the extended Press-Schechter formalism to compute the mass function of solitons, assuming various forms for the relationship between halo mass and soliton mass. We further provide a new calculation of the rate of soliton major mergers. Solitons composed of axion DM are unstable above a criti…
▽ More
Solitons are observed to form in simulations of dark matter (DM) halos consisting of bosonic fields. We use the extended Press-Schechter formalism to compute the mass function of solitons, assuming various forms for the relationship between halo mass and soliton mass. We further provide a new calculation of the rate of soliton major mergers. Solitons composed of axion DM are unstable above a critical mass, and decay to either relativistic axions or photons, depending on the values of the coupling constants. We use the computed soliton major merger rate to predict the enhanced DM decay rate due to soliton instability. For certain values of currently allowed axion parameters, the energy injection into the intergalactic medium from soliton decays to photons is comparable to or larger than the energy injection due to core collapse supernovae at $z>10$. A companion paper explores the phenomenology of such an energy injection.
△ Less
Submitted 12 February, 2024; v1 submitted 23 January, 2023;
originally announced January 2023.
-
Growing the First Galaxies' Merger Trees
Authors:
Ethan O. Nadler,
Andrew Benson,
Trey Driskell,
Xiaolong Du,
Vera Gluscevic
Abstract:
Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger tree…
▽ More
Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger trees with branches that are guaranteed to achieve a desired halo mass at a chosen redshift. This method augments extended Press Schechter theory solutions with constrained random processes known as Brownian bridges and is implemented in the open-source semi-analytic model $\texttt{Galacticus}$. We generate ensembles of constrained merger trees to predict the growth histories of seven high-redshift JWST galaxy candidates, finding that these systems most likely merge $\approx 2~\mathrm{Gyr}$ after the observation epoch and occupy haloes of mass $\gtrsim 10^{14}~M_{\mathrm{\odot}}$ today. These calculations are thousands of times more efficient than existing methods, are analytically controlled, and provide physical insights into the evolution of haloes with rapid early growth. Our constrained merger tree implementation is publicly available at http://github.com/galacticusorg/galacticus.
△ Less
Submitted 23 March, 2023; v1 submitted 16 December, 2022;
originally announced December 2022.
-
Inflationary and phase-transitional primordial magnetic fields in galaxy clusters
Authors:
Salome Mtchedlidze,
Paola Domínguez-Fernández,
Xiaolong Du,
Wolfram Schmidt,
Axel Brandenburg,
Jens Niemeyer,
Tina Kahniashvili
Abstract:
Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs…
▽ More
Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs, in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve the turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale inflationary models, while the phase transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as the amplitude and the shape of the magnetic power spectrum) and therefore on the final correlation length. We find a one order of magnitude difference between the final strengths of the inflation- and phase transition-generated magnetic fields, and a factor of 1.5 difference between their final coherence scales. Thus, the final configuration of the magnetic field retains information about the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios.
△ Less
Submitted 16 February, 2023; v1 submitted 18 October, 2022;
originally announced October 2022.
-
Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass
Authors:
Ethan O. Nadler,
Philip Mansfield,
Yunchong Wang,
Xiaolong Du,
Susmita Adhikari,
Arka Banerjee,
Andrew Benson,
Elise Darragh-Ford,
Yao-Yuan Mao,
Sebastian Wagner-Carena,
Risa H. Wechsler,
Hao-Yi Wu
Abstract:
We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-anal…
▽ More
We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-analog ($10^{13}~M_{\mathrm{\odot}}$) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of $\approx 3\times 10^{-7}$ times the host virial mass and a Plummer-equivalent gravitational softening length of $\approx 9\times 10^{-4}$ times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are $\approx 50\%$ higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts' underlying dark matter distributions. We compare our results to the semianalytic model $\mathrm{\texttt{Galacticus}}$, which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We use $\mathrm{\texttt{UniverseMachine}}$ to model halo and subhalo star formation histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available at http://web.stanford.edu/group/gfc/symphony.
△ Less
Submitted 16 March, 2023; v1 submitted 6 September, 2022;
originally announced September 2022.
-
Exploring the Correlation between $\rm{H}α$-to-UV Ratio and Burstiness for Typical Star-forming Galaxies at $z\sim2$
Authors:
Saeed Rezaee,
Naveen A. Reddy,
Michael W. Topping,
Irene Shivaei,
Alice E. Shapley,
Tara Fetherolf,
Mariska Kriek,
Alison Coil,
Bahram Mobasher,
Brian Siana,
Xinnan Du,
Ali Ahmad Khostovan,
Andrew Weldon,
Najmeh Emami,
Nima Chartab
Abstract:
The $\rm{H}α$-to-UV luminosity ratio ($L(\rm Hα)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $Σ_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm Hα)/L(\rm UV)$ for a sample o…
▽ More
The $\rm{H}α$-to-UV luminosity ratio ($L(\rm Hα)/L(\rm UV)$) is often used to probe SFHs of star-forming galaxies and it is important to validate it against other proxies for burstiness. To address this issue, we present a statistical analysis of the resolved distribution of $Σ_{\rm{SFR}}$ as well as stellar age and their correlations with the globally measured $L(\rm Hα)/L(\rm UV)$ for a sample of 310 star-forming galaxies in two redshift bins of $1.37 < z < 1.70$ and $ 2.09 < z < 2.61$ observed by the MOSDEF survey. We use the multi-waveband CANDELS/3D-HST imaging of MOSDEF galaxies to construct $Σ_{\rm{SFR}}$ and stellar age maps. We analyze the composite rest-frame far-UV spectra of a subsample of MOSDEF targets obtained by the Keck/LRIS, which includes 124 star-forming galaxies (MOSDEF-LRIS) at redshifts $1.4 < z < 2.6$, to examine the average stellar population properties, and the strength of age-sensitive FUV spectral features in bins of $L(\rm Hα)/L(\rm UV)$. Our results show no significant evidence that individual galaxies with higher $L(\rm Hα)/L(\rm UV)$ are undergoing a burst of star formation based on the resolved distribution of $Σ_{\rm{SFR}}$ of individual star-forming galaxies. We segregate the sample into subsets with low and high $L(\rm Hα)/L(\rm UV)$. The high-$L(\rm Hα)/L(\rm UV)$ subset exhibits, on average, an age of $\log[\rm{Age/yr}]$ = 8.0, compared to $\log[\rm{Age/yr}]$ = 8.4 for the low-$L(\rm Hα)/L(\rm UV)$ galaxies, though the difference in age is significant at only the $2σ$ level. Furthermore, we find no variation in the strengths of Siiv$λ\lambda1393, 1402$ and Civ$λ\lambda1548, 1550$ P-Cygni features from massive stars between the two subsamples.
△ Less
Submitted 19 September, 2023; v1 submitted 25 August, 2022;
originally announced August 2022.
-
Deep Near-Infrared Survey Toward the W40 and Serpens South Region in Aquila Rift: A Comprehensive Catalog of Young Stellar Objects
Authors:
Jia Sun,
Robert A. Gutermuth,
Hongchi Wang,
Miaomiao Zhang,
Shuinai Zhang,
Yuehui Ma,
Xinyu Du,
Min Long
Abstract:
Active star forming regions are excellent laboratories for studying the origins and evolution of young stellar object (YSO) clustering. The W40 - Serpens South region is such a region, and we compile a super near-and-mid-infrared catalog of point sources in it, based on deep NIR observations of CFHT in combination with 2MASS, UKIDSS, and Spitzer catalogs. From this catalog, we identify 832 YSOs, a…
▽ More
Active star forming regions are excellent laboratories for studying the origins and evolution of young stellar object (YSO) clustering. The W40 - Serpens South region is such a region, and we compile a super near-and-mid-infrared catalog of point sources in it, based on deep NIR observations of CFHT in combination with 2MASS, UKIDSS, and Spitzer catalogs. From this catalog, we identify 832 YSOs, and classify 15, 135, 647, and 35 of them to be the deeply embedded sources, Class I, Class II YSOs, and transition disk sources, respectively. In general, these YSOs are well correlated with the filamentary structures of molecular clouds, especially the deeply embedded sources and the Class I YSOs. The W40 central region is dominated by Class II YSOs, but in the Serpens South region, a half of the YSOs are Class I. We further generate a minimum spanning tree (MST) for all the YSOs. Around the W40 cluster, there are eight prominent MST branches that may trace vestigial molecular gas filaments that once fed gas to the central natal gas clump. Of the eight, only two now include detectable filamentary gas in Herschel data and corresponding Class I YSOs, while the other six are exclusively populated with Class II. Four MST branches overlap with the Serpens South main filament, and where they intersect, molecular gas "hubs" and more Class I YSOs are found. Our results imply a mixture of YSO distributions composed of both primordial and somewhat evolved YSOs in this star forming region.
△ Less
Submitted 18 July, 2022;
originally announced July 2022.
-
A semi-analytic study of self-interacting dark-matter haloes with baryons
Authors:
Fangzhou Jiang,
Andrew Benson,
Philip F. Hopkins,
Oren Slone,
Mariangela Lisanti,
Manoj Kaplinghat,
Annika H. G. Peter,
Zhichao Carton Zeng,
Xiaolong Du,
Shengqi Yang,
Xuejian Shen
Abstract:
We combine the isothermal Jeans model and the model of adiabatic halo contraction into a simple semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. We show that the model agrees well with cosmological SIDM simulations over the entire core-forming stage and up to the onset of gravother…
▽ More
We combine the isothermal Jeans model and the model of adiabatic halo contraction into a simple semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. We show that the model agrees well with cosmological SIDM simulations over the entire core-forming stage and up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable link in the context of the structural diversity of bright dwarf galaxies. The fast speed of the method facilitates analyses that would be challenging for numerical simulations -- notably, 1) we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, $M_{\rm b}/M_{\rm vir}$, and galaxy compactness, $r_{1/2}/R_{\rm vir}$; 2) we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and 3) we delineate the regime of gravothermal core-collapse in the $M_{\rm b}/M_{\rm vir}$-$r_{1/2}/R_{\rm vir}$ space, for a given cross section and a given halo concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible future improvements of the fluid model. We have made our programs for the model publicly available at https://github.com/JiangFangzhou/SIDM.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
Quantum fluctuations masquerade as halos: Bounds on ultra-light dark matter from quadruply-imaged quasars
Authors:
Alexander Laroche,
Daniel Gilman,
Xinyu Li,
Jo Bovy,
Xiaolong Du
Abstract:
Ultra-light dark matter (ULDM) refers to a class of theories, including ultra-light axions, in which particles with mass $m_ψ < 10^{-20}\, \rm{eV}$ comprise a significant fraction of the dark matter. A galactic scale de Broglie wavelength distinguishes these theories from cold dark matter (CDM), suppressing the overall abundance of structure on sub-galactic scales, and producing wave-like interfer…
▽ More
Ultra-light dark matter (ULDM) refers to a class of theories, including ultra-light axions, in which particles with mass $m_ψ < 10^{-20}\, \rm{eV}$ comprise a significant fraction of the dark matter. A galactic scale de Broglie wavelength distinguishes these theories from cold dark matter (CDM), suppressing the overall abundance of structure on sub-galactic scales, and producing wave-like interference phenomena in the density profiles of halos. With the aim of constraining the particle mass, we analyze the flux ratios in a sample of eleven quadruple-image strong gravitational lenses. We account for the suppression of the halo mass function and concentration-mass relation predicted by ULDM theories, and the wave-like fluctuations in the host halo density profile, calibrating the model for the wave interference against numerical simulations of galactic-scale halos. We show that the granular structure of halo density profiles, in particular, the amplitude of the fluctuations, significantly impacts image flux ratios, and therefore inferences on the particle mass derived from these data. We infer relative likelihoods of CDM to ULDM of 8:1, 7:1, 6:1, and 4:1 for particle masses $\log_{10}(m_ψ/\rm{eV})\in[-22.5,-22.25], [-22.25,-22.0],[-22.0,-21.75], [-21.75,-21.5]$, respectively. Repeating the analysis and omitting fluctuations associated with the wave interference effects, we obtain relative likelihoods of CDM to ULDM with a particle mass in the same ranges of 98:1, 48:1, 26:1 and 18:1, highlighting the significant perturbation to image flux ratios associated with the fluctuations. Nevertheless, our results disfavor the lightest particle masses with $m_ψ < 10^{-21.5}\,\rm{eV}$, adding to mounting pressure on ultra-light axions as a viable dark matter candidate.
△ Less
Submitted 15 September, 2022; v1 submitted 22 June, 2022;
originally announced June 2022.
-
Tidal Tracks and Artificial Disruption of Cold Dark Matter Halos
Authors:
Andrew J. Benson,
Xiaolong Du
Abstract:
We describe a simple extension to existing models for the tidal heating of dark matter subhalos which takes into account second order terms in the impulse approximation for tidal heating. We show that this revised model can accurately match the tidal tracks along which subhalos evolve as measured in high-resolution N-body simulations. We further demonstrate that, when a constant density core is in…
▽ More
We describe a simple extension to existing models for the tidal heating of dark matter subhalos which takes into account second order terms in the impulse approximation for tidal heating. We show that this revised model can accurately match the tidal tracks along which subhalos evolve as measured in high-resolution N-body simulations. We further demonstrate that, when a constant density core is introduced into a subhalo, this model is able to quantitatively reproduce the evolution and artificial disruption of N-body subhalos arising from finite resolution effects. Combining these results we confirm prior work indicating that artificial disruption in N-body simulations can result in a factor two underestimate of the subhalo mass function in the inner regions of host halos, and a 10--20% reduction over the entire virial volume.
△ Less
Submitted 22 September, 2022; v1 submitted 3 June, 2022;
originally announced June 2022.
-
Gravothermal solutions of SIDM halos: mapping from constant to velocity-dependent cross section
Authors:
Shengqi Yang,
Xiaolong Du,
Zhichao Carton Zeng,
Andrew Benson,
Fangzhou Jiang,
Ethan O. Nadler,
Annika H. G. Peter
Abstract:
The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the const…
▽ More
The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the constant-cross-section gravothermal solution to models with arbitrary velocity dependence in the cross section. We show that the gravothermal solutions after mapping are in good agreement with Arepo N-Body simulation results. We illustrate the power of this approach by applying this fast mapping method to a halo hosting a low surface brightness galaxy UGC 128. We show that this fast mapping method can be used to constrain free parameters in a physically motivated cross-section model and illustrate parameter space favored by the rotation curve measurement.
△ Less
Submitted 2 March, 2023; v1 submitted 5 May, 2022;
originally announced May 2022.
-
The MOSDEF-LRIS Survey: Connection between Galactic-scale Outflows and the Properties of $z$~2 star-forming Galaxies
Authors:
Andrew Weldon,
Naveen Reddy,
Michael Topping,
Alice Shapley,
Ryan Sanders,
Xinnan Du,
Sedona Price,
Alison Coil,
Brian Siana,
Bahram Mobasher,
Tara Fetherolf,
Irene Shivaei,
Saeed Rezaee
Abstract:
We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at $z$~2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionisation interstellar (LIS) metal absorption lin…
▽ More
We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at $z$~2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionisation interstellar (LIS) metal absorption lines and Ly$α$ emission. Outflow velocities are calculated from the centroids of the LIS absorption and/or Ly$α$ emission, as well as the highest-velocity component of the outflow from the blue wings of the LIS absorption lines. Outflow velocities are found to be marginally correlated or independent of galaxy properties, such as star-formation rate (SFR) and star-formation rate surface density ($Σ_{\rm SFR}$). Outflow velocity scales with SFR as a power-law with index 0.24, which suggests that the outflows may be primarily driven by mechanical energy generated by supernovae explosions, as opposed to radiation pressure acting on dusty material. On the other hand, outflow velocity and $Σ_{\rm SFR}$ are not significantly correlated, which may be due to the limited dynamic range of $Σ_{\rm SFR}$ probed by our sample. The relationship between outflow velocity and $Σ_{\rm SFR}$ normalised by stellar mass ($Σ_{\rm sSFR}$), as a proxy for gravitational potential, suggests that strong outflows (e.g., > 200 km s$^{-1}$) appear ubiquitous above a threshold of log($Σ_{\rm sSFR}/\rm{yr}^{-1}\ \rm{kpc}^{-2}$) ~ -11.3, and that above this threshold, outflow velocity uncouples from $Σ_{\rm sSFR}$. These results highlight the need for higher resolution spectroscopic data and spatially resolved imaging to test the driving mechanisms of outflows predicted by theory.
△ Less
Submitted 17 March, 2022;
originally announced March 2022.
-
Removing the divergence of Chandrasekhar limit caused by generalized uncertainty principle
Authors:
Xin-Dong Du,
Chao-Yun Long
Abstract:
The usual generalized uncertainty principle will lead to a divergent mass limit of white dwarf, and this divergence should be prevented for both scenarios including positive and negative parameters of generalized uncertainty principle. Although it has been shown that negative parameter can directly restore the mass limit, the underlying reason is not given to explain why the negative sign appears…
▽ More
The usual generalized uncertainty principle will lead to a divergent mass limit of white dwarf, and this divergence should be prevented for both scenarios including positive and negative parameters of generalized uncertainty principle. Although it has been shown that negative parameter can directly restore the mass limit, the underlying reason is not given to explain why the negative sign appears under the condition of white dwarf. In order to solve this problem, we derive a field-dependent parameter expression whose sign can change depending on the species of spin fields. Besides, we find that the actual physical effect of the negative sign is aimed at limiting the exorbitant uncertainty of momentum.
△ Less
Submitted 26 August, 2022; v1 submitted 12 January, 2022;
originally announced January 2022.
-
Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum: I. even-even nuclei
Authors:
DRHBc Mass Table Collaboration,
Kaiyuan Zhang,
Myung-Ki Cheoun,
Yong-Beom Choi,
Pooi Seong Chong,
Jianmin Dong,
Zihao Dong,
Xiaokai Du,
Lisheng Geng,
Eunja Ha,
Xiao-Tao He,
Chan Heo,
Meng Chit Ho,
Eun Jin In,
Seonghyun Kim,
Youngman Kim,
Chang-Hwan Lee,
Jenny Lee,
Hexuan Li,
Zhipan Li,
Tianpeng Luo,
Jie Meng,
Myeong-Hwan Mun,
Zhongming Niu,
Cong Pan
, et al. (22 additional authors not shown)
Abstract:
Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding ener…
▽ More
Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding energies, two-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. The rms deviation from the 637 mass data is 1.518 MeV, providing one of the best microscopic descriptions for nuclear masses. The drip lines obtained from DRHBc calculations are compared with other calculations, including the spherical relativistic continuum Hartree-Bogoliubov (RCHB) and triaxial relativistic Hartree-Bogoliubov (TRHB) calculations with PC-PK1. The deformation and continuum effects on the limits of the nuclear landscape are discussed. Possible peninsulas consisting of bound nuclei beyond the two-neutron drip line are predicted. The systematics of the two-nucleon separation energies, two-nucleon gaps, rms radii, quadrupole deformations, potential energy curves, neutron densities, neutron mean-field potentials, and pairing energies in the DRHBc calculations are also discussed. In addition, the $α$ decay energies extracted are in good agreement with available data.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo
Authors:
Zhichao Carton Zeng,
Annika H. G. Peter,
Xiaolong Du,
Andrew Benson,
Stacy Kim,
Fangzhou Jiang,
Francis-Yan Cyr-Racine,
Mark Vogelsberger
Abstract:
Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, follow…
▽ More
Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $\lesssim 10^8 M_\odot$ subhalos in lens-mass systems ($\sim 10^{13}M_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhalos with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhalos ($\lesssim 1/1000$ host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhalos to core-collapse, finding that it is nearly impossible to drive core-collapse in subhalos in SIDM models with constant cross sections. Any discovery of ultra-compact dark substructures with future substructure lensing observations favors additional degrees of freedom, such as velocity-dependence, in the cross section.
△ Less
Submitted 22 August, 2022; v1 submitted 1 October, 2021;
originally announced October 2021.
-
Evolution of primordial magnetic fields during large-scale structure formation
Authors:
Salome Mtchedlidze,
Paola Domínguez-Fernández,
Xiaolong Du,
Axel Brandenburg,
Tina Kahniashvili,
Shane O'Sullivan,
Wolfram Schmidt,
Marcus Brüggen
Abstract:
Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homo…
▽ More
Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homogeneous) or as (ii) scale-invariant stochastic fields, and phase transition-generated ones either as (iii) helical or as (iv) non-helical fields from the radiation-dominated epoch. We find that the final distribution of magnetic fields in the simulated cosmic web shows a dependence on the initial strength and the topology of the seed field. Thus, the observed field configuration retains information on the initial conditions at the moment of the field generation. If detected, primordial magnetic field observations would open a new window for indirect probes of the early universe. The differences between the competing models are revealed on the scale of galaxy clusters, bridges, as well as filaments and voids. The distinctive spectral evolution of different seed fields produces imprints on the correlation length today. We discuss how the differences between rotation measures from highly ionized regions can potentially be probed with forthcoming surveys.
△ Less
Submitted 20 April, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
Relaxation times for Bose-Einstein condensation by self-interaction and gravity
Authors:
Jiajun Chen,
Xiaolong Du,
Erik Lentz,
David J. E. Marsh
Abstract:
In this paper, we study the Bose-Einstein condensation of a scalar field with an attractive self-interaction, with or without gravitational interactions. We confirm through full dynamical simulation that the condensation timescale due to self-interaction is inversely proportional to the square of the number density $n$ and the self-coupling constant $g$ : $τ\propto n^{-2} g^{-2}$. We also investig…
▽ More
In this paper, we study the Bose-Einstein condensation of a scalar field with an attractive self-interaction, with or without gravitational interactions. We confirm through full dynamical simulation that the condensation timescale due to self-interaction is inversely proportional to the square of the number density $n$ and the self-coupling constant $g$ : $τ\propto n^{-2} g^{-2}$. We also investigate the condensation timescale when self-interaction and gravity are both important by solving the Gross-Pitaevskii-Poisson equations, and find that the condensation time scales according to an additive model for the cross section. We discuss the relevance of our results to theoretical models of boson star formation by condensation.
△ Less
Submitted 28 June, 2022; v1 submitted 23 September, 2021;
originally announced September 2021.
-
Diluted Axion Star Collisions with Neutron Stars
Authors:
Yang Bai,
Xiaolong Du,
Yuta Hamada
Abstract:
Diluted axion star, a self-gravitating object with the quantum pressure balancing gravity, has been predicted in many models with a QCD axion or axion-like particle. It can be formed in the early universe and composes a sizable fraction of dark matter. One could detect the transient radio signals when it passes by a magnetar with the axion particle converted into photon in the magnetic field. Usin…
▽ More
Diluted axion star, a self-gravitating object with the quantum pressure balancing gravity, has been predicted in many models with a QCD axion or axion-like particle. It can be formed in the early universe and composes a sizable fraction of dark matter. One could detect the transient radio signals when it passes by a magnetar with the axion particle converted into photon in the magnetic field. Using both numerical and semi-analytic approaches, we simulate the axion star's dynamic evolution and estimate the fraction of axion particles that can have a resonance conversion during such a collision event. We have found that both self-gravity and quantum pressure are not important after the diluted axion star enters the Roche radius. A free-fall approximate can capture individual particle trajectories very well. With some optimistic cosmological and astrophysical assumptions, the QCD axion parameter space can be probed from detecting such a collision event by radio telescopes.
△ Less
Submitted 2 September, 2021;
originally announced September 2021.
-
The Effects of Stellar Population and Gas Covering Fraction on the Emergent Lyman Alpha Emission of High-Redshift Galaxies
Authors:
Naveen A. Reddy,
Michael W. Topping,
Alice E. Shapley,
Charles C. Steidel,
Ryan L. Sanders,
Xinnan Du,
Alison L. Coil,
Bahram Mobasher,
Sedona H. Price,
Irene Shivaei
Abstract:
We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85<z<3.49 star-forming galaxies to deduce key properties of the massive stars, ionized ISM, and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Ly-alpha (Lya) photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra…
▽ More
We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85<z<3.49 star-forming galaxies to deduce key properties of the massive stars, ionized ISM, and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Ly-alpha (Lya) photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Ly-beta through CIII] 1907, 1909; and [OII], [NeIII], H-beta, [OIII], H-alpha, [NII], and [SII]. Spectral and photoionization modeling indicate that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Lya equivalent width, W(Lya), and only a marginal anti-correlation between age and W(Lya). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W(Lya). Rather, the covering fraction of optically-thick HI appears to be the principal factor modulating the escape of Lya, with most of the Lya photons in down-the-barrel observations of galaxies escaping through low-column-density or ionized channels in the ISM. Our analysis shows that a high star-formation-rate surface density, Sigma_SFR, particularly when coupled with a low galaxy potential (i.e., low stellar mass), can aid in reducing the covering fraction and ease the escape of Lya photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.
△ Less
Submitted 20 November, 2021; v1 submitted 11 August, 2021;
originally announced August 2021.
-
Hot and Dense Matter Equation of State Probability Distributions for Astrophysical Simulations
Authors:
Xingfu Du,
Andrew W. Steiner,
Jeremy W. Holt
Abstract:
We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microsco…
▽ More
We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. We compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.
△ Less
Submitted 14 July, 2021;
originally announced July 2021.
-
Light Bridges Can Suppress the Formation of Coronal Loops
Authors:
Yuhu Miao,
Libo Fu,
Xian Du,
Ding Yuan,
Chaowei Jiang,
Jiangtao Su,
Mingyu Zhao,
Sergey Anfinogentov
Abstract:
A light bridge is a magnetic intrusion into a sunspot, it interacts with the main magnetic field and excites a variety of dynamical processes. In the letter, we studied magnetic connectivity between a light bridge and coronal loops rooted at the sunspot. We used the data of the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory (SDO) to study the features of sunspots with light br…
▽ More
A light bridge is a magnetic intrusion into a sunspot, it interacts with the main magnetic field and excites a variety of dynamical processes. In the letter, we studied magnetic connectivity between a light bridge and coronal loops rooted at the sunspot. We used the data of the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory (SDO) to study the features of sunspots with light bridges. It is found that if a light bridge anchors at the umbra-penumbra boundary, the coronal loops could not be formed around the anchoring point. If the a light bridge become detached from the penumbra, the coronal loop starts to form again. The vector magnetogram provided by the Helioseismic Magnetic Imager onboard SDO shows that the anchoring region of a light bridge usually have an accompanying opposite minor-polarities. We conjugate that the magnetic field line could connect to these opposite polarities and form short-range magnetic loops, and therefore, coronal loops that extend to long-range could not be formed. A model of light bridge is proposed to explain the magnetic connectivity between a light bridge and the coronal loops. This model could explain many physical processes associated with light bridges.
△ Less
Submitted 24 June, 2021;
originally announced June 2021.
-
Core Formation in High-z Massive Haloes: Heating by Post Compaction Satellites and Response to AGN Outflows
Authors:
Avishai Dekel,
Jonathan Freundlich,
Fangzhou Jiang,
Sharon Lapiner,
Andreas Burkert,
Daniel Ceverino,
Xiaolong Du,
Reinhard Genzel,
Joel Primack
Abstract:
Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass $\geq\! 10^{12}M_\odot$ at $z\!\sim\! 2$. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-ma…
▽ More
Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass $\geq\! 10^{12}M_\odot$ at $z\!\sim\! 2$. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing AGN-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical-friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses $>\!10^{11.3}M_\odot$, that undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion onto haloes $\geq\!10^{12}M_\odot$, these satellites heat-up the cusps in half a virial time at $z\!\sim\! 2$. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGN indeed become effective in haloes $\geq\!10^{12}M_\odot$, where the black-hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot CGM. For simulations to reproduce the dynamical-friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black-hole accretion and clumpy response to AGN.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
The MOSDEF-LRIS Survey: Probing ISM/CGM Structure of Star-Forming Galaxies at z~2 Using Rest-UV Spectroscopy
Authors:
Xinnan Du,
Alice E. Shapley,
Michael W. Topping,
Naveen A. Reddy,
Ryan L. Sanders,
Alison L. Coil,
Mariska Kriek,
Bahram Mobasher,
Brian Siana
Abstract:
The complex structure of gas, metals, and dust in the interstellar and circumgalactic medium (ISM and CGM, respectively) in star-forming galaxies can be probed by Ly$α$ emission and absorption, low-ionization interstellar (LIS) metal absorption, and dust reddening E(B-V). We present a statistical analysis of the mutual correlations among Ly$α$ equivalent width (EW$_{Lyα}$), LIS equivalent width (E…
▽ More
The complex structure of gas, metals, and dust in the interstellar and circumgalactic medium (ISM and CGM, respectively) in star-forming galaxies can be probed by Ly$α$ emission and absorption, low-ionization interstellar (LIS) metal absorption, and dust reddening E(B-V). We present a statistical analysis of the mutual correlations among Ly$α$ equivalent width (EW$_{Lyα}$), LIS equivalent width (EW$_{LIS}$), and E(B-V) in a sample of 157 star-forming galaxies at $z\sim2.3$. With measurements obtained from individual, deep rest-UV spectra and spectral-energy distribution (SED) modeling, we find that the tightest correlation exists between EW$_{LIS}$ and E(B-V), although correlations among all three parameters are statistically significant. These results signal a direct connection between dust and metal-enriched HI gas, and that they are likely co-spatial. By comparing our results with the predictions of different ISM/CGM models, we favor a dusty ISM/CGM model where dust resides in HI gas clumps and Ly$α$ photons escape through the low HI covering fraction/column density intra-clump medium. Finally, we investigate the factors that potentially contribute to the intrinsic scatter in the correlations studied in this work, including metallicity, outflow kinematics, Ly$α$ production efficiency, and slit loss. Specifically, we find evidence that scatter in the relationship between EW$_{Lyα}$ and E(B-V) reflects the variation in metal-to-HI covering fraction ratio as a function of metallicity, and the effects of outflows on the porosity of the ISM/CGM. Future simulations incorporating star-formation feedback and the radiative transfer of Ly$α$ photons will provide key constraints on the spatial distributions of neutral hydrogen gas and dust in the ISM/CGM structure.
△ Less
Submitted 7 July, 2021; v1 submitted 29 March, 2021;
originally announced March 2021.
-
Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies
Authors:
Ethan O. Nadler,
Simon Birrer,
Daniel Gilman,
Risa H. Wechsler,
Xiaolong Du,
Andrew Benson,
Anna M. Nierenberg,
Tommaso Treu
Abstract:
Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and…
▽ More
Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model $\mathrm{\texttt{Galacticus}}$ to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.
△ Less
Submitted 9 August, 2021; v1 submitted 19 January, 2021;
originally announced January 2021.
-
New insights into the formation and growth of boson stars in dark matter halos
Authors:
Jiajun Chen,
Xiaolong Du,
Erik W. Lentz,
David J. E. Marsh,
Jens C. Niemeyer
Abstract:
This work studies the formation and growth of boson stars and their surrounding miniclusters by gravitational condensation using non-linear dynamical numerical methods. Fully dynamical attractive and repulsive self-interactions are also considered for the first time. In the case of pure gravity, we numerically prove that the growth of boson stars inside halos slows down and saturates as has been p…
▽ More
This work studies the formation and growth of boson stars and their surrounding miniclusters by gravitational condensation using non-linear dynamical numerical methods. Fully dynamical attractive and repulsive self-interactions are also considered for the first time. In the case of pure gravity, we numerically prove that the growth of boson stars inside halos slows down and saturates as has been previously conjectured, and detail its conditions. Self-interactions are included using the Gross-Pitaevskii-Poisson equations. We find that in the case of strong attractive self-interactions the boson stars can become unstable and collapse, in agreement with previous stationary computations. At even stronger coupling, the condensate fragments. Repulsive self-interactions, as expected, promote boson star formation, and lead to solutions with larger radii.
△ Less
Submitted 4 February, 2021; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Ensemble of meta-heuristic and exact algorithm based on the divide and conquer framework for multi-satellite observation scheduling
Authors:
Guohua Wu,
Qizhang Luo,
Xiao Du,
Xinwei Wang,
Yinguo Chen,
Ponnuthurai Nagaratnam Suganthan
Abstract:
Satellite observation scheduling plays a significant role in improving the efficiency of Earth observation systems. To solve the large-scale multi-satellite observation scheduling problem, this paper proposes an ensemble of meta-heuristic and exact algorithm based on a divide-and-conquer framework (EHE-DCF), including a task allocation phase and a task scheduling phase. In the task allocation phas…
▽ More
Satellite observation scheduling plays a significant role in improving the efficiency of Earth observation systems. To solve the large-scale multi-satellite observation scheduling problem, this paper proposes an ensemble of meta-heuristic and exact algorithm based on a divide-and-conquer framework (EHE-DCF), including a task allocation phase and a task scheduling phase. In the task allocation phase, each task is allocated to a proper orbit based on a meta-heuristic incorporated with a probabilistic selection and a tabu mechanism derived from ant colony optimization and tabu search respectively. In the task scheduling phase, we construct a task scheduling model for every single orbit, and use an exact method (i.e., branch and bound, B&B) to solve this model. The task allocation and task scheduling phases are performed iteratively to obtain a promising solution. To validate the performance of EHE-DCF, we compare it with B&B, three divide-and-conquer based meta-heuristics, and a state-of-the-art meta-heuristic. Experimental results show that EHE-DCF can obtain higher scheduling profits and complete more tasks compared with existing algorithms. EHE-DCF is especially efficient for large-scale satellite observation scheduling problems.
△ Less
Submitted 6 September, 2021; v1 submitted 7 July, 2020;
originally announced July 2020.
-
SatGen: a semi-analytical satellite galaxy generator -- I. The model and its application to Local-Group satellite statistics
Authors:
Fangzhou Jiang,
Avishai Dekel,
Jonathan Freundlich,
Frank C. van den Bosch,
Sheridan B. Green,
Philip F. Hopkins,
Andrew Benson,
Xiaolong Du
Abstract:
We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates c…
▽ More
We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates cosmological zoom-in hydro-simulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydro-simulations and physical modeling. The subhalo profile and the stellar mass and size of a satellite evolves depending on its tidal mass loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. Our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25% level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW disc depletes satellites by 20% and has a subtle effect of diversifying the internal structure of satellites, important for alleviating certain small-scale problems.
△ Less
Submitted 12 May, 2020;
originally announced May 2020.
-
A new calibration method of sub-halo orbital evolution for semi-analytic models
Authors:
Shengqi Yang,
Xiaolong Du,
Andrew J. Benson,
Anthony R. Pullen,
Annika H. G. Peter
Abstract:
Understanding the non-linear dynamics of satellite halos (a.k.a. "sub-halos") is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a mo…
▽ More
Understanding the non-linear dynamics of satellite halos (a.k.a. "sub-halos") is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a more efficient way to describe the key physical processes such as dynamical friction, tidal mass loss, and tidal heating, with only a few free parameters. In this work, we present a fast Monte Carlo Markov Chain fitting approach to explore the parameter space of such a sub-halo non-linear evolution model. We use the dynamical models described in an earlier work and calibrate the models to two sets of high-resolution cold dark matter N-body simulations, ELVIS and Caterpillar. Compared to previous calibrations that used manual parameter tuning, our approach provides a more robust way to determine the best-fit parameters and their posterior probabilities. We find that jointly fitting for the sub-halo mass and maximum velocity functions can break the degeneracy between tidal stripping and tidal heating parameters, as well as providing better constraints on the strength of dynamical friction. We show that our semi-analytic simulation can accurately reproduce N-body simulations statistics, and that the calibration results for the two sets of N-body simulations agree at 95% confidence level. Dynamical models calibrated in this work will be important for future dark matter substructure studies.
△ Less
Submitted 23 September, 2020; v1 submitted 23 March, 2020;
originally announced March 2020.
-
The DRAO Synthesis Telescope
Authors:
T. L. Landecker,
L. Belostotski,
J. C. Brown,
B. R. Carlson,
X. Du,
T. Foster,
A. S. Hill,
T. Johnson,
R. Kothes,
M. Rupen,
S. Safi-Harb,
K. Vanderlinde,
B. G. Veidt,
J. L. West
Abstract:
The DRAO Synthesis Telescope (ST) is a forefront telescope for imaging large-scale neutral hydrogen and polarized radio continuum emission at arcminute resolution. Equipped for observations at 1420 and 408 MHz, the ST completed the Canadian Galactic Plane Survey, providing pioneering measurements of arcminute-scale structure in HI emission and self-absorption and of the diffuse polarized emission,…
▽ More
The DRAO Synthesis Telescope (ST) is a forefront telescope for imaging large-scale neutral hydrogen and polarized radio continuum emission at arcminute resolution. Equipped for observations at 1420 and 408 MHz, the ST completed the Canadian Galactic Plane Survey, providing pioneering measurements of arcminute-scale structure in HI emission and self-absorption and of the diffuse polarized emission, using a fine grid of Rotation Measures to chart the large-scale Galactic magnetic field, and advancing the knowledge of the Galactic rotation curve. In this paper we describe a plan for renewal of the Synthesis Telescope that will create a forefront scientific instrument, a testbed for new radio astronomy technologies, and a training ground for the next generation of Canadian radio astronomers and radio telescope engineers. The renewed telescope will operate across the entire range 400 to 1800 MHz. Collaborations between DRAO and university partners have already demonstrated a novel feed antenna to cover this range, low-noise amplifiers, and a new GPU-based correlator of bandwidth 400 MHz. The renewed ST will provide excellent sensitivity to extended HI, covering the Galactic disk and halo, spectro-polarimetry with unprecedented resolution in angle and in Faraday depth, the ability to search for OH masers in all four 18-cm lines simultaneously, and sensitive recombination-line observations stacked over as many as forty transitions. As a testbed the renewed ST will evaluate low-cost digital clocking and sampling techniques of wide significance for the ngVLA, SKA, and other future telescopes, and a prototype of the digital correlator developed at DRAO for SKA-mid.
△ Less
Submitted 26 October, 2019;
originally announced October 2019.
-
Searching for z > 6.5 Analogs Near the Peak of Cosmic Star Formation
Authors:
Xinnan Du,
Alice E. Shapley,
Mengtao Tang,
Daniel P. Stark,
Crystal L. Martin,
Bahram Mobasher,
Michael W. Topping,
Jacopo Chevallard
Abstract:
Strong [OIII]$λλ$4959,5007+H$β$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$α$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the…
▽ More
Strong [OIII]$λλ$4959,5007+H$β$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$α$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$λλ$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$α$ and CIII]$λλ$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$α$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$α$ emission starts to emerge. Using stacked spectra, we find that Ly$α$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$λλ$1548,1550, HeII$λ$1640, and OIII]$λλ$1661,1665 emission in addition to intense Ly$α$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts.
△ Less
Submitted 4 January, 2020; v1 submitted 25 October, 2019;
originally announced October 2019.
-
The Redshift Evolution of Rest-UV Spectroscopic Properties to z~5
Authors:
Anthony J. Pahl,
Alice Shapley,
Andreas L. Faisst,
Peter L. Capak,
Xinnan Du,
Naveen A. Reddy,
Peter Laursen,
Michael W. Topping
Abstract:
We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$α$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$α$ equivalent width, p…
▽ More
We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$α$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$α$ equivalent width, performing corrections to Ly$α$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$α$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$α$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$α$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$α$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$α$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$α$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$α$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$α$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$α$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($ξ_{\rm ion}$).
△ Less
Submitted 7 February, 2020; v1 submitted 9 October, 2019;
originally announced October 2019.
-
Constraints on the mass-concentration relation of cold dark matter halos with 11 strong gravitational lenses
Authors:
Daniel Gilman,
Xiaolong Du,
Andrew Benson,
Simon Birrer,
Anna Nierenberg,
Tommaso Treu
Abstract:
The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the i…
▽ More
The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance.
△ Less
Submitted 9 September, 2019; v1 submitted 5 September, 2019;
originally announced September 2019.
-
Warm dark matter chills out: constraints on the halo mass function and the free-streaming length of dark matter with 8 quadruple-image strong gravitational lenses
Authors:
Daniel Gilman,
Simon Birrer,
Anna Nierenberg,
Tommaso Treu,
Xiaolong Du,
Andrew Benson
Abstract:
The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhal…
▽ More
The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhalos and halos along the line of sight, and account for warm dark matter (WDM) free-streaming effects on the mass function and mass-concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated halos, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parameterizing dark matter free-streaming through the half-mode mass $m_{\rm{hm}}$, we constrain the thermal relic particle mass $m_{\rm{DM}}$ corresponding to $m_{\rm{hm}}$. At $95 \%$ CI: $m_{\rm{hm}} < 10^{7.8} M_{\odot}$ ($m_{\rm{DM}} > 5.2 \ \rm{keV}$). We disfavor $m_{\rm{DM}} = 4.0 \rm{keV}$ and $ m_{\rm{DM}} = 3.0 \rm{keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between $10^6 - 10^{9} M_{\odot}$ near lensed images. At $68 \%$ CI, we infer $2.0 - 6.1 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm{sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \%$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to $\bar{f}_{\rm{sub}} > 0.005$. These results agree with the predictions of cold dark matter.
△ Less
Submitted 11 November, 2019; v1 submitted 19 August, 2019;
originally announced August 2019.
-
From the microscopic to the macroscopic world: from nucleons to neutron stars
Authors:
S. Gandolfi,
J. Lippuner,
A. W. Steiner,
I. Tews,
X. Du,
M. Al-Mamun
Abstract:
Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactio…
▽ More
Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactions among nucleons. We then review selected properties of neutron stars that can be directly compared with present and future observations.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Hyperons and quarks in proto-neutron stars
Authors:
J. Roark,
X. Du,
C. Constantinou,
V. Dexheimer,
A. W. Steiner,
J. R. Stone
Abstract:
In this work, we study matter in the cores of proto-neutron stars, focusing on the impact of their composition on the stellar structure. We begin by examining the effects of finite temperature (through a fixed entropy per baryon) and lepton fraction on purely nucleonic matter by making use of the DSH model . We then turn our attention to a relativistic mean-field model containing exotic degrees of…
▽ More
In this work, we study matter in the cores of proto-neutron stars, focusing on the impact of their composition on the stellar structure. We begin by examining the effects of finite temperature (through a fixed entropy per baryon) and lepton fraction on purely nucleonic matter by making use of the DSH model . We then turn our attention to a relativistic mean-field model containing exotic degrees of freedom, the Chiral Mean Field (CMF) model, again, under the conditions of finite temperature and trapped neutrinos. In the latter, since both hyperons and quarks are found in the cores of large-mass stars, their interplay and the possibility of mixtures of phases is taken into account and analyzed. Finally, we discuss how stellar rotation can affect our results.
△ Less
Submitted 13 May, 2019; v1 submitted 19 December, 2018;
originally announced December 2018.
-
The redshift evolution of rest-UV spectroscopic properties in Lyman Break Galaxies at z ~ 2-4
Authors:
Xinnan Du,
Alice E. Shapley,
Naveen A. Reddy,
Tucker Jones,
Daniel P. Stark,
Charles C. Steidel,
Allison L. Strom,
Gwen C. Rudie,
Dawn K. Erb,
Richard S. Ellis,
Max Pettini
Abstract:
We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$α$ equivalent width (EW), and various gal…
▽ More
We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$α$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$α$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$α$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$α$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$α$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$α$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$α$ strength with redshift, and the variation in Ly$α$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$α$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$α$ photons at the highest Ly$α$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6.
△ Less
Submitted 20 April, 2018; v1 submitted 15 March, 2018;
originally announced March 2018.
-
Hot and Dense Homogeneous Nucleonic Matter Constrained by Observations, Experiment, and Theory
Authors:
Xingfu Du,
Andrew W. Steiner,
Jeremy W. Holt
Abstract:
We construct a new class of phenomenological equations of state for homogeneous matter for use in simulations of hot and dense matter in local thermodynamic equilibrium. We construct a functional form which respects experimental, observational and theoretical constraints on the nature of matter in various density and temperature regimes. Our equation of state matches (i) the virial coefficients ex…
▽ More
We construct a new class of phenomenological equations of state for homogeneous matter for use in simulations of hot and dense matter in local thermodynamic equilibrium. We construct a functional form which respects experimental, observational and theoretical constraints on the nature of matter in various density and temperature regimes. Our equation of state matches (i) the virial coefficients expected from nucleon-nucleon scattering phase shifts, (ii) experimental measurements of nuclear masses and charge radii, (iii) observations of neutron star radii, (iv) theory results on the equation of state of neutron matter near the saturation density, and (v) theory results on the evolution of the EOS at finite temperatures near the saturation density. Our analytical model allows one to compute the variation in the thermodynamic quantities based on the uncertainties in the nature of the nucleon-nucleon interaction. Finally, we perform a correction to ensure the equation of state is causal at all densities, temperatures, and electron fractions.
△ Less
Submitted 26 February, 2018;
originally announced February 2018.
-
Molecular Gas Toward the Gemini OB1 Molecular Cloud Complex II: CO Outflow Candidates with Possible WISE Associations
Authors:
Yingjie Li,
Fa-Cheng Li,
Ye Xu,
Chen Wang,
Xin-Yu Du,
Wenjin Yang,
Ji Yang
Abstract:
We present a large scale survey of CO outflows in the Gem OB1 molecular cloud complex and its surroundings using the Purple Mountain Observatory Delingha 13.7 m telescope. A total of 198 outflow candidates were identified over a large area ($\sim$ 58.5 square degrees), of which 193 are newly detected. Approximately 68% (134/198) are associated with the Gem OB1 molecular cloud complex, including cl…
▽ More
We present a large scale survey of CO outflows in the Gem OB1 molecular cloud complex and its surroundings using the Purple Mountain Observatory Delingha 13.7 m telescope. A total of 198 outflow candidates were identified over a large area ($\sim$ 58.5 square degrees), of which 193 are newly detected. Approximately 68% (134/198) are associated with the Gem OB1 molecular cloud complex, including clouds GGMC 1, GGMC 2, BFS 52, GGMC 3 and GGMC 4. Other regions studied are: Local Arm (Local Lynds, West Front), Swallow, Horn, and Remote cloud. Outflow candidates in GGMC 1, BFS 52, and Swallow are mainly located at ring-like or filamentary structures. To avoid excessive uncertainty in distant regions ($\gtrsim$ 3.8 kpc), we only estimated the physical parameters for clouds in the Gem OB1 molecular cloud complex and in the Local arm. In those clouds, the total kinetic energy and the energy injection rate of the identified outflow candidates are $\lesssim$ 1% and $\lesssim$ 3% of the turbulent energy and the turbulent dissipation rate of each cloud, indicating that the identified outflow candidates cannot provide enough energy to balance turbulence of their host cloud at the scale of the entire cloud (several pc to dozens of pc). The gravitational binding energy of each cloud is $\gtrsim$ 135 times the total kinetic energy of the identified outflow candidates within the corresponding cloud, indicating that the identified outflow candidates cannot cause major disruptions to the integrity of their host cloud at the scale of the entire cloud.
△ Less
Submitted 31 January, 2018;
originally announced January 2018.