BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Ceverino

Daniel Ceverino appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.IM

Catalog footprint

What is connected

43works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2024arXiv

A merger-driven scenario for clumpy galaxy formation in the epoch of reionization: Physical properties of clumps in the FirstLight simulation

Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub-)galactic clumps that are bright in [OIII] 5007$\mathrm{\mathring{A}}$ line with flux greater than $\sim 10^{-18} \, {\rm erg\, s^{-1}\, cm^{-2}}$, to be detected by JWST. For 62 simulated galaxies that have stellar masses of $(0.5-6) \times 10^{10} \, M_\odot$ at $z=5$, we find clumps in 1828 snapshots in the redshift range $z = 9.5-5.5$. The clumps are identified by the surface density of star formation rate. About one-tenth of the snapshots show the existence of clumpy systems with two or more components. Most of the clumps are formed by mergers and can be characterized by their ages; central clumps dominated by stellar populations older than 50 Myr, and off-centered clumps dominated by younger stellar populations with specific star formation rates of $\sim 50 \, {\rm Gyr^{-1}}$. The latter type of young clumps is formed from gas debris in the tidal tails of major mergers with baryonic mass ratios of $1 \leq q < 4$. The merger-induced clumps are short-lived, and merge within a dynamical time of several tens million years. The number density of the clumpy systems is estimated to be $\sim 10^{-5}\, {\rm cMpc^{-3}}$, which is large enough to be detected in recent JWST surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Emission Line Galaxies in the SHARDS Hubble Frontier Fields II: Limits on Lyman-Continuum Escape Fractions of Lensed Emission Line Galaxies at Redshifts 2 < z < 3.5

We present an investigation on escape fractions of UV photons from a unique sample of lensed low-mass emission line selected galaxies at z < 3.5 found in the SHARDS Hubble Frontier Fields medium-band survey. We have used this deep imaging survey to locate 42 relatively low-mass galaxies, down to $log(M_{*}/M_{\odot}) = 7$, between redshifts 2.4 < z < 3.5 which are candidate line emitters. Using deep multi-band Hubble UVIS imaging we investigate the flux of escaping ionizing photons from these systems, obtaining 1$σ$ upper limits of $f^{rel}_{esc}$ ~7% for individual galaxies, and < 2% for stacked data. We measure potential escaping Lyman-continuum flux for two low-mass line emitters with values at $f^{\rm rel}_{\rm esc} = 0.032^{+0.081}_{-0.009}$ and $f^{\rm rel}_{\rm esc} = 0.021^{+0.101}_{-0.006}$, both detected at the ~3.2$σ$ level. A detailed analysis of possible contamination reveals a < 0.1% probability that these detections result from line-of-sight contamination. The relatively low Lyman-continuum escape fraction limit, and the low fraction of systems detected, is an indication that low-mass line emitting galaxies may not be as important a source of reionization as hoped if these are analogs of reionization sources. We also investigate the structures of our galaxy sample, finding no evidence for a correlation of escape fraction with asymmetric structure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Clump Survival and Migration in VDI Galaxies: an Analytic Model versus Simulations and Observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

OVI Traces Photoionized Streams With Collisionally Ionized Boundaries in Cosmological Simulations of $z \sim 1$ Massive Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A mass threshold for galactic gas discs by spin flips

We predict, analytically and by simulations, that gas discs tend to survive only in haloes above a threshold mass $\sim 2 \times 10^{11} M_\odot$ (stellar mass $\sim 10^9 M_\odot$), with only a weak redshift dependence. At lower masses, the disc spins typically flip in less than an orbital time due to mergers associated with a change in the pattern of the feeding cosmic-web streams. This threshold arises from the halo merger rate when accounting for the mass dependence of the ratio of galactic baryons and halo mass. Above the threshold, wet compactions lead to massive central nuggets that allow the longevity of extended clumpy gas rings. Supernova feedback has a major role in disrupting discs below the critical mass, by driving the stellar-to-halo mass ratio that affects the merger rate, by stirring up turbulence and suppressing high-angular-momentum gas supply, and by confining major compactions to the critical mass. Our predictions seem consistent with current observed fractions of gas discs, to be explored by future observations that will resolve galaxies below $10^9 M_\odot$ at high redshifts, e.g. by JWST.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Origin of Star-Forming Rings around Massive Centres in Massive Galaxies at $z\!<\!4$

