Source author record

Eric F. Bell

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph astro-ph.IM astro-ph.SR

Catalog footprint

What is connected

143works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

After The Fall: Resolving the Molecular Gas in Post-Starburst Galaxies

Post-starburst (PSB), or 'E+A', galaxies represent a rapid transitional phase between major, gas-rich mergers and gas-poor, quiescent early-type galaxies. Surprisingly, many PSBs have been shown to host a significant interstellar medium (ISM), despite theoretical predictions that the majority of star-forming gas should be expelled in AGN- or starburst-driven outflows. To-date, the resolved properties of this surviving ISM have remained unknown. We present high resolution ALMA continuum and CO(2$-$1) observations in six gas- and dust-rich PSBs, revealing for the first time the spatial and kinematic structure of their ISM on sub-kpc scales. We find extremely compact molecular reservoirs, with dust and gas surface densities rivaling those found in (ultra-)luminous infrared galaxies. We observe spatial and kinematic disturbances in all sources, with some also displaying disk-like kinematics. Estimates of the internal turbulent pressure in the gas exceed those of normal star-forming disks by at least 2 orders of magnitude, and rival the turbulent gas found in local interacting galaxies, such as the Antennae. Though the source of this high turbulent pressure remains uncertain, we suggest that the high incidence of tidal disruption events (TDEs) in PSBs could play a role. The star formation in these PSBs' turbulent central molecular reservoirs is suppressed, forming stars only 10% as efficiently as starburst galaxies with similar gas surface densities. "The fall" of star formation in these galaxies was not precipitated by complete gas expulsion or redistribution. Rather, this high-resolution view of PSBs' ISM indicates that star formation in their remaining compact gas reservoirs is suppressed by significant turbulent heating.

preprint2022arXiv

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics.

preprint2022arXiv

New Globular Cluster Candidates in the M81 group

The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of previously-discovered GCs to tune our selection criteria, finding that bright already-known GCs in the M81 group have sizes that are typically slightly larger than the Subaru PSF in our fields. In the optical bands, GCs appear to have colours that are only slightly different from stars. The inclusion of archival IRAC data yields dramatic improvements in colour separation, as the long wavelength baseline aids somewhat in the separation from stars and clearly separates GCs from many compact background galaxies. We show that some previously spectroscopically-identified GCs in the M81 group are instead foreground stars or background galaxies. GCs close to M82 have radial velocities suggesting that they fell into the M81 group along with M82. The overall M81 GC luminosity function is similar to the Milky Way and M31. M81's outer halo GCs are similar to the Milky Way in their metallicities and numbers, and much less numerous than M31's more metal-rich outer halo GC population. These properties reflect differences in the three galaxies' merger histories, highlighting the possibility of using outer halo GCs to trace merger history in larger samples of galaxies.

preprint2022arXiv

R2-D2: Roman and Rubin -- From Data to Discovery

The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community.

preprint2022arXiv

Relating the Diverse Merger Histories and Satellite Populations of Nearby Galaxies

We investigate whether the considerable diversity in the satellite populations of nearby Milky Way (MW)-mass galaxies is connected with the diversity in their host's merger histories. Analyzing 8 nearby galaxies with extensive observations of their satellite populations and stellar halos, we characterize each galaxy's merger history using the metric of its most dominant merger, $M_{\rm \star,Dom}$, defined as the greater of either its total accreted stellar mass or most massive current satellite. We find an unexpectedly tight relationship between these galaxies' number of $M_{V}\,{<}\,{-}9$ satellites within 150 kpc ($N_{\rm Sat}$) and $M_{\rm \star,Dom}$. This relationship remains even after accounting for differences in galaxy mass. Using the star formation and orbital histories of satellites around the MW and M81, we demonstrate that both likely evolved along the $M_{\rm\star,Dom}{-}N_{\rm Sat}$ relation during their current dominant mergers with the LMC and M82, respectively. We investigate the presence of this relation in galaxy formation models, including using the FIRE simulations to directly compare to the observations. We find no relation between $M_{\rm\star,Dom}$ and $N_{\rm Sat}$ in FIRE, and a universally large scatter in $N_{\rm Sat}$ with $M_{\rm \star,Dom}$ across models $-$ in direct contrast with the tightness of the empirical relation. This acute difference in the observed and predicted scaling relation between two fundamental galaxy properties signals that current simulations do not sufficiently reproduce diverse merger histories and their effects on satellite populations. Explaining the emergence of this relation is therefore essential for obtaining a complete understanding of galaxy formation.

preprint2022arXiv

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $σ'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $σ'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a Sérsic profile-based $R$ estimate, then a Sérsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex.

preprint2022arXiv

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center of M33 along the major and semi-major axes, respectively. We divide the PHATTER optical survey into 2005 regions that measure 24 arcsec, $\sim$100 pc, on a side and fit color magnitude diagrams for each region individually to measure the spatially resolved SFH of M33 within the PHATTER footprint. There are significant fluctuations in the SFH on small spatial scales and also galaxy-wide scales that we measure back to about 630 Myr ago. We observe a more flocculent spiral structure in stellar populations younger than about 80 Myr, while the structure of the older stellar populations is dominated by two spiral arms. We also observe a bar in the center of M33, which dominates at ages older than about 80 Myr. Finally, we find that the mean star formation rate (SFR) over the last 100 Myr within the PHATTER footprint is 0.32$\pm$0.02 M$_{\odot}$ yr$^{-1}$. We measure a current SFR (over the last 10 Myr) of 0.20$\pm$0.03 M$_{\odot}$ yr$^{-1}$. This SFR is slightly higher than previous measurements from broadband estimates, when scaled to account for the fraction of the D25 area covered by the PHATTER survey footprint.

preprint2022arXiv

Uncertainties associated with the backward integration of dwarf satellites using simple parametric potentials

In order to backward integrate the orbits of Milky Way (MW) dwarf galaxies, much effort has been invested in recent years to constrain their initial phase-space coordinates. Yet equally important are the assumptions on the potential that the dwarf galaxies experience over time, especially given the fact that the MW is currently accreting the Large Magellanic Cloud (LMC). In this work, using a dark matter-only zoom-in simulation, we test whether the use of common parametric forms of the potential is adequate to successfully backward integrate the orbits of the subhaloes from their present-day positions. We parametrise the recovered orbits and compare them with those from the simulations. We find that simple symmetric parametric forms of the potential fail to capture the complexities and the inhomogeneities of the true potential experienced by the subhaloes. More specifically, modelling a recent massive accretion like that of the LMC as a sum of two spherical parametric potentials leads to substantial errors in the recovered parameters of the orbits. These errors rival those caused due to a) a 30\% uncertainty in the virial mass of the MW and b) not modelling the potential of the recently accreted massive satellite. Our work suggests that i) the uncertainties in the parameters of the recovered orbits of some MW dwarfs may be under-estimated and that ii) researchers should characterise the uncertainties inherent to their choice of integration techniques and assumptions of the potential against cosmological zoom-in simulations of the MW, which include a recently-accreted LMC.

preprint2021arXiv

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra.

preprint2020arXiv

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $Hδ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $Hδ_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ Å/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ Å$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius.

preprint2020arXiv

Is NGC 300 a pure exponential disk galaxy?

NGC 300 is a low-mass disk galaxy in the Sculptor group. In the literature, it has been identified as a pure exponential disk galaxy, as its luminosity profile could be well fitted with a single exponential law over many disk scale lengths (Type I). We investigate the stellar luminosity distribution of NGC 300 using $Hubble$ $Space$ $Telescope$ (HST) archive data, reaching farther and deeper than any other previous studies. Color magnitude diagrams show a significant population of old red giant branch (RGB) stars in all fields out to $R\sim19$ kpc ($32'$), as well as younger populations in the inner regions. We construct the density profiles of the young, intermediate-aged, and old stellar populations. We find two clear breaks in the density profiles of the old RGB and intermediate-aged stars: one down-bending (Type II) at $R\sim5.9$ kpc, and another up-bending (Type III) at $R\sim8.3$ kpc. Moreover, the old RGB stars exhibit a negative radial color gradient with an up-bending at $R\sim8$~kpc, beyond which the stellar populations are uniformly old ($>$7~Gyr) and metal-poor ($\rm[Fe/H] = -1.6^{+0.2}_{-0.4}$ dex). The outer stellar component at $R\gtrapprox8$ kpc is, therefore, well separated from the inner disk in terms of the stellar density and stellar populations. While our results cast doubt on the currently established wisdom that NGC\,300 is a pure exponential disk galaxy, a more detailed survey should be carried out to identify the outskirts as either a disk or a stellar halo.

preprint2020arXiv

Mass-to-Light Ratios of Spatially Resolved Stellar Populations in M31

A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial variations in $M_\star/L$ and their dependence on color, star formation history, and dust across the disk of M31, using a map of $M^\mathrm{CMD}_\star$ derived from color-magnitude diagrams of resolved stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. First, we find comparable scatter in $M_\star/L$ for the optical and mid-IR, contrary to the common idea that $M_\star/L$ is less variable in the IR. Second, we confirm that $M_\star/L$ is correlated with color for both the optical and mid-IR and report color vs. $M_\star/L$ relations (CMLRs) in M31 for filters used in the Sloan Digital Sky Survey (SDSS) and Widefield Infrared Survey Explorer (WISE). Third, we show that the CMLR residuals correlate with recent SFH, such that quiescent regions are offset to higher $M_\star/L$ than star-forming regions at a fixed color. The mid-IR CMLR, however, is not linear due to the high scatter of $M_\star/L$ in star-forming regions. Finally, we find a flatter optical CMLR than any SPS-based CMLRs in the literature. We show this is an effect of dust geometry, which is typically neglected but should be accounted for when using optical data to map $M_\star/L$.

preprint2020arXiv

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

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

preprint2020arXiv

SMASHing the low surface brightness SMC

The periphery of the Small Magellanic Cloud (SMC) can unlock important information regarding galaxy formation and evolution in interacting systems. Here, we present a detailed study of the extended stellar structure of the SMC using deep colour-magnitude diagrams (CMDs), obtained as part of the Survey of the MAgellanic Stellar History (SMASH). Special care was taken in the decontamination of our data from MW foreground stars, including from foreground globular clusters NGC 362 and 47 Tuc. We derived the SMC surface brightness using a ``conservative'' approach from which we calculated the general parameters of the SMC, finding a staggered surface brightness profile. We also traced the fainter outskirts by constructing a stellar density profile. This approach, based on stellar counts of the oldest main sequence turn-off (MSTO) stars, uncovered a tidally disrupted stellar feature that reaches as far out as 12 degrees from the SMC centre. We also serendipitously found a faint feature of unknown origin located at $\sim 14$ degrees from the centre of the SMC and that we tentatively associated to a more distant structure. We compared our results to in-house simulations of a $1\times10^{9} M_\odot$ SMC, finding that its elliptical shape can be explained by its tidal disruption under the combined presence of the MW and the LMC. Finally, we found that the older stellar populations show a smooth profile while the younger component presents a jump in the density followed by a flat profile, confirming the heavily disturbed nature of the SMC.

preprint2020arXiv

The Second Data Release of the Survey of the MAgellanic Stellar History (SMASH)

The Large and Small Magellanic Clouds (LMC and SMC) are the largest satellite galaxies of the Milky Way and close enough to allow for a detailed exploration of their structure and formation history. The Survey of the MAgellanic Stellar History (SMASH) is a community Dark Energy Camera (DECam) survey of the Magellanic Clouds using $\sim$50 nights to sample over $\sim$2400 deg$^2$ centered on the Clouds at $\sim$20% filling factor (but with contiguous coverage in the central regions) and to depths of $\sim$24th mag in $ugriz$. The primary goals of SMASH are to map out the extended stellar peripheries of the Clouds and uncover their complicated interaction and accretion history as well as to derive spatially-resolved star formation histories of the central regions and create a "movie" of their past star formation. Here we announce the second SMASH public data release (DR2), which contains all 197 fully-calibrated DECam fields including the main body fields in the central regions. The DR2 data are available through the Astro Data Lab hosted by the NSF's National Optical-Infrared Astronomy Research Laboratory. We highlight three science cases that make use of the SMASH DR2 data and will be published in the future: (1) preliminary star formation histories of the LMC; (2) the search for Magellanic star clusters using citizen scientists; and, (3) photometric metallicities of Magellanic Cloud stars using the DECam $u$-band.

preprint2020arXiv

Tracing the anemic stellar halo of M101

Models of galaxy formation in a cosmological context predict that massive disk galaxies should have structured extended stellar halos. Recent studies in integrated light, however, report a few galaxies, including the nearby disk galaxy M101, that have no measurable stellar halos to the detection limit. We aim to quantify the stellar content and structure of M101's outskirts by resolving its stars. We present the photometry of its stars based on deep F606W and F814W images taken with Hubble Space Telescope as part of the GHOSTS survey. The constructed CMDs of stars reach down to two magnitudes below the tip of the red giant branch. We derived radial number density profiles of the bright red giant branch (RGB) stars. The mean color of the RGB stars at $R \sim$ 40 -- 60 kpc is similar to those of metal-poor globular clusters in the Milky Way. We also derived radial surface brightness profiles using the public image data provided by the Dragonfly team. Both the radial number density and surface brightness profiles were converted to radial mass density profiles and combined. We find that the mass density profiles show a weak upturn at the very outer region, where surface brightness is as faint as $μ_g\approx 34$ mag arcsec$^{-1}$. An exponential disk + power-law halo model on the mass density profiles finds the total stellar halo mass of $M_{halo}=8.2_{-2.2}^{+3.5}\times 10^7M_\odot$. The total stellar halo mass does not exceed $M_{halo} = 3.2 \times 10^8$ $M_{\odot}$ when strongly truncated disk models are considered. Combining the halo mass with the total stellar mass of M101, we obtain the stellar halo mass fraction of $M_{halo}/M_{gal} = 0.20_{-0.08}^{+0.10}\%$ with an upper limit of 0.78\%. We compare the halo properties of M101 with those of six GHOSTS survey galaxies as well as the Milky Way and M31 and find that M101 has an anemic stellar halo.

preprint2019arXiv

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

preprint2019arXiv

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(Hδ) \geq 4$Å) at $z\sim0.8$ drawn from DR2 of the LEGA-C survey to test the hypothesis that these galaxies experienced compact, central starbursts before quenching. These recently quenched galaxies, usually referred to as post-starburst galaxies, span a wide range of colors and we find a clear correlation between color and half-light radius, such that bluer galaxies are smaller. We build simple toy models to explain this correlation: a normal star-forming disk plus a central, compact starburst component. Bursts with exponential decay timescale of $\sim$~100 Myr that produce $\sim10\%$ to more than 100\% of the pre-existing masses can reproduce the observed correlation. More significant bursts also produce bluer and smaller descendants. Our findings imply that when galaxies shut down star formation rapidly, they generally had experienced compact, starburst events and that the large, observed spread in sizes and colors mostly reflects a variety of burst strengths. Recently quenched galaxies should have younger stellar ages in the centers; multi-wavelength data with high spatial resolution are required to reveal the age gradient. Highly dissipative processes should be responsible for this type of formation history. While determining the mechanisms for individual galaxies is challenging, some recently quenched galaxies show signs of gravitational interactions, suggesting that mergers are likely an important mechanism in triggering the rapid shut-down of star-formation activities at $z\sim0.8$.

preprint2019arXiv

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

