Violent Disk Instability Research Articles

Non-linear violent disc instability with high Toomre'sQin high-redshift clumpy disc galaxies

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.

Compaction and quenching of high-z galaxies in cosmological simulations: blue and red nuggets

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.

Four phases of angular-momentum buildup in high-z galaxies: from cosmic-web streams through an extended ring to disc and bulge

We study the angular-momentum (AM) buildup in high-$z$ massive galaxies using high-resolution cosmological simulations. The AM originates in co-planar streams of cold gas and merging galaxies tracing cosmic-web filaments, and it undergoes four phases of evolution. (I) Outside the halo virial radius ($R_{\rm v}\!\sim\!100\,{\rm kpc}$), the elongated streams gain AM by tidal torques with a specific AM (sAM) $\sim\!1.7$ times the dark-matter (DM) spin due to the gas' higher quadrupole moment. This AM is expressed as stream impact parameters, from $\sim\!0.3R_{\rm v}$ to counter rotation. (II) In the outer halo, while the incoming DM mixes with the existing halo of lower sAM to a spin $\lambda_{\rm dm}\!\sim\!0.04$, the cold streams transport the AM to the inner halo such that their spin in the halo is $\sim\!3\lambda_{\rm dm}$. (III) Near pericenter, the streams dissipate into an irregular rotating ring extending to $\sim\!0.3R_{\rm v}$ and tilted relative to the inner disc. Torques exerted partly by the disc make the ring gas lose AM, spiral in, and settle into the disc within one orbit. The ring is observable with 30\% probability as a damped Lyman-$\alpha$ absorber. (IV) Within the disc, $<\!0.1R_{\rm v}$, torques associated with violent disc instability drive AM out and baryons into a central bulge, while outflows remove low-spin gas, introducing certain sensitivity to feedback strength. Despite the different AM histories of gas and DM, the disc spin is comparable to the DM-halo spin. Counter rotation can strongly affect disc evolution.

CLUMPY GALAXIES IN CANDELS. I. THE DEFINITION OF UV CLUMPS AND THE FRACTION OF CLUMPY GALAXIES AT 0.5 &lt;z&lt; 3

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

Early formation of massive, compact, spheroidal galaxies with classical profiles by violent disc instability or mergers

We address the formation of massive stellar spheroids between redshifts $z=4$ and 1 using a suite of AMR hydro-cosmological simulations. The spheroids form as bulges, and the spheroid mass growth is partly driven by violent disc instability (VDI) and partly by mergers. A kinematic decomposition to disc and spheroid yields that the mass fraction in the spheroid is between 50\% and 90\% and is roughly constant in time, consistent with a cosmological steady state of VDI discs that are continuously fed from the cosmic web. The density profile of the spheroid is typically "classical", with a Sersic index $n = 4.5\pm 1$, independent of whether it grew by mergers or VDI and independent of the feedback strength. The disc is characterized by $n=1.5\pm 0.5$, and the whole galaxy by $n=3\pm 1$. The high-redshift spheroids are compact due to the dissipative inflow of gas and the high universal density. The stellar surface density within the effective radius of each galaxy as it evolves remains roughly constant in time after its first growth. For galaxies of a fixed stellar mass, the surface density is higher at higher redshifts.

NO MORE ACTIVE GALACTIC NUCLEI IN CLUMPY DISKS THAN IN SMOOTH GALAXIES ATz∼ 2 IN CANDELS/3D-HST

We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3 < z < 2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that despite being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z ∼ 2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z ∼ 2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile, the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z ∼ 1.85—whether violent disk instabilities or secular processes—are as efficient in smooth galaxies as they are in clumpy galaxies.

The population of giant clumps in simulated high-z galaxies: in situ and ex situ migration and survival