Using analytic modeling and simulations, we address the origin of an abundance of star-forming, clumpy, extended gas rings about massive central bodies in massive galaxies at $z \!<\! 4$. Rings form by high-angular-momentum streams and survive in galaxies of $M_{\rm star} \!>\! 10^{9.5-10} M_\odot$ where merger-driven spin flips and supernova feedback are ineffective. The rings survive after events of compaction to central nuggets. Ring longevity was unexpected based on inward mass transport driven by torques from violent disc instability. However, evaluating the torques from a tightly wound spiral structure, we find that the timescale for transport per orbital time is long and $\propto\! δ_{\rm d}^{-3}$, with $δ_{\rm d}$ the cold-to-total mass ratio interior to the ring. A long-lived ring forms when the ring transport is slower than its replenishment by accretion and the interior depletion by SFR, both valid for $δ_{\rm d} \!<\! 0.3$. The central mass that lowers $δ_{\rm d}$ is a compaction-driven bulge and/or dark matter, aided by the lower gas fraction at $z \!<\! 4$, provided that it is not too low. The ring is Toomre unstable for clump and star formation. The high-$z$ dynamic rings are not likely to arise form secular resonances or collisions. AGN feedback is not expected to affect the rings. Mock images of simulated rings through dust indicate qualitative consistency with observed rings about bulges in massive $z\!\sim\!0.5\!-\!3$ galaxies, in $H_α$ and deep HST imaging. ALMA mock images indicate that $z\!\sim\!0.5\!-\!1$ rings should be detectable. We quote expected observable properties of rings and their central nuggets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stellar Velocity Dispersion of a Massive Quenching Galaxy at z=4.01

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weak Evolution of the Mass-metallicity Relation at Cosmic Dawn in the FirstLight Simulations

Little is known about the mass-metallicity relation (MZR) in galaxies at cosmic dawn. Studying the first appearance of the MZR is one of the keys to understanding the formation and evolution of the first galaxies. In order to lay the groundwork for upcoming observational campaigns, we analyze 290 galaxies in halos spanning Mh = 10^9-10^11 Msun selected from the FirstLight (FL) cosmological zoom simulations to predict the MZR at z=5-8. Over this interval, the metallicity of FL galaxies with stellar mass Ms= 10^8 Msun declines by < 0.2 dex. This contrasts with the observed tendency for metallicities to increase at lower redshifts, and reflects weakly-evolving or even increasing gas fractions. We assess the use of the R3 strong-line diagnostic as a metallicity indicator, finding that it is informative for 12+log(O/H) < 8 but saturates to R3=3 at higher metallicities owing to a cancellation between enrichment and spectral softening. Nonetheless, campaigns with JWST should be able to detect a clear trend between R3 and stellar mass for Ms > 10^7.5 Msun. We caution that, at fixed metallicity, galaxies with higher specific star formation show higher R3 owing to their more intense radiation fields, indicating a potential for selection biases.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quiescent galaxies 1.5 billion years after the Big Bang and their progenitors

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

X-shooter Spectroscopy and HST Imaging of 15 Ultra Massive Quiescent Galaxies at $z\gtrsim2$

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Simulations at the Dwarf Scale: From Violent Dwarfs at Cosmic Dawn and Cosmic Noon to Quiet Discs today

