Cosmic Density Research Articles

The baryon census and the mass-density of stars, neutral gas, and hot gas as a function of halo mass

Abstract We study the stellar, neutral gas content within halos over a halo mass range 1010to1015.5M⊙ and hot X-ray gas content over a halo mass range 1012.8to1015.5M⊙ in the local universe. We combine various empirical datasets of stellar, H i and X-ray observations of galaxies, groups and clusters to establish fundamental baryonic mass vs halo mass scaling relations. These scaling relations are combined with halo mass function to obtain the baryon densities of stars, neutral gas and hot gas (T > 106K), as a function of halo mass. We calculate the contributions of the individual baryonic components to the cosmic baryon fraction. Cosmic stellar mass density ($\Omega _\text{star}=2.09^{+0.21}_{-0.18} \times 10^{-3}$), cosmic H i mass density ($\Omega _\rm{H\,{\small I}}=0.49^{+0.25}_{-0.12} \times 10^{-3}$) and cosmic neutral gas mass density ($\Omega _\text{neutral gas}=0.71^{+0.39}_{-0.18} \times 10^{-3}$) estimates are consistent with previous more direct method measurements of these values, thereby establishing the veracity of our method. We also give an estimate of the cosmic hot plasma density ($\Omega _\text{hot gas}=2.58^{+2.1}_{-0.66} \times 10^{-3}$).

Intertwined formation of H2, dust, and stars in cosmological simulations

Context. Molecular hydrogen (H2) plays a crucial role in the formation and evolution of galaxies, serving as the primary fuel reservoir for star formation. In a metal-enriched Universe, H2 forms mostly through catalysis on interstellar dust grain surfaces. However, due to the complexities of modelling this process, star formation in cosmological simulations often relies on empirical or theoretical frameworks that have only been validated in the local Universe to estimate the abundance of H2. Aims. The goal of this work is to model the connection between the processes of star, dust, and H2 formation in our cosmological simulations. Methods. Building upon our recent integration of a dust evolution model into the star formation and feedback model MUPPI, we included the formation of molecular hydrogen on the surfaces of dust grains. We also accounted for the destruction of molecules and their shielding from harmful radiation. Results. The model reproduces, reasonably well, the main statistical properties of the observed galaxy population for the stellar, dust, and H2 components. The evolution of the molecular hydrogen cosmic density (ρH2) in our simulated boxes peaks around redshift z = 1.5, consistent with observations. Following its peak, ρH2 decreases by a factor of two towards z = 0, which is a milder evolution than observed. Similarly, the evolution of the molecular hydrogen mass function since z = 2 displays a gentler evolution when compared to observations. Our model recovers satisfactorily the integrated molecular Kennicut-Schmidt (mKS) law between the surface star formation rate (ΣSFR) and surface H2 density (ΣH2) at z = 0. This relationship is already evident at z = 2, albeit with a higher normalization. We find hints of a broken power law with a steeper slope at higher ΣH2. We also study the H2-to-dust mass ratio in galaxies as a function of their gas metallicity and stellar mass, observing a decreasing trend with respect to both quantities. The H2-to-dust mass fraction for the global population of galaxies is higher at higher redshift. The analysis of the atomic-to-molecular transition on a particle-by-particle basis suggests that gas metallicity cannot reliably substitute the dust-to-gas ratio in models attempting to simulate dust-promoted H2.

ASW^2DF: Census of the obscured star formation in a galaxy cluster in formation at z=2.2

