Gaia Radial Velocity Research Articles

The Galactic Bulge Exploration. III. Calcium Triplet Metallicities for RR Lyrae Stars

RR Lyrae stars (RRLs) are excellent tracers of stellar populations for old, metal-poor components in the the Milky Way and the Local Group. Their luminosities have a metallicity dependence, but determining spectroscopic [Fe/H] metallicities for RRLs, especially at distances outside the solar neighborhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(ii) and Fe(i) abundances, we verify the calibration between the [Fe/H] of RRLs from the calcium triplet. Our calibration is applied to all RRLs with Gaia Radial Velocity Spectrometer (RVS) spectra in Gaia DR3 and to 80 stars in the inner Galaxy from the BRAVA-RR survey. The coadded Gaia RVS RRL spectra provide RRL metallicities with an uncertainty of 0.25 dex, which is a factor of two improvement over the Gaia photometric RRL metallicities. Within our Galactic bulge RRL sample, we find a dominant fraction with low energies without a prominent rotating component. Due to the large fraction of such stars, we interpret these stars as belonging to the in situ metal-poor Galactic bulge component, although we cannot rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles.

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Context. Gravitational waves from black-hole (BH) merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models – and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Aims. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. Methods. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. Results. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 ± 0.82 M⊙ BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. Conclusions. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way.

Adaptation of the phase distance correlation periodogram to account for measurement uncertainties

We present an improvement of the phase distance correlation (PDC) periodogram to account for uncertainties in the time-series data. The PDC periodogram introduced in our previous papers is based on the statistical concept of distance correlation. By viewing each measurement and its accompanying error estimate as a probability distribution, we are able to use the concept of energy distance to design a distance function (metric) between measurement-uncertainty pairs. We used this metric as the basis for the PDC periodogram, instead of the simple absolute difference. We demonstrate the periodogram's performance using both simulated and real-life data. This adaptation makes the PDC periodogram much more useful, demonstrating it can be helpful in the exploration of large time-resolved astronomical databases, ranging from Gaia radial velocity and photometry data releases to those of smaller surveys, such as APOGEE and LAMOST. We have made a public GitHub repository available, with a Python implementation of the new tools available to the community.

High mass function ellipsoidal variables in the Gaia Focused Product Release: searching for black hole candidates in the binary zoo

The recent Gaia Focused Product Release contains radial velocity time-series for more than 9,000 Gaia long-period photometric variables. Here we search for binary systems with large radial velocity amplitudes to identify candidates with massive, unseen companions. Eight targets have binary mass function f(M)>1M⊙, three of which are eclipsing binaries. The remaining five show evidence of ellipsoidal modulations. We fit spectroscopic orbit models to the Gaia radial velocities, and fit the spectral energy distributions of three targets. For the two systems most likely to host dark companions, J0946 and J1640, we use PHOEBE to fit the ASAS-SN light curves and Gaia radial velocities. The derived companion masses are >3M⊙, but the high Galactic dust extinctions towards these objects limit our ability to rule out main sequence companions or subgiants hotter than the photometric primaries. These systems are similar to other stellar-mass black hole impostors, notably the Unicorn (V723 Mon) and the Giraffe (2M04123153+6738486). While it is possible that J1640 and J0946 are similar examples of stripped giant star binaries, high-resolution spectra can be used to determine the nature of their companions.

The Sun Remains Relatively Refractory Depleted: Elemental Abundances for 17,412 Gaia RVS Solar Analogs and 50 Planet Hosts

