Source author record

Vasily Belokurov

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.HE astro-ph.IM gr-qc hep-ph physics.data-an

Catalog footprint

What is connected

62works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2026arXiv

JADES: the chemical enrichment pattern of distant galaxies -- $α$ enhancement, silicon depletion, and iron enhancement

We present gas-phase abundances of carbon (C), $α$-elements (O, Ne, Si, and Ar) and iron (Fe) obtained from stacked spectra of high-$z$ star-forming galaxies with the deep Near Infrared Spectrograph medium-resolution data from the James Webb Space Telescope Advanced Deep Extragalactic Survey. Our 564 sources at $z=4$--7 have a median stellar mass of $\log(M_{*}/M_{\odot})=8.46$ and a median star-formation rate of $\log(\mathrm{SFR}/M_{\odot}\,\mathrm{yr^{-1}})=0.30$, placing them close to the star-formation main sequence. We find that the stacked spectrum of all our 564 sources has relatively low [C/O]$=-0.70$, moderate [Ne/O]$=-0.09$, and low [Ar/O]$=-0.28$ values at a low gas-phase metallicity of $12+\log(\mathrm{O/H})=7.71$ ($Z\sim 0.1~Z_\odot$), suggesting dominant yields of core-collapse supernovae evolved from massive stars. The detection of a weak SiIII] emission line in our stacked spectrum provides a silicon-to-oxygen abundance ratio of [Si/O]$=-0.63$, which is lower than that of stars in the Milky Way disc and lower than expected by chemical evolution models, suggesting silicon depletion onto dust grains. Likewise, this Si/O value is lower than that we newly derive for two individual $z>6$ galaxies (GN-z11 and RXCJ2248) with negligible dust attenuation. By performing spectral stacking in bins of $M_{*}$, SFR, specific SFR (sSFR), and ultra-violet (UV) continuum slope $β_{\mathrm{UV}}$, we identify [FeIII] line detections in the high-sSFR bin and the blue-$β_{\mathrm{UV}}$ bin, both of which exhibit supersolar Fe/O ratios, while their C/O, Ar/O, and Si/O ratios are comparable to those of the all-sources stack. Our findings support a chemically young gas composition with rapid dust depletion in the general population of high-$z$ star-forming galaxies, while raising the possibility of anomalous, selective Fe/O enhancement at the very early epoch of star formation.

preprint2026arXiv

Multiple Populations and a CH Star Found in the 300S Globular Cluster Stellar Stream

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.

preprint2022arXiv

A gap in the double white dwarf separation distribution caused by the common-envelope evolution: astrometric evidence from Gaia

The trajectory of the center of light of an unresolved binary is different from that of its center of mass. Binary-induced stellar centroid wobbling can therefore be detected as an excess in the goodness-of-fit of the single-star astrometric model. We use reduced $χ^2$ of the astrometric fit in the {\it Gaia} Early Data Release 3 to detect the likely unresolved double white dwarfs (DWDs). Using parallax-based distances we convert the excess of reduced $χ^2$ into the amplitude of the centroid wobble $δa$, which is proportional to the binary separation $a$. The measured $δa$ distribution drops towards larger wobble amplitudes and shows a break around $δa \approx 0.2$ where it steepens. The integral of the distribution yields DWD fraction of $6.5 \pm 3.7$ per cent in the range $0.01 < a (\text{au}) < 2$. Using synthetic models of the Galactic DWDs we demonstrate that the break in the $δa$ distribution corresponds to one side of a deep gap in the DWD separation distribution at around $a\approx 1$ au. Model DWDs with separations less than several au shrink dramatically due to (al least one) common envelope phase, reshaping the original separation distribution, clearing a gap and creating a pile-up of systems with $a\approx 0.01$ au and $δa < 0.01$. Our models reproduce the overall shape of the observed $δa$ distribution and its normalisation, however the predicted drop in the numbers of DWDs beyond the break is steeper than in the data.

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

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

From dawn till disk: Milky Way's turbulent youth revealed by the APOGEE+Gaia data

We use accurate estimates of aluminium abundance provided as part of the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity $-1.5\lesssim {\rm [Fe/H]}\lesssim 0.5$ born in-situ in the Milky Way proper. We show that the low-metallicity ([Fe/H]$\lesssim -1.3$) in-situ component that we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in a number of element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with increasing metallicity between [Fe/H]$=-1.3$ and $-0.9$, the transition that we call the spin-up. The observed and theoretically expected age-metallicity correlations imply that this increase reflects a rapid formation of the Milky Way disk over $\approx 1-2$ Gyrs. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disk during later epochs corresponding to [Fe/H]$>-0.9$. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the Milky Way reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disk evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disk formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disk.

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

Merger-induced galaxy transformations in the ARTEMIS simulations

