Cosmology

Abstract:

Cosmic ray (CR) feedback has been proposed as a powerful mechanism for driving warm gas outflows in galaxies. We use cosmological magnetohydrodynamic simulations to investigate the impact of CR feedback on neutral hydrogen (H i) in a 1011 M⊙ dark matter halo at 2 < z < 4. To this end, we postprocess the simulations with ionizing radiative transfer and perform Monte Carlo Lyman-α (Lyα) transfer calculations. CR feedback reduces H i column densities around young stars, thereby allowing more Lyα photons to escape and consequently offering a better match to the Lyα luminosities of observed Lyα emitters. Although galaxies with CR-driven outflows have more extended H i in the circumgalactic medium, two Lyα line properties sensitive to optical depth and gas kinematic—the location of the red peak relative to the Lyα line center in velocity space (vred) and relative strength of the blue-to-red peaks (B/R)—cannot distinguish between the CR-driven and non-CR simulations. This is because Lyα photons propagate preferentially along low H i density channels created by the ionizing radiation, thereby limiting the scattering with volume-filling H i. In contrast, the observed low flux ratios between the valley and peak and the surface brightness profiles are better reproduced in the model with CR-driven outflows because the Lyα photons interact more before escaping, rather than being destroyed by dust as is the case in the non-CR simulation. We discuss the potential cause of the paucity of sightlines in simulations that exhibit prominent red peaks and large vred, which may require the presence of more volume-filling H i.

Abstract:

Abstract We investigate the clustering of Lyman-break galaxies at redshifts of 3 ≲ z ≲ 5 within the COSMOS field by measuring the angular two-point correlation function. Our robust sample of ~60,000 bright (mUV ≲ 27) Lyman-break galaxies was selected based on spectral energy distribution fitting across 14 photometric bands spanning optical and near-infrared wavelengths. We constrained both the 1- and 2-halo terms at separations up to 300 arcsec, finding an excess in the correlation function at scales corresponding to &lt;20 kpc, consistent with enhancement due to clumps in the same galaxy or interactions on this scale. We then performed Bayesian model fits on the correlation functions to infer the Halo Occupation Distribution parameters, star formation duty cycle, and galaxy bias in three redshift bins. We examined several cases where different combinations of parameters were varied, showing that our data can constrain the slope of the satellite occupation function, which previous studies have fixed. For an MUV-limited sub-sample, we found galaxy bias values of $b_g=3.18^{+0.14}_{-0.14}$ at z ≃ 3, $b_g=3.58^{+0.27}_{-0.29}$ at z ≃ 4, $b_g=4.27^{+0.25}_{-0.26}$ at z ≃ 5. The duty cycle values are $0.62^{+0.25}_{-0.26}$, $0.40^{+0.34}_{-0.22}$, and $0.39^{+0.31}_{-0.20}$,respectively. These results suggest that, as the redshift increases, there is a slight decrease in the host halo masses and a shorter timescale for star formation in bright galaxies, at a fixed rest-frame UV luminosity threshold.

Abstract:

Abstract This study investigates the estimation of the neutral hydrogen (H i) mass function (HiMF) using a Bayesian stacking approach with simulated data for the Five-hundred-meter Aperture Spherical radio Telescope (FAST) H i intensity mapping (HiIM) drift-scan surveys. Using data from the IllustrisTNG simulation, we construct H i sky cubes at redshift z ∼ 0.1 and the corresponding optical galaxy catalogs, simulating FAST observations under various survey strategies, including pilot, deep-field, and ultradeep-field surveys. The HiMF is measured for distinct galaxy populations—classified by optical properties into red, blue, and bluer galaxies—and injected with systematic effects such as observational noise and flux confusion caused by the FAST beam. The results show that Bayesian stacking significantly enhances HiMF measurements. For red and blue galaxies, the HiMF can be well constrained with pilot surveys, while deeper surveys are required for the bluer galaxy population. Our analysis also reveals that sample variance dominates over observational noise, emphasizing the importance of wide-field surveys to improve constraints. Furthermore, flux confusion shifts the HiMF toward higher masses, which we address using a transfer function for correction. Finally, we explore the effects of intrinsic sample incompleteness and propose a framework to quantify its impact. This work lays the groundwork for future HiMF studies with FAST HiIM, addressing key challenges and enabling robust analyses of H i content across galaxy populations.

Abstract:

Abstract Using deep observations across three of the LOFAR Two-metre Sky Survey Deep Fields, this work measures the angular clustering of star forming galaxies (SFGs) and low-excitation radio galaxies (LERGs) to z ≲1.5 for faint sources, S144 MHz ≥200 μJy. We measure the angular auto-correlation of LOFAR sources in redshift bins and their cross-correlation with multi-wavelength sources to measure the evolving galaxy bias for SFGs and LERGs. Our work shows the bias of the radio-selected SFGs increases from $b=0.90^{+0.11}_{-0.10}$ at z ∼0.2 to $b=2.94^{+0.36}_{-0.36}$ at z ∼1.2; faster than the assumed b($z$)∝1/D($z$) models adopted in previous LOFAR cosmology studies (at sensitivities where AGN dominate), but in broad agreement with previous work. We further study the luminosity dependence of bias for SFGs and find little evidence for any luminosity dependence at fixed redshift, although uncertainties remain large for the sample sizes available. The LERG population instead shows a weaker redshift evolution with $b=2.33^{+0.28}_{-0.27}$ at z ∼0.7 to $b=2.65^{+0.57}_{-0.55}$ at z ∼1.2, though it is also consistent with the assumed bias evolution model (b($z$)∝1/D($z$)) within the measured uncertainties. For those LERGs which reside in quiescent galaxies (QLERGs), there is weak evidence that they are more biased than the general LERG population and evolve from b=$2.62^{+0.33}_{-0.33}$ at z ∼0.7 to $b=3.08^{+0.85}_{-0.84}$ at z ∼1.2. This suggests the halo environment of radio sources may be related to their properties. These measurements can help constrain models for the bias evolution of these source populations, and can help inform multi-tracer analyses.

Abstract:

Understanding the origin and evolution of supermassive black holes (SMBHs) stands as one of the most important challenges in astrophysics and cosmology, with little current theoretical consensus. Improved observational constraints on the cosmological evolution of SMBH demographics are needed. Here we report results of a search via photometric variability for SMBHs appearing as active galactic nuclei (AGN) in the cosmological volume defined by the Hubble Ultra Deep Field. This work includes particular focus on a new observation carried out in 2023 with the Hubble Space Telescope using the WFC3/IR/F140W, which is compared directly to equivalent data taken 11 yr earlier in 2012. Two earlier pairs of observations from 2009 to 2012 with WFC3/IR/F105W and WFC3/IR/F160W are also analyzed. We identify 521, 188, and 109 AGN candidates as nuclear sources that exhibit photometric variability at a level of 2σ, 2.5σ, and 3σ, respectively, in at least one filter. This sample includes 13, 3, and 2 AGN candidates at redshifts z > 6, when the Universe was ≲900 Myr old. After variability and luminosity function (down to MUV = −17 mag) completeness corrections, we estimate the comoving number density of SMBHs, nSMBH(z). At z ≳ 6, nSMBH ≳ 6 × 10−3 cMpc−3. At low z our observations are sensitive to AGN fainter than MUV = −17 mag, and we estimate nSMBH ≳ 10−2 cMpc−3. We discuss how these results place strong constraints on a variety of SMBH seeding theories.