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 clu
mps 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.
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. T
he 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.
We analyse the distribution and origin of OVI in the Circumgalactic Medium (CGM) of dark-matter haloes of $sim 10^{12}$ M$_odot$ at $zsim1$ 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.
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.
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.
Cosmological simulations of galaxies have typically produced too many stars at early times. We study the global and morphological effects of radiation pressure (RP) in eight pairs of high-resolution cosmological galaxy formation simulations. We find
that the additional feedback suppresses star formation globally by a factor of ~2. Despite this reduction, the simulations still overproduce stars by a factor of ~2 with respect to the predictions provided by abundance matching methods for halos more massive than 5E11 Msun/h (Behroozi, Wechsler & Conroy 2013). We also study the morphological impact of radiation pressure on our simulations. In simulations with RP the average number of low mass clumps falls dramatically. Only clumps with stellar masses Mclump/Mdisk <= 5% are impacted by the inclusion of RP, and RP and no-RP clump counts above this range are comparable. The inclusion of RP depresses the contrast ratios of clumps by factors of a few for clump masses less than 5% of the disk masses. For more massive clumps, the differences between and RP and no-RP simulations diminish. We note however, that the simulations analyzed have disk stellar masses below about 2E10 Msun/h. By creating mock Hubble Space Telescope observations we find that the number of clumps is slightly reduced in simulations with RP. However, since massive clumps survive the inclusion of RP and are found in our mock observations, we do not find a disagreement between simulations of our clumpy galaxies and observations of clumpy galaxies. We demonstrate that clumps found in any single gas, stellar, or mock observation image are not necessarily clumps found in another map, and that there are few clumps common to multiple maps.
