Source author record

Hans-Walter Rix

Hans-Walter Rix appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

187works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

187 published item(s)

preprint2026arXiv

Gaia Barium Dwarfs and Their Ostensibly Ordinary Counterparts

The recently identified Gaia population of main-sequence--white dwarf (MS+WD) binaries at separations of ${\sim}\,1~{\rm AU}$, often with moderate eccentricities, is not readily reproduced by binary population synthesis models. Barium stars represent a closely related population whose enrichment in $s$-process elements both confirms the presence of a WD companion and attests to past binary interaction. It also indicates that mass transfer occurred at least during the late- and post-AGB phases of the WD progenitor, when $s$-process elements are dredged up. In this work, we further explore the connection between the astrometrically identified Gaia MS+WD binaries and the classical barium star population. To achieve this, we used high-resolution FEROS spectroscopy to measure abundances for 30 Gaia DR3 non-single-star binaries, identifying 10 as Ba-enriched. Together with our recent analysis of archival GALAH data, this yields a sample of 38 barium dwarfs with dynamically measured WD masses, compared to only 6 previously known systems with known WD masses at these separations. We find that, in cases where metallicity is sufficiently low to facilitate efficient $s$-process production, barium and yttrium enrichment is often detected. This enrichment is also identified in eccentric systems, suggesting that post-AGB mass transfer mechanisms are capable of pumping eccentricity into the orbit or occur without erasing it. Our results indicate that the Gaia MS+WD binaries trace the population from which barium stars emerge. Treating the large Gaia-discovered population as an extension of known $s$-process enriched dwarfs opens an avenue to empirically constrain their formation and evolution.

preprint2025arXiv

Millions of Main-Sequence Binary Stars from Gaia BP/RP Spectra

We present the main-sequence binary (MSMS) Catalog derived from Gaia Data Release 3 BP/RP (XP) spectra. Leveraging the vast sample of low-resolution Gaia XP spectra, we develop a forward modeling approach that maps stellar mass and photometric metallicity to XP spectra using a neural network. Our methodology identifies binary systems through statistical comparison of single- and binary-star model fits, enabling detection of binaries with mass ratios between 0.4 and 1.0 and flux ratios larger than 0.1. From an initial sample of 35 million stars within 1 kpc, we identify 14 million binary candidates and define a high-confidence "golden sample" of 1 million binary systems. This large, homogeneous sample enables detailed statistical analysis of binary properties across diverse Galactic environments, providing new insights into binary star formation and evolution. In addition, the $χ^2$ comparison allows us to distinguish stars with luminous companions from single stars or binaries with dark companions, such as white dwarfs, neutron stars and black hole candidates, improving our understanding of compact object populations.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2022arXiv

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc.

preprint2022arXiv

An empirical model of the Gaia DR3 selection function

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky. On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project: $\texttt{https://github.com/gaia-unlimited/gaiaunlimited}$

preprint2022arXiv

Analysis of high-resolution FEROS spectroscopy for a sample of variable B-type stars assembled from TESS photometry

Spectroscopic data are necessary to break degeneracies in asteroseismic modelling of the interior structure in high- and intermediate-mass stars. We derive precise photospheric stellar parameters for a sample of 166 B-type stars with TESS light curves, suitable for detailed asteroseismic studies, through a homogeneous spectroscopic analysis. The variability types of these stars are classified based on all currently available TESS sectors. We obtained high-resolution spectra for all 166 targets with the FEROS spectrograph in the context of a large program. The Least-Squares Deconvolution method is employed to investigate spectral line profile variability and to detect binary systems. We identify 26 spectroscopic double-lined binaries; the remainder of the sample are 42 supergiants in the LMC galaxy and 98 Galactic stars. The spectra of the Galactic stars are analysed with the zeta-Payne, a machine learning-based spectrum analysis algorithm. We determine the five main surface parameters with average formal precisions of 70 K (Teff), 0.03 dex (logg), 0.07 dex ([M/H]), 8 km/s (vsini), and 0.7 km/s (vmicro). The average internal uncertainties we find for FEROS spectra with our spectrum analysis method are 430 K (Teff), 0.12 dex (logg), 0.13 dex ([M/H]), 12 km/s (vsini), and 2 km/s (vmicro). We find spectroscopic evidence that eight of the 98 Galactic targets are fast rotating g-mode pulsators occurring in between the slowly pulsating B (SPB) stars and delta Scuti instability strips. The g-mode frequencies of these pulsators are shifted to relatively high frequency values due to their rotation. Their apparently too low Teff relative to the SPB instability region can in most cases be explained by the gravity darkening effect. We also discover 13 new HgMn stars of which only one is found in a spectroscopic binary, resulting in a biased and therefore unreliable low binary rate of only 8%.

preprint2022arXiv

Pearls on a String: Numerous Stellar Clusters Strung along the Same Orbit

Stars originate from the dense interstellar medium, which exhibits filamentary structure to scales of $\sim 1$ kpc in galaxies like our Milky Way. We explore quantitatively how much resulting large-scale correlation there is among different stellar clusters and associations in \emph{orbit phase space}, characterized here by actions and angles. As a starting point, we identified 55 prominent stellar overdensities in the 6D space of orbit (actions) and orbital phase (angles), among the $\sim$ 2.8 million stars with radial velocities from Gaia EDR3 and $d \leq 800$ pc. We then explored the orbital \emph{phase} distribution of all sample stars in the same \emph{orbit} patch as any one of these 55 overdensities. We find that very commonly numerous other distinct orbital phase overdensities exist along these same orbits, like pearls on a string. These `pearls' range from known stellar clusters to loose, unrecognized associations. Among orbit patches defined by one initial orbit-phase overdensity 50\% contain at least 8 additional orbital-phase pearls of 10 cataloged members; 20\% of them contain 20 additional pearls. This is in contrast to matching orbit patches sampled from a smooth mock catalog, or offset nearby orbit patches, where there are only 2 (or 5, respectively) comparable pearls. Our findings quantify for the first time how common it is for star clusters and associations to form at distinct orbital phases of nearly the same orbit. This may eventually offer a new way to probe the 6D orbit structure of the filamentary interstellar medium.

preprint2022arXiv

Stellar Abundance Maps of the Milky Way Disk

To understand the formation of the Milky Way's prominent bar it is important to know whether stars in the bar differ in the chemical element composition of their birth material as compared to disk stars. This requires stellar abundance measurements for large samples across the Milky Way's body. Such samples, e.g. luminous red giant stars observed by SDSS's Apogee survey, will inevitably span a range of stellar parameters; as a consequence, both modelling imperfections and stellar evolution may preclude consistent and precise estimates of their chemical composition at a level of purported bar signatures, which has left current analyses of a chemically distinct bar inconclusive. Here, we develop a new self-calibration approach to eliminate both modelling and astrophysical abundance systematics among red giant branch (RGB) stars of different luminosities (and hence surface gravity $\log g$). We apply our method to $48,853$ luminous Apogee DR16 RGB stars to construct spatial abundance maps of $20$ chemical elements near the Milky Way's mid-plane, covering Galactocentric radii of $0\,{\rm kpc}<R_{\rm GC}<20\,\rm kpc$. Our results indicate that there are no abundance variations whose geometry matches that of the bar, and that the mean abundance gradients vary smoothly and monotonically with Galactocentric radius. We confirm that the high-$α$ disk is chemically homogeneous, without spatial gradients. Furthermore, we present the most precise [Fe/H] vs. $R_{\rm GC}$ gradient to date with a slope of $-0.057\pm0.001\rm~dex\,kpc^{-1}$ out to approximately $15$ kpc.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2022arXiv

Unicorns and Giraffes in the binary zoo: stripped giants with subgiant companions

We analyze two binary systems containing giant stars, V723 Mon ("the Unicorn") and 2M04123153+6738486 ("the Giraffe"). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modeling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\,giant} = 3,800\,\rm K$ and $4,000\,\rm K$, $R_{\rm giant}= 22.5\,R_{\odot}$ and $25\,R_{\odot}$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\,M_{\odot}$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\,2} = 5,800\,\rm K$ and $5,150\,\rm K$, $R_2= 8.3\,R_{\odot}$ and $9\,R_{\odot}$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8\,M_{\odot}$ and $\approx 1.8\,M_{\odot}$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1-2\,M_{\odot}$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely-observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs.

preprint2022arXiv

What are the spectroscopic binaries with high mass functions near the Gaia DR3 main sequence?

The 3rd data release of the Gaia mission includes orbital solutions for $> 10^5$ single-lined spectroscopic binaries, representing more than an order of magnitude increase in sample size over all previous studies. This dataset is a treasure trove for searches for quiescent black hole + normal star binaries. We investigate one population of black hole candidate binaries highlighted in the data release: sources near the main sequence in the color-magnitude diagram (CMD) with dynamically-inferred companion masses $M_2$ larger than the CMD-inferred mass of the luminous star. We model light curves, spectral energy distributions, and archival spectra of the 14 such objects in DR3 with high-significance orbital solutions and inferred $M_2 > 3\,M_{\odot}$. We find that 100\% of these sources are mass-transfer binaries containing a highly stripped lower giant donor ($0.2 \lesssim M/M_{\odot} \lesssim 0.4$) and much more massive ($2 \lesssim M/M_{\odot} \lesssim 2.5$) main-sequence accretor. The Gaia orbital solutions are for the donors, which contribute about half the light in the Gaia RVS bandpass but only $\lesssim 20\%$ in the $g-$band. The accretors' broad spectral features likely prevented the sources from being classified as double-lined. The donors are all close to Roche lobe-filling ($R/R_{\rm Roche\,lobe}>0.8$), but modeling suggests that a majority are detached ($R/R_{\rm Roche\,lobe}<1$). Binary evolution models predict that these systems will soon become detached helium white dwarf + main sequence "EL CVn" binaries. Our investigation highlights both the power of Gaia data for selecting interesting sub-populations of binaries and the ways in which binary evolution can bamboozle standard CMD-based stellar mass estimates.

preprint2022arXiv

Wide twin binaries are extremely eccentric: evidence of twin binary formation in circumbinary disks

The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separations of the widest twins. Here, we infer the eccentricity distribution of wide twins from the distribution of their $v$-$r$ angles, i.e., the angle between the components' separation and relative velocity vectors. We find that wide twins must be on extremely eccentric orbits. For the excess-twin population at 400-1000 AU, we infer a near-delta function excess of high-eccentricity system, with eccentricity $0.95 \lesssim e \leq 1$. These high eccentricities for wide twins imply pericenter distances of order $10$ AU and suggest that their orbits were scattered via dynamical interactions in their birth environments, consistent with a scenario in which twins are born in circumbinary disks and subsequently widened. These results further establish twin wide binaries as a distinct population and imply that wide twins can be used as a probe of the dynamical history of stellar populations.

preprint2022arXiv

Zeta-Payne: a fully automated spectrum analysis algorithm for the Milky Way Mapper program of the SDSS-V survey

The Sloan Digital Sky Survey has recently initiated its 5th survey generation (SDSS-V), with a central focus on stellar spectroscopy. In particular, SDSS-V Milky Way Mapper program will deliver multi-epoch optical and near-infrared spectra for more than 5 million stars across the entire sky, covering a large range in stellar mass, surface temperature, evolutionary stage, and age. About 10% of those spectra will be of hot stars of OBAF spectral types, for whose analysis no established survey pipelines exist. Here we present the spectral analysis algorithm, Zeta-Payne, developed specifically to obtain stellar labels from SDSS-V spectra of stars with these spectral types and drawing on machine learning tools. We provide details of the algorithm training, its test on artificial spectra, and its validation on two control samples of real stars. Analysis with Zeta-Payne leads to only modest internal uncertainties in the near-IR with APOGEE (optical with BOSS): 3-10% (1-2%) for Teff, 5-30% (5-25%) for v*sin(i), 1.7-6.3 km/s(0.7-2.2 km/s) for RV, $<0.1$ dex ($<0.05$ dex) for log(g), and 0.4-0.5 dex (0.1 dex) for [M/H] of the star, respectively. We find a good agreement between atmospheric parameters of OBAF-type stars when inferred from their high- and low-resolution optical spectra. For most stellar labels the APOGEE spectra are (far) less informative than the BOSS spectra of these stars, while log(g), v*sin(i), and [M/H] are in most cases too uncertain for meaningful astrophysical interpretation. This makes BOSS low-resolution optical spectra better for stellar labels of OBAF-type stars, unless the latter are subject to high levels of extinction.

preprint2021arXiv

A million binaries from Gaia eDR3: sample selection and validation of Gaia parallax uncertainties

We construct from Gaia eDR3 an extensive catalog of spatially resolved binary stars within $\approx$ 1 kpc of the Sun, with projected separations ranging from a few au to 1 pc. We estimate the probability that each pair is a chance alignment empirically, using the Gaia catalog itself to calculate the rate of chance alignments as a function of observables. The catalog contains 1.3 (1.1) million binaries with >90% (>99%) probability of being bound, including 16,000 white dwarf -- main sequence (WD+MS) binaries and 1,400 WD+WD binaries. We make the full catalog publicly available, as well as the queries and code to produce it. We then use this sample to calibrate the published Gaia DR3 parallax uncertainties, making use of the binary components' near-identical parallaxes. We show that these uncertainties are generally reliable for faint stars ($G\gtrsim 18$), but are underestimated significantly for brighter stars. The underestimates are generally $\le 30\%$ for isolated sources with well-behaved astrometry, but are larger (up to 80%) for apparently well-behaved sources with a companion within $\lesssim 4$ arcsec, and much larger for sources with poor astrometric fits. We provide an empirical fitting function to inflate published $σ_{\varpi}$ values for isolated sources. The public catalog offers wide ranging follow-up opportunities: from calibrating spectroscopic surveys, to precisely constraining ages of field stars, to the masses and the initial-final mass relation of white dwarfs, to dynamically probing the Galactic tidal field.

preprint2021arXiv

Stellar labels for hot stars from low-resolution spectra - I. the HotPayne method and results for 330,000 stars from LAMOST DR6

We set out to determine stellar labels from low-resolution survey spectra of hot, OBA stars with effective temperature (Teff) higher than 7500K. This fills a gap in the scientific analysis of large spectroscopic stellar surveys such as LAMOST, which offers spectra for millions of stars at R=1800. We first explore the theoretical information content of such spectra for determining stellar labels, via the Cramér-Rao bound. We show that in the limit of perfect model spectra and observed spectra with S/N of 100, precise estimates are possible for a wide range of stellar labels: not only the effective temperature Teff, surface gravity logg, and projected rotation velocity vsini, but also the micro-turbulence velocity, Helium abundance and the elemental abundances [C/H], [N/H], [O/H], [Si/H], [S/H], and [Fe/H]. Our analysis illustrates that the temperature regime of around 9500K is challenging, as the dominant Balmer and Paschen line strength vary little with Teff. We implement the simultaneous fitting of these 11 stellar labels to LAMOST hot-star spectra using the Payne approach, drawing on Kurucz's ATLAS12/SYNTHE LTE spectra as the underlying models. We then obtain stellar parameter estimates for a sample of about 330,000 hot stars with LAMOST spectra, an increase by about two orders of magnitude in sample size. Among them, about 260,000 have good Gaia parallaxes (S/N>5), and more than 95 percent of them are luminous stars, mostly on the main sequence; the rest reflects lower luminosity evolved stars, such as hot subdwarfs and white dwarfs. We show that the fidelity of the abundance estimates is limited by the systematics of the underlying models, as they do not account for NLTE effects. Finally, we show the detailed distribution of vsini of stars with 8000-15,000K, illustrating that it extends to a sharp cut-off at the critical rotation velocity, across a wide range of temperatures.

preprint2020arXiv

A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

preprint2020arXiv

A Simple Unified Spectroscopic Indicator of Stellar Luminosity: the Extended Flux Weighted Gravity-Luminosity Relationship

We show that for a wide range of stellar masses, from 0.3 to 20 Msun, and for evolutionary phases from the main sequence to the beginning of the red giant stage, the stellar flux weighted gravity, g_F ~ g/Teff^4, is tightly correlated with absolute bolometric magnitude Mbol. Such a correlation is predicted by stellar evolution theory. We confirm this relation observationally, using a sample of 445 stars with precise stellar parameters. It holds over 17 stellar magnitudes from Mbol = 9.0 mag to -8.0 mag with a scatter of 0.17 mag above Mbol = -3.0 and 0.29 mag below this value. We then test the relation with 2.2 million stars with 6.5 mag > Mbol > 0.5 mag, where 'mass-produced' but robust log g, Teff and Mbol from LAMOST DR5 and Gaia DR2 are available. We find that the same relation holds with a scatter of ~0.2 mag for single stars offering a simple spectroscopic distance estimate good to ~10 percent.

preprint2020arXiv

Chemically peculiar A and F stars with enhanced s-process and iron-peak elements: stellar radiative acceleration at work

We present $\gtrsim 15,000$ metal-rich (${\rm [Fe/H]}>-0.2$dex) A and F stars whose surface abundances deviate strongly from Solar abundance ratios and cannot plausibly reflect their birth material composition. These stars are identified by their high [Ba/Fe] abundance ratios (${\rm [Ba/Fe]}>1.0$dex) in the LAMOST DR5 spectra analyzed by Xiang et al. (2019). They are almost exclusively main sequence and subgiant stars with $T_{\rm eff}\gtrsim6300$K. Their distribution in the Kiel diagram ($T_{\rm eff}$--$\log g$) traces a sharp border at low temperatures along a roughly fixed-mass trajectory (around $1.4M_\odot)$ that corresponds to an upper limit in convective envelope mass fraction of around $10^{-4}$. Most of these stars exhibit distinctly enhanced abundances of iron-peak elements (Cr, Mn, Fe, Ni) but depleted abundances of Mg and Ca. Rotational velocity measurements from GALAH DR2 show that the majority of these stars rotate slower than typical stars in an equivalent temperature range. These characteristics suggest that they are related to the so-called Am/Fm stars. Their abundance patterns are qualitatively consistent with the predictions of stellar evolution models that incorporate radiative acceleration, suggesting they are a consequence of stellar internal evolution particularly involving the competition between gravitational settling and radiative acceleration. These peculiar stars constitute 40% of the whole population of stars with mass above 1.5$M_\odot$, affirming that "peculiar" photospheric abundances due to stellar evolution effects are a ubiquitous phenomenon for these intermediate-mass stars. This large sample of Ba-enhanced chemically peculiar A/F stars with individual element abundances provides the statistics to test more stringently the mechanisms that alter the surface abundances in stars with radiative envelopes.

preprint2020arXiv

Close Binary Companions to APOGEE DR16 Stars: 20,000 Binary-star Systems Across the Color-Magnitude Diagram

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.

preprint2020arXiv

From birth associations to field stars: mapping the small-scale orbit distribution in the Galactic disc

Stars born at the same time in the same place should have formed from gas of the same element composition. But most stars subsequently disperse from their birth siblings, in orbit and orbital phase, becoming 'field stars'. Here we explore and provide direct observational evidence for this process in the Milky Way disc, by quantifying the probability that orbit-similarity among stars implies indistinguishable metallicity. We define the orbit similarity among stars through their distance in action-angle space, $Δ(J,θ)$, and their abundance similarity simply by $Δ$[Fe/H]. Analyzing a sample of main sequence stars from Gaia DR2 and LAMOST, we find an excess of pairs with the same metallicity ($Δ\mathrm{[Fe/H]}<0.1$) that extends to remarkably large separations in $Δ(J,θ)$ that correspond to nearly 1 kpc distances. We assess the significance of this effect through a mock sample, drawn from a smooth and phase-mixed orbit distribution. Through grouping such star pairs into associations with a friend-of-friends algorithm linked by $Δ(J,θ)$, we find 100s of mono-abundance groups with $\ge 3$ (to $\gtrsim 20$) members; these groups -- some clusters, some spread across the sky -- are over an order-of-magnitude more abundant than expected for a smooth phase-space distribution, suggesting that we are witnessing the 'dissolution' of stellar birth associations into the field.

preprint2020arXiv

Keeping it Cool: Much Orbit Migration, yet Little Heating, in the Galactic Disk

A star in the Milky Way's disk can now be at a Galactocentric radius quite distant from its birth radius for two reasons: either its orbit has become eccentric through radial heating, which increases its radial action $J_R$ (`blurring'); or merely its angular momentum $L_z$ has changed and thereby its guiding radius (`churning'). We know that radial orbit migration is strong in the Galactic low-$α$ disk and set out to quantify the relative importance of these two effects, by devising and applying a parameterized model for the distribution $p(L_z, J_R, τ, \mathrm[Fe/H])$ in the stellar disk. This model describes the orbit evolution for stars of age $τ$ and metallicity [Fe/H], presuming coeval stars were initially born on (near-)circular orbits, and with a unique [Fe/H] at a given birth angular momentum and age. We fit this model to APOGEE red clump stars, accounting for the complex selection function of the survey. The best fit model implies changes of angular momentum of $\sqrt{\langle ΔL_z \rangle^2} \approx 619\, \mathrm{kpc~km/s~}(τ/\mathrm{6~Gyr})^{0.5}$, and changes of radial action as $\sqrt{\langle ΔJ_R \rangle^2} \approx 63\, \mathrm{kpc~km/s~} (τ/\mathrm{6~Gyr})^{0.6}$ at 8 kpc. This suggests that the secular orbit evolution of the disk is dominated by diffusion in angular momentum, with radial heating being an order of magnitude lower.

