Showing 1–50 of 124 results for author: Bird, S

The Stellar "Snake" -- II: The Mass Function

Authors: Xiang-Ming Yang, Sarah A. Bird, Jiadong Li, Hai-Jun Tian, Dan Qiu, Jia-Peng Li, Cheng-Yuan Li, Gao-Chao Liu, Peng Zhang, Ju-Yong Zhang, Zhi-Ping Chen

Abstract: We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar ``Snake'' in the latest {\it Gaia} Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical para… ▽ More We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar ``Snake'' in the latest {\it Gaia} Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical parameters of the clusters within the Stellar Snake. In light of the high completeness of the member stars in the cluster regions, we develop a simulated model color-magnitude diagram-based inference method to derive the mass function, binary fraction, and mass-ratio distribution of the clusters in the Stellar Snake. Notably, despite their similar ages and metallicity, we discover systematic variations in the MFs along the elongation direction of the Snake in the mass range of 0.5 to 2.0 M$_\odot$. The ``head'' of the Snake conforms to a canonical initial mass function with a power-law slope of $α\sim-2.3$. Extending towards the ``tail,'' the MF becomes more top-light, indicating a deficiency of massive stars within these clusters. This result provides evidence for the delayed formation of massive stars in the clusters. Such clues give support to the hypothesis that the Stellar Snake constitutes as a hierarchically primordial structure. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 19 pages, 15 figures, Accepted for publication on MNRAS

The evolution of galaxy morphology from redshift z=6 to 3: Mock JWST observations of galaxies in the ASTRID simulation

Authors: Patrick LaChance, Rupert Croft, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Simeon Bird

Abstract: We present mock JWST observations for more than 215,000 different galaxies from the Astrid simulation with $3 \leq z \leq 6$. The mock observations are made using the BPASS stellar SED model, and a simple dust model. They are then viewed through NIRCam filters, convolved with a PSF, have noise added, and are drizzled together to emulate the Cosmic Evolution Early Release Science (CEERS) survey. We… ▽ More We present mock JWST observations for more than 215,000 different galaxies from the Astrid simulation with $3 \leq z \leq 6$. The mock observations are made using the BPASS stellar SED model, and a simple dust model. They are then viewed through NIRCam filters, convolved with a PSF, have noise added, and are drizzled together to emulate the Cosmic Evolution Early Release Science (CEERS) survey. We analyse this dataset by computing a number of morphological measures and find our catalog to have comparable statistics to similar mock catalogs, and the first release of CEERS data. We find that most of the Sersic indices of galaxies in our redshift range are lower than observed, with most having n less than one. Additionally, we observe the sizes of galaxies of all masses to increase from redshift z=6 to redshift z=3 consistent with other results. The number of galaxies in our catalog allows us to examine how relationships like the mass-size relation evolve with redshift, and compare the accuracy of a variety of traditional galaxy classification techniques (Sersic fit, Asymmetry-Concentration, and Gini-$M_{20}$) within our redshift range. We find the mass-size relation to be nearly flat at redshift z=6, and consistently increases as redshift decreases, and find the galaxy classification methods have minimal correlation with each other in our redshift range. We also investigate the impact that different stages of our imaging pipeline have on these morphological measures to determine how robust mock catalogs are to different choices at each step. Finally, we test the addition of incorporating light from AGNs into our pipeline and find that while the population of galaxies that have significant AGN luminosity is low, those galaxies do tend to have higher Sersic indices once the AGN luminosity is added, rectifying some of the systematic bias towards lower Sersic indices present in our dataset. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 17 pages, 14 figures

North-South asymmetries in the Galactic thin disk associated with the vertical phase spiral as seen using LAMOST-Gaia stars

Authors: Jun Lin, Rui Guo, Sarah A. Bird, Haijun Tian, Chao Liu, Chris Flynn, Gaochao Liu, Sheng Cui

