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

Avishai Dekel

Avishai Dekel 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-phastro-ph.IM

Catalog footprint

What is connected

106works
4topics
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

106 published item(s)

preprint2022arXiv

On the origin of surprisingly cold gas discs in galaxies at high redshift

We address the puzzling observational indications for very "cold" galactic discs at redshifts $z \gtrsim 3$, an epoch when discs are expected to be highly perturbed. Using a high-resolution cosmological zoom-in simulation, we identify such a cold disc at $z\sim 3.5$, with a rotation velocity to velocity dispersion ratio of $v_ϕ/σ_r \simeq 5$ for the total gas. It forms as a result of a period of intense accretion of co-planar, co-rotating gas via cold cosmic-web streams. This thin disc survives for $\sim 5$ orbital periods, after which it is disrupted by mergers and counter-rotating streams, longer but consistent with our estimate that a galaxy of this mass ($M_\star\sim10^{10}\mathrm{M_\odot}$) typically survives merger-driven spin flips for $\sim 2-3$ orbital periods. We find that $v_ϕ/σ_r$ is highly sensitive to the tracer used to perform the kinematic analysis. While it is $v_ϕ/σ_r \simeq 3.5$ for atomic HI gas, it is $v_ϕ/σ_r \simeq 8$ for molecular CO and H$_2$. This reflects the confinement of molecular gas to cold, dense clouds that reside near the disc mid-plane, while the atomic gas is spread into a turbulent and more extended thicker disc. The proposed mechanisms is a theoretical proposal that has not been validated yet with proper statistical measurements and it remains unclear whether it occurs frequently enough to explain the multiple discoveries of cold gas disks in high-z galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of turbulent galactic discs: gravitational instability, feedback and accretion

We study the driving of turbulence in star-froming disc galaxies of different masses at different epochs, using an analytic "bathtub" model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc instability and clumpy accretion via streams. The transport rate is computed with two different formalisms, with similar results. To achieve the energy balance, the disc self-regulates either the mass fraction in clumps or the turbulent viscous torque parameter. In this version of the model, the efficiency by which the stream kinetic energy is converted into turbulence is a free parameter, $ξ_a$. We find that the contributions of the three energy sources are in the same ball park, within a factor of $\sim\!2$ in all discs at all times. In haloes that evolve to a mass $\leq 10^{12}\,\Msun$ by $z=0$ ($\leq 10^{11.5}\,\Msun$ at $z\!\sim\!2$), feedback is the main driver throughout their lifetimes. Above this mass, the main driver is either transport or accretion for very low or very high values of $ξ_a$, respectively. For an assumed $ξ_a(t)$ that declines in time, galaxies in halos with present-day mass $>\!10^{12}$ M$_\odot$ make a transition from accretion to transport dominance at intermediate redshifts, $z\! \sim\!3$, when their mass was $\geq\!10^{11.5}\,\Msun$. The predicted relation between star-formation rate and gas velocity dispersion is consistent with observations.

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

Mock Lightcones and Theory Friendly Catalogs for the CANDELS Survey

We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass the footprints of observations associated with the five CANDELS fields. In order to allow field-to-field variance to be explored, we have created eight realizations of each field. In this paper, we present comparisons with observable global galaxy properties, including counts in observed frame bands, luminosity functions, color-magnitude distributions and color-color distributions. We additionally present comparisons with physical galaxy parameters derived from SED fitting for the CANDELS observations, such as stellar masses and star formation rates. We find relatively good agreement between the model predictions and CANDELS observations for luminosity and stellar mass functions. We find poorer agreement for colors and star formation rate distributions. All of the mock lightcones as well as curated "theory friendly" versions of the observational CANDELS catalogs are made available through a web-based data hub.

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

Instability of Supersonic Cold Streams Feeding Galaxies IV: Survival of Radiatively Cooling Streams

We study the effects of Kelvin Helmholtz Instability (KHI) on the cold streams that feed massive halos at high redshift, generalizing our earlier results to include the effects of radiative cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, $t_{\rm shear}$. This can be converted into a critical stream radius, $R_{\rm s,crit}$, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,mix}$. If $R_{\rm s}<R_{\rm s,crit}$, the non-linear evolution proceeds similarly to the non-radiative case studied by Mandelker et al. 2019a. If $R_{\rm s}>R_{\rm s,crit}$, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses onto the stream, and can increase the mass of cold gas by a factor of $\sim 3$ within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T<5\times 10^4{\rm K}$, suggesting much of it will be emitted in Ly$α$.

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

Quenching as a Contest between Galaxy Halos and their Central Black Holes

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The relationship between fine galaxy stellar morphology and star formation activity in cosmological simulations: a deep learning view

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SFR-radius connection: data and implications for wind strength and halo concentration

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A model for core formation in dark matter haloes and ultra diffuse galaxies by outflow episodes

We present a simple model for the response of a dissipationless spherical system to an instantaneous mass change at its center, describing the formation of flat cores in dark matter haloes and ultra-diffuse galaxies (UDGs) from feedback-driven outflow episodes in a specific mass range. This model generalizes an earlier simplified analysis of an isolated shell into a system with continuous density, velocity and potential profiles. The response is divided into an instantaneous change of potential at constant velocities due to a given mass loss or gain, followed by energy-conserving relaxation to a new Jeans equilibrium. The halo profile is modeled by a two-parameter function with a variable inner slope and an analytic potential profile (Dekel et al. 2017), which enables determining the associated kinetic energy at equilibrium. The model is tested against NIHAO cosmological zoom-in simulations, where it successfully predicts the evolution of the inner dark-matter profile between successive snapshots in about 75% of the cases, failing mainly in merger situations. This model provides a simple understanding of the formation of dark-matter halo cores and UDGs by supernova-driven outflows, and a useful analytic tool for studying such processes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Unveiling Sizes of Compact AGN Hosts with ALMA

We present rest-frame far-infrared (FIR) and optical size measurements of AGN hosts and star-forming galaxies in the COSMOS field, enabled by high-resolution ALMA/1 mm (0.1 arcsec - 0.4 arcsec) and HST/F814W imaging (~ 0.1 arcsec). Our sample includes 27 galaxies at z<2.5, classified as infrared-selected AGN (3 sources), X-ray selected AGN (4 sources), and non-AGN star-forming galaxies (20 sources), for which high-resolution Band 6/7 ALMA images are available at 1 mm from our own observing program as well as archival observations. The sizes and SFR surface densities measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies are more compact than non-AGN star-forming galaxies at similar redshift and stellar mass. This result suggests that the obscured accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could be associated with enhanced central star formation before a subsequent quenching driving the formation of compact passive galaxies. Moreover, most of the detected and stacked rest-frame FIR sizes of AGNs in our sample are similar or more compact than their rest-frame optical sizes, which is consistent with recent results of ALMA detected sources. This might be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those observed in star-forming galaxies observed with ALMA.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enhanced Momentum Feedback from Clustered Supernovae

Young stars typically form in star clusters, so the supernovae (SNe) they produce are clustered in space and time. This clustering of SNe may alter the momentum per SN deposited in the interstellar medium (ISM) by affecting the local ISM density, which in turn affects the cooling rate. We study the effect of multiple SNe using idealized 1D hydrodynamic simulations which explore a large parameter space of the number of SNe, and the background gas density and metallicity. The results are provided as a table and an analytic fitting formula. We find that for clusters with up to ~100 SNe the asymptotic momentum scales super-linearly with the number of SNe, resulting in a momentum per SN that can be an order of magnitude larger than for a single SN, with a maximum efficiency for clusters with 10-100 SNe. We argue that additional physical processes not included in our simulations -- self-gravity, breakout from a galactic disk, and galactic shear -- can slightly reduce the momentum enhancement from clustering, but the average momentum per SN still remains a factor of 4 larger than the isolated SN value when averaged over a realistic cluster mass function for a star-forming galaxy. We conclude with a discussion of the possible role of mixing between hot and cold gas, induced by multi-dimensional instabilities or preexisting density variations, as a limiting factor in the buildup of momentum by clustered SNe, and suggest future numerical experiments to explore these effects.

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