preprint2020arXiv

Not all stars form in clusters -- $Gaia$-DR2 uncovers the origin of OB associations

Historically, it has often been asserted that most stars form in compact clusters. In this scenario, present-day gravitationally-unbound OB associations are the result of the expansion of initially gravitationally-bound star clusters. However, this paradigm is inconsistent with recent results, both theoretical and observational, that instead favour a hierarchical picture of star formation in which stars are formed across a continuous distribution of gas densities and most OB associations never were bound clusters. Instead they are formed in-situ as the low-density side of this distribution, rather than as the remnants of expanding clusters. We utilise the second $Gaia$ data release to quantify the degree to which OB associations are undergoing expansion and, therefore, whether OB associations are the product of expanding clusters, or whether they were born in-situ, as the large-scale, globally-unbound associations that we see today. We find that the observed kinematic properties of associations are consistent with highly substructured velocity fields and additionally require some degree of localised expansion from sub-clusters within the association. While most present-day OB associations do exhibit low levels of expansion, there is no significant correlation between radial velocity and radius. Therefore, the large-scale structure of associations is not set by the expansion of clusters, rather it is a relic of the molecular gas cloud from which the association was formed. This finding is inconsistent with a monolithic model of association formation and instead favours a hierarchical model, in which OB associations form in-situ, following the fractal structure of the gas from which they form.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: A model to explain observed 1.1 and 0.85 millimeter dust continuum number counts

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $μ$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $μ$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $μ$m number counts is almost exhausted.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$. At all redshifts, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

preprint2020arXiv

The extended Gaia-PS1-SDSS (GPS1+) proper motion catalog

The GPS1 catalog was released in 2017. It delivered precise proper motions for around 350 million sources across three-fourths of the sky down to a magnitude of $r\sim20$\,mag. In this study, we present GPS1+ the extension GPS1 catalog down to $r\sim22.5$\,mag, based on {\it Gaia} DR2, PS1, SDSS and 2MASS astrometry. The GPS1+ totally provides proper motions for $\sim$400 million sources with a characteristic systematic error of less than 0.1\masyr. This catalog is divided into two sub-samples, i.e., the primary and secondary parts. The primary $\sim$264 million sources have either or both of the {\it Gaia} and SDSS astrometry, with a typical precision of 2.0-5.0 \masyr. In this part, $\sim$160 million sources have {\it Gaia} proper motions, we provide another new proper motion for each of them by building a Bayesian model. Relative to {\it Gaia}'s values, the precision is improved by $\sim$0.1\,dex on average at the faint end; $\sim$50 million sources are the objects whose proper motions are missing in {\it Gaia} DR2, we provide their proper motion with a precision of $\sim$4.5\masyr; the remaining $\sim$54 million faint sources are beyond {\it Gaia} detecting capability, we provide their proper motions for the first time with a precision of 7.0 \masyr. However, the secondary $\sim$136 million sources only have PS1 astrometry, the average precision is worse than 15.0 \masyr. All the proper motions have been validated using QSOs and the existing {\it Gaia} proper motions. The catalog will be released on-line and available via the VO-TAP Service, or via the National Astronomical Data Center serviced by China-VO: https://nadc.china-vo.org/data/data/gps1p/f.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

preprint2020arXiv

The Strength of the Dynamical Spiral Perturbation in the Galactic Disk

The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,φ)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $Σ_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $Ω_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$.

preprint2019arXiv

Abundance Estimates for 16 Elements in 6 Million Stars from LAMOST DR5 Low-Resolution Spectra

We present the determination of stellar parameters and individual elemental abundances for 6 million stars from $\sim$8 million low-resolution ($R\sim1800$) spectra from LAMOST DR5. This is based on a modeling approach that we dub $The$ $Data$--$Driven$ $Payne$ ($DD$--$Payne$), which inherits essential ingredients from both {\it The Payne} \citep{Ting2019} and $The$ $Cannon$ \citep{Ness2015}. It is a data-driven model that incorporates constraints from theoretical spectral models to ensure the derived abundance estimates are physically sensible. Stars in LAMOST DR5 that are in common with either GALAH DR2 or APOGEE DR14 are used to train a model that delivers stellar parameters ($T_{\rm eff}$, $\log g$, $V_{\rm mic}$) and abundances for 16 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Co, Ni, Cu, and Ba) when applied to LAMOST spectra. Cross-validation and repeat observations suggest that, for ${\rm S/N}_{\rm pix}\ge 50$, the typical internal abundance precision is 0.03--0.1\,dex for the majority of these elements, with 0.2--0.3\,dex for Cu and Ba, and the internal precision of $T_{\rm eff}$ and $\log g$ is better than 30\,K and 0.07\,dex, respectively. Abundance systematics at the $\sim$0.1\,dex level are present in these estimates, but are inherited from the high-resolution surveys' training labels. For some elements, GALAH provides more robust training labels, for others, APOGEE. We provide flags to guide the quality of the label determination and to identify binary/multiple stars in LAMOST DR5. The abundance catalogs are publicly accessible via \href{url}{http://dr5.lamost.org/doc/vac}.

preprint2019arXiv

The Separation Distribution of Ultra-Wide Binaries across Galactic Populations

We present an extensive sample of ultra-wide binary stars in the solar neighborhood, focusing on separations of $0.01<s/{pc}<1$. Using data from Gaia DR2, we define kinematic sub-populations via the systems' tangential velocities, i.e., disk-like ($v_{\perp, tot}<40$km/s), intermediate ($v_{\perp,tot}=(40-85)$km/s), and halo-like ($v_{\perp,tot}>85$km/s) samples, presuming that these velocity cuts represent a rough ordering in the binaries' age and metallicity. Through stringent cuts on astrometric precision, we can obtain pure binary samples with thousands of binaries in each sample. For all three populations, the distribution of binary separations extends smoothly to 1pc, displaying neither strong truncation nor bimodality. Fitting a smoothly-broken power law for the separation distribution, we find that its slope at separations $s=10^{2.5-4}$AU is the same for all sub-populations, $p(s)\propto s^γ$, with $γ\sim-1.54$. However, the logarithmic slope of $p(s)$ steepens at $s \gtrsim 10^4$AU. We find some evidences that the degree of steepening increases with the binaries' age, with a slope-change of only $Δγ\sim0.5$ for disk-like stars, but $Δγ>1$ for halo-like stars. This trend is contrary to what might be expected if steepening at wide separations were due to gravitational perturbations by molecular clouds or stars, which would preferentially disrupt disk binaries. If we were to interpret steepening at $s\gtrsim 10^4$AU as a consequence of disruption by massive halo objects, we would have to invoke a MACHO population inconsistent with other constraints. As a more plausible alternative, we propose a simple model to predict the separation distribution of wide binaries formed in dissolving star clusters. This model generically predicts $γ\sim-1.5$ as observed, with steepening at larger separations due to the finite size of binaries' birth clusters.

preprint2019arXiv

Tracing Kinematic and Chemical Properties of Sagittarius Stream by K-Giants, M-Giants, and BHB stars

We characterize the kinematic and chemical properties of $\sim$3,000 Sagittarius (Sgr) stream stars, including K-giants, M-giants, and BHBs, select from SEGUE-2, LAMOST, and SDSS separately in Integrals-of-Motion space. The orbit of Sgr stream is quite clear from the velocity vector in $X$-$Z$ plane. Stars traced by K-giants and M-giants present the apogalacticon of trailing steam is $\sim$ 100 kpc. The metallicity distributions of Sgr K-, M-giants, and BHBs present that the M-giants are on average the most metal-rich population, followed by K-giants and BHBs. All of the K-, M-giants, and BHBs indicate that the trailing arm is on average more metal-rich than leading arm, and the K-giants show that the Sgr debris is the most metal-poor part. The $α$-abundance of Sgr stars exhibits a similar trend with the Galactic halo stars at lower metallicity ([Fe/H] $<\sim$ $-$1.0 dex), and then evolve down to lower [$α$/Fe] than disk stars at higher metallicity, which is close to the evolution pattern of $α$-element of Milky Way dwarf galaxies. We find $V_Y$ and metallicity of K-giants have gradients along the direction of line-of-sight from the Galactic center in $X$-$Z$ plane, and the K-giants show that $V_Y$ increases with metallicity at [Fe/H] $>\sim-$1.5 dex. After dividing the Sgr stream into bright and faint stream according to their locations in equatorial coordinate, the K-giants and BHBs show that the bright and faint stream present different $V_Y$ and metallicities, the bright stream is on average higher in $V_Y$ and metallicity than the faint stream.

preprint2016arXiv

A radial age gradient in the geometrically thick disk of the Milky Way

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

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

Accelerated Fitting of Stellar Spectra

Stellar spectra are often modeled and fit by interpolating within a rectilinear grid of synthetic spectra to derive the stars' labels: stellar parameters and elemental abundances. However, the number of synthetic spectra needed for a rectilinear grid grows exponentially with the label space dimensions, precluding the simultaneous and self-consistent fitting of more than a few elemental abundances. Shortcuts such as fitting subsets of labels separately can introduce unknown systematics and do not produce correct error covariances in the derived labels. In this paper we present a new approach -- CHAT (Convex Hull Adaptive Tessellation) -- which includes several new ideas for inexpensively generating a sufficient stellar synthetic library, using linear algebra and the concept of an adaptive, data-driven grid. A convex hull approximates the region where the data lie in the label space. A variety of tests with mock datasets demonstrate that CHAT can reduce the number of required synthetic model calculations by three orders of magnitude in an 8D label space. The reduction will be even larger for higher-dimensional label spaces. In CHAT the computational effort increases only linearly with the number of labels that are fit simultaneously. Around each of these grid points in label space an approximate synthetic spectrum can be generated through linear expansion using a set of "gradient spectra" that represent flux derivatives at every wavelength point with respect to all labels. These techniques provide new opportunities to fit the full stellar spectra from large surveys with 15-30 labels simultaneously.

preprint2016arXiv

Action-based Dynamical Modelling for the Milky Way Disk