Abstract: We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($σ_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and… ▽ More We select 1,052,469 (754,635) thin disk stars from {\it Gaia} eDR3 and LAMOST DR7 in the range of Galactocentric radius $R$ (guiding center radius $R_\mathrm{g}$) from 8 to 11\,kpc to investigate the asymmetries between the North and South of the disk midplane. More specifically we analyze the vertical velocity dispersion profiles ($σ_{v_{z}}(z$)) in different bins of $R$ ($R_\mathrm{g}$) and $[\mathrm{Fe/H}]$. We find troughs in the profiles of $σ_{v_{z}}(z)$ located in both the North ($z \sim 0.7$\,kpc) and South ($z \sim -0.5$\,kpc) of the disk at all radial and chemical bins studied. The difference between the Northern and Southern vertical velocity dispersion profiles ($Δσ_{v_{z}}(|z|)$) shows a shift between curves of different $R$ and $R_\mathrm{g}$. A similar shift exists in these NS asymmetry profiles further divided into different $[\mathrm{Fe/H}]$ ranges. The sample binned with $R_\mathrm{g}$ more clearly displays the features in the velocity dispersion profiles. The shift in the peaks of the $Δσ_{v_{z}}$ profiles and the variation in the phase spiral shape binned by metallicity indicate the variation of the vertical potential profiles and the radial metallicity gradient. The wave-like signal in NS asymmetry of $σ_{v_{z}}(z)$ largely originates from phase spiral; while the NS asymmetry profiles of [Fe/H] only display a weak wave-like feature near solar radius. We perform a test particle simulation to qualitatively reproduce the observed results. A quantitative explanation of the NS asymmetry in the metallicity profile needs careful consideration of the spiral shape and the perturbation model, and we leave this for future work. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

z~2 dual AGN host galaxies are disky: stellar kinematics in the ASTRID Simulation

Authors: Ekaterine Dadiani, Tiziana Di Matteo, Nianyi Chen, Patrick Lachance, Yue Shen, Yu-Ching Chen, Rupert Croft, Yueying Ni, Simeon Bird

Abstract: We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δr<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'b… ▽ More We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δr<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'bulge' (bulge and halo) and define disk-dominated galaxies by the disk-to-total $D/T\geq0.5$. In ASTRID, $60.9\pm2.1\%$ of dual AGN hosts (independent of separation) are disk-dominated, with the $D/T$ distribution peaking at $\sim0.7$. Notably, hosts of BH pairs have similar morphologies (most either both disk or bulge-dominated). In dual-AGN hosts, the $D/T$ increases from $\sim17\% $ at $M_{\rm *}\sim 10^{9} M_{\odot}$ to $ 64\% $ for $M_{\rm *} \sim 10^{11.5} M_{\odot}$, and the pseudo-bulge is the dominant component of the disk fraction at the high mass end. Moreover, dual AGN hosts exhibit a higher fraction of disk/large pseudo-bulge than single-AGN hosts. The Disk-to-Total ratio is approximately constant with BH mass or AGN luminosity. We also create mock images of dual AGN host galaxies, employing morphological fitting software Statmorph to calculate morphological parameters and compare them with our kinematic decomposition results. Around $83.3\pm2.4\%$ of galaxies display disk-like profiles, of which $\sim60.7\pm2.2\%$ are kinematically confirmed as disks. Seŕsic indices and half-mass radii of dual AGN host galaxies align with observational measurements from HST at $z\sim2$. Around $34\%$ are identified as mergers from the $\text{Gini}-M_{20}$ relation. We find two dual AGN hosted by galaxies that exhibit disk-like seŕsic index $n_{12}<1$ and $(D/T)_{12}>0.5$, which are in remarkable agreement with properties of recently discovered dual quasars in disk galaxies at $z\sim 2$. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 15 pages, 12 figures, submitted to the Open Journal of Astrophysics

Astrometric Microlensing by Primordial Black Holes with The Roman Space Telescope

Authors: James Fardeen, Peter McGill, Scott E. Perkins, William A. Dawson, Natasha S. Abrams, Jessica R. Lu, Ming-Feng Ho, Simeon Bird

Abstract: Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing… ▽ More Primordial Black Holes (PBHs) could explain some fraction of dark matter and shed light on many areas of early-universe physics. Despite over half a century of research interest, a PBH population has so far eluded detection. The most competitive constraints on the fraction of dark matter comprised of PBHs ($f_{\rm DM}$) in the $(10^{-9}-10)M_{\odot}$ mass-ranges come from photometric microlensing and bound $f_{\rm DM}\lesssim10^{-2}-10^{-1}$. With the advent of the Roman Space Telescope with its sub-milliarcsecond (mas) astrometric capabilities and its planned Galactic Bulge Time Domain Survey (GBTDS), detecting astrometric microlensing signatures will become routine. Compared with photometric microlensing, astrometric microlensing signals are sensitive to different lens masses-distance configurations and contains different information, making it a complimentary lensing probe. At sub-mas astrometric precision, astrometric microlensing signals are typically detectable at larger lens-source separations than photometric signals, suggesting a microlensing detection channel of pure astrometric events. We use a Galactic simulation to predict the number of detectable microlensing events during the GBTDS via this pure astrometric microlensing channel. Assuming an absolute astrometric precision floor for bright stars of 0.1 mas for the GBTDS, we find that the number of detectable events peaks at $\approx 10^{3} f_{\rm DM}$ for a population of $ 1 M_{\odot}$ PBHs and tapers to $\approx 10f_{\rm DM}$ and $\approx 100f_{\rm DM}$ at $10^{-4}M_{\odot}$ and $10^{3}M_{\odot}$, respectively. Accounting for the distinguishability of PBHs from Stellar lenses, we conclude the GBTDS will be sensitive to a PBH population at $f_{\rm DM}$ down to $\approx10^{-1}-10^{-3}$ for $(10^{-1}-10^{2})M_{\odot}$ likely yielding novel PBH constraints. △ Less

Submitted 8 March, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 21 pages, 11 figures, accepted to AAS Journals

MAGICS I. The First Few Orbits Encode the Fate of Seed Massive Black Hole Pairs

Authors: Nianyi Chen, Diptajyoti Mukherjee, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft

Abstract: The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered fro… ▽ More The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered from large-volume cosmological simulation ASTRID. The simulations include subgrid star formation, supernovae feedback, and MBH accretion/feedback. In this first suite, we extract fifteen representative galaxy mergers with seed MBHs to examine their dynamics at an improved mass and spatial resolution (by $\sim2000$ and $\sim20$) and follow MBH orbits down to $\sim10\,\text{pc}$. We find that the seed MBH energy loss and orbital decay are largely governed by global torques induced by the galaxy merger process on scales resolvable by cosmological simulations. Specifically, pairs sink quickly if their orbits shrink rapidly below $1\,\text{kpc}$ during the first $\sim200\,\text{Myr}$ of pairing due to effective energy loss in major galaxy mergers, whereas MBHs gaining energy in minor galaxy mergers with head-on collisions are likely to stall. High initial eccentricities ($e_\text{init}>0.5$) and high stellar densities at kpc scales ($ρ_\text{star}>0.05\,M_\odot/\text{pc}^3$) also lead to most efficient decays. $\sim50\%$ high-redshift seed MBH pairs experience consecutive galaxy mergers and are more likely to stall at $\sim1\,\text{kpc}$. For a subset of systems, we carry out N-Body re-simulations until binary formation and find that some stalled systems merge at high-z when embedded in sufficient nuclear star clusters. △ Less

Submitted 25 April, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 21 pages, 20 Figures. Submitted to The Open Journal of Astrophysics. Comments welcome!

SDSS-IV from 2014 to 2016: A Detailed Demographic Comparison over Three Years

Authors: Amy M. Jones, Rachael L. Beaton, Brian A. Cherinka, Karen L. Masters, Sara Lucatello, Aleksandar M. Diamond-Stanic, Sarah A. Bird, Michael R. Blanton, Katia Cunha, Emily E. Farr, Diane Feuillet, Peter M. Frinchaboy, Alex Hagen, Karen Kinemuchi, Britt Lundgren, Mariarosa L. Marinelli, Adam D. Myers, Alexandre Roman-Lopes, Ashley J. Ross, Jose R. Sanchez-Gallego, Sarah J. Schmidt, Jennifer Sobeck, Keivan G. Stassun, Jamie Tayar, Mariana Vargas-Magana , et al. (2 additional authors not shown)

Abstract: The Sloan Digital Sky Survey (SDSS) is one of the largest international astronomy organizations. We present demographic data based on surveys of its members from 2014, 2015 and 2016, during the fourth phase of SDSS (SDSS-IV). We find about half of SDSS-IV collaboration members were based in North America, a quarter in Europe, and the remainder in Asia and Central and South America. Overall, 26-36%… ▽ More The Sloan Digital Sky Survey (SDSS) is one of the largest international astronomy organizations. We present demographic data based on surveys of its members from 2014, 2015 and 2016, during the fourth phase of SDSS (SDSS-IV). We find about half of SDSS-IV collaboration members were based in North America, a quarter in Europe, and the remainder in Asia and Central and South America. Overall, 26-36% are women (from 2014 to 2016), up to 2% report non-binary genders. 11-14% report that they are racial or ethnic minorities where they live. The fraction of women drops with seniority, and is also lower among collaboration leadership. Men in SDSS-IV were more likely to report being in a leadership role, and for the role to be funded and formally recognized. SDSS-IV collaboration members are twice as likely to have a parent with a college degree, than the general population, and are ten times more likely to have a parent with a PhD. This trend is slightly enhanced for female collaboration members. Despite this, the fraction of first generation college students (FGCS) is significant (31%). This fraction increased among collaboration members who are racial or ethnic minorities (40-50%), and decreased among women (15-25%). SDSS-IV implemented many inclusive policies and established a dedicated committee, the Committee on INclusiveness in SDSS (COINS). More than 60% of the collaboration agree that the collaboration is inclusive; however, collaboration leadership more strongly agree with this than the general membership. In this paper, we explain these results in full, including the history of inclusive efforts in SDSS-IV. We conclude with a list of suggested recommendations based on our findings, which can be used to improve equity and inclusion in large astronomical collaborations, which we argue is not only moral, but will also optimize their scientific output. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 30 pages, 9 figures, accepted in PASP

LATIS: Constraints on the Galaxy-halo Connection at $z \sim 2.5$ from Galaxy-galaxy and Galaxy-Ly$α$ Clustering

Authors: Andrew B. Newman, Mahdi Qezlou, Nima Chartab, Gwen C. Rudie, Guillermo A. Blanc, Simeon Bird, Andrew J. Benson, Daniel D. Kelson, Brian C. Lemaux

Abstract: The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the g… ▽ More The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Ly$α$ absorption measured in transverse sightlines; however, such studies are very few and have produced some unexpected or inconclusive results. We use $\sim$3000 spectra of $z\sim2.5$ galaxies from the Lyman-alpha Tomography IMACS Survey (LATIS) to measure the galaxy-galaxy and galaxy-Ly$α$ correlation functions in four bins of stellar mass spanning $10^{9.2} \lesssim M_* / M_{\odot} \lesssim 10^{10.5}$. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and are broadly in accord with standard SMHM relations. This consistency demonstrates that we are able to accurately measure and model Ly$α$ transmission fluctuations $δ_F$ in LATIS. We also show that the galaxy-Ly$α$ cross-correlation, a free byproduct of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy-galaxy clustering. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted to ApJ

Forecast Cosmological Constraints with the 1D Wavelet Scattering Transform and the Lyman-$α$ forest

Authors: Hurum Tohfa, Simeon Bird, Ming-Feng Ho, Mahdi Qezlou, Martin Fernandez

Abstract: We make forecasts for the constraining power of the 1D Wavelet Scattering Transform (WST) when used with a Lyman-$α$ forest cosmology survey. Using mock simulations and a Fisher matrix, we show that there is considerable cosmological information in the scattering transform coefficients not captured by the flux power spectrum. We estimate mock covariance matrices assuming uncorrelated Gaussian pixe… ▽ More We make forecasts for the constraining power of the 1D Wavelet Scattering Transform (WST) when used with a Lyman-$α$ forest cosmology survey. Using mock simulations and a Fisher matrix, we show that there is considerable cosmological information in the scattering transform coefficients not captured by the flux power spectrum. We estimate mock covariance matrices assuming uncorrelated Gaussian pixel noise for each quasar, at a level drawn from a simple lognormal model. The extra information comes from a smaller estimated covariance in the first-order wavelet power, and from second-order wavelet coefficients which probe non-Gaussian information in the forest. Forecast constraints on cosmological parameters from the WST are more than an order of magnitude tighter than for the power spectrum, shrinking a $4D$ parameter space by a factor of $10^6$. Should these improvements be realised with DESI, inflationary running would be constrained to test common inflationary models predicting $α_s = - 6\times 10^{-4}$ and neutrino mass constraints would be improved enough for a $5-σ$ detection of the minimal neutrino mass. △ Less

Submitted 3 June, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Journal ref: The 7 June 2024 issue of Physical Review Letters (Vol. 132, No. 23)

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

Authors: Scott Ellis Perkins, Peter McGill, William Dawson, Natasha S. Abrams, Casey Y. Lam, Ming-Feng Ho, Jessica R. Lu, Simeon Bird, Kerianne Pruett, Nathan Golovich, George Chapline

Abstract: From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However,… ▽ More From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens-class of any single event. We investigate this method's effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of $30 M_{\odot}$. On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ${\approx}25$\%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 31 pages, 18 figures, submitted to AAS

Modeling Quasar Proximity Zones in a Realistic Cosmological Environment with a Self-consistent Light Curve

Authors: Yihao Zhou, Huanqing Chen, Tiziana Di Matteo, Yueying Ni, Rupert A. C. Croft, Simeon Bird

Abstract: We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two or… ▽ More We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two orders of magnitude. Using this light curve to post-process the dense environment around the quasar, we find that the proximity zone size ($R_{p}$) ranges between $0.5-5$ pMpc. We show that the light curve variability causes a similar degree of scatter in $R_{p}$ as does the density fluctuation, both of which result in a standard deviation of $\sim 0.3$ pMpc). The $R_{p}$ traces the light curve fluctuations closely but with a time delay of $\sim 10^4\ \mathrm{yr}$, breaking the correspondence between the $R_{p}$ and the contemporaneous $M_{UV}$. This also indicates that we can only infer quasar activity within the past $\sim 10^4$ years instead of the integrated lifetime from $R_{p}$ in the later part of cosmic reionization. Compared with the variable light curve, a constant light curve underestimates the $R_{p}$ by 13% at the dim end ($M_{UV}\sim -23.5$), and overestimates the $R_{p}$ by 30% at the bright end ($M_{UV}\sim -26$). By calculating the $R_{p}$ generated by a number of quasars, we show that variable light curves predict a wider $R_{p}$ distribution than lightbulb models, and readily explain the extremely small $R_{p}$ values that have been observed. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 15 pages, 12 figures, comments welcome