The element abundances of stars, particularly the refractory elements (e.g., Fe, Si, and Mg), play an important role in connecting stars to their planets. Most Sun-like stars do not have refractory abundance measurements since obtaining a large sample of high-resolution spectra is difficult with oversubscribed observing resources. In this work we infer abundances for C, N, O, Na, Mn, Cr, Si, Fe, Ni, Mg, V, Ca, Ti, Al, and Y for solar analogs with Gaia Radial Velocity Spectrometer (RVS) spectra (R = 11,200) using TheCannon, a data-driven method. We train a linear model on a reference set of 34 stars observed by Gaia RVS with precise abundances measured from previous high-resolution spectroscopic efforts (R > 30,000–110,000). We then apply this model to several thousand Gaia RVS solar analogs. This yields abundances with average upper limit precisions of 0.04–0.1 dex for 17,412 stars, 50 of which are identified planet (candidate) hosts. We subsequently test the relative refractory depletion of these stars with increasing element condensation temperature compared to the Sun. The Sun remains refractory depleted compared to other Sun-like stars regardless of our current knowledge of the planets they host. This is inconsistent with theories of various types of planets locking up or sequestering refractories. Furthermore, we find no significant abundance differences between identified close-in giant planet hosts, giant planet hosts, and terrestrial/small planet hosts with the rest of the sample within our precision limits. This work demonstrates the utility of data-driven learning for future exoplanet composition and demographics studies.

Diffuse interstellar bands in Gaia DR3 RVS spectra

Diffuse interstellar bands (DIBs) are weak and broad interstellar absorption features in astronomical spectra that originate from unknown molecules. To measure DIBs in spectra of late-type stars more accurately and more efficiently, we developed a random forest model to isolate the DIB features from the stellar components. We applied this method to 780 thousand spectra collected by the Gaia Radial Velocity Spectrometer (RVS) that were published in the third data release (DR3). After subtracting the stellar components, we modeled the DIB at 8621 Å (λ8621) with a Gaussian function and the DIB around 8648 Å (λ8648) with a Lorentzian function. After quality control, we selected 7619 reliable measurements for DIB λ8621. The equivalent width (EW) of DIB λ8621 presented a moderate linear correlation with dust reddening, which was consistent with our previous measurements in Gaia DR3 and the newly focused product release. The rest-frame wavelength of DIB λ8621 was updated as λ0 = 8623.141 ± 0.030 Å in vacuum, corresponding to 8620.766 Å in air, which was determined by 77 DIB measurements toward the Galactic anticenter. The mean uncertainty of the fit central wave-length of these 77 measurements is 0.256 Å. With the peak-finding method and a coarse analysis, DIB λ8621 was found to correlate better with the neutral hydrogen than with the molecular hydrogen (represented by 12CO J = (1−0) emission). We also obtained 179 reliable measurements of DIB λ8648 in the RVS spectra of individual stars for the first time, further confirming this very broad DIB feature. Its EW and central wavelength presented a linear relation with those of DIB λ8621. A rough estimation of λ0 for DIB λ8648 was 8646.31 Å in vacuum, corresponding to 8643.93 Å in air, assuming that the carriers of λ8621 and λ8648 are comoving. Finally, we confirmed the impact of stellar residuals on the DIB measurements in Gaia DR3, which led to a distortion of the DIB profile and a shift of the center (≲0.5 Å), but the EW was consistent with our new measurements. With our measurements and analyses, we propose that the approach based on machine learning can be widely applied to measure DIBs in numerous spectra from spectroscopic surveys.

Cepheid Metallicity in the Leavitt Law (C- MetaLL) survey

