Search

Motivated by the early excess of bright galaxies seen by \textit{JWST}, we run zoom-in cosmological simulations of a massive galaxy at Cosmic Dawn (MDG), in a halo of $10^{11} M_\odot$ at $z = 9$, using the hydro-gravitational code RAMSES at an effective resolution $\sim 10~{\rm pc}$. We investigate physical mechanisms that enhance the star-formation efficiencies (SFEs) under the unique conditions of high gas density ($\sim 3\times 10^3~{\rm cm^{-3}}$, $\sim 10^4~M_\odot/{\rm pc^2}$). Our fiducial star formation recipe uses a physically-motivated, turbulence-based, multi-freefall model, avoiding ad-hoc extrapolation from lower redshifts. By $z = 9$, our simulated galaxy is a clumpy, thick, rotating disc with a high stellar mass of a few $10^9~M_\odot$ and high star formation rate of $\sim 100~M_\odot/{\rm yr}$. The high gas density makes supernova (SN) feedback less effective at suppressing star formation, producing a relatively high local SFE $\gtrsim 10\%$. The global SFE is dominated by feedback-driven outflows and is only weakly correlated with the local SFE. Photoionization heating can enhance the effects of SN feedback on the local SFE by making more SNe explode in diffuse environments, but the global SFE remains high even in our simulations with the strongest feedback. Intense accretion at Cosmic Dawn seeds strong turbulence, which reduces the local SFE for the same gas conditions due to turbulent pressure support. This prevents star-forming clouds from catastrophically collapsing, but this only weakly affects the global SFE. The star formation histories of our simulated galaxies are in the ballpark of the MDGs seen by \textit{JWST}, despite our limited resolution. They set the stage for future simulations which treat radiation self-consistently and use a higher effective resolution $\sim 1~{\rm pc}$ which captures the physics of star-forming clouds., Comment: 30 pages, 19 figures. Submitted to MNRAS on October 27th, 2024. Comments welcome. For associated movies, see https://www.youtube.com/playlist?list=PL7YbfRC6zxzAYgFEr5oefYb5dcv0dl7Ba. For associated data, see https://doi.org/10.34770/v56h-ps15

Theory and observations reveal that the circumgalactic medium (CGM) and the cosmic web at high redshifts are multiphase, with small clouds of cold gas embedded in a hot, diffuse medium. A proposed mechanism is `shattering' of large, thermally unstable clouds into tiny cloudlets of size lshatter~min(cs*tcool). We study these processes using idealized numerical simulations of thermally unstable gas clouds. We expand upon previous works by exploring the effects of cloud geometry (spheres, streams, and sheets), metallicity, and the inclusion of an ionizing UV background. We find that `shattering' is triggered by clouds losing sonic contact and rapidly imploding, leading to a reflected shock which causes the cloud to re-expand and induces Richtmyer-Meshkov instabilities at its interface. After fragmentation the cloudlets experience a drag force from the surrounding hot gas, leading to recoagulation into larger clouds. We distinguish between `fast' and `slow' coagulation regimes, showing that sheets are always in the `fast' coagulation regime while streams and spheres have a maximum overdensity for rapid coagulation. The critical overdensity for spheres is smaller than for streams, such that the coagulation efficiency increases from spheres to streams to sheets. Surprisingly, lshatter does not appear to be a characteristic clump size even if it is well resolved. Rather, fragmentation continues until the grid scale with a mass distribution of N(>m)~m^{-1}. We apply our results to the case of cold streams feeding massive (Mv>10^{12}Msun) high-z (z>2) galaxies from the cosmic web, finding that streams are likely to shatter upon entering the CGM through the virial shock. This could explain the large clumping factors and covering fractions of cold gas in the CGM around such galaxies, and may be related to galaxy quenching by preventing cold streams from reaching the central galaxy. [abridged], Comment: 36 pages, 20 figures, submitted to MNRAS. Comments are welcome!