Dwarf galaxies with stellar masses around 10^9 Msun can be explored at high and low redshifts and they give a glimpse of the different conditions of galaxy formation at different epochs. Using a large sample of about 300 zoom-in cosmological hydrodynamical simulations of galaxy formation I will briefly describe the formation of dwarfs at this mass scale at 3 different epochs: cosmic dawn (Ceverino, Klessen, Glover 2018), cosmic noon (Ceverino, Primack, Dekel 2015), and today (Ceverino et al. 2017). I will describe the FirstLight simulations of first galaxies at redshifts 5-15. These first dwarfs have extremely high star formation efficiencies due to high gas fractions and high gas accretion rates. These simulations will make predictions that will be tested for the first time with the James Webb Space Telescope (JWST). At cosmic noon, z = 2, galaxy formation is still a very violent and dynamic process. The VELA simulations have generated a set of dispersion-dominated dwarfs that show an elongated morphology due to their prolate dark-matter halos. Between z = 1 and 0, the AGORA simulation shows the formation of a low-mass disc due to slow gas accretion. The disc agrees with many local scaling relations, such as the stellar-mass-halo-mass and the baryonic Tully-Fisher relation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evolution of Density Profiles in High-z Galaxies: Compaction and Quenching Inside-Out

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evolution of Galaxy Shapes from Prolate to Oblate through Compaction Events

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gas inflow and metallicity drops in star-forming galaxies

Gas inflow feeds galaxies with low metallicity gas from the cosmic web, sustaining star formation across the Hubble time. We make a connection between these inflows and metallicity inhomogeneities in star-forming galaxies, by using synthetic narrow-band images of the Halpha emission line from zoom-in AMR cosmological simulations of galaxies with stellar masses of $M \simeq 10^9 $Msun at redshifts z=2-7. In $\sim$50\% of the cases at redshifts lower than 4, the gas inflow gives rise to star-forming, Halpha-bright, off-centre clumps. Most of these clumps have gas metallicities, weighted by Halpha luminosity, lower than the metallicity in the surrounding interstellar medium by $\sim$0.3 dex, consistent with observations of chemical inhomogeneities at high and low redshifts. Due to metal mixing by shear and turbulence, these metallicity drops are dissolved in a few disc dynamical times. Therefore, they can be considered as evidence for rapid gas accretion coming from cosmological inflow of pristine gas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Giant Clumps in Simulated High-z Galaxies: Properties, Evolution and Dependence on Feedback

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hubble Space Telescope Observations of Accretion-Induced Star Formation in the Tadpole Galaxy Kiso 5639

The tadpole galaxy Kiso 5639 has a slowly rotating disk with a drop in metallicity at its star-forming head, suggesting that star formation was triggered by the accretion of metal-poor gas. We present multi-wavelength HST WFC3 images of UV through I band plus Halpha to search for peripheral emission and determine the properties of various regions. The head has a mass in young stars of ~10^6 Mo and an ionization rate of 6.4x10^51 s^{-1}, equivalent to ~2100 O9-type stars. There are four older star-forming regions in the tail, and an underlying disk with a photometric age of ~1 Gyr. The mass distribution function of 61 star clusters is a power law with a slope of -1.73+-0.51. Fourteen young clusters in the head are more massive than 10^4 Mo, suggesting a clustering fraction of 30%-45%. Wispy filaments of Halpha emission and young stars extend away from the galaxy. Shells and holes in the head HII region could be from winds and supernovae. Gravity from the disk should limit the expansion of the HII region, although hot gas might escape through the holes. The star formation surface density determined from Halpha in the head is compared to that expected from likely pre-existing and accreted gas. Unless the surface density of the accreted gas is a factor of ~3 or more larger than what was in the galaxy before, the star formation rate has to exceed the usual Kennicutt-Schmidt rate by a factor of >5.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pathways to quiescence: SHARDS view on the Star Formation Histories of massive quiescent galaxies at 1.0 < z < 1.5

