Source author record

N. Wyn Evans

N. Wyn Evans 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

41works
11topics
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

41 published item(s)

preprint2022arXiv

Astrometric Identification of Nearby Binary Stars I: Predicted Astrometric Signals

We examine the capacity to identify binary systems from astrometric errors and deviations alone. Until the release of the fourth Gaia data release we lack the full astrometric time series that the satellite records, but as we show can still infer the presence of binaries from the best fit models, and their error, already available. We generate a broad catalog of simulated binary systems within 100 pc, and examine synthetic observations matching the Gaia survey's scanning law and astrometric data processing routine. We show how the Unit Weight Error (UWE) and Proper Motion Anomaly (PMA) vary as a function of period, and the properties of the binary. Both UWE and PMA peak for systems with a binary period close to the time baseline of the survey. Thus UWE can be expected to increase or remain roughly constant as we observe the same system over a longer baseline, and we suggest $UWE_{eDR3}>1.25$ and $ΔUWE/UWE_{eDR3}>-0.25$ as criteria to select astrometric binaries. For stellar binaries we find detectable significant astrometric deviations for 80-90\% of our simulated systems in a period range from months to decades. We confirm that for systems with periods less than the survey's baseline the observed $UWE$ scales $\propto \ \varpi$ (parallax), $a$ (semi-major axis) and $Δ=\frac{|q-l|}{(1+q)(1+l)}$ where $q$ and $l$ are the mass and light ratio respectively, with a modest dependence on viewing angle. For longer periods the signal is suppressed by a factor of roughly $\propto P^{-2}$ (period). PMA is largest in orbits with slightly longer periods but obeys the same approximate scaling relationships.

preprint2022arXiv

Astrometric identification of nearby binary stars II: Astrometric binaries in the Gaia Catalogue of Nearby Stars

We examine the capacity to identify binary systems from astrometric deviations alone. We apply our analysis to the Gaia eDR3 and DR2 data, specifically the Gaia Catalogue of Nearby Stars. We show we must renormalize (R)UWE over the local volume to avoid biasing local observations, giving a Local Unit Weight Error (LUWE). We use the simple criterion of LUWE>2, along with a handful of quality cuts to remove likely contaminants, to identify unresolved binary candidates. We identify 22,699 binary candidates within 100 pc of the Sun (just under 10% of sources in this volume). We find an astrometric binary candidate fraction of around 20% for giant stars, 10% on the Main Sequence and lower than 1% for White Dwarfs. We also look for Variability Induced Movers, by computing the correlation between photometric variability and astrometric noise -- and show that VIMs may dominate the binary population of sub-Solar mass MS stars. We discuss the possibility and limitations of identifying non-luminous massive companions from astrometry alone, but find that our method is insensitive to these. Finally, we compare the astrometric deviations of MS binaries to the simulated sample from paper I, which show excellent agreement, and compute the astrometric candidate binary fraction as a function of absolute magnitude.

preprint2022arXiv

Catalogues of Active Galactic Nuclei From Gaia and unWISE Data

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://zenodo.org/record/6837642#.YtEgjy8RrzI.

preprint2022arXiv

Dark lenses through the dust: parallax microlensing events in the VVV

We use near-infrared photometry and astrometry from the VISTA Variables in the Via Lactea (VVV) survey to analyse microlensing events containing annual microlensing parallax information. These events are located in highly extincted and low-latitude regions of the Galactic bulge typically off-limits to optical microlensing surveys. We fit a catalog of $1959$ events previously found in the VVV and extract $21$ microlensing parallax candidates. The fitting is done using nested sampling to automatically characterise the multi-modal and degenerate posterior distributions of the annual microlensing parallax signal. We compute the probability density in lens mass-distance using the source proper motion and a Galactic model of disc and bulge deflectors. By comparing the expected flux from a main sequence lens to the baseline magnitude and blending parameter, we identify 4 candidates which have probability $> 50$% that the lens is dark. The strongest candidate corresponds to a nearby ($\approx0.78$ kpc), medium-mass ($1.46^{+1.13}_{-0.71} \ M_{\odot}$) dark remnant as lens. In the next strongest, the lens is located at heliocentric distance $\approx5.3$ kpc. It is a dark remnant with a mass of $1.63^{+1.15}_{-0.70} \ M_{\odot}$. Both of those candidates are most likely neutron stars, though possibly high-mass white dwarfs. The last two events may also be caused by dark remnants, though we are unable to rule out other possibilities because of limitations in the data.

preprint2022arXiv

Information content of BP/RP spectra in Gaia DR3

