Search

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!

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

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

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

We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure identification framework. The MCPM-based reconstruction identifies filaments with higher fidelity, finding more low-prominence/diffuse filaments and better tracing the true underlying matter distribution than the DTFE-based reconstruction. Using our new filament catalogs, we find that most galaxies are located within 1.5-2.5 Mpc of a filamentary spine, with little change in the median specific star formation rate and the median galactic gas fraction with distance to the nearest filament. Instead, we introduce the filament line density, Sigma_fil(MCPM), as the total MCPM overdensity per unit length of a local filament segment, and find that this parameter is a superior predictor of galactic gas supply and quenching. Our results indicate that most galaxies are quenched and gas-poor near high-line density filaments at z<=1. At z=0, quenching in log(M*/Msun)>10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. In high-line density filaments, satellites are strongly quenched, whereas centrals have reduced star formation, but not gas fraction, at z<=0.5. We discuss the prospect of applying our new filament identification method to galaxy surveys with SDSS, DESI, Subaru PFS, etc. to elucidate the effect of large-scale structure on galaxy formation., Comment: Accepted for publication in ApJ. Data available at https://doi.org/10.5281/zenodo.10929036

We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays clear Lyman- and Balmer-breaks commensurate with this redshift. A fit to the spectrum suggests a B-type super-giant star of surface temperature $T_{\mathrm{eff,B}}\simeq15000$ K with either a redder F-type companion ($T_{\mathrm{eff,F}}\simeq6250$K) or significant dust attenuation ($A_V\simeq0.82$) along the line of sight. We also investigate the possibility that this object is a magnified young globular cluster rather than a single star. We show that the spectrum is in principle consistent with a star cluster, which could also accommodate the lack of flux variability between the two epochs. However, the lack of a counter image and the strong upper limit on the size of the object from lensing symmetry, $r\lesssim0.5$ pc, could indicate that this scenario is somewhat less likely -- albeit not completely ruled out by the current data. The presented spectrum seen at a time when the Universe was only $\sim1.2$ Gyr old showcases the ability of JWST to study early stars through extreme lensing., Comment: Accepted for publication in MNRAS letters. v2 updated to match the published version

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 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

In the Halo21 absorption modeling challenge we generated synthetic absorption spectra of the circumgalactic medium (CGM), and attempted to estimate the metallicity, temperature, and density (Z, T, and nH) of the underlying gas using observational methods. We iteratively generated and analyzed three increasingly-complex data samples: ion column densities of isolated uniform clouds, mock spectra of 1--3 uniform clouds, and mock spectra of high-resolution turbulent mixing zones. We found that the observational estimates were accurate for both uniform cloud samples, with Z, T, and nH retrieved within 0.1 dex of the source value for >90% of absorption systems. In the turbulent-mixing scenario, the mass, temperature, and metallicity of the strongest absorption components were also retrieved with high accuracy. However, the underlying properties of the subdominant components were poorly constrained because the corresponding simulated gas contributed only weakly to the H I absorption profiles. On the other hand, including additional components beyond the dominant ones did improve the fit, consistent with the true existence of complex cloud structures in the source data., Comment: 16 pages, 10 figures, submitted. Data and additional figures available upon request

The growth of galaxies in the early Universe is driven by accretion of circum- and inter-galactic gas. Simulations predict that steady streams of cold gas penetrate the dark matter halos of galaxies, providing the raw material necessary to sustain star formation. We report a filamentary stream of gas that extends for 100 kiloparsecs and connects to the massive radio galaxy 4C 41.17. The stream is detected using sub-millimeter observations of the [CI] line of atomic carbon, a tracer of neutral atomic or molecular hydrogen gas. The galaxy contains a central gas reservoir that is fueling a vigorous starburst. Our results show that the raw material for star formation can be present in cosmic streams outside galaxies., Comment: Published in Science, Volume 379, 31 March 2023 (accepted version, 31 pages)

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_{\ast}/\mathrm{M}_{\odot}) \geq 8$ to the nearest node ($d_{\mathrm{node}}$) and the nearest filament spine ($d_{\mathrm{fil}}$) to study the dependence of both median specific star formation rate () and median gas fraction (<$f_{\mathrm{gas}}$>) on these distances. We find that the of galaxies is only dependent on cosmic web environment at $z<2$, with the dependence increasing with time. At $z\leq0.5$, $8 \leq \log(M_{\ast}/\mathrm{M}_{\odot}) < 9$ galaxies are quenched at $d_{\mathrm{node}}\lesssim1$~Mpc, and have significantly-suppressed star formation at $d_{\mathrm{fil}}\lesssim1$~Mpc, trends driven mostly by satellite galaxies. At $z\leq1$, in contrast to the monotonic drop in of $\log(M_{\ast}/\mathrm{M}_{\odot}) <10$ galaxies with decreasing $d_{\mathrm{node}}$ and $d_{\mathrm{fil}}$, $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 10$ galaxies - both centrals and satellites - experience an upturn in at $d_{\mathrm{node}}\lesssim0.2$~Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in $$. Our results suggest that in the past $\sim$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, while at earlier times cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as SDSS and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST and Roman, are required to test our predictions against observations., Comment: Accepted for publication in ApJ. 5 Figures

