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.
We describe simple useful toy models for key processes of galaxy formation in its most active phase, at z > 1, and test the approximate expressions against the typical behaviour in a suite of high-resolution hydro-cosmological simulations of massive
galaxies at z = 4-1. We address in particular the evolution of (a) the total mass inflow rate from the cosmic web into galactic haloes based on the EPS approximation, (b) the penetration of baryonic streams into the inner galaxy, (c) the disc size, (d) the implied steady-state gas content and star-formation rate (SFR) in the galaxy subject to mass conservation and a universal star-formation law, (e) the inflow rate within the disc to a central bulge and black hole as derived using energy conservation and self-regulated Q ~ 1 violent disc instability (VDI), and (f) the implied steady state in the disc and bulge. The toy models provide useful approximations for the behaviour of the simulated galaxies. We find that (a) the inflow rate is proportional to mass and to (1+z)^5/2, (b) the penetration to the inner halo is ~50% at z = 4-2, (c) the disc radius is ~5% of the virial radius, (d) the galaxies reach a steady state with the SFR following the accretion rate into the galaxy, (e) there is an intense gas inflow through the disc, comparable to the SFR, following the predictions of VDI, and (f) the galaxies approach a steady state with the bulge mass comparable to the disc mass, where the draining of gas by SFR, outflows and disc inflows is replenished by fresh accretion. Given the agreement with simulations, these toy models are useful for understanding the complex phenomena in simple terms and for back-of-the-envelope predictions.
Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations res
olve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Ways dwarf satellites.
We develop a technique to investigate the possibility that some of the recently discovered ultra-faint dwarf satellites of the Milky Way might be cusp caustics rather than gravitationally self-bound systems. Such cusps can form when a stream of stars
folds, creating a region where the projected 2-D surface density is enhanced. In this work, we construct a Poisson maximum likelihood test to compare the cusp and exponential models of any substructure on an equal footing. We apply the test to the Hercules dwarf (d ~ 113 kpc, M_V ~ -6.2, e ~ 0.67). The flattened exponential model is strongly favored over the cusp model in the case of Hercules, ruling out at high confidence that Hercules is a cusp catastrophe. This test can be applied to any of the Milky Way dwarfs, and more generally to the entire stellar halo population, to search for the cusp catastrophes that might be expected in an accreted stellar halo.
Fully cosmological, high resolution N-Body + SPH simulations are used to investigate the chemical abundance trends of stars in simulated stellar halos as a function of their origin. These simulations employ a physically motivated supernova feedback r
ecipe, as well as metal enrichment, metal cooling and metal diffusion. As presented in an earlier paper, the simulated galaxies in this study are surrounded by stellar halos whose inner regions contain both stars accreted from satellite galaxies and stars formed in situ in the central regions of the main galaxies and later displaced by mergers into their inner halos. The abundance patterns ([Fe/H] and [O/Fe]) of halo stars located within 10 kpc of a solar-like observer are analyzed. We find that for galaxies which have not experienced a recent major merger, in situ stars at the high [Fe/H] end of the metallicity distribution function are more [alpha/Fe]-rich than accreted stars at similar [Fe/H]. This dichotomy in the [O/Fe] of halo stars at a given [Fe/H] results from the different potential wells within which in situ and accreted halo stars form. These results qualitatively match recent observations of local Milky Way halo stars. It may thus be possible for observers to uncover the relative contribution of different physical processes to the formation of stellar halos by observing such trends in the halo populations of the Milky Way, and other local L* galaxies.
We investigate the formation of the stellar halos of four simulated disk galaxies using high resolution, cosmological SPH + N-Body simulations. These simulations include a self-consistent treatment of all the major physical processes involved in gala
xy formation. The simulated galaxies presented here each have a total mass of ~10^12 M_sun, but span a range of merger histories. These simulations allow us to study the competing importance of in-situ star formation (stars formed in the primary galaxy) and accretion of stars from subhalos in the building of stellar halos in a LambdaCDM universe. All four simulated galaxies are surrounded by a stellar halo, whose inner regions (r < 20 kpc) contain both accreted stars, and an in-situ stellar population. The outer regions of the galaxies halos were assembled through pure accretion and disruption of satellites. Most of the in-situ halo stars formed at high redshift out of smoothly accreted cold gas in the inner 1 kpc of the galaxies potential wells, possibly as part of their primordial disks. These stars were displaced from their central locations into the halos through a succession of major mergers. We find that the two galaxies with recently quiescent merger histories have a higher fraction of in-situ stars (~20-50%) in their inner halos than the two galaxies with many recent mergers (~5-10% in-situ fraction). Observational studies concentrating on stellar populations in the inner halo of the Milky Way will be the most affected by the presence of in-situ stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic feature of galaxy formation.