preprint2016arXiv

A Synoptic Map of Halo Substructures from the Pan-STARRS1 3π Survey

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.

preprint2016arXiv

ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $μ$Jy beam$^{-1}$ and 12.7 $μ$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

preprint2016arXiv

Detailed Chemical Abundances in the r-Process-Rich Ultra-Faint Dwarf Galaxy Reticulum 2

The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of r-process material in Ret 2 (mean [Eu/Fe]=+1.69+/-0.05) found in three of these stars (mean [Fe/H]=-2.88+/-0.10). The abundances closely match the r-process pattern found in the well-studied metal-poor halo star CS22892-052. Such r-process-enhanced stars have not been found in any other ultra-faint dwarf galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 ([Fe/H]=-3.42+/-0.20) contains only trace amounts of Sr ([Sr/Fe]=-1.73+/-0.43) and no detectable heavier elements. One r-process enhanced star is also enhanced in C (natal [C/Fe]=+1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and r-process enhancements are decoupled. The r-process-deficient star is enhanced in Mg ([Mg/Fe]=+0.81+/-0.14), and the other three stars show normal levels of alpha-enhancement (mean [Mg/Fe]=+0.34+/-0.03). The abundances of other alpha and Fe-group elements closely resemble those in ultra-faint dwarf galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large r-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae.

preprint2016arXiv

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

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

preprint2016arXiv

First results from the MADCASH Survey: A Faint Dwarf Galaxy Companion to the Low Mass Spiral Galaxy NGC 2403 at 3.2 Mpc

We report the discovery of the faintest known dwarf galaxy satellite of an LMC stellar-mass host beyond the Local Group, based on deep imaging with Subaru/Hyper Suprime-Cam. MADCASH J074238+652501-dw lies $\sim$35 kpc in projection from NGC 2403, a dwarf spiral galaxy at $D$$\approx$3.2 Mpc. This new dwarf has $M_{g} = -7.4\pm0.4$ and a half-light radius of $168\pm70$ pc, at the calculated distance of $3.39\pm0.41$ Mpc. The color-magnitude diagram reveals no evidence of young stellar populations, suggesting that MADCASH J074238+652501-dw is an old, metal-poor dwarf similar to low luminosity dwarfs in the Local Group. The lack of either detected HI gas ($M_{\rm HI}/L_{V} < 0.69 M_\odot/L_\odot$, based on Green Bank Telescope observations) or $GALEX$ NUV/FUV flux enhancement is consistent with a lack of young stars. This is the first result from the MADCASH (Magellanic Analog Dwarf Companions And Stellar Halos) survey, which is conducting a census of the stellar substructure and faint satellites in the halos of Local Volume LMC analogs via resolved stellar populations. Models predict a total of $\sim$4-10 satellites at least as massive as MADCASH J074238+652501-dw around a host with the mass of NGC 2403, with 2-3 within our field of view, slightly more than the one such satellite observed in our footprint.

preprint2016arXiv

Is Draco II one of the faintest dwarf galaxies? First study from Keck/DEIMOS spectroscopy

We present the first spectroscopic analysis of the faint and compact stellar system Draco II (Dra II, M_V=-2.9+/-0.8, r_h=19^{+8}_{-6} pc), recently discovered in the Pan-STARRS1 3π survey. The observations, conducted with DEIMOS on the Keck II telescope, establish some of its basic characteristics: the velocity data reveal a narrow peak with 9 member stars at a systemic heliocentric velocity <v_r>=-347.6^{+1.7}_{-1.8} km/s, thereby confirming Dra II is a satellite of the Milky Way; we infer a velocity dispersion with σ_{vr}=2.9+/-2.1 km/s (<8.4 km/s at the 95% confidence level), which implies log_{10}(M_{1/2})=5.5^{+0.4}_{-0.6} and log_{10}((M/L)_{1/2})=2.7^{+0.5}_{-0.8}, in Solar units; furthermore, very weak Calcium triplet lines in the spectra of the high signal-to-noise member stars imply [Fe/H]<-2.1, whilst variations in the line strengths of two stars with similar colours and magnitudes suggest a metallicity spread in Dra II. These new data cannot clearly discriminate whether Draco II is a star cluster or amongst the faintest, most compact, and closest dwarf galaxies. However, the sum of the three --- individually inconclusive --- pieces of evidence presented here, seems to favour the dwarf galaxy interpretation.

preprint2016arXiv

Mapping the Monoceros Ring in 3D with Pan-STARRS1

Using the Pan-STARRS1 survey, we derive limiting magnitude, spatial completeness and density maps that we use to probe the three dimensional structure and estimate the stellar mass of the so-called Monoceros Ring. The Monoceros Ring is an enormous and complex stellar sub-structure in the outer Milky Way disk. It is most visible across the large Galactic Anticenter region, 120 < l < 240 degrees, -30 < b < +40 degrees. We estimate its stellar mass density profile along every line of sight in 2 X 2 degree pixels over the entire 30,000 square degree Pan-STARRS1 survey using the previously developed MATCH software. By parsing this distribution into a radially smooth component and the Monoceros Ring, we obtain its mass and distance from the Sun along each relevant line of sight. The Monoceros Ring is significantly closer to us in the South (6 kpc) than in the North (9 kpc). We also create 2D cross sections parallel to the Galactic plane that show 135 degrees of the Monoceros Ring in the South and 170 degrees of the Monoceros Ring in the North. We show that the Northern and Southern structures are also roughly concentric circles, suggesting that they may be a wave rippling from a common origin. Excluding the Galactic plane, we observe an excess stellar mass of 4 million solar masses across 120 < l < 240 degrees. If we interpolate across the Galactic plane, we estimate that this region contains 8 million solar masses. If we assume (somewhat boldly) that the Monoceros Ring is a set of two Galactocentric rings, its total stellar mass is 60 million solar masses. Finally, if we assume that it is a set of two circles centered at a point 4 kpc from the Galactic center in the anti-central direction, as our data suggests, we estimate its stellar mass to be 40 million solar masses.

preprint2016arXiv

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

preprint2016arXiv

Panchromatic Hubble Andromeda Treasury XVI. Star Cluster Formation Efficiency and the Clustered Fraction of Young Stars

We use the Panchromatic Hubble Andromeda Treasury (PHAT) survey dataset to perform spatially resolved measurements of star cluster formation efficiency ($Γ$), the fraction of stellar mass formed in long-lived star clusters. We use robust star formation history and cluster parameter constraints, obtained through color-magnitude diagram analysis of resolved stellar populations, to study Andromeda's cluster and field populations over the last $\sim$300 Myr. We measure $Γ$ of 4-8% for young, 10-100 Myr old populations in M31. We find that cluster formation efficiency varies systematically across the M31 disk, consistent with variations in mid-plane pressure. These $Γ$ measurements expand the range of well-studied galactic environments, providing precise constraints in an HI-dominated, low intensity star formation environment. Spatially resolved results from M31 are broadly consistent with previous trends observed on galaxy-integrated scales, where $Γ$ increases with increasing star formation rate surface density ($Σ_{\mathrm{SFR}}$). However, we can explain observed scatter in the relation and attain better agreement between observations and theoretical models if we account for environmental variations in gas depletion time ($τ_{\mathrm{dep}}$) when modeling $Γ$, accounting for the qualitative shift in star formation behavior when transitioning from a H$_2$-dominated to a HI-dominated interstellar medium. We also demonstrate that $Γ$ measurements in high $Σ_{\mathrm{SFR}}$ starburst systems are well-explained by $τ_{\mathrm{dep}}$-dependent fiducial $Γ$ models.

preprint2016arXiv

Quiescence correlates strongly with directly-measured black hole mass in central galaxies

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly-measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous AGN feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

preprint2016arXiv

The Bursty Star Formation Histories of Low-mass Galaxies at $0.4<z<1$ Revealed by Star Formation Rates Measured from H$β$ and FUV

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$β$ and FUV (1500 Å) (H$β$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$β$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$β$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$β$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$β$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The Diversity of Growth Histories of Milky Way-mass Galaxies

We use the semi-analytic model developed by Henriques et al. (2015) to explore the origin of star formation history diversity for galaxies that lie at the centre of their dark matter haloes and have present-day stellar masses in the range 5-8 $\times$ 10$^{10}$ M$_{\odot}$, similar to that of the Milky Way. In this model, quenching is the dominant physical mechanism for introducing scatter in the growth histories of these galaxies. We find that present-day quiescent galaxies have a larger variety of growth histories than star-formers since they underwent 'staggered quenching' - a term describing the correlation between the time of quenching and present-day halo mass. While halo mass correlates broadly with quiescence, we find that quiescence is primarily a function of black hole mass, where galaxies quench when heating from their active galactic nuclei becomes sufficient to offset the redshift-dependent cooling rate. In this model, the emergence of a prominent quiescent population is the main process that flattens the stellar mass-halo mass relation at mass scales at or above that of the Milky Way.

preprint2016arXiv

The GHOSTS survey. II. The diversity of Halo Color and Metallicity Profiles of Massive Disk Galaxies

We study the stellar halo color properties of six nearby massive highly inclined disk galaxies using Hubble Space Telescope Advanced Camera for Surveys and Wide Field Camera 3 observations in both F606W and F814W filters from the GHOSTS survey. The observed fields, placed both along the minor and major axis of each galaxy, probe the stellar outskirts out to projected distances of ~ 50-70 kpc from their galactic centre along the minor axis. The 50% completeness levels of the color magnitude diagrams are typically at two mag below the tip of the red giant branch. We find that all galaxies have extended stellar halos out to ~ 50 kpc and two out to ~ 70 kpc. We determined the halo color distribution and color profile for each galaxy using the median colors of stars in the RGB. Within each galaxy we find variations in the median colors as a function of radius which likely indicates population variations, reflecting that their outskirts were built from several small accreted objects. We find that half of the galaxies (NGC 0891, NGC 4565, and NGC 7814) present a clear negative color gradient, reflecting a declining metallicity in their halos; the other have no significant color or population gradient. In addition, notwithstanding the modest sample size of galaxies, there is no strong correlation between their halo color/metallicity or gradient with galaxy's properties such as rotational velocity or stellar mass. The diversity in halo color profiles observed in the GHOSTS galaxies qualitatively supports the predicted galaxy-to-galaxy scatter in halo stellar properties; a consequence of the stochasticity inherent in the assembling history of galaxies.

preprint2016arXiv

The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

preprint2016arXiv

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

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

preprint2016arXiv

The Stellar Density Profile of the Distant Galactic Halo

We use extensive gravity-sensitive DDO 51 photometry over 5100 square degrees, combined with SDSS broadband photometry, to select a catalog of $\sim 4,000$ giant stars covering a large fraction of the high Galactic latitude sky and reaching out to $\sim 80$ kpc in the Galactic halo. This sample of bright and unbiased tracers enables us to measure the radial profile and 3D structure of the stellar halo to large distance which had previously only been measured with sparse tracers or small samples. Using population synthesis models to reproduce the observed giant star luminosity function, we find that the halo maintains a $r^{-3.5}$ profile from $30$ to $80$ kpc with no signs of a truncation or sharp break over this range. The radial profile measurement is largely insensitive to individual halo substructure components, but we find that attempting to measure the shape of the halo is overwhelmed by the Sagittarius stream such that no ellipsoidal shape is a satisfactory description in this region. These measurements allow us to begin placing the Milky Way in context with the growing sample of external galaxies where similar halo profile measurements are available, with the goal of further linking the properties of stellar halos to the accretion histories that formed them.

preprint2016arXiv

Triangulum II: a very metal-poor and dynamically hot stellar system

We present a study of the recently discovered compact stellar system Triangulum II. From observations conducted with the DEIMOS spectrograph on Keck II, we obtained spectra for 13 member stars that follow the CMD features of this very faint stellar system and include two bright red giant branch stars. Tri II has a very negative radial velocity (<v_r>=-383.7^{+3.0}_{-3.3} km/s) that translates to <v_{r,gsr}> ~ -264 km/s and confirms it is a Milky Way satellite. We show that, despite the small data set, there is evidence that Tri II has complex internal kinematics. Its radial velocity dispersion increases from 4.4^{+2.8}_{-2.0} km/s in the central 2' to 14.1^{+5.8}_{-4.2} km/s outwards. The velocity dispersion of the full sample is inferred to be σ_{vr}=9.9^{+3.2}_{-2.2} km/s. From the two bright RGB member stars we measure an average metallicity <[Fe/H]>=-2.6+/-0.2, placing Tri II among the most metal-poor Milky Way dwarf galaxies. In addition, the spectra of the fainter member stars exhibit differences in their line widths that could be the indication of a metallicity dispersion in the system. All these properties paint a complex picture for Tri II, whose nature and current state are largely speculative. The inferred metallicity properties of the system however lead us to favor a scenario in which Tri II is a dwarf galaxy that is either disrupting or embedded in a stellar stream.

preprint2015arXiv

A clear age-velocity dispersion correlation in Andromeda's stellar disk

The stellar kinematics of galactic disks are key to constraining disk formation and evolution processes. In this paper, for the first time, we measure the stellar age-velocity dispersion correlation in the inner 20 kpc (3.5 disk scale lengths) of M31 and show that it is dramatically different from that in the Milky Way. We use optical Hubble Space Telescope/Advanced Camera for Surveys photometry of 5800 individual stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and Keck/DEIMOS radial velocity measurements of the same stars from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We show that the average line-of-sight velocity dispersion is a steadily increasing function of stellar age exterior to R=10 kpc, increasing from 30 km/s for the young upper main sequence stars to 90 km/s for the old red giant branch stars. This monotonic increase implies that a continuous or recurring process contributed to the evolution of the disk. Both the slope and normalization of the dispersion vs. age relation are significantly larger than in the Milky Way, allowing for the possibility that the disk of M31 has had a more violent history than the disk of the Milky Way, more in line with cosmological predictions. We also find evidence for an inhomogeneous distribution of stars from a second kinematical component in addition to the dominant disk component. One of the largest and hottest high-dispersion patches is present in all age bins, and may be the signature of the end of the long bar.

preprint2015arXiv

A Deep Study of the Dwarf Satellites Andromeda XXVIII & Andromeda XXIX

We present the results of a deep study of the isolated dwarf galaxies Andromeda XXVIII and Andromeda XXIX with Gemini/GMOS and Keck/DEIMOS. Both galaxies are shown to host old, metal-poor stellar populations with no detectable recent star formation, conclusively identifying both of them as dwarf spheroidal galaxies (dSphs). And XXVIII exhibits a complex horizontal branch morphology, which is suggestive of metallicity enrichment and thus an extended period of star formation in the past. Decomposing the horizontal branch into blue (metal poor, assumed to be older) and red (relatively more metal rich, assumed to be younger) populations shows that the metal rich are also more spatially concentrated in the center of the galaxy. We use spectroscopic measurements of the Calcium triplet, combined with the improved precision of the Gemini photometry, to measure the metallicity of the galaxies, confirming the metallicity spread and showing that they both lie on the luminosity-metallicity relation for dwarf satellites. Taken together, the galaxies exhibit largely typical properties for dSphs despite their significant distances from M31. These dwarfs thus place particularly significant constraints on models of dSph formation involving environmental processes such as tidal or ram pressure stripping. Such models must be able to completely transform the two galaxies into dSphs in no more than two pericentric passages around M31, while maintaining a significant stellar populations gradient. Reproducing these features is a prime requirement for models of dSph formation to demonstrate not just the plausibility of environmental transformation but the capability of accurately recreating real dSphs.

preprint2015arXiv

A Global Star Forming Episode in M31 2-4 Gyr Ago

We have identified a major global enhancement of star formation in the inner M31 disk that occurred between 2-4 Gyr ago, producing $\sim$60% of the stellar mass formed in the past 5 Gyr. The presence of this episode in the inner disk was discovered by modeling the optical resolved star color-magnitude diagrams of low extinction regions in the main disk of M31 (3$<$R$<$20 kpc) as part of the Panchromatic Hubble Andromeda Treasury. This measurement confirms and extends recent measurements of a widespread star formation enhancement of similar age in the outer disk, suggesting that this burst was both massive and global. Following the galaxy-wide burst, the star formation rate of M31 has significantly declined. We briefly discuss possible causes for these features of the M31 evolutionary history, including interactions with M32, M33 and/or a merger.

preprint2015arXiv

A New Faint Milky Way Satellite Discovered in the Pan-STARRS1 3 pi Survey

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Boötes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

preprint2015arXiv

And yet it moves: The dangers of artificially fixing the Milky Way center of mass in the presence of a massive Large Magellanic Cloud

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC's orbital period and apocentric distance derived assuming a fixed Milky Way are significantly shortened for models where the Milky Way is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the Milky Way center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km/s. Furthermore, we show that the gravitational pull of the LMC and response of the Milky Way are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the Milky Way around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the Milky Way and Local Group.

preprint2015arXiv

Are Compton-Thick AGN the Missing Link Between Mergers and Black Hole Growth?

We examine the host morphologies of heavily obscured active galactic nuclei (AGN) at $z\sim1$ to test whether obscured supermassive black hole growth at this epoch is preferentially linked to galaxy mergers. Our sample consists of 154 obscured AGN with $N_{\rm H}>10^{23.5}$ cm$^{-2}$ and $z<1.5$. Using visual classifications, we compare the morphologies of these AGN to control samples of moderately obscured ($10^{22}$ cm$^{-2}$ $<N_{\rm H}< 10^{23.5}$ cm$^{-2}$) and unobscured ($N_{\rm H}<10^{22}$ cm$^{-2}$) AGN. These control AGN are matched in redshift and intrinsic X-ray luminosity to our heavily obscured AGN. We find that heavily obscured AGN at z~1 are twice as likely to be hosted by late-type galaxies relative to unobscured AGN ($65.3^{+4.1}_{-4.6}\%$ vs $34.5^{+2.9}_{-2.7}\%$) and three times as likely to exhibit merger or interaction signatures ($21.5^{+4.2}_{-3.3}\%$ vs $7.8^{+1.9}_{-1.3}\%$). The increased merger fraction is significant at the 3.8$σ$ level. We also find that the incidence of point-like morphologies is inversely proportional to obscuration. If we exclude all point sources and consider only extended hosts, we find the correlation between merger fraction and obscuration is still evident, however at a reduced statistical significance ($2.5σ$). The fact that we observe a different disk/spheroid fraction versus obscuration indicates that viewing angle cannot be the only thing differentiating our three AGN samples, as a simple unification model would suggest. The increased fraction of disturbed morphologies with obscuration supports an evolutionary scenario, in which Compton-thick AGN are a distinct phase of obscured SMBH growth following a merger/interaction event. Our findings also suggest that some of the merger-triggered SMBH growth predicted by recent AGN fueling models may be hidden among the heavily obscured, Compton-thick population.

preprint2015arXiv

CANDELS Visual Classifications: Scheme, Data Release, and First Results

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

preprint2015arXiv

Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3

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

preprint2015arXiv

Finding, characterizing and classifying variable sources in multi-epoch sky surveys: QSOs and RR Lyrae in PS1 3$π$ data

In area and depth, the Pan-STARRS1 (PS1) 3$π$ survey is unique among many-epoch, multi-band surveys and has enormous potential for all-sky identification of variable sources. PS1 has observed the sky typically seven times in each of its five bands ($grizy$) over 3.5 years, but unlike SDSS not simultaneously across the bands. Here we develop a new approach for quantifying statistical properties of non-simultaneous, sparse, multi-color lightcurves through light-curve structure functions, effectively turning PS1 into a $\sim 35$-epoch survey. We use this approach to estimate variability amplitudes and timescales $(ω_r, τ)$ for all point-sources brighter than $r_{\mathrm{P1}}=21.5$ mag in the survey. With PS1 data on SDSS Stripe 82 as ``ground truth", we use a Random Forest Classifier to identify QSOs and RR Lyrae based on their variability and their mean PS1 and WISE colors. We find that, aside from the Galactic plane, QSO and RR Lyrae samples of purity $\sim$75\% and completeness $\sim$92\% can be selected. On this basis we have identified a sample of $\sim 1,000,000$ QSO candidates, as well as an unprecedentedly large and deep sample of $\sim$150,000 RR Lyrae candidates with distances from $\sim$10 kpc to $\sim$120 kpc. Within the Draco dwarf spheroidal, we demonstrate a distance precision of 6\% for RR Lyrae candidates. We provide a catalog of all likely variable point sources and likely QSOs in PS1, a total of $25.8\times 10^6$ sources.

preprint2015arXiv

Hydra II: a faint and compact Milky Way dwarf galaxy found in the Survey of the Magellanic Stellar History

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

preprint2015arXiv

Linking the Structural Properties of Galaxies and their Star Formation Histories with STAGES

We study the links between star formation history and structure for a large mass-selected galaxy sample at 0.05 < z_phot < 0.30. The galaxies inhabit a very broad range of environments, from cluster cores to the field. Using HST images, we quantify their structure following Hoyos et al. (2012), and divide them into disturbed and undisturbed. We also visually identify mergers. Additionally, we provide a quantitative measure of the degree of disturbance for each galaxy ("roughness"). The majority of elliptical and lenticular galaxies have relaxed structure, showing no signs of ongoing star formation. Structurally-disturbed galaxies, which tend to avoid the lowest-density regions, have higher star-formation activity and younger stellar populations than undisturbed systems. Cluster spirals with reduced/quenched star formation have somewhat less disturbed morphologies than spirals with "normal" star-formation activity, suggesting that these "passive" spirals have started their morphological transformation into S0s. Visually identified mergers and galaxies not identified as mergers but with similar roughness have similar specific star formation rates and stellar ages. The degree of enhanced star formation is thus linked to the degree of structural disturbance, regardless of whether it is caused by major mergers or not. This suggests that merging galaxies are not special in terms of their higher-than-normal star-formation activity. Any physical process that produces "roughness", or regions of enhanced luminosity density, will increase the star-formation activity in a galaxy with similar efficiency. An alternative explanation is that star formation episodes increase the galaxies' roughness similarly, regardless of whether they are merger-induced or not.

preprint2015arXiv

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud.

preprint2015arXiv

Quenching and Morphological Transformation in Semi-Analytic Models and CANDELS

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

preprint2015arXiv

Sagittarius II, Draco II and Laevens 3: three new Milky Way satellites discovered in the Pan-STARRS 1 3pi Survey

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3pi survey. The first satellite, Laevens 3, is located at a heliocentric distance of d=67+/-3 kpc. With a total magnitude of Mv=-4.4+/-0.3 and a half-light radius rh=7+/-2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4 (Dra II), is a closer and fainter satellite (d~20 kpc, Mv =-2.9+/-0.8), whose uncertain size (rh = 19 +8/-6 pc) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37 +9/-8 pc and Mv=-5.2+/-0.4). At a heliocentric distance of 67+/-5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.

preprint2015arXiv

Stellar Populations of Barred Quiescent Galaxies

Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to-noise (>100) stacked spectra with minimal emission line contamination, we accurately and precisely model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within galaxies. At 0.02<z<0.04, the SDSS fiber radius corresponds to ~1 kpc, which is the typical half-light radii of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the stellar populations at different redshift intervals from z=0.020 to z=0.085 at fixed stellar masses produces an estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming null evolution over this ~1 Gyr epoch. We find that there are no noticeable differences in the slopes of the azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars are not a strong influence on the chemical evolution of quiescent disk galaxies.

preprint2015arXiv

The High-Mass Stellar Initial Mass Function in M31 Clusters

We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope for stars >2 Msun. For the ensemble of clusters, the distribution of stellar MF slopes is best described by $Γ=+1.45^{+0.03}_{-0.06}$ with a very small intrinsic scatter. The data also imply no significant dependencies of the MF slope on cluster age, mass, and size, providing direct observational evidence that the measured MF represents the IMF. This analysis implies that the high-mass IMF slope in M31 clusters is universal with a slope ($Γ=+1.45^{+0.03}_{-0.06}$) that is steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values. Using our inference model on select Milky Way (MW) and LMC high-mass IMF studies from the literature, we find $Γ_{\rm MW} \sim+1.15\pm0.1$ and $Γ_{\rm LMC} \sim+1.3\pm0.1$, both with intrinsic scatter of ~0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be universal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study represents the most robust measurement of the high-mass IMF slope to date. We have grafted the M31 high-mass IMF slope onto widely used sub-solar mass Kroupa and Chabrier IMFs and show that commonly used UV- and Halpha-based star formation rates should be increased by a factor of ~1.3-1.5 and the number of stars with masses >8 Msun are ~25% fewer than expected for a Salpeter/Kroupa IMF. [abridged]

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]