Using the ARTEMIS set of 45 high-resolution cosmological simulations, we investigate a range of merger-induced dynamical transformations of Milky Way-like galaxies. We first identify populations of accreted stars on highly radial orbits, similar to the 'Gaia Sausage' in the Milky Way. We show that $\approx1/3$ of the ARTEMIS galaxies contain a similar feature, and confirm that they usually comprise stellar debris from the most massive accreted satellite. Selecting 15 galaxies with discs at the present-day, we study their changes around the times of the GS-like mergers. Dark matter haloes of many of these exhibit global changes in shape and orientation, with almost half becoming significantly more spherical when the mergers occur. Focusing on the galaxies themselves, we find that 4/15 have stellar discs which experience large changes in the orientation of their angular momentum (AM) axes, at rates of up to $\sim60$ degrees Gyr$^{-1}$. By calculating the orbital angular momentum axes of the satellites before they are accreted, we show that there is a tendency for the disc's AM to become more aligned with this axis after the merger. We also investigate the origin of in situ retrograde stars, analogous to the 'Splash' in the Milky Way. Tracing them back to earlier snapshots, we demonstrate that they were often disrupted onto their extreme orbits by multiple early mergers. We also find that the total mass of these stars outside the central regions positively correlates with the total accreted stellar mass.

preprint2022arXiv

Quenching of satellite galaxies of Milky Way analogues: reconciling theory and observations

The vast majority of low-mass satellite galaxies around the Milky Way and M31 appear virtually devoid of cool gas and show no signs of recent or ongoing star formation. Cosmological simulations demonstrate that such quenching is expected and is due to the harsh environmental conditions that satellites face when joining the Local Group (LG). However, recent observations of Milky Way analogues in the SAGA survey present a very different picture, showing the majority of observed satellites to be actively forming stars, calling into question the realism of current simulations and the typicality of the LG. Here we use the ARTEMIS suite of high-resolution cosmological hydrodynamical simulations to carry out a careful comparison with observations of dwarf satellites in the LG, SAGA, and the Local Volume (LV) survey. We show that differences between SAGA and the LG and LV surveys, as well as between SAGA and the ARTEMIS simulations, can be strongly reduced by considering differences in the host mass distributions and (more importantly) observational selection effects, specifically that low-mass satellites which have only recently been accreted are more likely to be star-forming, have a higher optical surface brightness, and are therefore more likely to be included in the SAGA survey. This picture is confirmed using data from the deeper LV survey, which shows pronounced quenching at low masses, in accordance with the predictions of LCDM-based simulations.

preprint2022arXiv

Radialization of satellite orbits in galaxy mergers

We consider the orbital evolution of satellites in galaxy mergers, focusing on the evolution of eccentricity. Using a large suite of N-body simulations, we study the phenomenon of satellite orbital radialization -- a profound increase in the eccentricity of its orbit as it decays under dynamical friction. While radialization is detected in a variety of different setups, it is most efficient in the cases of high satellite mass, not very steep host density profiles, and high initial eccentricity. To understand the origin of this phenomenon, we run additional simulations with various physical factors selectively turned off: satellite mass loss, reflex motion and distortion of the host, etc. We find that all these factors are important for radialization, since it does not occur for point-mass satellites or when the host potential is replaced with an unperturbed initial profile. The analysis of forces and torques acting on both galaxies confirms the major role of self-gravity of both host and satellite in the reduction of orbital angular momentum. The classical Chandrasekhar dynamical friction formula, which accounts only for the forces between the host and the satellite, but not for internal distortions of both galaxies, does not match the evolution of eccentricity observed in N-body simulations.

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.

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

Tango for three: Sagittarius, LMC, and the Milky Way

We assemble a catalogue of candidate Sagittarius stream members with 5d and 6d phase-space information, using astrometric data from Gaia DR2, distances estimated from RR Lyrae stars, and line-of-sight velocities from various spectroscopic surveys. We find a clear misalignment between the stream track and the direction of the reflex-corrected proper motions in the leading arm of the stream, which we interpret as a signature of a time-dependent perturbation of the gravitational potential. A likely cause of this perturbation is the recent passage of the most massive Milky Way satellite - the Large Magellanic Cloud (LMC). We develop novel methods for simulating the Sagittarius stream in the presence of the LMC, using specially tailored N-body simulations and a flexible parametrization of the Milky Way halo density profile. We find that while models without the LMC can fit most stream features rather well, they fail to reproduce the misalignment and overestimate the distance to the leading arm apocentre. On the other hand, models with an LMC mass in the range (1.3+-0.3)x10^11 Msun rectify these deficiencies. We demonstrate that the stream can not be modelled adequately in a static Milky Way. Instead, our Galaxy is required to lurch toward the massive in-falling Cloud, giving the Sgr stream its peculiar shape and kinematics. By exploring the parameter space of Milky Way potentials, we determine the enclosed mass within 100 kpc to be (5.6+-0.4)x10^11 Msun, and the virial mass to be (9.0+-1.3)x10^11 Msun, and find tentative evidence for a radially-varying shape and orientation of the Galactic halo.

preprint2021arXiv

The mass of the Milky Way out to 100 kpc using halo stars