Cosmological Constraints from the eBOSS Lyman-$α$ Forest using the PRIYA Simulations

Authors: M. A. Fernandez, Simeon Bird, Ming-Feng Ho

Abstract: We present new cosmological parameter constraints from the eBOSS Lyman-$α$ forest survey. We use a new theoretical model and likelihood based on the PRIYA simulation suite. PRIYA is the first suite to resolve the Lyman-$α$ forest in a ($120$~Mpc/h~)$^3$ volume, using a multi-fidelity emulation technique. We use PRIYA to predict Lyman-$α$ forest observables with $\lesssim 1\%$ interpolation error o… ▽ More We present new cosmological parameter constraints from the eBOSS Lyman-$α$ forest survey. We use a new theoretical model and likelihood based on the PRIYA simulation suite. PRIYA is the first suite to resolve the Lyman-$α$ forest in a ($120$~Mpc/h~)$^3$ volume, using a multi-fidelity emulation technique. We use PRIYA to predict Lyman-$α$ forest observables with $\lesssim 1\%$ interpolation error over an $11$ dimensional ($9$ simulated, $2$ in post-processing) parameter space. We identify an internal tension within the flux power spectrum data. Once the discrepant data is removed, we find the primeval scalar spectral index measured at a pivot scale of $k_0 = 0.78$ Mpc$^{-1}$ to be $n_P = 1.009^{+0.027}_{-0.018}$ at 68\% confidence. This measurement from the Lyman-$α$ forest flux power spectrum alone is in reasonable agreement with Planck, and in tension with earlier eBOSS analyses. The amplitude of matter fluctuations is $σ_8 = 0.733^{+0.026}_{-0.029}$ at 68\% confidence, in agreement with Dark Energy Survey weak lensing measurements and other small-scale structure probes and in tension with CMB measurements from Planck and ACT. The effective optical depth to Lyman-$α$ photons from our pipeline is in good agreement with earlier high resolution measurements. We find a linear power at $z=3$ and $k = 0.009$ s/km of $Δ_L^2 = 0.302^{+0.024}_{-0.027}$ with a slope $n_\mathrm{eff} = -2.264^{+0.026}_{-0.018}$. Our flux power spectrum only chains prefer a low level of heating during helium reionization. When we add IGM temperature data we find $n_P = 0.983\pm 0.020$ and $σ_8 = 0.703^{+0.023}_{-0.027}$. Our chains prefer an early and long helium reionization event, as suggested by measurements from the helium Lyman-$α$ forest. In the near future we will use our pipeline to infer cosmological parameters from the DESI Lyman-$α$ data. △ Less

Submitted 17 July, 2024; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 41 pages, 13 figures, published in JCAP https://iopscience.iop.org/article/10.1088/1475-7516/2024/07/029

Journal ref: JCAP07(2024)029

Probing the Galactic halo with RR Lyrae stars -- IV. On the Oosterhoff dichotomy of RR Lyrae stars

Authors: Shan Zhang, Gaochao Liu, Yang Huang, Zongfei Lv, Sarah Ann Bird, Bingqiu Chen, Huawei Zhang, Timothy C. Beers, Xinyi Li, Haijun Tian, Peng Zhang

Abstract: We use 3653 (2661 RRab, 992 RRc) RR Lyrae stars (RRLs) with 7D (3D position, 3D velocity, and metallicity) information selected from SDSS, LAMOST, and Gaia EDR3, and divide the sample into two Oosterhoff groups (Oo I and Oo II) according to their amplitude-period behaviour in the Bailey Diagram. We present a comparative study of these two groups based on chemistry, kinematics, and dynamics. We fin… ▽ More We use 3653 (2661 RRab, 992 RRc) RR Lyrae stars (RRLs) with 7D (3D position, 3D velocity, and metallicity) information selected from SDSS, LAMOST, and Gaia EDR3, and divide the sample into two Oosterhoff groups (Oo I and Oo II) according to their amplitude-period behaviour in the Bailey Diagram. We present a comparative study of these two groups based on chemistry, kinematics, and dynamics. We find that Oo I RRLs are relatively more metal rich, with predominately radially dominated orbits and large eccentricities, while Oo II RRLs are relatively more metal poor, and have mildly radially dominated orbits. The Oosterhoff dichotomy of the Milky Way's halo is more apparent for the inner-halo region than for the outer-halo region. Additionally, we also search for this phenomenon in the halos of the two largest satellite galaxies, the Large and Small Magellanic clouds (LMC, SMC), and compare over different bins in metallicity. We find that the Oosterhoff dichotomy is not immutable, and varies based on position in the Galaxy and from galaxy-to-galaxy. We conclude that the Oosterhoff dichotomy is the result of a combination of stellar and galactic evolution, and that it is much more complex than the dichotomy originally identified in Galactic globular clusters. △ Less

Submitted 12 September, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Exploring asymmetric substructures of the outer disk based on the conjugate angle of the radial action

Authors: Y. Xu, C. Liu, Z. Li, H. Tian, Sarah A. Bird, H. J. Newberg, S. Shao, L. C. Deng

Abstract: We use the conjugate angle of radial action ($θ_R$), the best representation of the orbital phase, to explore the "mid-plane, north branch, south branch" and "Monoceros area" disk structures that were previously revealed in the LAMOST K giants (Xu et al. 2020). The former three substructures, identified by their 3D kinematical distributions, have been shown to be projections of the phase space spi… ▽ More We use the conjugate angle of radial action ($θ_R$), the best representation of the orbital phase, to explore the "mid-plane, north branch, south branch" and "Monoceros area" disk structures that were previously revealed in the LAMOST K giants (Xu et al. 2020). The former three substructures, identified by their 3D kinematical distributions, have been shown to be projections of the phase space spiral (resulting from nonequilibrium phase mixing). In this work, we find that all of these substructures associated with the phase spiral show high aggregation in conjugate angle phase space, indicating that the clumping in conjugate angle space is a feature of ongoing, incomplete phase mixing. We do not find the $Z-V_Z$ phase spiral located in the "Monoceros area", but we do find a very highly concentrated substructure in the quadrant of conjugate angle space with the orbital phase from the apocenter to the guiding radius. The existence of the clump in conjugate angle space provides a complementary way to connect the "Monoceros area" with the direct response to a perturbation from a significant gravitationally interactive event. Using test particle simulations, we show that these features are analogous to disturbances caused by the impact of the last passage of the Sagittarius dwarf spheroidal galaxy. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: 53pages, 35 figures, 4 Tables, ApJ accepted

Enhanced Early Galaxy Formation in JWST from Axion Dark Matter?

Authors: Simeon Bird, Chia-Feng Chang, Yanou Cui, Daneng Yang

Abstract: We demonstrate that enhanced early galaxy formation can generically arise in axion-like particle (ALP) dark matter (DM) theories with the kinetic misalignment mechanism, potentially addressing the excess recently observed by the James Webb Space Telescope (JWST), while being consistent with existing constraints. We identify viable parameter regions with the ALP mass in the range of… ▽ More We demonstrate that enhanced early galaxy formation can generically arise in axion-like particle (ALP) dark matter (DM) theories with the kinetic misalignment mechanism, potentially addressing the excess recently observed by the James Webb Space Telescope (JWST), while being consistent with existing constraints. We identify viable parameter regions with the ALP mass in the range of $10^{-22}~{\rm eV}<m_a<10^{-19}~\rm eV$. In addition, we show that the ALP parameter regions of interest can lead to intriguing complementary signatures in the small-scale structure of DM halos and existing experimental searches for ALPs. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 6 pages, 5 figures