We study the properties of giant clumps and their radial gradients in high-$z$ disc galaxies using AMR cosmological simulations. Our sample consists of 770 snapshots in the redshift range $z=4-1$ from 29 galaxies that at $z=2$ span the stellar mass range $(0.2-3)\times 10^{11}M_{\odot}$. Extended gas discs exist in 83% of the snapshots. Clumps are identified by gas density in 3D and their stellar and dark matter components are considered thereafter. While most of the overdensities are diffuse and elongated, 91% of their mass and 83% of their star-fromation rate (SFR) are in compact round clumps. Nearly all galaxies have a central, massive bulge clump, while 70% of the discs show off-center clumps, 3-4 per galaxy. The fraction of clumpy discs peaks at intermediate disc masses. Clumps are divided based on dark-matter content into $\textit{in-situ}$ and $\textit{ex-situ}$, originating from violent disc instability (VDI) and minor mergers respectively. 60% of the discs are in a VDI phase showing off-center $\textit{in-situ}$ clumps, which contribute 1-7% of the disc mass and 5-45% of its SFR. The $\textit{in-situ}$ clumps constitute 75% of the off-center clumps in terms of number and SFR but only half the mass, each clump containing on average 1% of the disc mass and 6% of its SFR. They have young stellar ages, $100-400 {\rm Myr}$, and high specific SFR (sSFR), $1-10 {\rm Gyr}^{-1}$. They exhibit gradients resulting from inward clump migration, where the inner clumps are somewhat more massive and older, with lower gas fraction and sSFR and higher metallicity. Similar observed gradients indicate that clumps survive outflows. The \exsitu clumps have stellar ages $0.5-3 {\rm Gyr}$ and sSFR $\sim 0.1-2 {\rm Gyr}^{-1}$, and they exhibit weaker gradients. Massive clumps of old stars at large radii are likely \exsitu mergers, though half of them share the disc rotation.

Wet disc contraction to galactic blue nuggets and quenching to red nuggets

We study the origin of high-redshift, compact, quenched spheroids (red nuggets) through the dissipative shrinkage of gaseous discs into compact star-forming systems (blue nuggets). The discs, fed by cold streams, undergo violent disc instability (VDI) that drives gas into the centre (along with mergers). The inflow is dissipative when its timescale is shorter than the star formation timescale. This implies a threshold of 0.28 in the cold-to-total mass ratio within the disc radius. For the typical gas fraction 0.5 at z=2, this threshold is traced back to a maximum spin parameter of 0.05, implying that half the star-forming galaxies contract to blue nuggets, while the rest form extended stellar discs. Thus, the surface density of blue galaxies is expected to be bimodal about 10^9 Msun/kpc^2, slightly increasing with mass. Blue nuggets are expected to be rare at low z when the gas fraction is low. The blue nuggets quench to red nuggets by complementary internal and external mechanisms. Internal quenching by a compact bulge, operating as a fast mode and especially at high z, may involve starbursts, stellar and AGN feedback, or Q-quenching. Quenching due to hot-medium haloes above 10^{12} Msun provides maintenance and a slower mode at low redshift. These predictions are confirmed in simulations and are consistent with observations at z=0-3.

Balance among gravitational instability, star formation and accretion determines the structure and evolution of disc galaxies

Over the past 10 Gyr, star-forming galaxies have changed dramatically, from clumpy and gas rich, to rather quiescent stellar-dominated disks with specific star formation rates lower by factors of a few tens. We present a general theoretical model for how this transition occurs, and what physical processes drive it, making use of 1D axisymmetric thin disk simulations with an improved version of the Gravitational Instability-Dominated Galaxy Evolution Tool (GIDGET) code. We show that at every radius galaxies tend to be in a slowly evolving equilibrium state wherein new accretion is balanced by star formation, galactic winds, and radial transport of gas through the disk by gravitational instability (GI) -driven torques. The gas surface density profile is determined by which of these terms are in balance at a given radius - direct accretion is balanced by star formation and galactic winds near galactic centers, and by transport at larger radii. We predict that galaxies undergo a smooth transition from a violent disk instability phase to secular evolution. This model provides a natural explanation for the high velocity dispersions and large clumps in z~2 galaxies, the growth and subsequent quenching of bulges, and features of the neutral gas profiles of local spiral galaxies.

THE LONG LIVES OF GIANT CLUMPS AND THE BIRTH OF OUTFLOWS IN GAS-RICH GALAXIES AT HIGH REDSHIFT

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

Toy models for galaxy formation versus simulations

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.

Newborn spheroids at high redshift: when and how did the dominant, old stars in today's massive galaxies form?