The scenario of feedback-free starbursts (FFB), which predicts excessively bright galaxies at cosmic dawn as observed using JWST, may provide a natural setting for black hole (BH) growth. This involves the formation of intermediate-mass seed BHs and their runaway mergers into super-massive BHs with high BH-to-stellar mass ratios and low AGN luminosities. We present a scenario of merger-driven BH growth in FFB galaxies and study its feasibility. BH seeds form within the building blocks of the FFB galaxies, namely, thousands of compact star clusters, each starbursting in a free-fall time of a few Myr before the onset of stellar and supernova feedback. The BH seeds form by rapid core collapse in the FFB clusters, in a few free-fall times, sped up by the migration of massive stars due to the young, broad stellar mass function and stimulated by a `gravo-gyro' instability due to internal cluster rotation and flattening. BHs of $10^4 M_\odot$ are expected in $10^6 M_\odot$ FFB clusters within sub-kpc galactic disks at $z \sim 10$. The BHs then migrate to the galaxy center by dynamical friction, hastened by the compact FFB stellar galactic disk configuration. Efficient mergers of the BH seeds will produce $10^{6-8} M_\odot$ BHs with a BH-to-stellar mass ratio $\sim 0.01$ by $z \sim 4-7$, as observed. The growth of the central BH by mergers can overcome the bottleneck introduced by gravitational wave recoils if the BHs inspiral within a relatively cold disk or if the escape velocity from the galaxy is boosted by a wet compaction event. Such events, common in massive galaxies at high redshifts, can also help by speeding up the inward BH migration and by providing central gas to assist with the final parsec problem. The cold disk version of the FFB scenario provides a feasible route for the formation of supermassive BHs., Comment: 24 pages, 14 figures

We study the radial transport of cold gas within simulated disk galaxies at cosmic noon, aiming at distinguishing between disk instability and accretion along cold streams from the cosmic web as its driving mechanism. Disks are selected based on kinematics and flattening from the VELA zoom-in hydro-cosmological simulations. The radial velocity fields in the disks are mapped, their averages are computed as a function of radius and over the whole disk, and the radial mass flux in each disk as a function of radius is obtained. The transport directly associated with fresh incoming streams is identified by selecting cold gas cells that are either on incoming streamlines or have low metallicity. The radial velocity fields in VELA disks are found to be highly non-axisymmetric, showing both inflows and outflows. However, in most cases, the average radial velocities, both as a function of radius and over the whole disk, are directed inwards, with the disk-averaged radial velocities typically amounting to a few percent of the disk-averaged rotational velocities. This is significantly lower than the expectations from various models that analytically predict the inward mass transport as driven by torques associated with disk instability. Under certain simplifying assumptions, the latter typically predict average inflows of more than $10\%$ of the rotational velocities. Analyzing the radial motions of streams and off-stream material, we find that the radial inflow in VELA disks is dominated by the stream inflows themselves, especially in the outer disks. The high inward radial velocities inferred in observed disks at cosmic noon, at the level of $\sim \! 20\%$ of the rotational velocities, may reflect inflowing streams from the cosmic web rather than being generated by disk instability., Comment: 20 pages, 14 figures, submitted to A&A

We address the formation of giant clumps in violently unstable gas-rich disc galaxies at cosmic noon. While these are commonly thought to originate from gravitational Toomre instability, cosmological simulations have indicated that clumps form even in regions where the Toomre $Q$ parameter is well above unity, which should be stable according to linear Toomre theory (Inoue et al., 2016). Examining one of these cosmological simulations, we find that it exhibits an excess in compressive modes of turbulence with converging motions. The energy in converging motions within proto-clump regions is $\sim 70\%$ of the total turbulent energy, compared to $\sim 17\%$ expected in equipartition. When averaged over the whole disc, $\sim 32\%$ of the turbulent energy is in converging motions, with a further $\sim 8\%$ in diverging motions. Thus, a total of $\sim 40\%$ of the turbulent energy is in compressive modes, with the rest in solenoidal modes, compared to the $(1/3)-(2/3)$ division expected in equipartition. By contrast, we find that in an isolated-disc simulation with similar properties, resembling high-$z$ star-forming galaxies, the energy in the different turbulence modes are in equipartition, both in proto-clump regions and over the whole disc. We conclude that the origin of the excessive converging motions in proto-clump regions is external to the disc, and propose several mechanisms that can induce them. This is an additional mechanism for clump formation, complementary to and possibly preceding gravitational instability., Comment: 6 pages, 4 figures, Submitted to MNRAS Letters Comments welcome!

The observed rest-UV luminosity function at cosmic dawn ($z \sim 8-14$) measured by JWST revealed an excess of UV-luminous galaxies relative to many pre-launch theoretical predictions. A high star-formation efficiency (SFE) and a top-heavy initial mass function (IMF) are among the mechanisms proposed for explaining this excess. Although a top-heavy IMF has been proposed for its ability to increase the light-to-mass ratio (\(\Psi_{\mathrm{UV}}\)), the resulting enhanced radiative pressure from young stars could decrease the star formation efficiency (SFE), potentially driving galaxy luminosities back down. In this Letter, we use idealized radiation hydrodynamic simulations of star cluster formation to explore the effects of a top-heavy IMF on the SFE of clouds typical of the high pressure conditions found at these redshifts. We find that the SFE in star clusters with solar neighbourhood-like dust abundance decreases with increasingly top-heavy IMF's -- by $\sim 20 \%$ for an increase of factor 4 in $\Psi_{\mathrm{UV}}$, and by $50 \%$ for a factor $ \sim 10$ in $\Psi_{\mathrm{UV}}$. However, we find that an expected decrease in the dust-to-gas ratio ($\sim 0.01 \times \mathrm{Solar}$) at these redshifts can completely compensate for the enhanced light output. This leads to a (cloud-scale; $\sim 10 \, \mathrm{pc}$) SFE that is $\gtrsim 70\%$ even for a factor 10 increase in $\Psi_{\mathrm{UV}}$, implying that highly efficient star formation is unavoidable for high surface density and low metallicity conditions. Our results suggest that a top-heavy IMF, if present, likely coexists with efficient star formation in these galaxies., Comment: 12 pages, 4 figures. Accepted to the Astrophysical Journal Letters

Nearby dwarf galaxies display a variety of effective radii (sizes) at a given stellar mass, suggesting different evolution scenarios according to their final "stellar" size. The TNG hydrodynamical simulations present a bimodality in the z = 0 size - mass relation (SMRz0) of dwarf galaxies, at $r_{1/2,\star}$ ~ 450 pc. Using the TNG50 simulation, we explored the evolution of the most massive progenitors of dwarf galaxies (z=0 $\log( M_\star / \mathrm{M}_\odot)$ between 8.4 and 9.2) that end up as central galaxies of their groups. We split these dwarfs into three classes of the SMRz0: "Normals" from the central spine of the main branch, and "Compacts" from the secondary branch as well as the lower envelope of the main branch. Both classes of Compacts see their stellar sizes decrease from z ~ 1 onwards in contrast to Normals, while the sizes of the gas and dark matter (DM) components continue to increase (as for Normals). A detailed analysis reveals that Compacts live in poorer environments, and thus suffer fewer major mergers from z = 0.8 onwards, which otherwise would pump angular momentum into the gas, allowing strong gas inflows, producing inner star formation, and thus leading to the buildup of a stellar core. Compacts are predicted to be rounder and to have bluer cores. Compact dwarfs of similar sizes are observed in the GAMA survey, but the bimodality in size is less evident and the most compact dwarfs tend to be passive rather than star forming, as in TNG50. Our conclusions should therefore be confirmed with future cosmological hydrodynamical simulations., Comment: 15 pages, 15 figures, 3 tables, Accepted for Publication in A&A. Final version

We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their locations in the cosmic web, internal structures, and assembly histories. Our analysis of the TNG50 simulation reveals that the orbital-circularity threshold for disc differentiation varies across galaxies, with systematic trends in mass and redshift, so the widely used decomposition method with constant circularity cuts is oversimplified and underestimates thin disc at JWST redshifts. The energy threshold between the stellar halo and the inner galaxy is also a function of mass and redshift, minimizing at the sub-Galactic halo mass, where the circularity threshold peaks. Revisiting the issue of galaxy size predictor, we show that disc sizes in TNG50 exhibit correlations with three structural parameters besides virial mass and redshift: 1) a positive correlation with halo spin $\lambda$ across redshifts -- stronger than previously reported for zoom-in simulations but still weaker than the simple $r_{1/2}/R_{\rm vir} \propto \lambda$ scaling; 2) an anti-correlation with DM concentration $c$ that is well described by $r_{1/2}/R_{\rm vir} \propto c^{-0.7}$ even when $c$ is measured in the DM only run; 3) more actively accreting haloes having slightly larger discs, as well as more significant stellar haloes and lower thin-to-thick ratio. Disc mass fraction is higher in rounder haloes and in cosmic knots and filaments, implying that disc development needs both stable halo conditions and continuous material supply. Our methodology is public and adaptable to other simulations., Comment: 20 pages, 17 figures, submitted to MNRAS

One of the main challenges in galaxy formation that has emerged recently is the early assembly of massive galaxies. The observed number density and the maximum stellar mass ($M_{\star}$) of massive galaxies in the early Universe appear to be higher than model predictions, which may pose a serious problem to the LCDM cosmology. A major limitation in many previous studies is the large uncertainty in estimating $M_{\star}$ due to the lack of constraints in the rest-frame near-infrared part of the spectral energy distribution, which is critical to determining $M_{\star}$ accurately. Here we use data from a large JWST/MIRI survey in the PRIMER program to carry out a systematic analysis of massive galaxies at $z \sim 3-8$, leveraging photometric constraints at rest-frame $\gtrsim 1 \mu$m. We find a significant reduction in the number and mass densities of massive galaxies at $z > 5$ compared to earlier results that did not use the MIRI photometry. Within the standard $\Lambda$CDM cosmology, our results require a moderate increase in the baryon-to-star conversion efficiency ($\epsilon$) towards higher redshifts and higher $M_{\star}$. For the most massive galaxies at $z\sim 8$, the required $\epsilon$ is $\sim 0.3$, in comparison to $\epsilon \sim 0.14$ for typical low-redshift galaxies. Our findings are consistent with models assuming suppressed stellar feedback due to the high gas density and the associated short free-fall time expected for massive halos at high redshift., Comment: 23 pages, 10 figures, submitted

We explore the evolution of cold streams from the cosmic web that feed galaxies through their shock-heated circumgalactic medium (CGM) at cosmic noon, $z\simeq 1-5$. In addition to the hydrodynamical instabilities and radiative cooling that we have incorporated in earlier works, we embed the stream and the hot CGM in the gravitational potential of the host dark-matter halo, deriving equilibrium profiles for both. Self-gravity within the stream is tentatively ignored. We find that the cold streams gradually entrain a large mass of initially hot CGM gas that cools in the mixing layer and condenses onto the stream. This entrainment, combined with the acceleration down the gravitational potential well, typically triples the inward cold inflow rate into the central galaxy, compared to the original rate at the virial radius, which makes the entrained gas the dominant source of gas supply to the galaxy. The potential sources for the hot gas to be entrained are recycled enriched gas that has been previously ejected from the galaxy, and fresh virial-shock-heated gas that has accumulated in the CGM. This can naturally elevate the star formation rate in the galaxy by a factor of $\sim 3$ compared to the gas accretion rate onto the halo, thus explaining the otherwise puzzling observed excess of star formation at cosmic noon. When accounting for self-shielding of dense gas from the UV background, we find that the energy radiated from the streams, originating predominantly from the cooling of the entrained gas, is consistent with observed Lyman-$\alpha$ blobs around galaxies., Comment: 27 pages, 11 figures, accepted in MNRAS

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available., Comment: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org ) for more information

We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms., Comment: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature., Comment: 32 pages, 17 figures. Published in ApJ

Ultraviolet (UV; rest-frame $\sim1200-2000$ A) spectra provide a wealth of diagnostics to characterize fundamental galaxy properties, such as their chemical enrichment, the nature of their stellar populations, and their amount of Lyman-continuum (LyC) radiation. In this work, we leverage publicly released JWST data to construct the rest-frame UV-to-optical composite spectrum of a sample of 63 galaxies at $5.6

We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $\sigma(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $\sigma(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$., Comment: 28 pages, 12 figures, 2 Appendix figures

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 45000 AA. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: all spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the SED-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hb + [OIII] and Ha in the neighboring broad-band filters. We present examples of EELGs that could be incorrectly classified at ultra-high redshift (z>12) as a result of extreme Hb + [OIII] emission blended across the reddest photometric filters. We compare the EELGs to the broader (sub-extreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high equivalent width tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. Photometrically inferred emission-line ratios are consistent with ISM conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies., Comment: 24 pages, 25 figures, submitted to ApJ

We extend the analysis of a physical model within the standard cosmology that robustly predicts a high star-formation efficiency (SFE) in massive galaxies at cosmic dawn due to feedback-free starbursts (FFBs). It implies an excess of bright galaxies at z>~10 compared to the standard models based on the low SFE at later epochs, an excess indicated by JWST observations. Here we provide observable predictions based on the analytic FFB scenario. These can be compared with simulations and JWST observations. We approximate the SFE as a function of redshift and mass, assuming a maximum SFE of 0.2~1 in the FFB regime. From this, we derive the evolution of the galaxy mass and luminosity functions as well as the evolution of stellar and star-formation densities. We then predict the star-formation history (SFH), galaxy sizes, outflows, gas fractions, metallicities, and dust attenuation, all as functions of mass and redshift in the FFB regime. The major distinguishing feature is the occurrence of FFBs above a mass threshold that declines with redshift. The luminosities and star formation rates in bright galaxies are predicted to be in excess of extrapolations of standard empirical models and cosmological simulations, an excess that grows from z~9 to higher redshifts. The FFB phase of ~100 Myr is predicted to show a characteristic SFH that fluctuates on a timescale of ~10 Myr. The stellar systems are compact (Re~0.3 kpc at z~10 and declining with z). The galactic gas consists of a steady wind driven by supernovae from earlier generations, with high outflow velocities (FWHM~1400-6700km/s), low gas fractions (<0.1), low metallicities (<~0.1 solar), and low dust attenuation ($A_{UV}$~0.5 at z~10 and declining with z). We make tentative comparisons with current JWST observations for initial insights, anticipating more complete and reliable datasets for detailed quantitative comparisons in the future., Comment: 19 pages, 13 figures; accepted to A&A; code and tables avaible at https://github.com/syrte/ffb_predict

The wavelength-coverage and sensitivity of JWST now enables us to probe the rest-frame UV - optical spectral energy distributions (SEDs) of galaxies at high-redshift ($z>4$). From these SEDs it is, in principle, through SED fitting possible to infer key physical properties, including stellar masses, star formation rates, and dust attenuation. These in turn can be compared with the predictions of galaxy formation simulations allowing us to validate and refine the incorporated physics. However, the inference of physical properties, particularly from photometry alone, can lead to large uncertainties and potential biases. Instead, it is now possible, and common, for simulations to be \emph{forward-modelled} to yield synthetic observations that can be compared directly to real observations. In this work, we measure the JWST broadband fluxes and colours of a robust sample of $58$ the distributions differ somewhat, though our observed sample size is small and thus susceptible to statistical fluctuations. Likewise, the predicted and observed colour evolution show broad agreement, at least at $58$, though, again, the sample size is small here., Comment: 11 pages, 10 figures, submitted to MNRAS

The majority of massive disk galaxies in the local Universe show a stellar barred structure in their central regions, including our Milky Way. Bars are supposed to develop in dynamically cold stellar disks at low redshift, as the strong gas turbulence typical of disk galaxies at high redshift suppresses or delays bar formation. Moreover, simulations predict bars to be almost absent beyond $z = 1.5$ in the progenitors of Milky Way-like galaxies. Here we report observations of ceers-2112, a barred spiral galaxy at redshift $z_{\rm phot} \sim 3$, which was already mature when the Universe was only 2 Gyr old. The stellar mass ($M_{\star} = 3.9 \times 10^9 M_{\odot}$) and barred morphology mean that ceers-2112 can be considered a progenitor of the Milky Way, in terms of both structure and mass-assembly history in the first 2 Gyr of the Universe, and was the closest in mass in the first 4 Gyr. We infer that baryons in galaxies could have already dominated over dark matter at $z \sim 3$, that high-redshift bars could form in approximately 400 Myr and that dynamically cold stellar disks could have been in place by redshift $z = 4-5$ (more than 12 Gyrs ago)., Comment: Published in Nature. Accepted on September 12

We present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z~9, 11 and 14, and show that the space density of bright (M_UV=-20) galaxies changes only modestly from z~14 to z~9, compared to a steeper increase from z~8 to z~4. While our candidates are photometrically selected, spectroscopic followup has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z>10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models which naturally have enhanced star-formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic followup of a large sample of early galaxies can distinguish between these competing scenarios., Comment: Submitted to ApJL. Main paper is 33 pages, 13 figures, 5 tables. Two appendices with additional figures and tables

The 3D geometry of high-redshift galaxies remains poorly understood. We build a differentiable Bayesian model and use Hamiltonian Monte Carlo to efficiently and robustly infer the 3D shapes of star-forming galaxies in JWST-CEERS observations with $\log M_*/M_{\odot}=9.0-10.5$ at $z=0.5-8.0$. We reproduce previous results from HST-CANDELS in a fraction of the computing time and constrain the mean ellipticity, triaxiality, size and covariances with samples as small as $\sim50$ galaxies. We find high 3D ellipticities for all mass-redshift bins suggesting oblate (disky) or prolate (elongated) geometries. We break that degeneracy by constraining the mean triaxiality to be $\sim1$ for $\log M_*/M_{\odot}=9.0-9.5$ dwarfs at $z>1$ (favoring the prolate scenario), with significantly lower triaxialities for higher masses and lower redshifts indicating the emergence of disks. The prolate population traces out a ``banana'' in the projected $b/a-\log a$ diagram with an excess of low $b/a$, large $\log a$ galaxies. The dwarf prolate fraction rises from $\sim25\%$ at $z=0.5-1.0$ to $\sim50-80\%$ at $z=3-8$. If these are disks, they cannot be axisymmetric but instead must be unusually oval (triaxial) unlike local circular disks. We simultaneously constrain the 3D size-mass relation and its dependence on 3D geometry. High-probability prolate and oblate candidates show remarkably similar S\'ersic indices ($n\sim1$), non-parametric morphological properties and specific star formation rates. Both tend to be visually classified as disks or irregular but edge-on oblate candidates show more dust attenuation. We discuss selection effects, follow-up prospects and theoretical implications., Comment: Accepted version to appear in ApJ, main body is 36 pages of which ~half are full-page figures

Different star-formation models at Cosmic Dawn produce detectable signatures in the observables of upcoming 21-cm experiments. In this work, we consider the physical scenario of feedback-free starbursts (FFB), according to which the star-formation efficiency (SFE) is enhanced in sufficiently massive halos at early enough times, thus explaining the indication from the James Webb Space Telescope for an excess of bright galaxies at $z \geq 10$. We model the contribution of FFBs to popII SFE and compute the impact these have on the 21-cm global signal and power spectrum. We show that FFBs affect the evolution of the brightness temperature and the 21-cm power spectrum, but they only have a limited effect on the neutral hydrogen fraction. We investigate how the observables are affected by changes in the underlying star formation model and by contribution from popIII stars. Finally, we forecast the capability of next-generation Hydrogen Epoch of Reionization Array (HERA) to detect the existence of FFB galaxies via power spectrum measurements. Our results show the possibility of a significant detection, provided that popII stars are the main drivers of lowering the spin temperature. Efficient popIII star formation will make the detection more challenging., Comment: Version accepted for publication on MNRAS

Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 {\AA} luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar S\'ersic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center., Comment: 21 pages, 13 figures, accepted for publication in ApJ

We investigate the morphology and resolved physical properties of a sample of 22 IR-selected DSFG at cosmic noon using the JWST/NIRCam images obtained in the EGS field for the CEERS survey. The resolution of the NIRCam images allowed to spatially resolve these galaxies up to 4.4um and identify their bulge even when extinguished by dust. The goal of this study is to obtain a better understanding of the formation and evolution of FIR-bright galaxies by spatially resolving their properties using JWST in order to look through the dust and bridge the gap between the compact FIR sources and the larger optical SFG. Based on RGB images from the NIRCam filters, we divided each galaxy into several uniformly colored regions, fitted their respective SEDs, and measured physical properties. After classifying each region as SF or quiescent, we assigned galaxies to three classes, depending on whether active SF is located in the core, in the disk or in both. We find (i) that galaxies at a higher z tend to have a fragmented disk with a low core mass fraction. They are at an early stage of bulge formation. When moving toward a lower z, the core mass fraction increases, and the bulge growth is associated with a stabilization of the disk: the NIRCam data clearly point toward bulge formation in preexisting disks. (ii) Lopsidedness is a common feature of DSFGs. It could have a major impact on their evolution; (iii) 23% of galaxies have a SF core embedded in a quiescent disk. They seem to be undergoing outside-in quenching, often facilitated by their strong lopsidedness inducing instabilities. (iv) We show that half of our galaxies with SF concentrated in their core are good SMG counterpart candidates, demonstrating that compact SMGs are usually surrounded by a larger, less obscured disk. (v) Finally, we found surprising evidence for clump-like substructures being quiescent or residing in quiescent regions., Comment: 23 pages, 24 figures, Accepted by A&A

We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \,{\rm kpc}$. The post-shock gas is in virial equilibrium with the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\,\text{M}_\odot\,{\rm kpc}^{-1}$, the gas fraction within $r_{\rm shock}$ is the universal baryon fraction, and the virial temperature is $\sim 7\times 10^5 {\rm K}$. In the outer ''thermal'' (T) zone, $r \geq 0.65 \, r_{\rm shock}$, inward gravity and ram-pressure forces are over-balanced by outwards thermal pressure forces, decelerating the inflowing gas expanding the shock outward. In the intermediate ''vortex'' (V) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outwards force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. The shear and turbulent forces associated with the vortices act inward. The inner ''stream'' (S) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \,{\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outwards forces, though it exhibits both inflows and outflows. [abridged], Comment: 32 pages, 22 figures, submitted to MNRAS

The recent reported discovery of a low-mass $z=5.2$ and an intermediate-mass $z=7.3$ quenched galaxy with JWST/NIRSpec is the first evidence of halted star formation above $z\approx 5$. Here we show how bursty star formation at $z=4-8$ gives rise to temporarily quenched, or mini-quenched galaxies in the mass range $M_{\star} = 10^7-10^9 \ M_{\odot}$ using four models of galaxy formation: the periodic box simulation IllustrisTNG, the zoom-in simulations VELA and FirstLight and an empirical halo model. The main causes for mini-quenching are stellar feedback, lack of gas accretion onto galaxies and galaxy-galaxy interactions. The abundance of (mini-)quenched galaxies agrees across the models: the population first appears below $z\approx 8$, after which their proportion increases with cosmic time, from $\sim 0.5-1.0$% at $z=7$ to $\sim 2-4$% at $z=4$, corresponding to comoving number densities of $\sim 10^{-5}$ Mpc$^{-3}$ and $\sim 10^{-3}$ Mpc$^{-3}$, respectively. These numbers are consistent with star formation rate duty cycles inferred for VELA and FirstLight galaxies. Their star formation histories (SFHs) suggest that mini-quenching at $z=4-8$ is short-lived with a duration of $\sim 20-40$ Myr, which is close to the free-fall timescale of the inner halo. However, mock spectral energy distributions of mini-quenched galaxies in IllustrisTNG and VELA do not match JADES-GS-z7-01-QU photometry, unless their SFHs are artificially altered to be more bursty on timescales of $\sim 40$ Myr. Studying mini-quenched galaxies might aid in calibrating sub-grid models governing galaxy formation, as these may not generate sufficient burstiness at high redshift to explain the SFH inferred for JADES-GS-z7-01-QU., Comment: 13 pages, 5 figures, 5 tables

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

We present JWST/NIRSpec prism spectroscopy of seven galaxies selected from the Cosmic Evolution Early Release Science Survey (CEERS) NIRCam imaging with photometric redshifts z_phot>8. We measure emission line redshifts of z=7.65 and 8.64 for two galaxies, and z=9.77(+0.37,-0.29) and 10.01(+0.14,-0.19) for two others via the detection of continuum breaks consistent with Lyman-alpha opacity from a mostly neutral intergalactic medium. The presence (absense) of strong breaks (strong emission lines) give high confidence that these two galaxies are at z>9.6, but the break-derived redshifts have large uncertainties given the low spectral resolution and relatively low signal-to-noise of the CEERS NIRSpec prism data. The two z~10 sources are relatively luminous (M_UV<-20), with blue continua (-2.38 candidates with CEERS NIRSpec spectroscopy do not have secure redshifts, but the absence of emission lines in their spectra is consistent with redshifts z>9.6. We find that z>8 photometric redshifts are generally in agreement (within uncertainties) with the spectroscopic values. However, the photometric redshifts tend to be slightly overestimated (average Delta(z)=0.50+/-0.12), suggesting that current templates do not fully describe the spectra of very high-z sources. Overall, our results solidifies photometric evidence for a high space density of bright galaxies at z>8 compared to theoretical model predictions, and further disfavors an accelerated decline in the integrated UV luminosity density at z>8., Comment: Published in ApJL. 19 pages, 9 figures, 7 tables. File with Table 6 included in source .tar file

The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic winds in driving radial metallicity profiles and setting the mass-metallicity gradient relation (MZGR). We include a comprehensive treatment of these processes by including them as sources that supply mass, metals, and energy to marginally unstable galactic discs in pressure and energy balance. We show that both feedback and accretion that can drive turbulence and enhance metal-mixing via diffusion are crucial to reproduce the observed MZGR in local galaxies. Metal transport also contributes to setting metallicity profiles, but it is sensitive to the strength of radial gas flows in galaxies. While the mass loading of galactic winds is important to reproduce the mass metallicity relation (MZR), we find that metal mass loading is more important to reproducing the MZGR. Specifically, our model predicts preferential metal enrichment of galactic winds in low-mass galaxies. This conclusion is robust against our adopted scaling of the wind mass-loading factor, uncertainties in measured wind metallicities, and systematics due to metallicity calibrations. Overall, we find that at $z \sim 0$, galactic winds and metal transport are more important in setting metallicity gradients in low-mass galaxies whereas star formation feedback and gas accretion dominate setting metallicity gradients in massive galaxies., Comment: 21 pages (+ appendix), 11 figures, 4 tables. Accepted by MNRAS

JWST observations indicate a surprising excess of luminous galaxies at $z\sim 10$ and above, consistent with efficient conversion of the accreted gas into stars, unlike the suppression of star formation by feedback at later times. We show that the high densities and low metallicities at this epoch {\it guarantee} a high star-formation efficiency (SFE) in the most massive dark-matter haloes. Feedback-free starbursts (FFBs) occur when the free-fall time is shorter than $\sim 1$ Myr, below the time for low-metallicity massive stars to develop winds and supernovae. This corresponds to a characteristic density of $\sim 3\times 10^3$cm$^{-3}$. A comparable threshold density permits a starburst by allowing cooling to star-forming temperatures in a free-fall time. The galaxies within $\sim 10^{11} M_\odot$ haloes at $z \sim 10$ are expected to have FFB densities. The halo masses allow efficient gas supply by cold streams in a halo crossing time $\sim 80$ Myr. The FFBs gradually turn all the accreted gas into stars in clusters of $\sim 10^{4-7} M_\odot$ within galaxies that are rotating discs or shells. The starbursting clouds are insensitive to radiative feedback and are shielded against feedback from earlier stars. We predict high SFE above thresholds in redshift and halo mass, where the density is $10^{3-4}$cm$^{-3}$. The $z\sim 10$ haloes of $\sim 10^{10.8} M_\odot$ are predicted to host galaxies of $\sim 10^{10} M_\odot$ with SFR $\sim 65 M_\odot$ yr$^{-1}$ and sub-kpc sizes. The metallicity is $\leq 0.1 Z_\odot$ with little gas, dust, outflows and hot circumgalactic gas, allowing a top-heavy IMF but not requiring it. The compact galaxies with thousands of young FFB clusters may have implications on reionization, black-hole growth and globular clusters., Comment: 20 pages, 7 figures

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} 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., Comment: 33 pages, 26 figures in the main body (49 pages, 45 figures including appendix). Corresponds to accepted MNRAS version, minor revisions

We study the evolution of $L_{*}$ elliptical galaxies in the color-magnitude diagram in terms of their star-formation history and environment, in an attempt to learn about their quenching process. We have visually extracted 1109 $L_{*}$ galaxies from a sample of 36500 galaxies that were spectroscopically selected from Stripe82 of the Sloan Digital Sky Survey. From this sample we have selected 51 ellipticals based on their surface-brightness profile being well-fitted by a single S$\acute{e}$rsic profile with S$\acute{e}$rsic indices $3

We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z~0.5 to z>10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 sq. arcmin, to search for candidate galaxies at z>9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multi-band photometry across seven NIRCam broad and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z~9-16. These objects are compact with a median half-light radius of ~0.5 kpc. We present an early estimate of the z~11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M_UV ~ -20 appears to evolve very little from z~9 to z~11. We also find that the abundance (surface density [arcmin^-2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z>10 star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star-formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultra-high-redshift galaxies with which to further explore these conclusions., Comment: Replaced with published version

The new capabilities that JWST offers in the near- and mid-infrared (IR) are used to investigate in unprecedented detail the nature of optical/near-IR faint, mid-IR bright sources, HST-dark galaxies among them. We gather JWST data from the CEERS survey in the EGS, jointly with HST data, and analyze spatially resolved optical-to-mid-IR spectral energy distributions (SEDs) to estimate both photometric redshifts in 2 dimensions and stellar populations properties in a pixel-by-pixel basis. We select 138 galaxies with F150W-F356W>1.5 mag, F356W<27.5 mag. The nature of these sources is threefold: (1) 71% are dusty star-forming galaxies at 2100 Gyr^-1); (2) 18% are quiescent/dormant (i.e., subject to reignition and rejuvenation) galaxies at 3

Shell galaxies make a class of tidally distorted galaxies, characterised by wide concentric arc(s), extending out to large galactocentric distances with sharp outer edges. Recent observations of young massive star clusters in the prominent outer shell of NGC 474 suggest that such systems host extreme conditions of star formation. In this paper, we present a hydrodynamic simulation of a galaxy merger and its transformation into a shell galaxy. We analyse how the star formation activity evolves with time, location-wise within the system, and what are the physical conditions for star formation. During the interaction, an excess of dense gas appears, triggering a starburst, i.e. an enhanced star formation rate and a reduced depletion time. Star formation coincides with regions of high molecular gas fraction, such as the galactic nucleus, spiral arms, and occasionally the tidal debris during the early stages of the merger. Tidal interactions scatter stars into a stellar spheroid, while the gas cools down and reforms a disc. The morphological transformation after coalescence stabilises the gas and thus quenches star formation, without the need for feedback from an active galactic nucleus. This evolution shows similarities with a compaction scenario for compact quenched spheroids at high-redshift, yet without a long red nugget phase. Shells appear after coalescence, during the quenched phase, implying that they do not host the conditions necessary for in situ star formation. The results suggest that shell-forming mergers might be part of the process of turning blue late-type galaxies into red and dead early-types., Comment: Accepted for publication by MNRAS, 13 pages

The evolution of star-forming galaxies at high redshifts is very sensitive to the strength and nature of stellar feedback. Using two sets of cosmological, zoom-in simulations from the VELA suite, we compare the effects of two different models of feedback: with and without kinetic feedback from the expansion of supernovae shells and stellar winds. At a fixed halo mass and redshift, the stellar mass is reduced by a factor of 1-3 in the models with stronger feedback, so the stellar-mass-halo-mass relation is in better agreement with abundance matching results. On the other hand, the three-dimensional shape of low-mass galaxies is elongated along a major axis in both models. At a fixed stellar mass, Ms<10^10 Msun, galaxies are more elongated in the strong-feedback case. More massive, star-forming discs with high surface densities form giant clumps. However, the population of round, compact, old (age_c > 300 Myr), quenched, stellar (or gas-poor) clumps is absent in the model with strong feedback. On the other hand, giant star-forming clumps with intermediate ages (age_c= 100 - 300 Myr) can survive for several disc dynamical times, independently of feedback strength. The evolution through compaction followed by quenching in the plane of central surface density and specific star-formation rate is similar under the two feedback models., Comment: 14 pages, 12 figures, accepted version at MNRAS

Stellar bars are key drivers of secular evolution in galaxies and can be effectively studied using rest-frame near-infrared (NIR) images, which trace the underlying stellar mass and are less impacted by dust and star formation than rest-frame UV or optical images. We leverage the power of {\it{JWST}} CEERS NIRCam images to present the first quantitative identification and characterization of stellar bars at $z>1$ based on rest-frame NIR F444W images of high resolution (~1.3 kpc at z ~ 1-3). We identify stellar bars in these images using quantitative criteria based on ellipse fits. For this pilot study, we present six examples of robustly identified bars at $z>1$ with spectroscopic redshifts, including the two highest redshift bars at ~2.136 and 2.312 quantitatively identified and characterized to date. The stellar bars at $z$ ~ 1.1-2.3 presented in our study have projected semi-major axes of ~2.9-4.3 kpc and projected ellipticities of ~0.41-0.53 in the rest-frame NIR. The barred host galaxies have stellar masses ~ $ 1 \times 10^{10}$ to $2 \times 10^{11}$ $M_{\odot}$, star formation rates of ~ 21-295 $M_{\odot}$ yr$^{-1}$, and several have potential nearby companions. Our finding of bars at $z$ ~1.1-2.3 demonstrates the early onset of such instabilities and supports simulations where bars form early in massive dynamically cold disks. It also suggests that if these bars at lookback times of 8-10 Gyr survive out to present epochs, bar-driven secular processes may operate over a long time and have a significant impact on some galaxies by z ~ 0., Comment: 16 pages, 5 figures. Accepted for Publication in Astrophysical Journal Letters

We present the first comprehensive release of photometric redshifts (photo-z's) from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) team. We use statistics based upon the Quantile-Quantile (Q--Q) plot to identify biases and signatures of underestimated or overestimated errors in photo-z probability density functions (PDFs) produced by six groups in the collaboration; correcting for these effects makes the resulting PDFs better match the statistical definition of a PDF. After correcting each group's PDF, we explore three methods of combining the different groups' PDFs for a given object into a consensus curve. Two of these methods are based on identifying the minimum f-divergence curve, i.e., the PDF that is closest in aggregate to the other PDFs in a set (analogous to the median of an array of numbers). We demonstrate that these techniques yield improved results using sets of spectroscopic redshifts independent of those used to optimize PDF modifications. The best photo-z PDFs and point estimates are achieved with the minimum f-divergence using the best 4 PDFs for each object (mFDa4) and the Hierarchical Bayesian (HB4) methods, respectively. The HB4 photo-z point estimates produced $\sigma_{\rm NMAD} = 0.0227/0.0189$ and $|\Delta z/(1+z)| > 0.15$ outlier fraction = 0.067/0.019 for spectroscopic and 3D-HST redshifts, respectively. Finally, we describe the structure and provide guidance for the use of the CANDELS photo-z catalogs, which are available at https://archive.stsci.edu/hlsp/candels., Comment: 35 pages, 19 figures, published in ApJ, data available at https://archive.stsci.edu/hlsp/candels

We study the survival versus disruption of the giant clumps in high-redshift disc galaxies, short-lived (S) versus long-lived (L) clumps and two L sub-types, via analytic modeling tested against simulations. We develop a criterion for clump survival, with or without their gas, based on a predictive survivability parameter $S$. It compares the energy sources by supernova feedback and gravitational contraction to the clump binding energy and losses by outflows and turbulence dissipation. The clump properties are derived from Toomre instability, approaching virial and Jeans equilibrium, and the supernova energy deposit is based on an up-to-date bubble analysis. For moderate feedback levels, we find that L clumps exist with circular velocities $\sim\!50\, km\, s^{-1}$ and masses $\geq\!10^8\, M_\odot$. They are likely in galaxies with circular velocities $\geq \!200\, km\, s^{-1}$, consistent at $z \sim 2$ with the favored stellar mass for discs, $\geq\!10^{9.3}\, M_\odot$. L clumps favor disc gas fractions $\geq\!0.3$, low-mass bulges and redshifts $z\!\sim\! 2$. The likelihood of L clumps is reduced if the feedback is more ejective, e.g., if the supernovae are optimally clustered, if radiative feedback is very strong, if the stellar initial mass function is top-heavy, or if the star-formation-rate efficiency is particularly high. A sub-type of L clumps (LS), which lose their gas in several free-fall times but retain bound stellar components, may be explained by a smaller contraction factor and stronger external gravitational effects, where clump mergers increase the SFR efficiency. The more massive L clumps (LL) retain most of their baryons for tens of free-fall times with a roughly constant star-formation rate., Comment: 26 pages, 14 figures

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We propose an analytic model, CuspCore II, for the response of dark matter (DM) haloes to central gas ejection, as a mechanism for generating DM-deficient cores in dwarfs and high-z massive galaxies. We test this model and three other methods using idealized N-body simulations. The current model is physically justified and provides more accurate predictions than the earlier version, CuspCore I (Freundlich et al. 2020). The CuspCore model assumes an instantaneous change of potential, followed by a relaxation to a new Jeans equilibrium. The relaxation turns out to be violent relaxation during the first orbital period, followed by phase mixing. By tracing the energy diffusion dE=dU(r) iteratively, the model reproduces the simulated DM profiles with ~10% accuracy or better. A method based on adiabatic invariants shows similar precision for moderate mass change but underestimates the DM expansion for strong gas ejection. A method based on a simple empirical relation between DM and total mass ratios makes slightly inferior predictions. The crude assumption used in CuspCore I, of energy conservation for shells that encompass a fixed DM mass, turns out to underestimate the DM response, which can be partially remedied by introducing an alternative "energy" definition. Our model is being generalized to address the differential response of a multi-component system of stars and DM in the formation of DM-deficient galaxies., Comment: 20 pages, 15 figures, accepted version

We study the driving of turbulence in star-froming disc galaxies of different masses at different epochs, using an analytic "bathtub" model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc instability and clumpy accretion via streams. The transport rate is computed with two different formalisms, with similar results. To achieve the energy balance, the disc self-regulates either the mass fraction in clumps or the turbulent viscous torque parameter. In this version of the model, the efficiency by which the stream kinetic energy is converted into turbulence is a free parameter, $\xi_a$. We find that the contributions of the three energy sources are in the same ball park, within a factor of $\sim\!2$ in all discs at all times. In haloes that evolve to a mass $\leq 10^{12}\,\Msun$ by $z=0$ ($\leq 10^{11.5}\,\Msun$ at $z\!\sim\!2$), feedback is the main driver throughout their lifetimes. Above this mass, the main driver is either transport or accretion for very low or very high values of $\xi_a$, respectively. For an assumed $\xi_a(t)$ that declines in time, galaxies in halos with present-day mass $>\!10^{12}$ M$_\odot$ make a transition from accretion to transport dominance at intermediate redshifts, $z\! \sim\!3$, when their mass was $\geq\!10^{11.5}\,\Msun$. The predicted relation between star-formation rate and gas velocity dispersion is consistent with observations., Comment: 17 pages, 6 figures. Accepted for publication in MNRAS

We address the nature of the giant clumps in high-z galaxies that undergo Violent Disc Instability, attempting to distinguish between long-lived clumps that migrate inward and short-lived clumps that disrupt by feedback. We study the evolution of clumps as they migrate through the disc using an analytic model tested by simulations and confront theory with CANDELS observations. The clump ``bathtub" model, which considers gas and stellar gain and loss, is characterized by four parameters: the accretion efficiency, the star-formation-rate (SFR) efficiency, and the outflow mass-loading factors for gas and stars. The relevant timescales are all comparable to the migration time, two-three orbital times. A clump differs from a galaxy by the internal dependence of the accretion rate on the varying clump mass. The analytic solution, involving exponential growing and decaying modes, reveals a main evolution phase during the migration, where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse of the specific SFR an observable proxy for clump age. Later, the masses and SFR approach an exponential growth with a constant specific SFR, but this phase is hypothetical as the clump disappears in the galaxy center. The model matches simulations with different, moderate feedback, both in isolated and cosmological settings. The observed clumps agree with our theoretical predictions, indicating that the massive clumps are long-lived and migrating. A non-trivial challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio., Comment: 27 pages, 10 figures

We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration., Comment: Accepted for publication to the ApJ. The simulation snapshots used in this paper will be made publicly available in a user-friendly interface. The announcement of the data release will be posted on the AGORA Project webpage (www.agorasimulations.org). Early access to the data can be provided upon request to the authors, or to the project coordinator (sroca01@ucm.es, santacruzgalaxy@gmail.com)

Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass $\geq\! 10^{12}M_\odot$ at $z\!\sim\! 2$. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing AGN-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical-friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses $>\!10^{11.3}M_\odot$, that undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion onto haloes $\geq\!10^{12}M_\odot$, these satellites heat-up the cusps in half a virial time at $z\!\sim\! 2$. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGN indeed become effective in haloes $\geq\!10^{12}M_\odot$, where the black-hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot CGM. For simulations to reproduce the dynamical-friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black-hole accretion and clumpy response to AGN., Comment: 25+10 pages, 10+8 figures

We use the Cosmic Assembly Deep Near-infrared Extragalactic Legacy Survey (CANDELS) data to study the relationship between quenching and the stellar mass surface density within the central radius of 1 kpc ($\Sigma_1$) of low-mass galaxies (stellar mass $M_* \lesssim 10^{9.5} M_\odot$) at $0.5 \leq z < 1.5$. Our sample is mass complete down to $\sim 10^9 M_\odot$ at $0.5 \leq z < 1.0$. We compare the mean $\Sigma_1$ of star-forming galaxies (SFGs) and quenched galaxies (QGs) at the same redshift and $M_*$. We find that low-mass QGs have higher $\Sigma_1$ than low-mass SFGs, similar to galaxies above $10^{10} M_\odot$. The difference of $\Sigma_1$ between QGs and SFGs increases slightly with $M_*$ at $M_* \lesssim 10^{10} M_\odot$ and decreases with $M_*$ at $M_* \gtrsim 10^{10} M_\odot$. The turnover mass is consistent with the mass where quenching mechanisms transition from internal to environmental quenching. At $0.5 \leq z < 1.0$, we find that the $\Sigma_1$ of galaxies increases by about 0.25 dex in the green valley (i.e., the transitioning region from star forming to fully quenched), regardless of their $M_*$. Using the observed specific star formation rate (sSFR) gradient in the literature as a constraint, we estimate that the quenching timescale (i.e., time spent in the transition) of low-mass galaxies is a few ($\sim4$) Gyrs at $0.5 \leq z < 1.0$. The mechanisms responsible for quenching need to gradually quench star formation in an outside-in way, i.e., preferentially ceasing star formation in outskirts of galaxies while maintaining their central star formation to increase $\Sigma_1$. An interesting and intriguing result is the similarity of the growth of $\Sigma_1$ in the green valley between low-mass and massive galaxies, which suggests that the role of internal processes in quenching low-mass galaxies is a question worthy of further investigation., Comment: 12 pages, 11 figures. Accepted by ApJ

The phenomenological study of evolving galaxy populations has shown that star forming galaxies can be quenched by two distinct processes: mass quenching and environment quenching (Peng et al. 2010). To explore the mass quenching process in local galaxies, we study the massive central disk galaxies with stellar mass above the Schechter characteristic mass. In Zhang et al. (2019), we showed that during the quenching of the massive central disk galaxies as their star formation rate (SFR) decreases, their molecular gas mass and star formation efficiency drop rapidly, but their HI gas mass remains surprisingly constant. To identify the underlying physical mechanisms, in this work we analyze the change during quenching of various structure parameters, bar frequency, and active galactic nucleus (AGN) activity. We find three closely related facts. On average, as SFR decreases in these galaxies: (1) they become progressively more compact, indicated by their significantly increasing concentration index, bulge-to-total mass ratio, and central velocity dispersion, which are mainly driven by the growth and compaction of their bulge component; (2) the frequency of barred galaxies increases dramatically, and at a given concentration index the barred galaxies have a significantly higher quiescent fraction than unbarred galaxies, implying that the galactic bar may play an important role in mass quenching; and (3) the "AGN" frequency increases dramatically from 10% on the main sequence to almost 100% for the most quiescent galaxies, which is mainly driven by the sharp increase of LINERs. These observational results lead to a self-consistent picture of how mass quenching operates., Comment: 20 pages, 11 figures. Published in ApJ