Context. Classical Cepheids (DCEPs) play a fundamental role in the calibration of the extragalactic distance ladder, which eventually leads to the determination of the Hubble constant (H0) thanks to the period–luminosity (PL) and period–Wesenheit (PW) relations exhibited by these pulsating variables. Therefore, it is of great importance to establish the dependence of PL and PW relations on metallicity. Aims. We aim to quantify the metallicity dependence of the PL and PW relations of the Galactic DCEPs for a variety of photometric bands, ranging from optical to near-infrared. Methods. We gathered a literature sample of 910 DCEPs with available [Fe/H] values from high-resolution spectroscopy or metallicities from the Gaia Radial Velocity Spectrometer. For all these stars, we collected photometry in the GBP, GRP, G, I, V, J, H, and KS bands and astrometry from Gaia Data Release 3 (DR3). We used these data to investigate the metal dependence of both the intercepts and slopes of a variety of PL and PW relations at multiple wavelengths. Results. We find a large negative metallicity effect on the intercept (γ coefficient) of all the PL and PW relations investigated in this work, while present data still do not allow us to draw firm conclusions regarding the metal dependence of the slope (δ coefficient). The typical values of γ are around −0.4 : −0.5 mag dex−1, which is larger than most of the recent determinations present in the literature. We carried out several tests, which confirm the robustness of our results. As in our previous works, we find that the inclusion of a global zero point offset of Gaia parallaxes provides smaller values of γ (in an absolute sense). However, the assumption of the geometric distance of the Large Magellanic Cloud seems to indicate that larger values of γ (in an absolute sense) would be preferred.

Gaia Focused Product Release: Spatial distribution of two diffuse interstellar bands

Diffuse interstellar bands (DIBs) are absorption features seen in optical and infrared spectra of stars and extragalactic objects that are probably caused by large and complex molecules in the galactic interstellar medium (ISM). Here we investigate the Galactic distribution and properties of two DIBs identified in almost six million stellar spectra collected by the Gaia Radial Velocity Spectrometer. These measurements constitute a part of the Gaia Focused Product Release to be made public between the Gaia DR3 and DR4 data releases. In order to isolate the DIB signal from the stellar features in each individual spectrum, we identified a set of 160 000 spectra at high Galactic latitudes (|b| ⩾ 65°) covering a range of stellar parameters which we consider to be the DIB-free reference sample. Matching each target spectrum to its closest reference spectra in stellar parameter space allowed us to remove the stellar spectrum empirically, without reference to stellar models, leaving a set of six million ISM spectra. Using the star’s parallax and sky coordinates, we then allocated each ISM spectrum to a voxel (VOlume piXEL) on a contiguous three-dimensional grid with an angular size of 1.8° (level 5 HEALPix) and 29 unequally sized distance bins. Identifying the two DIBs at 862.1 nm (λ862.1) and 864.8 nm (λ864.8) in the stacked spectra, we modelled their shapes and report the depth, central wavelength, width, and equivalent width (EW) for each, along with confidence bounds on these measurements. We then explored the properties and distributions of these quantities and compared them with similar measurements from other surveys. Our main results are as follows: (1) the strength and spatial distribution of the DIB λ862.1 are very consistent with what was found in Gaia DR3, but for this work we attained a higher signal-to-noise ratio in the stacked spectra to larger distances, which allowed us to trace DIBs in the outer spiral arm and beyond the Scutum–Centaurus spiral arm; (2) we produced an all-sky map below ±65° of Galactic latitude to ~4000 pc of both DIB features and their correlations; (3) we detected the signals of DIB λ862.1 inside the Local Bubble (≲200 pc); and (4) there is a reasonable correlation with the dust reddening found from stellar absorption and EWs of both DIBs with a correlation coefficient of 0.90 for λ862.1 and 0.77 for λ864.8.

Verification of Gaia Data Release 3 Single-lined Spectroscopic Binary Solutions With Three Transiting Low-mass Secondaries

