Source author record

Aaron A. Dutton

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

49 published item(s)

preprint2020arXiv

Abundance matching tested on small scales with galaxy dynamics

We present a comprehensive test of the relation between stellar and total mass in galaxies as predicted by popular models based on abundance matching (AM) techniques. We use the ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey with photometric and dynamical profiles for 190 Virgo cluster galaxies to establish a relation between the stellar and dynamical masses measured within the isophotal radius $r_{23.5}$. Various dark matter and galaxy scaling relations are combined with results from the NIHAO suite of hydrodynamical simulations to recast AM predictions in terms of these observed quantities. Our results are quite insensitive to the exact choice of dark matter profile and halo response to baryon collapse. We find that theoretical models reproduce the slope and normalization of the observed stellar-to-halo mass relation (SHMR) over more than three orders of magnitude in stellar mass $(10^8 < M_*/M_{\odot} < 2 \times 10^{11})$. However, the scatter of the observed SHMR exceeds that of AM predictions by a factor $\sim$5. For systems with stellar masses exceeding $5 \times 10^{10}~M_{\odot}$, AM overpredicts the observed stellar masses for a given dynamical mass. The latter offset may support previous indications of a different stellar initial mass function in these massive galaxies. Overall our results support the validity of AM predictions on a wide dynamical range.

preprint2020arXiv

NIHAO-UHD: The properties of MW-like stellar disks in high resolution cosmological simulations

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

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.

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.

preprint2016arXiv

A spiral galaxy's mass distribution uncovered through lensing and dynamics

We investigate the matter distribution of a spiral galaxy with a counter-rotating stellar core, SDSS J1331+3628 (J1331), independently with gravitational lensing and stellar dynamical modelling. By fitting a gravitational potential model to a quadruplet of lensing images around J1331's bulge, we tightly constrain the mass inside the Einstein radius R_ein = (0.91 +/- 0.02)'' (~= 1.83 +/- 0.04 kpc) to within 4%: M_ein = (7.8 +/- 0.3) x 10^10 M_Sun. We model observed long-slit major axis stellar kinematics in J1331's central regions by finding Multi-Gaussian Expansion (MGE) models for the stellar and dark matter distribution that solve the axisymmetric Jeans equations. The lens and dynamical model are independently derived, but in very good agreement with each other around ~R_ein. We find that J1331's center requires a steep total mass-to-light ratio gradient. A dynamical model including an NFW halo (with virial velocity v_200 ~= 240 +/- 40 km/s and concentration c_200 ~= 8 +/- 2) and moderate tangential velocity anisotropy (beta_z ~= -0.4 +/- 0.1) can reproduce the signatures of J1331's counter-rotating core and predict the stellar and gas rotation curve at larger radii. However, our models do not agree with the observed velocity dispersion at large radii. We speculate that the reason could be a non-trivial change in structure and kinematics due to a possible merger event in J1331's recent past.

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.

preprint2016arXiv