We present results of [CII]$\,158\,\rm{\mu m}$ emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed ($\mu\sim20$) star-forming galaxy, MACS0308-zD1 at $z=6.2078\pm0.0002$. The [CII] emission line is detected with a signal-to-noise ratio $>6$ within the rest-frame UV bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the narrow [CII] has a velocity full-width-half-maximum (FWHM) of $110\pm20\,\rm{km/s}$, while broader [CII] is seen with an FWHM of $230\pm20\,\rm{km/s}$. The broader [CII] component is blueshifted ($-80\pm20\,\rm{km/s}$) with respect to the narrow [CII] component, and has a morphology which extends beyond the UV-bright clump. We find that while the narrow [CII] emission is most likely associated with zD1.1, the broader component is possibly associated with outflowing gas. Based on the non-detection of $\lambda_{\rm 158\,\mu m}$ dust continuum, we find that MACS0308-zD1's star-formation activity occurs in a dust-free environment with the stringent upper limit of infrared luminosity $\lesssim9\times10^{8}\,{\rm L_{\odot}}$. Targeting this strongly lensed faint galaxy for follow-up ALMA and JWST observations will be crucial to characterize the details of typical galaxy growth in the early Universe., Comment: 10 pages, 6 figures, 1 table; submitted to ApJ

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

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

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

The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$\mu$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $\mu > 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe., Comment: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

Detailed observations of star forming galaxies at high redshift are critical to understand the formation and evolution of the earliest galaxies. Gravitational lensing provides an important boost, allowing observations at physical scales unreachable in unlensed galaxies. We present three lensed galaxies from the RELICS survey at $z_{phot} = 6 - 10$, including the most highly magnified galaxy at $z_{phot} \sim 6$ (WHL0137-zD1, dubbed the Sunrise Arc), the brightest known lensed galaxy at $z_{phot} \sim 6$ (MACS0308-zD1), and the only spatially resolved galaxy currently known at $z_{phot} \sim 10$ (SPT0615-JD). The Sunrise Arc contains seven star-forming clumps with delensed radii as small as 3 pc, the smallest spatial scales yet observed in a $z>6$ galaxy, while SPT0615-JD contains features measuring a few tens of parsecs. MACS0308-zD1 contains a $r\sim 30$ pc clump with a star formation rate (SFR) of $\sim 3 M_{\odot} \textrm{ yr}^{-1}$, giving it a SFR surface density of $\Sigma_{SFR} \sim 10^3 M_{\odot}\textrm{ yr}^{-1}\textrm{ kpc}^{-2}$. These galaxies provide a unique window into small scale star formation during the Epoch of Reionization. They will be excellent targets for future observations with JWST, including one approved program targeting the Sunrise Arc., Comment: 15 pages, 7 figures, submitted to ApJ

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