An Exploration of AGN and Stellar Feedback Effects in the Intergalactic Medium via the Low Redshift Lyman-$α$ Forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Sultan Hassan, Rachel S. Somerville, Romeel Davé, Federico Marinacci, Lars Hernquist, Mark Vogelsberger

Abstract: We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both… ▽ More We explore the role of galactic feedback on the low redshift Lyman-$α$ (Ly$α$)~forest ($z \lesssim 2$) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba sub-grid models. We find that both AGN and stellar feedback in Simba play a role in setting the Ly$α$ forest column density distribution function (CDD) and the Doppler width ($b$-value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM causing changes in the shape and normalization of the CDD and a broadening of the $b$-value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Ly$α$ forest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low-$z$ Ly$α$ forest can be sensitive to changes in the ultraviolet background (UVB), stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM. △ Less

Submitted 1 November, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 27 pages, 12 figures, 2 tables, accepted to ApJ

Age of FGK Dwarfs Observed with LAMOST and GALAH: Considering the Oxygen Enhancement

Authors: Tiancheng Sun, Zhishuai Ge, Xunzhou Chen, Shaolan Bi, Tanda Li, Xianfei Zhang, Yaguang Li, Yaqian Wu, Sarah A. Bird, Ferguson J. W., Jianzhao Zhou, Lifei Ye, Liu Long, Jinghua Zhang

Abstract: Varying oxygen abundance could impact the modeling-inferred ages. This work aims to estimate the ages of dwarfs considering observed oxygen abundance. To characterize 67,503 LAMOST and 4,006 GALAH FGK-type dwarf stars, we construct a grid of stellar models which take into account oxygen abundance as an independent model input. Compared with ages determined with commonly-used $α$-enhanced models, w… ▽ More Varying oxygen abundance could impact the modeling-inferred ages. This work aims to estimate the ages of dwarfs considering observed oxygen abundance. To characterize 67,503 LAMOST and 4,006 GALAH FGK-type dwarf stars, we construct a grid of stellar models which take into account oxygen abundance as an independent model input. Compared with ages determined with commonly-used $α$-enhanced models, we find a difference of $\sim$9% on average when the observed oxygen abundance is considered. The age differences between the two types of models are correlated to [Fe/H] and [O/$α$], and they are relatively significant on stars with [Fe/H] $\lesssim$ -0.6 dex. Generally, varying 0.2 dex in [O/$α$] will alter the age estimates of metal-rich (-0.2 $<$ [Fe/H] $<$ 0.2) stars by $\sim$10%, and relatively metal-poor (-1 $<$ [Fe/H] $<$ -0.2) stars by $\sim$15%. Of the low-O stars with [Fe/H] $<$ 0.1 dex and [O/$α$] $\sim$ -0.2 dex, many have fractional age differences of $\geq$ 10%, and even reach up to 27%. The fractional age difference of high-O stars with [O/$α$] $\sim$ 0.4 dex reaches up to -33% to -42% at [Fe/H] $\lesssim$ -0.6 dex. We also analyze the chemical properties of these stars. We find a decreasing trend of [Fe/H] with age from 7.5-9 Gyr to 5-6.5 Gyr for the stars from the LAMOST and GALAH. The [O/Fe] of these stars increases with decreasing age from 7.5-9 Gyr to 3-4 Gyr, indicating that the younger population is more O-rich. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: 18 pages, 12 figures, accepted for publication in ApJS

Fly-by galaxy encounters with multiple black holes produce star-forming linear wakes

Authors: Nianyi Chen, Patrick LaChance, Yueying Ni, Tiziana Di Matteo, Rupert Croft, Priyamvada Natarajan, Simeon Bird

Abstract: We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of blac… ▽ More We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes) reproduce remarkably well many of the key properties (its length and linearity; recent star formation, etc.) of the observed star-forming linear feature. We predict the feature to persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (with $M_{\rm gal}=10^9 \sim 10^{10}M_\odot$) and a dual BH system ($M_{\rm BH}=10^5 \sim 10^7\,M_\odot$) in its brightest knot. X-ray emission from AGN in the knot should be detectable in such systems. After $100\sim 200\,{\rm Myrs}$ from the first fly-by, the galaxies merge leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass $\sim 5\times 10^{10}\,M_\odot$. Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a fly-by encounter from a massive BH wake to reveal the origin of such features. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 8 pages, 5 figures, comments welcome

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σ_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology

Authors: Simeon Bird, Martin Fernandez, Ming-Feng Ho, Mahdi Qezlou, Reza Monadi, Yueying Ni, Nianyi Chen, Rupert Croft, Tiziana Di Matteo

Abstract: We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a… ▽ More We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$α$ forest 1D flux power spectra from SDSS and future DESI data releases. △ Less

Submitted 2 January, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 25 pages, 11 figures, matches version published in JCAP

Journal ref: JCAP 10 (2023) 037

MF-Box: Multi-fidelity and multi-scale emulation for the matter power spectrum

Authors: Ming-Feng Ho, Simeon Bird, Martin A. Fernandez, Christian R. Shelton

Abstract: We introduce MF-Box, an extended version of MFEmulator, designed as a fast surrogate for power spectra, trained using N-body simulation suites from various box sizes and particle loads. To demonstrate MF-Box's effectiveness, we design simulation suites that include low-fidelity suites (L1 and L2) at $256 \,\mathrm{Mpc}/h$ and $100 \,\mathrm{Mpc}/h$, each with $128^3$ particles, and a high-fidelity… ▽ More We introduce MF-Box, an extended version of MFEmulator, designed as a fast surrogate for power spectra, trained using N-body simulation suites from various box sizes and particle loads. To demonstrate MF-Box's effectiveness, we design simulation suites that include low-fidelity suites (L1 and L2) at $256 \,\mathrm{Mpc}/h$ and $100 \,\mathrm{Mpc}/h$, each with $128^3$ particles, and a high-fidelity suite (HF) with $512^3$ particles at $256 \,\mathrm{Mpc}/h$, representing a higher particle load compared to the low-fidelity suites. MF-Box acts as a probabilistic resolution correction function, learning most of the cosmological dependencies from L1 and L2 simulations and rectifying resolution differences with just 3 HF simulations using a Gaussian process. MF-Box successfully emulates power spectra from our HF testing set with a relative error of $< 3\%$ up to $k \simeq 7 \,h/\mathrm{Mpc}$ at $z \in [0, 3]$, while maintaining a cost similar to our previous multi-fidelity approach, which was accurate only up to $z = 1$. The addition of an extra low-fidelity node in a smaller box significantly improves emulation accuracy for MF-Box at $k > 2 \,h/\mathrm{Mpc}$, increasing it by a factor of $10$. We conduct an error analysis of MF-Box based on computational budget, providing guidance for optimizing budget allocation per fidelity node. Our proposed MF-Box enables future surveys to efficiently combine simulation suites of varying quality, effectively expanding the range of emulation capabilities while ensuring cost efficiency. △ Less

Submitted 4 October, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 19 pages, 16 figures. Match the MNRAS version. Data available at: https://github.com/jibanCat/MFBoxData, Code and tutorials available at: https://github.com/jibanCat/matter_emu_mfbox

Discovery of two rotational modulation periods from a young hierarchical triple system

Authors: Yu-Tao Chen, Hai-Jun Tian, Min Fang, Xiao-Xiong Zuo, Sarah A. Bird, Di Liu, Xin-Yu Zhu, Peng Zhang, Gao-Chao Liu, Sheng Cui

Abstract: GW~Ori is a young hierarchical triple system located in $λ$ Orionis, consisting of a binary (GW~Ori\,A and B), a tertiary star (GW~Ori\,C) and a rare circumtriple disk. Due to the limited data with poor accuracy, several short-period signals were detected in this system, but the values from different studies are not fully consistent. As one of the most successful transiting surveys, the Transiting… ▽ More GW~Ori is a young hierarchical triple system located in $λ$ Orionis, consisting of a binary (GW~Ori\,A and B), a tertiary star (GW~Ori\,C) and a rare circumtriple disk. Due to the limited data with poor accuracy, several short-period signals were detected in this system, but the values from different studies are not fully consistent. As one of the most successful transiting surveys, the Transiting Exoplanet Survey Satellite (TESS) provides an unprecedented opportunity to make a comprehensive periodic analysis of GW~Ori. In this work we discover two significant modulation signals by analyzing the light curves of GW~Ori's four observations from TESS, i.e., 3.02 $\pm$ 0.15\,d and 1.92 $\pm$ 0.06\,d, which are very likely to be the rotational periods caused by starspot modulation on the primary and secondary components, respectively. We calculate the inclinations of GW~Ori\,A and B according to the two rotational periods. The results suggest that the rotational plane of GW~Ori\,A and B and the orbital plane of the binary are almost coplanar. We also discuss the aperiodic features in the light curves; these may be related to unstable accretion. The light curves of GW~Ori also include a third (possible) modulation signal with a period of 2.51$\pm$0.09\,d, but the third is neither quite stable nor statistically significant. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 13 pages, 5 figures, 3 tables, accepted for publication by SCIENCE CHINA Physics, Mechanics & Astronomy