We use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority (~98%) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar halos from cosmological simulations tend to underestimate the true halo mass. After applying our method to the Milky Way data we find a mass within 100 kpc of M(< 100 kpc) = 6.07 +/- 0.29 (stat.) +/- 1.21 (sys.) x 10^11 M_Sun. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5% systematic bias, and a 20% systematic uncertainty based on our tests with cosmological simulations. Assuming the mass-concentration relation for Navarro-Frenk-White haloes, our mass estimate favours a total (pre-LMC infall) Milky Way mass of M_200c = 1.01 +/- 0.24 x 10^12 M_Sun, or (post-LMC infall) mass of M_200c = 1.16 +/- 0.24 x 10^12 M_Sun when a 1.5 x 10^11 M_Sun mass of a rigid LMC is included.

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.

preprint2021arXiv

Variable star classification across the Galactic bulge and disc with the VISTA Variables in the Vía Láctea survey

We present VIVACE, the VIrac VAriable Classification Ensemble, a catalogue of variable stars extracted from an automated classification pipeline for the Vista Variables in the Vía Láctea (VVV) infrared survey of the Galactic bar/bulge and southern disc. Our procedure utilises a two-stage hierarchical classifier to first isolate likely variable sources using simple variability summary statistics and training sets of non-variable sources from the Gaia early third data release, and then classify candidate variables using more detailed light curve statistics and training labels primarily from OGLE and VSX. The methodology is applied to point-spread-function photometry for $\sim490$ million light curves from the VIRAC v2 astrometric and photometric catalogue resulting in a catalogue of $\sim1.4$ million likely variable stars, of which $\sim39,000$ are high-confidence (classification probability $>0.9$) RR Lyrae ab stars, $\sim8000$ RR Lyrae c/d stars, $\sim187,000$ detached/semi-detached eclipsing binaries, $\sim18,000$ contact eclipsing binaries, $\sim1400$ classical Cepheid variables and $\sim2200$ Type II Cepheid variables. Comparison with OGLE-4 suggests a completeness of around $90\,\%$ for RRab and $\lesssim60\%$ for RRc/d, and a misclassification rate for known RR Lyrae stars of around $1\%$ for the high confidence sample. We close with two science demonstrations of our new VIVACE catalogue: first, a brief investigation of the spatial and kinematic properties of the RR Lyrae stars within the disc/bulge, demonstrating the spatial elongation of bar-bulge RR Lyrae stars is in the same sense as the more metal-rich red giant population whilst having a slower rotation rate of $\sim40\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}$; and secondly, an investigation of the Gaia EDR3 parallax zeropoint using contact eclipsing binaries across the Galactic disc plane and bulge.

preprint2021arXiv

Weighing Milky Way Satellites with LISA

White dwarf stars are a well-established tool for studying Galactic stellar populations. Two white dwarfs in a tight binary system offer us an additional messenger - gravitational waves - for exploring the Milky Way and its immediate surroundings. Gravitational waves produced by double white dwarf (DWD) binaries can be detected by the future Laser Interferometer Space Antenna (LISA). Numerous and widespread DWDs have the potential to probe shapes, masses and formation histories of the stellar populations in the Galactic neighbourhood. In this work we outline a method for estimating the total stellar mass of Milky Way satellite galaxies based on the number of DWDs detected by LISA. To constrain the mass we perform a Bayesian inference using binary population synthesis models and considering the number of detected DWDs associated with the satellite and the measured distance to the satellite as the only inputs. Using a fiducial binary population synthesis model we find that for large satellites the stellar masses can be recovered to within 1) a factor two if the star formation history is known and 2) an order of magnitude when marginalising over different star formation history models. For smaller satellites we can place upper limits on their stellar mass. Gravitational wave observations can provide mass measurements for large satellites that are comparable, and in some cases more precise, than standard electromagnetic observations.

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

Discovery of a disrupting open cluster far into the Milky Way halo: a recent star formation event in the leading arm of the Magellanic stream?

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties.

preprint2020arXiv

Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$σ$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

preprint2020arXiv

Kinematics of the Palomar 5 stellar stream from RR Lyrae stars

Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work.

preprint2020arXiv

Sausage & Mash: The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus-Sausage merger

We analyse a set of cosmological magneto-hydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called "Gaia Sausage" found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus-Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50$\%$ of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars onto less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits.

preprint2020arXiv

SMASHing the low surface brightness SMC

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

preprint2020arXiv

The 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

Chemo-dynamical properties of the Anticenter Stream: a surviving disc fossil from a past satellite interaction

Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc.

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.

preprint2019arXiv

The little things matter: relating the abundance of ultrafaint satellites to the hosts' assembly history

Ultrafaint dwarf galaxies ($M_\star\lesssim10^{5-6}\,{\rm M}_\odot$) are relics of an early phase of galaxy formation. They contain some of the oldest and most metal-poor stars in the Universe which likely formed before the epoch of hydrogen reionisation. These galaxies are so faint that they can only be detected as satellites of the Milky Way. They are so small that they are just barely resolved in current cosmological hydrodynamics simulations. Here we combine very high resolution cosmological $N$-body simulations with a semi-analytic model of galaxy formation to study the demographics and spatial distribution of ultrafaint satellites in Milky Way-mass haloes. We show that the abundance of these galaxies is correlated with the assembly history of the host halo: at fixed mass, haloes assembled earlier contain, on average, more ultrafaint satellites today than haloes assembled later. We identify simulated galactic haloes that experience an ancient Gaia-Enceladus-Sausage-like and a recent LMC-like accretion event and find that the former occurs in 33% of the sample and the latter in 9%. Only 3% experience both events and these are especially rich in ultrafaint satellites, most acquired during the ancient accretion event. Our models predict that the radial distribution of satellites is more centrally concentrated in early-forming haloes. Accounting for the depletion of satellites by tidal interactions with the central disc, we find a very good match to the observed radial distribution of satellites in the Milky Way over the entire radial range. This agreement is mainly due to the ability of our model to track 'orphan' galaxies after their subhaloes fall below the resolution limit of the simulation.