Evolution of Intrinsic Scatter in the SFR-Stellar Mass Correlation at 0.5<z<3

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

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

Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MUSE searches for galaxies near very metal-poor gas clouds at z~3: new constraints for cold accretion models

We report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Lyα emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Lyα emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Lyα luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Lyα nebula at the quasar redshift and three Lyα emitters reminiscent of a "dark galaxy" population.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Ω_b/Ω_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

NIHAO XI: Formation of Ultra-Diffuse Galaxies by outflows

We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average Sérsic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stellar Mass--Gas-phase Metallicity Relation at $0.5\leq z\leq0.7$: A Power Law with Increasing Scatter toward the Low-mass Regime

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Suppression of star formation in dwarf galaxies by grain photoelectric feedback

Photoelectric heating has long been recognized as the primary source of heating for the neutral interstellar medium. Simulations of spiral galaxies found some indication that photoelectric heating could suppress star formation. However, simulations that include photoelectric heating have typically found that it has little effect on the rate of star formation in either spiral galaxies or dwarfs suggesting that supernovae and not photoelectric heating are responsible for setting the star formation law in galaxies. This result is in tension with recent work indicating that a star formation law that depends on galaxy metallicity, as expected for photoelectric heating but not for supernovae, reproduces the present-day galaxy population better than a metallicity-independent one. Here we report a series of simulations of dwarf galaxies, where the effects of both photoelectric heating and supernovae are expected to be strongest. We simultaneously include space- and time-dependent photoelectric heating, and we resolve the Sedov phase of every supernova blast wave, allowing us to make a direct measurement of the relative importance of momentum injection by supernovae and dust heating by far ultraviolet (FUV) photons in suppressing star formation. We find that supernovae are unable to account for the long observed gas depletion times in dwarf galaxies. Instead, ordinary photoelectric heating is the dominant means by which dwarf galaxies regulate their star formation rate at any given time, suppressing the star formation rate by more than an order of magnitude relative to simulations with only supernovae.

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

The Formation of Bulges, Discs and Two Component Galaxies in the CANDELS Survey at z < 3

We examine a sample of 1495 galaxies in the CANDELS fields to determine the evolution of two component galaxies, including bulges and discs, within massive galaxies at the epoch 1 < z < 3 when the Hubble sequence forms. We fit all of our galaxies' light profiles with a single Sérsic fit, as well as with a combination of exponential and Sérsic profiles. The latter is done in order to describe a galaxy with an inner and an outer component, or bulge and disc component. We develop and use three classification methods (visual, F-test and the RFF) to separate our sample into 1-component galaxies (disc/spheroids-like galaxies) and 2-component galaxies (galaxies formed by an 'inner part' or bulge and an 'outer part' or disc). We then compare the results from using these three different ways to classify our galaxies. We find that the fraction of galaxies selected as 2-component galaxies increases on average 50 per cent from the lowest mass bin to the most massive galaxies, and decreases with redshift by a factor of four from z = 1 to z = 3. We find that single Sérsic 'disc-like' galaxies have the highest relative number densities at all redshifts, and that 2-component galaxies have the greatest increase and become at par with Sérsic discs by z = 1. We also find that the systems we classify as 2-component galaxies have an increase in the sizes of their outer components, or 'discs', by about a factor of three from z = 3 to z = 1.5, while the inner components or 'bulges' stay roughly the same size. This suggests that these systems are growing from the inside out, whilst the bulges or protobulges are in place early in the history of these galaxies. This is also seen to a lesser degree in the growth of single 'disc-like' galaxies vs. 'spheroid-like' galaxies over the same epoch.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Relationship Between Star-formation Activity and Galaxy Structural Properties in CANDELS and a Semi-analytic Model

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS.

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

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

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

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

Quenching and Morphological Transformation in Semi-Analytic Models and CANDELS

We examine the spheroid growth and star formation quenching experienced by galaxies from z~3 to the present by studying the evolution with redshift of the quiescent and spheroid-dominated fractions of galaxies from the CANDELS and GAMA surveys. We compare the observed fractions with predictions from a semi-analytic model which includes prescriptions for bulge growth and AGN feedback due to mergers and disk instabilities. We facilitate direct morphological comparison by converting our model bulge-to-total stellar mass ratios to Sersic indices. We then subdivide our population into the four quadrants of the sSFR-Sersic index plane and study the buildup of each of these subpopulations. We find that the fraction of star forming disks declines steadily, while the fraction of quiescent spheroids builds up over cosmic time. The fractions of star forming spheroids and quiescent disks are both non-negligible, and stay nearly constant over the period we have studied, at about 10% and 15-20% respectively. Our model is qualitatively successful at reproducing the evolution of the two "main" populations (star forming disk-dominated galaxies and quiescent spheroid-dominated galaxies), and approximately reproduces the relative fractions of all four types, but predicts a stronger decline in star forming spheroids, and increase in quiescent disks, than seen in the observations. A model with an additional channel for bulge growth via disk instabilities agrees better overall with the observations than a model in which bulges may grow only through mergers. We study evolutionary tracks of some individual galaxies as they experience morphological transformation and quenching, and examine the importance of different physical drivers of this transformation (major and minor mergers and disk instabilities). We find that complex histories with multiple transformative events are the norm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Satellite Quenching and Galactic Conformity at 0.3 < z < 2.5

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3σ$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two Conditions for Galaxy Quenching: Compact Centres and Massive Haloes

We investigate the roles of two classes of quenching mechanisms for central and satellite galaxies in the SDSS ($z<0.075$): those involving the halo and those involving the formation of a compact centre. For central galaxies with inner compactness $Σ_{\rm 1kpc} \sim 10^{9-9.4}M_{\odot} {\rm kpc}^{-2}$, the quenched fraction $f_{q}$ is strongly correlated with $Σ_{\rm 1kpc}$ with only weak halo mass $M_{\rm h}$ dependence. However, at higher and lower $Σ_{\rm 1kpc}$, sSFR is a strong function of $M_{\rm h}$ and mostly independent of $Σ_{\rm 1kpc}$. In other words, $Σ_{\rm 1kpc} \sim 10^{9-9.4} M_{\odot} {\rm kpc}^{-2}$ divides galaxies into those with high sSFR below and low sSFR above this range. In both the upper and lower regimes, increasing $M_{\rm h}$ shifts the entire sSFR distribtuion to lower sSFR without a qualitative change in shape. This is true even at fixed $M_{*}$, but varying $M_{*}$ at fixed $M_{\rm h}$ adds no quenching information. Most of the quenched centrals with $M_{\rm h} > 10^{11.8}M_{\odot}$ are dense ($Σ_{\rm 1kpc} > 10^{9}~ M_{\odot} {\rm kpc}^{-2}$), suggesting compaction-related quenching maintained by halo-related quenching. However, 21% are diffuse, indicating only halo quenching. For satellite galaxies in the outskirts of halos, quenching is a strong function of compactness and a weak function of host $M_{\rm h}$. In the inner halo, $M_{\rm h}$ dominates quenching, with $\sim 90\%$ of the satellites being quenched once $M_{\rm h} > 10^{13}M_{\odot}$. This regional effect is greatest for the least massive satellites. As demonstrated via semi-analytic modelling with simple prescriptions for quenching, the observed correlations can be explained if quenching due to central compactness is rapid while quenching due to halo mass is slow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An analytic solution for the minimal bathtub toy model: challenges in the star-formation history of high-z galaxies