While secondary mass inferences based on single-lined spectroscopic binary (SB1) solutions are subject to degeneracies, this degeneracy can be lifted through the observations of eclipses. We combine the subset of Gaia Data Release 3 SB1 solutions consistent with brown dwarf-mass secondaries with the Transiting Exoplanet Survey Satellite (TESS) Object of Interest (TOI) list to identify three candidate transiting brown dwarf systems. Ground-based precision radial velocity follow-up observations confirm that TOI-2533.01 is a transiting brown dwarf with orbiting TYC 2010-124-1 and that TOI-5427.01 is a transiting very low-mass star with orbiting UCAC4 515-012898. We validate TOI-1712.01 as a very low-mass star with transiting the primary in the hierarchical triple system BD+45 1593. Even after accounting for third light, TOI-1712.01 has a radius nearly a factor of 2 larger than predicted for isolated stars with similar properties. We propose that the intense instellation experienced by TOI-1712.01 diminishes the temperature gradient near its surface, suppresses convection, and leads to its inflated radius. Our analyses verify Gaia DR3 SB1 solutions in the low Doppler semiamplitude limit, thereby providing the foundation for future joint analyses of Gaia radial velocities and Kepler, K2, TESS, and PLAnetary Transits and Oscillations light curves for the characterization of transiting massive brown dwarfs and very low-mass stars.

Mass Ratio of Single-line Spectroscopic Binaries with Visual Orbits Using Bayesian Inference and Suitable Priors* *Based in part on observations obtained at the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, e Inovaçãoes (MCTI) da República Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at

We present orbital elements for 22 single-line binaries, nine of them studied for the first time, determined from a joint spectroscopic and astrometric solution. The astrometry is based on interferometric measurements obtained with the HRCam Speckle camera on the SOAR 4.1 m telescope at Cerro Pachon, Chile, supplemented with historical data. The spectroscopic observations were secured using Echelle spectrographs (FEROS, FIDEOS, and HARPS) at La Silla, Chile. A comparison of our orbital elements and systemic velocities with previous studies, including Gaia radial velocities, shows the robustness of our estimations. By adopting suitable priors of the trigonometric parallax and spectral type of the primary component, and using a Bayesian inference methodology developed by our group, we were able to estimate mass ratios for these binaries. Combining the present results with a previous study of other single-line binaries from our team, we present a pseudo mass-to-luminosity relationship based on 23 systems (46 stars) in the mass range 0.6 ≤ M ⊙ ≤ 2.5. We find a reasonable correspondence with a fiducial mass-to-luminosity relationship. We conclude that our methodology does allow us to derive tentative mass ratios for these types of binaries.

Detection of Multiple Phase Space Overdensities of GSE Stars by Orbit Integration

In N-body simulations, nearly radial mergers can form shell-like overdensities in the sky position and phase space (r − v r ) due to the combination of dynamical friction and tidal stripping. The merger event of Gaia-Sausage-Enceladus (GSE) has provided a unique opportunity to study the shells in the phase space. To search for them, we integrate the orbits of 5949 GSE-related halo K giants from the LAMOST survey and record their positions at all time intervals in the r − v r diagram. After the subtraction of a smoothed background, we find six significant and complete thin chevron-like overdensities. The apocenters r apo of stars in the six chevrons are around 6.75, 12.75, 18.75, 25.25, 27.25, and 30.25 kpc. These chevrons reveal the multiple pileups of GSE stars at different apocenters. The application of a different Milky Way mass M vir will change the opening angles of these chevrons, while leaving their apocenters almost unchanged. By comparing with a recent study of the phase space overdensities of local halo stars from the Gaia Radial Velocity Spectrometer survey, our results are more inclined to a medium M vir of 1012 M ⊙. The application of a nonaxisymmetric Galactic potential with a steadily rotating bar has a blurring effect on the appearance of these chevron-like overdensities, especially for the chevrons with r apo > 20 kpc.

Gaia Data Release 3

