BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Vasily Belokurov

Vasily Belokurov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.SRastro-ph.COastro-ph.IMastro-ph.HEgr-qchep-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An empirical model of the Gaia DR3 selection function

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues&#39; 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&#39;&#39; exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia&#39;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&#39;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}$

0 citations0 reviewsNo code linkedOpen access
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&#39;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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From dawn till disk: Milky Way&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#39;Gaia Sausage&#39; 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&#39;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 &#39;Splash&#39; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;s stellar halo is imprinted with substructure from several accretion events. Evidence of these events is found in &#34;the Shards&#34;, 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 &#34;Dark Shards&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;Gaia Sausage&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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&#39;&#39; 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Halo&#39;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&#39;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 &#34;Gaia Sausage&#34;. 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 &#34;Sausage&#34; contributes at least 50 % of the Galactic halo. The fraction of &#34;Sausage&#34; stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;thin&#34; and &#34;thick&#34; 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 &#34;thick&#34;) 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&#39;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&#39;s formation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The little things matter: relating the abundance of ultrafaint satellites to the hosts&#39; 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 &#39;orphan&#39; galaxies after their subhaloes fall below the resolution limit of the simulation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access