Dwarf galaxies are thought to quench primarily due to environmental processes most typically occurring in galaxy groups and clusters or around single, massive galaxies. However, at earlier epochs, ($5 < z < 2$), the collapse of large scale structure (forming Zel'dovich sheets and subsequently filaments of the cosmic web) can produce volume-filling accretion shocks which elevate large swaths of the intergalactic medium (IGM) in these structures to a hot ($T>10^6$ K) phase. We study the impact of such an event on the evolution of central dwarf galaxies ($5.5 < \log M_* < 8.5$) in the field using a spatially large, high resolution cosmological zoom simulation which covers the cosmic web environment between two protoclusters. We find that the shock-heated sheet acts as an environmental quencher much like clusters and filaments at lower redshift, creating a population of quenched, central dwarf galaxies. Even massive dwarfs which do not quench are affected by the shock, with reductions to their sSFR and gas accretion. This process can potentially explain the presence of isolated quenched dwarf galaxies, and represents an avenue of pre-processing, via which quenched satellites of bound systems quench before infall., Comment: 15 pages, 10 figures. Submitted to MNRAS

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

Using a novel suite of cosmological simulations zooming in on a Mpc-scale intergalactic sheet or "pancake" at z~3-5, we conduct an in-depth study of the thermal properties and HI content of the warm-hot intergalactic medium (WHIM) at those redshifts. The simulations span nearly three orders of magnitude in gas-cell mass, from ~(7.7x10^6-1.5x10^4)Msun, one of the highest resolution simulations of such a large patch of the inter-galactic medium (IGM) to date. At z~5, a strong accretion shock develops around the main pancake following a collision between two smaller sheets. Gas in the post-shock region proceeds to cool rapidly, triggering thermal instabilities and the formation of a multiphase medium. We find neither the mass, nor the morphology, nor the distribution of HI in the WHIM to be converged at our highest resolution. Interestingly, the lack of convergence is more severe for the less dense, more metal-poor, intra-pancake medium (IPM) in between filaments and far from any star-forming galaxies. As the resolution increases, the IPM develops a shattered structure, with ~kpc scale clouds containing most of the HI. From our lowest to highest resolution, the covering fraction of metal-poor (Z<10^{-3}Zsun) Lyman-limit systems (NHI>10^{17.2}/cm^2) in the IPM at z~4 increases from (3-15)%, while that of Damped Lyman-alpha Absorbers (NHI>10^{20}/cm^2) with similar metallicity increases threefold, from (0.2-0.6)%, with no sign of convergence. We find that a necessary condition for the formation of a multiphase, shattered structure is resolving the cooling length, lcool=cs*tcool, at T~10^5K. If this scale is unresolved, gas "piles up" at these temperatures and cooling to lower temperatures becomes very inefficient. We conclude that state-of-the-art cosmological simulations are still unable to resolve the multi-phase structure of the low-density IGM, with potentially far-reaching implications., Comment: 24 pages, 15 (rather large) figures, 2 tables, 1 appendix (main body 21 pages) Comments welcome! Main observational results in sections 3-4 (Table 2) Main theory result in section 6.1

We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with $\sim \!\! 25\ \!\rm{pc}$ maximum resolution, targeting main sequence galaxies at $1\!<\!z\!<\!3$. The clumps are classified as long-lived clumps (LLCs) or short-lived clumps (SLCs) based on their longevity in the simulations. We then train neural networks to detect and classify the simulated clumps in mock, multi-color, dusty and noisy HST-like images. The clumps are detected using an encoder-decoder convolutional neural network (CNN), and are classified according to their longevity using a vanilla CNN. Tests using the simulations show our detector and classifier to be $\sim80\%$ complete and $\sim80\%$ pure for clumps more massive than $\sim10^{7.5}\rm{M_\odot}$. When applied to observed galaxies in the CANDELS/GOODS S+N fields, we find both types of clumps to appear in similar abundances in the simulations and the observations. LLCs are, on average, more massive than SLCs by $\sim 0.5\ \rm{dex}$, and they dominate the clump population above $M_{\rm c}\gtrsim 10^{7.6}\ \rm{M_\odot}$. LLCs tend to be found closer to the galactic centre, indicating clump migration to the centre or preferential formation at smaller radii. The LLCs are found to reside in high mass galaxies, indicating better clump survivability under supernova feedback there, due to clumps being more massive in these galaxies. We find the clump properties in the simulations and the observations to agree within a factor of $\sim\! 2$., Comment: 17 pages, 13 figures. Accepted for publication in MNRAS

We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range $\mathrm{M}_* = 10^{9-10} \ \mathrm{M}_\odot$, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below $\mathrm{M}_*=10^{9} \ \mathrm{M}_\odot$. We find that the low mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite apparent variation in the late time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path towards more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations., Comment: 19 pages, 8 figures. Figure 4 presents the main result. Code used in this work may be accessed at github.com/IQcollaboratory/orchard. Submitted to ApJ

We analyse the distribution and origin of OVI in the Circumgalactic Medium (CGM) of dark-matter haloes of $\sim 10^{12}$ M$_\odot$ at $z\sim1$ in the VELA cosmological zoom-in simulations. We find that the OVI in the inflowing cold streams is primarily photoionized, while in the bulk volume it is primarily collisionally ionized. The photoionized component dominates the observed column density at large impact parameters ($\gtrsim 0.3 R_{\rm vir}$), while the collisionally ionized component dominates closer in. We find that most of the collisional OVI, by mass, resides in the relatively thin boundaries of the photoionized streams. We discuss how the results are in agreement with analytic predictions of stream and boundary properties, and their compatibility with observations. This allows us to predict the profiles of OVI and other ions in future CGM observations and provides a toy model for interpreting them., Comment: 22 pages, 16 figures, submitted to MNRAS

We present an analytic toy model for the radiation produced by the interaction between the cold streams thought to feed massive halos at high redshift and their hot CGM. We begin by deriving cosmologically motivated parameters for the streams as they enter the halo virial radius, $R_{\rm v}$, as a function of halo mass and redshift. For $10^{12}M_{\odot}$ halos at $z=2$, we find the Hydrogen number density in streams to be $n_{\rm H,s}\sim (0.1-5)\times 10^{-2}{\rm cm}^{-3}$, a factor of $\delta \sim (30-300)$ times denser than the hot CGM density, while the stream radii are in the range $R_{\rm s}\sim (0.03-0.50)R_{\rm v}$. As the streams accelerate towards the halo centre, they become denser and narrower. The stream-CGM interaction induces Kelvin-Helmholtz Instability (KHI), which leads to entrainment of CGM mass by the stream and therefore to stream deceleration by momentum conservation. Assuming that the entrainment rates derived by Mandelker et al. 2019 in the absence of gravity can be applied locally at each halocentric radius, we derive equations of motion for the stream in the halo. Using these, we derive the net acceleration, mass growth, and energy dissipation induced by the stream-CGM interaction, as a function of halo mass and redshift, for different CGM density profiles. For the range of model parameters considered, we find that the interaction can induce dissipation luminosities $L_{\rm diss}>10^{42}~{\rm erg~s^{-1}}$ within $\le 0.6 R_{\rm v}$ of halos with $M_{\rm v}>10^{12}M_{\odot}$ at $z=2$, with the emission scaling with halo mass and redshift approximately as $\propto M_{\rm v}\,(1+z)^2$. The magnitude and spatial extent of the emission produced in massive halos at high redshift is consistent with observed Ly$\alpha$ blobs, though better treatment of the UV background and self-shielding is needed to solidify this conclusion., Comment: 16 pages, 6 figures, resubmitted to MNRAS with minor revisions following referee report. Model has been refined with fewer simplifying assumptions, but with no significant changes to our results. Several subtleties and caveats have been more clearly presented

We study the effects of Kelvin Helmholtz Instability (KHI) on the cold streams that feed massive halos at high redshift, generalizing our earlier results to include the effects of radiative cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, $t_{\rm shear}$. This can be converted into a critical stream radius, $R_{\rm s,crit}$, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,mix}$. If $R_{\rm s}R_{\rm s,crit}$, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses onto the stream, and can increase the mass of cold gas by a factor of $\sim 3$ within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T<5\times 10^4{\rm K}$, suggesting much of it will be emitted in Ly$\alpha$., Comment: 25 pages, 15 figures, 1 table, brief appendix. Revised version submitted to MNRAS. Some new simulations added to analysis, and mild changes to discussion of turbulence. Main conclusions unchanged. Comments still welcome!