NIHAO project II: Halo shape, phase-space density and velocity distribution of dark matter in galaxy formation simulations

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($ρ/σ^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

preprint2016arXiv

NIHAO VIII: Circum-galactic medium and outflows - The puzzles of HI and OVI gas distributions

We study the hot and cold circum-galactic medium (CGM) of 86 galaxies of the cosmological, hydrodynamical simulation suite NIHAO. NIHAO allows a study of how the $z=0$ CGM varies across 5 orders of magnitude of stellar mass using OVI and HI as proxies for hot and cold gas. The cool HI covering fraction and column density profiles match observations well, particularly in the inner CGM. OVI shows increasing column densities with mass, a trend seemingly echoed in the observations. As in multiple previous simulations, the OVI column densities in simulations are lower than observed and optically thick HI does not extend as far out as in observations. We take a look at the collisional ionisation fraction of OVI as a function of halo mass. We make observable predictions of the bipolarity of outflows and their effect on the general shape of the CGM. Bipolar outflows can be seen out to around 40 kpc in intermediate and low mass halos ($M_{\mathrm{halo}}<10^{11}M_{\mathrm{sun}}$), but outside that radius, the CGM is too well mixed to detect an elongated shape. Larger halos have extended gas discs beyond the stellar disc that dominate the shape of the inner CGM. The simulated CGM is remarkably spherical even in low mass simulations. The chemical enrichment of both halo and disc gas follow expected increasing trends as a function of halo mass that are well fit with power laws. These relations can be used in non-hydrodynamic models, such as semi-analytic models.

preprint2016arXiv

NIHAO X: Reconciling the local galaxy velocity function with Cold Dark Matter via mock HI observations

We used 87 high resolution hydrodynamical cosmological simulations from the NIHAO suite to investigate the relation between the maximum circular velocity (Vmax) of a dark matter halo in a collisionless simulation and the velocity width of the HI gas in the same halo in the hydrodynamical simulation. These two quantities are normally used to compare theoretical and observational velocity functions and have led to a possible discrepancy between observations and predictions based on the Cold Dark Matter (CDM) model. We show that below 100 km/s, there is clear bias between HI based velocities and Vmax, that leads to an underestimation of the actual circular velocity of the halo. When this bias is taken into account the CDM model has no trouble in reproducing the observed velocity function and no lack of low velocity galaxies is actually present. Our simulations also reproduce the linewidth - stellar mass (Tully-Fisher) relation and HI sizes, indicating that the HI gas in our simulations is as extended as observed. The physical reason for the lower than expected linewidths is that, in contrast to high mass galaxies, low mass galaxies no longer have extended thin HI rotating disks, as is commonly assumed.

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

preprint2016arXiv

Simulated ΛCDM analogues of the thin Plane of Satellites around the Andromeda galaxy are not kinematically coherent structures

A large fraction of the dwarf satellites orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and apparently kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model (ΛCDM). Using 21 high resolution cosmological simulations we analyse the kinematics of planes of satellites similar to the one around Andromeda. We find good agreement when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in agreement with the plane around our Galaxy. We find such planes to be common in our high concentration haloes. The number of co-rotating satellites obtained from the sign of the line-of-sight velocity shows large variations depending on the viewing angle and is consistent with that obtained from a sample with random velocities. We find that the clustering of angular momentum vectors of the satellites in the plane is a better measure of the kinematic coherence. Thus we conclude that the line-of- sight velocity is not well suited as a proxy for the kinematical coherence of the plane. Analysis of the kinematics of our planes shows a fraction of $\sim$30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that these planes are a transient feature and not kinematically coherent as would appear at first sight. Thus we expect some of the satellites in the plane around Andromeda to have high velocities perpendicular to the plane.

preprint2015arXiv

Evidence for Early Filamentary Accretion from the Andromeda Galaxy's Thin Plane of Satellites

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift.

preprint2015arXiv

NIHAO project I: Reproducing the inefficiency of galaxy formation across cosmic time with a large sample of cosmological hydrodynamical simulations

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

preprint2015arXiv

NIHAO V: Too big doesn't fail -- reconciling the conflict between LCDM predictions and the circular velocities of nearby field galaxies

We compare the half-light circular velocities, V_{1/2}, of dwarf galaxies in the Local Group to the predicted circular velocity curves of galaxies in the NIHAO suite of LCDM simulations. We use a subset of 34 simulations in which the central galaxy has a stellar luminosity in the range 0.5 x 10^5 < L_V < 2 x 10^8 L_{sun}. The NIHAO galaxy simulations reproduce the relation between stellar mass and halo mass from abundance matching, as well as the observed half-light size vs luminosity relation. The corresponding dissipationless simulations over-predict the V_{1/2}, recovering the problem known as too big to fail (TBTF). By contrast, the NIHAO simulations have expanded dark matter haloes, and provide an excellent match to the distribution of V_{1/2} for galaxies with L_V > 2 x 10^6 L_{sun}. For lower luminosities our simulations predict very little halo response, and tend to over predict the observed circular velocities. In the context of LCDM, this could signal the increased stochasticity of star formation in haloes below M_{halo} \sim 10^{10} M_{sun}, or the role of environmental effects. Thus, haloes that are "too big to fail", do not fail LCDM, but haloes that are "too small to pass" (the galaxy formation threshold) provide a future test of LCDM.

preprint2015arXiv

On the Global Mass Distribution in Disk Galaxies

We present a summary of recent assessments of the mass distribution in disk galaxies. Of issue in order to characterize galaxy formation models is to determine the relative fraction of baryons and dark matter at all radii in galaxies. For disk galaxies, various measurements of the mass distribution in galaxies based on vertical kinematics, strong lensing, residuals of scaling relations, fluid dynamical modeling, bar strength and pattern speed, warps, and others, have called for either a maximal or sub-maximal contribution of the baryons in the inner parts of the disk. We propose a global picture whereby all galaxies are typically baryon-dominated (maximal) at the center and dark-matter dominated (sub-maximal) in their outskirts. Using as a fiducial radius the peak of the rotation curve of a pure baryonic exponential disk at R_2.2=2.2 Rd, where Rd is the disk scale length, the transition from maximal to sub-maximal baryons occurs within or near R_2.2 for low-mass disk galaxies (with V_tot < 200 km/s) and beyond R_2.2 for more massive systems. The mass fraction's inverse dependence on circular velocity in disk galaxies is largely insensitive to the presence of a bar at R_2.2. The mean mass fractions of late- and early-type galaxies are shown to be different at the same fiducial radius and circular velocity, pointing to different galaxy formation mechanisms. Feedback, dynamical friction, size, and IMF variations are key ingredients for understanding these differences.

preprint2015arXiv

The response of dark matter haloes to elliptical galaxy formation: a new test for quenching scenarios

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\logρ/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms.

preprint2014arXiv

A mass-dependent density profile for dark matter haloes including the influence of galaxy formation

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($γ$) and outer ($β$) slopes of dark matter density, and the sharpness of transition between the slopes($α$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile.

preprint2014arXiv

Cold dark matter haloes in the Planck era: evolution of structural parameters for Einasto and NFW profiles

We present the evolution of the structure of relaxed cold dark matter haloes in the cosmology from the Planck satellite. Our simulations cover 5 decades in halo mass, from dwarf galaxies to galaxy clusters. Due to the increased matter density and power spectrum normalization the concentration mass relation in the Planck cosmology has a 20 percent higher normalization at redshift z=0 compared to WMAP cosmology. We confirm that CDM haloes are better described by the Einasto profile; for example, at scales near galaxy half-light radii CDM haloes have significantly steeper density profiles than implied by NFW fits. There is a scatter of 0.2 dex in the Einasto shape parameter at fixed halo mass, adding further to the diversity of CDM halo profiles. The evolution of the concentration mass relation in our simulations is not reproduced by any of the analytic models in the literature. We thus provide a simple fitting formula that accurately describes the evolution between redshifts z=5 to z=0 for both NFW and Einasto fits. Finally, the observed concentrations and halo masses of spiral galaxies, groups and clusters of galaxies at low redshifts are in good agreement with our simulations, suggesting only mild halo response to galaxy formation on these scales.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

preprint2014arXiv

The bulge-halo conspiracy in massive elliptical galaxies: implications for the stellar initial mass function and halo response to baryonic processes

Recent studies have shown that massive elliptical galaxies have total mass density profiles within an effective radius that can be approximated as ρ_{tot}\propto r^{-γ'}, with mean slope <γ'>=2.08 \pm 0.03 and scatter σ_γ'=0.16 \pm 0.02. The small scatter of the slope (known as the bulge-halo conspiracy) is not generic in LCDM based models and therefore contains information about the galaxy formation process. We compute the distribution of γ' for LCDM-based models that reproduce the observed correlations between stellar mass, velocity dispersion, and effective radius of early-type galaxies in the SDSS. The models have a range of stellar initial mass functions (IMFs) and dark halo responses to galaxy formation. The observed distribution of γ' is well reproduced by a model with cosmologically motivated but uncontracted dark matter haloes, and a Salpeter-type IMF. Other models are on average ruled out by the data, even though they may happen in individual cases. Models with adiabatic halo contraction (and lighter IMFs) predict too small values of γ'. Models with halo expansion, or mass-follows-light predict too high values of γ'. Our study shows that the non-homologous structure of massive early-type galaxies can be precisely reproduced by LCDM models if the IMF is not universal and if mechanisms such as feedback from active galactic nuclei, or dynamical friction, effectively on average counterbalance the contraction of the halo expected as a result of baryonic cooling.

preprint2014arXiv

The Distribution of Satellites Around Central Galaxies in a Cosmological Hydrodynamical Simulation

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

preprint2013arXiv

Dynamical Versus Stellar Masses in Compact Early-Type Galaxies: Further Evidence for Systematic Variation in the Stellar Initial Mass Function

Several independent lines of evidence suggest that the stellar initial mass function (IMF) in early-type galaxies becomes increasingly `bottom-heavy' with increasing galaxy mass and/or velocity dispersion, sigma. Here we consider evidence for IMF variation in a sample of relatively compact early-type galaxies drawn from the Sloan Digital Sky Survey. These galaxies are of sufficiently high stellar density that a dark halo likely makes a minor contribution to the total dynamical mass, Mdyn, within one effective radius. We fit our detailed stellar population synthesis models to the stacked absorption line spectra of these galaxies in bins of sigma and find evidence from IMF-sensitive spectral features for a bottom-heavy IMF at high sigma. We also apply simple `mass-follows-light' dynamical models to the same data and find that Mdyn is significantly higher than what would be expected if these galaxies were stellar dominated and had a universal Milky Way IMF. Adopting Mdyn M_* therefore implies that the IMF is `heavier' at high sigma. Most importantly, the quantitative amount of inferred IMF variation is very similar between the two techniques, agreeing to within 0.1 dex in mass. The agreement between two independent techniques, when applied to the same data, provides compelling evidence for systematic variation in the IMF as a function of early-type galaxy velocity dispersion. Any alternative explanations must reproduce both the results from dynamical and stellar population-based techniques.

preprint2013arXiv

The dependence of dark matter profiles on the stellar to halo mass ratio: a prediction for cusps vs cores

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, α in ρr^α, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of α on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how α varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with α decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes α<-0.8, approaching the NFW profile. This novel prediction for the dependence of α on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies.

preprint2013arXiv

The effect of Warm Dark Matter on galaxy properties: constraints from the stellar mass function and the Tully-Fisher relation

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the suppression of power on small scales, while no substantial difference is found at the high mass end. However these differences in the stellar mass function, vanish when different set of parameters are used to describe the (largely unknown) galaxy formation processes. We show that is possible to break this degeneracy between DM properties and the parameterization of baryonic physics by combining observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM models with a too warm candidate (m<0.75 keV) cannot simultaneously reproduce the stellar mass function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar mass function and the link between galaxies and dark matter haloes down to the very low-mass end can give very tight constraints on the nature of DM candidates.

preprint2013arXiv

The SWELLS Survey. VI. hierarchical inference of the initial mass functions of bulges and discs

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

preprint2013arXiv

Universal IMF vs dark halo response in early-type galaxies: breaking the degeneracy with the fundamental plane

We use the relations between aperture stellar velocity dispersion (σ_ap), stellar mass (M_sps), and galaxy size (R_e) for a sample of \sim 150,000 early-type galaxies from SDSS/DR7 to place constraints on the stellar initial mass function (IMF) and dark halo response to galaxy formation. We build LCDM based mass models that reproduce, by construction, the relations between galaxy size, light concentration and stellar mass, and use the spherical Jeans equations to predict σ_ap. Given our model assumptions (including those in the stellar population synthesis models), we find that reproducing the median σ_ap vs M_sps relation is not possible with {\it both} a universal IMF and a universal dark halo response. Significant departures from a universal IMF and/or dark halo response are required, but there is a degeneracy between these two solutions. We show that this degeneracy can be broken using the strength of the correlation between residuals of the velocity-mass (Δlog σ_ap) and size-mass (Δlog R_e) relations. The slope of this correlation, d_vr \equiv Δlog σ_ap/Δlog R_e, varies systematically with galaxy mass from d_vr \simeq -0.45 at M_sps \sim 10^{10}M_sun, to d_vr \simeq -0.15 at M_sps \sim 10^{11.6} M_sun. The virial fundamental plane (FP) has d_vr=-1/2, and thus we find the tilt of the observed FP is mass dependent. Reproducing this tilt requires {\it both} a non-universal IMF and a non-universal halo response. Our best model has mass-follows-light at low masses (Msps < 10^{11.2}M_sun) and unmodified NFW haloes at M_sps \sim 10^{11.5} M_sun. The stellar masses imply a mass dependent IMF which is "lighter" than Salpeter at low masses and "heavier" than Salpeter at high masses.

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.

preprint2012arXiv

Evidence for a non-universal stellar initial mass function in low-redshift high-density early-type galaxies

We determine an absolute calibration of stellar mass-to-light ratios for the densest \simeq 3% of early-type galaxies in the local universe (redshift z\simeq 0.08) from SDSS DR7. This sample of \sim 4000 galaxies has, assuming a Chabrier IMF, effective stellar surface densities, Sigma_e > 2500 M_sun/pc^2, stellar population synthesis (SPS) stellar masses log_10(M_sps/M_sun)<10.8, and aperture velocity dispersions of sigma_ap=168^{+37}_{-34} km/s (68% range). In contrast to typical early-type galaxies, we show that these dense early-type galaxies follow the virial fundamental plane, which suggests that mass-follows-light. With the additional assumption that any dark matter does not follow the light, the dynamical masses of dense galaxies provide a direct measurement of stellar masses. Our dynamical masses (M_dyn), obtained from the spherical Jeans equations, are only weakly sensitive to the choice of anisotropy (β) due to the relatively large aperture of the SDSS fiber for these galaxies: R_ap \simeq 1.5 R_e. Assuming isotropic orbits (β=0) we find a median log_{10} (M_dyn/M_sps) = 0.233 \pm 0.003, consistent with a Salpeter IMF, while more bottom heavy IMFs and standard Milky-Way IMFs are strongly disfavored. Our results are consistent with, but do not require, a dependence of the IMF on dynamical mass or velocity dispersion. We find evidence for a color dependence to the IMF such that redder galaxies have heavier IMFs with M_dyn/M_sps \propto (g-r)^{1.13\pm0.09}. This may reflect a more fundamental dependence of the IMF on the age or metallicity of a stellar population, or the density at which the stars formed.

preprint2012arXiv

The angular momentum of disc galaxies: implications for gas accretion, outflows, and dynamical friction

We combine constraints on the galaxy-dark matter connection with structural and dynamical scaling relations to investigate the angular momentum content of disc galaxies. For haloes with masses in the interval 10^{11.3} < M_vir/M_sun < 10^{12.7} we find that the galaxy spin parameters are independent of halo mass with <λ'_gal> = (J_gal/M_gal) / (\sqrt{2} R_vir V_vir) = 0.019^{+0.004}_{-0.003} (1sigma). This is significantly lower than for relaxed LCDM haloes, which have an average spin parameter <λ'_halo> = 0.031. The average ratio between the specific angular momentum of disk galaxies and their host dark matter haloes is therefore R_j = λ'_gal/λ'_halo = 0.61^{+0.13}_{-0.10}. This calls into question a standard assumption made in the majority of all (semi-analytical) models for (disc) galaxy formation, namely that R_j=1. Using simple disc formation models we show that it is particularly challenging to understand why R_j is independent of halo mass, while the galaxy formation efficiency (ε_GF, proportional to the ratio of galaxy mass to halo mass) reveals a strong halo mass dependence. We argue that the empirical scaling relations between ε_GF, R_j and halo mass require both feedback (i.e., galactic outflows) and angular momentum transfer from the baryons to the dark matter (i.e., dynamical friction). The efficiency of angular momentum loss need to decrease with increasing halo mass. Such a mass dependence may reflect a bias against forming stable discs in high mass, low spin haloes or a transition from cold-mode accretion in low mass haloes to hot-mode accretion at the massive end. However, current hydrodynamical simulations of galaxy formation, which should include these processes, seem unable to reproduce the empirical relation between ε_GF and R_j. We conclude that the angular momentum build-up of galactic discs remains poorly understood.

preprint2012arXiv

The baryonic Tully-Fisher relation and galactic outflows

Most of the baryons in the Universe are not in the form of stars and cold gas in galaxies. Galactic outflows driven by supernovae/stellar winds are the leading mechanism for explaining this fact. The scaling relation between galaxy mass and outer rotation velocity (also known as the baryonic Tully-Fisher relation, BTF) has recently been used as evidence against this viewpoint. We use a LCDM based semi-analytic disk galaxy formation model to investigate these claims. In our model, galaxies with less efficient star formation and higher gas fractions are more efficient at ejecting gas from galaxies. This is due to the fact that galaxies with less efficient star formation and higher gas fractions tend to live in dark matter haloes with lower circular velocities, from which less energy is required to escape the potential well. In our model the intrinsic scatter in the BTF is 0.15 dex, and mostly reflects scatter in dark halo concentration. The observed scatter, equal to 0.24 dex, is dominated by measurement errors. The best estimate for the intrinsic scatter is that it is less than 0.15 dex, and thus our LCDM based model (which does not include all possible sources of scatter) is only just consistent with this. In our model, gas rich galaxies, at fixed virial velocity (V_vir), with lower stellar masses have lower baryonic masses. This is consistent with the expectation that galaxies with lower stellar masses have had less energy available to drive an outflow. However, when the outer rotation velocity (V_flat) is used the correlation has the opposite sign, with a slope in agreement with observations. This is due to scatter in the relation between V_flat and V_vir. In summary, contrary to some previous claims, we show that basic features of the BTF are consistent with a LCDM based model in which the low efficiency of galaxy formation is determined by galactic outflows.

preprint2012arXiv

The DEEP2 Galaxy Redshift Survey: Design, Observations, Data Reduction, and Redshifts

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

preprint2012arXiv

The DEEP2 Galaxy Redshift Survey: The Voronoi-Delaunay Method catalog of galaxy groups

We present a public catalog of galaxy groups constructed from the spectroscopic sample of galaxies in the fourth data release from the DEEP2 Galaxy Redshift Survey, including the Extended Groth Strip (EGS). The catalog contains 1165 groups with two or more members in the EGS over the redshift range 0<z<1.5 and 1295 groups at z>0.6 in the rest of DEEP2. 25% of EGS galaxies and 14% of high-z DEEP2 galaxies are assigned to galaxy groups. The groups were detected using the Voronoi-Delaunay Method, after it has been optimized on mock DEEP2 catalogs following similar methods to those employed in Gerke et al. (2005). In the optimization effort, we have taken particular care to ensure that the mock catalogs resemble the data as closely as possible, and we have fine-tuned our methods separately on mocks constructed for the EGS and the rest of DEEP2. We have also probed the effect of the assumed cosmology on our inferred group-finding efficiency by performing our optimization on three different mock catalogs with different background cosmologies, finding large differences in the group-finding success we can achieve for these different mocks. Using the mock catalog whose background cosmology is most consistent with current data, we estimate that the DEEP2 group catalog is 72% complete and 61% pure (74% and 67% for the EGS) and that the group-finder correctly classifies 70% of galaxies that truly belong to groups, with an additional 46% of interloper galaxies contaminating the catalog (66% and 43% for the EGS). (Abridged)

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.

preprint2012arXiv

The Epoch of Disk Settling: z~1 to Now

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M* (M_Sun) < 10.7 over 0.2<z<1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar mass galaxies systematically decrease in disordered motions and increase in rotation velocity and potential well depth with time. Massive galaxies are the most well-ordered at all times examined, with higher rotation velocities and less disordered motions than less massive galaxies. We quantify disordered motions with an integrated gas velocity dispersion corrected for beam smearing (sigma_g). It is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of the galaxy sizes, sigma_g integrates over velocity gradients on large scales which can correspond to non-ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 1.2<z<3.8 and do not find any trends with redshift, likely for the most part because these datasets are biased toward the most highly star-forming systems. In summary, over the last ~8 billion years since z=1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being the most evolved at any time.

preprint2012arXiv

The SWELLS survey. III. Disfavouring "heavy" initial mass functions for spiral lens galaxies

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

preprint2012arXiv

The SWELLS survey. IV. Precision measurements of the stellar and dark matter distributions in a spiral lens galaxy

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

preprint2011arXiv

AGN Unification at z ~ 1: u - R Colors and Gradients in X-ray AGN Hosts

We present uncontaminated rest-frame u - R colors of 78 X-ray-selected AGN hosts at 0.5 < z < 1.5 in the Chandra Deep Fields measured with HST/ACS/NICMOS and VLT/ISAAC imaging. We also present spatially-resolved NUV - R color gradients for a subsample of AGN hosts imaged by HST/WFC3. Integrated, uncorrected photometry is not reliable for comparing the mean properties of soft and hard AGN host galaxies at z ~ 1 due to color contamination from point-source AGN emission. We use a cloning simulation to develop a calibration between concentration and this color contamination and use this to correct host galaxy colors. The mean u - R color of the unobscured/soft hosts beyond ~6 kpc is statistically equivalent to that of the obscured/hard hosts (the soft sources are 0.09 +/- 0.16 magnitudes bluer). Furthermore, the rest-frame V - J colors of the obscured and unobscured hosts beyond ~6 kpc are statistically equivalent, suggesting that the two populations have similar distributions of dust extinction. For the WFC3/IR sample, the mean NUV - R color gradients of unobscured and obscured sources differ by less than ~0.5 magnitudes for r > 1.1 kpc. These three observations imply that AGN obscuration is uncorrelated with the star formation rate beyond ~1 kpc. These observations favor a unification scenario for intermediate-luminosity AGNs in which obscuration is determined geometrically. Scenarios in which the majority of intermediate-luminosity AGN at z ~ 1 are undergoing rapid, galaxy-wide quenching due to AGN-driven feedback processes are disfavored.

preprint2011arXiv

An Investigation of Sloan Digital Sky Survey Imaging Data and Multi-Band Scaling Relations of Spiral Galaxies (with Dynamical Information)

We have compiled a sample of 3041 spiral galaxies with multi-band gri imaging from the Sloan Digital Sky Survey (SDSS) Data Release 7 and available galaxy rotational velocities derived from HI line widths. We compare the data products provided through the SDSS imaging pipeline with our own photometry of the SDSS images, and use the velocities (V) as an independent metric to determine ideal galaxy sizes (R) and luminosities (L). Our radial and luminosity parameters improve upon the SDSS DR7 Petrosian radii and luminosities through the use of isophotal fits to the galaxy images. This improvement is gauged via VL and RV relations whose respective scatters are reduced by ~8% and ~30% compared to similar relations built with SDSS parameters. The tightest VRL relations are obtained with the i-band radius, R235i, measured at 23.5 mag/arcsec^-2, and the luminosity L235i, measured within R235i. Our VRL scaling relations compare well, both in scatter and slope, with similar studies (such comparisons however depend sensitively on the nature and size of the compared samples). The typical slopes, b, and observed scatters, sigma, of the i-band VL, RL and RV relations are bVL=0.27+/-0.01, bRL=0.41+/-0.01, bRV=1.52+/-0.07, and sigmaVL=0.074, sigmaRL=0.071, sigmaRV=0.154 dex. Similar results for the SDSS g and r bands are also provided. Smaller scatters may be achieved for more pruned samples. We also compute scaling relations in terms of the baryonic mass (stars + gas), Mbar, ranging from 10^8.7 Msol to 10^11.6 Msol. Our baryonic velocity-mass (VM) relation has slope 0.29+/-0.01 and a measured scatter sigma_meas = 0.076 dex. While the observed VL and VM relations have comparable scatter, the stellar and baryonic VM relations may be intrinsically tighter, and thus potentially more fundamental, than other VL relations of spiral galaxies.

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.

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.

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.

preprint2011arXiv

The DEEP3 Galaxy Redshift Survey: The Impact of Environment on the Size Evolution of Massive Early-type Galaxies at Intermediate Redshift

Using data drawn from the DEEP2 and DEEP3 Galaxy Redshift Surveys, we investigate the relationship between the environment and the structure of galaxies residing on the red sequence at intermediate redshift. Within the massive (10 < log(M*/Msun) < 11) early-type population at 0.4 < z <1.2, we find a significant correlation between local galaxy overdensity (or environment) and galaxy size, such that early-type systems in higher-density regions tend to have larger effective radii (by ~0.5 kpc or 25% larger) than their counterparts of equal stellar mass and Sersic index in lower-density environments. This observed size-density relation is consistent with a model of galaxy formation in which the evolution of early-type systems at z < 2 is accelerated in high-density environments such as groups and clusters and in which dry, minor mergers (versus mechanisms such as quasar feedback) play a central role in the structural evolution of the massive, early-type galaxy population.

preprint2010arXiv

On the Evolution of the Velocity-Mass-Size Relations of Disk-Dominated Galaxies over the Past 10 Billion Years

We study the evolution of the scaling relations between maximum circular velocity, stellar mass and optical half-light radius of star-forming disk-dominated galaxies in the context of LCDM-based galaxy formation models. Using data from the literature combined with new data from the DEEP2 and AEGIS surveys we show that there is a consistent observational and theoretical picture for the evolution of these scaling relations from z\sim 2 to z=0. The evolution of the observed stellar scaling relations is weaker than that of the virial scaling relations of dark matter haloes, which can be reproduced, both qualitatively and quantitatively, with a simple, cosmologically-motivated model for disk evolution inside growing NFW dark matter haloes. In this model optical half-light radii are smaller, both at fixed stellar mass and maximum circular velocity, at higher redshifts. This model also predicts that the scaling relations between baryonic quantities evolve even more weakly than the corresponding stellar relations. We emphasize, though, that this weak evolution does not imply that individual galaxies evolve weakly. On the contrary, individual galaxies grow strongly in mass, size and velocity, but in such a way that they move largely along the scaling relations. Finally, recent observations have claimed surprisingly large sizes for a number of star-forming disk galaxies at z \sim 2, which has caused some authors to suggest that high redshift disk galaxies have abnormally high spin parameters. However, we argue that the disk scale lengths in question have been systematically overestimated by a factor \sim 2, and that there is an offset of a factor \sim 1.4 between Hαsizes and optical sizes. Taking these effects into account, there is no indication that star forming galaxies at high redshifts (z\sim 2) have abnormally high spin parameters.

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.

preprint2010arXiv

The Impact of Feedback on Disk Galaxy Scaling Relations

We use a disk formation model to study the effects of galactic outflows (a.k.a. feedback) on the rotation velocity - stellar mass - disk size, gas fraction - stellar mass, and gas phase metalicity - stellar mass scaling relations of disk galaxies. We show that models without outflows are unable to explain these scaling relations, having both the wrong slopes and normalization. The problem can be traced to the model galaxies having too many baryons. Models with outflows can solve this "over-cooling" problem by removing gas before it has time to turn into stars. Models with both momentum and energy driven winds can reproduce the observed scaling relations. However, these models predict different slopes which, with better observations, may be used to discriminate between these models.

preprint2010arXiv

The Kinematic Connection Between Galaxies and Dark Matter Haloes

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

preprint2009arXiv

Inspiralling Supermassive Black Holes: A New Signpost for Galaxy Mergers

We present a new technique for observationally identifying galaxy mergers spectroscopically rather than through host galaxy imaging. Our technique exploits the dynamics of supermassive black holes (SMBHs) powering active galactic nuclei (AGNs) in merger-remnant galaxies. Because structure in the universe is built up through galaxy mergers and nearly all galaxies host a central SMBH, some galaxies should possess two SMBHs near their centers as the result of a recent merger. These SMBHs spiral to the center of the resultant merger-remnant galaxy, and one or both of the SMBHs may power AGNs. Using the DEEP2 Galaxy Redshift Survey, we have examined 1881 red galaxies, of which 91 exhibit [O III] and Hbeta emission lines indicative of Seyfert 2 activity. Of these, 32 AGNs have [O III] emission-line redshifts significantly different from the redshifts of the host galaxies' stars, corresponding to velocity offsets of ~50 km/s to ~300 km/s. Two of these AGNs exhibit double-peaked [O III] emission lines, while the remaining 30 AGNs each exhibit a single set of velocity-offset [O III] emission lines. After exploring a variety of physical models for these velocity offsets, we argue that the most likely explanation is inspiralling SMBHs in merger-remnant galaxies. Based on this interpretation, we find that roughly half of the red galaxies hosting AGNs are also merger remnants, which implies that mergers may trigger AGN activity in red galaxies. The AGN velocity offsets we find imply a merger fraction of ~30% and a merger rate of ~3 mergers/Gyr for red galaxies at redshifts 0.34 < z < 0.82.

preprint2007arXiv

Concentration, Spin and Shape of Dark Matter Haloes: Scatter and the Dependence on Mass and Environment

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged)

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.