Context. The Gaia Radial Velocity Spectrometer (RVS) provides the unique opportunity of a spectroscopic analysis of millions of stars at medium resolution (λ/Δλ ∼ 11 500) in the near-infrared (845−872 nm). This wavelength range includes the Ca II infrared triplet (IRT) at 850.03, 854.44, and 866.45 nm, which is a good indicator of magnetic activity in the chromosphere of late–type stars. Aims. Here we present the method devised for inferring the Gaia stellar activity index from the analysis of the Ca II IRT in the RVS spectrum, together with its scientific validation. Methods. The Gaia stellar activity index is derived from the Ca II IRT excess equivalent width with respect to a reference spectrum, taking the projected rotational velocity (vsini) into account. We performed scientific validation of the Gaia stellar activity index by deriving a R′IRT index, which is largely independent of the photospheric parameters, and considering the correlation with the R′HK index for a sample of stars. A sample of well-studied pre-main-sequence (PMS) stars is considered to identify the regime in which the Gaia stellar activity index may be affected by mass accretion. The position of these stars in the colour–magnitude diagram and the correlation with the amplitude of the photometric rotational modulation is also scrutinised. Results.Gaia DR3 contains a stellar activity index derived from the Ca II IRT for some 2 × 106 stars in the Galaxy. This represents a ‘gold mine’ for studies on stellar magnetic activity and mass accretion in the solar vicinity. Three regimes of the chromospheric stellar activity are identified, confirming suggestions made by previous authors based on much smaller R′HK datasets. The highest stellar activity regime is associated with PMS stars and RS CVn systems, in which activity is enhanced by tidal interaction. Some evidence of a bimodal distribution in main sequence (MS) stars with Teff ≳ 5000 K is also found, which defines the two other regimes, without a clear gap in between. Stars with 3500 K ≲ Teff ≲ 5000 K are found to be either very active PMS stars or active MS stars with a unimodal distribution in chromospheric activity. A dramatic change in the activity distribution is found for Teff ≲ 3500 K, with a dominance of low activity stars close to the transition between partially- and fully convective stars and a rise in activity down into the fully convective regime.

The Period Distribution of Hot Jupiters Is Not Dependent on Host Star Metallicity

The probability that a Sun-like star has a close-orbiting giant planet (period ≲1 yr) increases with stellar metallicity. Previous work provided evidence that the period distribution of close-orbiting giant planets is also linked to metallicity, hinting that there two formation/evolution pathways for such objects, one of which is more probable in high-metallicity environments. Here, we check for differences in the period distribution of hot Jupiters (P < 10 days) as a function of host star metallicity, drawing on a sample of 232 transiting hot Jupiters and homogeneously derived metallicities from Gaia Data Release 3. We found no evidence for any metallicity dependence; the period distributions of hot Jupiters around metal-poor and metal-rich stars are indistinguishable. As a byproduct of this study, we provide transformations between metallicities from the Gaia Radial Velocity Spectrograph and from traditional high-resolution optical spectroscopy of main-sequence FGK stars.

Self-consistent models of our Galaxy

ABSTRACT A new class of models of stellar discs is introduced and used to build a self-consistent model of our Galaxy. The model is defined by the parameters that specify the action-based distribution functions (DFs) f(J) of four stellar discs (three thin-disc age cohorts and a thick disc), spheroidal bulge and spheroidal stellar and dark haloes. From these DFs plus a specified distribution of gas, we solve for the densities of stars and dark matter and the potential they generate. The principal observational constraints are the kinematics of stars with Gaia Radial Velocity Spectrometer (RVS) data and the density of stars in the column above the Sun. The model predicts the density and kinematics of stars and dark matter throughout the Galaxy, and suggests the structure of the dark halo prior to the infall of baryons. The code used to create the model is available on github.

The Poor Old Heart of the Milky Way

Our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population, which reflects the star formation and enrichment in the few most massive progenitors that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor “heart” of the Milky Way. We use information from Gaia Data Release 3, especially the XP spectra, to construct a sample of 2 million bright (G BP < 15.5 mag) giant stars within 30° of the Galactic center (GC) with robust [M/H] estimates, δ[M/H] ≲ 0.1. For ∼1.25 million stars we calculate orbits from Gaia Radial Velocity Spectrometer velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes ∼18,000 stars with −2.7 < [M/H] < −1.5, representing a stellar mass of ≳5 × 107 M ⊙. The spatial distribution of these [M/H] < −1.5 stars has a Gaussian extent of only around the GC, with most orbits confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Most stars show [α/Fe] enhancement and [Al/Fe]–[Mn/Fe] abundances expected for an origin in the more massive portions of the proto-Galaxy. Stars with [M/H] < −2 show no net rotation, whereas those with [M/H] ∼ −1 are rotation dominated. These central, metal-poor stars most likely predate the oldest disk population (τ age ≈ 12.5 Gyr), which implies that they formed at z ≳ 5, forging the proto-Milky Way.