We present RoadMapping, a full-likelihood dynamical modelling machinery that aims to recover the Milky Way's (MW) gravitational potential from large samples of stars in the Galactic disk. RoadMapping models the observed positions and velocities of stars with a parametrized, three-integral distribution function (DF) in a parametrized axisymmetric potential. We investigate through differential test cases with idealized mock data how the breakdown of model assumptions and data properties affect constraints on the potential and DF. Our key results are: (i) If the MW's true potential is not included in the assumed model potential family, we can - in the axisymmetric case - still find a robust estimate for the potential, with only <~ 10% difference in surface density within |z| <= 1.1 kpc inside the observed volume. (ii) Modest systematic differences between the true and model DF are inconsequential. E.g, when binning stars to define sub-populations with simple DFs, binning errors do not affect the modelling as long as the DF parameters of neighbouring bins differ by < 20%. In addition, RoadMapping ensures unbiased potential estimates for either (iii) small misjudgements of the spatial selection function (i.e., <~ 15% at the survey volume's edge), (iv) if distances are known to within 10%, or (v) if proper motion uncertainties are known within 10% or are smaller than delta mu <~ 1 mas/yr. Challenges are the rapidly increasing computational costs for large sample sizes. Overall, RoadMapping is well suited to making precise new measurements of the MW's potential with data from the upcoming Gaia releases.

preprint2016arXiv

Ages of massive galaxies at $0.5 < z < 2.0$ from 3D-HST rest-frame optical spectroscopy

We present low-resolution near-infrared stacked spectra from the 3D-HST survey up to $z=2.0$ and fit them with commonly used stellar population synthesis models: BC03 (Bruzual & Charlot, 2003), FSPS10 (Flexible Stellar Population Synthesis, Conroy & Gunn 2010), and FSPS-C3K (Conroy, Kurucz, Cargile, Castelli, in prep). The accuracy of the grism redshifts allows the unambiguous detection of many emission and absorption features, and thus a first systematic exploration of the rest-frame optical spectra of galaxies up to $z=2$. We select massive galaxies ($\rm log(M_{*} / M_{\odot}) > 10.8$), we divide them into quiescent and star-forming via a rest-frame color-color technique, and we median-stack the samples in 3 redshift bins between $z=0.5$ and $z=2.0$. We find that stellar population models fit the observations well at wavelengths below $\rm 6500 Å$ rest-frame, but show systematic residuals at redder wavelengths. The FSPS-C3K model generally provides the best fits (evaluated with a $χ^2_{red}$ statistics) for quiescent galaxies, while BC03 performs the best for star-forming galaxies. The stellar ages of quiescent galaxies implied by the models, assuming solar metallicity, vary from 4 Gyr at $z \sim 0.75$ to 1.5 Gyr at $z \sim 1.75$, with an uncertainty of a factor of 2 caused by the unknown metallicity. On average the stellar ages are half the age of the Universe at these redshifts. We show that the inferred evolution of ages of quiescent galaxies is in agreement with fundamental plane measurements, assuming an 8 Gyr age for local galaxies. For star-forming galaxies the inferred ages depend strongly on the stellar population model and the shape of the assumed star-formation history.

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

Chemical tagging can work: Identification of stellar phase-space structures purely by chemical-abundance similarity

Chemical tagging promises to use detailed abundance measurements to identify spatially separated stars that were in fact born together (in the same molecular cloud), long ago. This idea has not yielded much practical success, presumably because of the noise and incompleteness in chemical-abundance measurements. We have succeeded in substantially improving spectroscopic measurements with The Cannon, which has now delivered 15 individual abundances for ~100,000 stars observed as part of the APOGEE spectroscopic survey, with precisions around 0.04 dex. We test the chemical-tagging hypothesis by looking at clusters in abundance space and confirming that they are clustered in phase space. We identify (by the k-means algorithm) overdensities of stars in the 15-dimensional chemical-abundance space delivered by The Cannon, and plot the associated stars in phase space. We use only abundance-space information (no positional information) to identify stellar groups. We find that clusters in abundance space are indeed clusters in phase space. We recover some known phase-space clusters and find other interesting structures. This is the first-ever project to identify phase-space structures at survey-scale by blind search purely in abundance space; it verifies the precision of the abundance measurements delivered by The Cannon; the prospects for future data sets appear very good.

preprint2016arXiv

Constructing Polynomial Spectral Models for Stars

Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these $\mathcal{N}\sim 10-40$ model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with $\mathcal{N}$. We suggest instead the construction of a polynomial spectral model (PSM) of order $\mathcal{O}$ for the model flux at each wavelength. Building this approximation requires a minimum of only ${\mathcal{N}+\mathcal{O}\choose\mathcal{O}}$ calculations: e.g. a quadratic spectral model ($\mathcal{O}=2$) to fit $\mathcal{N}=20$ labels simultaneously, can be constructed from as few as $231$ ab initio spectral model calculations; in practice, a somewhat larger number ($\sim 300-1000$) of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation only over a portion of label space, which will vary case by case. Yet, taking the APOGEE survey as an example, a single quadratic PSM provides a remarkably good approximation to the exact ab initio spectral models across much of this survey: for random labels within that survey the PSM approximates the flux to within $10^{-3}$, and recovers the abundances to within $\sim 0.02$ dex rms of the exact models. This enormous speed-up enables the simultaneous many-label fitting of spectra with computationally expensive ab initio models for stellar spectra, such as non-LTE models. A PSM also enables the simultaneous fitting of observational parameters, such as the spectrum's continuum or line-spread function.

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

On Galactic density modeling in the presence of dust extinction

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic mid-plane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially-complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.

preprint2016arXiv

The Cannon 2: A data-driven model of stellar spectra for detailed chemical abundance analyses

We have shown that data-driven models are effective for inferring physical attributes of stars (labels; Teff, logg, [M/H]) from spectra, even when the signal-to-noise ratio is low. Here we explore whether this is possible when the dimensionality of the label space is large (Teff, logg, and 15 abundances: C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) and the model is non-linear in its response to abundance and parameter changes. We adopt ideas from compressed sensing to limit overall model complexity while retaining model freedom. The model is trained with a set of 12,681 red-giant stars with high signal-to-noise spectroscopic observations and stellar parameters and abundances taken from the APOGEE Survey. We find that we can successfully train and use a model with 17 stellar labels. Validation shows that the model does a good job of inferring all 17 labels (typical abundance precision is 0.04 dex), even when we degrade the signal-to-noise by discarding ~50% of the observing time. The model dependencies make sense: the spectral derivatives with respect to abundances correlate with known atomic lines, and we identify elements belonging to atomic lines that were previously unknown. We recover (anti-)correlations in abundance labels for globular cluster stars, consistent with the literature. However we find the intrinsic spread in globular cluster abundances is 3--4 times smaller than previously reported. We deliver 17 labels with associated errors for 87,563 red giant stars, as well as open-source code to extend this work to other spectroscopic surveys.

preprint2016arXiv

The Gaia-ESO survey: Hydrogen lines in red giants directly trace stellar mass

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

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 stellar population structure of the Galactic disk

The spatial structure of stellar populations with different chemical abundances in the Milky Way contains a wealth of information on Galactic evolution over cosmic time. We use data on 14,699 red-clump stars from the APOGEE survey, covering 4 kpc <~ R <~ 15 kpc, to determine the structure of mono-abundance populations (MAPs)---stars in narrow bins in [a/Fe] and [Fe/H]---accounting for the complex effects of the APOGEE selection function and the spatially-variable dust obscuration. We determine that all MAPs with enhanced [a/Fe] are centrally concentrated and are well-described as exponentials with a scale length of 2.2+/-0.2 kpc over the whole radial range of the disk. We discover that the surface-density profiles of low-[a/Fe] MAPs are complex: they do not monotonically decrease outwards, but rather display a peak radius ranging from ~5 kpc to ~13 kpc at low [Fe/H]. The extensive radial coverage of the data allows us to measure radial trends in the thickness of each MAP. While high-[a/Fe] MAPs have constant scale heights, low-[a/Fe] MAPs flare. We confirm, now with high-precision abundances, previous results that each MAP contains only a single vertical scale height and that low-[Fe/H], low-[a/Fe] and high-[Fe/H], high-[a/Fe] MAPs have intermediate (h_Z~300 to 600 pc) scale heights that smoothly bridge the traditional thin- and thick-disk divide. That the high-[a/Fe], thick disk components do not flare is strong evidence against their thickness being caused by radial migration. The correspondence between the radial structure and chemical-enrichment age of stellar populations is clear confirmation of the inside-out growth of galactic disks. The details of these relations will constrain the variety of physical conditions under which stars form throughout the MW disk.

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

A Three-Dimensional Map of Milky-Way Dust

We present a three-dimensional map of interstellar dust reddening, covering three-quarters of the sky out to a distance of several kiloparsecs, based on Pan-STARRS 1 and 2MASS photometry. The map reveals a wealth of detailed structure, from filaments to large cloud complexes. The map has a hybrid angular resolution, with most of the map at an angular resolution of 3.4' to 13.7', and a maximum distance resolution of ~25%. The three-dimensional distribution of dust is determined in a fully probabilistic framework, yielding the uncertainty in the reddening distribution along each line of sight, as well as stellar distances, reddenings and classifications for 800 million stars detected by Pan-STARRS 1. We demonstrate the consistency of our reddening estimates with those of two-dimensional emission-based maps of dust reddening. In particular, we find agreement with the Planck 353 GHz optical depth-based reddening map to within 0.05 mag in E(B-V) to a depth of 0.5 mag, and explore systematics at reddenings less than E(B-V) ~ 0.08 mag. We validate our per-star reddening estimates by comparison with reddening estimates for stars with both SDSS photometry and SEGUE spectral classifications, finding per-star agreement to within 0.1 mag out to a stellar E(B-V) of 1 mag. We compare our map to two existing three-dimensional dust maps, by Marshall et al. (2006) and Lallement et al. (2013), demonstrating our finer angular resolution, and better distance resolution compared to the former within ~3 kpc. The map can be queried or downloaded at http://argonaut.skymaps.info. We expect the three-dimensional reddening map presented here to find a wide range of uses, among them correcting for reddening and extinction for objects embedded in the plane of the Galaxy, studies of Galactic structure, calibration of future emission-based dust maps and determining distances to objects of known reddening.

preprint2015arXiv

APOGEE chemical tagging constraint on the maximum star cluster mass in the $α$-enhanced Galactic disk

Stars born from the same molecular cloud should be nearly homogeneous in their element abundances. The concept of chemical tagging is to identify members of disrupted clusters by their clustering in element abundance space. Chemical tagging requires large samples of stars with precise abundances for many individual elements. With uncertainties of $σ_{[X/{\rm Fe}]}$ and $σ_{\rm [Fe/H]} \simeq 0.05$ for 10 elements measured for $>10^4$ stars, the APOGEE DR12 spectra may be the first well-suited data set to put this idea into practice. We find that even APOGEE data offer only $\sim 500$ independent volume elements in the 10-dimensional abundance space, when we focus on the $α$-enhanced Galactic disk. We develop and apply a new algorithm to search for chemically homogeneous sets of stars against a dominant background. By injecting star clusters into the APOGEE data set we show that chemically homogeneous clusters with masses $\gtrsim 3 \times 10^7 \, {\rm M}_\odot$ would be easily detectable and yet no such signal is seen in the data. By generalizing this approach, we put a first abundance-based constraint on the cluster mass function for the old disk stars in the Milky Way.

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

Galaxy Structure as a Driver of the Star Formation Sequence Slope and Scatter

It is well established that (1) star-forming galaxies follow a relation between their star formation rate (SFR) and stellar mass (M$_{\star}$), the "star-formation sequence", and (2) the SFRs of galaxies correlate with their structure, where star-forming galaxies are less concentrated than quiescent galaxies at fixed mass. Here, we consider whether the scatter and slope of the star-formation sequence is correlated with systematic variations in the Sersic indices, $n$, of galaxies across the SFR-M$_{\star}$ plane. We use a mass-complete sample of 23,848 galaxies at $0.5<z<2.5$ selected from the 3D-HST photometric catalogs. Galaxy light profiles parameterized by $n$ are based on Hubble Space Telescope CANDELS near-infrared imaging. We use a single SFR indicator empirically-calibrated from stacks of Spitzer/MIPS 24$μ$m imaging, adding the unobscured and obscured star formation. We find that the scatter of the star-formation sequence is related in part to galaxy structure; the scatter due to variations in $n$ at fixed mass for star-forming galaxies ranges from 0.14$\pm$0.02 dex at $z\sim2$ to 0.30$\pm$0.04 dex at $z<1$. While the slope of the log(SFR)-log(M$_{\star}$) relation is of order unity for disk-like galaxies, galaxies with $n>2$ (implying more dominant bulges) have significantly lower SFR/M$_{\star}$ than the main ridgeline of the star-formation sequence. These results suggest that bulges in massive $z\sim2$ galaxies are actively building up, where the stars in the central concentration are relatively young. At $z<1$, the presence of older bulges within star-forming galaxies lowers global SFR/M$_{\star}$, decreasing the slope and contributing significantly to the scatter of the star-formation sequence.

preprint2015arXiv

Mass estimates from stellar proper motions: The mass of $ω$ Centauri

A projected mass estimator and a maximum likelihood estimator are developed to best determine the mass of stellar clusters from proper motion data. Both methods effectively account and correct for errors in the velocities, and provide an unbiased and robust estimate of the cluster mass. Using an extensive proper motion study of $ω$ Centauri by van Leeuwen et al. (2000), we estimate the mass of $ω$ Centauri using these two different methods to be $2.8 \times 10^{6} M_{\odot} \left[d/5.1 kpc\right]^{3}$. Within this modeling context, the statistical error is $3\%$ and the systematic error is $4\%$.

preprint2015arXiv

Red giant masses and ages derived from carbon and nitrogen abundances

We show that the masses of red giant stars can be well predicted from their photospheric carbon and nitrogen abundances, in conjunction with their spectroscopic stellar labels log g, Teff, and [Fe/H]. This is qualitatively expected from mass-dependent post main sequence evolution. We here establish an empirical relation between these quantities by drawing on 1,475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from APOGEE DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional r.m.s. errors of about 14% (typically 0.2 Msun). From these masses, we derive ages with r.m.s errors of 40%. This empirical model allows us for the first time to make age determinations (in the range 1-13 Gyr) for vast numbers of giant stars across the Galaxy. We apply our model to 52,000 stars in APOGEE DR12, for which no direct mass and age information was previously available. We find that these estimates highlight the vertical age structure of the Milky Way disk, and that the relation of age with [alpha/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

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 Masses and SFRs for 1M Galaxies from SDSS and WISE

We combine SDSS and WISE photometry for the full SDSS spectroscopic galaxy sample, creating SEDs that cover lambda=0.4-22 micron for an unprecedented large and comprehensive sample of 858,365 present-epoch galaxies. Using MAGPHYS we then model simultaneously and consistently both the attenuated stellar SED and the dust emission at 12 micron and 22 micron, producing robust new calibrations for monochromatic mid-IR star formation rate proxies. These modeling results provide the first mid-IR-based view of the bi-modality in star formation activity among galaxies, exhibiting the sequence of star-forming galaxies (main sequence) with a slope of dlogSFR/dlogM*=0.80 and a scatter of 0.39 dex. We find that these new star-formation rates along the SF main sequence are systematically lower by a factor of 1.4 than those derived from optical spectroscopy. We show that for most present-day galaxies the 0.4-22 micron SED fits can exquisitely predict the fluxes measured by Herschel at much longer wavelengths. Our analysis also illustrates that the majorities of stars in the present-day universe is formed in luminous galaxies (~L*) in and around the green valley of the color-luminosity plane. We make the matched photometry catalog and SED modeling results publicly available.

preprint2015arXiv

The 3D-HST Survey: Hubble Space Telescope WFC3/G141 grism spectra, redshifts, and emission line measurements for $\sim 100,000$ galaxies

We present reduced data and data products from the 3D-HST survey, a 248-orbit HST Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper (Skelton et al. 2014) we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N. The data analysis is complicated by the fact that no slits are used: all objects in the WFC3 field are dispersed, and many spectra overlap. We developed software to automatically and optimally extract interlaced 2D and 1D spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has 98,663 measured redshifts and line strengths down to $JH_{IR}\leq 26$ and 22,548 with $JH_{IR}\leq 24$, where we comfortably detect continuum emission. Of this sample 5,459 galaxies are at $z>1.5$ and 9,621 are at $0.7<z<1.5$, where H$α$ falls in the G141 wavelength coverage. Based on comparisons with ground-based spectroscopic redshifts, and on analyses of paired galaxies and repeat observations, the typical redshift error for $JH_{IR}\leq 24$ galaxies in our catalog is $σ_z \approx 0.003 \times (1+z)$, i.e., one native WFC3 pixel. The $3σ$ limit for emission line fluxes of point sources is $1.5\times10^{-17}$ ergs s$^{-1}$ cm$^{-2}$. We show various representations of the full dataset, as well as individual examples that highlight the range of spectra that we find in the survey.

preprint2015arXiv

The Cannon: A data-driven approach to stellar label determination

New spectroscopic surveys offer the promise of consistent stellar parameters and abundances ('stellar labels') for hundreds of thousands of stars in the Milky Way: this poses a formidable spectral modeling challenge. In many cases, there is a sub-set of reference objects for which the stellar labels are known with high(er) fidelity. We take advantage of this with The Cannon, a new data-driven approach for determining stellar labels from spectroscopic data. The Cannon learns from the 'known' labels of reference stars how the continuum-normalized spectra depend on these labels by fitting a flexible model at each wavelength; then, The Cannon uses this model to derive labels for the remaining survey stars. We illustrate The Cannon by training the model on only 542 stars in 19 clusters as reference objects, with Teff, log g and [Fe/H] as the labels, and then applying it to the spectra of 56,000 stars from APOGEE DR10. The Cannon is very accurate. Its stellar labels compare well to the stars for which APOGEE pipeline (ASPCAP) labels are provided in DR10, with rms differences that are basically identical to the stated ASPCAP uncertainties. Beyond the reference labels, The Cannon makes no use of stellar models nor any line-list, but needs a set of reference objects that span label-space. The Cannon performs well at lower signal-to-noise, as it delivers comparably good labels even at one ninth the APOGEE observing time. We discuss the limitations of The Cannon and its future potential, particularly, to bring different spectroscopic surveys onto a consistent scale of stellar labels.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

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 Magellanic Inter-Cloud Project (MAGIC) II: Slicing up the Bridge

The origin of the gas in between the Magellanic Clouds (MCs), known as the Magellanic Bridge (MB), has always been the subject of controversy. To shed light into this, we present the results from the MAGIC II project aimed at probing the stellar populations in ten large fields located perpendicular to the main ridge-line of HI in the Inter-Cloud region. We secured these observations of the stellar populations in between the MCs using the WFI camera on the 2.2 m telescope in La Silla. Using colour-magnitude diagrams (CMDs), we trace stellar populations across the Inter-Cloud region. In good agreement with MAGIC I, we find significant intermediate-age stars in the Inter-Cloud region as well as young stars of a similar age to the last pericentre passage in between the MCs (~200 Myr ago). We show here that the young, intermediate-age and old stars have distinct spatial distributions. The young stars correlate well with the HI gas suggesting that they were either recently stripped from the SMC or formed in-situ. The bulk of intermediate-age stars are located mainly in the bridge region where the HI column density is higher, but they are more spread out than the young stars. They have very similar properties to stars located ~2 Kpc from the SMC centre, suggesting that they were tidally stripped from this region. Finally, the old stars extend to some 8 Kpc from the SMC supporting the idea that all galaxies have a large extended metal poor stellar halo.

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]

preprint2015arXiv

The radial profile and flattening of the Milky Way's stellar halo to $\rm 80~$kpc from the SEGUE K-giant Survey

We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $ν_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $α_{in}=2.1\pm 0.3$ and $α_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the "outer" stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.

preprint2015arXiv

The SEGUE K Giant Survey. III. Quantifying Galactic Halo Substructure

We statistically quantify the amount of substructure in the Milky Way stellar halo using a sample of 4568 halo K giant stars at Galactocentric distances ranging over 5-125 kpc. These stars have been selected photometrically and confirmed spectroscopically as K giants from the Sloan Digital Sky Survey's SEGUE project. Using a position-velocity clustering estimator (the 4distance) and a model of a smooth stellar halo, we quantify the amount of substructure in the halo, divided by distance and metallicity. Overall, we find that the halo as a whole is highly structured. We also confirm earlier work using BHB stars which showed that there is an increasing amount of substructure with increasing Galactocentric radius, and additionally find that the amount of substructure in the halo increases with increasing metallicity. Comparing to resampled BHB stars, we find that K giants and BHBs have similar amounts of substructure over equivalent ranges of Galactocentric radius. Using a friends-of-friends algorithm to identify members of individual groups, we find that a large fraction (~33%) of grouped stars are associated with Sgr, and identify stars belonging to other halo star streams: the Orphan Stream, the Cetus Polar Stream, and others, including previously unknown substructures. A large fraction of sample K giants (more than 50%) are not grouped into any substructure. We find also that the Sgr stream strongly dominates groups in the outer halo for all except the most metal-poor stars, and suggest that this is the source of the increase of substructure with Galactocentric radius and metallicity.

preprint2015arXiv

Where stars form: inside-out growth and coherent star formation from HST Halpha maps of 2676 galaxies across the main sequence at z~1

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

preprint2015arXiv

Young alpha-enriched giant stars in the solar neighbourhood

We derive age constraints for 1639 red giants in the APOKASC sample for which seismic parameters from Kepler, as well as effective temperatures, metallicities and [alpha/Fe] values from APOGEE DR12 are available. We investigate the relation between age and chemical abundances for these stars, using a simple and robust approach to obtain ages. We first derive stellar masses using standard seismic scaling relations, then determine the maximum possible age for each star as function of its mass and metallicity, independently of its evolutionary stage. While the overall trend between maximum age and chemical abundances is a declining fraction of young stars with increasing [alpha/Fe], at least 14 out of 241 stars with [alpha/Fe]>0.13 are younger than 6 Gyr. Five stars with [alpha/Fe]>0.2 have ages below 4 Gyr. We examine the effect of modifications in the standard seismic scaling relations, as well as the effect of very low helium fractions, but these changes are not enough to make these stars as old as usually expected for alpha-rich stars (i.e., ages greater than 8-9 Gyr). Such unusual alpha-rich young stars have also been detected by other surveys, but defy simple explanations in a galaxy evolution context.

preprint2014arXiv

3D-HST WFC3-selected Photometric Catalogs in the Five CANDELS/3D-HST Fields: Photometry, Photometric Redshifts and Stellar Masses

The 3D-HST and CANDELS programs have provided WFC3 and ACS spectroscopy and photometry over ~900 square arcminutes in five fields: AEGIS, COSMOS, GOODS-North, GOODS-South, and the UKIDSS UDS field. All these fields have a wealth of publicly available imaging datasets in addition to the HST data, which makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range. In this paper we describe a photometric analysis of the CANDELS and 3D-HST HST imaging and the ancillary imaging data at wavelengths 0.3um to 8um. Objects were selected in the WFC3 near-IR bands, and their SEDs were determined by carefully taking the effects of the point spread function in each observation into account. A total of 147 distinct imaging datasets were used in the analysis. The photometry is made available in the form of six catalogs: one for each field, as well as a master catalog containing all objects in the entire survey. We also provide derived data products: photometric redshifts, determined with the EAZY code, and stellar population parameters determined with the FAST code. We make all the imaging data that were used in the analysis available, including our reductions of the WFC3 imaging in all five fields. 3D-HST is a spectroscopic survey with the WFC3 and ACS grisms, and the photometric catalogs presented here constitute a necessary first step in the analysis of these grism data. All the data presented in this paper are available through the 3D-HST website.

preprint2014arXiv

A Massive Galaxy in its Core Formation Phase Three Billion Years After the Big Bang

Most massive galaxies are thought to have formed their dense stellar cores at early cosmic epochs. However, cores in their formation phase have not yet been observed. Previous studies have found galaxies with high gas velocity dispersions or small apparent sizes but so far no objects have been identified with both the stellar structure and the gas dynamics of a forming core. Here we present a candidate core in formation 11 billion years ago, at z=2.3. GOODS-N-774 has a stellar mass of 1.0x10^11 Msun, a half-light radius of 1.0 kpc, and a star formation rate of 90[+45-20]Msun/yr. The star forming gas has a velocity dispersion 317+-30 km/s, amongst the highest ever measured. It is similar to the stellar velocity dispersions of the putative descendants of GOODS-N-774, compact quiescent galaxies at z~2 and giant elliptical galaxies in the nearby Universe. Galaxies such as GOODS-N-774 appear to be rare; however, from the star formation rate and size of the galaxy we infer that many star forming cores may be heavily obscured, and could be missed in optical and near-infrared surveys.

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

Andromeda's Dust

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

preprint2014arXiv

Estimating Black Hole Masses in Hundreds of Quasars

We explore the practical feasibility of AGN broad-band reverberation mapping and present first results. We lay out and apply a rigorous approach for stochastic reverberation mapping of unevenly sampled multi-broad-band flux measurements, assuming that the broad-line region (BLR) line flux is contributing up to 15 % in some bands, and is directly constrained by one spectroscopical epoch. The approach describes variations of the observed flux as the continuum, modeled as a stochastic Gaussian process, and emission line contribution, modeled as a scaled, smoothed and delayed version of the continuum. This approach is capable not only to interpolate in time between measurements, but also to determine confidence limits on continuum -- line emission delays. This approach is applied to SDSS observations in 'Stripe 82' (S82) providing flux measurements precise to 2 % at $\sim$ 60 epochs over $\sim$ 10 years. The strong annual variations in the epoch sampling prove a serious limitation in practice. Also, suitable redshift ranges must be identified, where strong broad emission line contribute to one filter, but not to another. Through generating and evaluating problem-specific mock data, we verify that S82-like data can constrain $τ_{\mathrm{delay}}$ for a simple transfer function model. In application to real data, we estimate $τ_{\mathrm{delay}}$ for 323 AGN with $0.225 < z < 0.846$, combining information for different objects through the ensemble-scaling relationships for BLR size and BH mass. Our analysis tentatively indicates a 1.7 times larger BLR size of H$α$ and MgII compared to Kaspi2000 and Vestergaard2002, but the seasonal data sampling casts doubt on the robustness of the inference.

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

How dead are dead galaxies? Mid-Infrared fluxes of quiescent galaxies at redshift 0.3 < z < 2.5: implications for star formation rates and dust heating

We investigate the star formation rates of quiescent galaxies at high redshift (0.3 < z < 2.5) using 3D-HST WFC3 grism spectroscopy and Spitzer mid-infrared data. We select quiescent galaxies on the basis of the widely used UVJ color-color criteria. Spectral energy distribution fitting (rest frame optical and near-IR) indicates very low star formation rates for quiescent galaxies (sSFR ~ 10^-12 yr^-1). However, SED fitting can miss star formation if it is hidden behind high dust obscuration and ionizing radiation is re-emitted in the mid-infrared. It is therefore fundamental to measure the dust-obscured SFRs with a mid-IR indicator. We stack the MIPS-24um images of quiescent objects in five redshift bins centered on z = 0.5, 0.9, 1.2, 1.7, 2.2 and perform aperture photometry. Including direct 24um detections, we find sSFR ~ 10^-11.9 * (1+z)^4 yr^-1. These values are higher than those indicated by SED fitting, but at each redshift they are 20-40 times lower than those of typical star forming galaxies. The true SFRs of quiescent galaxies might be even lower, as we show that the mid-IR fluxes can be due to processes unrelated to ongoing star formation, such as cirrus dust heated by old stellar populations and circumstellar dust. Our measurements show that star formation quenching is very efficient at every redshift. The measured SFR values are at z > 1.5 marginally consistent with the ones expected from gas recycling (assuming that mass loss from evolved stars refuels star formation) and well above that at lower redshifts.

preprint2014arXiv

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

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

preprint2014arXiv

On the importance of using appropriate spectral models to derive physical properties of galaxies at 0.7<z<2.8

Interpreting observations of distant galaxies in terms of constraints on physical parameters - such as stellar mass, star-formation rate (SFR) and dust optical depth - requires spectral synthesis modelling. We analyse the reliability of these physical parameters as determined under commonly adopted `classical' assumptions: star-formation histories assumed to be exponentially declining functions of time, a simple dust law and no emission-line contribution. Improved modelling techniques and data quality now allow us to use a more sophisticated approach, including realistic star-formation histories, combined with modern prescriptions for dust attenuation and nebular emission (Pacifici et al. 2012). We present a Bayesian analysis of the spectra and multi-wavelength photometry of 1048 galaxies from the 3D-HST survey in the redshift range 0.7<z<2.8 and in the stellar mass range 9<log(M/Mo)<12. We find that, using the classical spectral library, stellar masses are systematically overestimated (~0.1 dex) and SFRs are systematically underestimated (~0.6 dex) relative to our more sophisticated approach. We also find that the simultaneous fit of photometric fluxes and emission-line equivalent widths helps break a degeneracy between SFR and optical depth of the dust, reducing the uncertainties on these parameters. Finally, we show how the biases of classical approaches can affect the correlation between stellar mass and SFR for star-forming galaxies (the `Star-Formation Main Sequence'). We conclude that the normalization, slope and scatter of this relation strongly depend on the adopted approach and demonstrate that the classical, oversimplified approach cannot recover the true distribution of stellar mass and SFR.

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

The APOGEE red-clump catalog: Precise distances, velocities, and high-resolution elemental abundances over a large area of the Milky Way's disk

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar-evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5 to 10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc^3 over 5 kpc <~ R <~ 14 kpc, |Z| <~ 2 kpc, and -15 deg <~ Galactocentric azimuth <~ 30 deg. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from 2MASS, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ~45 degree-wide azimuthal region covered by the current sample to be <= 0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kpc from the Sun.

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

Tracing chemical evolution over the extent of the Milky Way's Disk with APOGEE Red Clump Stars

We employ the first two years of data from the near-infrared, high-resolution SDSS-III/APOGEE spectroscopic survey to investigate the distribution of metallicity and alpha-element abundances of stars over a large part of the Milky Way disk. Using a sample of ~10,000 kinematically-unbiased red-clump stars with ~5% distance accuracy as tracers, the [alpha/Fe] vs. [Fe/H] distribution of this sample exhibits a bimodality in [alpha/Fe] at intermediate metallicities, -0.9<[Fe/H]<-0.2, but at higher metallicities ([Fe/H]=+0.2) the two sequences smoothly merge. We investigate the effects of the APOGEE selection function and volume filling fraction and find that these have little qualitative impact on the alpha-element abundance patterns. The described abundance pattern is found throughout the range 5<R<11 kpc and 0<|Z|<2 kpc across the Galaxy. The [alpha/Fe] trend of the high-alpha sequence is surprisingly constant throughout the Galaxy, with little variation from region to region (~10%). Using simple galactic chemical evolution models we derive an average star formation efficiency (SFE) in the high-alpha sequence of ~4.5E-10 1/yr, which is quite close to the nearly-constant value found in molecular-gas-dominated regions of nearby spirals. This result suggests that the early evolution of the Milky Way disk was characterized by stars that shared a similar star formation history and were formed in a well-mixed, turbulent, and molecular-dominated ISM with a gas consumption timescale (1/SFE) of ~2 Gyr. Finally, while the two alpha-element sequences in the inner Galaxy can be explained by a single chemical evolutionary track this cannot hold in the outer Galaxy, requiring instead a mix of two or more populations with distinct enrichment histories.

preprint2013arXiv

A direct dynamical measurement of the Milky Way's disk surface density profile, disk scale length, and dark matter profile at 4 kpc < R < 9 kpc

We present and apply rigorous dynamical modeling with which we infer unprecedented constraints on the stellar and dark matter mass distribution within our Milky Way (MW), based on large sets of phase-space data on individual stars. Specifically, we model the dynamics of 16,269 G-type dwarfs from SEGUE, which sample 5 < R_GC/kpc < 12 and 0.3 < |Z|/kpc < 3. We independently fit a parameterized MW potential and a three-integral, action-based distribution function (DF) to the phase-space data of 43 separate abundance-selected sub-populations (MAPs), accounting for the complex selection effects affecting the data. We robustly measure the total surface density within 1.1 kpc of the mid-plane to 5% over 4.5 < R_GC/kpc < 9. Using metal-poor MAPs with small radial scale lengths as dynamical tracers probes 4.5 < R_GC/kpc < 7, while MAPs with longer radial scale lengths sample 7 < R_GC/kpc < 9. We measure the mass-weighted Galactic disk scale length to be R_d = 2.15+/-0.14 kpc, in agreement with the photometrically inferred spatial distribution of stellar mass. We thereby measure dynamically the mass of the Galactic stellar disk to unprecedented accuracy: M_* = 4.6+/-0.3+3.0x(R_0/kpc-8)x10^{10}Msun and a total local surface density of Σ_{R_0}(Z=1.1 kpc) = 68+/-4 Msun/pc^2 of which 38+/-4 Msun/pc^2 is contributed by stars and stellar remnants. By combining our surface density measurements with the terminal velocity curve, we find that the MW's disk is maximal in that V_{c,disk} / V_{c,total} = 0.83+/-0.04 at R=2.2 R_d. We also constrain for the first time the radial profile of the dark halo at such small Galactocentric radii, finding that ρ_{DM} (r;near R_0) \propto 1 / r^αwith α< 1.53 at 95% confidence. Our results show that action-based distribution-function modeling of complex stellar data sets is now a feasible approach that will be fruitful for interpreting Gaia data.

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

Constraining the Galactic potential via action-based distribution functions for mono-abundance stellar populations

We present a rigorous and practical way of constraining the Galactic potential based on the phase-space information for many individual stars. Such an approach is needed to dynamically model the data from ongoing spectroscopic surveys of the Galaxy and in the future Gaia. This approach describes the orbit distribution of stars by a family of parametrized distribution function (DF) proposed by McMillan and Binney, which are based on actions. We find that these parametrized DFs are flexible enough to capture well the observed phase-space distributions of individual abundance-selected Galactic subpopulations of stars (`mono-abundance populations') for a disc-like gravitational potential, which enables independent dynamical constraints from each of the Galactic mono-abundance populations. We lay out a statistically rigorous way to constrain the Galactic potential parameters by constructing the joint likelihood of potential and DF parameters, and subsequently marginalizing over the DF parameters. This approach explicitly incorporates the spatial selection function inherent to all Galactic surveys, and can account for the uncertainties of the individual position--velocity observations. On that basis, we study the precision of the parameters of the Galactic potential that can be reached with various sample sizes and realistic spatial selection functions. By creating mock samples from the DF, we show that, even under a restrictive and realistic spatial selection function, given a two-parameter gravitational potential, one can recover the true potential parameters to a few per cent with sample sizes of a few thousands. The assumptions of axisymmetry, of DFs that are smooth in the actions and of no time variation remain important limitations in our current study.

preprint2013arXiv

Empirical predictions for (sub-)millimeter line and continuum deep fields

[abridged] Modern (sub-)millimeter/radio interferometers will enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will provide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission. In this paper, we present empirical predictions for such (sub-)millimeter line and continuum deep fields. We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble Ultra Deep Field. Using a physically-motivated spectral energy distribution model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z=5, and obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminosity. We derive statistical constraints on the dust emission in the infrared and (sub-)millimeter which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA and JVLA bands. Using empirical relations between the observed CO/[CII] line luminosities and the infrared luminosity, we infer the flux of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our sample. We then predict the fluxes of higher CO transition lines CO(2-1) to CO(7-6) bracketing two extreme gas excitation scenarios. We use our predictions to discuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct benchmark for future deep field campaigns in the (sub-)millimeter regime.

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

Molecular Gas and Star Formation in Nearby Disk Galaxies

We compare molecular gas traced by 12CO(2-1) maps from the HERACLES survey, with tracers of the recent star formation rate (SFR) across 30 nearby disk galaxies. We demonstrate a first-order linear correspondence between Sig_mol and Sig_SFR but also find important second-order systematic variations in the apparent molecular gas depletion time, t_dep^mol = Sig_mol / Sig_SFR. At our 1 kpc common resolution, CO correlates closely with many tracers of the recent SFR. Weighting each line of sight equally and using a fixed, Milky Way alpha_CO, our data yield a molecular gas depletion time, t_dep^mol=Sig_mol/Sig_SFR ~ 2.2 Gyr with 0.3 dex scatter, in good agreement with literature data. We apply a forward-modeling approach to constrain the power-law index, N, that relates the SFR surface density and the molecular gas surface density and find N=1+/-0.15 for our full data set with some variation from galaxy to galaxy. However, we caution that a power law treatment oversimplifies the topic given that we observe correlations between t_dep^mol and other local and global quantities. The strongest of these are a decreased t_dep^mol in low-mass, low-metallicity galaxies and a correlation of the kpc-scale t_dep^mol with dust-to-gas ratio, D/G. These correlations can be explained by a CO-to-H2 conversion factor that depends on D/G in the theoretically expected way. This is not a unique interpretation, but external evidence of conversion factor variations makes it a conservative one. After applying a D/G-dependent alpha_CO, some weak correlations between t_dep^mol and local conditions persist. In particular, we observe lower t_dep^mol and enhanced CO excitation associated with some nuclear gas concentrations. These appear to reflect real enhancements in the SFR/H2 and t_dep appears multivalued at fixed Sig_mol, supporting the the idea of "disk" and "starburst" modes driven by environmental factors.

preprint2013arXiv

On the effect of the Cosmic Microwave Background in high-redshift (sub-)millimeter observations

Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.

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

Quiescent Galaxies in the 3D-HST Survey: Spectroscopic Confirmation of a Large Number of Galaxies with Relatively Old Stellar Populations at z~2

Quiescent galaxies at z~2 have been identified in large numbers based on rest-frame colors, but only a small number of these galaxies have been spectroscopically confirmed to show that their rest-frame optical spectra show either strong Balmer or metal absorption lines. Here, we median stack the rest-frame optical spectra for 171 photometrically-quiescent galaxies at 1.4 < z < 2.2 from the 3D-HST grism survey. In addition to Hbeta (4861A), we unambiguously identify metal absorption lines in the stacked spectrum, including the G-band (4304A), Mg I (5175A), and Na I (5894A). This finding demonstrates that galaxies with relatively old stellar populations already existed when the universe was ~3 Gyr old, and that rest-frame color selection techniques can efficiently select them. We find an average age of 1.3^0.1_0.3 Gyr when fitting a simple stellar population to the entire stack. We confirm our previous result from medium-band photometry that the stellar age varies with the colors of quiescent galaxies: the reddest 80% of galaxies are dominated by metal lines and have a relatively old mean age of 1.6^0.5_0.4 Gyr, whereas the bluest (and brightest) galaxies have strong Balmer lines and a spectroscopic age of 0.9^0.2_0.1 Gyr. Although the spectrum is dominated by an evolved stellar population, we also find [OIII] and Hbeta emission. Interestingly, this emission is more centrally concentrated than the continuum with L_[OIII] = 1.7 +/- 0.3 x 10^40 erg s^-1, indicating residual central star formation or nuclear activity.

preprint2013arXiv

Stellar Kinematics of z~2 Galaxies and the Inside-Out Growth of Quiescent Galaxies

Using stellar kinematics measurements, we investigate the growth of massive, quiescent galaxies from z~2 to today. We present X-Shooter spectra from the UV to NIR and dynamical mass measurements of 5 quiescent massive (>10^11 Msun) galaxies at z~2. This triples the sample of z>1.5 galaxies with well constrained (dsigma <100 km/s) velocity dispersion measurements. From spectral population synthesis modeling we find that these galaxies have stellar ages that range from 0.5-2 Gyr, with no signs of ongoing star formation. We measure velocity dispersions (290-450 km/s) from absorption lines and find that they are 1.6-2.1 times higher than those of galaxies in the SDSS at the same mass. Sizes are measured using GALFIT from HST-WFC3 H160 and UDS K-band images. The dynamical masses correspond well to the SED based stellar masses, with dynamical masses that are ~15% higher. We find that M_*/M_dyn may decrease slightly with time, which could reflect the increase of the dark matter fraction within an increasing effective radius. We combine different stellar kinematic studies from the literature, and examine the structural evolution from z~2 to z~0: we confirm that at fixed dynamical mass, the effective radius increases by a factor of ~2.8, and the velocity dispersion decreases by a factor of ~1.7. The mass density within one effective radius decreases by a factor of ~20, while within a fixed physical radius (1 kpc) it decreases only mildly (~2). When we allow for an evolving mass limit by selecting a population of galaxies at fixed number density, a stronger size growth with time is found (~4), velocity dispersion decreases by a factor of ~1.4, and interestingly, the mass density within 1 kpc is consistent with no evolution. This finding suggests that massive quiescent galaxies at ~2 grow inside-out, consistent with the expectations from minor mergers.

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

The Assembly of Milky Way-like Galaxies Since z~2.5

Galaxies with the mass of the Milky Way dominate the stellar mass density of the Universe but it is uncertain how and when they were assembled. Here we study progenitors of these galaxies out to z=2.5, using data from the 3D-HST and CANDELS Treasury surveys. We find that galaxies with present-day stellar masses of log(M)~10.7 built ~90% of their stellar mass since z=2.5, with most of the star formation occurring before z=1. In marked contrast to the assembly history of massive elliptical galaxies, mass growth is not limited to large radii: the mass in the central 2 kpc of the galaxies increased by a factor of 3.2+-0.8 between z=2.5 and z=1. We therefore rule out simple models in which bulges were fully assembled at high redshift and disks gradually formed around them. Instead, bulges (and black holes) likely formed in lockstep with disks, through bar instabilities, migration, or other processes. We find that after z=1 the growth in the central regions gradually stopped and the disk continued to build, consistent with recent studies of the gas distributions in z~1 galaxies and the properties of many spiral galaxies today.

preprint2013arXiv

The Gravitational Potential Near the Sun From SEGUE K-dwarf Kinematics

To constrain the Galactic gravitational potential near the Sun ($\sim$1.5 kpc), we derive and model the spatial and velocity distribution for a sample of 9000 K-dwarfs that have spectra from SDSS/SEGUE, which yield radial velocities and abundances ([Fe/H] & [$α$/Fe]). We first derive the spatial density distribution for stars of three abundance-selected sub-populations by accounting for the survey's selection function. The vertical profile of these sub-populations are simple exponentials and their vertical dispersion profile is nearly isothermal. To model these data, we apply the `vertical' Jeans Equation, which relates the observable tracer number density and vertical velocity dispersion to the gravitational potential or vertical force. We explore a number of functional forms for the vertical force law, and fit the dispersion and density profiles of all abundance selected sub-populations simultaneously in the same potential, and explore all parameter co-variances using MCMC. Our fits constrain a disk {\it mass} scale height $\lesssim$ 300 pc and the total surface mass density to be $67 \pm 6 M_{\odot} {\rm pc^{-2}}$ at $|z| = 1.0$ kpc of which the contribution from all stars is $42 \pm 5 M_{\odot} {\rm pc^{-2}}$ (presuming a contribution from cold gas of $13 M_{\odot} {\rm pc^{-2}}$). We find significant constraints on the local dark matter density of $0.0065\pm0.0023 M_{\odot} {\rm pc^{-3}}$ ($0.25\pm0.09 {\rm GeV cm^{-3}} $). Together with recent experiments this firms up the best estimate of $0.0075\pm0.0021 M_{\odot} {\rm pc^{-3}}$ ($0.28\pm0.08 {\rm GeV cm^{-3}} $), consistent with global fits of approximately round dark matter halos to kinematic data in the outskirts of the Galaxy.

preprint2013arXiv

The MAGellanic Inter-Cloud project (MAGIC) I: Evidence for intermediate-age stellar populations in between the Magellanic Clouds

The origin of the gas in between the Magellanic Clouds (MCs) -- known as the `Magellanic Bridge' (MB) -- is puzzling. Numerical simulations suggest that the MB formed from tidally stripped gas and stars in a recent interaction between the MCs. However, the apparent lack of stripped intermediate- or old-age stars associated with the MB is at odds with this picture. In this paper, we present the first results from the MAGellanic Inter-Cloud program (MAGIC) aimed at probing the stellar populations in the inter-cloud region. We present observations of the stellar populations in two large fields located in between the Large and Small Magellanic Clouds (LMC/SMC), secured using the WFI camera on the 2.2 m telescope in La Silla. Using a synthetic color-magnitude diagram (CMD) technique, we present the first quantitative evidence for the presence of intermediate-age and old stars in the inter-Cloud region. The intermediate-age stars -- that make up ~28% of all stars in the region -- are not present in fields at a similar distance from the SMC in a direction pointing away from the LMC. This provides potential evidence that these intermediate-age stars could have been tidally stripped from the SMC. However, spectroscopic studies will be needed to confirm or rule out the tidal origin for the inter-Cloud gas and stars.

preprint2013arXiv

The Milky Way's Stellar Disk

A suite of vast stellar surveys mapping the Milky Way, culminating in the Gaia mission, is revolutionizing the empirical information about the distribution and properties of stars in the Galactic stellar disk. We review and lay out what analysis and modeling machinery needs to be in place to test mechanisms of disk galaxy evolution and to stringently constrain the Galactic gravitational potential, using such Galactic star-by-star measurements. We stress the crucial role of stellar survey selection functions in any such modeling; and we advocate the utility of viewing the Galactic stellar disk as made up from `mono-abundance populations' (MAPs), both for dynamical modeling and for constraining the Milky Way's evolutionary processes. We review recent work on the spatial and kinematical distribution of MAPs, and lay out how further study of MAPs in the Gaia era should lead to a decisively clearer picture of the Milky Way's dark matter distribution and formation history.

preprint2013arXiv

The Radial Distribution of Star Formation in Galaxies at z~1 from the 3D-HST Survey

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic time. Thanks to the WFC3 grism on HST it is now possible to measure this beyond the local Universe. Here we present the spatial distribution of Halpha emission for a sample of 54 strongly star-forming galaxies at z~1 in the 3D-HST Treasury survey. By stacking the Halpha emission we find that star formation occurred in approximately exponential distributions at z~1, with median Sersic index of n=1.0+-0.2. The stacks are elongated with median axis ratios of b/a=0.58+-0.09 in Halpha, consistent with (possibly thick) disks at random orientation angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with inclination-corrected velocities of 90 to 330 km/s. The most straightforward interpretation of our results is that star formation in strongly star-forming galaxies at z~1 generally occurred in disks. The disks appear to be "scaled-up" versions of nearby spiral galaxies: they have EW(Halpha)~100 Angstroms out to the solar orbit and they have star formation surface densities above the threshold for driving galactic scale winds.

preprint2013arXiv

The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at z~1.5

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

preprint2013arXiv

The Structural Evolution of Milky Way-like Star Forming Galaxies since z~1.3

We follow the structural evolution of star forming galaxies (SFGs) like the Milky Way by selecting progenitors to z~1.3 based on the stellar mass growth inferred from the evolution of the star forming sequence. We select our sample from the 3D-HST survey, which utilizes spectroscopy from the HST WFC3 G141 near-IR grism and enables precise redshift measurements for our sample of SFGs. Structural properties are obtained from Sersic profile fits to CANDELS WFC3 imaging. The progenitors of z=0 SFGs with stellar mass M=10^{10.5} Msun are typically half as massive at z~1. This late-time stellar mass assembly is consistent with recent studies that employ abundance matching techniques. The descendant SFGs at z~0 have grown in half-light radius by a factor of ~1.4 since z~1. The half-light radius grows with stellar mass as r_e M^{0.29}. While most of the stellar mass is clearly assembling at large radii, the mass surface density profiles reveal ongoing mass growth also in the central regions where bulges and pseudobulges are common features in present day late-type galaxies. Some portion of this growth in the central regions is due to star formation as recent observations of H-alpha maps for SFGs at z~1 are found to be extended but centrally peaked. Connecting our lookback study with galactic archeology, we find the stellar mass surface density at R=8 kpc to have increased by a factor of ~2 since z~1, in good agreement with measurements derived for the solar neighborhood of the Milky Way.

preprint2012arXiv

3D-HST Grism Spectroscopy of a Gravitationally Lensed, Low-metallicity Starburst Galaxy at z=1.847

We present Hubble Space Telescope (HST) imaging and spectroscopy of the gravitational lens SL2SJ02176-0513, a cusp arc at z=1.847. The UV continuum of the lensed galaxy is very blue, which is seemingly at odds with its redder optical colors. The 3D-HST WFC3/G141 near-infrared spectrum of the lens reveals the source of this discrepancy to be extremely strong [OIII]5007 and H-beta emission lines with rest-frame equivalent widths of 2000 +/- 100 and 520 +/- 40 Angstroms, respectively. The source has a stellar mass ~10^8 Msun, sSFR \sim 100/Gyr, and detection of [OIII]4363 yields a metallicity of 12 + log(O/H) = 7.5 +/- 0.2. We identify local blue compact dwarf analogs to SL2SJ02176-0513, which are among the most metal-poor galaxies in the SDSS. The local analogs resemble the lensed galaxy in many ways, including UV/optical SED, spatial morphology and emission line equivalent widths and ratios. Common to SL2SJ02176-0513 and its local counterparts is an upturn at mid-IR wavelengths likely arising from hot dust heated by starbursts. The emission lines of SL2SJ02176-0513 are spatially resolved owing to the combination of the lens and the high spatial resolution of HST. The lensed galaxy is composed of two clumps with combined size r_e \sim 300 pc, and we resolve significant differences in UV color and emission line equivalent width between them. Though it has characteristics occasionally attributed to active galactic nuclei, we conclude that SL2SJ02176-0513 is a low-metallicity star-bursting dwarf galaxy. Such galaxies will be found in significant numbers in the full 3D-HST grism survey.

preprint2012arXiv

3D-HST: A wide-field grism spectroscopic survey with the Hubble Space Telescope

We present 3D-HST, a near-infrared spectroscopic Treasury program with the Hubble Space Telescope for studying the processes that shape galaxies in the distant Universe. 3D-HST provides rest-frame optical spectra for a sample of ~7000 galaxies at 1<z<3.5, the epoch when 60% of all star formation took place, the number density of quasars peaked, the first galaxies stopped forming stars, and the structural regularity that we see in galaxies today must have emerged. 3D-HST will cover 3/4 (625 sq.arcmin) of the CANDELS survey area with two orbits of primary WFC3/G141 grism coverage and two to four parallel orbits with the ACS/G800L grism. In the IR these exposure times yield a continuum signal-to-noise of ~5 per resolution element at H~23.1 and a 5sigma emission line sensitivity of 5x10-17 erg/s/cm2 for typical objects, improving by a factor of ~2 for compact sources in images with low sky background levels. The WFC3/G141 spectra provide continuous wavelength coverage from 1.1-1.6 um at a spatial resolution of ~0."13, which, combined with their depth, makes them a unique resource for studying galaxy evolution. We present the preliminary reduction and analysis of the grism observations, including emission line and redshift measurements from combined fits to the extracted grism spectra and photometry from ancillary multi-wavelength catalogs. The present analysis yields redshift estimates with a precision of sigma(z)=0.0034(1+z), or sigma(v)~1000 km/s. We illustrate how the generalized nature of the survey yields near-infrared spectra of remarkable quality for many different types of objects, including a quasar at z=4.7, quiescent galaxies at z~2, and the most distant T-type brown dwarf star known. The CANDELS and 3D-HST surveys combined will provide the definitive imaging and spectroscopic dataset for studies of the 1<z<3.5 Universe until the launch of the James Webb Space Telescope.

preprint2012arXiv

Chandra Observations of the Nuclear Star Cluster and Ultraluminous X-ray Sources in NGC 2139

We report Chandra observations of the Scd galaxy NGC 2139, which is known to host a recently formed (10^7.6 yrs) nuclear star cluster. The star cluster is undetected in X-rays, with an upper bound on 0.5-7 keV luminosity of L_X < 7.1 x 10^37 erg/s. This bound implies a bolometric accretion luminosity <0.3 percent of the Eddington luminosity for a black hole with the mass (approximately 3400 M_sun) expected from extrapolation of the M-sigma relation. The lack of X-ray emission indicates that a black hole, if present, is not undergoing significant accretion at the current time. While the central cluster is undetected, the data reveal a substantial population of bright X-ray point sources elsewhere in this galaxy, with eight qualifying as ultraluminous X-ray sources with L_X > 10^39 erg/s. We use archival Hubble Space Telescope images to identify candidate optical counterparts for seven Chandra sources, which in most cases have optical luminosities and spatial profiles consistent with star clusters. Compared with other galaxies, the number of luminous X-ray sources in NGC 2139 is larger by a factor of 4 - 10 than expected based on its present star formation rate and stellar mass. This finding can be understood if NGC 2139 has concluded a burst of star formation in the recent past, and suggests that this galaxy could be important for testing the use of X-ray source populations as a chronometer of star formation history.

preprint2012arXiv

Designing Imaging Surveys for a Retrospective Relative Photometric Calibration

In this paper, we investigate the impact of survey strategy on the performance of self-calibration when the goal is to produce accurate photometric catalogs from wide-field imaging surveys. This self-calibration technique utilizes multiple measurements of sources at different focal-plane positions to constrain instruments' large-scale response (flat-field) from survey science data alone. We create an artificial sky of sources and synthetically observe it under four basic survey strategies, creating an end-to-end simulation of an imaging survey for each. These catalog-level simulations include realistic measurement uncertainties and a complex focal-plane dependence of the instrument response. In the self-calibration step, we simultaneously fit for all the star fluxes and the parameters of a position-dependent flat-field. For realism, we deliberately fit with a wrong noise model and a flat-field functional basis that does not include the model that generated the synthetic data. We demonstrate that with a favorable survey strategy, a complex instrument response can be precisely self-calibrated. We show that returning the same sources to very different focal-plane positions is the key property of any survey strategy designed for accurate retrospective calibration of this type. The results of this work suggest the following advice for those considering the design of large-scale imaging surveys: Do not use a regular, repeated tiling of the sky; instead return the same sources to very different focal-plane positions.

preprint2012arXiv

Dwarfs Gobbling Dwarfs: A Stellar Tidal Stream Around NGC 4449 and Hierarchical Galaxy Formation on Small Scales

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf starburst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique combination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between luminosity or stellar-mass between the two galaxies is ~1:50, while the indirectly measured dynamical mass-ratio, when including dark matter, may be ~1:10-1:5. This system may thus represent a "stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and possibly in triggering their starbursts.

preprint2012arXiv

Halpha Equivalent Widths from the 3D-HST survey: evolution with redshift and dependence on stellar mass

We investigate the evolution of the Halpha equivalent width, EW(Halpha), with redshift and its dependence on stellar mass, taking advantage of the first data from the 3D-HST survey, a large spectroscopic Treasury program with the Hubble Space Telescope WFC3. Combining our Halpha measurements of 854 galaxies at 0.8<z<1.5 with those of ground based surveys at lower and higher redshift, we can consistently determine the evolution of the EW(Halpha) distribution from z=0 to z=2.2. We find that at all masses the characteristic EW(Halpha) is decreasing towards the present epoch, and that at each redshift the EW(Halpha) is lower for high-mass galaxies. We measure a slope of EW(Halpha) ~ (1+z)^(1.8) with little mass dependence. Qualitatively, this measurement is a model-independent confirmation of the evolution of star forming galaxies with redshift. A quantitative conversion of EW(Halpha) to sSFR is very model dependent, because of differential reddening corrections between the continuum SED and the Balmer lines. The observed EW(Halpha) can be reproduced with a simple model in which the SFR for galaxies rises to the epoch of z~2.5 and then decreases with time to z = 0. The model implies that the EW(Halpha) rises to 400 A at z=8. The sSFR evolves faster than EW(Halpha), as the mass-to-light ratio also evolves with redshift. In this context, we find that the sSFR evolves as (1+z)^(3.2), nearly independent of mass, consistent with previous reddening insensitive estimates. We confirm previous results that the observed slope of the sSFR-z relation is steeper than the one predicted by models, but models and observations agree in finding little mass dependence.

preprint2012arXiv

HST narrow-band search for extended Ly-alpha emission around two z>6 quasars

We search for extended Ly-alpha emission around two z>6 quasars, SDSS J1030+0524 (z=6.309) and SDSS J1148+5251 (z=6.419) using WFC3 narrow-band filters on board the Hubble Space Telescope. For each quasar, we collected two deep, narrow-band images, one sampling the Ly-alpha line+continuum at the quasar redshifts and one of the continuum emission redwards of the line. After carefully modeling the Point Spread Function, we find no evidence for extended Ly-alpha emission. These observations set 2-sigma limits of L(Ly-alpha, extended) < 3.2 x 10^{44} erg/s for J1030+0524 and L(Ly-alpha, extended) < 2.5 x 10^{44} erg/s for J1148+5251. Given the star formation rates typically inferred from (rest-frame) far-infrared measurements of z~6 quasars, these limits are well below the intrinsic bright Ly-alpha emission expected from the recombination of gas photoionized by the quasars or by the star formation in the host galaxies, and point towards significant Ly-alpha suppression or dust attenuation. However, small extinction values have been observed along the line of sight to the nuclei, thus reddening has to be coupled with other mechanisms for Ly-alpha suppression (e.g., resonance scattering). No Ly-alpha emitting companions are found, down to a 5-sigma sensitivity of ~ 1 x 10^{-17} erg/s/cm^2/arcsec^2 (surface brightness) and ~ 5 x 10^{-17} erg/s/cm^2 (assuming point sources).

preprint2012arXiv

Mass estimates from stellar proper motions: The mass of $ω$ Centauri

We lay out and apply methods to use proper motions of individual kinematic tracers for estimating the dynamical mass of star clusters. We first describe a simple projected mass estimator and then develop an approach that evaluates directly the likelihood of the discrete kinematic data given the model predictions. Those predictions may come from any dynamical modelling approach, and we implement an analytic King model, a spherical isotropic Jeans equation model and an axisymmetric, anisotropic Jeans equation model.We apply these approaches to the enigmatic globular cluster omega Centauri, combining the proper motion from van Leeuwen et al (2000) with improved photometric cluster membership probabilities. We show that all mass estimates based on spherical isotropic models yield $(4.55\pm 0.1) \times 10^6 M_{\odot} [D/5.5 \pm 0.2 kpc]^3$, where our modelling allows us to show how the statistical precision of this estimate improves as more proper motion data of lower signal-to-noise are included. MLM predictions, based on an anisotropic axisymmetric Jeans model, indicate for $ω$ Cen that the inclusion of anisotropies is not important for the mass estimates, but that accounting for the flattening is: flattened models imply $(4.05\pm 0.1) \times 10^6 M_{\odot} [D/5.5 \pm 0.2 kpc]^3$, 10% lower than when restricting the analysis to a spherical model. The best current distance estimates imply an additional uncertainty in the mass estimate of 12%.

preprint2012arXiv

On the nature of the stellar bridge between Leo IV and Leo V

We present a revised analysis of a speculated stellar bridge between the Milky Way dwarf galaxies Leo IV and Leo V. Using data acquired with Subaru/Suprime-Cam over a 1deg x 4deg field encompassing the two satellites and the region in between, we confirm our previous detection of a stellar overdensity between Leo IV and Leo V (de Jong et al. 2010). The larger area coverage and improved depth of our current dataset allow for an improved analysis of the stellar overdensity that had previously appeared to bridge the two galaxies. A main-sequence turn-off feature visible in the stacked colour-magnitude diagram of the contiguously observed Subaru fields reveals an extended stellar structure at a distance of approximately 20 kpc. Its angular proximity to the Virgo overdensity, as well as a good correspondence in distance and metallicity, suggests that the smaller structure we detect may be associated with the much larger Virgo stellar overdensity.

preprint2012arXiv

PHAT Stellar Cluster Survey I. Year 1 Catalog and Integrated Photometry

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an on-going Hubble Space Telescope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, representing more than a factor of four increase over previous catalogs within the current survey area (390 arcmin^2). This work presents results derived from the first \sim25% of the survey data; we estimate that the final sample will include \sim2500 clusters. For the Year 1 objects, we present a catalog with positions, radii, and six-band integrated photometry. Along with a general characterization of the cluster luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and highlight a number of individually interesting clusters found in the Year 1 search.

preprint2012arXiv

Shape Evolution of Massive Early-Type Galaxies: Confirmation of Increased Disk Prevalence at z>1

We use high-resolution K-band VLT/HAWK-I imaging over 0.25 square degrees to study the structural evolution of massive early-type galaxies since z~1. Mass-selected samples, complete down to log(M/M_sun)~10.7 such that `typical' L* galaxies are included at all redshifts, are drawn from pre-existing photometric redshift surveys. We then separated the samples into different redshift slices and classify them as late- or early-type galaxies on the basis of their specific star-formation rate. Axis-ratio measurements for the ~400 early-type galaxies in the redshift range 0.6<z<1.8 are accurate to 0.1 or better. The projected axis-ratio distributions are then compared with lower redshift samples. We find strong evidence for evolution of the population properties: early-type galaxies at z>1 are, on average, flatter than at z<1 and the median projected axis ratio at a fixed mass decreases with redshift. However, we also find that at all epochs z<~2 the very most massive early-type galaxies (log(M/M_sun)>11.3) are the roundest, with a pronounced lack among them of galaxies that are flat in projection. Merging is a plausible mechanism that can explain both results: at all epochs merging is required for early-type galaxies to grow beyond log(M/M_sun)~11.3, and all early types over time gradually and partially loose their disk-like characteristics.

preprint2012arXiv

Slicing The Monoceros Overdensity with Suprime-Cam

We derive distance, density and metallicity distribution of the stellar Monoceros Overdensity (MO) in the outer Milky Way, based on deep imaging with the Subaru Telescope. We applied CMD fitting techniques in three stripes at galactic longitudes: l=130 deg, 150 deg, 170 deg; and galactic latitudes: +15 < b [deg] < +25 . The MO appears as a wall of stars at a heliocentric distance of ~ 10.1\pm0.5 kpc across the observed longitude range with no distance change. The MO stars are more metal rich ([Fe/H] ~ -1.0) than the nearby stars at the same latitude. These data are used to test three different models for the origin of the MO: a perturbed disc model, which predicts a significant drop in density adjacent to the MO that is not seen; a basic flared disc model, which can give a good match to the density profile but the MO metallicity implies the disc is too metal rich to source the MO stars; and a tidal stream model, which bracket the distances and densities we derive for the MO, suggesting that a model can be found that would fully fit the MO data. Further data and modeling will be required to confirm or rule out the MO feature as a stream or as a flaring of the disc.

preprint2012arXiv

Spatially resolved Halpha maps and sizes of 57 strongly star-forming galaxies at z~1 from 3D-HST: evidence for rapid inside-out assembly of disk galaxies

We investigate the build-up of galaxies at z~1 using maps of Halpha and stellar continuum emission for a sample of 57 galaxies with rest-frame Halpha equivalent widths >100 Angstroms in the 3D-HST grism survey. We find that the Halpha emission broadly follows the rest-frame R-band light but that it is typically somewhat more extended and clumpy. We quantify the spatial distribution with the half-light radius. The median Halpha effective radius r_e(Halpha) is 4.2+-0.1 kpc but the sizes span a large range, from compact objects with r_e(Halpha) ~ 1.0 kpc to extended disks with r_e(Halpha) ~ 15 kpc. Comparing Halpha sizes to continuum sizes, we find <r_e(Halpha)/r_e(R)>=1.3+-0.1 for the full sample. That is, star formation, as traced by Halpha, typically occurs out to larger radii than the rest-frame R-band stellar continuum; galaxies are growing their radii and building up from the inside out. This effect appears to be somewhat more pronounced for the largest galaxies. Using the measured Halpha sizes, we derive star formation rate surface densities. We find that they range from ~0.05 Msun yr^{-1} kpc^{-2} for the largest galaxies to ~5 Msun yr^{-1} kpc^{-2} for the smallest galaxies, implying a large range in physical conditions in rapidly star-forming z~1 galaxies. Finally, we infer that all galaxies in the sample have very high gas mass fractions and stellar mass doubling times < 500 Myr. Although other explanations are also possible, a straightforward interpretation is that we are simultaneously witnessing the rapid formation of compact bulges and large disks at z~1.

preprint2012arXiv

The heating of dust by old stellar populations in the Bulge of M31

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust heating occurs in the central regions of galaxy spheroids that are essentially devoid of young stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean value in the disk of 17 pm 1K to \sim35K at the centre. UV to near-IR imaging of the central few kpc shows directly the absence of young stellar populations, delineates the radial profile of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, U\ast(r), as a function of radius. The increasing temperature profile in the centre matches that expected from the stellar heating, i.e. that the dust heating and cooling rates track each other over nearly two orders of magnitude in U\ast. The modelled dust heating is in excess of the observed dust temperatures, suggesting that it is more than sufficient to explain the observed IR emission. Together with the wavelength dependent absorption cross section of the dust, this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. This analysis shows that neither young stellar populations nor stellar near-UV radiation are necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat the dust. To the extent which these results pertain to the tenuous dust found in the centres of early-type galaxies remains yet to be explored.

preprint2012arXiv

The Milky Way has no thick disk

Different stellar sub-populations of the Milky Way's stellar disk are known to have different vertical scale heights, their thickness increasing with age. Using SEGUE spectroscopic survey data, we have recently shown that mono-abundance sub-populations, defined in the [α/Fe]-[Fe/H] space, are well described by single exponential spatial-density profiles in both the radial and the vertical direction; therefore any star of a given abundance is clearly associated with a sub-population of scale height h_z. Here, we work out how to determine the stellar surface-mass density contributions at the solar radius R_0 of each such sub-population, accounting for the survey selection function, and for the fraction of the stellar population mass that is reflected in the spectroscopic target stars given populations of different abundances and their presumed age distributions. Taken together, this enables us to derive Σ_{R_0}(h_z), the surface-mass contributions of stellar populations with scale height h_z. Surprisingly, we find no hint of a thin-thick disk bi-modality in this mass-weighted scale-height distribution, but a smoothly decreasing function, approximately Σ_{R_0}(h_z)\propto \exp(-h_z), from h_z ~ 200 pc to h_z ~ 1 kpc. As h_z is ultimately the structurally defining property of a thin or thick disk, this shows clearly that the Milky Way has a continuous and monotonic distribution of disk thicknesses: there is no 'thick disk' sensibly characterized as a distinct component. We discuss how our result is consistent with evidence for seeming bi-modality in purely geometric disk decompositions, or chemical abundances analyses. We constrain the total visible stellar surface-mass density at the Solar radius to be Σ^*_{R_0} = 30 +/- 1 M_\odot pc^{-2}.

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 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 spatial structure of mono-abundance sub-populations of the Milky Way disk

The spatial, kinematic, and elemental-abundance structure of the Milky Way's stellar disk is complex, and has been difficult to dissect with local spectroscopic or global photometric data. Here, we develop and apply a rigorous density modeling approach for Galactic spectroscopic surveys that enables investigation of the global spatial structure of stellar sub-populations in narrow bins of [α/Fe] and [Fe/H], using 23,767 G-type dwarfs from SDSS/SEGUE. We fit models for the number density of each such mono-abundance component, properly accounting for the complex spectroscopic SEGUE sampling of the underlying stellar population. We find that each mono-abundance sub-population has a simple spatial structure that can be described by a single exponential in both the vertical and radial direction, with continuously increasing scale heights (~200 pc to 1 kpc) and decreasing scale lengths (>4.5 kpc to 2 kpc) for increasingly older sub-populations, as indicated by their lower metallicities and [α/Fe] enhancements. That the abundance-selected sub-components with the largest scale heights have the shortest scale lengths is in sharp contrast with purely geometric `thick--thin disk' decompositions. To the extent that [α/Fe] is an adequate proxy for age, our results directly show that older disk sub-populations are more centrally concentrated, which implies inside-out formation of galactic disks. The fact that the largest scale-height sub-components are most centrally concentrated in the Milky Way is an almost inevitable consequence of explaining the vertical structure of the disk through internal evolution. Whether the simple spatial structure of the mono-abundance sub-components, and the striking correlations between age, scale length, and scale height can be plausibly explained by satellite accretion or other external heating remains to be seen.

preprint2012arXiv

The vertical motions of mono-abundance sub-populations in the Milky Way disk

We present the vertical kinematics of stars in the Milky Way's stellar disk inferred from SDSS/SEGUE G-dwarf data, deriving the vertical velocity dispersion, σ_z, as a function of vertical height |z| and Galactocentric radius R for a set of 'mono-abundance' sub-populations of stars with very similar elemental abundances [α/Fe] and [Fe/H]. We find that all components exhibit nearly isothermal kinematics in |z|, and a slow outward decrease of the vertical velocity dispersion: σ_z (z,R|[α/Fe],[Fe/H]) ~ σ_z ([α/Fe],[Fe/H]) x \exp (-(R-R_0)/7 kpc}). The characteristic velocity dispersions of these components vary from ~ 15 km/s for chemically young, metal-rich stars, to >~ 50 km/s for metal poor stars. The mean σ_z gradient away from the mid plane is only 0.3 +/- 0.2 km/s/kpc. We find a continuum of vertical kinetic temperatures (~σ^2_z) as function of ([α/Fe],[Fe/H]), which contribute to the stellar surface mass density as Σ_{R_0}(σ^2_z) ~ \exp(-σ^2_z). The existence of isothermal mono-abundance populations with intermediate dispersions reject the notion of a thin-thick disk dichotomy. This continuum of disks argues against models where the thicker disk portions arise from massive satellite infall or heating; scenarios where either the oldest disk portion was born hot, or where internal evolution plays a major role, seem the most viable. The wide range of σ_z ([α/Fe],[Fe/H]) combined with a constant σ_z(z) for each abundance bin provides an independent check on the precision of the SEGUE abundances: δ_[α/Fe] ~ 0.07 dex and δ_[Fe/H] ~ 0.15 dex. The radial decline of the vertical dispersion presumably reflects the decrease in disk surface-mass density. This measurement constitutes a first step toward a purely dynamical estimate of the mass profile the disk in our Galaxy. [abridged]

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 constant limiting mass scale for flat early-type galaxies from z=1 to z=0: density evolves but shapes do not

We measure the evolution in the intrinsic shape distribution of early-type galaxies from z~1 to z~0 by analyzing their projected axis-ratio distributions. We extract a low-redshift sample (0.04 < z < 0.08) of early-type galaxies with very low star-formation rates from the SDSS, based on a color-color selection scheme and verified through the absence of emission lines in the spectra. The inferred intrinsic shape distribution of these early-type galaxies is strongly mass dependent: the typical short-to-long intrinsic axis-ratio of high-mass early-type galaxies (>1e11 M_sun) is 2:3, where as at masses below 1e11 M_sun this ratio narrows to 1:3, or more flattened galaxies. In an entirely analogous manner we select a high-redshift sample (0.6 < z < 0.8) from two deep-field surveys: GEMS and COSMOS. We find a seemingly universal mass of ~1e11 M_sun for highly flatted early-type systems at all redshifts. This implies that the process that grows an early-type galaxy above this ceiling mass involves forming round systems. Using both parametric and non-parametric tests, we find no evolution in the projected axis-ratio distribution for galaxies with masses >3e10 M_sun with redshift. At the same time, our samples imply an increase of 2-3x in co-moving number density for early-type galaxies at masses >3e10 M_sun, in agreement with previous studies. Given the direct connection between the axis-ratio distribution and the underlying bulge-to-disk ratio distribution, our findings imply that the number density evolution of early-type galaxies is not exclusively driven by the emergence of either bulge- or disk-dominated galaxies, but rather by a balanced mix that depends only on the stellar mass of the galaxy. The challenge for galaxy formation models is to reproduce this overall non-evolving ratio of flattened to round early-type galaxies in the context of a continually growing population.

preprint2011arXiv

A Deep Photometric Look at Two of Andromeda's Dwarf Spheroidals: X and XVII

We use deep wide-field photometry from the Large Binocular Camera to study the stellar and structural properties of the recently discovered Andromeda X and Andromeda XVII (And X and And XVII) dwarf galaxies. Using the mean apparent magnitude of the horizontal branch (HB), we derive distances of 621 +- 20 kpc to And X and 734+- 23 kpc to And XVII, closer by >60 kpc than the previous estimates which were based on red giant branch (RGB) observations. Thus our results warrant against the use of the RGB tip method for determining distances to systems with sparsely populated RGBs, and show how crucial HB observations are in obtaining accurate distances in systems such as these. We find that And X is a relatively faint (MV = -7.36), highly elongated (e = 0.48) system at a distance of 174 +- 62 kpc from Andromeda. And XVII is brighter (MV = -8.61) with an M31-centric distance of 73 kpc which makes it one of the closest satellites to Andromeda. Both galaxies are metal-poor: we derive <[Fe/H]>=-2.2 for And X, while And XVII shows <[Fe/H]> = -2.0, consistent with the relation of higher luminosity dwarfs being more metal- rich. Additionally, both galaxies show considerable intrinsic spreads in metallicity (0.2 and 0.3 dex for And X and And XVII respectively), consistent with multiple stellar populations.

preprint2011arXiv

Astrometry with the MCAO instrument MAD - An analysis of single-epoch data obtained in the layer-oriented mode

Context: Current instrument developments at the largest telescopes worldwide have provisions for Multi-Conjugated Adaptive Optics (MCAO) modules. The large field of view and more uniform correction provided by these systems is not only highly beneficial for photometric studies but also for astrometric analysis of, e.g., large dense clusters and exoplanet detection and characterization. The Multi-conjugated Adaptive optics Demonstrator (MAD) is the first such instrument and was temporarily installed and tested at the ESO/VLT in 2007. We analyzed two globular cluster data sets in terms of achievable astrometric precision. Data were obtained in the layer-oriented correction mode, one in full MCAO correction mode with two layers corrected (NGC 6388) and the other applying ground-layer correction only (47 Tuc). Aims: We aim at analyzing the first available MCAO imaging data in the layer-oriented mode obtained with the MAD instrument in terms of astrometric precision and stability. Methods: We calculated Strehl maps for each frame in both data sets. Distortion corrections were performed and the astrometric precision was analyzed by calculating mean stellar positions over all frames and by investigation of the positional residuals present in each frame after transformation to a master-coordinate-frame. Results: The mean positional precision for stars between K = 14-18 mag is ~1.2 mas in the full MCAO correction mode data of the cluster NGC 6388. The precision measured in the GLAO data (47 Tuc) reaches ~1.0 mas for stars corresponding to 2MASS K magnitudes between 9 and 12. The observations were such that stars in these magnitude ranges correspond to the same detector flux range. The jitter movement used to scan a larger field of view introduced additional distortions in the frames, leading to a degradation of the achievable precision.

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

Evidence for non-evolving FeII/MgII ratios in rapidly accreting z~6 QSOs

[abridged] Quasars (QSOs) at the highest known redshift (z~6) are unique probes of the early growth of supermassive black holes (BHs). Until now, only the most luminous QSOs have been studied, often one object at a time. Here we present the most extensive consistent analysis to date of z>4 QSOs with observed NIR spectra, combining three new z~6 objects from our ongoing VLT-ISAAC program with nineteen 4<z<6.5 sources from the literature. The new sources extend the existing SDSS sample towards the faint end of the QSO luminosity function. Using a maximum likelihood fitting routine optimized for our spectral decomposition, we estimate the black hole mass (MBH), the Eddington ratio (defined as Lbol/LEdd) and the FeII/MgII line ratio, a proxy for the chemical abundance, to characterize both the central object and the broad line region gas. The QSOs in our sample host BHs with masses of ~10^9 M\odot that are accreting close to the Eddington luminosity, consistent with earlier results. We find that the distribution of observed Eddington ratios is significantly different than that of a luminosity-matched comparison sample of SDSS QSOs at lower redshift (0.35<z<2.25): the average <log(Lbol/LEdd)>=-0.37 (Lbol/LEdd~0.43) with a scatter of 0.20 dex for the z>4 sample and the <log(Lbol/LEdd)>=-0.80 (Lbol/LEdd~0.16) with a scatter of 0.24 dex for the 0.35<z<2.25 sample. This implies that, at a given luminosity, the MBH at high-z is typically lower than the average MBH of the lower-redshift population, i.e. the z>4 sources are accreting significantly faster than the lower-redshift ones. We show that the derived FeII/MgII ratios depend sensitively on the performed analysis: our self-consistent, homogeneous analysis significantly reduces the FeII/MgII scatter found in previous studies. The measured FeII/MgII line ratios show no sign of evolution with cosmic time in the redshift range 4<z<6.5 [...]

preprint2011arXiv

First Results from the 3D-HST Survey: The Striking Diversity of Massive Galaxies at z>1

We present first results from the 3D-HST program, a near-IR spectroscopic survey performed with the Wide Field Camera 3 on the Hubble Space Telescope. We have used 3D-HST spectra to measure redshifts and Halpha equivalent widths for a stellar mass-limited sample of 34 galaxies at 1<z<1.5 with M(stellar)>10^11 M(sun) in the COSMOS, GOODS, and AEGIS fields. We find that a substantial fraction of massive galaxies at this epoch are forming stars at a high rate: the fraction of galaxies with Halpha equivalent widths >10 A is 59%, compared to 10% among SDSS galaxies of similar masses at z=0.1. Galaxies with weak Halpha emission show absorption lines typical of 2-4 Gyr old stellar populations. The structural parameters of the galaxies, derived from the associated WFC3 F140W imaging data, correlate with the presence of Halpha: quiescent galaxies are compact with high Sersic index and high inferred velocity dispersion, whereas star-forming galaxies are typically large two-armed spiral galaxies, with low Sersic index. Some of these star forming galaxies might be progenitors of the most massive S0 and Sa galaxies. Our results challenge the idea that galaxies at fixed mass form a homogeneous population with small scatter in their properties. Instead we find that massive galaxies form a highly diverse population at z>1, in marked contrast to the local Universe.

preprint2011arXiv

Observational Evidence Against Long-Lived Spiral Arms in Galaxies

We test whether the spiral patterns apparent in many large disk galaxies should be thought of as dynamical features that are stationary in a co-rotating frame for > t_{dyn}, as implied by the density wave approach for explaining spiral arms. If such spiral arms have enhanced star formation (SF), observational tracers for different stages of the SF sequence should show a spatial ordering, from up-stream to downstream in the corotating frame: dense HI, CO, tracing molecular hydrogen gas, 24 micron emission tracing enshrouded SF and UV emission tracing unobscured young stars. We argue that such a spatial ordering should be reflected in the angular cross-correlation (CC, in polar coordinates) using all azimuthal positions among pairs of these tracers; the peak of the CC should be offset from zero, in different directions inside and outside the corotation radius. Recent spiral SF simulations by Dobbs & Pringle, show explicitly that for the case of a stationary spiral arm potential such angular offsets between gas and young stars of differing ages should be observable as cross-correlation offsets. We calculate the angular cross-correlations for different observational SF sequence tracers in 12 nearby spiral galaxies, drawing on a data set with high quality maps of the neutral gas HI, THINGS), molecular gas (CO, HERACLES) along with 24 micron emission (Spitzer, SINGS); we include FUV images (GALEX) and 3.6 $μ$m emission (Spitzer, IRAC) for some galaxies, tracing aging stars and longer timescales. In none of the resulting tracer cross-correlations for this sample do we find systematic angular offsets, which would be expected for a stationary dynamical spiral pattern of well-defined pattern speed. This result indicates that spiral density waves in their simplest form are not an important aspect of explaining spirals in large disk galaxies.

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

Redshift Evolution of the Galaxy Velocity Dispersion Function

We present a study of the evolution of the galaxy Velocity Dispersion Function (VDF) from z=0 to z=1.5 using photometric data from the UKIDSS Ultra Deep Survey (UDS) and Newfirm Medium Band Survey (NMBS) COSMOS surveys. The VDF has been measured locally using direct kinematic measurements from the Sloan Digital Sky Survey, but direct studies of the VDF at high redshift are difficult as they require velocity dispersion measurements of many thousands of galaxies. Taylor et al. (2010) demonstrated that dynamical and stellar mass are linearly related when the structure of the galaxy is accounted for. We show that the stellar mass, size and Sérsic index can reliably predict the velocity dispersions of SDSS galaxies. We apply this relation to galaxies at high redshift and determine the evolution of the inferred VDF. We find that the VDF at z~0.5 is very similar to the VDF at z=0. At higher redshifts, we find that the number density of galaxies with dispersions <~200 km/s is lower, but the number of high dispersion galaxies is constant or even higher. At fixed cumulative number density, the velocity dispersions of galaxies with logN[Mpc^{-3}]<-3.5 increase with time by a factor of ~1.4 from z~1.5-0, whereas the dispersions of galaxies with lower number density are approximately constant or decrease with time. The VDF appears to show less evolution than the stellar mass function, particularly at the lowest number densities. We note that these results are still somewhat uncertain and we suggest several avenues for further calibrating the inferred velocity dispersions.

preprint2011arXiv

Resolving the dynamical mass of a z~1.3 QSO host galaxy using SINFONI and Laser Guide Star assisted Adaptive Optics

Recent studies of the tight scaling relations between the masses of supermassive black holes and their host galaxies have suggested that in the past black holes constituted a larger fraction of their host galaxies' mass. However, these arguments are limited by selection effects and difficulties in determining robust host galaxy masses at high redshifts. Here we report the first results of a new, complementary diagnostic route: we directly determine a dynamical host galaxy mass for the z=1.3 luminous quasar J090543.56+043347.3 through high-spatial-resolution (0.47", 4kpc FWHM) observations of the host galaxy gas kinematics over 30x40 kpc using ESO/VLT/SINFONI with LGS/AO. Combining our result of M_dyn = 2.05+1.68_0.74 x 10^11 M_sun (within a radius 5.25 +- 1.05 kpc) with M_BH,MgII = 9.02 \pm 1.43 x 10^8 M_sun, M_BH,Halpha = 2.83 +1.93-1.13 x 10^8 M_sun, we find that the ratio of black hole mass to host galaxy dynamical mass for J090543.56+043347.3 matches the present-day relation for M_BH vs. M_Bulge,Dyn, well within the IR scatter, deviating at most a factor of two from the mean. J090543.56+043347.3 displays clear signs of an ongoing tidal interaction and of spatially extended star formation at a rate of 50-100 M_sun/yr, above the cosmic average for a galaxy of this mass and redshift. We argue that its subsequent evolution may move J090543.56+043347.3 even closer to the z=0 relation for M_BH vs. M_Bulge,Dyn. Our results support the picture where any substantive evolution in these relations must occur prior to z~1.3. Having demonstrated the power of this modelling approach we are currently analyzing similar data on seven further objects to better constrain such evolution.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Color Variability of Quasars

We quantify quasar color-variability using an unprecedented variability database - ugriz photometry of 9093 quasars from SDSS Stripe 82, observed over 8 years at ~60 epochs each. We confirm previous reports that quasars become bluer when brightening. We find a redshift dependence of this blueing in a given set of bands (e.g. g and r), but show that it is the result of the flux contribution from less-variable or delayed emission lines in the different SDSS bands at different redshifts. After correcting for this effect, quasar color-variability is remarkably uniform, and independent not only of redshift, but also of quasar luminosity and black hole mass. The color variations of individual quasars, as they vary in brightness on year timescales, are much more pronounced than the ranges in color seen in samples of quasars across many orders of magnitude in luminosity. This indicates distinct physical mechanisms behind quasar variability and the observed range of quasar luminosities at a given black hole mass - quasar variations cannot be explained by changes in the mean accretion rate. We do find some dependence of the color variability on the characteristics of the flux variations themselves, with fast, low-amplitude, brightness variations producing more color variability. The observed behavior could arise if quasar variability results from flares or ephemeral hot spots in an accretion disc.

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

preprint2011arXiv

The First High Redshift Quasar from Pan-STARRS

We present the discovery of the first high redshift (z > 5.7) quasar from the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1 or PS1). This quasar was initially detected as an i dropoutout in PS1, confirmed photometrically with the SAO Widefield InfraRed Camera (SWIRC) at Arizona's Multiple Mirror Telescope (MMT) and the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) at the MPG 2.2 m telescope in La Silla. The quasar was verified spectroscopically with the the MMT Spectrograph, Red Channel and the Cassegrain Twin Spectrograph (TWIN) at the Calar Alto 3.5 m telescope. It has a redshift of 5.73, an AB z magnitude of 19.4, a luminosity of 3.8 x 10^47 erg/s and a black hole mass of 6.9 x 10^9 solar masses. It is a Broad Absorption Line quasar with a prominent Ly-beta peak and a very blue continuum spectrum. This quasar is the first result from the PS1 high redshift quasar search that is projected to discover more than a hundred i dropout quasars, and could potentially find more than 10 z dropout (z > 6.8) quasars.

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