AI-assisted super-resolution cosmological simulations III: Time evolution

Authors: Xiaowen Zhang, Patrick Lachance, Yueying Ni, Yin Li, Rupert A. C. Croft, Tiziana Di Matteo, Simeon Bird, Yu Feng

Abstract: In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass fu… ▽ More In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass function agree well with results from high-resolution N-body simulations over the full trained redshift range ($10 \le z \le 0$). Furthermore, we assess the temporal consistency of our SR model by constructing halo merger trees. We examine progenitors, descendants and mass growth along the tree branches. All statistical indicators demonstrate the ability of our SR model to generate satisfactory high-resolution simulations based on low-resolution inputs. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Comments: 12 pages, 11 figures, code and movie available in https://github.com/sagasv5-xw/map2map on styled srsgan branch

Machine Learning Uncovers the Universe's Hidden Gems: A Comprehensive Catalogue of CIV Absorption Lines in SDSS DR12

Authors: Reza Monadi, Ming-Feng Ho, Kathy L. Cooksey, Simeon Bird

Abstract: We assemble the largest CIV absorption line catalogue to date, leveraging machine learning, specifically Gaussian processes, to remove the need for visual inspection for detecting CIV absorbers. The catalogue contains probabilities classifying the reliability of the absorption system within a quasar spectrum. Our training set was a sub-sample of DR7 spectra that had no detectable CIV absorption in… ▽ More We assemble the largest CIV absorption line catalogue to date, leveraging machine learning, specifically Gaussian processes, to remove the need for visual inspection for detecting CIV absorbers. The catalogue contains probabilities classifying the reliability of the absorption system within a quasar spectrum. Our training set was a sub-sample of DR7 spectra that had no detectable CIV absorption in a large visually inspected catalogue. We used Bayesian model selection to decide between our continuum model and our absorption-line models. Using a random hold-out sample of 1301 spectra from all of the 26,030 investigated spectra in DR7 CIV catalogue, we validated our pipeline and obtained an 87% classification performance score. We found good purity and completeness values, both ~80%, when a probability of ~95% is used as the threshold. Our pipeline obtained similar CIV redshifts and rest equivalent widths to our training set. Applying our algorithm to 185,425 selected quasar spectra from SDSS DR12, we produce a catalogue of 113,775 CIV doublets with at least 95% confidence. Our catalogue provides maximum a posteriori values and credible intervals for CIV redshift, column density, and Doppler velocity dispersion. We detect CIV absorption systems with a redshift range of 1.37 $\!-\!$ 5.1, including 33 systems with a redshift larger than 5 and 549 absorbers systems with a rest equivalent width greater than 2 A at more than 95% confidence. Our catalogue can be used to investigate the physical properties of the circumgalactic and intergalactic media. △ Less

Submitted 23 September, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 18 pages, 25 figures, 3 tables

The CAMELS project: Expanding the galaxy formation model space with new ASTRID and 28-parameter TNG and SIMBA suites

Authors: Yueying Ni, Shy Genel, Daniel Anglés-Alcázar, Francisco Villaescusa-Navarro, Yongseok Jo, Simeon Bird, Tiziana Di Matteo, Rupert Croft, Nianyi Chen, Natalí S. M. de Santi, Matthew Gebhardt, Helen Shao, Shivam Pandey, Lars Hernquist, Romeel Dave

Abstract: We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies.… ▽ More We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2,124 hydrodynamic simulation runs that vary 3 cosmological parameters ($Ω_m$, $σ_8$, $Ω_b$) and 4 parameters controlling stellar and AGN feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex non-linear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Boosting Line Intensity Map Signal-to-Noise with the Ly-$α$ Forest Cross-Correlation

Authors: Mahdi Qezlou, Simeon Bird, Adam Lidz, Guochao Sun, Andrew B. Newman, Gwen C. Rudie, Yueying Ni, Rupert Croft, Tiziana Di Matteo

Abstract: We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, an… ▽ More We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, and PFS. We show that CO$\times$Ly-$α$ forest can significantly enhance the detection signal-to-noise ratio of CO, with a $200$ to $300 \%$ improvement when cross-correlated with the forest observed in the Prime Focus Spectrograph (PFS) survey and a $50$ to $75\%$ enhancement for the currently available eBOSS or the upcoming DESI observations. We compare to the signal-to-noise improvements expected for a galaxy survey and show that CO$\times$Ly-$α$ is competitive with even a spectroscopic galaxy survey in raw signal-to-noise. Furthermore, our study suggests that the clustering of CO emission is tightly constrained by CO$\times$Ly-$α$ forest, due to the increased signal-to-noise ratio and the simplicity of Ly-$α$ absorption power spectrum modeling. Any foreground contamination or systematics are expected not to be shared between LIM surveys and Ly-$α$ forest observations; this provides an unbiased inference. Our findings highlight the potential benefits of utilizing the Ly-$α$ forest to aid in the initial detection of signals in line intensity experiments. For example, we also estimate that [CII]$\times$Ly-$α$ forest measurements from EXCLAIM and DESI/eBOSS, respectively, should have a larger signal-to-noise ratio than planned [CII]$\times$quasar observations by about an order of magnitude. Our results can be readily applied to actual data thanks to the observed quasar spectra in eBOSS Stripe 82, which overlaps with several LIM surveys. △ Less

Submitted 1 November, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: Codes and the produced data are available at https://github.com/qezlou/lali

Triple and Quadruple Black Holes in the ASTRID Simulation at $z \sim 2$

Authors: Calvin Hoffman, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft, Abraham Loeb