A machine learning approach to photometric metallicities of giant stars

ABSTRACT Despite the advances provided by large-scale photometric surveys, stellar features – such as metallicity – generally remain limited to spectroscopic observations often of bright, nearby low-extinction stars. To rectify this, we present a neural network approach for estimating the metallicities and distances of red giant stars with 8-band photometry and parallaxes from Gaia EDR3 and the 2MASS and WISE surveys. The algorithm accounts for uncertainties in the predictions arising from the range of possible outputs at each input and from the range of models compatible with the training set (through drop-out). A two-stage procedure is adopted where an initial network to estimate photoastrometric parallaxes is trained using a large sample of noisy parallax data from Gaia EDR3 and then a secondary network is trained using spectroscopic metallicities from the APOGEE and LAMOST surveys and an augmented feature space utilizing the first-stage parallax estimates. The algorithm produces metallicity predictions with an average uncertainty of $\pm 0.19\, \mathrm{dex}$. The methodology is applied to stars within the Galactic bar/bulge with particular focus on a sample of 1.69 million objects with Gaia radial velocities. We demonstrate the use and validity of our approach by inspecting both spatial and kinematic gradients with metallicity in the Galactic bar/bulge recovering previous results on the vertical metallicity gradient (−0.528 ± 0.002 dex kpc−1) and the vertex deviation of the bar ($-21.29\pm 2.74\, \mathrm{deg}$).

Stars That Approach within One Parsec of the Sun: New and More Accurate Encounters Identified in Gaia Data Release 3

Close encounters of stars to the Sun could affect life on Earth through gravitational perturbations of comets in the Oort cloud or exposure to ionizing radiation. By integrating orbits through the Galactic potential, I identify which of 33 million stars in Gaia DR3 with complete phase space information come close to the Sun. 61 stars formally approach within 1 pc, although there is high confidence in only 42 (two thirds) of these, the rest being spurious measurements or (in) binary systems. Most of the stars will encounter within the past or future 6 Myr; earlier/later encounters are less common due to the magnitude limit of the Gaia radial velocities (RVs). Several close encountering stars are identified for the first time, and the encounter times, distances, and velocities of previously known close encounters are determined more precisely on account of the significantly improved precision of Gaia DR3 over earlier releases. The K7 dwarf Gl 710 remains the closest known encounter, with an estimated (median) encounter distance of 0.0636 pc (90% confidence interval 0.0595–0.0678 pc) to take place in 1.3 Myr. The new second closest encounter took place 2.8 Myr ago: this was the G3 dwarf HD 7977, now 76 pc away, which approached within less than 0.05 pc of the Sun with a probability of one third. The apparent close encounter of the white dwarf UPM J0812-3529 is probably spurious due to an incorrect RV in Gaia DR3.

A measurement of stellar surface gravity hidden in radial velocity differences of comoving stars

ABSTRACT The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of the pairs of comoving stars and wide binaries, we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the 5σ level for our full data set. Additionally, we study the subdominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass–luminosity relation and stellar models that predict the magnitude of convective blueshift.

Gaia-predicted brown dwarf detection rates around FGK stars in astrometry, radial velocity, and photometric transits