We report the results of the deep and wide Atacama Large Millimeter/submillimeter Array (ALMA) 1.2\,mm mapping of the Spiderweb protocluster at $z=2.16$. The observations were divided into six contiguous fields covering a survey area of 19.3\,arcmin$^2$. With sim 13h of on-source time, the final maps in the six fields reach the 1sigma rms noise in a range of $40.3-57.1\ at a spatial resolution of $0 By using different source extraction codes and careful visual inspection, we detected 47 ALMA sources at a significance higher than 4sigma . We constructed the differential and cumulative number counts down to $ after the correction for purity and completeness obtained from Monte Carlo simulations. The ALMA 1.2\,mm number counts of dusty star-forming galaxies (DSFGs) in the Spiderweb protocluster are overall two times that of general fields, with some regions showing even higher overdensities (more than a factor of three). This is consistent with the results from previous studies over a larger scale using single-dish instruments. Comparison of the spatial distributions between different populations indicates that our ALMA sources are likely drawn from the same distribution as CO(1-0) emitters from the COALAS large program but are distinct from that of Halpha emitters. The cosmic star formation rate density of the ALMA sources is consistent with previous results (e.g., LABOCA 870\,mu m observations) after accounting for the difference in volume. We show that molecular gas masses estimates from dust measurements are not consistent with the ones derived from CO(1-0) and thus have to be taken with caution. The multiplicity fraction of single-dish DSFGs is higher than that of the field. Moreover, two extreme concentrations of ALMA sources were found on the outskirts of the Spiderweb protocluster, with an excess of more than 12 times that of the general fields. These results indicate that the ALMA-detected DSFGs are supplied through gas accretion along filaments and are triggered by intense star formation by accretion shocks before falling into the cluster center. The identified two galaxy groups are likely falling into the protocluster center and will trigger new merger events eventually, as indicated in simulations.

JWST FRESCO: a comprehensive census of H β + [O iii] emitters at 6.8 < z < 9.0 in the GOODS fields

ABSTRACT We present the census of H $\beta$ + [${\rm O\, {\small III}}$] $4960,5008\rm{\mathring{\rm\,A}}$ emitters at $6.8\lt z\lt 9.0$ from the JWST FRESCO survey over 124 arcmin$^2$ in the GOODS-North and GOODS-South fields. Our unbiased spectroscopic search results in 137 spectroscopically confirmed galaxies at $6.8\lt z\lt 9.0$ with observed [${\rm O\, {\small III}}$] fluxes $f_{\rm { [{\rm O\, {\small III}}]}}\gtrsim 1\times 10^{-18}\ \rm {ergs}\ \rm {s}^{-1} \ \rm {cm}^{-2}$. The rest-frame optical line ratios of the median stacked spectrum (median $M_{\rm {UV}}=-19.65^{+0.59}_{-1.05}$) indicate negligible dust attenuation, low metallicity ($12+\log (\rm {O/H})= 7.2-7.7$) and a high ionization parameter $\log _{10}U \simeq -2.5$. We find a factor $\times 1.3$ difference in the number density of $6.8\lt z\lt 9.0$ galaxies between GOODS-South and GOODS-North, which is caused by a single overdensity at $7.0\lt z\lt 7.2$ in GOODS-North. The bright end of the UV luminosity function of spectroscopically confirmed [${\rm O\, {\small III}}$] emitters is in good agreement with HST dropout-selected samples. Discrepancies between the observed [${\rm O\, {\small III}}$] LF, [${\rm O\, {\small III}}$]/UV ratio, and [${\rm O\, {\small III}}$] equivalent widths, and that predicted by theoretical models, suggest burstier star-formation histories and/or more heterogeneous metallicity and ionizing conditions in $z\gt 7$ galaxies. We report a rapid decline of the [${\rm O\, {\small III}}$] luminosity density at $z\gtrsim 6\!-\!7$ which cannot be explained by the evolution of the cosmic star-formation rate density. Finally we find that FRESCO detects in only 2h galaxies likely accounting for $\sim 10-20{{\ \rm per\ cent}}$ of the ionizing budget at $z=7\!-\!8$ (assuming an escape fraction of $10{{\ \rm per\ cent}}$), raising the prospect of directly detecting a significant fraction of the sources of reionization with JWST.

The Calibration of Polycyclic Aromatic Hydrocarbon Dust Emission as a Star Formation Rate Indicator in the AKARI NEP Survey

