Showing papers by "Joel R. Primack published in 2023"

Wet Compaction to a Blue Nugget: a Critical Phase in Galaxy Evolution

Abstract: We utilize high-resolution cosmological simulations to reveal that high-redshift galaxies tend to undergo a robust ‘wet compaction’ event when near a ‘golden’ stellar mass of $\sim \!\!10^{10}\, \rm M_\odot$. This is a gaseous shrinkage to a compact star-forming phase, a ‘blue nugget’ (BN), followed by central quenching of star formation to a compact passive stellar bulge, a ‘red nugget’ (RN), and a buildup of an extended gaseous disc and ring. Such nuggets are observed at cosmic noon and seed today’s early-type galaxies. The compaction is triggered by a drastic loss of angular momentum due to, e.g., wet mergers, counter-rotating cold streams, or violent disc instability. The BN phase marks drastic transitions in the galaxy structural, compositional and kinematic properties. The transitions are from star-forming to quenched inside-out, from diffuse to compact with an extended disc-ring and a stellar envelope, from dark matter to baryon central dominance, from prolate to oblate stellar shape, from pressure to rotation support, from low to high metallicity, and from supernova to AGN feedback. The central black hole growth, first suppressed by supernova feedback when below the golden mass, is boosted by the compaction, and the black hole keeps growing once the halo is massive enough to lock in the supernova ejecta.

On the nature of disks at high redshift seen by JWST/CEERS with contrastive learning and cosmological simulations

Abstract: Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $z \geq 3$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not populate the same manifold in the representation space from contrastive learning, partly because the observed galaxies tend to be more compact and more elongated than the simulated galaxies. We also find that about half the galaxies that were visually classified as disks based on their elongated images actually populate a similar region of the representation space than spheroids, which according to the TNG50 simulation is occupied by objects with low stellar specific angular momentum and non-oblate structure. This suggests that the disk fraction at $z>3$ as evaluated by visual classification may be severely overestimated by misclassifying compact, elongated galaxies as disks. Deeper imaging and/or spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies.

The Evolving Effect of Cosmic Web Environment on Galaxy Quenching

Abstract: We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar mass log(M*/M⊙)≥8 to the nearest node (d node) and the nearest filament spine (d fil) to study the dependence of both the median specific star formation rate (〈sSFR〉) and the median gas fraction (〈f gas〉) on these distances. We find that the 〈sSFR〉 of galaxies is only dependent on the cosmic web environment at z < 2, with the dependence increasing with time. At z ≤ 0.5, 8≤log(M*/M⊙)<9 galaxies are quenched at d node ≲ 1 Mpc, and have significantly suppressed star formation at d fil ≲ 1 Mpc, trends driven mostly by satellite galaxies. At z ≤ 1, in contrast to the monotonic drop in 〈sSFR〉 of log(M*/M⊙)<10 galaxies with decreasing d node and d fil, log(M*/M⊙)≥10 galaxies—both centrals and satellites—experience an upturn in 〈sSFR〉 at d node ≲ 0.2 Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in 〈f gas〉. Our results suggest that in the past ∼10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments. At earlier times, cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as the Sloan Digital Sky Survey and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST, and Roman, are required to test our predictions against observations.

A new derivation of the Hubble constant from $\gamma$-ray attenuation using improved optical depths for the Fermi and CTA era

Abstract: We present $\gamma$-ray optical-depth calculations from a recently published extragalactic background light (EBL) model built from multiwavelength galaxy data from the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (HST/CANDELS). CANDELS gathers one of the deepest and most complete observations of stellar and dust emissions in galaxies. This model resulted in a robust derivation of the evolving EBL spectral energy distribution up to $z\sim 6$, including the far-infrared peak. Therefore, the optical depths derived from this model will be useful for determining the attenuation of $\gamma$-ray photons coming from high-redshift sources, such as those detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and for multi-TeV photons that will be detected from nearby sources by the future Cherenkov Telescope Array. From these newly calculated optical depths, we derive the cosmic $\gamma$-ray horizon and also measure the expansion rate and matter content of the Universe including an assessment of the impact of the EBL uncertainties. We find $H_{0}=61.9$ $^{+2.9}_{-2.4}$ km s$^{-1}$ Mpc$^{-1}$ when fixing $\Omega_{m}=0.32$, and $H_{0}=65.6$ $^{+5.6}_{-5.0}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m}=0.19\pm 0.07$, when exploring these two parameters simultaneously.