preprint2011arXiv

Two-dimensional H_alpha kinematics of bulgeless disk galaxies

We present two-dimensional H_alpha velocity fields for 20 late-type, disk-dominated spiral galaxies, the largest sample to date with high-resolution H_alpha velocity fields for bulgeless disks. From these data we derive rotation curves and the location of the kinematic centers. The galaxy sample was selected to contain nucleated and non-nucleated galaxies, which allows us to investigate what impact the gas kinematics in the host disk have on the presence (or absence) of a nuclear star cluster. In general, we find that the velocity fields span a broad range of morphologies. While some galaxies show regular rotation, most have some degree of irregular gas motions. There appears to be no systematic difference in the kinematics of nucleated and non-nucleated disks. Due to the large fields of view of the integral field units we use, we are able to observe the flattening of the rotation curve in almost all of our sample galaxies. This makes modeling of the velocity fields relatively straight-forward. Due to the complexities of the velocity fields, we obtain reliable determinations of the kinematic center for only 6 of our 20 sample galaxies. For all of these the locations of the nuclear star cluster/photometric center and the kinematic center agree within the uncertainties. If we disregard all kinematically irregular galaxies, our study concludes that nuclear star clusters truly occupy the nuclei, or dynamical centers, of their hosts. Our results are thus consistent with in-situ formation of nuclear star clusters. Yet, many well-motivated formation scenarios for nuclear clusters invoke off-center cluster formation and subsequent sinking of clusters due to dynamical friction. In that case, our results imply that dynamical friction in the centers of bulgeless galaxies must be very effective in driving massive clusters to the kinematic center. (abridged)