We study ∼330 massive (M* > 109.5 M⊙), newborn spheroidal galaxies (SGs) around the epoch of peak star formation (1 < z < 3) to explore the high-redshift origin of SGs and gain insight into when and how the old stellar populations that dominate today's Universe formed. The sample is drawn from the Hubble Space Telescope (HST)/WFC3 Early-Release Science programme, which provides deep 10-filter (0.2–1.7 μm) HST imaging over one-third of the GOODS-South field. We find that the star formation episodes that built our SGs likely peaked in the redshift range 2 < z < 5 (with a median of z ∼ 3) and have decay time-scales shorter than ∼1.5 Gyr. Starburst time-scales and ages show no trend with stellar mass in the range 109.5 < M* < 1010.5 M⊙. However, the time-scales show increased scatter towards lower values (<0.3 Gyr) for M* > 1010.5 M⊙, and an age trend becomes evident in this mass regime: SGs with M* > 1011.5 M⊙ are ∼2 Gyr older than their counterparts with M* < 1010.5 M⊙. Nevertheless, a smooth downsizing trend with galaxy mass is not observed, and the large scatter in starburst ages indicates that SGs are not a particularly coeval population. Around half of the blue SGs appear not to drive their star formation via major mergers, and those that have experienced a recent major merger show only modest enhancements (∼40 per cent) in their specific star formation rates. Our empirical study indicates that processes other than major mergers (e.g. violent disc instability driven by cold streams and/or minor mergers) likely play a dominant role in building SGs, and creating a significant fraction of the old stellar populations that dominate today's Universe.

Detectability of cold streams into high-redshift galaxies by absorption lines

Cold gas streaming along the dark-matter filaments of the cosmic web is predicted to be the major provider of resources for disc buildup, violent disk instability and star formation in massive galaxies in the early universe. We study to what extent these cold streams are traceable in the extended circum-galactic environment of galaxies via Ly alpha absorption and selected low ionisation metal absorption lines. We model the expected absorption signatures using high resolution zoom-in AMR cosmological simulations. In the postprocessing, we distinguish between self-shielded gas and unshielded gas. In the self-shielded gas, which is optically thick to Lyman continuum radiation, we assume pure collisional ionisation for species with an ionisation potential greater than 13.6 eV. In the optically thin, unshielded gas these species are also photoionised by the meta-galactic radiation. In addition to absorption of radiation from background quasars, we compute the absorption line profiles of radiation emitted by the galaxy at the centre of the same halo. We predict the strength of the absorption signal for individual galaxies without stacking. We find that the Ly alpha absorption profiles produced by the streams are consistent with observations of absorption and emission Ly alpha profiles in high redshift galaxies. Due to the low metallicities in the streams, and their low covering factors, the metal absorption features are weak and difficult to detect.

Evolution of violent gravitational disc instability in galaxies: late stabilization by transition from gas to stellar dominance

We address the cosmological evolution of violent gravitational instability in high-redshift, massive, star-forming galactic discs. To this aim, we integrate in time the equations of mass and energy conservation under self-regulated instability of a two-component disc of gas and stars. The disc is assumed to be continuously fed by cold gas at the average cosmological rate. The gas forms stars and is partly driven away by stellar feedback. The gas and stars flow inward through the disc to a central bulge due to torques that drive angular momentum outwards. The gravitational energy released by the mass inflow down the gravitational potential gradient drives the disc turbulence that maintains the disc unstable with a Toomre instability parameter Q∼ 1, compensating for the dissipative losses of the gas turbulence and raising the stellar velocity dispersion. We follow the velocity dispersion of stars and gas as they heat and cool, respectively, and search for disc ‘stabilization’, to be marked by a low gas velocity dispersion comparable to the speed of sound ∼10 km s−1. We vary the model parameters that characterize the accreted gas fraction, turbulence dissipation rate, star formation rate and stellar feedback. We find that as long as the gas input roughly follows the average cosmological rate, the disc instability is a robust phenomenon at high redshift till z∼ 1, driven by the high surface density and high gas fraction due to the intense cosmological accretion. For a broad range of model parameter values, the discs tend to ‘stabilize’ at z∼ 0–0.5 as they become dominated by hot stars. When the model parameters are pushed to extreme values, the discs may stabilize as early as z∼ 2, with the gas loss by strong outflows serving as the dominant stabilizing factor.

Investigating the Merger Origin of Early-type Galaxies using Ultra-deep Optical Images

AbstractThe mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry/wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra–deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the ATLAS3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec−2 in the g-band. Initial results for two galaxies are presented here.