Abstract: We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the… ▽ More We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at $Δr \sim 50-100$ kpc. They include some of the most massive BHs (up to $10^{10} \,M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \,M_\odot$. Tuples' host galaxies are typically massive with $M_* \sim 10^{10-11} \,M_\odot$. We find that $>10\%$ massive halos with $M_{\rm halo} > 10^{13} M_\odot$ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5\%$ of the triplets all three MBHs merge within a Gyr, and $15\%$ go through one merger. As a by-product of the complex multi-galaxy interaction of these systems, we also find that up to $\sim 5\%$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($>10^{10} \,M_\odot $) have undergone a triple quasar phase while for BHs with $M_{\rm BH} \sim 10^9 \,M_\odot$ this fraction drops to $50\%$. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 10 pages, 9 figures; comments welcome

High-redshift supermassive black hole mergers in simulations with dynamical friction modelling

Authors: Colin DeGraf, Nianyi Chen, Yueying Ni, Tiziana Di Matteo, Simeon Bird, Michael Tremmel, Rupert Croft

Abstract: In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological si… ▽ More In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass BHs, and directly incorporates dynamical friction so as to accurately model the dynamical processes which bring black holes to the galaxy center where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to 10$^{4.5}$ M$_\odot$, and remain broadly consistent with the TNG simulations at scales > 10$^6$ M$_\odot$ (the seed mass used in TNG). By resolving lower-mass black holes, the overall merger rate is ~5x higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentering scheme, reducing the high-z merger rate mass-matched mergers by a factor of ~2x. We also calculate the expected LISA Signal-to-Noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA's peak sensitivity (10$^6$ M$_\odot$) to resolve the majority of LISA's GW detections. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 11 pages, 10 figures, submitted to MNRAS

A vast population of wandering and merging IMBHs at cosmic noon

Authors: Tiziana Di Matteo, Yueying Ni, Nianyi Chen, Rupert Croft, Simeon Bird, Fabio Pacucci, Angelo Ricarte, Michael Tremmel

Abstract: Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate… ▽ More Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate the population of IMBH seeds. Dynamical friction is largely inefficient at sinking and merging seed IMBHs at high-z. This leads to an extensive population (several hundred per galaxy) of wandering IMBHs in large halos at z~2. A small fraction of these IMBHs are detectable as HLXs, Hyper Luminous X-ray sources. Importantly, at z ~ 2, IMBHs mergers produce the peak of GW events. We find close to a million GW events in Astrid between z=2-3 involving seed IMBH mergers. These GW events (almost all detectable by LISA) at cosmic noon should provide strong constraints on IMBH seed models and their formation mechanisms. At the center of massive galaxies, where the number of IMBHs can be as high as 10-100, SMBH-IMBH pairs can form. These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio inspirals (EMRIs), will require the next generation of milli-muHz space-based GW interferometers to be detected. Large populations of IMBHs around massive black holes will probe their environments and MBH causal structure. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Efficient long-range AGN feedback affects the low redshift Lyman-$α$ forest

Authors: Megan Taylor Tillman, Blakesley Burkhart, Stephanie Tonnesen, Simeon Bird, Greg L. Bryan, Daniel Anglés-Alcázar, Romeel Davé, Shy Genel

Abstract: Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedbac… ▽ More Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$α$ forest can provide constraints on the impact of AGN feedback. We show that TNG over-predicts the number density of absorbers at column densities $N_{\rm HI} < 10^{14}$ cm$^{-2}$ compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model which includes AGN jet feedback is the preferred fit to the observed CDD of the $z = 0.1$ Lyman-$α$ forest across five orders of magnitude in column density. We show that the Simba results with jets produce a quantitatively better fit to the observational data than the Simba results without jets, even when the UVB is left as a free parameter. AGN jets in Simba are high speed, collimated, weakly-interacting with the interstellar medium (via brief hydrodynamic decoupling) and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the IGM when compared to TNG's kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low redshift Lyman-$α$ forest provides plausible evidence for long-range AGN jet feedback. △ Less

Submitted 7 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 15 pages, 5 figures, 1 table, accepted to ApJL

Journal ref: ApJL 945 L17 (2023)

A Forecast for Large Scale Structure Constraints on Horndeski Gravity with Line Intensity Mapping

Authors: Bryan R. Scott, Kirit S. Karkare, Simeon Bird

Abstract: We consider the potential for line intensity mapping (LIM) of the rotational CO(1-0), CO(2-1) and CO(3-2) transitions to detect deviations from General Relativity from $0 < z < 3$ within the framework of a very general class of modified gravity models, called Horndeski theories. Our forecast assumes a multi-tracer analysis separately obtaining information from the matter power spectrum and the fir… ▽ More We consider the potential for line intensity mapping (LIM) of the rotational CO(1-0), CO(2-1) and CO(3-2) transitions to detect deviations from General Relativity from $0 < z < 3$ within the framework of a very general class of modified gravity models, called Horndeski theories. Our forecast assumes a multi-tracer analysis separately obtaining information from the matter power spectrum and the first two multipoles of the redshift space distortion power spectrum. To achieve $\pm 0.1$ level constraints on the slope of the kinetic gravity braiding and Planck mass evolution parameters, a mm-wave LIM experiment would need to accumulate $\approx 10^8-10^9$ spectrometer hours, feasible with instruments that could be deployed in the 2030s. Such a measurement would constrain large portions of the remaining parameter space available to Scalar-Tensor modified gravity theories. Our modeling code is publicly available. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 13 pages, 5 figures; to be submitted to Monthly Notices of the Royal Astronomical Society

Report of the Topical Group on Cosmic Probes of Dark Matter for Snowmass 2021

Authors: Alex Drlica-Wagner, Chanda Prescod-Weinstein, Hai-Bo Yu, Andrea Albert, Mustafa Amin, Arka Banerjee, Masha Baryakhtar, Keith Bechtol, Simeon Bird, Simon Birrer, Torsten Bringmann, Regina Caputo, Sukanya Chakrabarti, Thomas Y. Chen, Djuna Croon, Francis-Yan Cyr-Racine, William A. Dawson, Cora Dvorkin, Vera Gluscevic, Daniel Gilman, Daniel Grin, Renée Hložek, Rebecca K. Leane, Ting S. Li, Yao-Yuan Mao , et al. (15 additional authors not shown)

Abstract: Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to… ▽ More Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to inform our understanding of the fundamental nature of dark matter in the coming decade. △ Less

Submitted 13 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Report of the CF3 Topical Group for Snowmass 2021; 35 pages, 10 figures, many references. V3 updates Fig 3-2 and the author list

Probing the Galactic halo with RR Lyrae stars $-$ III. The chemical and kinematic properties of the stellar halo

Authors: Gaochao Liu, Yang Huang, Sarah Ann Bird, Huawei Zhang, Fei Wang, Haijun Tian

Abstract: Based on a large spectroscopic sample of $\sim$ 4,300 RR Lyrae stars with metallicity, systemic radial velocity and distance measurements, we present a detailed analysis of the chemical and kinematic properties of the Galactic halo. Using this sample, the metallicity distribution function (MDF) as a function of $r$ and the velocity anisotropy parameter $β$ profiles (for different metallicity popul… ▽ More Based on a large spectroscopic sample of $\sim$ 4,300 RR Lyrae stars with metallicity, systemic radial velocity and distance measurements, we present a detailed analysis of the chemical and kinematic properties of the Galactic halo. Using this sample, the metallicity distribution function (MDF) as a function of $r$ and the velocity anisotropy parameter $β$ profiles (for different metallicity populations) are derived for the stellar halo. Both the chemical and kinematic results suggest that the Galactic halo is composed of two distinct parts, the inner halo and outer halo. The cutoff radius ($\sim$ 30 kpc) is similar to the previous break radius found in the density distribution of the stellar halo. We find that the inner part is dominated by a metal-rich population with extremely radial anisotropy ($β\sim 0.9$). These features are in accordance with those of ``{\it Gaia}-Enceladus-Sausage'' (GES) and we attribute this inner halo component as being dominantly composed of stars deposited from this ancient merged satellite. We find that GES probably has a slightly negative metallicity gradient. The metal-poor populations in the inner halo are characterized as a long-tail in MDF with an anisotropy of $β\sim 0.5$, which is similar to that of the outer part. The MDF for the outer halo is very broad with several weak peaks and the value of $β$ is around 0.5 for all metallicities. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 15 pages, 14 figures, MNRAS accepted

Ultramassive black holes formed by triple quasar mergers at $z\sim 2$

Authors: Yueying Ni, Tiziana DiMatteo, Nianyi Chen, Rupert Croft, Simeon Bird

Abstract: The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a br… ▽ More The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a bright, rare triple quasar system powered by three ~10^9 Msun black holes. The second merger of supermassive black holes (SMBHs) follows the first after 150 Myrs. The merger events lead to sustained Eddington accretion onto the central SMBH, forming an UMBH in the center of a massive compact stellar core with M* > 2e12 Msun. The strong feedback of the UMBH quenches the surrounding star formation to < 10 Msun/yr in the inner 50 kpc/h region. There are two more UMBHs with MBH > 5e10 Msun at z>2 in ASTRID which are also produced by major mergers of galaxies, and their progenitors can be observed as quasar triplets of lower luminosity. The rarely observed quasar multiples can be the cradle of UMBHs at high redshift, and likely end up in the center of the most massive clusters. △ Less

Submitted 6 September, 2022; v1 submitted 2 September, 2022; originally announced September 2022.

Comments: 7 pages, 4 figures, comments welcome

Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon

Authors: Yue Shen, Hsiang-Chih Hwang, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Nadia Zakamska, Xin Liu, Yu-Ching Chen, Kaitlin M. Kratter

Abstract: The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of… ▽ More The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5<z<3.5$, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases towards smaller separations by a factor of $\sim 5$ over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, sub-arcsec-resolution IR imaging for the closest pairs). Intriguingly, at this point the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 19 pages, 9 figures; submitted to ApJ

Journal ref: ApJ, 943, 38 (2023)

Properties and Evolution of Dual and Offset AGN in the ASTRID Simulation at $z \sim 2$

Authors: Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Michael Tremmel, Colin DeGraf, Yue Shen, A. Miguel Holgado, Simeon Bird, Rupert Croft, Yu Feng

Abstract: We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Com… ▽ More We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Compared with the full AGN population, duals are characterized by a low $M_\text{BH}/M_*$ ratio, a high specific star-formation rate (sSFR) of $\sim 1\,\text{Gyr}^{-1}$, and a high Eddington ratio ($\sim 0.05$, double that of single AGN). The dual AGN are formed in major galaxy mergers (typically involving $M_\text{halo}<10^{13}\,M_\odot$), with BHs that have similar masses. At small separations (when their host galaxies are in the late phase of the merger) duals become $2\sim8$ times brighter (albeit more obscured) than at larger separations. $80\%$ of these bright, close duals merge in the simulation within $\sim500\,\text{Myrs}$. Notably, the initially less-massive BH in duals frequently becomes the brighter AGN during the galaxy merger. In offset AGN, the active BH is typically $\gtrsim 10$ times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the center of massive halos ($ M_\text{ halo}\sim 10^{13-14}\,M_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $Δr\sim5\,\text{kpc}$ for a few hundred Myrs. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 17 Pages, 18 Figures; comments welcome!

Milky Way Mass with K Giants and BHB Stars Using LAMOST, SDSS/SEGUE, and Gaia: 3D Spherical Jeans Equation and Tracer Mass Estimator

Authors: Sarah A. Bird, Xiang-Xiang Xue, Chao Liu, Chris Flynn, Juntai Shen, Jie Wang, Chengqun Yang, Meng Zhai, Ling Zhu, Gang Zhao, Hai-Jun Tian

Abstract: We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it… ▽ More We measure the enclosed Milky Way mass profile to Galactocentric distances of $\sim70$ and $\sim50$ kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are LAMOST and SDSS/SEGUE K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and {\it Gaia} DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to $\sim70$ and $\sim50$ kpc, the Jeans method yields total enclosed masses of $4.3\pm0.95$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ and $4.1\pm1.2$ (random) $\pm0.6$ (systematic) $\times10^{11}$ M$_\odot$ for the KG and BHB stars, respectively. For the KG and BHB samples we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) $\pm0.083$ (systematic) $\times10^{12}$ M$_\odot$ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) $\pm0.15$ (systematic) $\times10^{12}$ M$_\odot$, respectively. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS. Comments welcomed