preprint2014arXiv

A new distant Milky Way globular cluster in the Pan-STARRS1 3π survey

We present a new satellite in the outer halo of the Galaxy, the first Milky Way satellite found in the stacked photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1) Survey. From follow-up photometry obtained with WFI on the MPG/ESO 2.2m telescope, we argue that the object, located at a heliocentric distance of 145+/-17 kpc, is the most distant Milky Way globular cluster yet known. With a total magnitude of M_V=-4.3+/-0.2 and a half-light radius of 20+/-2 pc, it shares the properties of extended globular clusters found in the outer halo of our Galaxy and the Andromeda galaxy. The discovery of this distant cluster shows that the full spatial extent of the Milky Way globular cluster system has not yet been fully explored.

preprint2014arXiv

A new population of recently quenched elliptical galaxies in the SDSS

We use the Sloan Digital Sky Survey to investigate the properties of massive elliptical galaxies in the local Universe (z\leq0.08) that have unusually blue optical colors. Through careful inspection, we distinguish elliptical from non-elliptical morphologies among a large sample of similarly blue galaxies with high central light concentrations (c_r\geq2.6). These blue ellipticals comprise 3.7 per cent of all c_r\geq2.6 galaxies with stellar masses between 10^10 and 10^11 h^{-2} {\rm M}_{\sun}. Using published fiber spectra diagnostics, we identify a unique subset of 172 non-star-forming ellipticals with distinctly blue urz colors and young (< 3 Gyr) light-weighted stellar ages. These recently quenched ellipticals (RQEs) have a number density of 2.7-4.7\times 10^{-5}\,h^3\,{\rm Mpc}^{-3} and sufficient numbers above 2.5\times10^{10} h^{-2} {\rm M}_{\sun} to account for more than half of the expected quiescent growth at late cosmic time assuming this phase lasts 0.5 Gyr. RQEs have properties that are consistent with a recent merger origin (i.e., they are strong `first-generation' elliptical candidates), yet few involved a starburst strong enough to produce an E+A signature. The preferred environment of RQEs (90 per cent reside at the centers of < 3\times 10^{12}\,h^{-1}{\rm M}_{\sun} groups) agrees well with the `small group scale' predicted for maximally efficient spiral merging onto their halo center and rules out satellite-specific quenching processes. The high incidence of Seyfert and LINER activity in RQEs and their plausible descendents may heat the atmospheres of small host halos sufficiently to maintain quenching.

preprint2014arXiv

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

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

preprint2014arXiv

Bulgeless Galaxies at Intermediate Redshift: Sample Selection, Colour Properties, and the Existence of Powerful AGN

We present a catalogue of bulgeless galaxies, which includes 19225 objects selected in four of the deepest, largest multi-wavelength datasets available -- COSMOS, AEGIS, GEMS and GOODS -- at intermediate redshift ($0.4 \leq z \leq 1.0$). The morphological classification was provided by the Advanced Camera for Surveys General Catalogue (ACS-GC), which used publicly available data obtained with the ACS instrument on the Hubble Space Telescope. Rest-frame photometric quantities were derived using kcorrect. We analyse the properties of the sample and the evolution of pure-disc systems with redshift. Very massive [$\log (M_\star/M_{\odot}) > 10.5$] bulgeless galaxies contribute to ~30% of the total galaxy population number density at $z \geq 0.7$, but their number density drops substantially with decreasing redshift. We show that only a negligible fraction of pure discs appear to be quiescent systems, and red sequence bulgeless galaxies show indications of dust-obscured star formation. X-ray catalogues were used to search for X-ray emission within our sample. After visual inspection and detailed parametric morphological fitting we identify 30 AGN that reside in galaxies without a classical bulge. The finding of such peculiar objects at intermediate redshift shows that while AGN growth in merger-free systems is a rare event (0.2% AGN hosts in this sample of bulgeless galaxies), it can indeed happen relatively early in the Universe history.

preprint2014arXiv

Galactic Globular and Open Cluster Fiducial Sequences in the Pan-STARRS1 Photometric System

We present the fiducial sequences of a sample of Galactic star clusters in the five bands of the Pan-STARRS1 (PS1) photometric system (g_P1, r_P1, i_P1, z_P1, y_P1). These empirical sequences -- which include the red giant and sub-giant branches, the main sequence, and the horizontal branch -- were defined from deep colour-magnitude diagrams reaching below the oldest main-sequence turn-offs of 13 globular and 3 old open clusters covering a wide range of metallicities (-2.4 < [Fe/H] < +0.4). We find excellent agreement for the nine clusters in common with previous studies in similar photometric systems when transformed to the PS1 system. Because the photometric and spectroscopic properties of these stellar populations are accurately known, the fiducials provide a solid basis for the interpretation of observations in the PS1 system, as well as valuable constraints to improve the empirical colour--$T_{eff}$ relations.

preprint2014arXiv

Galaxy Zoo: Are Bars Responsible for the Feeding of Active Galactic Nuclei at 0.2 < z < 1.0?

We present a new study investigating whether active galactic nuclei (AGN) beyond the local universe are preferentially fed via large-scale bars. Our investigation combines data from Chandra and Galaxy Zoo: Hubble (GZH) in the AEGIS, COSMOS, and GOODS-S surveys to create samples of face-on, disc galaxies at 0.2 < z < 1.0. We use a novel method to robustly compare a sample of 120 AGN host galaxies, defined to have 10^42 erg/s < L_X < 10^44 erg/s, with inactive control galaxies matched in stellar mass, rest-frame colour, size, Sersic index, and redshift. Using the GZH bar classifications of each sample, we demonstrate that AGN hosts show no statistically significant enhancement in bar fraction or average bar likelihood compared to closely-matched inactive galaxies. In detail, we find that the AGN bar fraction cannot be enhanced above the control bar fraction by more than a factor of two, at 99.7% confidence. We similarly find no significant difference in the AGN fraction among barred and non-barred galaxies. Thus we find no compelling evidence that large-scale bars directly fuel AGN at 0.2<z<1.0. This result, coupled with previous results at z=0, implies that moderate-luminosity AGN have not been preferentially fed by large-scale bars since z=1. Furthermore, given the low bar fractions at z>1, our findings suggest that large-scale bars have likely never directly been a dominant fueling mechanism for supermassive black hole growth.

preprint2014arXiv

Galaxy Zoo: CANDELS Barred Disks and Bar Fractions

The formation of bars in disk galaxies is a tracer of the dynamical maturity of the population. Previous studies have found that the incidence of bars in disks decreases from the local Universe to z ~ 1, and by z > 1 simulations predict that bar features in dynamically mature disks should be extremely rare. Here we report the discovery of strong barred structures in massive disk galaxies at z ~ 1.5 in deep rest-frame optical images from CANDELS. From within a sample of 876 disk galaxies identified by visual classification in Galaxy Zoo, we identify 123 barred galaxies. Selecting a sub-sample within the same region of the evolving galaxy luminosity function (brighter than L*), we find that the bar fraction across the redshift range 0.5< z < 2 (f_bar = 10.7 +6.3 -3.5% after correcting for incompleteness) does not significantly evolve. We discuss the implications of this discovery in the context of existing simulations and our current understanding of the way disk galaxies have evolved over the last 11 billion years.

preprint2014arXiv

No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~2 in CANDELS / 3D-HST

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

preprint2014arXiv

Serendipitous Discovery of a Thin Stellar Stream near the Galactic Bulge in the Pan-STARRS1 3Pi Survey

We report the discovery of a thin stellar stream found in Pan-STARRS1 photometry near the Galactic bulge in the constellation of Ophiuchus. It appears as a coherent structure in the colour-selected stellar density maps produced to search for tidal debris around nearby globular clusters. The stream is exceptionally short and narrow; it is about 2.5° long and 6' wide in projection. The colour-magnitude diagram of this object, which harbours a blue horizontal-branch, is consistent with an old and relatively metal-poor population ([Fe/H]~-1.3) located 9.5 +/- 0.9 kpc away at (l,b) ~ (5°,+32°), and 5.0 +/- 1.0 kpc from the Galactic centre. These properties argue for a globular cluster as progenitor. The finding of such a prominent, nearby stream suggests that many streams could await discovery in the more densely populated regions of our Galaxy.

preprint2014arXiv

Spectroscopy of the three distant Andromedan satellites Cassiopeia III, Lacerta I, and Perseus I

We present Keck II/DEIMOS spectroscopy of the three distant dwarf galaxies of M31 Lacerta I, Cassiopeia III, and Perseus I, recently discovered within the Pan-STARRS1 3πimaging survey. The systemic velocities of the three systems (v_{r,helio} = -198.4 +/- 1.1 km/s, -371.6 +/- 0.7 km/s, and -326 +/- 3 km/s, respectively) confirm that they are satellites of M31. In the case of Lacerta I and Cassiopeia III, the high quality of the data obtained for 126 and 212 member stars, respectively, yields reliable constraints on their global velocity dispersions (σ_{vr} = 10.3 +/- 0.9 km/s and 8.4 +/- 0.6 km/s, respectively), leading to dynamical-mass estimates for both of ~4x10^7 Msun within their half-light radius. These translate to V-band mass-to-light ratios of 15^{+12}_{-9} and 8^{+9}_{-5} in solar units. We also use our spectroscopic data to determine the average metallicity of the 3 dwarf galaxies ([Fe/H] = -2.0 +/- 0.1, -1.7 +/- 0.1, and -2.0 +/- 0.2, respectively). All these properties are typical of dwarf galaxy satellites of Andromeda with their luminosity and size.

preprint2014arXiv

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

preprint2014arXiv

The Complex Structure of Stars in the Outer Galactic Disk as revealed by Pan-STARRS1

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement.

preprint2014arXiv

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

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

preprint2014arXiv

The Mass Dependence of Dwarf Satellite Galaxy Quenching

We combine observations of the Local Group with data from the NASA-Sloan Atlas to show the variation in the quenched fraction of satellite galaxies from low mass dwarf spheroidals and dwarf irregulars to more massive dwarfs similar to the Magellanic clouds. While almost all of the low mass ($M_\star \lesssim 10^7$ $M_\odot$) dwarfs are quenched, at higher masses the quenched fraction decreases to approximately 40-50%. This change in the quenched fraction is large, and suggests a sudden change in the effectiveness of quenching that correlates with satellite mass. We combine this observation with models of satellite infall and ram pressure stripping to show that the low mass satellites must quench within 1-2 Gyr of pericenter passage to maintain a high quenched fraction, but that many more massive dwarfs must continue to form stars today even though they likely fell in to their host >5 Gyr ago. We also characterize how the susceptibility of dwarfs to ram pressure must vary as a function of mass if it is to account for the change in quenched fractions. Though neither model predicts the quenching effectiveness a priori, this modeling illustrates the physical requirements that the observed quenched fractions place on possible quenching mechanisms.

preprint2014arXiv

The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy

We present near-infrared spectroscopy of a sample of 22 Extreme Emission Line Galaxies at redshifts 1.3 < z < 2.3, confirming that these are low-mass (M* = 10^8 - 10^9 M_sun) galaxies undergoing intense starburst episodes (M*/SFR ~ 10-100 Myr). The sample is selected by [O III] or Hα emission line flux and equivalent width using near-infrared grism spectroscopy from the 3D-HST survey. High-resolution NIR spectroscopy is obtained with LBT/LUCI and VLT/X-SHOOTER. The [O III]/Hβ line ratio is high (> 5) and [N II]/Hα is always significantly below unity, which suggests a low gas-phase metallicity. We are able to determine gas-phase metallicities for 7 of our objects using various strong-line methods, with values in the range 0.05-0.30 Z_sun and with a median of 0.15 Z_sun; for 3 of these objects we detect [O III]λ4363 which allows for a direct constraint on the metallicity. The velocity dispersion, as measured from the nebular emission lines, is typically ~50 km/s. Combined with the observed star-forming activity, the Jeans and Toomre stability criteria imply that the gas fraction must be large (> 2/3), consistent with the difference between our dynamical and stellar mass estimates. The implied gas depletion time scale (several hundred Myr) is substantially longer than the inferred mass-weighted ages (~50 Myr), which further supports the emerging picture that most stars in low-mass galaxies form in short, intense bursts of star formation.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury V: Ages and Masses of the Year 1 Stellar Clusters

We present ages and masses for 601 star clusters in M31 from the analysis of the six filter integrated light measurements from near ultraviolet to near infrared wavelengths, made as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We derive the ages and masses using a probabilistic technique, which accounts for the effects of stochastic sampling of the stellar initial mass function. Tests on synthetic data show that this method, in conjunction with the exquisite sensitivity of the PHAT observations and their broad wavelength baseline, provides robust age and mass recovery for clusters ranging from $\sim 10^2 - 2 \times 10^6 M_\odot$. We find that the cluster age distribution is consistent with being uniform over the past $100$ Myr, which suggests a weak effect of cluster disruption within M31. The age distribution of older ($>100$ Myr) clusters fall towards old ages, consistent with a power-law decline of index $-1$, likely from a combination of fading and disruption of the clusters. We find that the mass distribution of the whole sample can be well-described by a single power-law with a spectral index of $-1.9 \pm 0.1$ over the range of $10^3-3 \times 10^5 M_\odot$. However, if we subdivide the sample by galactocentric radius, we find that the age distributions remain unchanged. However, the mass spectral index varies significantly, showing best fit values between $-2.2$ and $-1.8$, with the shallower slope in the highest star formation intensity regions. We explore the robustness of our study to potential systematics and conclude that the cluster mass function may vary with respect to environment.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury X. Ultraviolet to Infrared Photometry of 117 Million Equidistant Stars

We have measured stellar photometry with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) in near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) bands for 117 million resolved stars in M31. As part of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, we measured photometry with simultaneous point spread function fitting across all bands and at all source positions after precise astrometric image alignment (<5-10 milliarcsecond accuracy). In the outer disk, the photometry reaches a completeness-limited depth of F475W~28, while in the crowded, high surface brightness bulge, the photometry reaches F475W~25. We find that simultaneous photometry and optimized measurement parameters significantly increase the detection limit of the lowest resolution filters (WFC3/IR) providing color-magnitude diagrams that are up to 2.5 magnitudes deeper when compared with color-magnitude diagrams from WFC3/IR photometry alone. We present extensive analysis of the data quality including comparisons of luminosity functions and repeat measurements, and we use artificial star tests to quantify photometric completeness, uncertainties and biases. We find that largest sources of systematic error in the photometry are due to spatial variations in the point spread function models and charge transfer efficiency corrections. This stellar catalog is the largest ever produced for equidistant sources, and is publicly available for download by the community.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury. VI. The reliability of far-ultraviolet flux as a star formation tracer on sub-kpc scales

