Source author record

Douglas Boubert

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

10 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'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.

preprint2016arXiv

A Dipole on the Sky: Predictions for Hypervelocity Stars from the Large Magellanic Cloud

We predict the distribution of hypervelocity stars (HVSs) ejected from the Large Magellanic Cloud (LMC), under the assumption that the dwarf galaxy hosts a central massive black hole (MBH). For the majority of stars ejected from the LMC the orbital velocity of the LMC has contributed a significant fraction of their galactic rest frame velocity, leading to a dipole density distribution on the sky. We quantify the dipole using spherical harmonic analysis and contrast with the monopole expected for HVSs ejected from the Galactic Center. There is a tendril in the density distribution that leads the LMC which is coincident with the well-known and unexplained clustering of HVSs in the constellations of Leo and Sextans. Our model is falsifiable, since it predicts that Gaia will reveal a large density of HVSs in the southern hemisphere.

preprint2016arXiv

Photochemical-dynamical models of externally FUV irradiated protoplanetary discs

There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their "critical radius". This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.