A Multi-Fidelity Emulator for the Lyman-$α$ Forest Flux Power Spectrum

Authors: M. A. Fernandez, Ming-Feng Ho, Simeon Bird

Abstract: In this work we extend our recently developed multi-fidelity emulation technique to the simulated Lyman-$α$ forest flux power spectrum. Multi-fidelity emulation allows interpolation of simulation outputs between cosmological parameters using many cheap low-fidelity simulations and a few expensive high-fidelity simulations. Using a test suite of small box (30 Mpc/h) simulations, we show that multi-… ▽ More In this work we extend our recently developed multi-fidelity emulation technique to the simulated Lyman-$α$ forest flux power spectrum. Multi-fidelity emulation allows interpolation of simulation outputs between cosmological parameters using many cheap low-fidelity simulations and a few expensive high-fidelity simulations. Using a test suite of small box (30 Mpc/h) simulations, we show that multi-fidelity emulation is able to reproduce the Lyman-$α$ forest flux power spectrum well, achieving an average accuracy when compared to a test suite of $0.8\%$. We further show that it has a substantially increased accuracy over single-fidelity emulators, constructed using either the high or low-fidelity simulations only. In particular, it allows the extension of an existing simulation suite to smaller scales and higher redshifts. △ Less

Submitted 24 October, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 12 pages, 10 figures

Journal ref: MNRAS 517, 3, 3200 (2022)

A population of ultraviolet-dim protoclusters detected in absorption

Authors: Andrew B. Newman, Gwen C. Rudie, Guillermo A. Blanc, Mahdi Qezlou, Simeon Bird, Daniel D. Kelson, Victoria Pérez, Enrico Congiu, Brian C. Lemaux, Alan Dressler, John S. Mulchaey

Abstract: Galaxy protoclusters, which will eventually grow into the massive clusters we see in the local universe, are usually traced by locating overdensities of galaxies. Large spectroscopic surveys of distant galaxies now exist, but their sensitivity depends mainly on a galaxy's star formation activity and dust content rather than its mass. Tracers of massive protoclusters that do not rely on their galax… ▽ More Galaxy protoclusters, which will eventually grow into the massive clusters we see in the local universe, are usually traced by locating overdensities of galaxies. Large spectroscopic surveys of distant galaxies now exist, but their sensitivity depends mainly on a galaxy's star formation activity and dust content rather than its mass. Tracers of massive protoclusters that do not rely on their galaxy constituents are therefore needed. Here we report observations of Lyman-$α$ absorption in the spectra of a dense grid of background galaxies, which we use to locate a substantial number of candidate protoclusters at redshifts 2.2-2.8 via their intergalactic gas. We find that the structures producing the most absorption, most of which were previously unknown, contain surprisingly few galaxies compared to the dark matter content of their analogs in cosmological simulations. Nearly all are expected to be protoclusters, and we infer that half of their expected galaxy members are missing from our survey because they are unusually dim at rest-frame ultraviolet wavelengths. We attribute this to an unexpectedly strong and early influence of the protocluster environment on the evolution of these galaxies that reduced their star formation or increased their dust content. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: Published in Nature. The version of record is available at https://www.nature.com/articles/s41586-022-04681-6 or view-only at https://rdcu.be/cPIIX

Journal ref: Nature 606, 475-478 (2022)

The low redshift Lyman-$α$ Forest as a constraint for models of AGN feedback

Authors: Blakesley Burkhart, Megan Tillman, Alexander B. Gurvich, Simeon Bird, Stephanie Tonnesen, Greg L. Bryan, Lars E. Hernquist, Rachel S. Somerville

Abstract: We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and… ▽ More We study the sensitivity of the $z=0.1$ Lyman-$α$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and similar cosmological parameters, but TNG features an entirely reworked AGN feedback model. Due to changes in the AGN radio mode model, the original Illustris simulations have a factor of 2-3 fewer Lyman-$α$ absorbers than TNG at column densities $N_{\rm HI}< 10^{15.5}$ cm$^{-2}$. We compare the simulated forest statistics to UV data from the Cosmic Origins Spectrograph (COS) and find that neither simulation can reproduce the slope of the absorber distribution. Both Illustris and TNG also produce significantly smaller line width distributions than observed in the COS data. We show that TNG is in much better agreement with the observed $z=0.1$ flux power spectrum than Illustris. We explore which statistics can disentangle the effects of AGN feedback from alternative UVB models by rescaling the UVB of Illustris to produce a CDD match to TNG. While this UVB rescaling is degenerate with the effect of AGN feedback on the CDD, the amplitude and shape of the flux PDF and 1D flux power spectrum change in a way distinct from a scaling of the UVB. Our study suggests that the $z=0.1$ Lyman-$α$ forest observables can be used as a diagnostic of AGN feedback models. △ Less

Submitted 18 July, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: ApJL

Snowmass2021 Cosmic Frontier White Paper:Primordial Black Hole Dark Matter

Authors: Simeon Bird, Andrea Albert, Will Dawson, Yacine Ali-Haimoud, Adam Coogan, Alex Drlica-Wagner, Qi Feng, Derek Inman, Keisuke Inomata, Ely Kovetz, Alexander Kusenko, Benjamin V. Lehmann, Julian B. Munoz, Rajeev Singh, Volodymyr Takhistov, Yu-Dai Tsai

Abstract: Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivit… ▽ More Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivity to PBHs. However, devoted analysis pipelines and theoretical modeling are required to fully leverage these novel data. The search for PBHs constitutes a low-cost, high-reward science case with significant impact on the high energy physics community. △ Less

Submitted 1 July, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 22 pages, 4 figures, submission to the Snowmass 2021 process. Update to add references

Snowmass2021 Cosmic Frontier White Paper: Dark Matter Physics from Halo Measurements

Authors: Keith Bechtol, Simon Birrer, Francis-Yan Cyr-Racine, Katelin Schutz, Susmita Adhikari, Mustafa Amin, Arka Banerjee, Simeon Bird, Nikita Blinov, Kimberly K. Boddy, Celine Boehm, Kevin Bundy, Malte Buschmann, Sukanya Chakrabarti, David Curtin, Liang Dai, Alex Drlica-Wagner, Cora Dvorkin, Adrienne L. Erickcek, Daniel Gilman, Saniya Heeba, Stacy Kim, Vid Iršič, Alexie Leauthaud, Mark Lovell , et al. (19 additional authors not shown)

Abstract: The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determin… ▽ More The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics. △ Less

Submitted 24 April, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). 88 pages, 9 figures. Comments welcome

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Authors: Yao-Yuan Mao, Annika H. G. Peter, Susmita Adhikari, Keith Bechtol, Simeon Bird, Simon Birrer, Jonathan Blazek, Jeffrey L. Carlin, Nushkia Chamba, Johann Cohen-Tanugi, Francis-Yan Cyr-Racine, Tansu Daylan, Birendra Dhanasingham, Alex Drlica-Wagner, Cora Dvorkin, Christopher Fassnacht, Eric Gawiser, Maurizio Giannotti, Vera Gluscevic, Alma Gonzalez-Morales, Renee Hlozek, M. James Jee, Stacy Kim, Akhtar Mahmood, Rachel Mandelbaum , et al. (8 additional authors not shown)

Abstract: Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requir… ▽ More Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation

Authors: Nianyi Chen, Yueying Ni, A. Miguel Holgado, Tiziana Di Matteo, Michael Tremmel, Colin DeGraf, Simeon Bird, Rupert Croft, Yu Feng

Abstract: We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to… ▽ More We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 17 pages, 13 Figures; comments are welcome

Characterizing Protoclusters and Protogroups at z $\sim$ 2.5 Using Ly$α$ Tomography

Authors: Mahdi Qezlou, Andrew B. Newman, Gwen C. Rudie, Simeon Bird