preprint2016arXiv

A 10 kpc stellar substructure at the edge of the Large Magellanic Cloud: perturbed outer disk or evidence for tidal stripping?

We report the discovery of a substantial stellar overdensity in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The structure appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the feature are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like structure between $\sim9-12$ degrees from the LMC centre. Based on all these observations, we conclude that our remote, stream-like feature is likely comprised of material that has been stripped from the outskirts of the LMC disk, although we cannot rule out that it represents a transient overdensity in the disk itself. We conduct a simple $N$-body simulation to show that either type of structure could plausibly arise due to the tidal force of the Milky Way; however we also recognize that a recent close interaction between the LMC and the SMC may be the source of the stripping or perturbation. Finally, we observe evidence for extremely diffuse LMC populations extending to radii of $\sim 18.5$ kpc in the disk plane ($\approx 20$ degrees on the sky), corroborating previous spectroscopic detections at comparable distances.

preprint2016arXiv

A Magellanic Origin of the DES Dwarfs

We establish the connection between the Magellanic Clouds (MCs) and the dwarf galaxy candidates discovered in the Dark Energy Survey (DES) by building a dynamical model of the MC satellite populations, based on an extensive suite of tailor-made numerical simulations. Our model takes into account the response of the Galaxy to the MCs infall, the dynamical friction experienced by the MCs and the disruption of the MC satellites by their hosts. The simulation suite samples over the uncertainties in the MC's proper motions, the masses of the MW and the Clouds themselves and allows for flexibility in the intrinsic volume density distribution of the MC satellites. As a result, we can accurately reproduce the DES satellites' observed positions and kinematics. Assuming that Milky Way (MW) dwarfs follow the distribution of subhaloes in $Λ$CDM, we further demonstrate that, of 14 observed satellites, the MW halo contributes fewer than 4 (8) of these with 68% (95%) confidence and that 7 (12) DES dwarfs have probabilities greater than 0.7 (0.5) of belonging to the LMC. Marginalising over the entire suite, we constrain the total number of the Magellanic satellites at ~70, the mass of the LMC around $10^{11}M_\odot$ and show that the Clouds have likely endured only one Galactic pericentric passage so far. Finally, we give predictions for the line-of-sight velocities and the proper motions of the satellites discovered in the vicinity of the LMC.

preprint2016arXiv

Stray, swing and scatter: angular momentum evolution of orbits and streams in aspherical potentials

In aspherical potentials orbital planes continuously evolve. The gravitational torques impel the angular momentum vector to precess, that is to slowly stray around the symmetry axis, and nutate, i.e. swing up and down periodically in the perpendicular direction. This familiar orbital pole motion - if detected and measured - can reveal the shape of the underlying gravitational potential, the quantity only crudely gauged in the Galaxy so far. Here we demonstrate that the debris poles of stellar tidal streams show a very similar straying and swinging behavior, and give analytic expressions to link the amplitude and the frequency of the pole evolution to the flattening of the dark matter distribution. While these results are derived for near-circular orbits, we show they are also valid for eccentric orbits. Most importantly, we explain how the differential orbital plane precession leads to the broadening of the stream and show that streams on polar orbits ought to scatter faster. We provide expressions for the stream width evolution as a function of the axisymmetric potential flattening and the angle from the symmetry plane and prove that our models are in good agreement with streams produced in N-body simulations. Interestingly, the same intuition applies to streams whose progenitors are on short or long-axis loops in a triaxial potential. Finally, we present a compilation of the Galactic cold stream data, and discuss how the simple picture developed here, along with stream modelling, can be used to constrain the symmetry axes and flattening of the Milky Way.

preprint2016arXiv

The number and size of subhalo-induced gaps in stellar streams

Ample observational capabilities exist today to detect the small density perturbations that low-mass dark matter subhaloes impart on stellar streams from disrupting Galactic satellites. In anticipation of these observations, we investigate the expected number and size of gaps by combining an analytic prescription for gap evolution on circular orbits with the flux of subhaloes near the stream. We explore the distribution of gap sizes and depths for a typical cold stream around the Milky Way and find that for a given stream age and gap depth, each subhalo mass produces a characteristic gap size. For a stream with an age of a few Gyr, orbiting at a distance of 10-20 kpc from the Galactic center, even modest subhaloes with a mass of $10^6-10^7 M_\odot$ produce gaps with sizes that are on the order of several degrees. We consider the number and distribution of gap sizes created by subhaloes with masses $10^5-10^9 M_\odot$, accounting for the expected depletion of subhaloes by the Milky Way disk, and present predictions for six cold streams around the Milky Way. For Pal 5, we forecast 0.7 gaps with a density depletion of at least 25\% and a typical gap size of $8^\circ$. Thus, there appears to be {\it no tension} between the recent non-detection of density depletions in the Pal 5 tidal tails and $Λ$CDM expectations. These predictions can be used to guide the scale of future gap searches.