We present a new cosmological zoom-in simulation, where the zoom region consists of two halos with virial mass M_v~5x10^{12}M_{sun} and a ~Mpc long cosmic filament connecting them at z~2. Using this simulation, we study the evolution of the intergalactic medium in between these two halos at unprecedented resolution. At 5>z>3, the two halos are found to lie in a large intergalactic sheet, or "pancake", consisting of multiple co-planar dense filaments along which nearly all halos with M_v>10^9M_{sun} are located. This sheet collapses at z~5 from the merger of two smaller sheets. The strong shock generated by this merger leads to thermal instabilities in the post-shock region, and to a shattering of the sheet resulting in <~kpc scale clouds with temperatures of T>~2x10^4K and densities of n>~10^{-3}cm^{-3}, which are pressure confined in a hot medium with T~10^6K and n>~10^{-5}cm^{-3}. When the sheet is viewed face on, these cold clouds have neutral hydrogen column densities of N_{HI}>10^{17.2}cm^{-2}, making them detectable as Lyman limit systems, though they lie well outside the virial radius of any halo and even well outside the dense filaments. Their chemical composition is pristine, having zero metalicity, similar to several recently observed systems. Since these systems form far from any galaxies, these results are robust to galaxy formation physics, resulting purely from the collapse of large scale structure and radiative cooling, provided sufficient spatial resolution is available., Comment: 7 pages, 5 figures. Resubmitted to ApJL following first referee report. Movie generalizing figure 2 available with ancillary material. Additional comments welcome!