Abstract: Ly-$α$ tomography surveys have begun to produce three-dimensional (3D) maps of the intergalactic medium (IGM) opacity at $z \sim 2.5$ with Mpc resolution. These surveys provide an exciting new way to discover and characterize high-redshift overdensities, including the progenitors of today's massive groups and clusters of galaxies, known as protogroups and protoclusters. We use the IllustrisTNG-… ▽ More Ly-$α$ tomography surveys have begun to produce three-dimensional (3D) maps of the intergalactic medium (IGM) opacity at $z \sim 2.5$ with Mpc resolution. These surveys provide an exciting new way to discover and characterize high-redshift overdensities, including the progenitors of today's massive groups and clusters of galaxies, known as protogroups and protoclusters. We use the IllustrisTNG-$300$ hydrodynamical simulation to build mock maps that realistically mimic those observed in the Ly-$α$ Tomographic IMACS Survey (LATIS). We introduce a novel method for delineating the boundaries of structures detected in 3D Ly-$α$ flux maps by applying the watershed algorithm. We provide estimators for the dark matter masses of these structures (at $z\sim2.5$), their descendant halo masses at $z=0$, and the corresponding uncertainties. We also investigate the completeness of this method for the detection of protogroups and protoclusters. Compared to earlier work, we apply and characterize our method over a wider mass range that extends to massive protogroups. We also assess the widely used fluctuating Gunn-Peterson approximation (FGPA) applied to dark-matter-only simulations; we conclude while it is adequate for estimating the Ly-$α$ absorption signal from moderate-to-massive protoclusters ($ \gtrsim 10^{14.2} M_{\odot}$), it artificially merges a minority of lower-mass structures with more massive neighbors. Our methods will be applied to current and future Ly-$α$~tomography surveys to create catalogs of overdensities and study environment-dependent galaxy evolution in the Cosmic~Noon~era. △ Less

Submitted 11 May, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 18 pages, 18 figures, codes available at https://github.com/mahdiqezlou/LyTomo-Watershed

Journal ref: The Astrophysical Journal (2022), 930-109

Densified pupil spectrograph as high-precision radial velocimetry: From direct measurement of the Universe's expansion history to characterization of nearby habitable planet candidates

Authors: Taro Matsuo, Thomas P. Greene, Mahdi Qezlou, Simeon Bird, Kiyotomo Ichiki, Yuka Fujii, Tomoyasu Yamamuro

Abstract: The direct measurement of the Universe's expansion history and the search for terrestrial planets in habitable zones around solar-type stars require extremely high-precision radial velocity measures over a decade. This study proposes an approach for enabling high-precision radial velocity measurements from space. The concept presents a combination of a high-dispersion densified pupil spectrograph… ▽ More The direct measurement of the Universe's expansion history and the search for terrestrial planets in habitable zones around solar-type stars require extremely high-precision radial velocity measures over a decade. This study proposes an approach for enabling high-precision radial velocity measurements from space. The concept presents a combination of a high-dispersion densified pupil spectrograph and a novel telescope line-of-sight monitor. The precision of the radial velocity measurements is determined by combining the spectrophotometric accuracy and the quality of the absorption lines in the recorded spectrum. Therefore, a highly dispersive densified pupil spectrograph proposed to perform stable spectroscopy can be utilized for high-precision radial velocity measures. A concept involving the telescope line-of-sight monitor is developed to minimize the change of the telescope line-of-sight over a decade. This monitor allows the precise measurement of a long-term telescope drift without any significant impact on the Airy disk when the densified pupil spectra are recorded. We analytically derive the uncertainty of the radial velocity measurements, which is caused by the residual offset of the line-of-sights at two epochs. We find that the error could be reduced down to approximately 1 $cm/s$, and the precision will be limited by another factor (e.g., wavelength calibration uncertainty). A combination of the high precision spectrophotometry and the high spectral resolving power could open a new path toward the characterization of nearby non-transiting habitable planet candidates orbiting late-type stars. We present two simple and compact high-dispersed densified pupil spectrograph designs for the cosmology and exoplanet sciences. △ Less

Submitted 12 November, 2021; originally announced November 2021.

Comments: 37 pages, 18 figures, Accepted for publication in the Astronomical Journal

The ASTRID Simulation: Galaxy Formation and Reionization

Authors: Simeon Bird, Yueying Ni, Tiziana Di Matteo, Rupert Croft, Yu Feng, Nianyi Chen

Abstract: We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=3$. As a particular focus is modelling the hi… ▽ More We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=3$. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarise the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed UV luminosity functions, galaxy stellar mass functions and specific star formation rates. We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at $z=6$, halos with $M \sim 2\times 10^9 M_\odot$ which have been reionized have a star formation rate $1.5$ times greater than those which have not yet been reionized. △ Less

Submitted 23 March, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures. Accepted by MNRAS, added star formation rate density

Contribution of Gaia Sausage to the Galactic Stellar Halo Revealed by K Giants and Blue Horizontal Branch Stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, Sloan Digital Sky Survey, and Gaia

Authors: Wenbo Wu, Gang Zhao, Xiang-Xiang Xue, Sarah A. Bird, Chengqun Yang

Abstract: We explore the contribution of the Gaia-Sausage to the stellar halo of the Milky Way by making use of a Gaussian Mixture model (GMM) and applying it to halo star samples of LAMOST K giants, SEGUE K giants, and SDSS blue horizontal branch stars. The GMM divides the stellar halo into two parts, of which one represents a more metal-rich and highly radially biased component associated with an ancient,… ▽ More We explore the contribution of the Gaia-Sausage to the stellar halo of the Milky Way by making use of a Gaussian Mixture model (GMM) and applying it to halo star samples of LAMOST K giants, SEGUE K giants, and SDSS blue horizontal branch stars. The GMM divides the stellar halo into two parts, of which one represents a more metal-rich and highly radially biased component associated with an ancient, head-on collision referred to as the Gaia-Sausage, and the other one is a more metal-poor and isotropic halo. A symmetric bimodal Gaussian is used to describe the distribution of spherical velocity of the Gaia-Sausage, and we find that the mean absolute radial velocity of the two lobes decreases with Galactocentric radius. We find that the Gaia-Sausage contributes about $41\%-74\%$ of the inner (Galactocentric radius $r_\mathrm{gc} < 30$ kpc) stellar halo. The fraction of stars of the Gaia-Sausage starts to decline beyond $r_\mathrm{gc} \sim$ $25-30$ kpc, and the outer halo is found to be significantly less influenced by the Gaia-Sausage than the inner halo. After the removal of halo substructures found by integrals of motion, the contribution of the \textit{Gaia}-Sausage falls slightly within $r_\mathrm{gc} \sim$ 25 kpc, but is still as high as $30\%-63\%$. Finally, we select several possible Sausage-related substructures consisting of stars on highly eccentric orbits. The GMM/Sausage component agrees well with the selected substructure stars in their chemodynamical properties, which increases our confidence in the reliability of the GMM fits. △ Less

Submitted 12 January, 2022; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: accepted for publication in ApJ

The ASTRID simulation: the evolution of Supermassive Black Holes

Authors: Yueying Ni, Tiziana Di Matteo, Simeon Bird, Rupert Croft, Yu Feng, Nianyi Chen, Michael Tremmel, Colin DeGraf, Yin Li

Abstract: We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses… ▽ More We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers. △ Less

Submitted 20 March, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: Match the published version

The Stellar "Snake" I: Whole Structure and Properties

Authors: Fan Wang, Hai-Jun Tian, Dan Qiu, Qi Xu, Min Fang, Hao Tian, Di Li, Sarah Bird, Jian-Rong Shi, Xiao-Ting Fu, Gao-Chao Liu, Sheng Cui, Yong Zhang

Abstract: To complement our previous discovery of the young snake-like structure in the solar neighborhood and reveal the structure's full extent, we build two samples of stars within the Snake and its surrounding territory from {\tt Gaia EDR3}. With the friends-of-friends algorithm, we identify 2694 and 9615 Snake member candidates from the two samples. Thirteen open clusters are embedded in these member c… ▽ More To complement our previous discovery of the young snake-like structure in the solar neighborhood and reveal the structure's full extent, we build two samples of stars within the Snake and its surrounding territory from {\tt Gaia EDR3}. With the friends-of-friends algorithm, we identify 2694 and 9615 Snake member candidates from the two samples. Thirteen open clusters are embedded in these member candidates. By combining the spectroscopic data from multiple surveys, we investigate the comprehensive properties of the candidates and find that they \thj{are very likely to} belong to one sizable structure, since most of the components are well bridged in their spatial distributions, and follow a single stellar population with an age of $30-40$\,Myr and solar metallicity. This sizable structure is best explained as hierarchically primordial, and probably formed from a filamentary giant molecular cloud with unique formation history in localized regions. To analyze the dynamics of the Snake, we divide the structure into five groups according to their tangential velocities; we find that the groups are expanding at a coherent rate ($κ_X\sim3.0\,\times10^{-2}\,\rm km\,s^{-1}\,pc^{-1}$) along the length of the structure ($X$-direction). \thj{The corresponding expansion age ($τ\sim33$\,Myr) is highly consistent with the age of the Snake}. With over ten thousand member stars, the Snake is an ideal laboratory to study nearby coeval stellar formation, stellar physics, and environmental evolution over a large spatial extent. △ Less

Submitted 4 December, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 13 figures, 1 table, submitted to MNRAS (revised according to the referee's comments)