Gaia Data Release 3 has provided the astronomical community with the largest stellar spectroscopic survey to date ($>$ 220 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar atmospheric parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fisher information matrices to probe the resolution limit of stellar parameter measurements using BP/RP spectra. The best-case scenario uncertainties that this analysis provides are then used to produce a mock-observed stellar population in order to evaluate the false positive rate (FPR) of identifying extremely metal-poor (EMP) stars. We conclude that the community will be able to confidently identify metal-poor stars at magnitudes brighter than $G = 16$ using BP/RP spectra. At fainter magnitudes true detections will start to be overwhelmed by false positives. When adopting the commonly-used $G < 14$ limit for metal-poor star searches, we find a FPR for the low-metallicity regimes [Fe/H] $<$ -2, -2.5 and -3 of just 14$\%$, 33$\%$ and 56$\%$ respectively, offering the potential for significant improvements on previous targeting campaigns. Additionally, we explore the chemical sensitivity obtainable directly from BP/RP spectra for Carbon and $α$-elements. We find an absolute Carbon abundance uncertainty of $σ_{A(C)} < 1$ dex for Carbon-enriched metal-poor (CEMP) stars, indicating the potential to identify a CEMP stellar population for follow-up confirmation with higher resolution spectroscopy. Finally, we find that large uncertainties in $α$-element abundance measurements using BP/RP spectra means that efficiently obtaining these abundances will be challenging.

preprint2022arXiv

Mass Segregation in the Hyades Cluster

Using the Gaia colour-magnitude diagram, we assign masses to a catalogue of 979 confirmed members of the Hyades cluster and tails. By fitting the cumulative mass profile, stars within the tidal radius have a Plummer-like profile with half-mass radius $r_{\rm h}$ of 5.75 pc. The tails are extended with $r_{\rm h} = 69.35$ pc and fall off more slowly than Plummer with density proportional to distance$^{-1.36}$. The cluster stars are separated into two groups at BP-RP $=2$ or $0.56 M_\odot$ to give a high mass (${\bar M} = 0.95 M_\odot$) and a low mass (${\bar M} = 0.32 M_\odot$) population. We show that: (i) the high mass population has a half-mass radius $r_{\rm h}$ of 4.88 pc, whilst the low mass population has $r_{\rm h} = 8.10$ pc; (ii) despite the differences in spatial extent, the kinematics and binarity properties of the high and low mass populations are similar. They have isotropic velocity ellipsoids with mean 1d velocity dispersions $σ$ of 0.427 and 0.415 km s$^{-1}$ respectively. The dynamical state of the Hyades is far from energy equipartition ($σ\propto {\bar M}^{-1/2}$). We identify a new mass segregation instability for clusters with escape speed $V$. Populations with $V/σ\lesssim 2\sqrt{2}$ can never attain thermal equilibrium and equipartition. This regime encompasses many Galactic open and globular clusters. For the Hyades, there must be an outward energy flux of at least $9.5 \times 10^{-4} M_\odot\,{\rm km^2\, s^{-2} Myr^{-1}}$ to maintain its current configuration. The present mass loss of $0.26 M_\odot {\rm Myr}^{-1}$ due to tidal stripping by itself implies a substantial energy flow beyond the required magnitude.

preprint2022arXiv

The impact of a massive Sagittarius dSph on GD-1-like streams

We investigate the effects of a massive ($\gtrsim4\times10^{10}M_\odot$) Sagittarius dwarf spheroidal galaxy (Sgr) on stellar streams using test particle simulations in a realistic Milky Way potential. We find that Sgr can easily disrupt streams formed more than $\sim3$ Gyr ago, while stars stripped more recently are generally unaffected. In certain realizations, Sgr is able to produce asymmetry between the leading and trailing tails of Pal 5, qualitatively similar to observations. Using data from the Gaia space telescope and elsewhere, we fit models to the GD-1 stream in the presence of a Sgr with various initial masses. While the best-fitting models do show perturbations resulting from interactions with Sgr, we find that the level of disruption is not significantly greater than in the observed stream. To investigate the general effects of Sgr on a population of streams, we generate 1000 mock streams on GD-1-like orbits with randomized orientations. Some streams show clear evidence of disruption, becoming folded on the sky or developing asymmetry betweeen their two tails. However, many survive unaffected and the peak surface brightness of stars is decreased by no more than $\sim0.3$ mag/arcsec$^2$ on average. We conclude that Sgr having an initial mass of $\gtrsim4\times10^{10}M_\odot$ is compatible with the survival and detection of streams formed more than 3 Gyr ago.

preprint2022arXiv

The Local Group Mass in the light of Gaia

High accuracy proper motions (PMs) of M31 and other Local Group satellites have now been provided by the {\it Gaia} satellite. We revisit the Timing Argument to compute the total mass $M$ of the Local Group from the orbit of the Milky Way and M31, allowing for the Cosmological Constant. We rectify for a systematic effect caused by the presence of the Large Magellanic Cloud (LMC). The interaction of the LMC with the Milky Way induces a motion towards the LMC. This contribution to the measured velocity of approach of the Milky Way and M31 must be removed. We allow for cosmic bias and scatter by extracting correction factors tailored to the accretion history of the Local Group. The distribution of correction factors is centered around $0.63$ with a scatter $\pm 0.2$, indicating that the Timing Argument significantly overestimates the true mass. Adjusting for all these effects, the estimated mass of the Local Group is $ M = 3.4^{+1.4}_{-1.1} \times 10^{12} M_{\odot}$ (68 % CL) when using the M31 tangential velocity $ 82^{+38}_{-35}$ km/s. Lower tangential velocity models with $59^{+42}_{-38}$ km/s (derived from the same PM data with a flat prior on the tangential velocity) lead to an estimated mass of $ M = 3.1^{+1.3}_{-1.0} \times 10^{12} M_{\odot}$ (68 % CL). By making an inventory of the total mass associated with the 4 most substantial LG members (the Milky Way, M31, M33 and the LMC), we estimate the known mass is in the range $3.7^{+0.5}_{-0.5} \times 10^{12} \, M_{\odot}$.

preprint2021arXiv

Charting galactic accelerations II: how to 'learn' accelerations in the solar neighbourhood

Gravitational acceleration fields can be deduced from the collisionless Boltzmann equation, once the distribution function is known. This can be constructed via the method of normalizing flows from datasets of the positions and velocities of stars. Here, we consider application of this technique to the solar neighbourhood. We construct mock data from a linear superposition of multiple `quasi-isothermal' distribution functions, representing stellar populations in the equilibrium Milky Way disc. We show that given a mock dataset comprising a million stars within 1 kpc of the Sun, the underlying acceleration field can be measured with excellent, sub-percent level accuracy, even in the face of realistic errors and missing line-of-sight velocities. The effects of disequilibrium can lead to bias in the inferred acceleration field. This can be diagnosed by the presence of a phase space spiral, which can be extracted simply and cleanly from the learned distribution function. We carry out a comparison with two other popular methods of finding the local acceleration field (Jeans analysis and 1D distribution function fitting). We show our method most accurately measures accelerations from a given mock dataset, particularly in the presence of disequilibria.

preprint2021arXiv

Discovering Strongly-lensed QSOs From Unresolved Light Curves

We present a new method of discovering galaxy-scale, strongly-lensed QSO systems from unresolved light curves using the autocorrelation function. The method is tested on five rungs of simulated light curves from the Time Delay Challenge 1 that were designed to match the light-curve qualities from existing, ongoing, and forthcoming time-domain surveys such as the Medium Deep Survey of the Panoramic Survey Telescope And Rapid Response System 1, the Zwicky Transient Facility, and the Rubin Observatory Legacy Survey of Space and Time. Among simulated lens systems for which time delays can be successfully measured by current best algorithms, our method achieves an overall true positive rate of 28--58% for doubly-imaged QSOs (doubles) and 36--60% for quadruply-imaged QSOs (quads) while maintains $\lesssim$10% false positive rates. We also apply the method to observed light curves of 22 known strongly-lensed QSOs, and recover 20% of doubles and 25% of quads. The tests demonstrate the capability of our method for discovering strongly-lensed QSOs from major time domain surveys. The performance of our method can be further improved by analysing multi-filter light curves and supplementing with morphological, colour, and/or astrometric constraints. More importantly, our method is particularly useful for discovering small-separation strongly-lensed QSOs, complementary to traditional imaging-based methods.

preprint2021arXiv

Elevated r-process enrichment in Gaia Sausage and Sequoia

The Gaia Sausage (GS) and the Sequoia represent the major accretion events that formed the stellar halo of the Milky Way. A detailed chemical study of these main building blocks provides a pristine view of the early steps of the Galaxy's assembly. We present the results of the analysis of the UVES high-resolution spectroscopic observations at the 8.2m VLT of 9 Sausage/Sequoia members selected kinematically using Gaia DR2. We season this set of measurements with archival data from Nissen & Schuster (2011) and GALAH DR3 (2020). Here, we focus on the neutron-capture process by analysing Sr, Y, Ba and Eu behavior. We detect clear enhancement in Eu abundance ([Eu/Fe]~0.6-0.7) indicative of large prevalence of r-process in the stellar n-capture makeup. We are also able to trace the evolution of the heavy element production across a wide range of metallicity. The barium to europium changes from a tight, flat sequence with [Ba/Eu]=-0.7 reflecting dominant contribution from exploding massive stars, to a clear upturn at higher iron abundances, betraying the onset of contamination from asymptotic giant branch (AGB) ejecta. Additionally, we discover two clear sequences in [Fe/H]-[Ba/Fe] plane likely caused by distinct levels of s-process pollution and mixing within the GS progenitor.

preprint2021arXiv

The Photo-Astrometric Vertical Tracer Density of the Milky Way II: Results from Gaia

We use Gaia photometry and astrometry to estimate the vertical spatial structure of the Milky Way at the Solar radius, formally accounting for sample incompleteness (the selection function) and parallax measurement uncertainty. Our results show impressive precision demonstrating the power of the Gaia data. However, systematic errors dominate the parameter value uncertainties. We thoroughly test and quantify the impacts of all systematic uncertainties. The vertical tracer density is modelled as a sum of two exponential profiles for the thin and thick discs, together with a spherically symmetric power-law for the stellar halo. We constrain the thin disc scale height as ${h_\mathrm{Tn}=260 \pm 3\, (\mathrm{stat}) \pm 26\,\mathrm{pc}\, (\mathrm{sys})}$ and thick disc ${h_\mathrm{Tk}=693 \pm 7 \,(\mathrm{stat}) \pm 121\,\mathrm{pc}\, (\mathrm{sys})}$. For the halo, we obtain a power law profile with $n_\mathrm{H}=3.543\pm0.023 \,(\mathrm{stat}) \pm0.259\, (\mathrm{sys})$. We infer a local stellar mass density for non-compact object stars of ${ρ_\mathrm{local}^* = 3.66\pm0.03\,(\mathrm{stat})\pm0.52 \times10^{-2}\,\mathrm{M}_\odot/\mathrm{pc}^3\,(\mathrm{sys})}$ and surface density of ${Σ_\mathrm{local}^* = 23.17\pm0.08\,(\mathrm{stat})\pm2.43\,\mathrm{M}_\odot/\mathrm{pc}^2\,(\mathrm{sys})}$. We find asymmetries above and below the disc with longer disc scale heights in the north but a flatter halo in the south at the $\lesssim 10$ per cent level.

preprint2020arXiv

Binary deviations from single object astrometry

Most binaries are undetected. Astrometric reductions of a system using the assumption that the object moves like a single point mass can be biased by unresolved binary stars. The discrepancy between the centre of mass of the system (which moves like a point mass) and the centre of light (which is what we observe) introduces additional motion. We explore the extent to which binary systems affect single object models fit to astrometric data. This tells us how observations are diluted by binaries and which systems cause the largest discrepancies - but also allows us to make inferences about the binarity of populations based on observed astrometric error. By examining a sample of mock observations, we show that binaries with periods close to one year can mimic parallax and thus bias distance measurements, whilst long period binaries can introduce significant apparent proper motion. Whilst these changes can soak up some of the error introduced by the binary, the total deviation from the best fitting model can be translated into a lower limit on the on-sky separation of the pair. Throughout we link these predictions to data from the Gaia satellite, whilst leaving the conclusions generalizable to other surveys.

preprint2020arXiv

Dark Shards: velocity substructure from Gaia and direct searches for dark matter

Data from the Gaia satellite show that the solar neighbourhood of the Milky Way's stellar halo is imprinted with substructure from several accretion events. Evidence of these events is found in "the Shards", stars clustering with high significance in both action space and metallicity. Stars in the Shards share a common origin, likely as ancient satellite galaxies of the Milky Way, so will be embedded in dark matter (DM) counterparts. These "Dark Shards" contain two substantial streams (S1 and S2), as well as several retrograde, prograde and lower energy substructures. The retrograde stream S1 has a very high Earth-frame speed of $\sim 550$ km s$^{-1}$ while S2 moves on a prograde, but highly polar orbit and enhances the peak of the speed distribution at around $300$ km s$^{-1}$. The presence of the Dark Shards locally leads to modifications of many to the fundamental properties of experimental DM signals. The~S2 stream in particular gives rise to an array of effects in searches for axions and in the time dependence of nuclear recoils: shifting the peak day, inducing nonsinusoidal distortions, and increasing the importance of the gravitational focusing of DM by the Sun. Dark Shards additionally bring new features for directional signals, while also enhancing the DM flux towards Cygnus.

preprint2020arXiv

Kinematic modelling of clusters with Gaia: the Death Throes of the Hyades

The precision of the Gaia data offers a unique opportunity to study the internal velocity field of star clusters. We develop and validate a forward-modelling method for the internal motions of stars in a cluster. The model allows an anisotropic velocity dispersion matrix and linear velocity gradient describing rotation and shear, combines radial velocities available for a subset of stars, and accounts for contamination from background sources via a mixture model. We apply the method to Gaia DR2 data of the Hyades cluster and its tidal tails, dividing and comparing the kinematics of stars within and beyond 10 pc, which is roughly the tidal radius of the cluster. While the velocity dispersion for the cluster is nearly isotropic, the velocity ellipsoid for the tails is clearly elongated with the major axis pointing towards the Galactic centre. We find positive and negative expansion at $\approx 2σ$ significance in Galactic azimuthal and vertical direction for the cluster but no rotation. The tidal tails are stretching in a direction tilted from the Galactic centre while equally contracting as the cluster in Galactic vertical direction. The tails have a shear (A) of $16.90\pm 0.92$ m/s/pc and a vorticity (B) of $-6.48\pm1.15$ m/s/pc, values distinct from the local Oort constants.

preprint2020arXiv

The Halo's Ancient Metal-Rich Progenitor Revealed with BHB Stars

Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

preprint2020arXiv

Unresolved stellar companions with Gaia DR2 astrometry

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

preprint2019arXiv

Dynamical Relics of the Ancient Galactic Halo

We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing substructures in the nearby halo, such as the $Gaia$ Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric $Gaia$ Sausage groups contain representatives both of debris from the satellite itself (which is $α$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $α$-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the $Gaia$ Sausage host two relatively metal-rich $r$-II stars and six CEMP stars in different sub-classes, consistent with the idea that the $Gaia$ Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly $r$-process-enhanced stars with [Fe/H] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers.

preprint2019arXiv

Spectroscopy of the Young Stellar Association Price-Whelan 1: Origin in the Magellanic Leading Arm and Constraints on the Milky Way Hot Halo

We report spectroscopic measurements of stars in the recently discovered young stellar association Price-Whelan 1 (PW 1), which was found in the vicinity of the Leading Arm (LA) of the Magellanic Stream. We obtained Magellan+MIKE high-resolution spectra of the 28 brightest stars in PW 1 and used The Cannon to determine their stellar parameters. We find that the mean metallicity of PW 1 is [Fe/H]=-1.23 with a small scatter of 0.06 dex and the mean radial velocity is Vhelio=276.7 km/s with a dispersion of 11.0 km/s. Our results are consistent in Teff, logg, and [Fe/H] with the young and metal-poor characteristics (116 Myr and [Fe/H]=-1.1) determined for PW 1 from our discovery paper. We find a strong correlation between the spatial pattern of the PW 1 stars and the LA II gas with an offset of -10.15 deg in L_MS and +1.55 deg in B_MS. The similarity in metallicity, velocity, and spatial patterns indicates that PW 1 likely originated in LA II. We find that the spatial and kinematic separation between LA II and PW 1 can be explained by ram pressure from Milky Way gas. Using orbit integrations that account for the LMC and MW halo and outer disk gas, we constrain the halo gas density at the orbital pericenter of PW 1 to be n_halo (17 kpc) = 2.7 (3) x 10^-3 atoms/cm^3 and the disk gas density at the midplane at 20 kpc to be n_disk (20 kpc,0) = 6.0 (1.8) x 10^-2 atoms/cm^3. We, therefore, conclude that PW 1 formed from the LA II of the Magellanic Stream, making it a powerful constraint on the Milky Way-Magellanic interaction.

preprint2019arXiv

The Biggest Splash

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation.

preprint2016arXiv

Indirect Dark Matter Detection for Flattened Dwarf Galaxies

Gamma-ray experiments seeking to detect evidence of dark matter annihilation in dwarf spheroidal galaxies require knowledge of the distribution of dark matter within these systems. We analyze the effects of flattening on the annihilation (J) and decay (D) factors of dwarf spheroidal galaxies with both analytic and numerical methods. Flattening has two consequences: first, there is a geometric effect as the squeezing (or stretching) of the dark matter distribution enhances (or diminishes) the J-factor; second, the line of sight velocity dispersion of stars must hold up the flattened baryonic component in the flattened dark matter halo. We provide analytic formulae and a simple numerical approach to estimate the correction to the J- and D-factors required over simple spherical modeling. The formulae are validated with a series of equilibrium models of flattened stellar distributions embedded in flattened dark-matter distributions. We compute corrections to the J- and D-factors for the Milky Way dwarf spheroidal galaxies under the assumption that they are all prolate or all oblate and find that the hierarchy of J-factors for the dwarf spheroidals is slightly altered (typical correction factors for an ellipticity of $0.4$ are $0.75$ for the oblate case and $1.6$ for the prolate case). We demonstrate that spherical estimates of the D-factors are very insensitive to the flattening and introduce uncertainties significantly less than the uncertainties in the D-factors from the other observables for all the dwarf spheroidals (for example, ${}^{+10\mathrm{\,percent}}_{-3\mathrm{\,percent}}$ for a typical ellipticity of $0.4$). We conclude by investigating the spread in correction factors produced by triaxial figures and provide uncertainties in the J-factors for the dwarf spheroidals using different physically-motivated assumptions for their intrinsic shape and axis alignments. (abridged)

preprint2016arXiv

Mass estimators for flattened dispersion-supported galaxies

We investigate the reliability of mass estimators based on the observable velocity dispersion and half-light radius $R_\mathrm{h}$ for dispersion-supported galaxies. We show how to extend them to flattened systems and provide simple formulae for the mass within an ellipsoid under the assumption the dark matter density and the stellar density are stratified on the same self-similar ellipsoids. We demonstrate explicitly that the spherical mass estimators (Walker et al. 2009, Wolf et al. 2010) give accurate values for the mass within the half-light ellipsoid, provided $R_\mathrm{h}$ is replaced by its `circularized' analogue $R_\mathrm{h}\sqrt{1-ε}$. We provide a mathematical justification for this surprisingly simple and effective workaround. It means, for example, that the mass-to-light ratios are valid not just when the light and dark matter are spherically distributed, but also when they are flattened on ellipsoids of the same constant shape.

preprint2015arXiv

Beasts of the Southern Wild: Discovery of nine Ultra Faint satellites in the vicinity of the Magellanic Clouds

We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be established with the DES data alone. Based on the morphological properties, three of the new satellites are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue~1 satellite and can not be classified straightforwardly without follow-up spectroscopic observations. The satellites we have discovered cluster around the LMC and the SMC. We show that such spatial distribution is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new satellites must have been associated with the Magellanic Clouds in the past.

preprint2015arXiv

Indication of Gamma-ray Emission from the Newly Discovered Dwarf Galaxy Reticulum II

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

preprint2015arXiv

Magellan/M2FS Spectroscopy of the Reticulum 2 Dwarf Spheroidal Galaxy

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 182 stellar targets along the line of sight to the newly-discovered `ultrafaint' object Reticulum 2 (Ret 2). For 38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity ($v_{\rm los}$, median random error $δ_{v_{\rm los}}=1.3$ km s$^{-1}$), effective temperature ($T_{\rm eff}$, $δ_{T_{\rm eff}}=464$ K), surface gravity ($\log g$, $δ_{\rm logg}=0.54$ dex) and iron abundance ([Fe/H], $δ_{\mathrm{[Fe/H]}}=0.45$ dex). We use these results to confirm 18 stars as members of Ret 2. From the member sample we estimate a velocity dispersion of $σ_{v_{\rm los}}=3.6_{-0.6}^{+0.9}$ km s$^{-1}$ about a mean of $\langle v_{\rm los}\rangle =64.8_{-1.0}^{+1.1}$ km s$^{-1}$ in the solar rest frame ($\sim -90.9$ km s$^{-1}$ in the Galactic rest frame), and a metallicity dispersion of $σ_{\rm [Fe/H]}=0.50_{-0.13}^{+0.17}$ dex about a mean of $\langle \mathrm{[Fe/H]} \rangle =-2.67_{-0.34}^{+0.34}$. These estimates marginalize over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on chemodynamical scaling relations followed by the Milky Way's dwarf-galactic satellites. Under assumptions of dynamic equilibrium and negligible contamination from binary stars---both of which must be checked with deeper imaging and repeat spectroscopic observations---the estimated velocity dispersion suggests a dynamical mass of $M(R_{\rm h})\approx 5R_{\rm h}σ_{v_{\rm los}}^2/(2G)=2.4_{-0.8}^{+1.3}\times 10^5$ $M_{\odot}$ enclosed within projected halflight radius $R_{\rm h}\sim 32$ pc, with mass-to-light ratio $\approx 2M(R_{\rm h})/L_{\rm V}=462_{-157}^{+264}$ in solar units.

preprint2015arXiv

Resonant Orbits and the High Velocity Peaks Towards the Bulge

We extract the resonant orbits from an N-body bar that is a good representation of the Milky Way, using the method recently introduced by Molloy et al. (2015). By decomposing the bar into its constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards the Galactic centre. The resonant orbits are shown to have distinct kinematic features and may be used to explain the cold velocity peak seen in the APOGEE commissioning data (Nidever at al., 2012). We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit family and that stars on other higher order resonances can contribute to the peaks. The locations of the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 10 < theta_bar < 25 (deg). However, some important questions about the nature of the peaks remain, such as their apparent absence in other surveys of the Bulge and the deviations from symmetry between equivalent fields in the north and south. We show that the absence of a peak in surveys at higher latitudes is likely due to the combination of a less prominent peak and a lower number density of bar supporting orbits at these latitudes.

preprint2015arXiv

Self-consistent triaxial models

We present self-consistent triaxial stellar systems that have analytic distribution functions (DFs) expressed in terms of the actions. These provide triaxial density profiles with cores or cusps at the centre. They are the first self-consistent triaxial models with analytic DFs suitable for modelling giant ellipticals and dark haloes. Specifically, we study triaxial models that reproduce the Hernquist profile from Williams & Evans (2015), as well as flattened isochrones of the form proposed by Binney (2014). We explore the kinematics and orbital structure of these models in some detail. The models typically become more radially anisotropic on moving outwards, have velocity ellipsoids aligned in Cartesian coordinates in the centre and aligned in spherical polar coordinates in the outer parts. In projection, the ellipticity of the isophotes and the position angle of the major axis of our models generally changes with radius. So, a natural application is to elliptical galaxies that exhibit isophote twisting. As triaxial Stäckel models do not show isophote twists, our DFs are the first to generate mass density distributions that do exhibit this phenomenon, typically with a gradient of $\approx 10^\circ$/effective radius, which is comparable to the data. Triaxiality is a natural consequence of models that are susceptible to the radial orbit instability. We show how a family of spherical models with anisotropy profiles that transition from isotropic at the centre to radially anisotropic becomes unstable when the outer anisotropy is made sufficiently radial. Models with a larger outer anisotropy can be constructed but are found to be triaxial. We argue that the onset of the radial orbit instability can be identified with the transition point when adiabatic relaxation yields strongly triaxial rather than weakly spherical endpoints.

preprint2013arXiv

Sagittarius stream 3-d kinematics from SDSS Stripe 82

Using multi-epoch observations of the Stripe 82 region done by Sloan Digital Sky Survey, we measure precise statistical proper motions of the stars in the Sagittarius stellar stream. The multi-band photometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the proper motion precision to ~0.1 mas/yr. We measure separately the proper motion of a photometrically selected sample of the main sequence turn-off stars, as well as of a spectroscopically selected Sgr giants. The data allow us to determine the proper motion separately for the two Sgr streams in the South found in Koposov et al.(2012). Together with the precise velocities from SDSS, our proper motion provide exquisite constraints of the 3-D motions of the stars in the Sgr streams.

preprint2013arXiv

The Great Escape III: Placing post-main-sequence evolution of planetary and binary systems in a Galactic context

Our improving understanding of the life cycle of planetary systems prompts investigations of the role of the Galactic environment before, during and after Asymptotic Giant Branch (AGB) stellar evolution. Here, we investigate the interplay between stellar mass loss, Galactic tidal perturbations, and stellar flybys for evolving stars which host one planet, smaller body or stellar binary companion and reside in the Milky Way's bulge or disc. We find that the potential evolutionary pathways from a main sequence (MS) to a white dwarf (WD) planetary system are a strong function of Galactocentric distance only with respect to the prevalence of stellar flybys. Planetary ejection and collision with the parent star should be more common towards the bulge. At a given location anywhere in the Galaxy, if the mass loss is adiabatic, then the secondary is likely to avoid close flybys during AGB evolution, and cannot eventually escape the resulting WD because of Galactic tides alone. Partly because AGB mass loss will shrink a planetary system's Hill ellipsoid axes by about 20 to 40 per cent, Oort clouds orbiting WDs are likely to be more depleted and dynamically excited than on the MS.

preprint2012arXiv

A census of orbital properties of the M31 satellites

We present an analysis of the dynamics of the M31 satellite system. Proper motion data exist for only two of the M31 satellites. We account for this incompleteness in velocity data by a statistical analysis using a combination of the timing argument and phase-space distribution functions. The bulk of the M31 satellites are well fit by these models and we offer a table of orbital properties, including period, eccentricity and semi-major axis. This enables us to search for evidence of group infall based on orbital similarity rather than propinquity on the sky. Our results favour an association between Cass II and NGC 185, as the orbital parameters are in close agreement, but not for NGC 185 and NGC 147, which have often been associated in the past. Other possible satellite groupings include the pair And I and And XVII; the pair And IX and And X; and the triple And V, And XXV and NGC 147. And XXII has been claimed as a satellite of M33; we find that they are not moving independently along the same orbit, but cannot determine whether they are orbiting each other or are unrelated. Two satellites, And XII and And XIV, have high line-of-sight velocities, consistent with very recent infall from the edge of the Local Group. They are not well described by our underlying smooth phase space distribution function, and are reanalysed without priors on their orbital parameters. For And XIV, multiple pericentric passages are possible and improved distance information is needed to draw further conclusions. For And XII, orbits which assume at least one pericentric passage can be ruled out and it must be on its first infall into the M31 system.

preprint2012arXiv

Exoplanets Beyond the Solar Neighbourhood: Galactic Tidal Perturbations

The majority of Milky Way extrasolar planets likely reside within a few kpc of the Galactic centre. The Galactic tidal forces acting on planets scale inversely with radius in the Galaxy and so are much greater in the inner Galaxy than in the Solar neighbourhood. Within a range of 3.5 to 10 kpc, the vertical tide from the Galactic disc is predominant. Interior to 3.5 kpc, the effects of the Galactic bulge cannot be neglected and the in-plane tidal components are as important as the vertical ones. Here, we quantify the orbital changes induced by these tides. We find that the greatest perturbations occur when the planetary orbit is severely misaligned to the parent star's orbit. When both planes are perpendicular, the eccentricity of the planet is driven to unity, although the semimajor axis is secularly unaffected. When both planes are coincident, the effect from Galactic tides is minimized, but remains non-zero. In these cases, we provide estimates for the survival times, as well as the minimum baseline eccentricity variation for all Milky Way exoplanets as a function of Galactic parameters. Inclinations similar to the Solar System's (about 60 degrees) can easily cause eccentric Neptunes (at about 30 AU) around host stars deep within the Galactic bulge (within 50 pc) to experience eccentricity variations of several tenths, and cause the exoplanets with the widest-known separations (at about 1000 AU) to experience similar variations in the Galactic disc. These variations occur on timescales of a few Gyr, a fraction of a typical main sequence lifetime.

preprint2012arXiv

Planetary Orbital Equations in Externally-Perturbed Systems: Position and Velocity-Dependent Forces

The increasing number and variety of extrasolar planets illustrates the importance of characterizing planetary perturbations. Planetary orbits are typically described by physically intuitive orbital elements. Here, we explicitly express the equations of motion of the unaveraged perturbed two-body problem in terms of planetary orbital elements by using a generalized form of Gauss' equations. We consider a varied set of position and velocity-dependent perturbations, and also derive relevant specific cases of the equations: when they are averaged over fast variables (the "adiabatic" approximation), and in the prograde and retrograde planar cases. In each instance, we delineate the properties of the equations. As brief demonstrations of potential applications, we consider the effect of Galactic tides. We measure the effect on the widest-known exoplanet orbit, Sedna-like objects, and distant scattered disk objects, particularly with regard to where the adiabatic approximation breaks down. The Mathematica code which can help derive the equations is freely available upon request.

preprint2012arXiv

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

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

preprint2011arXiv

Accurate Stellar Kinematics at Faint Magnitudes: application to the Bootes~I dwarf spheroidal galaxy

We develop, implement and characterise an enhanced data reduction approach which delivers precise, accurate, radial velocities from moderate resolution spectroscopy with the fibre-fed VLT/FLAMES+GIRAFFE facility. This facility, with appropriate care, delivers radial velocities adequate to resolve the intrinsic velocity dispersions of the very faint dSph dwarf galaxies. Importantly, repeated measurements let us reliably calibrate our individual velocity errors ($0.2 \leq δ_V\leq 5$ km s$^{-1}$) and directly detect stars with variable radial velocities. We show, by application to the Bootes-1 dwarf spheroidal, that the intrinsic velocity dispersion of this system is significantly below 6.5\,km/s reported by previous studies. Our data favor a two-population model of Bootes-1, consisting of a majority `cold' stellar component, with velocity dispersion $2.4^{+0.9}_{-0.5}$\,km/s, and a minority `hot' stellar component, with velocity dispersion $\sim 9$\,km/s, although we can not completely rule out a single component distribution with velocity dispersion $4.6^{0.8}_{-0.6}$\,km/s. We speculate this complex velocity distribution actually reflects the distribution of velocity anisotropy in Bootes-1, which is a measure of its formation processes.

preprint2010arXiv

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

We present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used to reconstruct the unlensed source profile and the lens galaxy mass profiles. We introduce a prior distribution over multi-wavelength sources that is realistic as a representation of our knowledge about the surface brightness profiles of galaxies and groups of galaxies. To carry out the inference computationally, we use Diffusive Nested Sampling, an efficient variant of Nested Sampling that uses Markov Chain Monte Carlo (MCMC) to sample the complex posterior distributions and compute the normalising constant. We demonstrate the efficacy of this approach with the reconstruction of the group-group gravitational lens system CSWA 31, finding the source to be composed of five merging spiral galaxies magnified by a factor of 13.

preprint2010arXiv

No evidence for internal rotation in the remnant core of the Sagittarius dwarf

We have conducted a spectroscopic survey of the inner regions of the Sagittarius (Sgr) dwarf galaxy using the AAOmega spectrograph on the Anglo-Australian Telescope. We determine radial velocities for over 1800 Sgr star members in 6 fields that cover an area 18.84 deg^2, with a typical accuracy of ~2 km/s. Motivated by recent numerical models of the Sgr tidal stream that predict a substantial amount of rotation in the dwarf remnant core, we compare the kinematic data against N-body models that simulate the stream progenitor as (i) a pressure-supported, mass-follows-light system, and (ii) a late-type, rotating disc galaxy embedded in an extended dark matter halo. We find that the models with little, or no intrinsic rotation clearly yield a better match to the mean line-of-sight velocity in all surveyed fields, but fail to reproduce the shape of the line-of-sight velocity distribution. This result rules out models wherein the prominent bifurcation observed in the leading tail of the Sgr stream was caused by a transfer from intrinsic angular momentum from the progenitor satellite into the tidal stream. It also implies that the trajectory of the young tidal tails has not been affected by internal rotation in the progenitor system. Our finding indicates that new, more elaborate dynamical models, in which the dark and luminous components are treated independently, are necessary for simultaneously reproducing both the internal kinematics of the Sgr dwarf and the available data for the associated tidal stream.

preprint2010arXiv

X-shooter Observations of the Gravitational Lens System CASSOWARY 5

We confirm an eighth gravitational lens system in the CASSOWARY catalogue. Exploratory observations with the X-shooter spectrograph on the VLT show the system CSWA5 to consist of at least three images of a blue star-forming galaxy at z = 1.0686, lensed by an apparent foreground group of red galaxies one of which is at z = 0.3877. The lensed galaxy exhibits a rich spectrum with broad interstellar absorption lines and a wealth of nebular emission lines. Preliminary analysis of these features shows the galaxy to be young, with an age of 25-50 Myr. With a star-formation rate of approximately 20 solar masses/yr, the galaxy has already assembled a stellar mass of 3 x 10^9 solar masses and reached half-solar metallicity. Its blue spectral energy distribution and Balmer line ratios suggest negligible internal dust extinction. A more in-depth analysis of the properties of this system is currently hampered by the lack of a viable lensing model. However, it is already clear that CSWA5 shares many of its physical characteristics with the general population of UV-selected galaxies at redshifts z = 1-3, motivating further study of both the source and the foreground mass concentration responsible for the gravitational lensing.

preprint2009arXiv

A theorem on central velocity dispersions

It is shown that, if the tracer population is supported by a spherical dark halo with a core or a cusp diverging more slowly than that of a singular isothermal sphere, the logarithmic cusp slope 'g' of the tracers must be given exactly by g=2b where b is their velocity anisotropy parameter at the center unless the same tracers are dynamically cold at the center. If the halo cusp diverges faster than that of the singular isothermal sphere, the velocity dispersion of the tracers must diverge at the center too. In particular, if the logarithmic halo cusp slope is larger than two, the diverging velocity dispersion also traces the behavior of the potential. The implication of our theorem on projected quantities is also discussed. We argue that our theorem should be understood as a warning against interpreting results based on simplifying assumptions such as isotropy and spherical symmetry.

preprint2009arXiv

A Universal Mass Profile for Dwarf Spheroidal Galaxies?

(abridged) We apply the Jeans equation to estimate masses for eight of the brightest dSph galaxies. For Fornax we obtain a model-independent constraint on the maximum-circular velocity, Vmax=18_{-3}^{+5} km/s. Although we obtain only lower-limits of Vmax > 10 km/s for the remaining dSphs, we find that in all cases the enclosed mass at the projected half-light radius is well constrained and robust to a wide range of halo models and velocity anisotropies. We derive a simple analytic formula that estimates M(rhalf) accurately with respect to results from the full Jeans analysis. Applying this formula to the entire population of Local Group dSphs, we demonstrate a correlation such that M(rhalf)\propto rhalf^{1.4\pm 0.4}, or in terms of the mean density interior to the half-light radius, <rho> \propto rhalf^{-1.6\pm 0.4}. This relation is driven by the fact that the dSph data exhibit a correlation between global velocity dispersion and half-light radius. We argue that tidal forces are unlikely to have introduced this relation, but tides may have increased the scatter and/or altered the slope. While the data are well described by mass profiles ranging over a factor of < 2 in normalization (Vmax ~ 10-20 km/s), we consider the hypothesis that all dSphs are embedded within a "universal" dark matter halo. We show that in addition to the power law M\propto r^{1.4}, viable candidates include a cuspy "NFW" halo with Vmax ~ 15 km/s and scale radius r_0 ~ 800 pc, as well as a cored halo with Vmax ~ 13 km/s and r_0 ~ 150 pc. Finally, assuming that their measured velocity dispersions accurately reflect their masses, the smallest dSphs now allow us to resolve dSph densities at radii as small as a few tens of pc. At these small scales we find mean densities as large as <ρ> ~ 5 Msun/pc^3 (200 GeV/cm^3).

preprint2009arXiv

CASSOWARY 20: a Wide Separation Einstein Cross Identified with the X-shooter Spectrograph

We have used spectra obtained with X-shooter, the triple arm optical-infrared spectrograph recently commissioned on the Very Large Telescope of the European Southern Observatory, to confirm the gravitational lens nature of the CASSOWARY candidate CSWA 20. This system consists of a luminous red galaxy at redshift z(abs) = 0.741, with a very high velocity dispersion sigma(lens) = 500 km/s, which lenses a blue star-forming galaxy at z(em) = 1.433 into four images with mean separation of about 6 arcseconds. The source shares many of its properties with those of UV-selected galaxies z = 2-3: it is forming stars at a rate of 25 solar masses per year, has a metallicity of about 1/4 solar, and shows nebular emission from two components separated by 0.4 arcseconds in the image plane, possibly indicating a merger. It appears that foreground interstellar material within the galaxy has been evacuated from the sight-line along which we observe the starburst, giving an unextinguished view of its stars and H II regions. CSWA 20, with its massive lensing galaxy producing a high magnification of an intrinsically luminous background galaxy, is a promising target for future studies at a variety of wavelengths.

preprint2009arXiv

The abundance spread in the Bootes I dwarf spheroidal galaxy

We present medium-resolution spectra of 16 radial velocity red-giant members of the low-luminosity Bootes I dwarf spheroidal (dSph) galaxy, that have sufficient S/N for abundance determination, based on the strength of the Ca II K line. Assuming [Ca/Fe] ~ +0.3, the abundance range in the sample is Delta [Fe/H] ~ 1.7 dex, with one star having [Fe/H] = -3.4. The dispersion is sigma([Fe/H]) = 0.45 +/- 0.08 -- similar to those of the Galaxy's more luminous dSph systems and Omega Centauri. This suggests that the large mass (greater than approximately 10 million solar masses) normally assumed to foster self-enrichment and the production of chemical abundance spreads was provided by the non-baryonic material in Bootes I.