preprint2015arXiv

A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

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

Forensics of Subhalo-Stream Encounters: The Three Phases of Gap Growth

There is hope to discover dark matter subhalos free of stars (predicted by the current theory of structure formation) by observing gaps they produce in tidal streams. In fact, this is the most promising technique for dark substructure detection and characterization as such gaps grow with time, magnifying small perturbations into clear signatures observable by ongoing and planned Galaxy surveys. To facilitate such future inference, we develop a comprehensive framework for studies of the growth of the stream density perturbations. Starting with simple assumptions and restricting to streams on circular orbits, we derive analytic formulae that describe the evolution of all gap properties (size, density contrast etc) at all times. We uncover complex, previously unnoticed behavior, with the stream initially forming a density enhancement near the subhalo impact point. Shortly after, a gap forms due to the relative change in period induced by the subhalo's passage. There is an intermediate regime where the gap grows linearly in time. At late times, the particles in the stream overtake each other, forming caustics, and the gap grows like $\sqrt{t}$. In addition to the secular growth, we find that the gap oscillates as it grows due to epicyclic motion. We compare this analytic model to N-body simulations and find an impressive level of agreement. Importantly, when analyzing the observation of a single gap we find a large degeneracy between the subhalo mass, the impact geometry and kinematics, the host potential and the time since flyby.

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

Kinematics and chemistry of recently discovered Reticulum 2 and Horologium 1 dwarf galaxies

Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters and dwarf galaxies because of their overlap in morphological properties. Here we report on VLT/GIRAFFE spectra of candidate member stars in two recently discovered ultra-faint satellites Reticulum 2 and Horologium 1, obtained as part of the ongoing Gaia-ESO Survey. We identify 18 members in Reticulum 2 and 5 in Horologium 1. We find Reticulum 2 to have a velocity dispersion of ~3.22 km/s, implying a M/L ratio of ~ 500. We have inferred stellar parameters for all candidates and we find Reticulum 2 to have a mean metallicity of [Fe/H] = -2.46+/-0.1, with an intrinsic dispersion of ~ 0.29, and is alpha-enhanced to the level of [alpha/Fe]~0.4. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Fe/H] ~ -2 and V_hel ~ 220 km/s, far from the systemic velocity of Reticulum 2. For Horologium 1 we infer a velocity dispersion of 4.9^{+2.8}_{-0.9} km/s and a consequent M/L ratio of ~ 600, leading us to conclude that Horologium 1 is also a dwarf galaxy. Horologium 1 is slightly more metal-poor than Reticulum 2 [Fe/H] = -2.76 +/- 0.1 and is similarly alpha-enhanced: [alpha/Fe] ~ 0.3. Despite a large error-bar, we also measure a significant spread of metallicities of 0.17 dex which strengthen the evidence that Horologium 1 is indeed a dwarf galaxy. The line-of-sight velocity of Reticulum 2 is offset by some 100 km/s from the prediction of the orbital velocity of the LMC, thus making its association with the Cloud uncertain. However, at the location of Horologium 1, both the backward integrated LMC's orbit and the LMC's halo are predicted to have radial velocities similar to that of the dwarf. Therefore, it is likely that Horologium 1 is or once was a member of the Magellanic Family.

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

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

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

preprint2015arXiv

Properties of Dark Subhaloes from Gaps in Tidal Streams

Cold or Warm, the Dark Matter substructure spectrum must extend to objects with masses as low as $10^7 M_\odot$, according to the most recent Lyman-$α$ measurements. Around a Milky Way-like galaxy, more than a thousand of these subhaloes will not be able to form stars but are dense enough to survive even deep down in the potential well of their host. There, within the stellar halo, these dark pellets will bombard tidal streams as they travel around the Galaxy, causing small but recognizable damage to the stream density distribution. The detection and characterization of these stream ruptures will allow us to constrain the details of the subhalo-stream interaction. In this work, for the first time, we will demonstrate how the properties of a subhalo, most importantly its mass and size, can be reliably inferred from the gap it produces in a tidal stream. For a range of realistic observational setups, mimicking e.g. SDSS, DES, Gaia and LSST data, we find that it is possible to measure the {\it complete set} of properties (including the phase-space coordinates during the flyby) of dark perturbers with $M>10^7 M_\odot$, up to a 1d degeneracy between the mass and velocity.

preprint2015arXiv

Satellites of LMC-Mass Dwarfs: Close Friendships Ruined by Milky Way Mass Halos

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhalos, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.

preprint2015arXiv

Uncovering Blue Diffuse Dwarf Galaxies