Polycyclic aromatic hydrocarbon (PAH) dust emission has been proposed as an effective extinction-independent star formation rate (SFR) indicator in the mid-infrared, but this may depend on conditions in the interstellar medium. The coverage of the AKARI/Infrared Camera (IRC) allows us to study the effects of metallicity, starburst intensity, and active galactic nuclei on PAH emission in galaxies with f ν (L18W) ≲ 19 AB mag. Observations include follow-up, rest-frame optical spectra of 443 galaxies within the AKARI North Ecliptic Pole survey that have IRC detections from 7 to 24 μm. We use optical emission line diagnostics to infer SFR based on Hα and [O ii]λ λ3726, 3729 emission line luminosities. The PAH 6.2 μm and PAH 7.7 μm luminosities (L(PAH 6.2 μm) and L(PAH 7.7 μm), respectively) derived using multiwavelength model fits are consistent with those derived from slitless spectroscopy within 0.2 dex. L(PAH 6.2 μm) and L(PAH 7.7 μm) correlate linearly with the 24 μm dust-corrected Hα luminosity only for normal, star-forming “main-sequence” galaxies. Assuming multilinear correlations, we quantify the additional dependencies on metallicity and starburst intensity, which we use to correct our PAH SFR calibrations at 0 < z < 1.2 for the first time. We derive the cosmic star formation rate density (SFRD) per comoving volume from 0.15 ≲ z ≲ 1. The PAH SFRD is consistent with that of the far-infrared and reaches an order of magnitude higher than that of uncorrected UV observations at z ∼ 1. Starburst galaxies contribute ≳0.7 of the total SFRD at z ∼ 1 compared to main-sequence galaxies.

Constraining the Initial Mass Function via Stellar Transients

The stellar initial mass function (IMF) represents a fundamental quantity in astrophysics and cosmology describing the mass distribution of stars from low mass all the way up to massive and very massive stars. It is intimately linked to a wide variety of topics, including stellar and binary evolution, galaxy evolution, chemical enrichment, and cosmological reionization. Nonetheless, the IMF still remains highly uncertain. In this work, we aim to determine the IMF with a novel approach based on the observed rates of transients of stellar origin. We parametrize the IMF with a simple but flexible Larson shape, and insert it into a parametric model for the cosmic UV luminosity density, local stellar mass density, type Ia supernova (SN Ia), core-collapse supernova (CCSN), and long gamma-ray burst (LGRB) rates as a function of redshift. We constrain our free parameters by matching the model predictions to a set of empirical determinations for the corresponding quantities via a Bayesian Markov Chain Monte Carlo method. Remarkably, we are able to provide an independent IMF determination with a characteristic mass mc=0.10−0.08+0.24M⊙ and high-mass slope ξ=−2.53−0.27+0.24 that are in accordance with the widely used IMF parameterizations (e.g., Salpeter, Kroupa, Chabrier). Moreover, the adoption of an up-to-date recipe for the cosmic metallicity evolution allows us to constrain the maximum metallicity of LGRB progenitors to Zmax=0.12−0.05+0.29Z⊙. We also find which progenitor fraction actually leads to SN Ia or LGRB emission (e.g., due to binary interaction or jet-launching conditions), put constraints on the CCSN and LGRB progenitor mass ranges, and test the IMF universality. These results show the potential of this kind of approach for studying the IMF, its putative evolution with the galactic environment and cosmic history, and the properties of SN Ia, CCSN, and LGRB progenitors, especially considering the wealth of data incoming in the future.

Dispersion and rotation measures from fast radio burst (FRB) host galaxies based on the TNG50 simulation

Context. Fast radio bursts (FRBs) are poised to become important cosmological tools in the near future, as the number of observed FRBs is increasing rapidly with multiple surveys underway. A large sample of FRBs will soon have available dispersion measures (DMs) and rotation measures (RMs), which can be used to study the cosmic baryon density and the intergalactic magnetic field. However, the observed DM and RM of FRBs consists of multiple contributions that must be quantified to estimate the DM and RM of the intergalactic medium (IGM). Aims. In this paper, we estimate one such contribution to DM and RM, namely, of FRB host galaxies. We show how this contribution changes with redshift, galaxy type, and the stellar mass of the galaxies. We also investigate its dependence on galaxy inclination and on an FRB’s offset from the center of the galaxy. Methods. Using the TNG50 simulation of the IllustrisTNG project, we selected 16 500 galaxies at redshifts of 0≤ ɀ ≤2, with stellar masses in the range of 9 ≤ log(M*/M⊙) ≤ 12. In each galaxy, we calculated the DM and RM contributions of 1000 sightlines; from these, we constructed the DM and RM probability density functions (PDFs). Results. We find that the rest frame DM distributions of all galaxies at a given redshift can be fitted by a log normal function and its median and width increase as a function of redshift. The rest-frame RM distribution is symmetric, with a median RMhost,rf=0 rad m–2 and it can be fitted by a combination of a Lorentzian and two Gaussian functions. The redshift evolution of the distribution width can be fitted by a curved power law. The parameters of these functions change for different subsets of galaxies with different stellar mass, inclination, and FRB offset. These changes are due to an increasing ne with redshift, SFR, and stellar mass. We do find a more ordered B field at lower ɀ compared to higher ɀ, as suggested by the presence of more galaxies with B field reversals and B fields dominated by random B field at higher ɀ. Conclusions. We estimated the FRB host DM and RM contributions, which can be used in the future to isolate the IGM contribution from the observed DM and RM of FRBs. We predict that to constrain a σRM,IGM of 2 rad m–2 to the 95% confidence level, we would need to observe 95 000 FRBs at ɀ = 0.5, but only 9 500 FRBs at ɀ = 2.

