We provide a purely dynamical global map of the non-axisymmetric structure of the Milky Way disk. For this, we exploited the information contained within the in-plane motions of disk stars from Gaia DR3 to adjust a model of the Galactic potential, including a detailed parametric form for the bar and spiral arms. We explored the parameter space of the non-axisymmetric components with the backward integration method, first adjusting the bar model to selected peaks of the stellar velocity distribution in the solar neighborhood, and then adjusting the amplitude, phase, pitch angle, and pattern speed of spiral arms to the median radial velocity as a function of position within the disk. We checked a posteriori that our solution also qualitatively reproduces various other features of the global non-axisymmetric phase-space distribution, including most moving groups and phase-space ridges, despite those not being primarily used in the adjustment. This fiducial model has a bar with a pattern speed of 37 km s −1 kpc −1 and two spiral modes that are twoarmed and three-armed, respectively. The two-armed spiral mode has a ~25% local contrast surface density and a low pattern speed of 13.1 km s −1 kpc −1 , and matches the location of the Crux-Scutum, Local, and Outer arm segments. The three-armed spiral mode has a ~9% local contrast density, a slightly higher pattern speed of 16.4 km s −1 kpc −1 , and matches the location of the Carina-Sagittarius and Perseus arm segments. The Galactic bar, with a higher pattern speed than both spiral modes, has recently disconnected from those two arms. The fiducial non-axisymmetric potential presented in this paper, reproducing most non-axisymmetric signatures detected in the stellar kinematics of the Milky Way disk, can henceforth be used to confidently integrate orbits within the Galactic plane.

The most metal-poor stars (e.g. [Fe/H] ≤-2.5) are the ancient fossils from the early assembly epoch of our Galaxy, very likely before the formation of the thick disc. Recent studies have shown that a non-negligible fraction of them have prograde planar orbits, which makes their origin a puzzle. It has been suggested that a later-formed rotating bar could have driven these old stars from the inner Galaxy outward, and transformed their orbits to be more rotation-dominated. However, it is not clear if this mechanism can explain these stars as observed in the solar neighborhood. In this paper, we explore the possibility of this scenario by tracing these stars backwards in an axisymmetric Milky Way potential with a bar perturber. We integrate their orbits backward for 6 Gyr under two bar models: one with a constant pattern speed and another one with a decelerating speed. Our experiments show that, under the constantly-rotating bar model, the stars of interest are little affected by the bar and cannot have been shepherded from a spheroidal inner Milky Way to their current orbits. In the extreme case of a rapidly decelerating bar, some of the very metal-poor stars on planar and prograde orbits can be brought from the inner Milky Way, but ∼90% of them were nevertheless already rotation-dominated (J_ϕ ≥ 1000 km s^-1 kpc) 6 Gyr ago. The chance of these stars having started with spheroid-like orbits with small rotation (J_ϕ ≲ 600 km s^-1 kpc) is very low (< 3%). We therefore conclude that, within the solar neighborhood, the bar is unlikely to have shepherded a significant fraction of inner Galaxy spheroid stars to produce the overdensity of stars on prograde, planar orbits that is observed today.

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further.

Gravitational waves from black-hole 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. 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. 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. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot 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. 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.

The most metal-poor tail of the Milky Way ([Fe/H] $\leq$ $-$2.5) contains a population of stars with very prograde planar orbits, which is puzzling in both their origin and evolution. A possible scenario is that they are shepherded by the bar from the inner Galaxy, where many of the old and low-metallicity stars in the Galaxy are located. To investigate this scenario, we use test-particle simulations with an axisymmetric background potential plus a central bar model. The test particles are generated by an extended distribution function (EDF) model based on the observational constraints of bulge stars. According to the simulation results, a bar with constant pattern speed cannot help bring stars from the bulge to the solar vicinity. In contrast, when the model includes a rapidly decelerating bar, some bulge stars can gain rotation and move outwards as they are trapped in the co-rotation regions of the bar. The resulting distribution of shepherded stars heavily depends on the present-day azimuthal angle between the bar and the Sun. The majority of the low-metallicity bulge stars driven outwards are distributed in the fourth quadrant of the Galaxy with respect to the Sun, and about 10$\%$ of them are within 6 kpc from us. Our experiments indicate that the decelerating bar perturbation can be a contributing process to explain part of the most metal-poor stars with prograde planar orbits seen in the Solar neighborhood but is unlikely to be the dominant one.