Context. After more than two decades of relevant radial velocity surveys, the current sample of known brown dwarfs (BDs) around FGK stars is only of the order of 100, limiting our understanding of their occurrence rate, properties, and formation. The ongoing ESA mission Gaia has already collected more than its nominal 5 years of mission data, and is expected to operate for up to 10 years in total. Its exquisite astrometric precision allows for the detection of (unseen) companions down to the Jupiter-mass level, allowing the efficient detection of large numbers of BDs. Additionally, its low-accuracy multi-epoch radial velocity measurements for GRVS &lt; 12 can provide additional detections or constraints for the more massive BDs, while a further small sample will have detectable transits in Gaia photometry. Aims. Using detailed simulations, we provide an assessment of the number of BDs that could be discovered using Gaia astrometry, radial velocity, and photometric transits around main sequence (V) and subgiant (IV) FGK host stars for the nominal five-year and extended ten-year mission. Methods. Using a robust Δχ2 statistic we analyse the BD companion detectability from the Besançon Galaxy population synthesis model for G = 10.5 − 17.5 mag, complemented by Gaia DR2 data for G &lt; 10.5, using the latest Gaia performance and scanning law, and literature-based BD-parameter distributions. Results. We report here only reliable detection numbers with Δχ2 &gt; 50 for a five-year mission, and those in square brackets are for a ten-year mission. For astrometry alone, we expect 28 000–42 000 [45 000–55 000] detections out to several hundred parsecs [up to more than a kiloparsec]. The majority of these have G ∼ 14 − 15 [14 − 16] and periods of greater than 200 d, extending up to the longest simulated periods of 5 yr. Gaia radial velocity time-series for GRVS &lt; 12 (G ≲ 12.7), should allow the detection of 830–1100 [1500–1900] BDs, most having orbital periods of &lt; 10 days, and being amongst the most massive BDs (55 − 80MJ), though several tens will extend down to the ‘desert’ and lowest BD masses. Systems with at least three photometric transits with S/N &gt; 3 are expected for 720–1100 [1400–2300] BDs, averaging at 4–5 [5–6] transits per source. The combination of astrometric and radial velocity detections should yield some 370–410 [870–950] detections. Perhaps 17–27 [35–56] BDs will have both transit and radial velocity detections, while both transits and astrometric detection will lead to a meagre 1–3 [4–6] detection(s). Conclusions. Though the above numbers are affected by ±50% uncertainty due to the uncertain occurrence rate and period distribution of BDs around FGK host stars, detections of BDs with Gaia will number in the tens of thousands, enlarging the current BD sample by at least two orders of magnitude, allowing us to investigate the BD fraction and orbital architectures as a function of host stellar parameters in greater detail than every before.

Kinematical Analysis of Substructure in the Southern Periphery of the Large Magellanic Cloud

Abstract We report the first 3D kinematical measurements of 88 stars in the direction of several recently discovered substructures in the southern periphery of the Large Magellanic Cloud (LMC) using a combination of Gaia proper motions and radial velocities from the APOGEE-2 survey. More specifically, we explore stars in assorted APOGEE-2 pointings in a region of the LMC periphery where various overdensities of stars have previously been identified in maps of stars from Gaia and DECam. By using a model of the LMC disk rotation, we find that a sizable fraction of the APOGEE-2 stars have extreme space velocities that are distinct from, and not a simple extension of, the LMC disk. Using N-body hydrodynamical simulations of the past dynamical evolution and interaction of the LMC and Small Magellanic Cloud (SMC), we explore whether the extreme-velocity stars may be accounted for as tidal debris created in the course of that interaction. We conclude that the combination of LMC and SMC debris produced from their interaction is a promising explanation, although we cannot rule out other possible origins, and that these new data should be used to constrain future simulations of the LMC–SMC interaction. We also conclude that many of the stars in the southern periphery of the LMC lie outside of the LMC plane by several kiloparsecs. Given that the metallicity of these stars suggests that they are likely of Magellanic origin, our results suggest that a wider exploration of the past interaction history of the Magellanic Clouds is needed.