preprint2011arXiv

What Sets the Sizes of the Faintest Galaxies?

We provide a comprehensive description and offer an explanation for the sizes of the faintest known galaxies in the universe, the dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda. After compiling a consistent data set of half-light radii (r_{1/2}) and luminosities, we describe the size-luminosity relation of dSphs by a log-normal distribution in r_{1/2} with a mean size that varies as a function of luminosity. Accounting for modest number statistics, measurement uncertainties and surface brightness limitations, we find that the size-luminosity relations of the Milky Way and Andromeda dSph populations are statistically indistinguishable, and also very similar: their mean sizes at a given stellar luminosity differ by no more than 30%. In addition, we find that the mean size, slope and scatter of this log-normal size description of Local Group dSphs matches onto the relation of more massive low-concentration galaxies. This suggests that the stellar sizes of dSphs are ultimately related to their overall initial baryonic angular momentum. To test this hypothesis we perform a series of high resolution N-body simulations that we couple with a semi-analytic model of galaxy formation. These predict the same mean size and slope as observed in dSph satellites. At the same time, these models predict that the size-luminosity distributions for satellite galaxies around similar host-halos must be similar providing a natural explanation as to why the size distributions of Milky Way and Andromeda satellites are similar. Although strong rotation is currently not observed in dSphs, this may well be consistent with our angular-momentum-based explanation for their sizes if the disks of these galaxies have become sufficiently stirred through tidal interaction.

