Showing papers by "Joel R. Primack published in 2017"

Structural and Star-forming Relations since z similar to 3: Connecting Compact Star-forming and Quiescent Galaxies

Abstract: We study the evolution of the scaling relations that compare the effective density (Sigma(e), r 9.6 -9.3M(circle dot) kpc(-2), allowing the most efficient identification of compact SFGs and quiescent galaxies at every redshift.

Constraining the galaxy–halo connection over the last 13.3 Gyr: star formation histories, galaxy mergers and structural properties

Abstract: We present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star formation histories, merger rates, and structural properties for average galaxies, combining star-forming and quenched galaxies. Our main findings: (1) The halo mass $M_{50}$ above which 50\% of galaxies are quenched coincides with sSFR/sMAR$\sim1$, where sMAR is the specific halo mass accretion rate. (2) $M_{50}$ increases with redshift, presumably due to cold streams being more efficient at high redshift while virial shocks and AGN feedback become more relevant at lower redshifts. (3) The ratio sSFR/sMAR has a peak value, which occurs around $M_{\rm vir}\sim2\times10^{11}M_{\odot}$. (4) The stellar mass density within 1 kpc, $\Sigma_1$, is a good indicator of the galactic global sSFR. (5) Galaxies are statistically quenched after they reach a maximum in $\Sigma_1$, consistent with theoretical expectations of the gas compaction model; this maximum depends on redshift. (6) In-situ star formation is responsible for most galactic stellar mass growth, especially for lower-mass galaxies. (7) Galaxies grow inside out. The marked change in the slope of the size--mass relation when galaxies became quenched, from $d\log R_{\rm eff}/d\log M_*\sim0.35$ to $\sim2.5$, could be the result of dry minor mergers.

Giant clumps in simulated high-z Galaxies: properties, evolution and dependence on feedback

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

The nature of massive transition galaxies in CANDELS, GAMA and cosmological simulations

Abstract: We explore observational and theoretical constraints on how galaxies might transition between the 'star-forming main sequence' (SFMS) and varying 'degrees of quiescence' out to z = 3. Our analysis is focused on galaxies with stellar mass M* > 1010 M⊙, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical 'zoom in' simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate Sersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to z = 3. We place an observational upper limit on the average population transition time-scale as a function of redshift, finding that the average high-redshift galaxy is on a 'fast track' for quenching whereas the average low-redshift galaxy is on a 'slow track' for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching time-scales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at z ∼ 3. In the SAM, we do not find a clear-cut morphological dependence of quenching time-scales, but we do predict that the mean stellar ages, cold gas fractions, SMBH (supermassive black hole) masses and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at z < 3.

CANDELS: elevated black hole growth in the progenitors of compact quiescent galaxies at z ~ 2

Abstract: We examine the fraction of massive (M* > 1010 M) compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at z ~ 2. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that - 39.2+ % 3.6 3.9 (65/166) of cSFGs at 1.4 < z < 3.0 host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the 6.2s level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

z ∼ 2: An Epoch of Disk Assembly

Abstract: We explore the evolution of the internal gas kinematics of star-forming galaxies from the peak of cosmic star-formation at $z\sim2$ to today. Measurements of galaxy rotation velocity $V_{rot}$, which quantify ordered motions, and gas velocity dispersion $\sigma_g$, which quantify disordered motions, are adopted from the DEEP2 and SIGMA surveys. This sample covers a continuous baseline in redshift from $z=2.5$ to $z=0.1$, spanning 10 Gyrs. At low redshift, nearly all sufficiently massive star-forming galaxies are rotationally supported ($V_{rot}>\sigma_g$). By $z=2$, the percentage of galaxies with rotational support has declined to 50$\%$ at low stellar mass ($10^{9}-10^{10}\,M_{\odot}$) and 70$\%$ at high stellar mass ($10^{10}-10^{11}M_{\odot}$). For $V_{rot}\,>\,3\,\sigma_g$, the percentage drops below 35$\%$ for all masses. From $z\,=\,2$ to now, galaxies exhibit remarkably smooth kinematic evolution on average. All galaxies tend towards rotational support with time, and it is reached earlier in higher mass systems. This is mostly due to an average decline in $\sigma_g$ by a factor of 3 since a redshift of 2, which is independent of mass. Over the same time period, $V_{rot}$ increases by a factor of 1.5 for low mass systems, but does not evolve for high mass systems. These trends in $V_{rot}$ and $\sigma_g$ with time are at a fixed stellar mass and should not be interpreted as evolutionary tracks for galaxy populations. When galaxy populations are linked in time with abundance matching, not only does $\sigma_g$ decline with time as before, but $V_{rot}$ strongly increases with time for all galaxy masses. This enhances the evolution in $V_{rot}/\sigma_g$. These results indicate that $z\,=\,2$ is a period of disk assembly, during which the strong rotational support present in today's massive disk galaxies is only just beginning to emerge.

Clumpy Galaxies in CANDELS. II. Physical Properties of UV-bright Clumps at $0.5\leq z<3$

Abstract: Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged]