We study the minimal ``bathtub" toy model as an analytic tool for capturing key processes of galaxy evolution and identifying robust successes and challenges in reproducing observations at high redshift. The source and sink terms of the continuity equations for gas and stars are expressed in simple terms from first principles. The assumed dependence of star-formation rate (SFR) on gas mass self-regulates the system into a unique asymptotic behavior, which is approximated by an analytic quasi-steady-state solution (QSS). We address the validity of the QSS at different epochs independent of earlier conditions. At high z, where the accretion is assumed to consist of gas only, the specific SFR is robustly predicted to be sSFR = [(1+z)/3]^{5/2} Gyr^{-1}, slightly higher than the cosmological specific accretion rate, in agreement with observations at z=3-8. The gas fraction is expected to decline slowly, and the observations constrain the SFR efficiency per dynamical time to epsilon=0.02. The stellar-to-virial mass ratio f_sv is predicted to be constant in time, and the observed value requires an outflow mass-lading factor of eta=1-3, depending on the penetration efficiency of gas into the galaxy. However, at z=2, where stars are also accreted through mergers, the simplest model has an apparent difficulty in matching observations. The model that maximizes the sSFR, with the outflows fully recycled, falls short by a factor 3 in sSFR, and overestimates f_sv. With strong outflows, the model can reproduce the observed f_sv but at the expense of underestimating the sSFR by an order of magnitude. We discuss potential remedies including a bias due to the exclusion of quenched galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bulge Growth and Quenching since z = 2.5 in CANDELS/3D-HST

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deconstructing the Galaxy Stellar Mass Function with UKIDSS and CANDELS: the Impact of Colour, Structure and Environment

We combine photometry from the UDS, and CANDELS UDS and CANDELS GOODS-S surveys to construct the galaxy stellar mass function probing both the low and high mass end accurately in the redshift range 0.3<z<3. The advantages of using a homogeneous concatenation of these datasets include meaningful measures of environment in the UDS, due to its large area (0.88 deg^2), and the high resolution deep imaging in CANDELS (H_160 > 26.0), affording us robust measures of structural parameters. We construct stellar mass functions for the entire sample as parameterised by the Schechter function, and find that there is a decline in the values of phi and of alpha with higher redshifts, and a nearly constant M* up to z~3. We divide the galaxy stellar mass function by colour, structure, and environment and explore the links between environmental over-density, morphology, and the quenching of star formation. We find that a double Schechter function describes galaxies with high Sersic index (n>2.5), similar to galaxies which are red or passive. The low-mass end of the n>2.5 stellar mass function is dominated by blue galaxies, whereas the high-mass end is dominated by red galaxies. This hints that possible links between morphological evolution and star formation quenching are only present in high-mass galaxies. This is turn suggests that there are strong mass dependent quenching mechanisms. In addition, we find that the number density of high mass systems is elevated in dense environments, suggesting that an environmental process is building up massive galaxies quicker in over densities than in lower densities.

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

No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Observations of environmental quenching in groups in the 11 Gyr since z=2.5: different quenching for central and satellite galaxies