We present Star Formation Histories (SFHs) for a sample of 104 massive (stellar mass M > 10$^{10}$ M$_{\odot}$) quiescent galaxies (MQGs) at $z$=1.0-1.5 from the analysis of spectro-photometric data from the SHARDS and HST/WFC3 G102 and G141 surveys of the GOODS-N field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (Far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFR=SFR/Mass). The Spectral Energy Distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85%) by galaxies with young mass-weighted ages, tM $<$ 2 Gyr, short star formation timescales, $\langle τ\rangle$ $\sim$ 60-200 Myr, and masses log(M/M$_{\odot}$) $\sim$ 10.5. There is an older population (15%) with tM$=$2 - 4 Gyr, longer star formation timescales, $\langle τ\rangle$$\sim$ 400 Myr, and larger masses, log(M/M$_{\odot}$) $\sim$ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the Main Sequence (MS) of star-forming galaxies at $z$ $>$ 1.0, when our galaxies were 0.5-1.0 Gyr old. According to these SFHs, all the MQGs experienced a Luminous Infrared Galaxy (LIRG) phase that lasts for $\sim$ 500 Myr, and half of them an Ultra Luminous Infrared Galaxy (ULIRG) phase for $\sim$ 100 Myr. We find that the MQG population is almost assembled at $z$ $\sim$ 1, and continues evolving passively with few additions to the population.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Confinement of Star-Forming Galaxies into a Main Sequence through Episodes of Gas Compaction, Depletion, and Replenishment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Velocities of warm galactic outflows from synthetic Hα observations of star-forming galaxies

The velocity structure imprinted in the Hα emission line profiles contains valuable information about galactic outflows. Using a set of high-resolution zoom-in cosmological simulations of galaxies at z=2, we generate Hα emission line profiles, taking into account the temperature-dependent Hα emissivity, as well as dust extinction. The Hα line can be described as a sum of two gaussians, as typically done with observations. In general, its properties are in good agreement with those observed in local isolated galaxies with similar masses and star formation rates, assuming a spatially constant clumping factor of c=24. Blueshifted outflows are very common in the sample. They extend several kpc above the galaxy discs. They are also spread over the full extent of the discs. However, at small radii, the material with high velocities tends to remain confined within a thick disc, as part of galactic fountains or a turbulent medium, most probably due to the deeper gravitational potential at the galaxy centre.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<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)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 ($Σ_1$). Massive galaxies quench earlier, faster, and at a higher $Σ_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 $Σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Distribution of streaming rates into high-redshift galaxies

We study the accretion along streams from the cosmic web into high-redshift massive galaxies using three sets of AMR hydro-cosmological simulations. We find that the streams keep a roughly constant accretion rate as they penetrate into the halo centre. The mean accretion rate follows the mass and redshift dependence predicted for haloes by the EPS approximation, dM / dt is proportional to Mvir^{1.25} (1 + z)^{2.5}. The distribution of the accretion rates can well be described by a sum of two Gaussians, the primary corresponding to "smooth inflow" and the secondary to "mergers". The same functional form was already found for the distributions of specific star formation rates in observations. The mass fraction in the smooth component is 60 - 90 %, insensitive to redshift or halo mass. The simulations with strong feedback show clear signs of reaccretion due to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2 - 3 larger than that of the smooth component. The standard deviation of the accretion rate is 0.2 - 0.3 dex, showing no trend with mass or redshift. For the smooth component it is 0.12 - 0.24 dex.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diverse Structural Evolution at z > 1 in Cosmologically Simulated Galaxies

From mock Hubble Space Telescope images, we quantify non-parametric statistics of galaxy morphology, thereby predicting the emergence of relationships among stellar mass, star formation, and observed rest-frame optical structure at 1 < z < 3. We measure automated diagnostics of galaxy morphology in cosmological simulations of the formation of 22 central galaxies with 9.3 < log10 M_*/M_sun < 10.7. These high-spatial-resolution zoom-in calculations enable accurate modeling of the rest-frame UV and optical morphology. Even with small numbers of galaxies, we find that structural evolution is neither universal nor monotonic: galaxy interactions can trigger either bulge or disc formation, and optically bulge-dominated galaxies at this mass may not remain so forever. Simulated galaxies with M_* > 10^10 M_sun contain relatively more disc-dominated light profiles than those with lower mass, reflecting significant disc brightening in some haloes at 1 < z < 2. By this epoch, simulated galaxies with specific star formation rates below 10^-9.7 yr^-1 are more likely than normal star-formers to have a broader mix of structural types, especially at M_* > 10^10 M_sun. We analyze a cosmological major merger at z ~ 1.5 and find that the newly proposed MID morphology diagnostics trace later merger stages while G-M20 trace earlier ones. MID is sensitive also to clumpy star-forming discs. The observability time of typical MID-enhanced events in our simulation sample is less than 100 Myr. A larger sample of cosmological assembly histories may be required to calibrate such diagnostics in the face of their sensitivity to viewing angle, segmentation algorithm, and various phenomena such as clumpy star formation and minor mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of elongated galaxies with low masses at high redshift