preprint2010arXiv

A Cosmic Variance Cookbook

Deep pencil beam surveys (<1 deg^2) are of fundamental importance for studying the high-redshift universe. However, inferences about galaxy population properties are in practice limited by 'cosmic variance'. This is the uncertainty in observational estimates of the number density of galaxies arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Dz. Using a halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting function to compute cosmic variance as a function of the angular dimensions of the field, z, Dz and stellar mass m*. We also provide tabulated values and a software tool. We find that for GOODS at z=2 and with Dz=0.5 the relative cosmic variance of galaxies with m*>10^11 Msun is ~38%, while it is ~27% for GEMS and ~12% for COSMOS. For galaxies of m*~10^10 Msun the relative cosmic variance is ~19% for GOODS, ~13% for GEMS and ~6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z=2 for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is less serious.

preprint2010arXiv

A Giant Planet Around a Metal-poor Star of Extragalactic Origin

Stars in their late stage of evolution, such as Horizontal Branch stars, are still largely unexplored for planets. We report the detection of a planetary companion around HIP 13044, a very metal-poor star on the red Horizontal Branch, based on radial velocity observations with a high-resolution spectrograph at the 2.2-m MPG/ESO telescope. The star's periodic radial velocity variation of P=16.2 days caused by the planet can be distinguished from the periods of the stellar activity indicators. The minimum mass of the planet is 1.25 Jupiter masses and its orbital semi-major axis 0.116 AU. Because HIP 13044 belongs to a group of stars that have been accreted from a disrupted satellite galaxy of the Milky Way, the planet most likely has an extragalactic origin.