Deep learning the intergalactic medium using lyman-alpha forest at 4 ≤ z ≤ 5

Abstract Unveiling the thermal history of the intergalactic medium (IGM) at 4 ≤ z ≤ 5 holds the potential to reveal early onset He ii reionization or lingering thermal fluctuations from H i reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep neural networks. Our approach leverages the Sherwood-Relics simulation suite, consisting of diverse thermal histories, to generate mock spectra. Our convolutional and residual networks with likelihood metric predicts the Lyα optical depth-weighted density or temperature for each pixel in the Lyα forest skewer. We find that our network can successfully reproduce IGM conditions with high fidelity across range of instrumental signal-to-noise. These predictions are subsequently translated into the temperature-density plane, facilitating the derivation of reliable constraints on thermal parameters. This allows us to estimate temperature at mean cosmic density, T0 , with one sigma confidence, $\delta {\rm T_{\rm 0}}\,$≲ 1000K, using only one 20h−1cMpc sightline (Δz ≃ 0.04) with a typical reionization history. Existing studies utilize redshift pathlength comparable to Δz ≃ 4 for similar constraints. We can also provide more stringent constraints on the slope (1σ confidence interval, δγ ≲ 0.1) of the IGM temperature-density relation as compared to other traditional approaches. We test the reconstruction on a single high signal-to-noise observed spectrum (20h−1cMpc segment), and recover thermal parameters consistent with current measurements. This machine learning approach has the potential to provide accurate yet robust measurements of IGM thermal history at the redshifts in question.

Redshift-dependent galaxy formation efficiency at z = 5 − 13 in the FirstLight Simulations

Context. Some models of the formation of first galaxies predict low masses and faint objects at extremely high redshifts, z ≃ 9 − 15. However, the first observations of this epoch indicate a higher-than-expected number of bright (sometimes massive) galaxies. Aims. Numerical simulations can help to elucidate the mild evolution of the bright end of the UV luminosity function and they can provide the link between the evolution of bright galaxies and variations of the galaxy formation efficiency across different redshifts. Methods. We use the FirstLight database of 377 zoom-in cosmological simulations of a volume- and mass-complete sample of galaxies. Mock luminosities are estimated by a dust model constrained by observations of the β–MUV relation at z = 6 − 9. Results. FirstLight contains a high number of bright galaxies, MUV ≤ −20, consistent with current data at z = 6 − 13. The evolution of the UV cosmic density is driven by the evolution of the galaxy efficiency and the relation between MUV and halo mass. The efficiency of galaxy formation increases significantly with mass and redshift. At a fixed mass, galactic halos at extremely high redshifts convert gas into stars at a higher rate than at lower redshifts. The high gas densities in these galaxies enable high efficiencies. Our simulations predict higher number densities of massive galaxies, M* ≃ 109 M⊙, than other models with constant efficiency. Conclusions. Cosmological simulations of galaxy formation with detailed models of star formation and feedback can reproduce the different regimes of galaxy formation across cosmic history.

JWST PRIMER: a new multifield determination of the evolving galaxy UV luminosity function at redshifts z ≃ 9 – 15