Extremely metal poor (XMP) galaxies are known to be very rare, despite the large numbers of low-mass galaxies predicted by the local galaxy luminosity function. This paper presents a sub-sample of galaxies that were selected via a morphology-based search on SDSS images with the aim of finding these elusive XMP galaxies. By using the recently discovered extremely metal-poor galaxy, Leo P, as a guide, we obtained a collection of faint, blue systems, each with isolated HII regions embedded in a diffuse continuum, that have remained optically undetected until now. Here we show the first results from optical spectroscopic follow-up observations of 12 of ~100 of these blue, diffuse dwarf (BDD) galaxies yielded by our search algorithm. Oxygen abundances were obtained via the direct method for eight galaxies, and found to be in the range 7.45<12+log(O/H)<8.0, with two galaxies being classified as XMPs. All BDDs were found to currently have a young star-forming population (<10 Myr) and relatively high ionisation parameters of their HII regions. Despite their low luminosities (-11< M_B < -18) and low surface brightnesses (~23-25 mag arcsec^-2), the galaxies were found to be actively star-forming, with current star-formation rates between 0.0003 and 0.078 Msol/yr. From our current subsample, BDD galaxies appear to be a population of non-quiescent dwarf irregular (dIrr) galaxies, or the diffuse counterparts to blue compact galaxies (BCDs) and as such may bridge the gap between these two populations. Our search algorithm demonstrates that morphology-based searches are successful in uncovering more diffuse metal-poor star-forming galaxies, which traditional emission-line based searches overlook.

preprint2014arXiv

Balancing mass and momentum in the Local Group

In the rest frame of the Local Group (LG), the total momentum of the Milky Way (MW) and Andromeda (M31) should balance to zero. We use this fact to constrain new solutions for the solar motion with respect to the LG centre-of-mass, the total mass of the LG, and the individual masses of M31 and the MW. Using the set of remote LG galaxies at $>350$ kpc from the MW and M31, we find that the solar motion has amplitude $V_{\odot}=299\pm 15 {\rm ~km~s^{-1}}$ in a direction pointing toward galactic longitude $l_{\odot}=98.4^{\circ}\pm 3.6^{\circ}$ and galactic latitude $b_{\odot}=-5.9^{\circ}\pm 3.0^{\circ}$. The velocities of M31 and the MW in this rest frame give a direct measurement of their mass ratio, for which we find $\log_{10} (M_{\rm M31}/M_{\rm MW})=0.36 \pm 0.29$. We combine these measurements with the virial theorem to estimate the total mass within the LG as $M_{\rm LG}=(2.5\pm 0.4)\times 10^{12}~{\rm M}_{\odot}$. Our value for $M_{\rm LG}$ is consistent with the sum of literature values for $M_{\rm MW}$ and $M_{\rm M31}$. This suggests that the mass of the LG is almost entirely located within the two largest galaxies rather than being dispersed on larger scales or in a background medium. The outskirts of the LG are seemingly rather empty. Combining our measurement for $M_{\rm LG}$ and the mass ratio, we estimate the individual masses of the MW and M31 to be $M_{\rm MW}=(0.8\pm 0.5)\times 10^{12}~{\rm M}_{\odot}$ and $M_{\rm M31}=(1.7\pm 0.3)\times 10^{12}~{\rm M}_{\odot}$, respectively. Our analysis favours M31 being more massive than the MW by a factor of $\sim$2.3, and the uncertainties allow only a small probability (9.8%) that the MW is more massive. This is consistent with other properties such as the maximum rotational velocities, total stellar content, and numbers of globular clusters and dwarf satellites, which all suggest that $M_{\rm M31}/M_{\rm MW}>1$.

preprint2014arXiv

Strong RR Lyrae excess in the Hercules-Aquila Cloud

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30° < l < 55° and -45° < b < -25°, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.

preprint2014arXiv

The evolution of late-type galaxies from CASSOWARY lensing systems

We explore the properties of lensing galaxies and lensed faint sources at redshifts between 1.5 and 3.0. Our sample consists of 9 carefully selected strongly-lensed galaxies discovered by the CASSOWARY survey in the Sloan Digital Sky Survey (SDSS) data. We find that, despite some limitations of the original SDSS data, the homogeneous sample of lensing systems can provide a useful insight into lens and source properties. We also explore the limitations of using low-resolution data to model and analyse galaxy-galaxy lensing. We derive the relative alignment of mass and light in fitted lens profiles. The range of magnification extends above 5, hence we are able to analyse potentially small and low-mass galaxies at high redshifts. We confirm the likely evolution of the size-luminosity relation for blue star-forming galaxies as a function of redshift

preprint2014arXiv

Touching The Void: A Striking Drop in Stellar Halo Density Beyond 50 kpc