preprint2010arXiv

Arm & Interarm Star Formation in Spiral Galaxies

We investigate the relationship between spiral arms and star formation in the grand-design spirals NGC 5194 and NGC 628 and in the flocculent spiral NGC 6946. Filtered maps of near-IR (3.6 micron) emission allow us to identify "arm regions" that should correspond to regions of stellar mass density enhancements. The two grand-design spirals show a clear two-armed structure, while NGC 6946 is more complex. We examine these arm and interarm regions, looking at maps that trace recent star formation - far-ultraviolet (GALEX NGS) and 24 micron emission (Spitzer, SINGS) - and cold gas - CO (Heracles) and HI (Things). We find the star formation tracers and CO more concentrated in the spiral arms than the stellar 3.6 micron flux. If we define the spiral arms as the 25% highest pixels in the filtered 3.6 micron images, we find that the majority (60%) of star formation tracers occurs in the interarm regions; this result persists qualitatively even when considering the potential impact of finite data resolution and diffuse interarm 24 micron emission. Even with a generous definition of the arms (45% highest pixels), interarm regions still contribute at least 30% to the integrated star formation rate tracers. We look for evidence that spiral arms trigger star or cloud formation using the ratios of star formation rate (SFR, traced by a combination of FUV and 24 micron emission) to H_2 (traced by CO) and H_2 to HI. Any enhancement of SFR / M(H_2) in the arm region is very small (less than 10%) and the grand design spirals show no enhancement compared to the flocculent target. Arm regions do show a weak enhancement in H_2/HI compared to the interarm regions, but at a fixed gas surface density there is little clear enhancement in the H_2/HI ratio in the arm regions. Thus, it seems that spiral arms may only act to concentrate the gas to higher densities in the arms.

preprint2010arXiv

Classification of field dwarfs and giants in RAVE and its use in stellar stream detection

An efficient separation between dwarfs and giants in surveys of bright stars is important, especially for studies in which distances are estimated through photometric parallax relations. We use the available spectroscopic log g estimates from the second RAVE data release (DR2) to assign each star a probability for being a dwarf or subgiant/giant based on mixture model fits to the log g distribution in different color bins. We further attempt to use these stars as a labeled training set in order to classify stars which lack log g estimates into dwarfs and giants with a SVM algorithm. We assess the performance of this classification against different choices of the input feature vector. In particular, we use different combinations of reduced proper motions, 2MASS JHK, DENIS IJK and USNO-B B2R2 apparent magnitudes. Our study shows that -- for our color ranges -- the infrared bands alone provide no relevant information to separate dwarfs and giants. Even when optical bands and reduced proper motions are added, the fraction of true giants classified as dwarfs (the contamination) remains above 20%. Using only the dwarfs with available spectroscopic log g and distance estimates (the latter from Breddels et al. 2010), we then repeat the stream search by Klement, Fuchs & Rix (2008, KFR08), which assumed all stars were dwarfs and claimed the discovery of a new stellar stream at V = -160 km/s in a sample of 7015 stars from RAVE DR1. Our re-analysis of the pure DR2 dwarf sample exhibits an overdensity of 5 stars at the phase-space position of the KFR08 stream, with a metallicity distribution that appears inconsistent with that of stars at comparably low rotational velocities. Compared to several smooth Milky Way models, the mean standardized deviation of the KFR08 stream is only marginal at 1.6$\pm$0.4... (abbreviated)

preprint2010arXiv

Constraining the Milky Way potential with a 6-D phase-space map of the GD-1 stellar stream

The narrow GD-1 stream of stars, spanning 60 deg on the sky at a distance of ~10 kpc from the Sun and ~15 kpc from the Galactic center, is presumed to be debris from a tidally disrupted star cluster that traces out a test-particle orbit in the Milky Way halo. We combine SDSS photometry, USNO-B astrometry, and SDSS and Calar Alto spectroscopy to construct a complete, empirical 6-dimensional phase-space map of the stream. We find that an eccentric orbit in a flattened isothermal potential describes this phase-space map well. Even after marginalizing over the stream orbital parameters and the distance from the Sun to the Galactic center, the orbital fit to GD-1 places strong constraints on the circular velocity at the Sun's radius V_c=224 \pm 13 km/s and total potential flattening q_Φ=0.87^{+0.07}_{-0.04}. When we drop any informative priors on V_c the GD-1 constraint becomes V_c=221 \pm 18 km/s. Our 6-D map of GD-1 therefore yields the best current constraint on V_c and the only strong constraint on q_Φat Galactocentric radii near R~15 kpc. Much, if not all, of the total potential flattening may be attributed to the mass in the stellar disk, so the GD-1 constraints on the flattening of the halo itself are weak: q_{Φ,halo}>0.89 at 90% confidence. The greatest uncertainty in the 6-D map and the orbital analysis stems from the photometric distances, which will be obviated by Gaia.

preprint2010arXiv

Detailed Decomposition of Galaxy Images. II. Beyond Axisymmetric Models

We present a two-dimensional (2-D) fitting algorithm (GALFIT, Version 3) with new capabilities to study the structural components of galaxies and other astronomical objects in digital images. Our technique improves on previous 2-D fitting algorithms by allowing for irregular, curved, logarithmic and power-law spirals, ring and truncated shapes in otherwise traditional parametric functions like the Sersic, Moffat, King, Ferrer, etc., profiles. One can mix and match these new shape features freely, with or without constraints, apply them to an arbitrary number of model components and of numerous profile types, so as to produce realistic-looking galaxy model images. Yet, despite the potential for extreme complexity, the meaning of the key parameters like the Sersic index, effective radius or luminosity remain intuitive and essentially unchanged. The new features have an interesting potential for use to quantify the degree of asymmetry of galaxies, to quantify low surface brightness tidal features beneath and beyond luminous galaxies, to allow more realistic decompositions of galaxy subcomponents in the presence of strong rings and spiral arms, and to enable ways to gauge the uncertainties when decomposing galaxy subcomponents. We illustrate these new features by way of several case studies that display various levels of complexity.

preprint2010arXiv

Mapping the stellar structure of the Milky Way thick disk and halo using SEGUE photometry

We map the stellar structure of the Galactic thick disk and halo by applying color-magnitude diagram (CMD) fitting to photometric data from the SEGUE survey, allowing, for the first time, a comprehensive analysis of their structure at both high and low latitudes using uniform SDSS photometry. Incorporating photometry of all relevant stars simultaneously, CMD fitting bypasses the need to choose single tracer populations. Using old stellar populations of differing metallicities as templates we obtain a sparse 3D map of the stellar mass distribution at |Z|>1 kpc. Fitting a smooth Milky Way model comprising exponential thin and thick disks and an axisymmetric power-law halo allows us to constrain the structural parameters of the thick disk and halo. The thick-disk scale height and length are well constrained at 0.75+-0.07 kpc and 4.1+-0.4 kpc, respectively. We find a stellar halo flattening within ~25 kpc of c/a=0.88+-0.03 and a power-law index of 2.75+-0.07 (for 7<R_{GC}<~30 kpc). The model fits yield thick-disk and stellar halo densities at the solar location of rho_{thick,sun}=10^{-2.3+-0.1} M_sun pc^{-3} and rho_{halo,sun}=10^{-4.20+-0.05} M_sun pc^{-3}, averaging over any substructures. Our analysis provides the first clear in situ evidence for a radial metallicity gradient in the Milky Way's stellar halo: within R<~15 kpc the stellar halo has a mean metallicity of [Fe/H]=-1.6, which shifts to [Fe/H]=-2.2 at larger radii. Subtraction of the best-fit smooth and symmetric model from the overall density maps reveals a wealth of substructures at all latitudes, some attributable to known streams and overdensities, and some new. A simple warp cannot account for the low latitude substructure, as overdensities occur simultaneously above and below the Galactic plane. (abridged)

preprint2010arXiv

Observational Evidence from SDSS for a Merger Origin of the Milky Way's Thick Disk

We test competing models that aim at explaining the nature of stars in the Milky Way that are well away (|z|$\gtrsim$ 1kpc) from the midplane, the so-called thick disk: the stars may have gotten there through orbital migration, through satellite mergers and accretion, or through heating of pre-existing thin disk stars. Sales et al. (2009) proposed the eccentricity distribution of thick disk stars as a diagnostic to differentiate between these mechanisms. Drawing on SDSS DR7, we have assembled a sample of 34,223 G-dwarfs with 6-D phase-space information and metallicities, and have derived orbital eccentricities for them. Comparing the resulting eccentricity distributions, p(e|z), with the models, we find that: a) the observed p(e|z) is inconsistent with that predicted by orbital migration only, as there are more observed stars of high and of very low eccentricity; b) scenarios where the thick disk is made predominantly through abrupt heating of a pre-existing thin disk are also inconsistent, as they predict more high-eccentricity stars than observed; c) the observed p(e|z) fits well with a "gas-rich merger" scenario, where most thick disk stars were born from unsettled gas in situ.

preprint2010arXiv

Selecting Quasars by their Intrinsic Variability

We present a new and simple technique for selecting extensive, complete and pure quasar samples, based on their intrinsic variability. We parametrize the single-band variability by a power-law model for the light-curve structure function, with amplitude A and power-law index gamma. We show that quasars can be efficiently separated from other non-variable and variable sources by the location of the individual sources in the A-gamma plane. We use ~60 epochs of imaging data, taken over ~5 years, from the SDSS stripe 82 (S82) survey, where extensive spectroscopy provides a reference sample of quasars, to demonstrate the power of variability as a quasar classifier in multi-epoch surveys. For UV-excess selected objects, variability performs just as well as the standard SDSS color selection, identifying quasars with a completeness of 90% and a purity of 95%. In the redshift range 2.5<z<3, where color selection is known to be problematic, variability can select quasars with a completeness of 90% and a purity of 96%. This is a factor of 5-10 times more pure than existing color-selection of quasars in this redshift range. Selecting objects from a broad griz color box without u-band information, variability selection in S82 can afford completeness and purity of 92%, despite a factor of 30 more contaminants than quasars in the color-selected feeder sample. This confirms that the fraction of quasars hidden in the 'stellar locus' of color-space is small. To test variability selection in the context of Pan-STARRS 1 (PS1) we created mock PS1 data by down-sampling the S82 data to just 6 epochs over 3 years. Even with this much sparser time sampling, variability is an encouragingly efficient classifier. For instance, a 92% pure and 44% complete quasar candidate sample is attainable from the above $griz$-selected catalog.

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

Stellar Tidal Streams in Spiral Galaxies of the Local Volume: A Pilot Survey with Modest Aperture Telescopes

[Abridged] Within the hierarchical framework for galaxy formation, minor merging and tidal interactions are expected to shape all large galaxies to the present day. As a consequence, most seemingly normal disk galaxies should be surrounded by spatially extended stellar 'tidal features' of low surface brightness. As part of a pilot survey for such interaction signatures, we have carried out ultra deep, wide field imaging of 8 isolated spiral galaxies in the Local Volume, with data taken at small (D=0.1-0.5m) robotic telescopes that provide exquisite surface brightness sensitivity (mu_V)~28.5$ mag/arcsec^2). This initial observational effort has led to the discovery of six previously undetected extensive (to ~30 kpc) stellar structures in the halos surrounding these galaxies, likely debris from tidally disrupted satellites. In addition, we confirm and clarify several enormous stellar over-densities previously reported in the literature, but never before interpreted as tidal streams. Even this pilot sample of galaxies exhibits strikingly diverse morphological characteristics of these extended stellar features: great circle-like features that resemble the Sagittarius stream surrounding the Milky Way, remote shells and giant clouds of presumed tidal debris far beyond the main stelar body, as well as jet-like features emerging from galactic disks. A qualitative comparison with available simulations set in a Lambda-Cold Dark Matter cosmology shows that the extraordinary variety of stellar morphologies detected in this pilot survey matches that seen in those simulations. The common existence of these tidal features around 'normal' disk galaxies and the morphological match to the simulations constitutes new evidence that these theoretical models also apply to a large number of other Milky Way-mass disk galaxies in the Local Volume.

preprint2010arXiv

The enigmatic pair of dwarf galaxies Leo IV and Leo V: coincidence or common origin?

We have obtained deep photometry in two 1x1 degree fields covering the close pair of dwarf spheroidal galaxies (dSph) Leo IV and Leo V and part of the area in between. We find that both systems are significantly larger than indicated by previous measurements based on shallower data and also significantly elongated. With half-light radii of r_h=4'.6 +- 0'.8 (206 +- 36 pc) and r_h=2'.6 +- 0'.6 (133 +- 31 pc), respectively, they are now well within the physical size bracket of typical Milky Way dSph satellites. Their ellipticities of epsilon ~0.5 are shared by many faint (M_V>-8) Milky Way dSphs. The large spatial extent of our survey allows us to search for extra-tidal features with unprecedented sensitivity. The spatial distribution of candidate red giant branch and horizontal branch stars is found to be non-uniform at the ~3 sigma level. This substructure is aligned along the direction connecting the two systems, indicative of a possible `bridge' of extra-tidal material. Fitting the stellar distribution with a linear Gaussian model yields a significance of 4 sigma for this overdensity, a most likely FWHM of ~16 arcmin and a central surface brightness of ~32 mag arcsec^{-2}. We investigate different scenarios to explain the close proximity of Leo IV and Leo V and the possible tidal bridge between them. Orbit calculations demonstrate that they are unlikely to be remnants of a single disrupted progenitor, while a comparison with cosmological simulations shows that a chance collision between unrelated subhalos is negligibly small. Leo IV and Leo V could, however, be a bound `tumbling pair' if their combined mass exceeds 8 +- 4 x 10^9 M_sun. The scenario of an internally interacting pair appears to be the most viable explanation for this close celestial companionship. (abridged)

preprint2010arXiv

The visitor from an ancient galaxy: A planetary companion around an old, metal-poor red horizontal branch star

We report the detection of a planetary companion around HIP 13044, a metal-poor red horizontal branch star belonging to a stellar halo stream that results from the disruption of an ancient Milky Way satellite galaxy. The detection is based on radial velocity observations with FEROS at the 2.2-m MPG/ESO telescope. The periodic radial velocity variation of P=16.2 days can be distinguished from the periods of the stellar activity indicators. We computed a minimum planetary mass of 1.25 Jupiter masses and an orbital semimajor axis of 0.116 AU for the planet. This discovery is unique in three aspects: First, it is the first planet detection around a star with a metallicity much lower than few percent of the solar value; second, the planet host star resides in a stellar evolutionary stage that is still unexplored in the exoplanet surveys; third, the planetary system HIP 13044 most likely has an extragalactic origin in a disrupted former satellite of the Milky Way.

preprint2009arXiv

A quantitative explanation of the observed population of Milky Way satellite galaxies

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.

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

Measuring the undetectable: Proper motions and parallaxes of very faint sources

The near future of astrophysics involves many large solid-angle, multi-epoch, multi-band imaging surveys. These surveys will, at their faint limits, have data on large numbers of sources that are too faint to be detected at any individual epoch. Here we show that it is possible to measure in multi-epoch data not only the fluxes and positions, but also the parallaxes and proper motions of sources that are too faint to be detected at any individual epoch. The method involves fitting a model of a moving point source simultaneously to all imaging, taking account of the noise and point-spread function in each image. By this method it is possible to measure the proper motion of a point source with an uncertainty close to the minimum possible uncertainty given the information in the data, which is limited by the point-spread function, the distribution of observation times (epochs), and the total signal-to-noise in the combined data. We demonstrate our technique on multi-epoch Sloan Digital Sky Survey imaging of the SDSS Southern Stripe. We show that we can distinguish very red brown dwarfs by their proper motions from very high-redshift quasars more than $1.6\mag$ fainter than with traditional technique on these SDSS data, and with better better fidelity than by multi-band imaging alone. We re-discover all 10 known brown dwarfs in our sample and present 9 new candidate brown dwarfs, identified on the basis of high proper motion.

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

Resolved maps of stellar mass and SED of galaxies from optical/NIR imaging and SPS models