Self-gravitating gaseous filaments exist on many astrophysical scales, from sub-pc filaments in the interstellar medium to Mpc scale streams feeding galaxies from the cosmic web. These filaments are often subject to Kelvin-Helmotz Instability (KHI) due to shearing against a confining background medium. We study the nonlinear evolution of KHI in pressure-confined self-gravitating gas streams initially in hydrostatic equilibrium, using analytic models and hydrodynamic simulations, not including radiative cooling. We derive a critical line-mass as a function of the stream Mach number and density contrast with respect to the background, $\mu_{cr}(M_b,\delta_c)\le 1$, where $\mu=1$ is normalized to the maximal line mass for which initial hydrostatic equilibrium is possible. For $\mu<\mu_{cr}$, KHI dominates the stream evolution. A turbulent shear layer expands into the background and leads to stream deceleration at a similar rate to the non-gravitating case. However, with gravity, penetration of the shear layer into the stream is halted at roughly half the initial stream radius by stabilizing buoyancy forces, significantly delaying total stream disruption. Streams with $\mu_{cr}<\mu\le 1$ fragment and form round, long-lived clumps by gravitational instability (GI), with typical separations roughly 8 times the stream radius, similar to the case without KHI. When KHI is still somewhat effective, these clumps are below the spherical Jeans mass and are partially confined by external pressure, but they approach the Jeans mass as $\mu\rightarrow 1$ and GI dominates. We discuss potential applications of our results to filaments in the ISM and dense streams feeding galaxies at high redshift., Comment: 23 pages, 13 figures, accepted in MNRAS. Minor updates to text and additional discussion on tidal streams

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise., Comment: White Paper submitted to the Astro2020 Decadal Survey

The majority of massive star-forming galaxies at $z\sim2$ have velocity gradients suggestive of rotation, in addition to large amounts of disordered motions. In this paper, we demonstrate that it is challenging to distinguish the regular rotation of a disk galaxy from the orbital motions of merging galaxies with seeing-limited data. However, the merger fractions at $z\sim2$ are likely too low for this to have a large effect on measurements of disk fractions. To determine how often mergers pass for disks, we look to galaxy formation simulations. We analyze $\sim$24000 synthetic images and kinematic maps of 31 high-resolution simulations of isolated galaxies and mergers at $z\sim2$. We determine if the synthetic observations pass criteria commonly used to identify disk galaxies, and whether the results are consistent with their intrinsic dynamical states. Galaxies that are intrinsically mergers pass the disk criteria for anywhere from 0 to 100$\%$ of sightlines. The exact percentage depends strongly on the specific disk criteria adopted, and weakly on the separation of the merging galaxies. Therefore, one cannot tell with certainty whether observations of an individual galaxy indicate a merger or a disk. To estimate the fraction of mergers passing as disks in current kinematics samples, we combine the probability that a merger will pass as a disk with theoretical merger fractions from a cosmological simulation. Taking the latter at face-value, the observed disk fractions are overestimated by small amounts: at most by $5\%$ at high stellar mass ($10^{10-11}$ M$_{\odot}$) and $15\%$ at low stellar mass ($10^{9-10}$ M$_{\odot}$)., Comment: 15 pages, 10 figures, accepted for publication in ApJ

