Context.Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for half a billion stars. In this work, we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. Aims. We produce homogeneous samples of stars with high-quality astrophysical parameters by exploiting Gaia DR3, while focusing on many regimes across the Hertzsprung-Russell (HR) diagram; spectral types OBA, FGKM, and ultracool dwarfs (UCDs). We also focus on specific subsamples of particular interest to the community: solar analogues, carbon stars, and the Gaia spectrophotometric standard stars (SPSS). Methods. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples of the stars of interest. We validate our results using the Gaia catalogue itself and by comparison with external data. Results. We produced six homogeneous samples of stars with high-quality astrophysical parameters across the HR diagram for the community to exploit. We first focus on three samples that span a large parameter space: young massive disc stars (OBA; about 3 Million), FGKM spectral type stars (about 3 Million), and UCDs (about 20 000). We provide these sources along with additional information (either a flag or complementary parameters) as tables that are made available in the Gaia archive. We also identify 15 740 bone fide carbon stars and 5863 solar analogues, and provide the first homogeneous set of stellar parameters of the SPSS sample. We demonstrate some applications of these samples in different astrophysical contexts. We use a subset of the OBA sample to illustrate its usefulness in analysing the Milky Way rotation curve. We then use the properties of the FGKM stars to analyse known exoplanet systems. We also analyse the ages of some unseen UCD-companions to the FGKM stars. We additionally predict the colours of the Sun in various passbands (Gaia, 2MASS, WISE) using the solar-analogue sample. Conclusions.Gaia DR3 contains a wealth of new high-quality astrophysical parameters for the community to exploit.@en

Context. The Gaia DR3 catalogue contains, for the first time, about 800 000 solutions with either orbital elements or trend parameters for astrometric, spectroscopic, and eclipsing binaries, and combinations of these three. Aims. With this paper, we aim to illustrate the huge potential of this large non-single-star catalogue. Methods. Using the orbital solutions and models of the binaries, we have built a catalogue of tens of thousands of stellar masses or lower limits thereof, some with consistent flux ratios. Properties concerning the completeness of the binary catalogues are discussed, statistical features of the orbital elements are explained, and a comparison with other catalogues is performed. Results. Illustrative applications are proposed for binaries across the Hertzsprung-Russell Diagram (HRD). Binarity is studied in the giant branch and a search for genuine spectroscopic binaries among long-period variables is performed. The discovery of new EL CVn systems illustrates the potential of combining variability and binarity catalogues. Potential compact object companions are presented, mainly white dwarf companions or double degenerates, but one candidate neutron star is also found. Towards the bottom of the main sequence, the orbits of previously suspected binary ultracool dwarfs are determined and new candidate binaries are discovered. The long awaited contribution of Gaia to the analysis of the substellar regime shows the brown dwarf desert around solar-type stars using true rather than minimum masses, and provides new important constraints on the occurrence rates of substellar companions to M dwarfs. Several dozen new exoplanets are proposed, including two with validated orbital solutions and one super-Jupiter orbiting a white dwarf, all being candidates requiring confirmation. Besides binarity, higher order multiple systems are also found. Conclusions. By increasing the number of known binary orbits by more than one order of magnitude, Gaia DR3 will provide a rich reservoir of dynamical masses and an important contribution to the analysis of stellar multiplicity.@en