We use A-type stars selected from Sloan Digital Sky Survey data release 9 photometry to measure the outer slope of the Milky Way stellar halo density profile beyond 50 kpc. A likelihood-based analysis is employed that models the ugr photometry distribution of blue horizontal branch (BHB) and blue straggler (BS) stars. In the magnitude range, 18.5 < g < 20.5, these stellar populations span a heliocentric distance range of: 10 kpc < D_BS < 75 kpc, 40 kpc < D_BHB < 100 kpc. Contributions from contaminants, such as QSOs, and the effect of photometric uncertainties, are also included in our modeling procedure. We find evidence for a very steep outer halo profile, with power-law index alpha ~ 6 beyond Galactocentric radii r=50 kpc, and even steeper slopes favored (alpha ~ 6-10) at larger radii. This result holds true when stars belonging to known overdensities, such as the Sagittarius stream, are included or excluded. We show that, by comparison to numerical simulations, stellar halos with shallower slopes at large distances tend to have more recent accretion activity. Thus, it is likely that the Milky Way has undergone a relatively quiet accretion history over the past several Gyr. Our measurement of the outer stellar halo profile may have important implications for dynamical mass models of the Milky Way, where the tracer density profile is strongly degenerate with total mass-estimates.

preprint2013arXiv

A 10,000 star spectroscopic survey of the thick disk-halo interface : Phase-space sub-structure in the thick disk

We analyse a 10,000 star spectroscopic survey, focused on Galactic thick disk stars typically 2-5 kpc from the Sun, carried out using the AAOmega Spectrograph on the AAT. We develop methods for completeness-correction of the survey based on SDSS photometry, and we derive star distances using an improved isochrone-fitting method with accuracies better than 10%. We determine the large-scale kinematic (Vphi=172 km/s and sigma(phi,r,z) = (49,51,40) km/s), and abundance properties of the thick disk, showing these representative values are a fair description within about 3 kpc of the sun, and in the range 1-3 kpc from the Galactic Plane. We identify a substantial overdensity in lines of sight towards the inner Galaxy, with metallicity [Fe/H] -1 dex, and higher line of sight velocities than the thick disk, localised along the direction (l,b)=(48,-26). This overdensity appears to be towards, but closer than, the known Hercules-Aquila halo overdensity and may be related to the Humphreys-Larsen inner galaxy thick disk asymmetry.

preprint2013arXiv

Galactic Archaeology. The dwarfs that survived and perished

From the archaeological point of view, the local dwarf galaxies are unique objects in which the imprint of the conditions that shaped the early structure formation can be studied today at high resolution. Over the last decade, this new window into the high redshift Universe has started to be exploited using deep wide-field imaging, high resolution spectroscopy and cutting edge N-body and hydro-dynamical simulations. We review the recent advances in the observational studies of the Milky Way dwarf galaxies, with the aim to understand the properties of the population as a whole and to assist an objective comparison between the models and the data.

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 CASSOWARY spectroscopy survey: A new sample of gravitationally lensed galaxies in SDSS

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

preprint2012arXiv

A Tale Twice Told: The Luminosity Profiles of the Sagittarius Tails

The Sagittarius dwarf galaxy is the archetype of a tidally disrupting system. Both leading and trailing tails can be observed across at least 180 degrees of the sky and measurements of their luminosity density profiles have recently become available. Using numerical simulations, we explore the factors that control the appearance of such profiles. We use two possible models for the Sgr progenitor. The first is a one-component Plummer model, which may represent either a dark matter free progenitor, or one in which pre-existing dark matter has already been largely stripped. The second is a two-component model in which the stars are represented by a Hernquist sphere embedded in a cosmologically modish Navarro-Frenk-White dark halo. Disruption of the models in the Milky Way galaxy provides us with two tellings of the tale of the formation of the Sgr stream. The initial disintegration of the baryons proceeds more slowly for the two-component models because of the protective cocoon of dark matter. Once this has been stripped, though, matters proceed apace. In both cases, the profiles after ~6 pericentric passages provide good matches to the observational data, but the tails are more extended for the two-component models. The leading and trailing tails are symmetric at apocentre or pericentre. At other orbital phases, asymmetries are present, as tails are compressed as they approach apocentre and stretched out as they approach pericentre. There may exist density enhancements corresponding to such pile-ups which may be observable in current survey data. We re-visit the calculation of Niederste-Ostholt et al. (2010) and slightly revise upwards the luminosity of the Sgr progenitor to 9.9-14.4 x 10^7 solar luminosities based on insights from the simulations.

preprint2012arXiv

Constraining the dark energy equation of state with double source plane strong lenses

We investigate the possibility of constraining the dark energy equation of state by measuring the ratio of Einstein radii in a strong gravitational lens system with two source planes. This quantity is independent of the Hubble parameter and directly measures the growth of angular diameter distances as a function of redshift. We investigate the prospects for a single double source plane system and for a forecast population of systems discovered by re-observing a population of single source lenses already known from a photometrically selected catalogue such as CASSOWARY or from a spectroscopically selected catalogue such as SLACS. We find that constraints comparable to current data-sets (15% uncertainty on the dark equation of state at 68%CL) are possible with a handful of double source plane systems. We also find that the method's degeneracy between Omega_M and w is almost orthogonal to that of CMB and BAO measurements, making this method highly complimentary to current probes.

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

preprint2010arXiv

A strong redshift dependence of the broad absorption line quasar fraction