Properties of dark matter haloes as a function of local environment density

Abstract: We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and rounder shapes than halos in median density environments. Halos in median and low-density environments have similar accretion rates and concentrations, but halos in low density environments have lower spins and are more elongated. Halos of a given mass in high-density regions accrete material earlier than halos of the same mass in lower-density regions. All but the most massive halos in high-density regions are losing mass (i.e., being stripped) at low redshifts, which causes artificially lowered NFW scale radii and increased concentrations. Tidal effects are also responsible for the decreasing spins of low mass halos in high density regions at low redshifts $z < 1$, by preferentially removing higher angular momentum material from halos. Halos in low-density regions have lower than average spins because they lack nearby halos whose tidal fields can spin them up. We also show that the simulation density distribution is well fit by an Extreme Value Distribution, and that the density distribution becomes broader with cosmic time.

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

Abstract: We study the correlation of galaxy structural properties with their location relative to the SFR–M* correlation, also known as the star formation ‘star-forming main sequence’ (SFMS), in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey and Galaxy and Mass Assembly Survey 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 SAM, bulge growth is driven by mergers and disc instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via active galactic nucleus feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and active galactic nucleus feedback plays a critical role in moving galaxies off of the SFMS.

Spatially Resolved Kinematics in the Central 1 kpc of a Compact Star-forming Galaxy at z ∼ 2.3 from ALMA CO Observations