We have used optical observations of resolved stars from the Panchromatic Hubble Andromeda Treasury (PHAT) to measure the recent (< 500 Myr) star formation histories (SFHs) of 33 FUV-bright regions in M31. The region areas ranged from ~$10^4$ to $10^6$ pc$^2$, which allowed us to test the reliability of FUV flux as a tracer of recent star formation on sub-kpc scales. The star formation rates (SFRs) derived from the extinction-corrected observed FUV fluxes were, on average, consistent with the 100-Myr mean SFRs of the SFHs to within the 1$σ$ scatter. Overall, the scatter was larger than the uncertainties in the SFRs and particularly evident among the smallest regions. The scatter was consistent with an even combination of discrete sampling of the initial mass function and high variability in the SFHs. This result demonstrates the importance of satisfying both the full-IMF and the constant-SFR assumptions for obtaining precise SFR estimates from FUV flux. Assuming a robust FUV extinction correction, we estimate that a factor of 2.5 uncertainty can be expected in FUV-based SFRs for regions smaller than $10^5$ pc$^2$, or a few hundred pc. We also examined ages and masses derived from UV flux under the common assumption that the regions are simple stellar populations (SSPs). The SFHs showed that most of the regions are not SSPs, and the age and mass estimates were correspondingly discrepant from the SFHs. For those regions with SSP-like SFHs, we found mean discrepancies of 10 Myr in age and a factor of 3 to 4 in mass. It was not possible to distinguish the SSP-like regions from the others based on integrated FUV flux.

preprint2013arXiv

A Critical Assessment of Photometric Redshift Methods: A CANDELS Investigation

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

preprint2013arXiv

A New Approach to Detailed Structural Decomposition from the SPLASH and PHAT Surveys: Kicked-up Disk Stars in the Andromeda Galaxy?

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 4 < R_proj < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ~20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < m_814W < 22) stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sersic, exponential, and cored power-law, respectively, and the LFs are modeled as broken power-laws. We present probability distributions for each of 32 parameters describing the SB profiles and LFs of the three subcomponents. We find that the number of stars with a disk-like LF is ~5% larger than the the number with disk-like (dynamically cold) kinematics, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a "kicked-up disk" halo population in M31. The fraction of kicked-up disk stars is consistent with that found in simulations. We also find evidence for a radially varying disk LF, consistent with a negative metallicity gradient in the stellar disk.

preprint2013arXiv

A Tidally-Stripped Stellar Component of the Magellanic Bridge

Deep photometry of the Small Magellanic Cloud (SMC) stellar periphery (R=4 deg, 4.2 kpc) is used to study its line-of-sight depth with red clump (RC) stars. The RC luminosity function is affected little by young (<1 Gyr) blue-loop stars in these regions because their main-sequence counterparts are not observed in the color magnitude diagrams. The SMC's eastern side is found to have a large line-of-sight depth (~23 kpc) while the western side has a much shallower depth (~10 kpc), consistent with previous photographic plate photometry results. We use a model SMC RC luminosity function to deconvolve the observed RC magnitudes and construct the density function in distance for our fields. Three of the eastern fields show a distance bimodality with one component at the "systemic" ~67 kpc SMC distance and a second component at ~55 kpc. Our data are not reproduced well by the various extant Magellanic Cloud and Stream simulations. However, the models predict that the known HI Magellanic Bridge (stretching from the SMC eastward towards the LMC) has a decreasing distance with angle from the SMC and should be seen in both the gaseous and stellar components. From comparison with these models we conclude that the most likely explanation for our newly identified ~55 kpc stellar structure in the eastern SMC is a stellar counterpart of the HI Magellanic Bridge that was tidally stripped from the SMC ~200 Myr ago during a close encounter with the LMC. This discovery has important implications for microlensing surveys of the SMC.

preprint2013arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

Confirmation of Small Dynamical and Stellar Masses for Extreme Emission Line Galaxies at z~2

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal kinematically narrow lines (~50 km/s) for a sample of 14 Extreme Emission Line Galaxies (EELGs) at redshifts 1.4 < z < 2.3. These measurements imply that the total dynamical masses of these systems are low (< 3x10^9 M_sun). Their large [O III] 5007 equivalent widths (500-1100 Angstroms) and faint blue continuum emission imply young ages of 10-100 Myr and stellar masses of 10^8-10^9 M_sun, confirming the presence of a violent starburst. The dynamical masses represent the first such determinations for low-mass galaxies at z > 1. The stellar mass formed in this vigorous starburst phase represents a large fraction of the total (dynamical) mass, without a significantly massive underlying population of older stars. The occurrence of such intense events in shallow potentials strongly suggests that supernova-driven winds must be of critical importance in the subsequent evolution of these systems.

preprint2013arXiv

Confronting Models of Dwarf Galaxy Quenching with Observations of the Local Group

A number of mechanisms have been proposed to connect star-forming dwarf irregular galaxies with the formation of non-star-forming dwarf spheroidal galaxies, but distinguishing between these mechanisms has been difficult. We use the Via Lactea dark matter only cosmological simulations to test two well-motivated simple hypotheses---transformation of irregulars into dwarf spheroidal galaxies by tidal stirring and ram pressure stripping following a close passage to the host galaxy, and transformation via mergers between dwarfs---and predict the radial distribution and inferred formation times of the resulting dwarf spheroidal galaxies. We compare this to the observed distribution in the Local Group and show that 1) the observed dSph distribution far from the Galaxy or M31 can be matched by the VL halos that have passed near the host galaxy at least once, though significant halo-to-halo scatter exists, 2) models that require two or more pericenter passages for dSph-formation cannot account for the dSphs beyond 500 kpc such as Cetus and Tucana, and 3) mergers predict a flat radial distribution of dSphs and cannot account for the high dSph fraction near the Galaxy, but are not ruled out at large distances. The models also suggest that for dSphs found today beyond 500 kpc, mergers tend to occur significantly earlier than dwarf--host encounters, thus leading to a potentially observable difference in stellar populations. We argue that tidal interactions are sufficient to reproduce the observed distribution of dSphs if and only if a single pericenter passage is sufficient to form a dSph.

preprint2013arXiv

Constraining the age of the NGC 4565 HI Disk Warp: Determining the Origin of Gas Warps

We have mapped the distribution of young and old stars in the gaseous HI warp of NGC 4565. We find a clear correlation of young stars (<600 Myr) with the warp, but no coincident old stars (>1 Gyr), which places an upper limit on the age of the structure. The formation rate of the young stars, which increased ~300 Myr ago relative to the surrounding regions, is (6.3 +2.5/-1.5) x 10^-5 M_sol/yr/kpc^2. This implies a ~60+/-20 Gyr depletion time of the HI warp, similar to the timescales calculated for the outer HI disks of nearby spiral galaxies. While some stars associated with the warp fall into the asymptotic giant branch (AGB) region of the color magnitude diagram, where stars could be as old as 1 Gyr, further investigation suggests that they may be interlopers rather than real AGB stars. We discuss the implications of these age constraints for the formation of HI warps, and the gas fueling of disk galaxies.