Context. With the most recent Gaia data release, the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, and more than 11 million variable stars are identified. Aims. Using the astrophysical parameters and variability classifications provided in Gaia DR3, we selected various stellar populations to explore and identify non-axisymmetric features in the disc of the Milky Way in configuration and velocity space. Methods. Using more about 580 000 sources identified as hot OB stars, together with 988 known open clusters younger than 100 Myr, we mapped the spiral structure associated with star formation 4'5 kpc from the Sun. We selected over 2800 Classical Cepheids younger than 200 Myr that show spiral features extending as far as 10 kpc from the Sun in the outer disc. We also identified more than 8.7 million sources on the red giant branch (RGB), of which 5.7 million have line-of-sight velocities. This later sample allows the velocity field of the Milky Way to be mapped as far as 8 kpc from the Sun, including the inner disc. Results. The spiral structure revealed by the young populations is consistent with recent results using Gaia EDR3 astrometry and source lists based on near-infrared photometry, showing the Local (Orion) Arm to be at least 8 kpc long, and an outer arm consistent with what is seen in HI surveys, which seems to be a continuation of the Perseus arm into the third quadrant. The subset of RGB stars with velocities clearly reveals the large-scale kinematic signature of the bar in the inner disc, as well as evidence of streaming motions in the outer disc that might be associated with spiral arms or bar resonances. A local comparison of the velocity field of the OB stars reveals similarities and differences with the RGB sample. Conclusions. This cursory study of Gaia DR3 data shows there is a rich bounty of kinematic information to be explored more deeply, which will undoubtedly lead us to a clearer understanding of the dynamical nature of the non-axisymmetric structures of the Milky Way.@en

Context. Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame isdefined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. Aims. We describe the construction of Gaia-CRF3 and its properties in terms of the distributions in magnitude, colour, and astrometricquality. Methods. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasi-stellar objects (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio–optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). Results. Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13–21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 μas yr−1 on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but it exceeds 4 mas in either coordinate for 127 sources. We outline the future of Gaia-CRF in the next Gaia data releases. Appendices give further details on the external catalogues used, how to extract information about the Gaia-CRF3 sources, potential (Galactic) confusion sources, and the estimation of the spin and orientation of an astrometric solution.@en