The circumgalactic medium (CGM), i.e. the gaseous haloes around galaxies, is both the reservoir of gas that fuels galaxy growth and the repository of gas expelled by galactic winds. Most cosmological, hydrodynamical simulations focus their computational effort on the galaxies themselves and treat the CGM more coarsely, which means small-scale structure cannot be resolved. We get around this issue by running zoom-in simulations of a Milky Way-mass galaxy with standard mass refinement and additional uniform spatial refinement within the virial radius. This results in a detailed view of its gaseous halo at unprecedented (1 kpc) uniform resolution with only a moderate increase in computational time. The improved spatial resolution does not impact the central galaxy or the average density of the CGM. However, it drastically changes the radial profile of the neutral hydrogen column density, which is enhanced at galactocentric radii larger than 40 kpc. The covering fraction of Lyman-Limit Systems within 150 kpc is almost doubled. We therefore conclude that some of the observational properties of the CGM are strongly resolution dependent. Increasing the resolution in the CGM, without increasing the resolution of the galaxies, is a promising and computationally efficient method to push the boundaries of state-of-the-art simulations., Comment: Accepted for publication in MNRAS Letters. Revised version: minor changes

We study the effects of Kelvin-Helmholtz instability (KHI) on the cold streams that feed high-redshift galaxies through their hot haloes, generalizing our earlier analyses of a 2D slab to a 3D cylinder, but still limiting our analysis to the adiabatic case with no gravity. We combine analytic modeling and numerical simulations in the linear and non-linear regimes. For subsonic or transonic streams with respect to the halo sound speed, the instability in 3D is qualitatively similar to 2D, but progresses at a faster pace. For supersonic streams, the instability grows much faster in 3D and can be qualitatively different due to azimuthal modes, which introduce a strong dependence on the initial width of the stream-background interface. Using analytic toy models and approximations supported by high-resolution simulations, we apply our idealized hydrodynamical analysis to the astrophysical scenario. The upper limit for the radius of a stream that disintegrates prior to reaching the central galaxy is ~70% larger than the 2D estimate; it is in the range (0.5-5)% of the halo virial radius, decreasing with increasing stream density and velocity. Stream disruption generates a turbulent mixing zone around the stream with velocities at the level of ~20% of the initial stream velocity. KHI can cause significant stream deceleration and energy dissipation in 3d, contrary to 2D estimates. For typical streams, up to (10-50)% of the gravitational energy gained by inflow down the dark-matter halo potential can be dissipated, capable of powering Lyman-alpha blobs if most of it is dissipated into radiation., Comment: Accepted to MNRAS, updated to accepted version. Very minor changes to astrophysical conclusions, but significant additions to hydrodynamical analysis, in particular shear layer growth and bulk kinetic energy dissipation. Main astrophysical applications in sections 1 and 5-7, totaling 11 pages and 3 figures

As part of our long-term campaign to understand how cold streams feed massive galaxies at high redshift, we study the Kelvin-Helmholtz instability (KHI) of a supersonic, cold, dense gas stream as it penetrates through a hot, dilute circumgalactic medium (CGM). A linear analysis (Paper I) showed that, for realistic conditions, KHI may produce nonlinear perturbations to the stream during infall. Therefore, we proceed here to study the nonlinear stage of KHI, still limited to a two-dimensional slab with no radiative cooling or gravity. Using analytic models and numerical simulations, we examine stream breakup, deceleration and heating via surface modes and body modes. The relevant parameters are the density contrast between stream and CGM ($\delta$), the Mach number of the stream velocity with respect to the CGM ($M_{\rm b}$) and the stream radius relative to the halo virial radius ($R_{\rm s}/R_{\rm v}$). We find that sufficiently thin streams disintegrate prior to reaching the central galaxy. The condition for breakup ranges from $R_{\rm s} < 0.03 R_{\rm v}$ for $(M_{\rm b} \sim 0.75, \delta \sim 10)$ to $R_{\rm s} < 0.003 R_{\rm v}$ for $(M_{\rm b} \sim 2.25, \delta \sim 100)$. However, due to the large stream inertia, KHI has only a small effect on the stream inflow rate and a small contribution to heating and subsequent Lyman-$\alpha$ cooling emission., Comment: The main astrophysical results are Figure 22 and Figure 23. Final 7 pages are appendices. Accepted to MNRAS

Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged], Comment: 23 pages, 14 figures. Accepted by ApJ. The clump catalog will be published online in ApJ. It can also be requested through emails to the first author or this link http://faculty.missouri.edu/guoyic/data/clump/