preprint2013arXiv

Evidence for an Interaction in the Nearest Starbursting Dwarf Irregular Galaxy IC 10

Using deep 21-cm HI data from the Green Bank Telescope we have detected an ~18.3 kpc-long gaseous extension associated with the starbursting dwarf galaxy IC 10. The newly-found feature stretches 1.3 deg to the northwest and has a large radial velocity gradient reaching to ~65 km/s lower than the IC 10 systemic velocity. A region of higher column density at the end of the extension that possesses a coherent velocity gradient (~10 km/s across ~26 arcmin) transverse to the extension suggests rotation and may be a satellite galaxy of IC 10. The HI mass of IC 10 is 9.5x10^7 (d/805 kpc)^2 Msun and the mass of the new extension is 7.1x10^5 (d/805 kpc)^2 Msun. An IC 10-M31 orbit using known radial velocity and proper motion values for IC 10 show that the HI extension is inconsistent with the trailing portion of the orbit so that an M31-tidal or ram pressure origin seems unlikely. We argue that the most plausible explanation for the new feature is that it is the result of a recent interaction (and possible late merger) with another dwarf galaxy. This interaction could not only have triggered the origin of the recent starburst in IC 10, but could also explain the existence of previously-found counter-rotating HI gas in the periphery of the IC 10 which was interpreted as originating from primordial gas infall.

preprint2013arXiv

Galaxy Zoo: Observing Secular Evolution Through Bars

In this paper, we use the Galaxy Zoo 2 dataset to study the behavior of bars in disk galaxies as a function of specific star formation rate (SSFR), and bulge prominence. Our sample consists of 13,295 disk galaxies, with an overall (strong) bar fraction of $23.6\pm 0.4\%$, of which 1,154 barred galaxies also have bar length measurements. These samples are the largest ever used to study the role of bars in galaxy evolution. We find that the likelihood of a galaxy hosting a bar is anti-correlated with SSFR, regardless of stellar mass or bulge prominence. We find that the trends of bar likelihood and bar length with bulge prominence are bimodal with SSFR. We interpret these observations using state-of-the-art simulations of bar evolution which include live halos and the effects of gas and star formation. We suggest our observed trends of bar likelihood with SSFR are driven by the gas fraction of the disks; a factor demonstrated to significantly retard both bar formation and evolution in models. We interpret the bimodal relationship between bulge prominence and bar properties as due to the complicated effects of classical bulges and central mass concentrations on bar evolution, and also to the growth of disky pseudobulges by bar evolution. These results represent empirical evidence for secular evolution driven by bars in disk galaxies. This work suggests that bars are not stagnant structures within disk galaxies, but are a critical evolutionary driver of their host galaxies in the local universe ($z<1$).

preprint2013arXiv

GHOSTS I: A New Faint very Isolated Dwarf Galaxy at D = 12 +/- 2 Mpc

We report the discovery of a new faint dwarf galaxy, GHOSTS I, using HST/ACS data from one of our GHOSTS (Galaxy Halos, Outer disks, Substructure, Thick disk, and Star clusters) fields. Its detected individual stars populate an approximately one magnitude range of its luminosity function (LF). Using synthetic color-magnitude diagrams (CMDs) to compare with the galaxy's CMD, we find that the colors and magnitudes of GHOSTS I's individual stars are most consistent with being young helium-burning and asymptotic giant branch stars at a distance of 12 +/- 2 Mpc. Morphologically, GHOSTS I appears to be actively forming stars, so we tentatively classify it as a dwarf irregular (dIrr) galaxy, although future HST observations deep enough to resolve a larger magnitude range in its LF are required to make a more secure classification. GHOSTS I's absolute magnitude is $M_V = -9.85^{+ 0.40}_{- 0.33}$, making it one of the least luminous dIrr galaxies known, and its metallicity is lower than [Fe/H] =-1.5 dex. The half-light radius of GHOSTS I is 226 +/- 38 pc and its ellipticity is 0.47 +/- 0.07, similar to Milky Way and M31 dwarf satellites at comparable luminosity. There are no luminous massive galaxies or galaxy clusters within ~ 4 Mpc from GHOSTS I that could be considered as its host, making it a very isolated dwarf galaxy in the Local Universe.

preprint2013arXiv

Lacerta I and Cassiopeia III: Two luminous and distant Andromeda satellite dwarf galaxies found in the 3π Pan-STARRS1 survey

We report the discovery of two new dwarf galaxies, Lacerta I/Andromeda XXXI (Lac I/And XXXI) and Cassiopeia III/Andromeda XXXII (Cas III/And XXXII), in stacked Pan-STARRS1 r_P1- and i_P1-band imaging data. Both are luminous systems (M_V ~ -12) located at projected distances of 20.3° and 10.5° from M31. Lac I and Cas III are likely satellites of the Andromeda galaxy with heliocentric distances of 756^{+44}_{-28} kpc and 772^{+61}_{-56} kpc, respectively, and corresponding M31-centric distances of 275+/-7 kpc and 144^{+6}_{-4} kpc . The brightest of recent Local Group member discoveries, these two new dwarf galaxies owe their late discovery to their large sizes (r_h = 4.2^{+0.4}_{-0.5} arcmin or 912^{+124}_{-93} pc for Lac I; r_h = 6.5^{+1.2}_{-1.0} arcmin or 1456+/-267 pc for Cas III), and consequently low surface brightness (μ_0 ~ 26.0 mag/arcsec^2), as well as to the lack of a systematic survey of regions at large radii from M31, close to the Galactic plane. This latter limitation is now alleviated by the 3π Pan-STARRS1 survey, which could lead to the discovery of other distant Andromeda satellite dwarf galaxies.

preprint2013arXiv

Perseus I: A distant satellite dwarf galaxy of Andromeda

We present the discovery of a new dwarf galaxy, Perseus I/Andromeda XXXIII, found in the vicinity of Andromeda (M31) in stacked imaging data from the Pan-STARRS1 3π survey. Located 27.9° away from M31, Perseus I has a heliocentric distance of 785 +/- 65 kpc, compatible with it being a satellite of M31 at 374 +14/-10 kpc from its host. The properties of Perseus I are typical for a reasonably bright dwarf galaxy (M_V = -10.3 +/- 0.7), with an exponential half-light radius of r_h = 1.7 +/- 0.4 arcminutes or r_h = 400 +105/-85 pc at this distance, and a moderate ellipticity (ε= 0.43 +0.15/-0.17). The late discovery of Perseus I is due to its fairly low surface brightness (μ_0=25.7 +1.0/-0.9 mag/arcsec^2), and to the previous lack of deep, high quality photometric data in this region. If confirmed to be a companion of M31, the location of Perseus I, far east from its host, could place interesting constraints on the bulk motion of the satellite system of M31.

preprint2013arXiv

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

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

preprint2013arXiv

Testing galaxy formation models with the GHOSTS survey: The color profile of M81's stellar halo

We study the properties of the stellar populations in M81's outermost part, which hereafter we will term the stellar halo, using HST ACS/WFC observations of 19 fields from the GHOSTS survey. The observed fields probe the stellar halo out to a projected distance of ~ 50 kpc from the galactic center. Each field was observed in both F606W and F814W filters. The 50% completeness levels of the color magnitude diagrams (CMDs) are typically at 2 mag below the tip of the red giant branch (TRGB). Fields at distances closer than 15 kpc show evidence of disk-dominated populations whereas fields at larger distances are mostly populated by halo stars. The RGB of the M81's halo CMDs is well matched with isochrones of ~ 10 Gyr and metallicities [Fe/H] ~ -1.2 dex, suggesting that the dominant stellar population of M81's halo has a similar age and metallicity. The halo of M81 is characterized by a color distribution of width ~ 0.4 mag and an approximately constant median value of (F606W - F814W) ~ 1 mag measured using stars within the magnitude range 23.7 < F814W < 25.5. When considering only fields located at galactocentric radius R > 15 kpc, we detect no color gradient in the stellar halo of M81. We place a limit of 0.03+/-0.11 mag difference between the median color of RGB M81 halo stars at ~ 15 and at 50 kpc, corresponding to a metallicity difference of 0.08+/-0.35 dex over that radial range for an assumed constant age of 10 Gyr. We compare these results with model predictions for the colors of stellar halos formed purely via accretion of satellite galaxies. When we analyze the cosmologically motivated models in the same way as the HST data, we find that they predict no color gradient for the stellar halos, in good agreement with the observations.

preprint2013arXiv

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

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

preprint2012arXiv

A CANDELS WFC3 Grism Study of Emission-Line Galaxies at z~2: A Mix of Nuclear Activity and Low-Metallicity Star Formation

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

preprint2012arXiv

AGN host galaxies at redshift z~0.7: peculiar or not?

We perform a quantitative morphological comparison between the hosts of Active Galactic Nuclei (AGN) and quiescent galaxies at intermediate redshifts (z~0.7). The imaging data are taken from the large HST/ACS mosaics of the GEMS and STAGES surveys. Our main aim is to test whether nuclear activity at this cosmic epoch is triggered by major mergers. Using images of quiescent galaxies and stars, we create synthetic AGN images to investigate the impact of an optical nucleus on the morphological analysis of AGN hosts. Galaxy morphologies are parameterized using the asymmetry index A, concentration index C, Gini coefficient G and M20 index. A sample of ~200 synthetic AGN is matched to 21 real AGN in terms of redshift, host brightness and host-to-nucleus ratio to ensure a reliable comparison between active and quiescent galaxies. The optical nuclei strongly affect the morphological parameters of the underlying host galaxy. Taking these effects into account, we find that the morphologies of the AGN hosts are clearly distinct from galaxies undergoing violent gravitational interactions. In fact, the host galaxies' distributions in morphological descriptor space are more similar to undisturbed galaxies than major mergers. Intermediate-luminosity (Lx < 10^44 erg/s) AGN hosts at z~0.7 show morphologies similar to the general population of massive galaxies with significant bulges at the same redshifts. If major mergers are the driver of nuclear activity at this epoch, the signatures of gravitational interactions fade rapidly before the optical AGN phase starts, making them undetectable on single-orbit HST images, at least with usual morphological descriptors. This could be investigated in future synthetic observations created from numerical simulations of galaxy-galaxy interactions.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Modeling the red sequence: Hierarchical growth yet slow luminosity evolution

We explore the effects of mergers on the evolution of massive early-type galaxies by modeling the evolution of their stellar populations in a hierarchical context. We investigate how a realistic red sequence population set up by z~1 evolves under different assumptions for the merger and star formation histories, comparing changes in color, luminosity and mass. The purely passive fading of existing red sequence galaxies, with no further mergers or star formation, results in dramatic changes at the bright end of the luminosity function and color-magnitude relation. Without mergers there is too much evolution in luminosity at a fixed space density compared to observations. The change in color and magnitude at a fixed mass resemble that of a passively evolving population that formed relatively recently, at z~2. Mergers amongst the red sequence population ("dry mergers") occurring after z=1 build up mass, counteracting the fading of the existing stellar populations to give smaller changes in both color and luminosity for massive galaxies. By allowing some galaxies to migrate from the blue cloud onto the red sequence after z=1 through gas-rich mergers, younger stellar populations are added to the red sequence. This manifestation of the progenitor bias increases the scatter in age and results in even smaller changes in color and luminosity between z=1 and z=0 at a fixed mass. The resultant evolution appears much slower, resembling the passive evolution of a population that formed at high redshift (z~3-5) and is in closer agreement with observations. Measurements of the luminosity and color evolution alone are not sufficient to distinguish between the purely passive evolution of an old population and cosmologically motivated hierarchical growth, although these scenarios have very different implications for the mass growth of early-type galaxies over the last half of cosmic history.

preprint2012arXiv

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

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

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 Nature of Extremely Red H-[4.5]>4 Galaxies revealed with SEDS and CANDELS

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury I: Bright UV Stars in the Bulge of M31

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqué stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqué (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or α abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury II. Tracing the Inner M31 Halo with Blue Horizontal Branch Stars

We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galactocentric distances. If we assume a power-law profile r^(-α) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of α=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to α<1.2 inside of 3 kpc. This slope is consistent with previous measurements but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass of 2.1^{+1.7}_{-0.4} x 10^9 M_sun. These properties are comparable with both simulations of stellar halo formation formed by satellite disruption alone, and with simulations that include some in situ formation of halo stars.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury III. Measuring Ages and Masses of Partially Resolved Stellar Clusters

The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the stellar initial mass function, when inferred from the integrated color of low mass clusters (less than ~10^4 solar masses). We use simulated star clusters to show that these effects are minimized when the brightest, rapidly evolving stars in a cluster can be resolved, and the light of the fainter, more numerous unresolved stars can be analyzed separately. When comparing the light from the less luminous cluster members to models of unresolved light, more accurate age estimates can be obtained than when analyzing the integrated light from the entire cluster under the assumption that the initial mass function is fully populated. We show the success of this technique first using simulated clusters, and then with a stellar cluster in M31. This method represents one way of accounting for the discrete, stochastic sampling of the stellar initial mass function in less massive clusters and can be leveraged in studies of clusters throughout the Local Group and other nearby galaxies.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury IV. A Probabilistic Approach to Inferring the High Mass Stellar Initial Mass Function and Other Power-law Functions

We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, α, are unbiased and that the uncertainty, Δα, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on α and provide an analytic approximation for Δα as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <α>=2.46 with a 1-σ dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged)

preprint2012arXiv

The Star Formation History of Leo T from Hubble Space Telescope Imaging

We present the star formation history (SFH) of the faintest known star-forming galaxy, Leo T, based on imaging taken with the Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC2). The HST/WFPC2 color-magnitude diagram (CMD) of Leo T is exquisitely deep, extending ~ 2 magnitudes below the oldest main sequence turnoff, permitting excellent constraints on star formation at all ages. We use a maximum likelihood CMD fitting technique to measure the SFH of Leo T assuming three different sets of stellar evolution models: Padova (solar-scaled metallicity) and BaSTI (both solar-scaled and alpha-enhanced metallicities). The resulting SFHs are remarkably consistent at all ages, indicating that our derived SFH is robust to the choice of stellar evolution model. From the lifetime SFH of Leo T, we find that 50% of the total stellar mass formed prior to z ~ 1 (7.6 Gyr ago). Subsequent to this epoch, the SFH of Leo T is roughly constant until the most recent ~ 25 Myr, where the SFH shows an abrupt drop. This decrease could be due to a cessation of star formation or stellar initial mass function sampling effects, but we are unable to distinguish between the two scenarios. Overall, our measured SFH is consistent with previously derived SFHs of Leo T. However, the HST-based solution provides improved age resolution and reduced uncertainties at all epochs. The SFH, baryonic gas fraction, and location of Leo T are unlike any of the other recently discovered faint dwarf galaxies in the Local Group, and instead bear strong resemblance to gas-rich dwarf galaxies (irregular or transition), suggesting that gas-rich dwarf galaxies may share common modes of star formation over a large range of stellar mass (~ 10^5-10^9 Msun).

preprint2012arXiv

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

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

preprint2011arXiv

A New Automatic Method to Identify Galaxy Mergers I. Description and Application to the STAGES Survey