We report the identification of elongated (triaxial or prolate) galaxies in cosmological simulations at $z\simeq2$. These are preferentially low-mass galaxies ($M_s \le 10^{9.5} \ M_\odot$), residing in dark-matter (DM) haloes with strongly elongated inner parts, a common feature of high-redshift DM haloes in the $Λ$CDM cosmology. Feedback slows formation of stars at the centres of these halos, so that a dominant and prolate DM distribution gives rise to galaxies elongated along the DM major axis. As galaxies grow in stellar mass, stars dominate the total mass within the galaxy half-mass radius, making stars and DM rounder and more oblate. A large population of elongated galaxies produces a very asymmetric distribution of projected axis ratios, as observed in high-z galaxy surveys. This indicates that the majority of the galaxies at high redshifts are not discs or spheroids but rather galaxies with elongated morphologies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 $λ_{\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λ_{\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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inflow velocities of cold flows streaming into massive galaxies at high redshifts

We study the velocities of the accretion along streams from the cosmic web into massive galaxies at high redshift with the help of three different suites of AMR hydrodynamical cosmological simulations. The results are compared to free-fall velocities and to the sound speeds of the hot ambient medium. The sound speed of the hot ambient medium is calculated using two different methods to determine the medium's temperature. We find that the simulated cold stream velocities are in violent disagreement with the corresponding free-fall profiles. The sound speed is a better albeit not always correct description of the cold flows' velocity. Using these calculations as a first order approximation for the gas inflow velocities vinflow = 0.9 vvir is given. We conclude from the hydrodynamical simulations as our main result that the velocity profiles for the cold streams are constant with radius. These constant inflow velocities seem to have a "parabola-like" dependency on the host halo mass in units of the virial velocity that peaks at Mvir = 10^{12} Msun and we also propose that the best fit functional form for the dependency of the inflow velocity on the redshift is a square root power law relation: vinflow is proportional to sqrt(z + 1) vvir.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-mass galaxy assembly in simulations: regulation of early star formation by radiation from massive stars

Despite recent success in forming realistic present-day galaxies, simulations still form the bulk of their stars earlier than observations indicate. We investigate the process of stellar mass assembly in low-mass field galaxies, a dwarf and a typical spiral, focusing on the effects of radiation from young stellar clusters on the star formation (SF) histories. We implement a novel model of SF with a deterministic low efficiency per free-fall time, as observed in molecular clouds. Stellar feedback is based on observations of star-forming regions, and includes radiation pressure from massive stars, photoheating in H II regions, supernovae and stellar winds. We find that stellar radiation has a strong effect on the formation of low-mass galaxies, especially at z > 1, where it efficiently suppresses SF by dispersing cold and dense gas, preventing runaway growth of the stellar component. This behaviour is evident in a variety of observations but had so far eluded analytical and numerical models without radiation feedback. Compared to supernovae alone, radiation feedback reduces the SF rate by a factor of ~100 at z < 2, yielding rising SF histories which reproduce recent observations of Local Group dwarfs. Stellar radiation also produces bulgeless spiral galaxies and may be responsible for excess thickening of the stellar disc. The galaxies also feature rotation curves and baryon fractions in excellent agreement with current data. Lastly, the dwarf galaxy shows a very slow reduction of the central dark matter density caused by radiation feedback over the last ~7 Gyr of cosmic evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-linear violent disc instability with high Toomre's Q in 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The co-planarity of satellite galaxies delivered by randomly aligned cold mode accretion streams

Recent observations have shown that the majority of the Andromeda galaxy's satellites are aligned in a thin plane. On the theoretical side it has been proposed that galaxies acquire their gas via cold streams. In addition, numerical simulations show that the same streams also deliver gas clumps which could potentially develop into satellite galaxies. Assuming that cold streams are a major source of satellite systems around galaxies we calculate the probabilities in two different models to find a certain fraction of satellites within a thin plane around the central galaxy of the host halo with and without having the same sense of rotation. Using simple geometrical considerations and adopting a random orientation of the streams we demonstrate that the vast thin disk of satellites detected around Andromeda can naturally be explained within this framework. In fact, without any satellite scattering, two streams would lead to too many satellites in the thin plane, compared with the observations. Three streams reproduce the observations very well. Natural implications from our model are that all massive galaxies should have a thin plane of satellites and that the satellites should naturally distribute themselves not only into a single plane but into several inclined ones. We estimate the effect of additional satellites accreted from random directions and find it to be of minor relevance for a mild inflow of satellites from random directions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiative feedback and the low efficiency of galaxy formation in low-mass haloes at high redshift

Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modeled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the $10^{4-4.5}$ K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of $\sim 1 \ {\rm M_\odot} \ {\rm yr}^{-1}$ in halos with masses about $10^{11} \ {M_\odot}$ at $z\simeq3$. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass halos contain only $\sim0.6\%$ of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Star Formation and Clumps in Cosmological Galaxy Simulations with Radiation Pressure Feedback

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Confronting Simulations of Optically Thick Gas in Massive Halos with Observations at z=2-3

Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M_vir~ 2x10^11 - 4x10^12 M_sun at z~2. We find that both the mass fraction of cool (T <= 3x10^4 K) gas and the covering fraction of optically-thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f_c~0.05 - 0.2, with significant scatter among halos. Our simulations of massive halos (M_vir >= 10^12 M_sun) underpredict the covering fraction of optically-thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M_vir >= 5x10^11 M_sun) hosting star-forming galaxies, the predicted covering factor agrees with observations, however current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale auto-correlation function of optically-thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z~2-3, as well as the characteristic size and covering factor of the CGM.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Absorption line systems in simulated galaxies fed by cold streams

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long narrow streams of merging galaxies and a smoother component of cold gas. We post-process seven high-resolution simulated galaxies with radiative transfer to study the absorption characteristics of the gas in galaxies and streams, in comparison with the statistics of observed absorption-line systems. We find that much of the stream gas is ionized by UV radiation from background and local stellar sources, but still optically thick (N_HI > 10^17 cm^-2) so that the streams appear as Lyman-limit systems (LLSs). At z>3, the fraction of neutral gas in streams becomes non-negligible, giving rise to damped Lyman-alpha (DLA) absorbers as well. The gas in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within one (two) virial radii, the covering factor of optically thick gas is <25% (10%) for LLSs and <5% (1%) for DLAs, slowly declining with time following the universal expansion. Nevertheless, galaxies and their cold streams in the studied mass range, M_vir = 10^10 - 10^12 Msun, account for >30% of the observed absorbers in the foreground of quasars, the rest possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in the streams is ~1% solar, much lower than in the galaxies. The simulated galaxies reproduce the Lyalpha-absorption equivalent widths observed around Lyman-break galaxies, but they severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise from outflows. We conclude that the observed metal-poor LLSs are likely detections of the predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations with more massive outflows in conjunction with the inflows may enable a clearer distinction between the signatures of the various gas modes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotational Support of Giant Clumps in High-z Disc Galaxies

We address the internal support against total free-fall collapse of the giant clumps that form by violent gravitational instability in high-z disc galaxies. Guidance is provided by an analytic model, where the proto-clumps are cut from a rotating disc and collapse to equilibrium while preserving angular momentum. This model predicts prograde clump rotation. This is confirmed in hydro-AMR zoom-in simulations of galaxies in a cosmological context. In most high-z clumps, the centrifugal force dominates the support, R=Vrot^2/Vcirc^2 > 0.5, where Vrot is the rotation velocity and Vcirc is the circular velocity. The clump spin indeed tends to be in the sense of the global disc angular momentum, but substantial tilts are frequent. Most clumps are in Jeans equilibrium, with the rest of the support provided by turbulence. Simulations of isolated gas-rich discs that resolve the clump substructure reveal that the cosmological simulations may overestimate R by ~30%, but the dominance of rotational support at high-z is not a resolution artifact. In turn, isolated gas-poor disc simulations produce at z=0 smaller gaseous non-rotating transient clouds, indicating that the difference in rotational support is associated with the fraction of cold baryons in the disc. In our current cosmological simulations, the clump rotation velocity is typically Vrot~100 km/s, but when beam smearing of \geq 0.1 arcsec is imposed, the rotation signal is reduced to a small gradient of \leq 30 km/s/kpc across the clump. The velocity dispersion in the simulated clumps is comparable to the disc dispersion so it is expected to leave only a marginal signal. Retrograde minor-merging galaxies could lead to massive clumps that do not show rotation.Testable predictions of the scenario as simulated are that the mean stellar age of the clumps, and the stellar fraction, are declining linearly with distance from the disc center.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-redshift clumpy discs and bulges in cosmological simulations

We analyze the first cosmological simulations that recover the fragmentation of high-redshift galactic discs driven by cold streams. The fragmentation is recovered owing to an AMR resolution better than 70 pc with cooling below 10^4 K. We study three typical star-forming galaxies in haloes of approx. 5 10^11 Msun at z=2.3, when they were not undergoing a major merger. The steady gas supply by cold streams leads to gravitationally unstable, turbulent discs, which fragment into giant clumps and transient features on a dynamical timescale. The disc clumps are not associated with dark-matter haloes. The clumpy discs are self-regulated by gravity in a marginaly unstable state. Clump migration and angular-momentum transfer on an orbital timescale help the growth of a central bulge with a mass comparable to the disc. The continuous gas input keeps the system of clumpy disc and bulge in a near "steady state", for several Gyr. The average star-formation rate, much of which occurs in the clumps, follows the gas accretion rate of approx. 45 Msun/yr. The simulated galaxies resemble in many ways the observed star-forming galaxies at high redshift. Their properties are consistent with the simple theoretical framework presented in Dekel, Sari & Ceverino (2009). In particular, a two-component analysis reveals that the simulated discs are indeed marginally unstable, and the time evolution confirms the robustness of the clumpy configuration in a cosmological steady state. By z=1 the simulated systems are stabilized by a dominant stellar spheroid, demonstrating the process of "morphological quenching" of star formation (Martig et al. 2009) . We demonstrate that the disc fragmentation is not a numerical artifact once the Jeans length is kept larger than 7 resolution elements, i.e. beyond the standard Truelove criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Low-mass galaxy formation in cosmological AMR simulations: the effects of varying the sub-grid physics parameters

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter halo (~7 times 10^10 Msun/h at redshift z=0). The simulations are run within a cosmological setting with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. In all of our simulations, an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed. We found that a non negligible fraction of the halo stars are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold densities nsf produce larger stellar effective radii Re, less peaked circular velocity curves V_c(R), and greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled by temporarily switching off radiative cooling in the star forming regions, Re increases (by a factor of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi-phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, measured by a stellar particle mass distribution biased toward larger values, increases the strength of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback is too strong, gas loss by galactic outflows -which are easier to produce in low-mass galaxies- interrupts SF, whose history becomes episodic. The simulations exhibit two important shortcomings: the baryon fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for redshifts ~ 0.4-1.

0 citations0 reviewsNo code linkedOpen access