Researcher profile

Douglas Boubert

Douglas Boubert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
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

8 published item(s)

preprint2021arXiv

A selection function toolbox for subsets of astronomical catalogues

Large catalogues are ubiquitous throughout astronomy, but most scientific analyses are carried out on smaller samples selected from these catalogues by chosen cuts on catalogued quantities. The selection function of that scientific sample - the probability that a star in the catalogue will satisfy these cuts and so make it into the sample - is thus unique to each scientific analysis. We have created a general framework that can flexibly estimate the selection function of a sample drawn from a catalogue in terms of position, magnitude and colour. Our method is unique in using the binomial likelihood and accounting for correlations in the selection function across position, magnitude and colour using Gaussian processes and spherical harmonics. We have created a new open-source Python package selectionfunctiontoolbox that implements this framework and used it to make three different estimates of the APOGEE DR16 red giant sample selection function, as a subset of 2MASS, with each estimate using an increasing amount of technical knowledge of the APOGEE targeting. In a companion paper we applied our methodology to derive estimates of the astrometric and spectroscopic selection functions of Gaia EDR3. Our framework will make it trivial for astrophysicists to estimate the selection function that they should be using with the custom sample of stars that they have chosen to answer their scientific question.

preprint2021arXiv

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

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

preprint2020arXiv

Completeness of the Gaia-verse I: when and where were Gaia's eyes on the sky during DR2?

The Gaia space mission is crafting revolutionary astrometric, photometric and spectroscopic catalogues that will allow us to map our Galaxy, but only if we know the completeness of this Gaia-verse of catalogues: what stars does it contain and what stars is it missing? We argue that the completeness is driven by Gaia's spinning-and-precessing scanning law and will apply this principle to the Gaia-verse over this series. We take a first step by identifying the periods in time that did not contribute any measurements to Gaia DR2; these gaps create ribbons of incompleteness across the sky that will bias any study that ignores them, although some of these gaps may be filled in future data releases. Our first approach was to use the variable star photometry to identify the 94 gaps longer than 1% of a day. Our second approach was to predict the number of observations of every point on the sky, which in comparison to the reported number of detections revealed additional gaps in the astrometry and spectroscopy. Making these predictions required us to make the most precise, publicly-available determination of the Gaia scanning law. Using this scanning law, we further identified that most stars fainter than $G=22$ in DR2 have spurious magnitudes due to a miscalibration resulting from a thunderstorm over Madrid. Our list of gaps and precision scanning law will allow astronomers to know when Gaia's eye was truly on their binary star, exoplanet or microlensing event during the time period of the second data release.

preprint2020arXiv

Completeness of the Gaia-verse II: what are the odds that a star is missing from Gaia DR2?

The second data release of the Gaia mission contained astrometry and photometry for an incredible 1,692,919,135 sources, but how many sources did Gaia miss and where do they lie on the sky? The answer to this question will be crucial for any astronomer attempting to map the Milky Way with Gaia DR2. We infer the completeness of Gaia DR2 by exploiting the fact that it only contains sources with at least five astrometric detections. The odds that a source achieves those five detections depends on both the number of observations and the probability that an observation of that source results in a detection. We predict the number of times that each source was observed by Gaia and assume that the probability of detection is either a function of magnitude or a distribution as a function of magnitude. We fit both these models to the 1.7 billion stars of Gaia DR2, and thus are able to robustly predict the completeness of Gaia across the sky as a function of magnitude. We extend our selection function to account for crowding in dense regions of the sky, and show that this is vitally important, particularly in the Galactic bulge and the Large and Small Magellanic Clouds. We find that the magnitude limit at which Gaia is still 99% complete varies over the sky from $G=18.9$ to $21.3$. We have created a new Python package selectionfunctions (https://github.com/gaiaverse/selectionfunctions) which provides easy access to our selection functions.

preprint2020arXiv

Deflection of the hypervelocity stars by the dance of the Milky Way and Large Magellanic Cloud

Stars slingshotted by the supermassive black hole at the Galactic centre will escape the Milky Way so quickly that their trajectories will be almost straight lines. Previous works have shown how these `hypervelocity stars' are subsequently deflected by the gravitational field of the Milky Way and the Large Magellanic Cloud (LMC), but have neglected to account for the reflex motion of the Milky Way in response to the fly by of the LMC. A consequence of this motion is that the hypervelocity stars we see on the outskirts of the Milky Way today were ejected from where the Milky Way centre was hundreds of millions of years ago. This change in perspective causes large apparent deflections in the trajectories of the hypervelocity stars, which are of the same order as the deflections caused by the gravitational force of the Milky Way and LMC. We quantify these deflections by simulating the production of hypervelocity stars in an isolated Milky Way (with a spherical or flattened dark matter halo), in a fixed-in-place Milky Way with a passing LMC, and in a Milky Way which responds to the passage of the LMC. The proper motion precision necessary to measure these deflections will be possible with the combination of Gaia with the proposed GaiaNIR successor mission, and these measurements will unlock the hypervelocity stars as probes of the shape of the Milky Way, the mass of the LMC, and of the dance of these two galaxies.

preprint2020arXiv

Kinematic study of the association Cyg OB3 with Gaia DR2

We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is ~22 km/s, dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds <20 km/s with respect to the association for a majority of the sample stars.

preprint2020arXiv

Predictions of Gaia&#39;s prize microlensing events are flawed

Precision astrometry from the second Gaia data release has allowed astronomers to predict 5,787 microlensing events, with 528 of these having maximums within the extended Gaia mission (J2014.5 - J2026.5). Future analysis of the Gaia time-series astrometry of these events will, in some cases, lead to precise gravitational mass measurements of the lens. We find that 61% of events predicted during the extended Gaia mission with sources brighter than G = 18 are likely to be spurious, with the background source in these cases commonly being either a duplicate detection or a binary companion of the lens. We present quality cuts to identify these spurious events and a revised list of microlensing event candidates. Our findings imply that half of the predictable astrometric microlensing events during the Gaia mission have yet to be identified.

preprint2019arXiv

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

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