We present an automatic method to identify galaxy mergers using the morphological information contained in the residual images of galaxies after the subtraction of a Sersic model. The removal of the bulk signal from the host galaxy light is done with the aim of detecting the fainter minor mergers. The specific morphological parameters that are used in the merger diagnostic suggested here are the Residual Flux Fraction and the asymmetry of the residuals. The new diagnostic has been calibrated and optimized so that the resulting merger sample is very complete. However, the contamination by non-mergers is also high. If the same optimization method is adopted for combinations of other structural parameters such as the CAS system, the merger indicator we introduce yields merger samples of equal or higher statistical quality than the samples obtained through the use of other structural parameters. We explore the ability of the method presented here to select minor mergers by identifying a sample of visually classified mergers that would not have been picked up by the use of the CAS system, when using its usual limits. Given the low prevalence of mergers among the general population of galaxies and the optimization used here, we find that the merger diagnostic introduced in this work is best used as a negative merger test, i.e., it is very effective at selecting non-merging galaxies. As with all the currently available automatic methods, the sample of merger candidates selected is contaminated by non-mergers, and further steps are needed to produce a clean sample. This merger diagnostic has been developed using the HST/ACS F606W images of the A901/02 cluster (z=0.165) obtained by the STAGES team. In particular, we have focused on a mass and magnitude limited sample (log M/M_{O}>9.0, R_{Vega}<23.5mag)) which includes 905 cluster galaxies and 655 field galaxies of all morphological types.

preprint2011arXiv

Andromeda XXIX: A New Dwarf Spheroidal Galaxy 200kpc from Andromeda