ABSTRACT We present a new determination of the evolving galaxy ultraviolet (UV) luminosity function (LF) over the redshift range $8.5&amp;lt; z&amp;lt; 15.5$ using a combination of several major Cycle-1 JWST imaging programmes – Public Release IMaging for Extragalactic Research, JWST Advanced Deep Extragalactic Survey, and Next Generation Deep Extragalactic Exploratory Public Survey. This multifield approach yields a total of $\simeq 370$ arcmin2 of JWST/NIRCam imaging, reaching (5-$\sigma$) depths of $\simeq 30$ AB mag in the deepest regions. We select a sample of 2548 galaxies with a significant probability of lying at high redshift ($p(z&amp;gt; 8.5)&amp;gt; 0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq 4$ mag in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($\rho _{\rm {SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the observed early evolution of the galaxy UV LF (and $\rho _{\rm {SFR}}$) can be reproduced in a ${\rm \Lambda }$cold dark matter Universe, with no change in dust properties or star-formation efficiency required out to $z \simeq 12$. Instead, a progressive trend towards younger stellar population ages can reproduce the observations, and the typical ages required at $z \simeq$ 8, 9, 10, and 11 all converge on $\simeq 380-330$ Myr after the big bang, indicative of a rapid emergence of early galaxies at $z \simeq 12 - 13$. This is consistent with the first indications of a steeper drop-off in $\rho _{\rm {SFR}}$ we find beyond $z \simeq 13$, possibly reflecting the rapid evolution of the halo mass function at earlier times.

Gravitational Particle Production and the Hubble Tension

The effect of gravitational particle production of scalar particles on the total effective cosmic energy density (in the era after photon decoupling till the present) is considered. The effect is significant for heavy particles. It is found that gravitational particle production results in an effective increase in the directly measured value of the Hubble constant H0, while it does not affect the value of the Hubble constant in the calculation of the number density of baryons at the present time that is used to calculate recombination redshift. This may explain why the Hubble constants determined by local measurements and non-local measurements (such as CMB) are different. This suggests that gravitational particle production may have a non-negligible impact on H0 tension.

The detection and characterization of highly magnified stars with JWST: Prospects of finding Population III

Abstract Gravitational lensing may render individual high-mass stars detectable out to cosmological distances, and several extremely magnified stars have in recent years been detected out to redshifts z ≈ 6. Here, we present Muspelheim, a model for the evolving spectral energy distributions of both metal-enriched and metal-free stars at high redshifts. Using this model, we argue that lensed stars will form a highly biased sample of the intrinsic distribution of stars across the Hertzsprung-Russell diagram, and that this bias will typically tend to favour the detection of lensed stars in evolved stages characterized by low effective temperatures, even though stars only spend a minor fraction of their lifetimes in such states. We also explore the prospects of detecting individual, lensed metal-free (Population III) stars at high redshifts using the James Webb Space Telescope (JWST). We find that very massive (≳ 100 M⊙) Population III stars at z ≳ 6 may potentially be detected by JWST in surveys covering large numbers of strong lensing clusters, provided that the Population III stellar initial mass function is sufficiently top-heavy, that these stars evolve to effective temperatures ≤15000 K, and that the cosmic star formation rate density of Pop III stars reaches ≳ 10−4 M⊙ cMpc−3 yr−1 at z ≈ 6–10. Various ways to distinguish metal-free lensed stars from metal-enriched ones are also discussed.

Filaments of the Slime Mold Cosmic Web and How They Affect Galaxy Evolution

We present a novel approach for identifying cosmic web filaments within the DisPerSE structure identification framework, using cosmic density field estimates from the Monte Carlo Physarum Machine (MCPM), inspired by the slime mold organism. We apply our method to the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. The MCPM density field is superior to the Delaunay tessellation field estimator in tracing the true underlying matter distribution and allows filaments to be identified with higher fidelity, finding more low-prominence/diffuse filaments. Using our new filament catalogs, we find that ≳90% of galaxies are located within ∼1.5 Mpc of a filamentary spine, with little change in the median star formation activity with distance to the nearest filament. Instead, we uncover a differential effect of the local filament line density, Σfil (MCPM)—the total MCPM overdensity per unit length along a filament segment—on galaxy formation: most galaxies are quenched and gas-poor near high-line density filaments at z ≤ 1. At earlier times, the filamentary environment appears to have no effect on galactic gas supply and quenching. At z = 0, quenching in log(M*/M⊙)≳10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. Satellites are far more susceptible to filaments than centrals. The local environments of massive halos are not sufficient to account for the effect of filament line density on gas removal and quenching. Our new approach holds great promise for observationally identifying filaments from galaxy surveys such as SDSS and DESI.

The FLAMINGO simulation view of cluster progenitors observed in the epoch of reionization with JWST

ABSTRACT Motivated by the recent JWST discovery of galaxy overdensities during the Epoch of Reionzation, we examine the physical properties of high-z protoclusters and their evolution using the Full-hydro Large-scale structure simulations with All-sky Mapping for the Interpretation of Next Generation Observations (FLAMINGO) simulation suite. We investigate the impact of the apertures used to define protoclusters, because the heterogeneous apertures used in the literature have limited our understanding of the population. Our results are insensitive to the uncertainties of the subgrid models at a given resolution, whereas further investigation into the dependence on numerical resolution is needed. When considering galaxies more massive than $M_\ast \, {\simeq }\, 10^8\, {\rm M_\odot }$, the FLAMINGO simulations predict a dominant contribution from progenitors similar to those of the Coma cluster to the cosmic star formation rate density during the reionization epoch. Our results indicate the onset of suppression of star formation in the protocluster environments as early as $z\, {\simeq }\, 5$. The galaxy number density profiles are similar to NFW (Navarro–Frenk–White profile) at $z\, {\lesssim }\, 1$ while showing a steeper slope at earlier times before the formation of the core. Different from most previous simulations, the predicted star formation history for individual protoclusters is in good agreement with observations. We demonstrate that, depending on the aperture, the integrated physical properties including the total (dark matter and baryonic) mass can be biased by a factor of 2 to 5 at $z\, {=}\, 5.5$–7, and by an order of magnitude at $z\, {\lesssim }\, 4$. This correction suffices to remove the ${\simeq }\, 3\, \sigma$ tensions with the number density of structures found in recent JWST observations.

Limits on non-relativistic matter during Big-bang nucleosynthesis

Big-bang nucleosynthesis (BBN) probes the cosmic mass-energy density at temperatures ∼ 10 MeV to ∼ 100 keV. Here, we consider the effect of a cosmic matter-like species that is non-relativistic and pressureless during BBN. Such a component must decay; doing so during BBN can alter the baryon-to-photon ratio, η, and the effective number of neutrino species. We use light element abundances and the cosmic microwave background (CMB) constraints on η and Nν to place constraints on such a matter component. We find that electromagnetic decays heat the photons relative to neutrinos, and thus dilute the effective number of relativistic species to N eff < 3 for the case of three Standard Model neutrino species. Intriguingly, likelihood results based on Planck CMB data alone find Nν = 2.800 ± 0.294, and when combined with standard BBN and the observations of D and 4He give Nν = 2.898 ± 0.141. While both results are consistent with the Standard Model, we find that a nonzero abundance of electromagnetically decaying matter gives a better fit to these results. Our best-fit results are for a matter species that decays entirely electromagnetically with a lifetime τX = 0.89 sec and pre-decay density that is a fraction ξ = (ρX /ρ rad|10 MeV = 0.0026 of the radiation energy density at 10 MeV; similarly good fits are found over a range where ξτX 1/2 is constant. On the other hand, decaying matter often spoils the BBN+CMB concordance, and we present limits in the (τX ,ξ) plane for both electromagnetic and invisible decays. For dark (invisible) decays, standard BBN (i.e. ξ = 0) supplies the best fit. We end with a brief discussion of the impact of future measurements including CMB-S4.

Earliest Galaxies in the JADES Origins Field: Luminosity Function and Cosmic Star Formation Rate Density 300 Myr after the Big Bang

We characterize the earliest galaxy population in the JADES Origins Field, the deepest imaging field observed with JWST. We make use of ancillary Hubble Space Telescope optical images (five filters spanning 0.4–0.9 μm) and novel JWST images with 14 filters spanning 0.8−5 μm, including seven medium-band filters, and reaching total exposure times of up to 46 hr per filter. We combine all our data at >2.3 μm to construct an ultradeep image, reaching as deep as ≈31.4 AB mag in the stack and 30.3–31.0 AB mag (5σ, r = 0.″1 circular aperture) in individual filters. We measure photometric redshifts and use robust selection criteria to identify a sample of eight galaxy candidates at redshifts z = 11.5−15. These objects show compact half-light radii of R 1/2 ∼ 50−200 pc, stellar masses of M ⋆ ∼ 107−108 M ☉, and star formation rates ∼ 0.1−1 M ☉ yr−1. Our search finds no candidates at 15 < z < 20, placing upper limits at these redshifts. We develop a forward-modeling approach to infer the properties of the evolving luminosity function without binning in redshift or luminosity that marginalizes over the photometric redshift uncertainty of our candidate galaxies and incorporates the impact of nondetections. We find a z = 12 luminosity function in good agreement with prior results, and that the luminosity function normalization and UV luminosity density decline by a factor of ∼2.5 from z = 12 to z = 14. We discuss the possible implications of our results in the context of theoretical models for evolution of the dark matter halo mass function.

Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Resolving the hot and ionized Universe through the Sunyaev-Zeldovich effect

An omnipresent feature of the multi-phase “cosmic web” — the large-scale filamentary backbone of the Universe — is that warm/hot (≳ 105 K) ionized gas pervades it. This gas constitutes a relevant contribution to the overall universal matter budget across multiple scales, from the several tens of Mpc-scale intergalactic filaments, to the Mpc intracluster medium (ICM), all the way down to the circumgalactic medium (CGM) surrounding individual galaxies, on scales from ∼ 1 kpc up to their respective virial radii (∼ 100 kpc). The study of the hot baryonic component of cosmic matter density represents a powerful means for constraining the intertwined evolution of galactic populations and large-scale cosmological structures, for tracing the matter assembly in the Universe and its thermal history. To this end, the Sunyaev-Zeldovich (SZ) effect provides the ideal observational tool for measurements out to the beginnings of structure formation. The SZ effect is caused by the scattering of the photons from the cosmic microwave background off the hot electrons embedded within cosmic structures, and provides a redshift-independent perspective on the thermal and kinematic properties of the warm/hot gas. Still, current and next-generation (sub)millimeter facilities have been providing only a partial view of the SZ Universe due to any combination of: limited angular resolution, spectral coverage, field of view, spatial dynamic range, sensitivity, or all of the above. In this paper, we motivate the development of a wide-field, broad-band, multi-chroic continuum instrument for the Atacama Large Aperture Submillimeter Telescope (AtLAST) by identifying the scientific drivers that will deepen our understanding of the complex thermal evolution of cosmic structures. On a technical side, this will necessarily require efficient multi-wavelength mapping of the SZ signal with an unprecedented spatial dynamic range (from arcsecond to degree scales) and we employ detailed theoretical forecasts to determine the key instrumental constraints for achieving our goals.

Luminous giants populate the dense Cosmic Web

Context. Giant radio galaxies (GRGs, giant RGs, or giants) are megaparsec-scale, jet-driven outflows from accretion disks of supermassive black holes, and represent the most extreme pathway by which galaxies can impact the Cosmic Web around them. A long-standing but unresolved question is why giants are so much larger than other radio galaxies. Aims. It has been proposed that, in addition to having higher jet powers than most RGs, giants might live in especially low-density Cosmic Web environments. In this work, we aim to test this hypothesis by pinpointing Local Universe giants and other RGs in physically principled, Bayesian large-scale structure reconstructions. Methods. More specifically, we localised a LOFAR Two-metre Sky Survey (LoTSS) DR2–dominated sample of luminous (lν(ν = 150 MHz)≥1024 W Hz−1) giants and a control sample of LoTSS DR1 RGs, both with spectroscopic redshifts up to zmax = 0.16, in the BORG SDSS Cosmic Web reconstructions. We measured the Cosmic Web density on a smoothing scale of ∼2.9 Mpc h−1 for each RG; for the control sample, we then quantified the relation between RG radio luminosity and Cosmic Web density. With the BORG SDSS tidal tensor, we also measured for each RG whether the gravitational dynamics of its Cosmic Web environment resemble those of clusters, filaments, sheets, or voids. Results. For both luminous giants and general RGs, the Cosmic Web density distribution is gamma distribution–like. Luminous giants populate large-scale environments that tend to be denser than those of general RGs. This result is corroborated by gravitational dynamics classification and a cluster catalogue crossmatching analysis. We find that the Cosmic Web density around RGs with 150 MHz radio luminosity lν is distributed as 1 + ΔRG | Lν = lν ∼ Γ(k, θ), where k = 4.8 + 0.2 · √, θ = 1.4 + 0.02 · √, and √:= log10(lν (1023 W Hz−1)−1). Conclusions. This work presents more than a thousand inferred megaparsec-scale densities around radio galaxies, which may be correct up to a factor of order unity – except in clusters of galaxies, where the densities can be more than an order of magnitude too low. We pave the way to a future in which megaparsec-scale densities around RGs are common inferred quantities, which help to better understand their dynamics, morphology, and interaction with the enveloping Cosmic Web. Our data demonstrate that luminous giants inhabit denser environments than general RGs. This shows that – at least at high jet powers – low-density environments are no prerequisite for giant growth. Using general RGs, we quantified the relation between radio luminosity at 150 MHz and Cosmic Web density on a smoothing scale of ∼2.9 Mpc h−1. This positive relation, combined with the discrepancy in radio luminosity between known giants and general RGs, reproduces the discrepancy in Cosmic Web density between known giants and general RGs. Our findings are consistent with the view that giants are regular, rather than mechanistically special, members of the radio galaxy population.

Neural style transfer of weak lensing mass maps

We propose a new generative model of projected cosmic mass density maps inferred from weak gravitational lensing observations of distant galaxies (weak lensing mass maps). We construct the model based on a neural style transfer so that it can transform Gaussian weak lensing mass maps into deeply non-Gaussian counterparts as predicted in ray-tracing lensing simulations. We develop an unpaired image-to-image translation method with Cycle-Consistent Generative Adversarial Networks (Cycle GAN), which learn efficient mapping from an input domain to a target domain. Our model is designed to enjoy important advantages; it is trainable with no need for paired simulation data, flexible to make the input domain visually meaningful, and expandable to rapidly-produce a map with a larger sky coverage than training data without additional learning. Using 10,000 lensing simulations, we find that appropriate labeling of training data based on field variance allows the model to reproduce a correct scatter in summary statistics for weak lensing mass maps. Compared with a popular log-normal model, our model improves in predicting the statistical natures of three-point correlations and local properties of rare high-density regions. We also demonstrate that our model enables us to produce a continuous map with a sky coverage of ∼166deg2 but similar non-Gaussian features to training data covering ∼12deg2 in a GPU minute. Hence, our model can be beneficial to massive productions of synthetic weak lensing mass maps, which is of great importance in future precise real-world analyses.

Gauge invariance, gauge field quantization, and renormalization in autoregularization

Autoregularization, a new divergence-free framework for calculating scattering amplitudes, uses a Lorentz-invariant scale harvested from the kinematics of a scattering process to regularize the amplitude of the process [N. Prabhu, .]. Preliminary validation studies show that autoregularization’s predictions are in good agreement with experimental data—across several scattering processes and a wide range of energy scales. Further, tree-level calculation of the vacuum energy density of the free fields in the Standard Model, using autoregularization, is shown to yield a value that is smaller than the current estimate of the cosmic critical density. In this paper, we prove that the scattering amplitudes in QED, calculated using autoregularization, are gauge invariant. Our proof, which is valid both for autoregularization and current theory, is stronger in that it shows the amplitude of every Feynman diagram is gauge invariant in contrast to previous proofs, which establish gauge invariance only for sum of amplitudes of Feynman diagrams of a process. Next, we show that—unlike in the standard quantization framework, which requires modification of both the quantization framework itself as well as the Lagrangian in order to quantize gauge fields in covariant gauge—in autoregularization the gauge field in QED can be quantized, in covariant gauge, without modifying the standard quantization procedure or the Lagrangian and without introducing the ghost field. Finally, we illustrate renormalization based on autoregularization up to 1-loop in φ4 theory. Since perturbative corrections are finite in autoregularization, the counterterms are not designed to remove divergences but to implement renormalization prescriptions at every order of perturbation. We also derive the renormalization group equation (RGE). Unlike in some regularization schemes (such as dimensional regularization), in which the physical meaning of the fictitious scale introduced by regularization is unclear, in autoregularization the scale in RGE has a transparent physical meaning—it is the Lorentz-invariant kinematic scale of the scattering process of interest. The increasing simplifications resulting within autoregularization and the agreement between its predictions and experimental data, together with the underlying thermodynamic argument, which shows that the framework is essential for a complete description of quantum fields, all converge to suggest that autoregularization provides the proper framework for the description of quantum fields. Published by the American Physical Society 2024