We propose that cold filamentary accretion in massive galaxies at high redshift can lead to the formation of star-forming clumps in the halos of these galaxies without the presence of dark matter sub-structure. In certain cases, these clumps can be the birth places of metal poor globular-clusters (MP GCs). Using cosmological simulations, we show that narrow streams of dense gas feeding massive galaxies from the cosmic web can fragment, producing star-forming clumps. We then derive an analytical model for the properties of streams as a function of halo mass and redshift, and assess when these are gravitationally unstable, when this can lead to collapse and star-formation in the halo, and when it may result in the formation of MP GCs. For stream metalicities >~1% solar, this is likely to occur at z>4.5. At z~6, the collapsing clouds have masses of (5-10)x10^7 M_{sun} and the average stream pressure is ~10^6 K cm^{-3}. The conditions for GC formation are met in the extremely turbulent "eyewall" at ~0.3 R_{vir}, where counter-rotating streams can collide, driving very large densities. Our scenario can account for the observed kinematics and spatial distribution of MP GCs, the correlation between their mass and metalicity, and the mass ratio between the GC system and the host halo. For MW mass halos, we infer that ~30% of MP GCs could have formed in this way, the remainder likely accreted in mergers. Our predictions for GC formation along circumgalactic filaments at high-redshift are testable with JWST., Comment: Published in ApJ. Updated to published version. Only minor changes from previous version. Additional references and discussion, new table summarizing model parameters

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes., Comment: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields., Comment: Accepted to MNRAS. Minor changes from previous version

We study the evolution of giant clumps in high-z disc galaxies using AMR cosmological simulations at redshifts z=6-1. Our sample consists of 34 galaxies, of halo masses 10^{11}-10^{12}M_s at z=2, run with and without radiation pressure (RP) feedback from young stars. While RP has little effect on the sizes and global stability of discs, it reduces the amount of star-forming gas by a factor of ~2, leading to a decrease in stellar mass by a similar factor by z~2. Both samples undergo violent disc instability (VDI) and form giant clumps of masses 10^7-10^9M_s at a similar rate, though RP significantly reduces the number of long-lived clumps. When RP is (not) included, clumps with circular velocity <40(20)km/s, baryonic surface density <200(100)M_s/pc^2 and baryonic mass <10^{8.2}(10^{7.3})M_s are short-lived, disrupted in a few free-fall times. The more massive and dense clumps survive and migrate toward the disc centre over a few disc orbital times. In the RP simulations, the distribution of clump masses and star-formation rates (SFRs) normalized to their host disc is very similar at all redshifts. They exhibit a truncated power-law with a slope slightly shallower than -2. Short-lived clumps preferentially have young stellar ages, low masses, high gas fractions and specific SFRs (sSFR), and they tend to populate the outer disc. The sSFR of massive, long-lived clumps declines with age as they migrate towards the disc centre, producing gradients in mass, stellar age, gas fraction, sSFR and metallicity that distinguish them from short-lived clumps. Ex situ mergers make up ~37% of the mass in clumps and ~29% of the SFR. They are more massive and with older stellar ages than the in situ clumps, especially near the disc edge. Roughly half the galaxies at redshifts z=4-1 are clumpy over a wide range of stellar mass, with clumps accounting for ~3-30% of the SFR but ~0.1-3% of the stellar mass., Comment: Accepted to MNRAS. Updated measurements of outflows from clumps and new study of clump virial parameters. Other small changes

We study the evolution of global shapes of galaxies using cosmological simulations. The shapes refer to the components of dark matter (DM), stars and gas at the stellar half-mass radius. Most galaxies undergo a characteristic compaction event into a blue nugget at $z\sim2-4$, which marks the transition from a DM-dominated central body to a self-gravitating baryonic core. We find that in the high-$z$, DM-dominated phase, the stellar and DM systems tend to be triaxial, preferentially prolate and mutually aligned. The elongation is supported by an anisotropic velocity dispersion that originates from the assembly of the galaxy along a dominant large-scale filament. We estimate that torques by the dominant halo are capable of inducing the elongation of the stellar system and its alignment with the halo. Then, in association with the transition to self-gravity, small-pericenter orbits puff up and the DM and stellar systems evolve into a more spherical and oblate configuration, aligned with the gas disc and associated with rotation. This transition typically occurs when the stellar mass is $\sim 10^9$ M$_\odot$ and the escape velocity in the core is $\sim 100$ km s$^{-1}$, indicating that supernova feedback may be effective in keeping the core DM-dominated and the system prolate. The early elongated phase itself may be responsible for the compaction event, and the transition to the oblate phase may be associated with the subsequent quenching in the core., Comment: 23 pages, 20 figures; accepted for publication in MNRAS

