Gaia Catalogue Research Articles

Gaia Focused Product Release: Sources from Service Interface Function image analysis

Context. Gaia’s readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri (ω Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Aims. Our aim is to improve the completeness of the Gaia source inventory in a very dense region in the sky, ω Cen. Methods. An adapted version of Gaia’s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. These source detections were clustered and assigned to new SIF CF or existing Gaia sources by Gaia s cross-match software. For the new sources, astrometry was calculated using the Astrometric Global Iterative Solution software, and photometry was obtained in the Gaia DR3 reference system. We validated the results by comparing them to the public Gaia DR3 catalogue and external Hubble Space Telescope data. Results. With this Focused Product Release, 526 587 new sources have been added to the Gaia catalogue in ω Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal Gaia sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. Conclusions. This first SIF CF data publication already adds great value to the Gaia catalogue. It demonstrates what to expect for the fourth Gaia catalogue, which will contain additional sources for all nine SIF CF regions.

Gaia Focused Product Release: Radial velocity time series of long-period variables

Context. The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity data planned with Data Release 4, this Focused Product Release (FPR) provides radial-velocity time series for a selection of LPV candidates with high-quality observations. Aims. We describe the production and content of the Gaia catalog of LPV radial-velocity time series, and the methods used to compute the variability parameters published as part of the Gaia FPR. Methods. Starting from the DR3 catalog of LPV candidates, we applied several filters to construct a sample of sources with high-quality radial-velocity measurements. We modeled their radial-velocity and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the radial-velocity period and at least one of the G, GBP, or GRP photometric periods. Results. The catalog includes radial-velocity time series and variability parameters for 9614 sources in the magnitude range 6 ≲ G/mag ≲ 14, including a flagged top-quality subsample of 6093 stars whose radial-velocity periods are fully compatible with the values derived from the G, GBP, and GRP photometric time series. The radial-velocity time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great majority of the sources (88%) as genuine LPV candidates, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% of the catalog consists of candidate ellipsoidal binaries. Quality checks against radial velocities available in the literature show excellent agreement. We provide some illustrative examples and cautionary remarks. Conclusions. The publication of radial-velocity time series for almost ten thousand LPV candidates constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog.

Methodology for obtaining the relative orbit and individual masses of Gaia astrometric binaries

The recent Gaia Data Release 3 has revealed a catalogue of more than eight hundred thousand binary systems. The release provides orbital solutions for half of the systems, with the majority of them being unresolved astrometric binaries. However, some astrophysical parameters are still unknown for most of them, such as the spectral type and the mass of each companion, since they can only be derived from the relative orbit and spectroscopic data. The purpose of this work is to develop a methodology that would allow us to obtain those fundamental stellar parameters, along with those related to the geometry and the ephemeris of the system, to find out whether it can be optically resolved. To obtain precise values for each component, we proposed an analytic algorithm to estimate the only two possible relative orbits and pairs of masses of main sequence (MS) astrometric binaries using all the available astrometric, photometric, and spectroscopic data from Gaia DR3. In some cases, it is possible to select the solution that is more aligned with the rest of the data. We deduced two possible values for the individual absolute magnitudes, masses, and effective temperatures for each binary, as well as the size of the telescope necessary to resolve their components. We present the workflow of our algorithm applied to the Ephemeris, Stellar Masses, and relative ORbits from GAia (ESMORGA) catalogue, along with the individual masses, absolute magnitudes, and effective temperatures derived for 49\,530 binaries.

Estimating Distances from Parallaxes. VI. A Method for Inferring Distances and Transverse Velocities from Parallaxes and Proper Motions Demonstrated on Gaia Data Release 3

The accuracy of stellar distances inferred purely from parallaxes degrades rapidly with distance. Proper motion measurements, when combined with some idea of typical velocities, provide independent information on stellar distances. Here, I build a direction- and distance-dependent model of the distribution of stellar velocities in the Galaxy, then use this together with parallaxes and proper motions to infer kinegeometric distances and transverse velocities for stars in Gaia DR3. Using noisy simulations, I assess the performance of the method and compare its accuracy to purely parallax-based (geometric) distances. Over the whole Gaia catalog, kinegeometric distances are on average 1.25 times more accurate than geometric ones. This average masks a large variation in the relative performance, however. Kinegeometric distances are considerably better than geometric ones beyond several kpc, for example. On average, kinegeometric distances can be measured to an accuracy of 19% and velocities () to 16 km s−1 (median absolute deviations). In Gaia DR3, kinegeometric distances are smaller than geometric ones on average for distant stars, but the pattern is more complex in the bulge and disk. With the much more accurate proper motions expected in Gaia DR5, a further improvement in the distance accuracy by a factor of (only) 1.35 on average is predicted (with kinegeometric distances still 1.25 times more accurate than geometric ones). The improvement attained from proper motions is limited by the width of the velocity prior, in a way that the improvement from better parallaxes is not limited by the width of the distance prior.

Statistical analysis of the gravitational anomaly in Gaia wide binaries

ABSTRACT The exploration of the low acceleration a < a0 regime, where a0 = 1.2 × 10−10 m s−2 is the acceleration scale of MOND around which gravitational anomalies at galactic scale appear, has recently been extended to the much smaller mass and length scales of local wide binaries thanks to the availability of the Gaia catalogue. Statistical methods to test the underlying structure of gravity using large samples of such binary stars and dealing with the necessary presence of kinematic contaminants in such samples have also been presented. However, an alternative approach using binary samples carefully selected to avoid any such contaminants, and consequently much smaller samples, has been lacking a formal statistical development. In the interest of having independent high-quality checks on the results of wide binary gravity tests, we here develop a formal statistical framework for treating small, clean, wide binary samples in the context of testing modifications to gravity of the form G → γG. The method is validated through extensive tests with synthetic data samples, and applied to recent Gaia DR3 binary star observational samples of relative velocities and internal separations on the plane of the sky, v2D and r2D, respectively. Our final results for a high acceleration r2D < 0.01 pc region are of γ = 1.000 ± 0.096, in full accordance with Newtonian expectations. For a low acceleration r2D > 0.01 pc region, however, we obtain γ = 1.5 ± 0.2, inconsistent with the Newtonian value of γ = 1 at a 2.6σ level, and much more indicative of MOND AQUAL predictions of close to γ = 1.4.

Blending from binarity in microlensing searches toward the Large Magellanic Cloud

Context. Studies of gravitational microlensing effects require the estimation of their detection efficiency as soon as one wants to quantify the massive compact objects along the line of sight of source targets. This is particularly important for setting limits on the contribution of massive compact objects to the Galactic halo. These estimates of detection efficiency must not only account for the blending effects of accidentally superimposed sources in crowded fields, but also for possible mixing of light from stars belonging to multiple gravitationally bound stellar systems. Aims. Until now, only blending due to accidental alignment of stars had been studied, in particular as a result of high-resolution space images. In this paper, we address the impact of unresolved binary sources that are physically gravitationally bound and not accidentally aligned, in the case of microlensing detection efficiencies toward the Large Magellanic Cloud (LMC). Methods. We used the Gaia catalog of nearby stars to constrain the local binarity rate, which we extrapolated to the distance of the LMC. Then we estimated an upper limit to the impact of this binarity on the detection efficiency of microlensing effects, as a function of lens mass. Results. We find that a maximum of 6.2% of microlensing events on LMC sources due to halo lenses heavier than 30 M⊙ could be affected as a result of the sources belonging to unresolved binary systems. This number is the maximum fraction of events for which the source is a binary system separated by about one angular Einstein radius or more in a configuration where light-curve distortion could affect the efficiency of some detection algorithms. For events caused by lighter lenses on LMC sources, our study shows that the chances of blending effects by binary systems is likely to be higher and should be studied in more detail to improve the accuracy of efficiency calculations.

Dynamical Evolution of White Dwarfs in Triples in the Era of Gaia

The Gaia mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple-stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orbit about the inner binary. In such a system, the two orbits torque each other via the eccentric Kozai–Lidov mechanism, which can alter the orbital configuration of the inner binary. We simulate thousands of triple-stellar systems for over 10 Gyr, tracking gravitational interactions, tides, general relativity, and stellar evolution up to their WD fate. As demonstrated here, three-body dynamics coupled with stellar evolution is a critical channel to form tight WD binaries or merge a WD binary. Among these triples, we explore their manifestations as cataclysmic variables, Type Ia supernovae, and gravitational-wave events. The simulated systems are then compared to a sample of WD triples selected from the Gaia catalog. We find that including the effect of mass-loss-induced kicks is crucial for producing a distribution of the inner binary–tertiary separations that is consistent with Gaia observations. Lastly, we leverage this consistency to estimate that, at minimum, 30% of solar-type stars in the local 200 pc were born in triples.

The Potential of Detecting Nearby Terrestrial Planets in the HZ with Different Methods

Terrestrial planets in the habitable zone (hereafter HZ) around nearby stars are of great interest and provide a good sample to characterize their habitability. In this paper, we collect a nearby star catalog within 20 pc according to the Gaia Catalog of Nearby Stars, complete the physical parameters of the stars, and select stars that are not brown or white dwarfs. After selection, a sample of 2234 main-sequence stars is used to estimate the extended HZ. Then we inject Earth-like planets into the extended HZ around each star and calculate the signals with four methods, i.e., velocity amplitude for radial velocity, transit probability and depth for transit, stellar displacements for astrometry, and contrast and angular separation for imaging. Considering a typical noise model based on classic instruments, e.g., ESPRESSO, Kepler, Gaia, HabEx, and LIFE, we predict the highest possible detection number of Earth-like planets via different methods in the best-case hypothetical scenario. According to this, we conclude that both astrometry and imaging have the potential to detect nearby Earth-like planets around G-type stars, while radial velocity has the potential to detect 2% of nearby Earth-like planets around M stars with a precision of 0.2 m s−1. Our work also provides the precision requirements for future missions to reveal the nearby Earth-like planet in the HZ.

Estimating the selection function of Gaia DR3 subsamples

Context. Understanding the intricacies behind the presence and absence of sources in an astronomical catalogue is crucial for the accurate interpretation of astronomical data. In particular, for the multi-dimensional Gaia data, filters and cuts on different parameters or measurements introduce a selection function that may unintentionally alter scientific conclusions in subtle ways. Aims. We aim to develop a methodology to estimate the selection function for different subsamples of stars in the Gaia catalogue. Methods. Comparing the number of stars in a given subsample to that in the overall Gaia catalogue provides an estimate of the subsample membership probability as a function of sky position, magnitude, and colour. The method used to make this estimate must differentiate the stochastic absence of subsample stars from selection effects. When multiplied with the overall Gaia catalogue selection function, this provides the total selection function of the subsample. Results. We present our new method for estimating the selection function by applying it to the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. Conclusions. The method we have developed estimates the selection function of the stars present in a subsample of Gaia data, given that the subsample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package.

Astrometry via close approach events: applications to main-belt asteroid (702) Alauda

ABSTRACT The release of Gaia catalogue is revolutionary to the astronomy of Solar system objects. After some effects such as atmospheric refraction and CCD geometric distortion have been taken into account, the astrometric precision for ground-based telescopes can reach the level of tens of milli-arcseconds (mas). If an object approaches a reference star in a small relative angular distance (less than 100 arcsec), which is called close approach event in this work, the relative positional precision between the object and reference star will be further improved since the systematic effects of atmospheric turbulence and local telescope optics can be reduced. To obtain the precise position of a main-belt asteroid in an close approach event, a second-order angular velocity model with time is supposed in the sky plane. By fitting the relationship between the relative angular distance and observed time, we can derive the time of maximum approximation and calculate the corresponding position of the asteroid. In practice, five nights’ CCD observations including 15 close approach events of main-belt asteroid (702) Alauda are taken for testing by the 1 m telescope at Yunnan Observatory, China. Compared with conventional solutions, our results show that the positional precision significantly improves, which reaches better than 4 mas, and 1 mas in the best case when referenced for the Jet Propulsion Laboratory (JPL) ephemeris in both right ascension and declination.

Internal kinematics of Gaia DR3 wide binaries: anomalous behaviour in the low acceleration regime

ABSTRACT The Gaia eDR3 catalogue has recently been used to study statistically the internal kinematics of wide binary populations using relative velocities of the two component stars, ΔV, total binary masses, mB, and separations, s. For s ≳ 0.01 pc, these binaries probe the low-acceleration a ≲ 2a0 regime where gravitational anomalies usually attributed to dark matter are observed in the flat rotation curves of spiral galaxies, where a0 ≈ 1.2 × 10−10 m s−2 is the acceleration scale of MOND. Such experiments test the degree of generality of these anomalies, by exploring the same acceleration regime using independent astronomical systems of vastly smaller mass and size. A signal above Newtonian expectations has been observed when a ≲ 2a0, alternatively interpreted as evidence of a modification of gravity, or as due to kinematic contaminants; undetected stellar components, unbound encounters, or spurious projection effects. Here I take advantage of the enhanced DR3 Gaia catalogue to perform a more rigorous study of the internal kinematics of wide binaries than what has previously been possible. Internally determined Gaia stellar masses and estimates of binary probabilities for each star using spectroscopic information, together with a larger sample of radial velocities, allow for a significant improvement in the analysis and careful exclusion of possible kinematic contaminants. Resulting ΔV scalings accurately tracing Newtonian expectations for the high acceleration regime, but markedly inconsistent with these expectations in the low acceleration one, are obtained. A non-Newtonian low acceleration phenomenology is thus confirmed.

Revisiting the mass of open clusters with Gaia data

ABSTRACT The publication of the Gaia catalogue and improvements in methods to determine memberships and fundamental parameters of open clusters has led to major advances in recent years. However, important parameters such as the masses of these objects, although being studied mostly in some isolated cases, have not been addressed in large homogeneous samples based on Gaia data, taking into account details such as binary fractions. Consequently, relevant aspects such as the existence of mass segregation were not adequately studied. Within this context, in this work, we introduce a new method to determine individual stellar masses, including an estimation for the ones in binary systems. This method allows us to study the mass of open clusters, as well as the mass functions of the binary star populations. We validate the method and its efficiency and characterize uncertainties using a grid of synthetic clusters with predetermined parameters. We highlight the application of the method to the Pleiades cluster, showing that the results obtained agree with the current consensus in the literature as well as recent Gaia data. We then applied the procedure to a sample of 773 open clusters with fundamental parameters determined using Gaia Early Data Release 3 data, obtaining their masses. Subsequently, we investigated the relation between the masses and other fundamental parameters of the clusters. Among the results, we found no significant evidence that clusters in our sample lose and segregate mass with age.

A Hunting Expedition For High-Order Hierarchies

ABSTRACT Stellar hierarchical systems of high order containing more than three stars are rare and fascinating objects; their discovery and study highlight still unknown aspects of star formation and early evolution. We matched eclipsing binaries discovered by Transiting Exoplanet Survey Satellite with a Gaia catalogue of wide binaries and selected candidate quadruple (or higher order) systems based on excessive astrometric noise. A subset of 192 southern candidates located within 500 pc was observed by speckle interferometry, and we resolved for the first time 50 close pairs, confirming their high-order (from four to five components) multiplicity. These observations are reported, and some remarkable hierarchical systems are discussed.

The magnetic fields around the cometary globules L328, L323, and L331

ABSTRACT This work presents the magnetic field geometry in a complex of three cometary (with head–tail morphology) globules, namely LDN 323, LDN 328, and LDN 331, using R-band polarization measurements of background stars. These observations were combined with a Planck sky survey to study the large-scale morphology of the magnetic fields in the region. The distances of the target stars were adopted from the Gaia catalogue. The variation of the degree of polarization and polarization position angle with distances of stars is analysed. The field geometry is mostly found to follow the cometary shape of the cloud, with some randomness at certain locations. For studying the correlation between cloud morphology and magnetic field orientations, a modified version of the histogram of relative orientation (HRO) analysis was employed.

Gaia Data Release 3

Context. One of the rarest types of variability is the phenomenon of gravitational microlensing, a transient brightening of a background star due to an intervening lensing object. Microlensing is a powerful tool for studying the invisible or otherwise undetectable populations in the Milky Way, including planets and black holes. Aims. We describe the first Gaia catalogue of candidate microlensing events, give an overview of its content, and discuss its validation. Methods. The catalogue of Gaia microlensing events was composed by analysing the light curves of around 2 billion sources of Gaia DR3 from all over the sky covering 34 months, between 2014 and 2017. Results. We present 363 Gaia microlensing events and discuss their properties. Of these, 90 have never been reported before and have not been discovered by other surveys. The contamination of the catalogue is assessed to 0.6%–1.7%.

Gaia Data Release 3

Context. The third Gaia Data Release covers 34 months of data and includes the second Gaia catalogue of long-period variables (LPVs), with G variability amplitudes larger than 0.1 mag (5–95% quantile range). Aims. The paper describes the production and content of the second Gaia catalogue of LPVs and the methods we used to compute the published variability parameters and identify C-star candidates. Methods. We applied various filtering criteria to minimise contamination from variable star types other than LPVs. The period and amplitude of the detected variability were derived from model fits to the G-band light curve wherever possible. C stars were identified using their molecular signature in the low-resolution RP spectra. Results. The catalogue contains 1 720 558 LPV candidates, including 392 240 stars with published periods (ranging from 35 to ∼1000 days) and 546 468 stars classified as C-star candidates. Comparison with literature data (OGLE and ASAS-SN) leads to an estimated completeness of 80%. The recovery rate is about 90% for the most regular stars (typically miras) and 60% for SRVs and irregular stars. At the same time, the number of known LPVs is increased by a factor of 6 with respect to literature data for amplitudes larger than 0.1 mag in G, and the contamination is estimated to be below 2%. Our C-star classification, based on solid theoretical arguments, is consistent with spectroscopically identified C stars in the literature. Caution must be taken in crowded regions, however, where the signal-ro-noise ratio of the RP spectra can become very low, or if the source is reddened by some kind of extinction. The quality and potential of the catalogue are illustrated by presenting and discussing LPVs in the solar neighbourhood, in globular clusters, and in galaxies of the Local Group. Conclusions. This is the largest all-sky LPVs catalogue to date. The photometric depth reaches G = 20 mag. This is a unique dataset for research into the late stages of stellar evolution.

An interferometric SETI observation of Kepler-111 b

ABSTRACT The application of very long baseline interferometry (VLBI) to the search for extraterrestrial intelligence (SETI) has been limited to date, despite the technique offering many advantages over traditional single-dish SETI observations. In order to further develop interferometry for SETI, we used the European VLBI network (EVN) at 21 cm to observe potential secondary phase calibrators in the Kepler field. Unfortunately, no secondary calibrators were detected. However, a VLBA primary calibrator in the field, J1926+4441, offset only ∼1.88 arcmin from a nearby exoplanet Kepler-111 b, was correlated with high temporal $\left(0.25 \, \rm {s}\right)$ and spectral $\left(16384 \times 488\, \rm {Hz \ channels}\right)$ resolution. During the analysis of the high-resolution data, we identified a spectral feature that was present in both the auto and cross-correlation data with a central frequency of 1420.424 ± 0.0002 MHz and a width of 0.25 MHz. We demonstrate that the feature in the cross-correlations is an artefact in the data, associated with a significant increase in each telescope’s noise figure due to the presence of H i in the beam. This would typically go unnoticed in data correlated with standard spectral resolution. We flag (excluded from the subsequent analysis) these channels and phase rotate the data to the location of Kepler-111 b aided by the Gaia catalogue and search for signals with $\rm {SNR}\gt 7$ . At the time of our observations, we detect no transmitters with an equivalent isotropically radiated power ≳4 × 1015 W.

A Catalog of Nearby Accelerating Star Candidates in Gaia DR3

We describe a new catalog of accelerating star candidates with Gaia G ≤ 17.5 mag and distances d ≤ 100 pc. Designated as the Gaia Nearby Accelerating Star Catalog (GNASC), it contains 29,684 members identified using a supervised machine-learning algorithm trained on the Hipparcos–Gaia Catalog of Accelerations (HGCA), Gaia Data Release 2, and Gaia Early Data Release 3. We take advantage of the difference in observation timelines between the two Gaia catalogs and information about the quality of the astrometric modeling based on the premise that acceleration will correlate with astrometric uncertainties. Catalog membership is based on whether constant proper motion over three decades can be ruled out at high confidence (greater than 99.9%). Test data suggest that catalog members each have a 68% likelihood of true astrometric acceleration; subsets of the catalog perform even better, with the likelihood exceeding 85%. We compare the GNASC with Gaia Data Release 3 and its table of stars for which acceleration is detected at high confidence based on precise astrometric fits. Our catalog, derived without this information, captures over 96% of the sources in the table that meet our selection criteria. In addition, the GNASC contains bright, nearby candidates that were not in the original Hipparcos survey, including members of known binary systems as well as stars with companions yet to be identified. It thus extends the HGCA and demonstrates the potential of the machine-learning approach for discovering hidden partners of nearby stars in future astrometric surveys.

Exploring Thousands of Nearby Hierarchical Systems with Gaia and Speckle Interferometry

There should be about 10,000 stellar hierarchical systems within 100 pc with primary stars more massive than 0.5 M ⊙, and a similar amount of less-massive hierarchies. A list of 8000 candidate multiples is derived from wide binaries found in the Gaia Catalog of Nearby Stars where one or both components have excessive astrometric noise or other indicators of inner subsystems. A subset of 1243 southern candidates were observed with high angular resolution at the 4.1 m Southern Astrophysical Research Telescope, and 503 new pairs with separations from 0.″03 to 1″ were resolved. These data allow estimation of the inner mass ratios and periods, and help to quantify the ability of Gaia to detect close pairs. Another 621 hierarchies with known inner periods come from the Gaia catalog of astrometric and spectroscopic orbits. These two nonoverlapping groups, combined with existing ground-based data, bring the total number of known nearby hierarchies to 2754, reaching a completeness of ∼22% for stars above 0.5 M ⊙. Distributions of their periods and mass ratios are briefly discussed, and the prospects of further observations are outlined.

Photometric Catalogue for Space and Ground Night-Time Remote-Sensing Calibration: RGB Synthetic Photometry from Gaia DR3 Spectrophotometry

Recent works have made strong efforts to produce standardised photometry in RGB bands. For this purpose, we carefully defined the transmissivity curves of RGB bands and defined a set of standard sources using the photometric information present in Gaia EDR3. This work aims not only to significantly increase the number and accuracy of RGB standards but also to provide, for the first time, reliable uncertainty estimates using the BP and RP spectrophotometry published in Gaia DR3 instead of their integrated photometry to predict RGB photometry. Furthermore, this method allows including calibrated sources regardless of how they are affected by extinction, which was a major shortcoming of previous work. The RGB photometry is synthesised from the Gaia BP and RP low-resolution spectra by directly using their set of coefficients multiplied with some basis functions provided in the Gaia catalogue for all sources published in Gaia DR3. The output synthetic magnitudes are compared with the previous catalogue of RGB standards available.