Abstract: We present high spatial resolution (FWHM similar to 0."14) observations of the CO(8-7) line in GDS-14876, a compact star-forming galaxy at z = 2.3 with a total stellar mass of log(M-*/M-circle dot) = 10.9. The spatially resolved velocity map of the inner r less than or similar to 1 kpc reveals a continuous velocity gradient consistent with the kinematics of a rotating disk with v(rot)(r = 1 kpc) = 163. +/-. 5 km.s(-1) and v(rot)/s similar to 2.5. The gas-to-stellar ratios estimated from CO(8-7) and the dust continuum emission span a broad range, f(gas)(CO) = M-gas/M-* = 13%-45% and f(gas)(CO) = 50%-67% but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of log(M-dyn/M-circle dot.) = 10.5(-0.2)(+0.5), which is lower, but still consistent with the baryonic mass, log(M-bar = M-circle dot + M-gas(CO)/M-circle dot)= 11.0 if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO(8-7), similar to 3x smaller than the stellar size, implies a strong relative concentration that increases the gas fraction up to f(gas)(CO) similar to 85% , 1 kpc in the central 1 kpc. Such a gas-rich center, coupled with a high star formation rate (SFR) similar to 500. M-circle dot yr(-1), suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly (t(depl) similar to 27 Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.

z~2: An Epoch of Disk Assembly

Abstract: We explore the evolution of the internal gas kinematics of star-forming galaxies from the peak of cosmic star-formation at $z\sim2$ to today. Measurements of galaxy rotation velocity $V_{rot}$, which quantify ordered motions, and gas velocity dispersion $\sigma_g$, which quantify disordered motions, are adopted from the DEEP2 and SIGMA surveys. This sample covers a continuous baseline in redshift from $z=2.5$ to $z=0.1$, spanning 10 Gyrs. At low redshift, nearly all sufficiently massive star-forming galaxies are rotationally supported ($V_{rot}>\sigma_g$). By $z=2$, the percentage of galaxies with rotational support has declined to 50$\%$ at low stellar mass ($10^{9}-10^{10}\,M_{\odot}$) and 70$\%$ at high stellar mass ($10^{10}-10^{11}M_{\odot}$). For $V_{rot}\,>\,3\,\sigma_g$, the percentage drops below 35$\%$ for all masses. From $z\,=\,2$ to now, galaxies exhibit remarkably smooth kinematic evolution on average. All galaxies tend towards rotational support with time, and it is reached earlier in higher mass systems. This is mostly due to an average decline in $\sigma_g$ by a factor of 3 since a redshift of 2, which is independent of mass. Over the same time period, $V_{rot}$ increases by a factor of 1.5 for low mass systems, but does not evolve for high mass systems. These trends in $V_{rot}$ and $\sigma_g$ with time are at a fixed stellar mass and should not be interpreted as evolutionary tracks for galaxy populations. When galaxy populations are linked in time with abundance matching, not only does $\sigma_g$ decline with time as before, but $V_{rot}$ strongly increases with time for all galaxy masses. This enhances the evolution in $V_{rot}/\sigma_g$. These results indicate that $z\,=\,2$ is a period of disk assembly, during which the strong rotational support present in today's massive disk galaxies is only just beginning to emerge.

CANDELS: Elevated Black Hole Growth in the Progenitors of Compact Quiescent Galaxies at z~2

Abstract: We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4

Formation and settling of a disc galaxy during the last 8 billion years in a cosmological simulation

Abstract: We present results of a high-resolution zoom cosmological simulation of the evolution of a low-mass galaxy with a maximum velocity of V=100 km/s at z=0, using the initial conditions from the AGORA project (Kim et al. 2014). The final disc-dominated galaxy is consistent with local disc scaling relations, such as the stellar-to-halo mass relation and the baryonic Tully-Fisher. The galaxy evolves from a compact, dispersion-dominated galaxy into a rotation-dominated but dynamically hot disc in about 0.5 Gyr (from z=1.4 to z=1.2). The disc dynamically cools down for the following 7 Gyr, as the gas velocity dispersion decreases over time, in agreement with observations. The primary cause of this slow evolution of velocity dispersion in this low-mass galaxy is stellar feedback. It is related to the decline in gas fraction, and to the associated gravitational disk instability, as the disc slowly settles from a global Toomre Q>1 turbulent disc to a marginally unstable disc (Q=1).

Spatially resolved kinematics in the central 1 kpc of a compact star-forming galaxy at z=2.3 from ALMA CO observations

Abstract: We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/\sigma\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.

Understanding large-scale structure in the SSA22 protocluster region using cosmological simulations

Abstract: We investigate the nature and evolution of large-scale structure within the SSA22 protocluster region at $z=3.09$ using cosmological simulations. A redshift histogram constructed from current spectroscopic observations of the SSA22 protocluster reveals two separate peaks at $z = 3.065$ (blue) and $z = 3.095$ (red). Based on these data, we report updated overdensity and mass calculations for the SSA22 protocluster. We find $\delta_{b,gal}=4.8 \pm 1.8$, $\delta_{r,gal}=9.5 \pm 2.0$ for the blue and red peaks, respectively, and $\delta_{t,gal}=7.6\pm 1.4$ for the entire region. These overdensities correspond to masses of $M_b = (0.76 \pm 0.17) \times 10^{15} h^{-1} M_{\odot}$, $M_r = (2.15 \pm 0.32) \times 10^{15} h^{-1} M_{\odot}$, and $M_t=(3.19 \pm 0.40) \times 10^{15} h^{-1} M_{\odot}$ for the red, blue, and total peaks, respectively. We use the Small MultiDark Planck (SMDPL) simulation to identify comparably massive $z\sim 3$ protoclusters, and uncover the underlying structure and ultimate fate of the SSA22 protocluster. For this analysis, we construct mock redshift histograms for each simulated $z\sim 3$ protocluster, quantitatively comparing them with the observed SSA22 data. We find that the observed double-peaked structure in the SSA22 redshift histogram corresponds not to a single coalescing cluster, but rather the proximity of a $\sim 10^{15}h^{-1} M_{\odot}$ protocluster and at least one $>10^{14} h^{-1} M_{\odot}$ cluster progenitor. Such associations in the SMDPL simulation are easily understood within the framework of hierarchical clustering of dark matter halos. We finally find that the opportunity to observe such a phenomenon is incredibly rare, with an occurrence rate of $7.4h^3 \mbox{ Gpc}^{-3}$.

Understanding substructure in the SSA22 protocluster region using cosmological simulations

Abstract: We investigate the nature and evolution of substructure within the SSA22 protocluster region at $z=3.09$ using cosmological simulations. A redshift histogram constructed from current spectroscopic observations of the SSA22 protocluster reveals two separate peaks at $z = 3.065$ (blue) and $z = 3.095$ (red). Based on these data, we report updated overdensity and mass calculations for the SSA22 protocluster. We find $\delta_{b,gal}=4.8 \pm 1.8$, $\delta_{r,gal}=9.5 \pm 2.0$ for the blue and red peaks, respectively, and $\delta_{t,gal}=7.6\pm 1.4$ for the entire region. These overdensities correspond to masses of $M_b = (0.76 \pm 0.17) \times 10^{15} h^{-1} M_{\odot}$, $M_r = (2.15 \pm 0.32) \times 10^{15} h^{-1} M_{\odot}$, and $M_t=(3.19 \pm 0.40) \times 10^{15} h^{-1} M_{\odot}$ for the red, blue, and total peaks, respectively. We use the Small MultiDark Planck (SMDPL) simulation to identify comparably massive $z\sim 3$ protoclusters, and uncover the underlying structure and ultimate fate of the SSA22 protocluster. For this analysis, we construct mock redshift histograms for each simulated $z\sim 3$ protocluster, quantitatively comparing them with the observed SSA22 data. We find that the observed double-peaked substructure in the SSA22 redshift histogram corresponds not to a single coalescing cluster, but rather the proximity of a $\sim 10^{15}h^{-1} M_{\odot}$ protocluster and at least one $>10^{14} h^{-1} M_{\odot}$ cluster progenitor. Such associations in the SMDPL simulation are easily understood within the framework of hierarchical clustering of dark matter halos. We finally find that the opportunity to observe such a phenomenon is incredibly rare, with an occurrence rate of $8h^3 \mbox{ Gpc}^{-3}$.