Source author record

Emily C. Cunningham

Emily C. Cunningham 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

13works
3topics
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

13 published item(s)

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

Implications of the Milky Way travel velocity for dynamical mass estimates of the Local Group

The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history including the recent pericentric passage of the LMC, and recent work has found that the MW disk is moving with a lower bound "travel velocity" of $\sim 32~{\rm km}~{\rm s}^{-1}$ with respect to the outer stellar halo. Previous Timing Argument measurements attempt to account for this non-equilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10--12\% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume.

preprint2022arXiv

The Baltimore Oriole's Nest: Cool Winds from the Inner and Outer Parts of a Star-Forming Galaxy at z=1.3

Strong galactic winds are ubiquitous at $z\gtrsim 1$. However, it is not well known where inside galaxies these winds are launched from. We study the cool winds ($\sim 10^4$\,K) in two spatial regions of a massive galaxy at $z=1.3$, which we nickname the "Baltimore Oriole's Nest." The galaxy has a stellar mass of $10^{10.3\pm 0.3} M_\odot$, is located on the star-forming main sequence, and has a morphology indicative of a recent merger. Gas kinematics indicate a dynamically complex system with velocity gradients ranging from 0 to 60 $\mathrm{km}\cdot\mathrm{s}^{-1}$. The two regions studied are: a dust-reddened center (Central region), and a blue arc at 7 kpc from the center (Arc region). We measure the \ion{Fe}{2} and \ion{Mg}{2} absorption line profiles from deep Keck/DEIMOS spectra. Blueshifted wings up to 450 km$\cdot$s$^{-1}$ are found for both regions. The \ion{Fe}{2} column densities of winds are $10^{14.7\pm 0.2}\,\mathrm{cm}^{-2}$ and $10^{14.6\pm 0.2}\,\mathrm{cm}^{-2}$ toward the Central and Arc regions, respectively. Our measurements suggest that the winds are most likely launched from both regions. The winds may be driven by the spatially extended star formation, the surface density of which is around 0.2 $M_\odot\,\mathrm{yr}^{-1}\cdot \mathrm{kpc}^{-2}$ in both regions. The mass outflow rates are estimated to be $4\,M_\odot\,\mathrm{yr}^{-1}$ and $3\,M_\odot\,\mathrm{yr}^{-1}$ for the Central and Arc regions, with uncertainties of one order-of-magnitude or more. Findings of this work and a few previous studies suggest that the cool galactic winds at $z\gtrsim 1$ might be commonly launched from the entire spatial extents of their host galaxies due to extended galaxy star formation.

preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2021arXiv

Reading the CARDs: the Imprint of Accretion History in the Chemical Abundances of the Milky Way's Stellar Halo

In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies' assembly histories by modelling the CARDs of the stellar halos of the Latte galaxies as a linear combination of template CARDs from disrupted dwarfs, with different stellar masses $M_{\star}$ and quenching times $t_{100}$. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of $<10\%$. For the fraction of mass accreted as a function of $t_{100}$, we find residuals of $20-30\%$ for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Though the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos.

preprint2021arXiv

The homogeneity of the star forming environment of the Milky Way disk over time

Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$α$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using measurements from $\sim$27,000 APOGEE stars (R=22,500, SNR$>$200), we build simple local linear models to predict a sample of elements (X = Si, O, Ca, Ti, Ni, Al, Mn, Cr) using (Fe, Mg) abundances alone, as fiducial tracers of supernovae production channels. Given [Fe/H] and [Mg/H], we predict these elements, [X/H], to about double the uncertainty of their measurements. The intrinsic dispersion, after subtracting measurement errors in quadrature is $\approx 0.015-0.04$~dex. The residuals of the prediction (measurement $-$ model) for each element demonstrate that each element has an individual link to birth properties at fixed (Fe, Mg). Residuals from primarily massive-star supernovae (i.e. Si, O, Al) partially correlate with guiding radius. Residuals from primarily supernovae Ia (i.e. Mn, Ni) partially correlate with age. A fraction of the intrinsic scatter that persists at fixed (Fe, Mg), however, after accounting for correlations, does not appear to further discriminate between birth properties that can be traced with present-day measurements. Presumably, this is because the residuals are also, in part, a measure of the typical (in)-homogeneity of the disk's stellar birth environments, previously inferred only using open-cluster systems. Our study implies at fixed birth radius and time, there is a median scatter of $\approx 0.01-0.015$ dex in elements generated in supernovae sources.

preprint2020arXiv

Elemental Abundances in M31: A Comparative Analysis of Iron and Alpha Element Abundances in the Outer Disk, Giant Stellar Stream, and Inner Halo of M31