We utilize zoom-in cosmological simulations to study the nature of violent disc instability (VDI) in clumpy galaxies at high redshift, $z=1$--$5$. Our simulated galaxies are not in the ideal state assumed in Toomre instability, of linear fluctuations in an isolated, uniform, rotating disk. There, instability is characterised by a $Q$ parameter below unity, and lower when the disk is thick. Instead, the high-redshift discs are highly perturbed. Over long periods they consist of non-linear perturbations, compact massive clumps and extended structures, with new clumps forming in inter-clump regions. This is while the galaxy is subject to frequent external perturbances. We compute the local, two-component $Q$ parameter for gas and stars, smoothed on a $\sim1~{\rm kpc}$ scale to capture clumps of $10^{8-9}~{\rm M}_\odot$. The $Q<1$ regions are confined to collapsed clumps due to the high surface density there, while the inter-clump regions show $Q$ significantly higher than unity. Tracing the clumps back to their relatively smooth Lagrangian patches, we find that $Q$ prior to clump formation typically ranges from unity to a few. This is unlike the expectations from standard Toomre instability. We discuss possible mechanisms for high-$Q$ clump formation, e.g. rapid turbulence decay leading to small clumps that grow by mergers, non-axisymmetric instability, or clump formation induced by non-linear perturbations in the disk. Alternatively, the high-$Q$ non-linear VDI may be stimulated by the external perturbations such as mergers and counter-rotating streams. The high $Q$ may represent excessive compressive modes of turbulence, possibly induced by tidal interactions., Comment: Accepted for publication in MNRAS. 20 pages, 21 figures

Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end., Comment: 27 pages, 17 figures. Accepted for publication in MNRAS

Using cosmological simulations, we address the interplay between structure and star formation in high-redshift galaxies via the evolution of surface density profiles. Our sample consists of 26 galaxies evolving in the redshift range $z=7-1$, spanning the stellar mass range $(0.2-6.4)\times 10^{10}M_\odot$ at $z=2$. We recover the main trends by stacking the profiles in accordance to their evolution phases. Following a wet compaction event that typically occurs when the stellar mass is $\sim10^{9.5}~M_{\odot}$ at $z\sim2-4$, the gas develops a cusp inside the effective radius, associated with a peak in star-formation rate (SFR). The SFR peak and the associated feedback, in the absence of further gas inflow to the centre, marks the onset of gas depletion from the central 1 kpc, leading to quenching of the central SFR. An extended, star-forming ring that forms by fresh gas during the central quenching process shows as a rising specific SFR (sSFR) profile, which is interpreted as inside-out quenching. Before quenching, the stellar density profile grows self-similarly, maintaining its log-log shape because the sSFR is similar at all radii. During the quenching process, the stellar density saturates to a constant value, especially in the inner 1 kpc. The stellar mass and SFR profiles deduced from observations show very similar shapes, consistent with the scenario of wet compaction leading to inside-out quenching and the subsequent saturation of a dense stellar core. We predict a cuspy gas profile during the blue nugget phase, and a gas-depleted core, sometimes surrounded by a ring, in the post-blue nugget phase., Comment: 24 pages, 15 figures. Accepted for publication in MNRAS

We use cosmological simulations to study a characteristic evolution pattern of high redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (blue nuggets) at z~4-2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z~2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star-formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc ($\Sigma_1$). Massive galaxies quench earlier, faster, and at a higher $\Sigma_1$ than lower-mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-radius-mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum $\Sigma_1$ anti-correlated with the initial spin parameter, as predicted by Dekel & Burkert (2014). The central quenching is triggered by the high SFR and stellar/supernova feedback (possibly also AGN feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing., Comment: Resubmitted to MNRAS after responding to referee's comments; Updated and added two figures