Search

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., Codes and the produced data are available at https://github.com/qezlou/lali

Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, $k$, we find that the baryon fraction in halos contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel'dovich (tSZ) signal from group-scale halos contains similar information. Using recent Dark Energy Survey (DES) weak lensing and Atacama Cosmology Telescope (ACT) tSZ cross-correlation measurements and models trained on CAMELS, we obtain $10\%$ constraints on feedback effects on the power spectrum at $k \sim 5\, h/{\rm Mpc}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal{O}(, 18 pages, 15 figures. Comments are welcome

We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $Ω_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $Ω_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $Ω_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $Ω_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$., 32 pages, 13 figures, summary video: https://youtu.be/STZHvDHkVgo

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 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., 11 pages, 10 figures. Sumbitted to Monthly Notices of the Royal Astronomical Society

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $\Omega_{\rm m}$ and $\sigma_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $\Omega_{\rm m}$ and $\sigma_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $\Omega_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters., Comment: 25 pages, 11 figures, summary video: https://youtu.be/qkw92Z6owJU

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.

Recently, a family of hairy black holes in 4-dimensional Einstein gravity minimally coupled to a complex, massive scalar field was discovered~\cite{hbh}. Besides the mass $M$ and spin angular momentum $J$, these objects are characterized by a Noether charge $Q$, measuring the amount of scalar hair, which is not associated to a Gauss law and cannot be measured at spatial infinity. Introducing a dimensionless scalar hair parameter $q$, ranging from 0 to 1, we recover (a subset of) Kerr black holes for $q=0$ and a family of rotating boson stars for $q=1$. In the present paper, we explore the possibility of measuring $q$ for astrophysical black holes with current and future X-ray missions. We study the iron K$��$ line expected in the reflection spectrum of such hairy black holes and we simulate observations with Suzaku and eXTP. As a proof of concept, we point out, by analyzing a sample of hairy black holes, that current observations can already constrain the scalar hair parameter $q$, because black holes with $q$ close to 1 would have iron lines definitively different from those we observe in the available data. We conclude that a detailed scanning of the full space of solutions, together with data from the future X-ray missions, like eXTP, will be able to put relevant constraints on the astrophysical realization of Kerr black holes with scalar hair., 1+18 pages, 9 figures. v2: refereed version