We measured [Fe/H] and [$α$/Fe] using spectral synthesis of low-resolution stellar spectroscopy for 70 individual red giant branch stars across four fields spanning the outer disk, Giant Stellar Stream (GSS), and inner halo of M31. Fields at M31-centric projected distances of 23 kpc in the halo, 12 kpc in the halo, 22 kpc in the GSS, and 26 kpc in the outer disk are $α$-enhanced, with $\langle$[$α$/Fe]$\rangle$ = 0.43, 0.50, 0.41, and 0.58, respectively. The 23 kpc and 12 kpc halo fields are relatively metal-poor, with $\langle$[Fe/H]$\rangle$ = $-$1.54 and $-$1.30, whereas the 22 kpc GSS and 26 kpc outer disk fields are relatively metal-rich with $\langle$[Fe/H]$\rangle$ = $-$0.84 and $-$0.92, respectively. For fields with substructure, we separated the stellar populations into kinematically hot stellar halo components and kinematically cold components. We did not find any evidence of an [$α$/Fe] gradient along the high surface brightness core of the GSS between $\sim$17$-$22 kpc. However, we found tentative suggestions of a negative [$α$/Fe] gradient in the stellar halo, which may indicate that different progenitor(s) or formation mechanisms contributed to the build up of the inner versus outer halo. Additionally, the [$α$/Fe] distribution of the metal-rich ([Fe/H] $>$ $-$1.5), smooth inner stellar halo (r$_{\rm{proj}}$ $\lesssim$ 26 kpc) is inconsistent with having formed from the disruption of progenitor(s) similar to present-day M31 satellite galaxies. The 26 kpc outer disk is most likely associated with the extended disk of M31, where its high $α$-enhancement provides support for an episode of rapid star formation in M31's disk, possibly induced by a major merger.

preprint2020arXiv

Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.

preprint2020arXiv

The power of co-ordinate transformations in dynamical interpretations of Galactic structure

$Gaia$ DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the $Gaia$ data \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-ϕ$ and the $v_{\mathrm{R}}-v_ϕ$ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in $Gaia$ and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure.

preprint2016arXiv

Isotropic at the Break? 3D Kinematics of Milky Way Halo Stars in the Foreground of M31

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $β$ using 3D kinematics outside of the solar neighborhood. We find $β=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $β$ measurements by 2$σ$ ($β_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $β=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1σ$. The potential decrease in $β$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

preprint2016arXiv

Runaway Dwarf Carbon Stars as Candidate Supernova Ejecta

The dwarf carbon (dC) star SDSS J112801.67+004034.6 has an unusually high radial velocity, 531$\pm 4$ km s$^{-1}$. We present proper motion and new spectroscopic observations which imply a large Galactic rest frame velocity, 425$\pm 9$ km s$^{-1}$. Several other SDSS dC stars are also inferred to have very high galactocentric velocities, again each based on both high heliocentric radial velocity and also confidently detected proper motions. Extreme velocities and the presence of $C_2$ bands in the spectra of dwarf stars are both rare. Passage near the Galactic center can accelerate stars to such extreme velocities, but the large orbital angular momentum of SDSS J1128 precludes this explanation. Ejection from a supernova in a binary system or disruption of a binary by other stars are possibilities, particularly as dC stars are thought to obtain their photospheric $C_2$ via mass transfer from an evolved companion.

preprint2015arXiv

The Most Massive Ultra-Compact Dwarf Galaxy in the Virgo Cluster

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$α$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

preprint2013arXiv

A Search for RR Lyrae Stars in Segue 2 and Segue 3

We present an extensive search for RR Lyrae stars in and around the ultra-faint Milky Way companions Segue 2 and Segue 3. The former (M_V = -2.5, Belokurov et al. 2009) appears to be an extremely faint dwarf galaxy companion of the Milky Way. The latter (M_V = 0.0, Fadely et al. 2011) is among the faintest star clusters known. We use B and V band time-series imaging obtained at the WIYN 0.9 meter telescope at Kitt Peak National Observatory to search for RR Lyrae in these objects. In our Segue 2 observations, we present a previously unknown fundamental mode (RRab) RR Lyrae star with a period of P_ab = 0.748 days. With this measurement, we revisit the inverse correlation between <P_ab> and <[Fe/H]> established in the literature for Milky Way dwarf galaxies and their RR Lyrae. In this context, the long period of Segue 2's RRab star as well as the known significant spread in metallicity in this dwarf galaxy are consistent with the observed trend in <P_ab> and <[Fe/H]>. We derive the first robust distance to Segue 2, using both its RRab star and spectroscopically confirmed blue horizontal branch stars. Using [Fe/H] = -2.16 and -2.44 dex, we find d_RRL = 36.6 +2.5/-2.4 and 37.7 +/- 2.7 kpc; assuming [Fe/H] = -2.257 dex, we find d_BHB = 34.4 +/- 2.6 kpc. Although no RR Lyrae were present in the Segue 3 field, we found a candidate eclipsing binary star system.