ABSTRACT WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, saw first light in late 2022. WEAVE comprises a new 2-deg field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable ‘mini’ integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366–959 nm at R ∼ 5000, or two shorter ranges at $R\sim 20\, 000$. After summarizing the design and implementation of WEAVE and its data systems, we present the organization, science drivers, and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy’s origins by completing Gaia’s phase-space information, providing metallicities to its limiting magnitude for ∼3 million stars and detailed abundances for ∼1.5 million brighter field and open-cluster stars; (ii) survey ∼0.4 million Galactic-plane OBA stars, young stellar objects, and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey ∼400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionized gas in z < 0.5 cluster galaxies; (vi) survey stellar populations and kinematics in ${\sim} 25\, 000$ field galaxies at 0.3 ≲ z ≲ 0.7; (vii) study the cosmic evolution of accretion and star formation using >1 million spectra of LOFAR-selected radio sources; and (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.

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 similar to 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 similar to 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.

We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 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. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ 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 GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 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 BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged)

Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. 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 by using the Gaia catalogue itself and by comparison with external data. We have 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 disk stars (~3M), FGKM spectral type stars (~3M), and UCDs (~20K). 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 furthermore identify 15740 bone fide carbon stars, 5863 solar-analogues, and provide the first homogeneous set of stellar parameters of the Spectro Photometric Standard Stars. We use a subset of the OBA sample to illustrate its usefulness to analyse 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.

Disc stars from the Gaia DR3 RVS catalogue are selected to explore the phase spiral as a function of position in the Galaxy. The data reveal a two-armed phase spiral pattern in the local $z-v_z$ plane inside the solar radius, which appears clearly when colour-coded by $\langle v_R \rangle (z,v_z)$: this is characteristic of a breathing mode that can in principle be produced by in-plane non-axisymmetric perturbations. We note the phase spiral pattern becomes single armed outside the solar radius. When a realistic analytic model with an axisymmetric background potential plus a steadily rotating bar and 2-armed spiral arms as perturbation is used to perform particle test integrations, the pseudo stars get a prominent spiral pattern in the $\langle v_R \rangle$ map in the $x-y$ plane. Additionally, clear breathing mode evidence at a few $\rm{kms}^{-1}$ level can be seen in the $\langle v_z \rangle$ map on the $x-z$ plane, confirming that such breathing modes are non-negligible in the joint presence of a bar and spiral arms. However, no phase-spiral is perceptible in the $(z, v_z)$ plane. When an initial vertical perturbation is added to all pseudo stars to carry out the simulation, the one-armed phase spirals can clearly be seen 500~Myr after the perturbation and gradually disappear inside-out. Finally, we show as a proof of concept how a toy model of a time-varying non-axisymmetric in-plane perturbation with varying amplitude and pattern speed can produce a strong two-armed phase-spiral. We conclude a time-varying strong internal perturbation together with an external vertical perturbation could perhaps explain the transition between the two-armed and one-armed phase-spirals around the Solar radius.

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 >= 1.3 M-circle dot) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), delta Sct, and gamma 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 delta 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 delta 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 HertzsprungRussell (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 ffi Sct stars and provide evidence that the relation features di fferent regimes depending on the oscillation period. We demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of ffi 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.

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 ($\omega$ 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. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $\omega$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $\omega$ 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 {\it 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. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions.

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 (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs 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 (abridged)

The Gaia DR3 Catalogue contains for the first time about eight hundred thousand solutions with either orbital elements or trend parameters for astrometric, spectroscopic and eclipsing binaries, and combinations of them. This paper aims to illustrate the huge potential of this large non-single star catalogue. Using the orbital solutions together with models of the binaries, a catalogue of tens of thousands of stellar masses, or lower limits, partly together with consistent flux ratios, has been built. 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. Illustrative applications are proposed for binaries across the H-R diagram. The binarity is studied in the RGB/AGB and a search for genuine SB1 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 presented. 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. Beside binarity, higher order multiple systems are also found.

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.

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, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provided in Gaia DR3, we select various stellar populations to explore and identify non-axisymmetric features in the disc of the Milky Way in both configuration and velocity space. Using more about 580 thousand sources identified as hot OB stars, together with 988 known open clusters younger than 100 million years, we map the spiral structure associated with star formation 4-5 kpc from the Sun. We select over 2800 Classical Cepheids younger than 200 million years, which show spiral features extending as far as 10 kpc from the Sun in the outer disc. We also identify more than 8.7 million sources on the red giant branch (RGB), of which 5.7 million have line-of-sight velocities, allowing the velocity field of the Milky Way to be mapped as far as 8 kpc from the Sun, including the inner disc. 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. Meanwhile, 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. (abridged)

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.

Gaia Data Release 3 provides novel flux-calibrated low-resolution spectrophotometry for about 220 million sources in the wavelength range 330nm - 1050nm (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 ~220 million sources in several wide bands of widely used photometric systems is provided (Gaia Synthetic Photometry Catalogue; GSPC) as well as a catalogue of $\simeq 10^5$ white dwarfs with DA/non-DA classification obtained with a Random Forest algorithm (Gaia Synthetic Photometry Catalogue for White Dwarfs; GSPC-WD).

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.

Context. Di ffuse 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 angstrom kpc(-1) 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 (lambda(0) = 8620:86 +/- 0.019 angstrom 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.