We report the discovery of a new dwarf galaxy, Andromeda XXIX, using data from the recently-released Sloan Digital Sky Survey DR8, and confirmed by Gemini North telescope Multi-Object Spectrograph imaging data. And XXIX appears to be a dwarf spheroidal galaxy, separated on the sky by a little more than 15 degrees from M31, with a distance inferred from the tip of the red giant branch of 730kpc+/-75kpc, corresponding to a three dimensional separation from M31 of between 205kpc and 227kpc (close to M31's virial radius). Its absolute magnitude, as determined by comparison to the red giant branch luminosity function of the Draco dwarf spheroidal, is M_V = -8.3+/-0.4. And XXIX's stellar populations appear very similar to Draco's; consequently, we estimate a metallicity for And XXIX of [Fe/H]\sim-1.8. The half-light radius of And XXIX is 360pc+/-60pc and its ellipticity is 0.35+/-0.06, typical of dwarf satellites of the Milky Way and M31 at this absolute magnitude range.

preprint2011arXiv

Andromeda XXVIII: A Dwarf Galaxy More Than 350 kpc from Andromeda

We report the discovery of a new dwarf galaxy, Andromeda XXVIII, using data from the recently-released SDSS DR8. The galaxy is a likely satellite of Andromeda, and, at a separation of $365^{+17}_{-1}$ kpc, would be one of the most distant of Andromeda's satellites. Its heliocentric distance is $650^{+150}_{-80}$ kpc, and analysis of its structure and luminosity show that it has an absolute magnitude of $M_V = -8.5^{+0.4}_{-1.0}$ and half-light radius of $r_h = 210^{+60}_{-50}$ pc, similar to many other faint Local Group dwarfs. With presently-available imaging we are unable to determine if there is ongoing or recent star formation, which prevents us from classifying it as a dwarf spheroidal or dwarf irregular.

preprint2011arXiv

CANDELS: Constraining the AGN-Merger Connection with Host Morphologies at z~2

Using HST/WFC3 imaging taken as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we examine the role that major galaxy mergers play in triggering active galactic nuclei (AGN) activity at z~2. Our sample consists of 72 moderate-luminosity (Lx ~ 1E42-1E44 erg/s) AGN at 1.5<z<2.5 that are selected using the 4 Msec Chandra observations in the Chandra Deep Field South, the deepest X-ray observations to date. Employing visual classifications, we have analyzed the rest-frame optical morphologies of the AGN host galaxies and compared them to a mass-matched control sample of 216 non-active galaxies at the same redshift. We find that most of the AGN reside in disk galaxies (51.4%), while a smaller percentage are found in spheroids (27.8%). Roughly 16.7% of the AGN hosts have highly disturbed morphologies and appear to be involved in a major merger or interaction, while most of the hosts (55.6%) appear relatively relaxed and undisturbed. These fractions are statistically consistent with the fraction of control galaxies that show similar morphological disturbances. These results suggest that the hosts of moderate-luminosity AGN are no more likely to be involved in an ongoing merger or interaction relative to non-active galaxies of similar mass at z~2. The high disk fraction observed among the AGN hosts also appears to be at odds with predictions that merger-driven accretion should be the dominant AGN fueling mode at z~2, even at moderate X-ray luminosities. Although we cannot rule out that minor mergers are responsible for triggering these systems, the presence of a large population of relatively undisturbed disk-like hosts suggests that secular processes play a greater role in fueling AGN activity at z~2 than previously thought.

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

Quantifying Kinematic Substructure in the Milky Way's Stellar Halo

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

preprint2011arXiv

The environmental dependence of the structure of outer galactic discs in STAGES spiral galaxies

We present an analysis of V-band radial surface brightness profiles for spiral galaxies from the field and cluster environments using Hubble Space Telescope/Advanced Camera for Surveys imaging and data from the Space Telescope A901/2 Galaxy Evolution Survey (STAGES). We use a large sample of ~330 face-on to intermediately inclined spiral galaxies and assess the effect of the galaxy environment on the azimuthally averaged radial surface brightness mu profiles for each galaxy in the outer stellar disc (24 < mu < 26.5 mag per sq arcsec). For galaxies with a purely exponential outer disc (~50 per cent), we determine the significance of an environmental dependence on the outer disc scalelength h_out. For galaxies with a broken exponential in their outer disc, either down-bending (truncation, ~10 per cent) or up-bending (anti-truncation, ~40 per cent), we measure the strength T (outer-to-inner scalelength ratio, log_10(h_out/h_in) of the mu breaks and determine the significance of an environmental dependence on break strength T. Surprisingly, we find no evidence to suggest any such environmental dependence on either outer disc scalelength h_out or break strength T, implying that the galaxy environment is not affecting the stellar distribution in the outer stellar disc. We also find that for galaxies with small effective radii (r_e < 3 kpc) there is a lack of outer disc truncations in both the field and cluster environments. Our results suggest that the stellar distribution in the outer disc of spiral galaxies is not significantly affected by the galaxy environment.

preprint2011arXiv

The Majority of Compact Massive Galaxies at z~2 are Disk Dominated

We investigate the stellar structure of massive, quiescent galaxies at z~2, based on Hubble Space Telescope/WFC3 imaging from the Early Release Science program. Our sample of 14 galaxies has stellar masses of M* > 10^{10.8} Msol and photometric redshifts of 1.5 < z < 2.5. In agreement with previous work, their half-light radii are <2 kpc, much smaller than equally massive galaxies in the present-day universe. A significant subset of the sample appears highly flattened in projection, which implies, considering viewing angle statistics, that a significant fraction of the galaxies in our sample have pronounced disks. This is corroborated by two-dimensional surface brightness profile fits. We estimate that 65% +/- 15% of the population of massive, quiescent z~2 galaxies are disk-dominated. The median disk scale length is 1.5 kpc, substantially smaller than the disks of equally massive galaxies in the present-day universe. Our results provide strong observational evidence that the much-discussed ultra-dense high-redshift galaxies should generally be thought of as disk-like stellar systems with the majority of stars formed from gas that had time to settle into a disk.

preprint2011arXiv

The Resolved Stellar Halo of NGC 253

We have obtained Magellan/IMACS and HST/ACS imaging data that resolve red giant branch stars in the stellar halo of the starburst galaxy NGC 253. The HST data cover a small area, and allow us to accurately interpret the ground-based data, which cover 30% of the halo to a distance of 30 kpc, allowing us to make detailed quantitative measurements of the global properties and structure of a stellar halo outside of the Local Group. The geometry of the halo is significantly flattened in the same sense as the disk, with a projected axis ratio of b/a ~ 0.35 +/- 0.1. The total stellar mass of the halo is estimated to be M_halo ~ 2.5 +/- 1.5 x 10^9 M_sun, or 6% of the total stellar mass of the galaxy, and has a projected radial dependence that follows a power law of index -2.8 +/- 0.6, corresponding to a three-dimensional power law index of ~ -4. The total luminosity and profile shape that we measure for NGC 253 are somewhat larger and steeper than the equivalent values for the Milky Way and M31, but are well within the scatter of model predictions for the properties of stellar halos built up in a cosmological context. Structure within the halo is seen at a variety of scales: there is small kpc-scale density variation and a large shelf-like feature near the middle of the field. The techniques that have been developed will be essential for quantitatively comparing our upcoming larger sample of observed stellar halos to models of halo formation.

preprint2011arXiv

The star formation history of mass-selected galaxies in the COSMOS field

We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.2<z<3 both populations show a strong and M_*-independent decrease in their SSFR towards z=0.2, best described by a power- law (1+z)^n, where n~4.3 for all galaxies and n~3.5 for SF sources. The decrease appears to have started at z>2, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples.

preprint2011arXiv

The Structure of the Sagittarius Stellar Stream as Traced by Blue Horizontal Branch Stars

We use a sample of blue horizontal branch (BHB) stars from the Sloan Digital Sky Survey Data Release 7 to explore the structure of the tidal tails from the Sagittarius Dwarf Galaxy. We use a method yielding BHB star candidates with up to ~70% purity from photometry alone. The resulting sample has a distance precision of roughly 5% and can probe distances in excess of 100 kpc. Using this sample, we identify a possible extension to the trailing arm at distances of 60-80 kpc from the Sun with an estimated significance of at least 3.8 sigma. Current models predict that a distant `returning' segment of the debris stream should exist, but place it substantially closer to the Sun where no debris is observed in our data. Exploiting the distance precision of our tracers, we estimate the mean line-of-sight thickness of the leading arm to be ~3 kpc, and show that the two `bifurcated' branches of the debris stream differ by only 1-2 kpc in distance. With a spectroscopic very pure BHB star subsample, we estimate the velocity dispersion in the leading arm, 37 km s^-1, which is in reasonable agreement with models of Sgr disruption. We finally present a sample of high-probability Sgr BHB stars in the leading arm of Sgr, selected to have distances and velocities consistent with Sgr membership, to allow further study.

preprint2010arXiv

Dust-Obscured Star-Formation in Intermediate Redshift Galaxy Clusters

We present Spitzer MIPS 24-micron observations of 16 0.4<z<0.8 galaxy clusters drawn from the ESO Distant Cluster Survey (EDisCS). This is the first large 24-micron survey of clusters at intermediate redshift. The depth of our imaging corresponds to a total IR luminosity of 8x10^10 Lsun, just below the luminosity of luminous infrared galaxies (LIRGs), and 6^{+1}_{-1}% of M_V < -19 cluster members show 24-micron emission at or above this level. We compare with a large sample of coeval field galaxies and find that while the fraction of cluster LIRGs lies significantly below that of the field, the IR luminosities of the field and cluster galaxies are consistent. However, the stellar masses of the EDisCS LIRGs are systematically higher than those of the field LIRGs. A comparison with optical data reveals that ~80% of cluster LIRGs are blue and the remaining 20% lie on the red sequence. Of LIRGs with optical spectra, 88^{+4}_{-5}% show [O II] emission with EW([O II])>5A, and ~75% exhibit optical signatures of dusty starbursts. On average, the fraction of cluster LIRGs increases with projected cluster-centric radius but remains systematically lower than the field fraction over the area probed (< 1.5xR200). The amount of obscured star formation declines significantly over the 2.4 Gyr interval spanned by the EDisCS sample, and the rate of decline is the same for the cluster and field populations. Our results are consistent with an exponentially declining LIRG fraction, with the decline in the field delayed by ~1 Gyr relative to the clusters.

preprint2010arXiv

Stellar Population Variations in the Milky Way's Stellar Halo

If the stellar halos of disk galaxies are built up from the disruption of dwarf galaxies, models predict highly structured variations in the stellar populations within these halos. We test this prediction by studying the ratio of blue horizontal branch stars (BHB stars; more abundant in old, metal-poor populations) to main-sequence turn-off stars (MSTO stars; a feature of all populations) in the stellar halo of the Milky Way using data from the Sloan Digital Sky Survey. We develop and apply an improved technique to select BHB stars using ugr color information alone, yielding a sample of ~9000 g<18 candidates where ~70% of them are BHB stars. We map the BHB/MSTO ratio across ~1/4 of the sky at the distance resolution permitted by the absolute magnitude distribution of MSTO stars. We find large variations of BHB/MSTO star ratio in the stellar halo. Previously identified, stream-like halo structures have distinctive BHB/MSTO ratios, indicating different ages/metallicities. Some halo features, e.g., the low-latitude structure, appear to be almost completely devoid of BHB stars, whereas other structures appear to be rich in BHB stars. The Sagittarius tidal stream shows an apparent variation in BHB/MSTO ratio along its extent, which we interpret in terms of population gradients within the progenitor dwarf galaxy. Our detection of coherent stellar population variations between different stellar halo substructures provides yet more support to cosmologically motivated models for stellar halo growth.

preprint2010arXiv

Systematic errors in weighted 2-point correlation functions: An application to interaction-induced star formation

Weighted correlation functions are an increasingly important tool for understanding how galaxy properties depend on their separation from each other. We use a mock galaxy sample drawn from the Millenium simulation, assigning weights using a simple prescription to illustrate and explore how well a weighted correlation function recovers the true radial dependence of the input weights. We find that the use of a weighted correlation function results in a dilution of the magnitude of any radial dependence of properties and a smearing out of that radial dependence in radius, compared to the input behavior. We present a quantitative discussion of the dilution in the magnitude of radial dependence in properties in the special case of a constant enhancement at r < rc. In this particular case where there was a SFR enhancement at small radius r < rc = 35 kpc, the matching of one member of an enhanced pair with a non-enhanced galaxy in the same group gives an artificial enhancement out to large radius ~ 200 kpc. We compare this with observations of SFR enhancement from the SDSS (Li et al. 2008; MNRAS, 385, 1903) finding very similar behavior - a significant enhancement at radii < 40 kpc and a weak enhancement out to more than 150 kpc. While we explore a particular case in this Letter, it is easy to see that the phenomenon is general, and precision analyses of weighted correlation functions will need to account carefully for this effect using simulated mock catalogs.

preprint2010arXiv

Testing a prediction of the merger origin of early-type galaxies: a correlation between stellar populations and asymmetry

One of the key predictions of the merger hypothesis for the origin of early-type (elliptical and lenticular) galaxies is that tidally-induced asymmetric structure should correlate with signatures of a relatively young stellar population. Such a signature was found by Schweizer and Seitzer (1992; AJ, 104, 1039) at roughly 4sigma confidence. In this paper, we revisit this issue with a nearly ten-fold larger sample of 0.01<z<0.03 galaxies selected from the Two Micron All-Sky Survey and the Sloan Digital Sky Survey. We parameterize tidal structure using a repeatable algorithmic measure of asymmetry, and correlate this with color offset from the early-type galaxy color-magnitude relation. We recover the color offset-asymmetry correlation; furthermore, we demonstrate observationally for the first time that this effect is driven by a highly-significant trend towards younger ages at higher asymmetry values. We present a simple model for the evolution of early-type galaxies through gas-rich major and minor mergers that reproduces their observed build-up from z=1 to the present day and the distribution of present-day colors and ages. We show using this model that if both stellar populations and asymmetry were ideal `clocks' measuring the time since last major or minor gas-rich interaction, then we would expect a rather tight correlation between age and asymmetry. We suggest that the source of extra scatter is natural diversity in progenitor star formation history, gas content, and merger mass ratio, but quantitative confirmation of this conjecture will require sophisticated modeling. We conclude that the asymmetry-age correlation is in basic accord with the merger hypothesis, and indicates that an important fraction of the early-type galaxy population is affected by major or minor mergers at cosmologically-recent times.

preprint2010arXiv

The merger-driven evolution of massive galaxies

We explore the rate and impact of galaxy mergers on the massive galaxy population using the amplitude of the two-point correlation function on small scales for M > 5e10 M_sun galaxies from the COSMOS and COMBO-17 surveys. Using a pair fraction derived from the Sloan Digital Sky Survey as a low-redshift benchmark, the large survey area at intermediate redshifts allows us to determine the evolution of the close pair fraction with unprecedented accuracy for a mass-selected sample: we find that the fraction of galaxies more massive than 5e10M_sun in pairs separated by less than 30 kpc in 3D space evolves as F(z) = (0.0130+/-0.0019)x(1+z)^1.21+/-0.25 between z = 0 and z = 1.2. Assuming a merger time scale of 0.5 Gyrs, the inferred merger rate is such that galaxies with mass in excess of 1e11 M_sun have undergone, on average, 0.5 (0.7) mergers involving progenitor galaxies both more massive than 5e10 M_sun since z = 0.6 (1.2). We also study the number density evolution of massive red sequence galaxies using published luminosity functions and constraints on the M/L evolution from the fundamental plane. Moreover, we demonstrate that the measured merger rate of massive galaxies is sufficient to explain this observed number density evolution in massive red sequence galaxies since z = 1.

preprint2010arXiv

The Physical Origins of The Morphology-Density Relation: Evidence for Gas Stripping from the SDSS

We provide a physical interpretation and explanation of the morphology-density relation for galaxies, drawing on stellar masses, star formation rates, axis ratios and group halo masses from the Sloan Digital Sky Survey. We first re-cast the classical morphology-density relation in more quantitative terms, using low star formation rate (quiescence) as a proxy for early-type morphology and dark matter halo mass from a group catalog as a proxy for environmental density: for galaxies of a given stellar mass the quiescent fraction is found to increase with increasing dark matter halo mass. Our novel result is that - at a given stellar mass - quiescent galaxies are significantly flatter in dense environments, implying a higher fraction of disk galaxies. Supposing that the denser environments differ simply by a higher incidence of quiescent disk galaxies that are structurally similar to star-forming disk galaxies of similar mass, explains simultaneously and quantitatively these quiescence -nvironment and shape-environment relations. Our findings add considerable weight to the slow removal of gas as the main physical driver of the morphology-density relation, at the expense of other explanations.

preprint2009arXiv

History of Galaxy Interactions and their Impact on Star Formation over the Last 7 Gyr from GEMS

We perform a comprehensive estimate of the frequency of galaxy mergers and their impact on star formation over z~0.24--0.80 (lookback time T_b~3--7 Gyr) using 3698 (M*>=1e9 Msun) galaxies with GEMS HST, COMBO-17, and Spitzer data. Our results are: (1) Among 790 high mass (M*>=2.5e10 Msun) galaxies, the visually-based merger fraction over z~0.24--0.80, ranges from 9%+-5% to 8%+-2%. Lower limits on the major and minor merger fractions over this interval range from 1.1% to 3.5%, and 3.6% to 7.5%, respectively. This is the first approximate empirical estimate of the frequency of minor mergers at z<1. For a visibility timescale of ~0.5 Gyr, it follows that over T_b~3--7 Gyr, ~68% of high mass systems have undergone a merger of mass ratio >1/10, with ~16%, 45%, and 7% of these corresponding respectively to major, minor, and ambiguous `major or minor' mergers. The mean merger rate is a few x 1e-4 Gyr-1 Mpc-3. (2) We compare the empirical merger fraction and rate for high mass galaxies to a suite of Lambda CDM-based models: halo occupation distribution models, semi-analytic models, and hydrodynamic SPH simulations. We find qualitative agreement between observations and models such that the (major+minor) merger fraction or rate from different models bracket the observations, and show a factor of five dispersion. Near-future improvements can now start to rule out certain merger scenarios. (3) Among ~3698 M*>=1e9 Msun galaxies, we find that the mean SFR of visibly merging systems is only modestly enhanced compared to non-interacting galaxies over z~0.24--0.80. Visibly merging systems only account for less than 30% of the cosmic SFR density over T_b~3--7 Gyr. This suggests that the behavior of the cosmic SFR density over the last 7 Gyr is predominantly shaped by non-interacting galaxies.

preprint2009arXiv

Less than 10 percent of star formation in z=0.6 massive galaxies is triggered by major interactions

Both observations and simulations show that major tidal interactions or mergers between gas-rich galaxies can lead to intense bursts of starformation. Yet, the average enhancement in star formation rate (SFR) in major mergers and the contribution of such events to the cosmic SFR are not well estimated. Here we use photometric redshifts, stellar masses and UV SFRs from COMBO-17, 24 micron SFRs from Spitzer and morphologies from two deep HST cosmological survey fields (ECDFS/GEMS and A901/STAGES) to study the enhancement in SFR as a function of projected galaxy separation. We apply two-point projected correlation function techniques, which we augment with morphologically-selected very close pairs (separation <2 arcsec) and merger remnants from the HST imaging. Our analysis confirms that the most intensely star-forming systems are indeed interacting or merging. Yet, for massive (M* > 10^10 Msun) star-forming galaxies at 0.4<z<0.8, we find that the SFRs of galaxies undergoing a major interaction (mass ratios <1:4 and separations < 40 kpc) are only 1.80 +/- 0.30 times higher than the SFRs of non-interacting galaxies when averaged over all interactions and all stages of the interaction, in good agreement with other observational works. We demonstrate that these results imply that <10% of star formation at 0.4 < z < 0.8 is triggered directly by major mergers and interactions; these events are not important factors in the build-up of stellar mass since z=1.

preprint2009arXiv

Major Merging: The Way to Make a Massive, Passive Galaxy

We analyze the projected axial ratio distribution, p(b/a), of galaxies that were spectroscopically selected from the Sloan Digital Sky Survey (DR6) to have low star-formation rates. For these quiescent galaxies we find a rather abrupt change in p(b/a) at a stellar mass of ~10^{11} M_sol: at higher masses there are hardly any galaxies with b/a<0.6, implying that essentially none of them have disk-like intrinsic shapes and must be spheroidal. This transition mass is ~3-4 times higher than the threshold mass above which quiescent galaxies dominate in number over star-forming galaxies, which suggests these mass scales are unrelated. At masses lower than ~10^{11} M_sol, quiescent galaxies show a large range in axial ratios, implying a mix of bulge- and disk-dominated galaxies. Our result strongly suggests that major merging is the most important, and perhaps only relevant, evolutionary channel to produce massive (>10^{11} M_sol), quiescent galaxies, as it inevitably results in spheroids.

preprint2009arXiv

On the Size and Comoving Mass Density Evolution of Early-Type Galaxies

We present a simple, empirically motivated model that simultaneously predicts the evolution of the mean size and the comoving mass density of massive early-type galaxies from z=2 to the present. First we demonstrate that some size evolution of the population can be expected simply due to the continuous emergence of early-type galaxies. SDSS data reveal that in the present-day universe more compact early-type galaxies with a given dynamical mass have older stellar populations. In contrast, at a given stellar velocity dispersion, SDSS data show that there is no relation between size and age, which implies that the velocity dispersion can be used to estimate the epoch at which galaxies stopped forming stars, turning into early-type galaxies. Applying such a 'formation' criterion to a large sample of nearby early-type galaxies, we predict the redshift evolution in the size distribution and the comoving mass density. The resulting evolution in the mean size is roughly half of the observed evolution. Then we include a prescription for the merger histories of galaxies between the 'formation' redshift and the present, based on cosmological simulations of the assembly of dark matter halos. Such mergers after the transformation into an early-type galaxy are presumably dissipationless ('dry'), where the increase in size is expected to be approximately proportional to the increase in mass. This model successfully reproduces the observed evolution since z~2 in the mean size and in the comoving mass density of massive early-type galaxies. We conclude that the recently measured, substantial size evolution of early-type galaxies can be explained by the combined effect of the continuous emergence of galaxies as early types and their subsequent growth through dry merging.

preprint2009arXiv

Optically-passive spirals: The missing link in gradual star formation suppression upon cluster infall

Galaxies migrate from the blue cloud to the red sequence when their star formation is quenched. Here, we report on galaxies quenched by environmental effects and not by mergers or strong AGN as often invoked: They form stars at a reduced rate which is optically even less conspicuous, and manifest a transition population of blue spirals evolving into S0 galaxies. These 'optically passive' or 'red spirals' are found in large numbers in the STAGES project (and by Galaxy Zoo) in the infall region of clusters and groups.

preprint2009arXiv

Stellar mass-to-light ratios from galaxy spectra: how accurate can they be?

Stellar masses play a crucial role in the exploration of galaxy properties and the evolution of the galaxy population. In this paper, we explore the minimum possible uncertainties in stellar mass-to-light (M/L) ratios from the assumed star formation history (SFH) and metallicity distribution, with the goals of providing a minimum set of requirements for observational studies. We use a large Monte Carlo library of SFHs to study as a function of galaxy spectral type and signal-to-noise ratio (S/N) the statistical uncertainties of M/L values using either absorption-line data or broad band colors. The accuracy of M/L estimates can be significantly improved by using metal-sensitive indices in combination with age-sensitive indices, in particular for galaxies with intermediate-age or young stellar populations. While M/L accuracy clearly depends on the spectral S/N ratio, there is no significant gain in improving the S/N much above 50/pix and limiting uncertainties of 0.03 dex are reached. Assuming that dust is accurately corrected or absent and that the redshift is known, color-based M/L estimates are only slightly more uncertain than spectroscopic estimates (at comparable spectroscopic and photometric quality), but are more easily affected by systematic biases. This is the case in particular for galaxies with bursty SFHs (high Hdelta at fixed D4000), the M/L of which cannot be constrained any better than 0.15 dex with any indicators explored here. Finally, we explore the effects of the assumed prior distribution in SFHs and metallicity, finding them to be higher for color-based estimates.

preprint2009arXiv

The effect of dry mergers on the color-magnitude relation of early-type galaxies

We investigate the effect of dry merging on the color-magnitude relation (CMR) of galaxies and find that the amount of merging predicted by a hierarchical model results in a red sequence that compares well with the observed low-redshift relation. A sample of 29,000 early-type galaxies selected from the Sloan Digital Sky Survey Data Release 6 shows that the bright end of the CMR has a shallower slope and smaller scatter than the faint end. This magnitude dependence is predicted by a simple toy model in which gas-rich mergers move galaxies onto a "creation red sequence" (CRS) by quenching their star formation, and subsequent mergers between red, gas-poor galaxies (so-called "dry" mergers) move galaxies along the relation. We use galaxy merger trees from a semianalytic model of galaxy formation to test the amplitude of this effect and find a change in slope at the bright end that brackets the observations, using gas fraction thresholds of 10-30% to separate wet and dry mergers. A more realistic model that includes scatter in the CRS shows that dry merging decreases the scatter at the bright end. Contrary to previous claims, the small scatter in the observed CMR thus cannot be used to constrain the amount of dry merging.

preprint2009arXiv

The evolution of the scatter of the cosmic average color-magnitude relation: Demonstrating consistency with the ongoing formation of elliptical galaxies

We present first measurements of the evolution of the scatter of the cosmic average early-type galaxy color-magnitude relation (CMR) from z=1 to the present day, finding that it is consistent with models in which galaxies are constantly being added to the red sequence through truncation of star formation in blue cloud galaxies. We used a sample of over 700 red sequence, structurally-selected early-type galaxies (defined to have Sersic index >2.5) with redshifts 0<z<1 taken from the Extended Chandra Deep Field South (173 galaxies) and the Sloan Digital Sky Survey (550 galaxies), constructing rest-frame U-V colors accurate to <0.04mag. We find that the scatter of the CMR of cosmic average early-type galaxies is ~0.1mag in rest-frame U-V color at 0.05<z<0.75, and somewhat higher at z=1. We compared these observations with a model in which new red sequence galaxies are being constantly added at the rate required to match the observed number density evolution, and found that this model predicts the correct CMR scatter and its evolution. Furthermore, this model predicts approximately the correct number density of `blue spheroids' - structurally early-type galaxies with blue colors - albeit with considerable model dependence. Thus, we conclude that both the evolution of the number density and colors of the early-type galaxy population paint a consistent picture in which the early-type galaxy population grows significantly between z=1 and the present day through the quenching of star formation in blue cloud galaxies.

preprint2009arXiv

The SPLASH Survey: A Spectroscopic Analysis of the Metal-Poor, Low-Luminosity M31 dSph Satellite Andromeda X

Andromeda X (And X) is a newly discovered low-luminosity M31 dwarf spheroidal galaxy (dSph) found by Zucker et al. (2007) in the Sloan Digital Sky Survey (SDSS - York et al. 2000). In this paper, we present the first spectroscopic study of individual red giant branch stars in And X, as a part of the SPLASH Survey (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo). Using the Keck II telescope and multiobject DEIMOS spectrograph, we target two spectroscopic masks over the face of the galaxy and measure radial velocities for ~100 stars with a median accuracy of sigma_v ~ 3 km/s. The velocity histogram for this field confirms three populations of stars along the sight line: foreground Milky Way dwarfs at small negative velocities, M31 halo red giants over a broad range of velocities, and a very cold velocity ``spike'' consisting of 22 stars belonging to And X with v_rad = -163.8 +/- 1.2 km/s. By carefully considering both the random and systematic velocity errors of these stars (e.g., through duplicate star measurements), we derive an intrinsic velocity dispersion of just sigma_v = 3.9 +/- 1.2 km/s for And X, which for its size, implies a minimum mass-to-light ratio of M/L =37^{+26}_{-19} assuming the mass traces the light. Based on the clean sample of member stars, we measure the median metallicity of And X to be [Fe/H] = -1.93 +/- 0.11, with a slight radial metallicity gradient. The dispersion in metallicity is large, sigma([Fe/H]) = 0.48, possibly hinting that the galaxy retained much of its chemical enrichment products. We discuss the potential for better understanding the formation and evolution mechanisms for M31's system of dSphs through (current) kinematic and chemical abundance studies, especially in relation to the Milky Way sample. (abridged version)

preprint2008arXiv

Mapping low-latitude stellar substructure with SEGUE photometry

Encircling the Milky Way at low latitudes, the Low Latitude Stream is a large stellar structure, the origin of which is as yet unknown. As part of the SEGUE survey, several photometric scans have been obtained that cross the Galactic plane, spread over a longitude range of 50 to 203 degrees. These data allow a systematic study of the structure of the Galaxy at low latitudes, where the Low Latitude Stream resides. We apply colour-magnitude diagram fitting techniques to map the stellar (sub)structure in these regions, enabling the detection of overdensities with respect to smooth models. These detections can be used to distinguish between different models of the Low Latitude Stream, and help to shed light on the nature of the system.

preprint2007arXiv

An Explanation for the Observed Weak Size Evolution of Disk Galaxies

Surveys of distant galaxies with the Hubble Space Telescope and from the ground have shown that there is only mild evolution in the relationship between radial size and stellar mass for galactic disks from z~1 to the present day. Using a sample of nearby disk-dominated galaxies from the Sloan Digital Sky Survey (SDSS), and high redshift data from the GEMS (Galaxy Evolution from Morphology and SEDs) survey, we investigate whether this result is consistent with theoretical expectations within the hierarchical paradigm of structure formation. The relationship between virial radius and mass for dark matter halos in the LCDM model evolves by about a factor of two over this interval. However, N-body simulations have shown that halos of a given mass have less centrally concentrated mass profiles at high redshift. When we compute the expected disk size-stellar mass distribution, accounting for this evolution in the internal structure of dark matter halos and the adiabatic contraction of the dark matter by the self-gravity of the collapsing baryons, we find that the predicted evolution in the mean size at fixed stellar mass since z~1 is about 15-20 percent, in good agreement with the observational constraints from GEMS. At redshift z~2, the model predicts that disks at fixed stellar mass were on average only 60% as large as they are today. Similarly, we predict that the rotation velocity at a given stellar mass (essentially the zero-point of the Tully-Fisher relation) is only about 10 percent larger at z~1 (20 percent at z~2) than at the present day.

preprint2007arXiv

GEMS: Galaxy fitting catalogues and testing parametric galaxy fitting codes

In the context of measuring structure and morphology of intermediate redshift galaxies with recent HST/ACS surveys, we tune, test, and compare two widely used fitting codes (GALFIT and GIM2D) for fitting single-component Sersic models to the light profiles of both simulated and real galaxy data. We find that fitting accuracy depends sensitively on galaxy profile shape. Exponential disks are well fit with Sersic models and have small measurement errors, whereas fits to de Vaucouleurs profiles show larger uncertainties owing to the large amount of light at large radii. We find that both codes provide reliable fits and little systematic error, when the effective surface brightness is above that of the sky. Moreover, both codes return errors that significantly underestimate the true fitting uncertainties, which are best estimated with simulations. We find that GIM2D suffers significant systematic errors for spheroids with close companions owing to the difficulty of effectively masking out neighboring galaxy light; there appears to be no work around to this important systematic in GIM2D's current implementation. While this crowding error affects only a small fraction of galaxies in GEMS, it must be accounted for in the analysis of deeper cosmological images or of more crowded fields with GIM2D. In contrast, GALFIT results are robust to the presence of neighbors because it can simultaneously fit the profiles of multiple companions thereby deblending their effect on the fit to the galaxy of interest. We find GALFIT's robustness to nearby companions and factor of >~20 faster runtime speed are important advantages over GIM2D for analyzing large HST/ACS datasets. Finally we include our final catalog of fit results for all 41,495 objects detected in GEMS.

preprint2006arXiv

Andromeda X, A New Dwarf Spheroidal Satellite of M31: Photometry

We report the discovery of Andromeda X, a new dwarf spheroidal satellite of M31, based on stellar photometry from the Sloan Digital Sky Survey (SDSS). Using follow-up imaging data we have estimated its distance and other physical properties. We find that Andromeda X has a dereddened central surface brightness of mu_V,0 ~ 26.7 mag arcsec^-2 and a total apparent magnitude of V_tot ~ 16.1, which at the derived distance modulus, (m - M)_0 ~ 24.12 - 24.34, yields an absolute magnitude of M_V ~ -8.1 +/- 0.5; these values are quite comparable to those of Andromeda IX, a previously-discovered low luminosity M31 satellite. The discoveries of Andromeda IX and Andromeda X suggest that such extremely faint satellites may be plentiful in the Local Group.

preprint2006arXiv

The merger rate of massive galaxies

We calculate the projected two point correlation function for samples of luminous and massive galaxies in the COMBO-17 photometric redshift survey, focusing particularly on the amplitude of the correlation function at small projected radii and exploring the constraints such measurements can place on the galaxy merger rate. For nearly volume-limited samples with 0.4<z<0.8, we find that 4+/-1% of luminous M_B<-20 galaxies are in close physical pairs (with real space separation of <30 proper kpc). The corresponding fraction for massive galaxies with M_*>2.5e10 M_sun is 5+/-1%. Incorporating close pair fractions from the literature, the 2dFGRS and the SDSS, we find a fairly rapid evolution of the merger fraction of massive galaxies between z=0.8 and the present day. Assuming that the major merger timescale is of order the dynamical timescale for close massive galaxy pairs, we tentatively infer that ~50% (70%) of all galaxies with present-day masses M_*>5e10 M_sun (remnants of mergers between galaxies with M_*>2.5e10 M_sun) have undergone a major merger since z=0.8(1): major mergers between massive galaxies are a significant driver of galaxy evolution over the last eight billion years.

preprint2005arXiv

GEMS: The Size Evolution of Disk Galaxies

We combine HST imaging from the GEMS survey with photometric redshifts from COMBO-17 to explore the evolution of disk-dominated galaxies since z<1.1. The sample is comprised of all GEMS galaxies with Sersic indices n<2.5, derived from fits to the galaxy images. We account fully for selection effects through careful analysis of image simulations; we are limited by the depth of the redshift and HST data to the study of galaxies with absolute magnitudes M(V)<-20, or equivalently stellar masses log(M)>10. We find strong evolution in the magnitude-size scaling relation for galaxies with M(V)<-20, corresponding to a brightening of 1 mag per sqarcsec in rest-frame V-band by z=1. Yet, disks at a given absolute magnitude are bluer and have lower stellar mass-to-light ratios at z=1 than at the present day. As a result, our findings indicate weak or no evolution in the relation between stellar mass and effective disk size for galaxies with log(M)>10 over the same time interval. This is strongly inconsistent with the most naive theoretical expectation, in which disk size scales in proportion to the halo virial radius, which would predict that disks are a factor of two denser at fixed mass at z=1. The lack of evolution in the stellar mass-size relation is consistent with an ``inside-out'' growth of galaxy disks on average (galaxies increasing in size as they grow more massive), although we cannot rule out more complex evolutionary scenarios.

preprint2004arXiv

Andromeda IX: A New Dwarf Spheroidal Satellite of M31

We report the discovery of a new dwarf spheroidal satellite of M31, Andromeda IX, based on resolved stellar photometry from the Sloan Digital Sky Survey (SDSS). Using both SDSS and public archival data we have estimated its distance and other physical properties, and compared these to the properties of a previously known dwarf spheroidal companion, Andromeda V, also observed by SDSS. Andromeda IX is the lowest surface brightness galaxy found to date (mu_{V} ~ 26.8 mag arcsec^-2), and at the distance we estimate from the position of the tip of Andromeda IX's red giant branch, (m - M)_{0} ~ 24.5 (805 kpc), Andromeda IX would also be the faintest galaxy known (M_{V} ~ -8.3).

preprint2004arXiv

GEMS: Galaxy Evolution from Morphologies and SEDs

GEMS, Galaxy Evolution from Morphologies and SEDs, is a large-area (800 arcmin2) two-color (F606W and F850LP) imaging survey with the Advanced Camera for Surveys on HST. Centered on the Chandra Deep Field South, it covers an area of ~28'x28', or about 120 Hubble Deep Field areas, to a depth of m_AB(F606W)=28.3 (5sigma and m_AB(F850LP)=27.1 (5sigma) for compact sources. In its central ~1/4, GEMS incorporates ACS imaging from the GOODS project. Focusing on the redshift range 0.2<=z<=1.1, GEMS provides morphologies and structural parameters for nearly 10,000 galaxies where redshift estimates, luminosities and SEDs exist from COMBO-17. At the same time, GEMS contains detectable host galaxy images for several hundred faint AGN. This paper provides an overview of the science goals, the experiment design, the data reduction and the science analysis plan for GEMS.

preprint2003arXiv

The Optical and Near-Infrared Properties of Galaxies: I. Luminosity and Stellar Mass Functions

We use a large sample of galaxies from the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) to calculate galaxy luminosity and stellar mass functions in the local Universe. We estimate k-corrections, evolution corrections, and stellar mass-to-light ratios (M/Ls) by fitting the galaxy colors with simple models. Our optical and near-infrared luminosity functions agree with most recent literature optical and near-infrared determinations within the uncertainties. We argue that 2MASS is biased against low surface brightness galaxies, and use SDSS plus our knowledge of stellar populations to estimate that the true K-band luminosity function has a steeper faint end slope than the direct estimate. Assuming a universally-applicable stellar initial mass function (IMF), we derive the stellar mass function of galaxies. The faint end slope slope for the stellar mass function is steeper than -1.1, reflecting the low stellar M/Ls characteristic of low-mass galaxies. We estimate an upper limit to the stellar mass density in the local Universe Omega* h = 0.0020+/-0.0006 by assuming an IMF as rich in low-mass stars as allowed by observations of galaxy dynamics. Finally, we find that the characteristic luminosity or mass of early-type galaxies is larger than for later types, and the faint end slope is steeper for later types than for earlier types. Accounting for typing uncertainties, we estimate that at least half, and perhaps as much as 3/4, of the stellar mass in the Universe is in early-type galaxies. We present also SDSS/2MASS color-M/L correlations, an updated discussion of near-infrared stellar M/L estimates, and the volume-corrected distribution of g and K-band stellar M/Ls as a function of stellar mass. [Abridged]

preprint2000arXiv

Stellar mass-to-light ratios and the Tully-Fisher relation

We have used a suite of simplified spectrophotometric spiral galaxy evolution models to argue that there are substantial variations in stellar mass-to-light ratios (M/Ls) within and among galaxies, amounting to factors of between 3 and 7 in the optical, and factors of 2 in the near-infrared. Under the assumption of a universal spiral galaxy IMF, our models show a strong and robust correlation between stellar M/L and the optical colors of the integrated stellar populations. Using observed maximum disk M/Ls we conclude that a Salpeter Initial Mass Function (IMF) has too many low-mass stars per unit luminosity, but that an IMF similar to the Salpeter IMF at the high-mass end with less low-mass stars (giving stellar M/Ls 30% or more lower than the Salpeter value) is consistent with the maximum disk constraints. We apply the model trends in stellar M/L with color to the Tully-Fisher (TF) relation. We find that the stellar mass TF relation is relatively steep and has modest scatter, and is independent of the passband and color used to derive the stellar masses. The difference in slope between the optical (especially blue) and near-infrared TF relations is due to the combined effects of dust attenuation and stellar M/L variations with galaxy mass. We find that the baryonic TF relation has a slope of 3.5 (with random and systematic errors of +/- 0.2 each). Since we have normalized the stellar M/L to be as high as can possibly be allowed by maximum disk constraints, the slope of the baryonic TF relation will be somewhat shallower than 3.5 if all disks are substantially sub-maximal. [Abridged]

Eric F. Bell

What is connected

Connect this record

See the researcher in context

Building this map preview

143 published item(s)

After The Fall: Resolving the Molecular Gas in Post-Starburst Galaxies

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

New Globular Cluster Candidates in the M81 group

R2-D2: Roman and Rubin -- From Data to Discovery

Relating the Diverse Merger Histories and Satellite Populations of Nearby Galaxies

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

Uncertainties associated with the backward integration of dwarf satellites using simple parametric potentials

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

Is NGC 300 a pure exponential disk galaxy?

Mass-to-Light Ratios of Spatially Resolved Stellar Populations in M31

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

SMASHing the low surface brightness SMC

The Second Data Release of the Survey of the MAgellanic Stellar History (SMASH)

Tracing the anemic stellar halo of M101

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

A Synoptic Map of Halo Substructures from the Pan-STARRS1 3π Survey

ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description

Detailed Chemical Abundances in the r-Process-Rich Ultra-Faint Dwarf Galaxy Reticulum 2

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

First results from the MADCASH Survey: A Faint Dwarf Galaxy Companion to the Low Mass Spiral Galaxy NGC 2403 at 3.2 Mpc

Is Draco II one of the faintest dwarf galaxies? First study from Keck/DEIMOS spectroscopy

Mapping the Monoceros Ring in 3D with Pan-STARRS1

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

Panchromatic Hubble Andromeda Treasury XVI. Star Cluster Formation Efficiency and the Clustered Fraction of Young Stars

Quiescence correlates strongly with directly-measured black hole mass in central galaxies

The Bursty Star Formation Histories of Low-mass Galaxies at $0.4<z<1$ Revealed by Star Formation Rates Measured from H$β$ and FUV

The DESI Experiment Part I: Science,Targeting, and Survey Design

The DESI Experiment Part II: Instrument Design

The Diversity of Growth Histories of Milky Way-mass Galaxies

The GHOSTS survey. II. The diversity of Halo Color and Metallicity Profiles of Massive Disk Galaxies

The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

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

The Stellar Density Profile of the Distant Galactic Halo

Triangulum II: a very metal-poor and dynamically hot stellar system

A clear age-velocity dispersion correlation in Andromeda's stellar disk

A Deep Study of the Dwarf Satellites Andromeda XXVIII & Andromeda XXIX

A Global Star Forming Episode in M31 2-4 Gyr Ago

A New Faint Milky Way Satellite Discovered in the Pan-STARRS1 3 pi Survey

And yet it moves: The dangers of artificially fixing the Milky Way center of mass in the presence of a massive Large Magellanic Cloud

Are Compton-Thick AGN the Missing Link Between Mergers and Black Hole Growth?

CANDELS Visual Classifications: Scheme, Data Release, and First Results

Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3

Finding, characterizing and classifying variable sources in multi-epoch sky surveys: QSOs and RR Lyrae in PS1 3$π$ data

Hydra II: a faint and compact Milky Way dwarf galaxy found in the Survey of the Magellanic Stellar History

Linking the Structural Properties of Galaxies and their Star Formation Histories with STAGES

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

Quenching and Morphological Transformation in Semi-Analytic Models and CANDELS

Sagittarius II, Draco II and Laevens 3: three new Milky Way satellites discovered in the Pan-STARRS 1 3pi Survey

Stellar Populations of Barred Quiescent Galaxies

The High-Mass Stellar Initial Mass Function in M31 Clusters

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

A new distant Milky Way globular cluster in the Pan-STARRS1 3π survey

A new population of recently quenched elliptical galaxies in the SDSS

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

Bulgeless Galaxies at Intermediate Redshift: Sample Selection, Colour Properties, and the Existence of Powerful AGN

Galactic Globular and Open Cluster Fiducial Sequences in the Pan-STARRS1 Photometric System

Galaxy Zoo: Are Bars Responsible for the Feeding of Active Galactic Nuclei at 0.2 < z < 1.0?

Galaxy Zoo: CANDELS Barred Disks and Bar Fractions

No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~2 in CANDELS / 3D-HST

Serendipitous Discovery of a Thin Stellar Stream near the Galactic Bulge in the Pan-STARRS1 3Pi Survey

Spectroscopy of the three distant Andromedan satellites Cassiopeia III, Lacerta I, and Perseus I

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

The Complex Structure of Stars in the Outer Galactic Disk as revealed by Pan-STARRS1

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

The Mass Dependence of Dwarf Satellite Galaxy Quenching

The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy

The Panchromatic Hubble Andromeda Treasury V: Ages and Masses of the Year 1 Stellar Clusters