We present direct observational evidence for star formation quenching in galaxy groups in the redshift range 0<z<2.5. We utilize a large sample of nearly 6000 groups, selected by fixed cumulative number density from three photometric catalogs, to follow the evolving quiescent fractions of central and satellite galaxies over roughly 11 Gyr. At z~0, central galaxies in our sample range in stellar mass from Milky Way/M31 analogs (M=6.5x10^10 M\solar) to nearby massive ellipticals (M=1.5x10^11 M\solar). Satellite galaxies in the same groups reach masses as low as twice that of the Large Magellanic Cloud (M=6.5x10^9 M\solar). Using statistical background subtraction, we measure the average rest-frame colors of galaxies in our groups and calculate the evolving quiescent fractions of centrals and satellites over seven redshift bins. Our analysis shows clear evidence for star formation quenching in group halos, with a different quenching onset for centrals and their satellite galaxies. Using halo mass estimates for our central galaxies, we find that star formation shuts off in centrals when typical halo masses reach between 10^12 and 10^13 M\solar, consistent with predictions from the halo quenching model. In contrast, satellite galaxies in the same groups most likely undergo quenching by environmental processes, whose onset is delayed with respect to their central galaxy. Although star formation is suppressed in all galaxies over time, the processes that govern quenching are different for centrals and satellites. While mass plays an important role in determining the star formation activity of central galaxies, quenching in satellite galaxies is dominated by the environment in which they reside.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Properties of Submillimeter Galaxies in the CANDELS GOODS-S Field

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

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 distribution of satellites around massive galaxies at 1<z<3 in ZFOURGE/CANDELS: dependence on star formation activity

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$σ$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The interstellar medium and feedback in the progenitors of the compact passive galaxies at z~2

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The mass evolution of the first galaxies: stellar mass functions and star formation rates at $4 < z < 7$ in the CANDELS GOODS-South field

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($α\approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} ρ_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($β$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

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

A Link Between Star Formation Quenching and Inner Stellar Mass Density in SDSS Central Galaxies

We study the correlation between galaxy structure and the quenching of star formation using a sample of SDSS central galaxies with stellar masses 9.75< log M_*/M_sun<11.25 and redshifts z<0.075. GALEX UV data are used to cleanly divide the sample into star-forming and quenched galaxies, and to identify galaxies in transition (the green valley). Despite a stark difference in visual appearance between blue and red galaxies, their average radial stellar mass density profiles are remarkably similar (especially in the outer regions) at fixed mass. The inner stellar mass surface density within a radius of 1 kpc, Σ_1, is used to quantify the growth of the bulge as galaxies evolve. When galaxies are divided into narrow mass bins, their distribution in the color-Σ_1 plane at fixed mass forms plausible evolutionary tracks. Σ_1 seems to grow as galaxies evolve through the blue cloud, and once it crosses a threshold value, galaxies are seen to quench at fixed Σ_1. The Σ_1 threshold for quenching grows with stellar mass, Σ_1 ~ M_*^{0.64}. However, the existence of some star-forming galaxies above the threshold Σ_1 implies that a dense bulge is necessary but not sufficient to quench a galaxy fully. This would be consistent with a two-step quenching process in which gas within a galaxy is removed or stabilized against star formation by bulge-driven processes (such as a starburst, AGN feedback, or morphological quenching), whereas external gas accretion is suppressed by separate halo-driven processes (such as halo gas shock heating). Quenching thus depends on an interplay between the inner structure of a galaxy and its surrounding dark matter halo, and lack of perfect synchrony between the two could produce the observed scatter in color vs. Σ_1. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Balance among gravitational instability, star formation, and accretion determines the structure and evolution of disk 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CANDELS Observations of the Environmental Dependence of the Color-Mass-Morphology Relation at z = 1.6

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and Sérsic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller Sérsic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low Sérsic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CANDELS: The correlation between galaxy morphology and star formation activity at z~2

We discuss the state of the assembly of the Hubble Sequence in the mix of bright galaxies at redshift 1.4< z \le 2.5 with a large sample of 1,671 galaxies down to H_{AB}~26, selected from the HST/ACS and WFC3 images of the GOODS--South field obtained as part of the GOODS and CANDELS observations. We investigate the relationship between the star formation properties and morphology using various parametric diagnostics, such as the Sersic light profile, Gini (G), M_{20}, Concentration (C), Asymmetry (A) and multiplicity parameters. Our sample clearly separates into massive, red and passive galaxies versus less massive, blue and star forming ones, and this dichotomy correlates very well with the galaxies' morphological properties. Star--forming galaxies show a broad variety of morphological features, including clumpy structures and bulges mixed with faint low surface brightness features, generally characterized by disky-type light profiles. Passively evolving galaxies, on the other hand, very often have compact light distribution and morphology typical of today's spheroidal systems. We also find that artificially redshifted local galaxies have a similar distribution with z~2galaxies in a G-M_{20} plane. Visual inspection between the rest-frame optical and UV images show that there is a generally weak morphological k-correction for galaxies at z~2, but the comparison with non-parametric measures show that galaxies in the rest-frame UV are somewhat clumpier than rest-frame optical. Similar general trends are observed in the local universe among massive galaxies, suggesting that the backbone of the Hubble sequence was already in place at z~2.

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
preprint2013arXiv

On the origin of the fundamental metallicity relation and the scatter in galaxy scaling relations

We present a simple toy model to understand what sets the scatter in star formation and metallicity of galaxies at fixed mass. The scatter ultimately arises from the intrinsic scatter in the accretion rate, but may be substantially reduced depending on the timescale on which the accretion varies compared to the timescale on which the galaxy loses gas mass. This model naturally produces an anti-correlation between star formation and metallicity at a fixed mass, the basis of the fundamental metallicity relation. We show that observational constraints on the scatter in galaxy scaling relations can be translated into constraints on the galaxy-to-galaxy variation in the mass loading factor, and the timescales and magnitude of stochastic accretion onto star-forming galaxies. We find a remarkably small scatter in the mass loading factor, < 0.1 dex, and that the scatter in accretion rates is smaller than expected from N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Steady Outflows in Giant Clumps of High-z Disk Galaxies During Migration and Growth by Accretion

We predict the evolution of giant clumps undergoing star-driven outflows in high-z gravitationally unstable disk galaxies. We find that the mass loss is expected to occur through a steady wind over many tens of free-fall times (t_ff ~ 10 Myr) rather than by an explosive disruption in one or a few t_ff. Our analysis is based on the finding from simulations that radiation trapping is negligible because it destabilizes the wind (Krumholz & Thompson 2012, 2013). Each photon can therefore contribute to the wind momentum only once, so the radiative force is limited to L/c. When combining radiation, protostellar and main-sequence winds, and supernovae, we estimate the total direct injection rate of momentum into the outflow to be 2.5 L/c. The adiabatic phase of supernovae and main-sequence winds can double this rate. The resulting outflow mass-loading factor is of order unity, and if the clumps were to deplete their gas the timescale would have been a few disk orbital times, to end with half the original clump mass in stars. However, the clump migration time to the disk center is on the order of an orbital time, about 250 Myr, so the clumps are expected to complete their migration prior to depletion. Furthermore, the clumps are expected to double their mass in a disk orbital time by accretion from the disk and clump-clump mergers, so their mass actually grows in time and with decreasing radius. From the 6-7 giant clumps with observed outflows, 5 are consistent with these predictions, and one has a much higher mass-loading factor and momentum injection rate. The latter either indicates that the estimated outflow is an overestimate (within the 1-sigma error), that the SFR has dropped since the time when the outflow was launched, or that the driving mechanism is different, e.g. supernova feedback in a cavity generated by the other feedbacks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structural Evolution of Early-type Galaxies to z=2.5 in CANDELS

Projected axis ratio measurements of 880 early-type galaxies at redshifts 1<z<2.5 selected from CANDELS are used to reconstruct and model their intrinsic shapes. The sample is selected on the basis of multiple rest-frame colors to reflect low star-formation activity. We demonstrate that these galaxies as an ensemble are dust-poor and transparent and therefore likely have smooth light profiles, similar to visually classified early-type galaxies. Similar to their present-day counterparts, the z>1 early-type galaxies show a variety of intrinsic shapes; even at a fixed mass, the projected axis ratio distributions cannot be explained by the random projection of a set of galaxies with very similar intrinsic shapes. However, a two-population model for the intrinsic shapes, consisting of a triaxial, fairly round population, combined with a flat (c/a~0.3) oblate population, adequately describes the projected axis ratio distributions of both present-day and z>1 early-type galaxies. We find that the proportion of oblate versus triaxial galaxies depends both on the galaxies' stellar mass, and - at a given mass - on redshift. For present-day and z<1 early-type galaxies the oblate fraction strongly depends on galaxy mass. At z>1 this trend is much weaker over the mass range explored here (10^10<M*/M_sun<10^11), because the oblate fraction among massive (M*~10^11 M_sun) was much higher in the past: 0.59+-0.10 at z>1, compared to 0.20+-0.02 at z~0.1. In contrast, the oblate fraction among low-mass early-type galaxies (log(M*/M_sun)<10.5) increased toward the present, from 0.38+-0.11 at z>1 to 0.72+-0.06 at z=0. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Progenitors of the Compact Early-Type Galaxies at High-Redshift

We use GOODS and CANDELS images to identify progenitors of massive (log M > 10 Msun) compact "early-type" galaxies (ETGs) at z~1.6. Since merging and accretion increase the size of the stellar component of galaxies, if the progenitors are among known star-forming galaxies, these must be compact themselves. We select candidate progenitors among compact Lyman-break galaxies at z~3 based on their mass, SFR and central stellar density and find that these account for a large fraction of, and possibly all, compact ETGs at z~1.6. We find that the average far-UV SED of the candidates is redder than that of the non-candidates, but the optical and mid-IR SED are the same, implying that the redder UV of the candidates is inconsistent with larger dust obscuration, and consistent with more evolved (aging) star-formation. This is in line with other evidence that compactness is a sensitive predictor of passivity among high-redshift massive galaxies. We also find that the light distribution of both the compact ETGs and their candidate progenitors does not show any extended "halos" surrounding the compact "core", both in individual images and in stacks. We argue that this is generally inconsistent with the morphology of merger remnants, even if gas-rich, as predicted by N-body simulations. This suggests that the compact ETGs formed via highly dissipative, mostly gaseous accretion of units whose stellar components are very small and undetected in the HST images, with their stellar mass assembling in-situ, and that they have not experienced any major merging until the epoch of observations at z~1.6.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Redshift and Mass Dependence on the Formation of The Hubble Sequence at z>1 from CANDELS/UDS

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The star-forming progenitors of massive red galaxies

The link between massive red galaxies in the local Universe and star-forming galaxies at high redshift is investigated with a semi-analytic model that has proven successful in many ways, e.g. explaining the galaxy colour-magnitude bimodality and the stellar mass-age relation for red-sequence galaxies. The model is used to explore the processes that drive star formation in different types of galaxies as a function of stellar mass and redshift. We find that most z=2-4 star-forming galaxies with M_*>10^10 M_Sun evolve into red-sequence galaxies. Also, most of the massive galaxies on the red-sequence today have passed through a phase of intense star formation at z>2. Specifically, ~ 90% of today's red galaxies with M_*>10^11 M_Sun were fed during this phase by cold streams including minor mergers. Gas-rich major mergers are rare and the effects of merger-driven starbursts are ephemeral. On the other hand, major mergers are important in powering the most extreme starbursts. Gas-rich mergers also explain the tail of intermediate-mass red galaxies that form relatively late, after the epoch of peak star formation. In two thirds of the currently red galaxies that had an intense star formation event at $z<1$, this event was triggered by a merger.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A diversity of progenitors and histories for isolated spiral galaxies

We analyze a suite of 33 cosmological simulations of the evolution of Milky Way-mass galaxies in low-density environments. Our sample spans a broad range of Hubble types at z=0, from nearly bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is typically in place by z=1, we find no significant correlation between the morphology at z=1 and at z=0. The z=1 progenitors of disk galaxies span a range of morphologies, including smooth disks, unstable disks, interacting galaxies and bulge-dominated systems. By z=0.5, spiral arms and bars are largely in place and the progenitor morphology is correlated with the final morphology. We next focus on late-type galaxies with a bulge-to-total ratio B/T<0.3 at z=0. These show a correlation between B/T at z=0 and the mass ratio of the largest merger at z<2, as well as with the gas accretion rate at z>1. We find that the galaxies with the lowest B/T tend to have a quiet baryon input history, with no major mergers at z<2, and with a low and constant gas accretion rate that keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n<2. The galaxies with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability rather than active mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Observed Link between Active Galactic Nuclei and Violent Disk Instabilities in High-Redshift Galaxies

We provide evidence for a correlation between the presence of giant clumps and the occurrence of active galactic nuclei (AGN) in disk galaxies. Giant clumps of 10^8-9 Msun arise from violent gravitational instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks in our sample have a high probability of containing AGN. Their [OIII] emission line is strongly excited, inconsistent with low-metallicity star formation alone. [NeIII] excitation is also higher. Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals an X-ray excess in clumpy galaxies, which cannot be solely due to star formation and confirms the presence of AGN. The clumpy galaxies in our intermediate-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, being in particular on the Main Sequence of star formation. This suggests that our findings apply to the physically-similar and numerous gas-rich unstable disks at z>1. Using the observed [OIII] and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical bolometric luminosities of the order of a few 10^43 erg/s, BH growth rates ~10^-2 Msun/yr, and that these AGN are substantially obscured in X-rays. This moderate-luminosity mode could be sufficient to provide a large fraction of today's BH mass over a couple of Gyr given that our observations suggest a high duty cycle (>10%), accretion bursts with higher luminosities being possible over shorter phases. The observed evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing of star formation and of BH growth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Baryons Matter: Why Luminous Satellite Galaxies Have Reduced Central Masses

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CANDELS: Correlations of SEDs and Morphologies with Star-formation Status for Massive Galaxies at z ~ 2

We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CANDELS: The progenitors of compact quiescent galaxies at z~2

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Co-planar streams, pancakes, and angular-momentum exchange in high-z disc galaxies

We study the feeding of massive galaxies at high redshift through streams from the cosmic web using the Mare Nostrum hydro-cosmological simulation. Our statistical sample consists of 350 dark-matter haloes of ~10^12 Msun at z = 2.5. We find that ~70% of the influx into the virial radius Rv is in narrow streams covering 10% of the virial shell. On average 64% of the stream influx is in one stream, and 95% is in three dominant streams. The streams that feed a massive halo tend to lie in a plane that extends from half to a few Rv, hereafter "the stream plane" (SP). The streams are typically embedded in a thin sheet of low-entropy gas, a Zel'dovich pancake, which carries ~20% of the influx into Rv. The filaments-in-a-plane configuration about the massive haloes at the nodes of the cosmic web differs from the large- scale structure of the web where the filaments mark the intersections of slanted sheets. The stream plane is only weakly aligned with the angular momentum (AM) near Rv, consistent with the fact that typically 80% of the AM is carried by one dominant stream. The galactic disc plane shows a weak tendency to be perpendicular to the large-scale SP, consistent with tidal-torque theory. Most interesting, the direction of the disc AM is only weakly correlated with the AM direction at Rv. This indicates a significant AM exchange at the interphase between streams and disc in the greater environment of the disc inside an "AM sphere" of radius ~0.3Rv . The required large torques are expected based on the perturbed morphology and kinematics within this interaction sphere. This AM exchange may or may not require a major modification of the standard disc modeling based on AM conservation, depending on the extent to which the amplitude of the disc AM is affected, which is yet to be studied.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dependence of Galaxy Quenching on Halo Mass and Distance from its Centre

We study the dependence of star-formation quenching on galaxy mass and environment, in the SDSS (z~0.1) and the AEGIS (z~1). It is crucial that we define quenching by low star-formation rate rather than by red colour, given that one third of the red galaxies are star forming. We address stellar mass M*, halo mass Mh, density over the nearest N neighbours deltaN, and distance to the halo centre D. The fraction of quenched galaxies appears more strongly correlated with Mh at fixed M* than with M* at fixed Mh, while for satellites quenching also depends on D. We present the M*-Mh relation for centrals at z~1. At z~1, the dependence of quenching on M* at fixed Mh is somewhat more pronounced than at z~0, but the quenched fraction is low (10%) and the haloes are less massive. For satellites, M*-dependent quenching is noticeable at high D, suggesting a quenching dependence on sub-halo mass for recently captured satellites. At small D, where satellites likely fell in more than a few Gyr ago, quenching strongly depends on Mh, and not on M*. The Mh-dependence of quenching is consistent with theoretical wisdom where virial shock heating in massive haloes shuts down accretion and triggers ram-pressure stripping, causing quenching. The interpretation of deltaN is complicated by the fact that it depends on the number of observed group members compared to N, motivating the use of D as a better measure of local environment.

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
preprint2012arXiv

GOODS-Herschel and CANDELS: The Morphologies of Ultraluminous Infrared Galaxies at z~2

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Luminous and High Stellar Mass Candidate Galaxies at z ~ 8 Discovered in CANDELS

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Metallicity-Dependent quenching of Star Formation at High Redshift in Small Galaxies

[abridged] The star formation rates (SFR) of low-metallicity galaxies depend sensitively on the gas metallicity, because metals are crucial to mediating the transition from intermediate-temperature atomic gas to cold molecular gas, a necessary precursor to star formation. We study the impact of this effect on the star formation history of galaxies. We incorporate metallicity-dependent star formation and metal enrichment in a simple model that follows the evolution of a halo main progenitor. Our model shows that including the effect of metallicity leads to suppression of star formation at redshift z>2 in dark halos with masses <~ 10^11 Msun, with the suppression becoming near total for halos below ~10^9.5-10 Msun. We find that at high redshift the SFR cannot catch up with the gas inflow rate (IR), because the SFR is limited by the free-fall time, and because it is suppressed further by a lack of metals. As a result, in each galaxy the SFR is growing in time faster than the IR, and the integrated cosmic SFR density is rising with time. The suppressed in situ SFR at high z makes the growth of stellar mass dominated by ex situ SFR which implies that the specific SFR (sSFR) remains constant with time. The intensely accreted gas at high z is accumulating as an atomic gas reservoir. This provides additional fuel for star formation in 10^10 - 10^12 Msun halos at z ~ 1-3, which allows the SFR to exceed the instantaneous IR, and may enable an even higher outflow rate. At z<1, following the natural decline in IR with time due to the universal expansion, the SFR and sSFR are expected to drop. We specify the expected dependence of sSFR and metallicity on stellar mass and redshift. At a given z, and below a critical mass, these relations are predicted to be flat and rising respectively. Our model predictions qualitatively match some of the puzzling features in the observed star formation history.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Combining Galaxy Clustering and Weak Lensing to Unveil Galaxy Biasing via the Halo Model

We formulate the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics. Using two-point statistics in projection, we define the galaxy bias function, b_g(r_p), and the galaxy-dark matter cross-correlation function, R_{gm}(r_p), where r_p is the projected distance. We use the analytical halo model to predict how the scale dependence of b_g and R_{gm}, over the range 0.1 Mpc/h < r_p < 30 Mpc/h, depends on the non-linearity and stochasticity in halo occupation models. In particular we quantify the effect due to the presence of central galaxies, the assumption for the radial distribution of satellite galaxies, the richness of the halo, and the Poisson character of the probability to have a certain number of satellite galaxies in a halo of a certain mass. Overall, brighter galaxies reveal a stronger scale dependence, and out to a larger radius. In real-space, we find that galaxy bias becomes scale independent, with R_{gm} = 1, for radii r > 1 - 5 Mpc/h, depending on luminosity. However, galaxy bias is scale-dependent out to much larger radii when one uses the projected quantities defined in this paper. These projected bias functions have the advantage that they are more easily accessible observationally and that their scale dependence carries a wealth of information regarding the properties of galaxy biasing. To observationally constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing we propose the use of the bias function Gamma_{gm}(r_p)=b_g(r_p)/R_{gm}(r_p). This function is obtained via a combination of weak gravitational lensing and galaxy clustering, and it can be measured using existing and forthcoming imaging and spectroscopic galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Galaxies and Homology

The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to a set of three dimensional scaling constants which may differ from one galaxy to the next. Noting that we are free to choose the dimensional constants, we find that a set of hydrodynamic simulated galaxy merger remnants is significantly closer to homologous when the dimensional length constant is taken to be the radius containing equal amounts of dark and baryonic matter rather than the usual observationally motivated choice of the baryonic half-mass radius. Once the correct dimensional scaling constants are used, the stellar velocity dispersion anisotropy is essentially the sole source of the variation in the kinematic structure of these simulated merger remnants. In order to facilitate the use of these scaling constants to analyse observed galaxies, we calculated the relationship between our preferred dimensional scaling constants and the typical observationally accessible quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the effect of cosmological inflow on turbulence and instability in galactic discs

We analyse the evolution of turbulence and gravitational instability of a galactic disc in a quasi-steady state governed by cosmological inflow. We focus on the possibility that the coupling between the in-streaming gas and the disc is maximal, e.g., via dense clumps, and ask whether the streams could be the driver of turbulence in an unstable disc with a Toomre parameter Q~1. Our fiducial model assumes an efficiency of ~0.5 per dynamical time for the decay of turbulence energy, and ~0.02 for each of the processes that deplete the disc gas, i.e., star formation, outflow, and inflow within the disc into a central bulge. In this case, the in-streaming drives a ratio of turbulent to rotation velocity sigma/V~0.2-0.3, which at z~2 induces an instability with Q~1, both as observed. However, in conflict with observations, this model predicts that sigma/V remains constant with time, independent of the cosmological accretion rate, because mass and turbulence have the same external source. Such strongly coupled cosmological inflow tends to stabilize the disc at low z, with Q ~ a few, which may be consistent with observations. The instability could instead be maintained for longer, with a properly declining sigma/V, if it is self-regulated to oscillations about Q~1 by a duty cycle for disc depletion. However, the 'off' phases of this duty cycle become long at low z, which may be hard to reconcile with observations. Alternatively, the coupling between the in-streaming gas and the disc may weaken in time, reflecting an evolving nature of the accretion. If, instead, that coupling is weak at all times, the likely energy source for self-regulated stirring up of the turbulence is the inflow within the disc down the potential gradient (studied in a companion paper).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Smooth(er) Stellar Mass Maps in CANDELS: Constraints on the Longevity of Clumps in High-redshift Star-forming Galaxies

We perform a detailed analysis of the resolved colors and stellar populations of a complete sample of 323 star-forming galaxies at 0.5 < z < 1.5, and 326 star-forming galaxies at 1.5 < z < 2.5 in the ERS and CANDELS-Deep region of GOODS-South. Galaxies were selected to be more massive than 10^10 Msun and have specific star formation rates above 1/t_H. We model the 7-band optical ACS + near-IR WFC3 spectral energy distributions of individual bins of pixels, accounting simultaneously for the galaxy-integrated photometric constraints available over a longer wavelength range. We analyze variations in rest-frame color, stellar surface mass density, age, and extinction as a function of galactocentric radius and local surface brightness/density, and measure structural parameters on luminosity and stellar mass maps. We find evidence for redder colors, older stellar ages, and increased dust extinction in the nuclei of galaxies. Big star-forming clumps seen in star formation tracers are less prominent or even invisible on the inferred stellar mass distributions. Off-center clumps contribute up to ~20% to the integrated SFR, but only 7% or less to the integrated mass of all massive star-forming galaxies at z ~ 1 and z ~ 2, with the fractional contributions being a decreasing function of wavelength used to select the clumps. The stellar mass profiles tend to have smaller sizes and M20 coefficients, and higher concentration and Gini coefficients than the light distribution. Our results are consistent with an inside-out disk growth scenario with brief (100 - 200 Myr) episodic local enhancements in star formation superposed on the underlying disk. Alternatively, the young ages of off-center clumps may signal inward clump migration, provided this happens efficiently on the order of an orbital timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Dependence of Quenching upon the Inner Structure of Galaxies at 0.5<z< 0.8 in the DEEP2/AEGIS Survey

The shutdown of star formation in galaxies is generally termed `quenching'. Although quenching may occur through a variety of processes, the exact mechanism(s) that is in fact responsible for quenching is still in question. This paper addresses quenching by searching for traces of possible quenching processes through their effects on galaxy structural parameters such as surface stellar mass density and Sersic index (n). We analyze the rest-frame U-B color correlations versus these structural parameters using a sample of galaxies in the redshift range 0.5< z<0.8 from the DEEP2/AEGIS survey. We find that Sersic index (n) has the smallest overlap region among all tested parameters and resembles a step-function with a threshold value of n=2.3. There exists, however, a significant population of outliers with blue colors yet high n values that seem to contradict this behavior. We hypothesize that their Sersic values may be distorted by bursts of star formation, AGNs, and/or poor fits, leading us to consider central surface stellar mass density as an alternative to Sersic index. Not only does it correct the outliers, it also forms a tight relationship with color, suggesting that the innermost structure of galaxies is most physically linked with quenching. Furthermore, at z~0.65, the majority of the blue cloud galaxies cannot simply fade onto the red sequence since their GIM2D bulge masses are only half as large on average as the bulge masses of similar red sequence galaxies, thus demonstrating that stellar mass must absolutely increase at the centers of galaxies as they quench. We discuss a two-stage model for quenching in which galaxy star formation rates are controlled by their dark halos while they are still in the blue cloud and a second quenching process sets in later, associated with the central stellar mass build-up.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

What turns galaxies off? The different morphologies of star-forming and quiescent galaxies since z~2 from CANDELS

We use HST/WFC3 imaging from the CANDELS Multicycle Treasury Survey, in conjunction with the Sloan Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses >3e10M_sun from z=2.2 to the present epoch, a time span of 10Gyr. We explore the relationship between rest-frame optical color, stellar mass, star formation activity and galaxy structure. We confirm the dramatic increase from z=2.2 to the present day in the number density of non-star-forming galaxies above 3e10M_sun reported by others. We further find that the vast majority of these quiescent systems have concentrated light profiles, as parametrized by the Sersic index, and the population of concentrated galaxies grows similarly rapidly. We examine the joint distribution of star formation activity, Sersic index, stellar mass, inferred velocity dispersion, and stellar surface density. Quiescence correlates poorly with stellar mass at all z<2.2. Quiescence correlates well with Sersic index at all redshifts. Quiescence correlates well with `velocity dispersion' and stellar surface density at z>1.3, and somewhat less well at lower redshifts. Yet, there is significant scatter between quiescence and galaxy structure: while the vast majority of quiescent galaxies have prominent bulges, many of them have significant disks, and a number of bulge-dominated galaxies have significant star formation. Noting the rarity of quiescent galaxies without prominent bulges, we argue that a prominent bulge (and perhaps, by association, a supermassive black hole) is an important condition for quenching star formation on galactic scales over the last 10Gyr, in qualitative agreement with the AGN feedback paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Universal, Local Star Formation Law in Galactic Clouds, Nearby Galaxies, High-Redshift Disks, and Starbursts

[abridged] While observations of Local Group galaxies show a very simple, local star formation law in which the star formation rate per unit area in each patch of a galaxy scales linearly with the molecular gas surface density, recent observations of both Milky Way molecular clouds and high redshift galaxies apparently show a more complicated relationship, in which regions of equal surface density can form stars at quite different rates. These data have been interpreted as implying either that different star formation laws apply in different circumstances, that the star formation law is sensitive to large-scale galaxy properties rather than local properties, or that there are high density thresholds for star formation. Here we collate resolved observations of Milky Way molecular clouds, kpc-scale observations of Local Group galaxies, and unresolved observations of both disk and starburst galaxies in the local universe and at high redshift. We show that all of these data are in fact consistent with a simple, local, volumetric star formation law. The apparent variations stem from the fact that the observed objects have a wide variety of 3D size scales and degrees of internal clumping, so even at fixed gas column density the regions being observed can have wildly varying volume densities. We provide a simple theoretical framework to remove this projection effect, and we show that all the data, from small Solar neighborhood clouds with masses ~10^3 Msun to sub-mm galaxies with masses ~10^11 Msun, fall on a single star formation law in which the SFR is simply ~1% of the molecular gas mass per local free-fall time. In contrast, proposed star formation laws in which the star formation timescale is set by the galactic rotation period or the SFR is linearly proportional to the gas mass above some density threshold fail to match at least some of the data.

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

Black Hole growth and AGN obscuration by instability-driven inflows in high-redshift disk galaxies fed by cold streams

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dark halo response and the stellar initial mass function in early-type and late-type galaxies

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Quenching by Clumpy Accretion in Cool Core Clusters: Convective Dynamical Response to Overheating

Many galaxy clusters pose a "cooling-flow problem", where the observed X-ray emission from their cores is not accompanied by enough cold gas or star formation. A continuous energy source is required to balance the cooling rate over the whole core volume. We address the feasibility of a gravitational heating mechanism, utilizing the gravitational energy released by the gas that streams into the potential well of the cluster dark-matter halo. We focus here on a specific form of gravitational heating in which the energy is transferred to the medium thorough the drag exerted on inflowing gas clumps. Using spheri-symmetric hydro simulations with a subgrid representation of these clumps, we confirm our earlier estimates that in haloes >=10^13 solar masses the gravitational heating is more efficient than the cooling everywhere. The worry was that this could overheat the core and generate an instability that might push it away from equilibrium. However, we find that the overheating does not change the global halo properties, and that convection can stabilize the cluster by carrying energy away from the overheated core. In a typical rich cluster of 10^{14-15}solar masses, with ~5% of the accreted baryons in gas clumps of ~10^8 solar masses, we derive upper and lower limits for the temperature and entropy profiles and show that they are consistent with those observed in cool-core clusters. We predict the density and mass of cold gas and the level of turbulence driven by the clump accretion. We conclude that gravitational heating is a feasible mechanism for preventing cooling flows in clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the puzzling plateau in the specific star formation rate at z=2-7

The observational indications for a constant specific star-formation rate (sSFR) in the redshift range z=2-7 are puzzling in the context of current galaxy-formation models. Despite the tentative nature of the data, their marked conflict with theory motivates a study of the possible implications. The plateau at sSFR ~ 2 Gyr^-1 is hard to reproduce because (a) its level is low compared to the cosmological specific accretion rate at z > 6, (b) it is higher than the latter at z ~ 2, (c) the natural correlation between SFR and stellar mass makes it difficult to manipulate their ratio, and (d) a low SFR at high z makes it hard to produce enough massive galaxies by z ~ 2. Using a flexible semi-analytic model, we explore ad-hoc modifications to the standard physical recipes trying to obey the puzzling observational constraints. Successful models involve non-trivial modifications, such as (a) a suppressed SFR at z > 4 in galaxies of all masses, by enhanced feedback or reduced SFR efficiency, following an initial active phase at z > 7, (b) a delayed gas consumption into stars, allowing the gas that was prohibited from forming stars or ejected at high z to form stars later in more massive galaxies, and (c) enhanced growth of massive galaxies, in terms of either faster assembly or more efficient starbursts in mergers, or by efficient star formation in massive haloes.

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

On the Origin of the Galaxy Star-Formation-Rate Sequence: Evolution and Scatter

We use a semi-analytic model for disk galaxies to explore the origin of the time evolution and small scatter of the galaxy SFR sequence -- the tight correlation between star-formation rate (SFR) and stellar mass (M_star). The steep decline of SFR from z~2 to the present, at fixed M_star, is a consequence of the following: First, disk galaxies are in a steady state with the SFR following the net (i.e., inflow minus outflow) gas accretion rate. The evolution of the SFR sequence is determined by evolution in the cosmological specific accretion rates, \propto (1+z)^{2.25}, but is found to be independent of feedback. Although feedback determines the outflow rates, it shifts galaxies along the SFR sequence, leaving its zero point invariant. Second, the conversion of accretion rate to SFR is materialized through gas density, not gas mass. Although the model SFR is an increasing function of both gas mass fraction and gas density, only the gas densities are predicted to evolve significantly with redshift. Third, star formation is fueled by molecular gas. Since the molecular gas fraction increases monotonically with increasing gas density, the model predicts strong evolution in the molecular gas fractions, increasing by an order of magnitude from z=0 to z~2. On the other hand, the model predicts that the effective surface density of atomic gas is ~10 M_sun pc^{-2}, independent of redshift, stellar mass or feedback. Our model suggests that the scatter in the SFR sequence reflects variations in the gas accretion history, and thus is insensitive to stellar mass, redshift or feedback. The large scatter in halo spin contributes negligibly, because it scatters galaxies along the SFR sequence. An observational consequence of this is that the scatter in the SFR sequence is independent of the size (both stellar and gaseous) of galaxy disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Survival of Star-Forming Giant Clumps in High-Redshift Galaxies

We investigate the effects of radiation pressure from stars on the survival of the star-forming giant clumps in high-redshift massive disc galaxies, during the most active phase of galaxy formation. The clumps, typically of mass ~10^8-10^9 Msun and radius ~0.5-1, are formed in the turbulent gas-rich discs by violent gravitational instability and then migrate into a central bulge in ~10 dynamical times. We show that the survival or disruption of these clumps under the influence of stellar feedback depends critically on the rate at which they form stars. If they convert a few percent of their gas mass to stars per free-fall time, as observed for all local star-forming systems and implied by the Kennicutt-Schmidt law, they cannot be disrupted. Only if clumps convert most of their mass to stars in a few free-fall times can feedback produce significant gas expulsion. We consider whether such rapid star formation is likely in high-redshift giant clumps.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Morphological quenching of star formation: making early-type galaxies red

We point out a natural mechanism for quenching of star formation in early-type galaxies. It automatically links the color of a galaxy with its morphology and does not require gas consumption, removal or termination of gas supply. Given that star formation takes place in gravitationally unstable gas disks, it can be quenched when a disk becomes stable against fragmentation to bound clumps. This can result from the growth of a stellar spheroid, for instance by mergers. We present the concept of morphological quenching (MQ) using standard disk instability analysis, and demonstrate its natural occurrence in a cosmological simulation using an efficient zoom-in technique. We show that the transition from a stellar disk to a spheroid can be sufficient to stabilize the gas disk, quench star formation, and turn an early-type galaxy red and dead while gas accretion continues. The turbulence necessary for disk stability can be stirred up by sheared perturbations within the disk in the absence of bound star-forming clumps. While gas stripping processes are limited to dense groups and clusters, and other quenching mechanisms like AGN feedback, virial shock heating and gravitational heating, are limited to halos more massive than 10^12 Mo, the MQ can explain the appearance of red ellipticals even in less massive halos and in the field. The dense gas disks observed in some of today's red ellipticals may be the relics of this mechanism, whereas red galaxies with quenched gas disks are expected to be more frequent at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tidal Effects and the Environment Dependence of Halo Assembly

We explore a possible origin for the puzzling anti-correlation between the formation epoch of galactic dark-matter haloes and their environment density. This correlation has been revealed from cosmological N-body simulations and is in conflict with the Extended Press-Schechter model of halo clustering. Using similar simulations, we first quantify the straightforward association of an early formation epoch with a reduced mass growth rate at late times. We then find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the assembly of haloes due to the shear along the large-scale filaments that feed the massive halo to actual mass loss in haloes that pass through the massive halo. Using the restricted three-body problem, we show that haloes are prone to tidal mass loss within 1.5 virial radii of a larger halo. Our results suggest that the dependence of formation epoch on environment density is a secondary effect induced by the enhanced density of haloes in filaments near massive haloes where the tides are strong. Our measures of assembly rate are particularly correlated with the tidal field at high redshifts z~1.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Universal Merger Histories of Dark-Matter Haloes

We study merger histories of dark-matter haloes in a suite of N-body simulations that span different cosmological models. The simulated cases include the up-to-date WMAP5 cosmology and other test cases based on the Einstein-deSitter cosmology with different power spectra. We provide a robust fitting function for the conditional mass function (CMF) of progenitor haloes of a given halo. This fit is valid for the different cosmological models and for different halo masses and redshifts, and it is a significant improvement over earlier estimates. Based on this fit, we develop a simple and accurate technique for transforming the merger history of a given simulated halo into haloes of different mass, redshift and cosmology. Other statistics such as main-progenitor history and merger rates are accurately transformed as well. This method can serve as a useful tool for studying galaxy formation. It is less sensitive to the low accuracy of the fit at small time-steps, and it can thus replace the more elaborate task of construction Monte-Carlo realizations. As an alternative approach, we confirm the earlier finding by Neistein & Dekel that the main-progenitor follows a log-normal distribution. This property of merger trees allows us to better capture their behaviour as a function of time and descendant mass, but a broader suite of simulations is required for evaluating the dependence of the log-normal parameters on the cosmological model.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Origin of the Joint Distribution of Structural Parameters in Disc Galaxies

We use a simple model of disc-galaxy formation in the LCDM cosmology to simultaneously reproduce the slopes, zero-points, scatter and uncorrelated residuals in the observed velocity-luminosity (VL) and radius-luminosity (RL) relations. Observed I-band luminosities are converted to stellar masses using the IMF-dependent relation between stellar mass-to-light ratio, M/L, and color. The model treats halo concentration, spin parameter, and disc mass fraction as independent log-normal random variables. Our main conclusion is that the VL and RL zero-points and the uncorrelated residuals can only be reproduced simultaneously if adiabatic contraction is avoided. One or the other could be fixed by appealing to unrealistic values for M/L, halo concentration c, and spin parameter lambda, but not all together. The small VL scatter is naturally determined by the predicted scatter in c and in M/L, quite independent of the large scatter in lambda. However, the RL scatter, driven by the scatter in lambda, can be as low as observed only if sigma_lambda simeq 0.25, about half the value predicted for CDM haloes. This may indicate that discs form in a special subset of haloes. The VL slope is reproduced once star formation occurs only above a threshold surface density, and the RL slope implies that the disc mass fraction is increasing with halo mass, consistent with feedback effects. A model that incorporates the above ingredients provides a simultaneous fit to all the observed features. In particular, the elimination of adiabatic contraction allows 80% disc contribution to the observed rotation velocity of bright discs at 2.2 disc scale lengths. The lack of halo contraction may indicate that disc formation is not as smooth as typically envisioned, but instead may involve clumpy, cold streams.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Virial shocks in galactic haloes?

We investigate the conditions for the existence of an expanding virial shock in the gas falling within a spherical dark-matter halo. The shock relies on pressure support by the shock-heated gas behind it. When the radiative cooling is efficient compared to the infall rate the post-shock gas becomes unstable; it collapses inwards and cannot support the shock. We find for a monoatomic gas that the shock is stable when the post-shock pressure and density obey gamma effective>10/7, with gamma effective begin the time depended equivalent to the adiabatic index. We express the effective gamma in terms of r, u and rho at the shock to obtain a simple condition for shock stability. This result is confirmed by hydrodynamical simulations, using an accurate spheri-symmetric Lagrangian code. When the stability analysis is applied in cosmology, we find that a virial shock does not develop in most haloes that form before z ~ 2, and it never forms in haloes less massive than a few 10^11 solar masses. In such haloes the infalling gas is never heated to the virial temperature, and it does not need to cool radiatively before it drops into a disc. Instead, the gas collapses at T ~ 10^4K directly into the disc. This should have nontrivial effects on the star-formation rate and on the gas removal by supernova-driven winds. Instead of radiating soft x rays, this gas would emit lyman alpha thus helping explain the low flux of soft x-ray background and the lyman alpha emitters observed at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

A universal angular momentum profile for galactic halos

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Profiles of dark haloes: evolution, scatter, and environment

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Large-Scale Power-Spectrum from Peculiar Velocities

The power spectrum (PS) of mass density fluctuations, in the range $0.05 \leq k \leq 0.2 \ihmpc$, is derived from the Mark III catalog of peculiar velocities of galaxies, independent of ``biasing". It is computed from the density field as recovered by POTENT with Gaussian smoothing of $12\hmpc$, within a sphere of radius $\sim 60\hmpc$ about the Local Group. The density is weighted inversely by the square of the errors. The PS is corrected for the effects of smoothing, random errors, sparse sampling and finite volume using mock catalogs that mimic in detail the Mark III catalog and the dynamics of our cosmological neighborhood. The mock catalogs are also used for error analysis. The PS at $k = 0.1 \ihmpc$ (for example) is $(4.6 \pm 1.4) \times 10^3 Ω^{-1.2} \3hmpc$, with a local logarithmic slope of $-1.45 \pm 0.5$. An itegration yields $σ_8 Ω^{0.6} \simeq 0.7-0.8$, depending on where the PS peak is. Direct comparisons of the mass PS with the galaxy PS derived from sky and redshift surveys show a similarity in shape and yield for $β\equiv Ω^{0.6}/b$ values in the range $0.77-1.21$, with a typical error of $\pm 0.1$ per galaxy sample. A comparison of the mass PS at $\sim 100\hmpc$ with the large-angle CMB fluctuations by COBE provides constraints on cosmological parameters and on the slope of the initial PS ($n$). The ``standard" CDM model is marginally rejected at the $\sim 2σ$ level, while each of the following modifications leads to a good fit: $n \lsim 1$, $Ω_ν\sim 0.3$, or $Ω\lsim 1$. Values of $Ω$ as low as $\sim 0.2$ are ruled out with high confidence (independent of $Λ$).

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Large-Scale Mass Distribution Behind the Galactic Plane

We map the smoothed mass-density distribution in the Galactic zone of avoidance (zoa), within $6000\kms$ of the Local Group, using potent reconstruction from peculiar velocities of galaxies. The interpolation into the zoa is based on the assumed gravitational nature of the velocity field implying a potential flow. The main dynamical features at a distance r=4000 km/s are (a) the peak of the Great Attractor connecting Centaurus and Pavo at l=330, (b) a moderate bridge connecting Perseus-Pisces and Cepheus at l= 140, and (c) an extension of a large void from the southern Galactic hemisphere into the zoa near the direction of Puppis, l=220-270. We find a strong correlation between the mass density and the IRAS and optical galaxy density at $b=+-20, which indicates that the main dynamical features in the zoa should also be seen in galaxy surveys through the Galactic plane. The gravitational acceleration at the Local Group, based on the mass distribution out to 6000 km/s, is strongly affected by the mass distribution in the zoa: its direction changes by 31 when the |b| < 20 zoa is included, bringing it to within 4 +- 19 of the CMB dipole.

0 citations0 reviewsNo code linkedOpen access