We report on the method developed by Zibetti, Charlot and Rix (2009) to construct resolved stellar mass maps of galaxies from optical and NIR imaging. Accurate pixel-by-pixel colour information (specifically g-i and i-H) is converted into stellar mass-to-light ratios with typical accuracy of 30%, based on median likelihoods derived from a Monte Carlo library of 50,000 stellar population synthesis models that include dust and updated TP-AGB phase prescriptions. Hence, surface mass densities are computed. In a pilot study, we analyze 9 galaxies spanning a broad range of morphologies. Among the main results, we find that: i) galaxies appear much smoother in stellar mass maps than at any optical or NIR wavelength; ii) total stellar mass estimates based on unresolved photometry are biased low with respect to the integral of resolved stellar mass maps, by up to 40%, due to dust obscured regions being under-represented in global colours; iii) within a galaxy, on local scales colours correlate with surface stellar mass density; iv) the slope and tightness of this correlation reflect/depend on the morphology of the galaxy.

preprint2009arXiv

Resolved stellar mass maps of galaxies. I: method and implications for global mass estimates

(Abridged) We introduce a novel technique to construct spatially resolved maps of stellar mass surface density in galaxies based on optical and near IR imaging. We use optical/NIR colour(s) to infer effective stellar mass-to-light ratios (M/L) at each pixel, which are then multiplied by the surface brightness to obtain the local stellar surface mass density. We build look-up tables to express M/L as a function of colour(s) by marginalizing over a Monte Carlo library of 50,000 stellar population synthesis (SPS) models by Charlot & Bruzual (2007), which include a revised prescription for the TP-AGB stellar evolutionary phase, with a wide range of dust exinctions. In order to extract reliable flux and colour information at any position in the galaxy, we perform a median adaptive smoothing of the images that preserves the highest possible spatial resolution. As the most practical and robust, and hence fiducial method, we express the M/L in the H band as a function of (g-i) and (i-H). Stellar mass maps computed in this way have a typical accuracy of 30 per cent or less at any given pixel, determined from the scatter in the models. We compare maps obtained with our fiducial method with those derived using other combinations of bandpasses and the old BC03 TP-AGB prescription. Finally, we compare total stellar mass estimates obtained by integrating resolved mass maps with those obtained with unresolved photometry. In galaxies with evident dust lanes, unresolved estimates may miss up to 40 per cent of the total stellar mass because dusty regions are strongly under-represented in the luminous fluxes.

preprint2008arXiv

A comprehensive Maximum Likelihood analysis of the structural properties of faint Milky Way satellites

We derive the structural parameters of the recently discovered very low luminosity Milky Way satellites through a Maximum Likelihood algorithm applied to SDSS data. For each satellite, even when only a few tens of stars are available down to the SDSS flux limit, the algorithm yields robust estimates and errors for the centroid, position angle, ellipticity, exponential half-light radius and number of member stars. This latter parameter is then used in conjunction with stellar population models of the satellites to derive their absolute magnitudes and stellar masses, accounting for `CMD shot-noise'. We find that faint systems are somewhat more elliptical than initially found and ascribe that to the previous use of smoothed maps which can be dominated by the smoothing kernel. As a result, the faintest half of the Milky Way dwarf galaxies (M_V>-7.5) is significantly (4-sigma) flatter (e=0.47+/-0.03) than its brightest half (M_V<-7.5, e=0.32+/-0.02). From our best models, we also investigate whether the seemingly distorted shape of the satellites, often taken to be a sign of tidal distortion, can be quantified. We find that, except for tentative evidence of distortion in CVnI and UMaII, these can be completely accounted for by Poisson scatter in the sparsely sampled systems. We consider three scenarios that could explain the rather elongated shape of faint satellites: rotation supported systems, stars following the shape of more triaxial dark matter subhalos, or elongation due to tidal interaction with the Milky Way. Although none of these is entirely satisfactory, the last one appears the least problematic, but warrants much deeper observations to track evidence of such tidal interaction.

preprint2008arXiv

A pilot survey of stellar tidal streams in nearby spiral galaxies

Within the hierarchical framework for galaxy formation, merging and tidal interactions are expected to shape large galaxies to this day. While major mergers are quite rare at present, minor mergers and satellite disruptions - which result in stellar streams - should be common, and are indeed seen in both the Milky Way and the Andromeda Galaxy. As a pilot study, we have carried out ultra-deep, wide-field imaging of some spiral galaxies in the Local Volume, which has revealed external views of such stellar tidal streams at unprecedented detail, with data taken at small robotic telescopes (0.1-0.5-meter) that provide exquisite surface brightness sensitivity. The goal of this project is to undertake the first systematic and comprehensive imaging survey of stellar tidal streams, from a sample of ~50 nearby Milky-Way-like spiral galaxies within 15 Mpc, that features a surface brightness sensitivity of ~ 30 mag/arcsec^2 The survey will result in estimates of the incidence, size/geometry and stellar luminosity/mass distribution of such streams. This will not only put our Milky Way and M31 in context but, for the first time, also provide an extensive statistical basis for comparison with state-of-the-art, self-consistent cosmological simulations of this phenomenon.

preprint2008arXiv

Dynamical Evolution of AGN Host Galaxies - Gas In/Out-Flow Rates in 7 NUGA Galaxies

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUGA sample (NUclei of GAlaxies). We calculated gravitational torques from near-IR images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using PdBI) and atomic (HI using the VLA) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the Corotation Resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other non-axisymmmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 solar masses per year and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for 3 galaxies. Non-circular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g. 5 times 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas.

preprint2008arXiv

Identifying stellar streams in the 1st RAVE public data release

We searched for and detected stellar streams or moving groups in the solar neighbourhood, using the data provided by the 1st RAVE public data release. This analysis is based on distances to RAVE stars estimated from a color-magnitude relation that was calibrated on Hipparcos stars. Our final sample consists of 7015 stars selected to be within 500 pc of the Sun and to have distance errors better than 25%. Together with radial velocities from RAVE and proper motions from various data bases, there are estimates for all 6 phase-space coordinates of the stars in the sample. We characterize the orbits of these stars through suitable proxies for their angular momentum and eccentricity, and compare the observed distribution to the expectations from a smooth distribution. On this basis we identify at least four "phase space overdensities" of stars on very similar orbits in the Solar neighbourhood. We estimate the statistical significance of these overdensities by Monte Carlo simulations. Three of them have been identified previously: the Sirius and Hercules moving group and a stream found independently in 2006 by Arifyanto and Fuchs and Helmi et al. In addition, we have found a new stream candidate on a quite radial orbit, suggesting an origin external to the Milky Way's disk. Also, there is evidence for the Arcturus stream and the Hyades-Pleiades moving group in the sample. This analysis, using only a minute fraction of the final RAVE data set, shows the power of this experiment to probe the phase-space substructure of stars around the Sun.

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.

preprint2008arXiv

The Rise of Massive Red Galaxies: the color-magnitude and color-stellar mass diagrams for z < ~2 from the MUltiwavelength Survey by Yale-Chile (MUSYC)

We present the color-magnitude and color-stellar mass diagrams for galaxies with z_phot < ~2, based on a K < 22 (AB) catalog of the Extended Chandra Deep Field South (ECDFS) from the MUltiwavelength Survey by Yale-Chile (MUSYC). Our main sample of 7840 galaxies contains 1297 M_* > 10^11 M_Sol galaxies in the range 0.2 < z_phot < 1.8. We show empirically that this catalog is approximately complete for M_* > 10^11 M_Sol galaxies for z_phot < 1.8. For this mass-limited sample, we show that the locus of the red sequence color-stellar mass relation evolves as Del(u-r) ~ (-0.44+/-0.02) z_phot for z_phot < ~1.2. For z_phot > ~1.3, however, we are no longer able to reliably distinguish red and blue subpopulations from the observed color distribution; we show that this would require much deeper near infrared data. At 1.5 < z_phot <1.8, the comoving number density of M_* > 10^11 M_Sol galaxies is ~50% of the local value, with a red fraction of ~33%. Making a parametric fit to the observed evolution, we find n_tot(z) ~ (1+z_phot)^(-0.52+/-0.12(+/-0.20)). We find stronger evolution in the red fraction: f_red(z) ~ (1+z_phot)^(-1.17+/-0.18(+/-0.21)). Through a series of sensitivity analyses, we show that the most important sources of systematic error are: 1. systematic differences in the analysis of the z~0 and z>>0 samples; 2. systematic effects associated with details of the photometric redshift calculation; and 3. uncertainties in the photometric calibration. With this in mind, we show that our results based on photometric redshifts are consistent with a completely independent analysis which does not require redshift information for individual galaxies. Our results suggest that, at most, 1/5 of local red sequence galaxies with M_* >10^11 M_Sol were already in place at z ~ 2.

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

Comparing Dynamical and Photometric Mass Estimates of Low- and High-Redshift Galaxies: Random and Systematic Uncertainties

We determine the importance of redshift-dependent systematic effects in the determination of stellar masses from broad band spectral energy distributions (SEDs), using high quality kinematic and photometric data of early-type galaxies at z~1 and z~0. We find that photometric masses of z~1 galaxies can be systematically different, by up to a factor of 2, from photometric masses of z~0 galaxies with the same dynamical mass. The magnitude of this bias depends on the choice of stellar population synthesis model and the rest-frame wavelength range used in the fits. The best result, i.e., without significant bias, is obtained when rest-frame optical SEDs are fitted with models from Bruzual&Charlot (2003). When the SEDs are extended to the rest-frame near-IR, a bias is introduced: photometric masses of the z~1 galaxies increase by a factor of 2 relative to the photometric masses of the z~0 galaxies. When we use the Maraston (2005) models, the photometric masses of the z~1 galaxies are low relative to the photometric masses of the z~0 galaxies by a factor of ~1.8. This offset occurs both for fits based on rest-frame optical SEDs, and fits based on rest-frame optical+near-IR SEDs. The results indicate that model uncertainties produce uncertainties as high as a factor of 2.5 in mass estimates from rest-frame near-IR photometry, independent of uncertainties due to unknown star formation histories.

preprint2006arXiv

HST/STIS Spectra of Nuclear Star Clusters in Spiral Galaxies: Dependence of Age and Mass On Hubble Type

(Abridged) We study the nuclear star clusters in spiral galaxies of various Hubble types using spectra obtained with STIS on-board HST. We observed the nuclear clusters in 40 galaxies, selected from two previous HST/WFPC2 imaging surveys. The spectra provide a better separation of cluster light from underlying galaxy light than is possible with ground-based spectra. To infer the star formation history, metallicity and dust extinction, we fit weighted superpositions of single-age stellar population templates to the spectra. The luminosity-weighted age ranges from 10 Myrs to 10 Gyrs. The stellar populations of NCs are generally best fit as a mixture of populations of different ages. This indicates that NCs did not form in a single event, but instead they had additional star formation long after the oldest stars formed. On average, the sample clusters in late-type spirals have a younger luminosity-weighted mean age than those in early-type spirals (log(age/yr) = 8.37+/-0.25 vs. 9.23+/-0.21). The average cluster masses are smaller in late-type spirals than in early-type spirals (log(M/Msun) = 6.25+/-0.21 vs. 7.63+/-0.24), and exceed the masses typical of globular clusters. The cluster mass correlates strongly with both the Hubble type of the host galaxy and the luminosity of its bulge. The latter correlation has the same slope as the well-known correlation between supermassive black hole mass and bulge luminosity. The properties of both nuclear clusters and black holes are therefore intimately connected to the properties of the host galaxy.

preprint2006arXiv

Spectroscopic Identification of Massive Galaxies at z~2.3 with Strongly Suppressed Star Formation

We present first results of a spectroscopic survey targeting K-selected galaxies at z=2.0-2.7 using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength coverage of 1.0-2.5 micron for 26 K-bright galaxies (K<19.7) selected from the MUSYC survey using photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest-frame optical emission lines or the redshifted Balmer/4000 Angstrom break. Twenty galaxies have spectroscopic redshifts in the range 2.0<z<2.7, for which bright emission lines like Halpha and [OIII] fall in atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The median 2 sigma upper limit on the rest-frame equivalent width of Halpha for these nine galaxies is ~10 Angstrom. The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with zero within 1 sigma for the remaining four. Thus, both the Halpha measurements and the independent stellar continuum modeling imply that 45% of our K-selected galaxies are not forming stars intensely. This high fraction of galaxies without detected line emission and low SFRs may imply that the suppression of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy formation models. However, obscured star formation may have been missed, and deep mid-infrared imaging is needed to clarify this situation.

preprint2006arXiv

Stellar populations in the nuclei of late-type spiral galaxies

(Abridged) As part of an ongoing effort to study the stellar nuclei of very late-type, bulge-less spirals, we present results from a high-resolution spectroscopic survey of nine such nuclear star clusters, undertaken with VLT/UVES. We fit the spectra with population synthesis models and measure Lick-type indices to determine mean luminosity-weighted ages, which range from 4.1*10^7 to 1.1*10^10 years and are insensitive to assumed metallicity or internal extinction. The average metallicity of nuclear clusters in late-type spirals is slightly sub-solar (<Z> = 0.015) but shows significant scatter. The fits also show that the nuclear cluster spectra are best described by a mix of several generations of stars. This is supported by the fact that only models with composite stellar populations yield mass-to-light ratios that match those obtained from dynamical measurements. The last star formation episode was on average 34 Myr ago, while all clusters experienced some star formation in the last 100 Myr. We thus conclude that the nuclear clusters undergo repeated episodes of star formation. The robustness with respect to possible contamination from the underlying galaxy disk is demonstrated by comparison to spectra obtained with HST/STIS. Combining these results with those from Walcher et al. (2005), we have thus shown that the stellar nuclei of these bulge-less galaxies are massive and dense star clusters that form stars recurrently until the present day. This unique set of properties is likely due to the central location of these clusters in their host galaxies.

preprint2005arXiv

Discovery of a Fifth Image of the Large Separation Gravitationally Lensed Quasar SDSS J1004+4112

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

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.

preprint2004arXiv

Stellar Populations and Kinematics of Red Galaxies at z>2: Implications for the Formation of Massive Galaxies

We recently identified a substantial population of galaxies at z>2 with red rest-frame optical colors. These distant red galaxies (DRGs) are efficiently selected by the simple observed color criterion J-K>2.3. In this paper we present NIR spectroscopy with Keck/NIRSPEC of six DRGs at 2.4<z<3.2. We detect continuum emission and emission lines of all observed galaxies. Equivalent widths of H alpha are 20-30 Ang, smaller than measured for LBGs and nearby LIRGs, and comparable to normal nearby galaxies. The modest equivalent widths imply that the galaxies either have a decreasing star formation rate, or that they are very dusty. Fitting both the photometry and the H alpha lines, we find continuum extinction A_V=1-2 mag, ages 1-2.5 Gyr, star formation rates 200-400 solar masses/yr, and stellar masses 1-5x10^11 solar masses for models with constant star formation rates. From [NII]/H alpha ratios we infer that the metallicities are high, 1-1.5 x Solar. For four galaxies we can determine line widths from the optical emission lines. The widths are high, ranging from 130-240 km/s, and by combining data for LBGs and DRGs we find significant correlations between linewidth and restframe U-V color, and between linewidth and stellar mass. The latter correlation has a similar slope and offset as the ``baryonic Tully-Fisher relation'' for nearby galaxies. The median dynamical mass is ~2x10^11 solar masses, supporting the high stellar masses inferred from the photometry. We find that the median M/L_V ~ 0.8, a factor of ~5 higher than measured for LBGs. We infer from our small sample that DRGs are dustier, more metal rich, more massive, and have higher ages than z=3 LBGs of the same rest-frame V-band luminosity. Their high M/L ratios imply that they contribute significantly to the stellar mass density at z~2.5. [ABRIDGED]

preprint2002arXiv

Detailed Structural Decomposition of Galaxy Images

We present a two-dimensional (2-D) fitting algorithm (GALFIT) designed to extract structural components from galaxy images, with emphasis on closely modeling light profiles of spatially well-resolved, nearby galaxies observed with the Hubble Space Telescope. Our algorithm improves on previous techniques in two areas, by being able to simultaneously fit a galaxy with an arbitrary number of components, and with optimization in computation speed, suited for working on large galaxy images. We use 2-D models such as the ``Nuker'' law, the Sersic (de Vaucouleurs) profile, an exponential disk, and Gaussian or Moffat functions. The azimuthal shapes are generalized ellipses that can fit disky and boxy components. Many galaxies with complex isophotes, ellipticity changes, and position-angle twists can be modeled accurately in 2-D. When examined in detail, we find that even simple-looking galaxies generally require at least three components to be modeled accurately, rather than the one or two components more often employed. We illustrate this by way of 7 case studies, which include regular and barred spiral galaxies, highly disky lenticular galaxies, and elliptical galaxies displaying various levels of complexities. A useful extension of this algorithm is to accurately extract nuclear point sources in galaxies. We compare 2-D and 1-D extraction techniques on simulated images of galaxies having nuclear slopes with different degrees of cuspiness, and we then illustrate the application of the program to several examples of nearby galaxies with weak nuclei.

preprint2001arXiv

The Ghost of Sagittarius and Lumps in the Halo of the Milky Way

We identify new structures in the halo of the Milky Way Galaxy from positions, colors and magnitudes of five million stars detected in the Sloan Digital Sky Survey. Most of these stars are within 1.26 degrees of the celestial equator. We present color-magnitude diagrams (CMDs) for stars in two previously discovered, tidally disrupted structures. The CMDs and turnoff colors are consistent with those of the Sagittarius dwarf galaxy, as had been predicted. In one direction, we are even able to detect a clump of red stars, similar to that of the Sagittarius dwarf, from stars spread across 110 square degrees of sky. Focusing on stars with the colors of F turnoff objects, we identify at least five additional overdensities of stars. Four of these may be pieces of the same halo structure, which would cover a region of the sky at least 40 degrees in diameter, at a distance of 11 kpc from the Sun (18 kpc from the center of the Galaxy). The turnoff is significantly bluer than that of thick disk stars, and closer to the Galactic plane than a power-law spheroid. We suggest two models to explain this new structure. One possibility is that this new structure could be a new dwarf satellite of the Milky Way, hidden in the Galactic plane, and in the process of being tidally disrupted. The other possibility is that it could be part of a disk-like distribution of stars which is metal-poor, with a scale height of approximately 2 kpc and a scale length of approximately 10 kpc. The fifth overdensity, which is 20 kpc away, is some distance from the Sagittarius dwarf streamer orbit and is not associated with any known structure in the Galactic plane. It is likely that there are many smaller streams of stars in the Galactic halo.

preprint1995arXiv

CONSTRAINTS ON THE OPACITY OF SPIRAL DISKS FROM NEAR-INFRARED OBSERVATIONS

In this paper I review how near-infrared (NIR) observations can constrain the opacity of spiral disks. Basic considerations show that NIR photometry provides a powerful probe of the optical depths in spiral galaxy disks in the regime of interest, where the optical depth in the V-band is near unity. I review the existing opacity constraints from the analysis of dust lanes in edge-on and face-on galaxies. The ``internal extinction correction" in the NIR-Tully-Fisher relation deserves particular attention as the most powerful constraint on the impact of dust on the total luminosity of spiral galaxies. All observations for luminous spirals point towards an effective, face-on optical depth of tau_V=0.5-1.

preprint1995arXiv

THE SHAPES OF DARK HALOS

The shapes of dark halos (e.g. spherical, oblate or triaxial) provide important clues to the formation history of the halos, to their baryon content and to the nature of the dark matter. The halo shapes also influence the dynamics of the galaxy they harbor. This review discusses the current status of theoretical predictions (halo formation in a cosmological context with and without a dissipative component) and of the observational constraints on halo shapes. Halos are neither predicted nor observed to be spherical.