We describe the application of non-negative matrix factorisation to generate compact reconstructions of quasar spectra from the Sloan Digital Sky Survey (SDSS), with particular reference to broad absorption line quasars (BALQSOs). BAL properties are measured for SiIV lambda1400, CIV lambda1550, AlIII lambda1860 and MgII lambda2800, resulting in a catalogue of 3547 BALQSOs. Two corrections, based on extensive testing of synthetic BALQSO spectra, are applied in order to estimate the intrinsic fraction of CIV BALQSOs. First, the probability of an observed BALQSO spectrum being identified as such by our algorithm is calculated as a function of redshift, signal-to-noise ratio and BAL properties. Second, the different completenesses of the SDSS target selection algorithm for BALQSOs and non-BAL quasars are quantified. Accounting for these selection effects the intrinsic CIV BALQSO fraction is 41+/-5 per cent. Our analysis of the selection effects allows us to measure the dependence of the intrinsic CIV BALQSO fraction on luminosity and redshift. We find a factor of 3.5+/-0.4 decrease in the intrinsic fraction from the highest redshifts, z~4.0, down to z~2.0. The redshift dependence implies that an orientation effect alone is not sufficient to explain the presence of BAL troughs in some but not all quasar spectra. Our results are consistent with the intrinsic BALQSO fraction having no strong luminosity dependence, although with 3-sigma limits on the rate of change of the intrinsic fraction with luminosity of -6.9 and 7.0 per cent dex^-1 we are unable to rule out such a dependence.

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

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.

Vasily Belokurov

What is connected

Connect this record

See the researcher in context

Building this map preview

62 published item(s)

JADES: the chemical enrichment pattern of distant galaxies -- $α$ enhancement, silicon depletion, and iron enhancement

Multiple Populations and a CH Star Found in the 300S Globular Cluster Stellar Stream

A gap in the double white dwarf separation distribution caused by the common-envelope evolution: astrometric evidence from Gaia

An empirical model of the Gaia DR3 selection function

Astrometric Identification of Nearby Binary Stars I: Predicted Astrometric Signals

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

Catalogues of Active Galactic Nuclei From Gaia and unWISE Data

From dawn till disk: Milky Way's turbulent youth revealed by the APOGEE+Gaia data

Information content of BP/RP spectra in Gaia DR3

Merger-induced galaxy transformations in the ARTEMIS simulations

Quenching of satellite galaxies of Milky Way analogues: reconciling theory and observations

Radialization of satellite orbits in galaxy mergers

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

Discovering Strongly-lensed QSOs From Unresolved Light Curves

Elevated r-process enrichment in Gaia Sausage and Sequoia

Tango for three: Sagittarius, LMC, and the Milky Way

The mass of the Milky Way out to 100 kpc using halo stars

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

Variable star classification across the Galactic bulge and disc with the VISTA Variables in the Vía Láctea survey

Weighing Milky Way Satellites with LISA

Binary deviations from single object astrometry

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

Discovery of a disrupting open cluster far into the Milky Way halo: a recent star formation event in the leading arm of the Magellanic stream?

Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

Kinematics of the Palomar 5 stellar stream from RR Lyrae stars

Sausage & Mash: The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus-Sausage merger

SMASHing the low surface brightness SMC

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

Unresolved stellar companions with Gaia DR2 astrometry

Chemo-dynamical properties of the Anticenter Stream: a surviving disc fossil from a past satellite interaction

The Biggest Splash

The little things matter: relating the abundance of ultrafaint satellites to the hosts' assembly history

A 10 kpc stellar substructure at the edge of the Large Magellanic Cloud: perturbed outer disk or evidence for tidal stripping?

A Magellanic Origin of the DES Dwarfs

Stray, swing and scatter: angular momentum evolution of orbits and streams in aspherical potentials

The number and size of subhalo-induced gaps in stellar streams

A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

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

Forensics of Subhalo-Stream Encounters: The Three Phases of Gap Growth

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

Kinematics and chemistry of recently discovered Reticulum 2 and Horologium 1 dwarf galaxies

Magellan/M2FS Spectroscopy of the Reticulum 2 Dwarf Spheroidal Galaxy

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

Properties of Dark Subhaloes from Gaps in Tidal Streams

Satellites of LMC-Mass Dwarfs: Close Friendships Ruined by Milky Way Mass Halos

Uncovering Blue Diffuse Dwarf Galaxies

Balancing mass and momentum in the Local Group

Strong RR Lyrae excess in the Hercules-Aquila Cloud

The evolution of late-type galaxies from CASSOWARY lensing systems

Touching The Void: A Striking Drop in Stellar Halo Density Beyond 50 kpc

A 10,000 star spectroscopic survey of the thick disk-halo interface : Phase-space sub-structure in the thick disk

Galactic Archaeology. The dwarfs that survived and perished

Sagittarius stream 3-d kinematics from SDSS Stripe 82

The CASSOWARY spectroscopy survey: A new sample of gravitationally lensed galaxies in SDSS

A Tale Twice Told: The Luminosity Profiles of the Sagittarius Tails

Constraining the dark energy equation of state with double source plane strong lenses

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

A strong redshift dependence of the broad absorption line quasar fraction

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

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

X-shooter Observations of the Gravitational Lens System CASSOWARY 5

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