Context. The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators (M'a ¥''1.3'M') of spectral types O, B, A, or F, known as β Cep, slowly pulsating B (SPB), δ Sct, and γ Dor stars. These stars are often multi-periodic and display low amplitudes, making them challenging targets to analyse with sparse time series. Aims. We investigate the extent to which the sparse Gaia DR3 data can be used to detect OBAF-type pulsators and discriminate them from other types of variables. We aim to probe the empirical instability strips and compare them with theoretical predictions. The most populated variability class is that of the δ Sct variables. For these stars, we aim to confirm their empirical period-luminosity (PL) relation, and verify the relation between their oscillation amplitude and rotation. Methods. All datasets used in this analysis are part of the Gaia DR3 data release. The photometric time series were used to perform a Fourier analysis, while the global astrophysical parameters necessary for the empirical instability strips were taken from the Gaia DR3 gspphot tables, and the v sin i data were taken from the Gaia DR3 esphs tables. The δ'Sct PL relation was derived using the same photometric parallax method as the one recently used to establish the PL relation for classical Cepheids using Gaia data. Results. We show that for nearby OBAF-type pulsators, the Gaia DR3 data are precise and accurate enough to pinpoint them in the Hertzsprung-Russell (HR) diagram. We find empirical instability strips covering broader regions than theoretically predicted. In particular, our study reveals the presence of fast rotating gravity-mode pulsators outside the strips, as well as the co-existence of rotationally modulated variables inside the strips as reported before in the literature. We derive an extensive period'luminosity relation for δ Sct stars and provide evidence that the relation features different regimes depending on the oscillation period. We demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of δ Sct stars. Conclusions. The Gaia DR3 time-series photometry already allows for the detection of the dominant (non-)radial oscillation mode in about 100 000 intermediate- and high-mass dwarfs across the entire sky. This detection capability will increase as the time series becomes longer, allowing the additional delivery of frequencies and amplitudes of secondary pulsation modes.@en

Context. We present the third data release of the European Space Agency's Gaia mission, Gaia DR3. This release includes a large variety of new data products, notably a much expanded radial velocity survey and a very extensive astrophysical characterisation of Gaia sources. Aims. We outline the content and the properties of Gaia DR3, providing an overview of the main improvements in the data processing in comparison with previous data releases (where applicable) and a brief discussion of the limitations of the data in this release. Methods. The Gaia DR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. Results. The Gaia DR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, GBP, and GRP pass-bands already present in the Early Third Data Release, Gaia EDR3. Gaia DR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges GRVS'<'14 and 3100'<'Teff'<'14'500, have new determinations of their mean radial velocities based on data collected by Gaia. We provide GRVS magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The Gaia DR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BP/RP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. Gaia DR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some 800'000 astrometric, spectroscopic and eclipsing binaries. More than 150'000 Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BP/RP spectral data are published for about 60 000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey, consisting of the photometric time series for all sources located in a 5.5 degree radius field centred on the Andromeda galaxy. Conclusions. This data release represents a major advance with respect to Gaia DR2 and Gaia EDR3 because of the unprecedented quantity, quality, and variety of source astrophysical data. To date this is the largest collection of all-sky spectrophotometry, radial velocities, variables, and astrophysical parameters derived from both low- and high-resolution spectra and includes a spectrophotometric and dynamical survey of SSOs of the highest accuracy. The non-single star content surpasses the existing data by orders of magnitude. The quasar host and galaxy light profile collection is the first such survey that is all sky and space based. The astrophysical information provided in Gaia DR3 will unleash the full potential of Gaia's exquisite astrometric, photometric, and radial velocity surveys.@en

Context. The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was derived from measurements obtained by means of the Blue and Red photometers (BP/RP), which were binned in 16 discrete wavelength bands. For every spectrum, the DR3 also contains additional information about the data quality for each band. Aims. We describe the processing of the Gaia spectral data of SSOs, explaining both the criteria used to select the subset of asteroid spectra published in Gaia DR3, and the different steps of our internal validation procedures. In order to further assess the quality of Gaia SSO reflectance spectra, we carried out external validation against SSO reflectance spectra obtained from ground-based and space-borne telescopes and available in the literature; we present our validation approach. Methods. For each selected SSO, an epoch reflectance was computed by dividing the calibrated spectrum observed by the BP/RP at each transit on the focal plane by the mean spectrum of a solar analogue. The latter was obtained by averaging the Gaia spectral measurements of a selected sample of stars known to have very similar spectra to that of the Sun. Finally, a mean of the epoch reflectance spectra was calculated in 16 spectral bands for each SSO. Results.Gaia SSO reflectance spectra are in general agreement with those obtained from a ground-based spectroscopic campaign specifically designed to cover the same spectral interval as Gaia and mimic the illumination and observing geometry characterising Gaia SSO observations. In addition, the agreement between Gaia mean reflectance spectra and those available in the literature is good for bright SSOs, regardless of their taxonomic spectral class. We identify an increase in the spectral slope of S-type SSOs with increasing phase angle. Moreover, we show that the spectral slope increases and the depth of the 1 μm absorption band decreases for increasing ages of S-type asteroid families. The latter can be interpreted as proof of progressive ageing of S-type asteroid surfaces due to their exposure to space weathering effects.@en

Gaia Data Release 3 provides novel flux-calibrated low-resolution spectrophotometry for '220 million sources in the wavelength range 330 nm ≤ λ ≤ 1050 nm (XP spectra). Synthetic photometry directly tied to a flux in physical units can be obtained from these spectra for any passband fully enclosed in this wavelength range. We describe how synthetic photometry can be obtained from XP spectra, illustrating the performance that can be achieved under a range of different conditions - for example passband width and wavelength range - as well as the limits and the problems affecting it. Existing top-quality photometry can be reproduced within a few per cent over a wide range of magnitudes and colour, for wide and medium bands, and with up to millimag accuracy when synthetic photometry is standardised with respect to these external sources. Some examples of potential scientific application are presented, including the detection of multiple populations in globular clusters, the estimation of metallicity extended to the very metal-poor regime, and the classification of white dwarfs. A catalogue providing standardised photometry for 2.2×108sources in several wide bands of widely used photometric systems is provided (Gaia Synthetic Photometry Catalogue; GSPC) as well as a catalogue of '105 white dwarfs with DA/non-DA classification obtained with a Random Forest algorithm (Gaia Synthetic Photometry Catalogue for White Dwarfs; GSPC-WD).@en

The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G 21 mag. Due to the nature of the Gaia onboard-selection algorithms, these are mostly point-source-like objects. Using data provided by the satellite, we have identified quasar and galaxy candidates via supervised machine learning methods, and estimate their redshifts using the low resolution BP/RP spectra. We further characterise the surface brightness profiles of host galaxies of quasars and of galaxies from pre-defined input lists. Here we give an overview of the processing of extragalactic objects, describe the data products in Gaia DR3, and analyse their properties. Two integrated tables contain the main results for a high completeness, but low purity (50-70%), set of 6.6 million candidate quasars and 4.8 million candidate galaxies. We provide queries that select purer sub-samples of these containing 1.9 million probable quasars and 2.9 million probable galaxies (both 95% purity). We also use high quality BP/RP spectra of 43 thousand high probability quasars over the redshift range 0.05-4.36 to construct a composite quasar spectrum spanning restframe wavelengths from 72 1000 nm.@en

Context. Diffuse interstellar bands (DIBs) are common interstellar absorption features in spectroscopic observations but their origins remain unclear. DIBs play an important role in the life cycle of the interstellar medium (ISM) and can also be used to trace Galactic structure. Aims. Here, we demonstrate the capacity of the Gaia-Radial Velocity Spectrometer (RVS) in Gaia DR3 to reveal the spatial distribution of the unknown molecular species responsible for the most prominent DIB at 862 nm in the RVS passband, exploring the Galactic ISM within a few kiloparsecs from the Sun. Methods. The DIBs are measured within the GSP-Spec module using a Gaussian profile fit for cool stars and a Gaussian process for hot stars. In addition to the equivalent widths and their uncertainties, Gaia DR3 provides their characteristic central wavelength, width, and quality flags. Results. We present an extensive sample of 476 117 individual DIB measurements obtained in a homogeneous way covering the entire sky. We compare spatial distributions of the DIB carrier with interstellar reddening and find evidence that DIB carriers are present in a local bubble around the Sun which contains nearly no dust. We characterised the DIB equivalent width with a local density of 0.19 ± 0.04 kpc1 and a scale height of 98.60 8.46+11.10 pc. The latter is smaller than the dust scale height, indicating that DIBs are more concentrated towards the Galactic plane. We determine the rest-frame wavelength with unprecedented precision (?0 = 8620.86 ± 0.019 in air) and reveal a remarkable correspondence between the DIB velocities and the CO gas velocities, suggesting that the 862 nm DIB carrier is related to macro-molecules. Conclusions. We demonstrate the unique capacity of Gaia to trace the spatial structure of the Galactic ISM using the 862 nm DIB.@en

Context. The motion of stars has been used to reveal details of the complex history of the Milky Way, in constant interaction with its environment. Nevertheless, to reconstruct the Galactic history puzzle in its entirety, the chemo-physical characterisation of stars is essential. Previous Gaia data releases were supported by a smaller, heterogeneous, and spatially biased mixture of chemical data from ground-based observations. Aims. Gaia Data Release 3 opens a new era of all-sky spectral analysis of stellar populations thanks to the nearly 5.6 million stars observed by the Radial Velocity Spectrometer (RVS) and parametrised by the GSP-Spec module. In this work, we aim to demonstrate the scientific quality of Gaia s Milky Way chemical cartography through a chemo-dynamical analysis of disc and halo populations. Methods. Stellar atmospheric parameters and chemical abundances provided by Gaia DR3 spectroscopy are combined with DR3 radial velocities and EDR3 astrometry to analyse the relationships between chemistry and Milky Way structure, stellar kinematics, and orbital parameters. Results. The all-sky Gaia chemical cartography allows a powerful and precise chemo-dynamical view of the Milky Way with unprecedented spatial coverage and statistical robustness. First, it reveals the strong vertical symmetry of the Galaxy and the flared structure of the disc. Second, the observed kinematic disturbances of the disc seen as phase space correlations and kinematic or orbital substructures are associated with chemical patterns that favour stars with enhanced metallicities and lower [&alpha;/Fe] abundance ratios compared to the median values in the radial distributions. This is detected both for young objects that trace the spiral arms and older populations. Several &alpha;, iron-peak elements and at least one heavy element trace the thin and thick disc properties in the solar cylinder. Third, young disc stars show a recent chemical impoverishment in several elements. Fourth, the largest chemo-dynamical sample of open clusters analysed so far shows a steepening of the radial metallicity gradient with age, which is also observed in the young field population. Finally, the Gaia chemical data have the required coverage and precision to unveil galaxy accretion debris and heated disc stars on halo orbits through their [&alpha;/Fe] ratio, and to allow the study of the chemo-dynamical properties of globular clusters. Conclusions. Gaia DR3 chemo-dynamical diagnostics open new horizons before the era of ground-based wide-field spectroscopic surveys. They unveil a complex Milky Way that is the outcome of an eventful evolution, shaping it to the present day.@en

Uncertainties in radiation pressure modelling play a significant role in the thermospheric density and crosswind observations derived from the GRACE-FO accelerometer, especially during low solar activity. Under such conditions, the radiation pressure acceleration matches the magnitude of the aerodynamic acceleration along the track and exceeds it in the cross-track direction. The GRACE-FO mission has been operating for several years at such high altitudes during both low and rising solar activity, providing a perfect opportunity to study the effects of radiation pressure. This research uses ray tracing based on a high-fidelity satellite geometry model to calculate the radiation pressure acceleration. We numerically fine-tuned the coefficients describing the thermo-optical surface properties to obtain more accurate radiation pressure accelerations than those specified in the GRACE-FO mission manual. We also used in situ temperature measurements from thermistors on the solar arrays to model the satellite's thermal emission. These temperature measurements allowed a realistic setup of the thermal model, extended by the parameter describing the efficiency of the solar cells, and reproduced the acceleration of the thermal emission with an accuracy of RMS 0.148 nms−2. The combination of the updated thermal model and the fine-tuning of the surface coefficients improved the accuracy of the crosswind acceleration to an RMS of 0.55 nms−2, compared to an RMS of 4.22 nms−2 when using panel models and instantaneous thermal radiation. We compared the observed crosswind with two models: HWM14 and TIE-GCM. While both models capture most of the salient features of the observed crosswind, HWM14 shows particularly good agreement at high latitudes. Compared to the previously employed radiation pressure model, the crosswind observations have been improved in low and mid-latitudes, especially during periods of higher solar activity. Since the effect of radiation pressure is most significant in the crosswind direction, the effect on density was small compared to previously published datasets.@en

We present new neutral mass density and crosswind observations for the CHAMP, GRACE, and GRACE-FO missions, filling the last gaps in our database of accelerometer-derived thermosphere observations. For consistency, we processed the data over the entire lifetime of these missions, noting that the results for GRACE in 2011- 2017 and GRACE-FO are entirely new. All accelerometer data are newly calibrated. We modeled the temperature-induced bias variations for the GRACE accelerometer data to counter the detrimental effects of the accelerometer thermal control deactivation in April 2011. Further, we developed a new radiation pressure model, which uses ray tracing to account for shadowing and multiple reflections and calculates the satellitea's thermal emissions based on the illumination history. The advances in calibration and radiation pressure modeling are essential when the radiation pressure acceleration is significant compared to the aerodynamic one above 450 km altitude during low solar activity, where the GRACE and GRACE-FO satellites spent a considerable fraction of their mission lifetime. The mean of the new density observations changes only marginally, but their standard deviation shows a substantial reduction compared to thermosphere models, up to 15% for GRACE in 2009. The mean and standard deviation of the new GRACE-FO density observations are in good agreement with the GRACE observations. The GRACE and CHAMP crosswind observations agree well with the physics-based TIE-GCM winds, particularly the polar wind patterns. The mean observed crosswind is a few tens of m·s-1 larger than the model one, which we attribute primarily to the crosswind errors being positive by the definition of the retrieval algorithm. The correlation between observed and model crosswind is about 60%, except for GRACE in 2004- 2011 when the signal was too small to retrieve crosswinds reliably.@en