Francesca Figueras

Distribution of known fundamental‐mode classical Cepheids in the Milky Way (data from Pietrukowicz et al. 2021). Top: positions in Galactocentric Cartesian coordinates, colour‐coded by age; bottom: Age versus Galactocentric distance, colour‐coded by metallicity (for objects with spectroscopic [Fe/H] measurements). The age distribution of Milky Way Cepheids is a particularly interesting and puzzling observable in the Galactic astrophysics context. Finding the youngest Cepheids in the inner Galaxy seems to be in conflict with the widely accepted inside‐out formation scenario where by star formation progressively reached the outer regions of the Galaxy over the past billion years. A solution of this puzzle is offered in the paper by F. Anders et al. published in this issue (e70009). image

Context. Despite the fundamental importance of the star formation history (SFH) and the initial mass function (IMF) in the description of the Milky Way, their consistent and robust derivation is still elusive. Recent and accurate astrometry and photometry collected by the Gaia satellite provide the natural framework to consolidate these ingredients in our local Galactic environment. Aims. We aim to simultaneously infer the IMF and the SFH of the Galactic disc by comparing Gaia data with the mock catalogue resulting from the Besançon population synthesis model (BGM). Our goal is also to estimate the impact of the systematics present in current stellar evolutionary models (SEMs) on this inference. Methods. We used a new implementation of the BGM fast approximate simulations (BGM FASt) framework to fit the seven-million-star Gaia DR3 all-sky G < 13 colour–magnitude diagram (CMD) to the most up-to-date dynamically self-consistent BGM. Results. Our derived SFH supports an abrupt decrease in star formation approximately 1–1.5 Gyr ago followed by a significant enhancement with a wide plateau in the range of 2–6 Gyr ago. A remarkable hiatus appears around 5–7 Gyr ago, with a ∼ 1 Gyr shift depending on the set of stellar models. A complex evolution at ages older than 8 Gyr deserves further investigation. Precise but discrepant values using different SEMs are found for the power-law indices of the IMF. In our fiducial execution with PARSEC SEM, the slope takes a value of α₂ = 1.45_-0.12^+0.19 for the range [0.5–1.53] M_⊙, while for masses larger than 1.53 M_⊙ we obtain α₃ = 1.98_-0.05^+0.13 Using STAREVOL SEM, the inferred values are α₂ = 2.48_-0.11^+0.09 and α₃ = 1.64_-0.02^+0.15. We find the solution with PARSEC to have a significantly higher likelihood than that obtained with STAREVOL. Conclusions. The current implementation of the BGM FASt framework is ready to address executions fitting all-sky Gaia data up to an apparent limiting magnitude of 14–17. This will naturally allow us to derive both a reliable SFH for the early epochs of the Galactic disc evolution and a precise slope for the IMF at low masses.

We revisit the problem of the positive correlation between age and Galactocentric distance seen in Galactic Classical Cepheids, which at first sight may seem counter‐intuitive in the context of inside‐out galaxy formation. To explain it, we use the Besançon Galaxy Model and a simulation of star particles in the Galactic disc coupled with stellar evolutionary models. We then select Classical Cepheids from this simulation and test in qualitative terms which ingredients are necessary to find agreement with the observational data. We show that the interplay of the Galactic disc's metallicity gradient and the metallicity dependence of the Cepheids' life‐time in the instability strip results in a pronounced positive age‐Galactocentric distance relation. This renders a reconstruction of the recent star‐formation history based on Classical Cepheids unrealistic. It also has important consequences on our interpretation of the observed scatter about the radial metallicity gradient measured with Galactic Classical Cepheids.

The TW Hydrae Association (TWA) is a young local association (YLA) about 50 pc from the Sun, offering a unique opportunity to study star and planet formation processes in detail. We characterized TWA’s location, kinematics, and age, investigating its origin within the Scorpius-Centaurus (Sco-Cen) OB association. Using Gaia DR3 astrometric data and precise ground-based radial velocities, we identified substructures within TWA, tentatively dividing them into TWA-a and TWA-b. Sco-Cen’s massive cluster σ Cen (15 Myr, 1805 members) may have influenced TWA’s formation. The alignment of σ Cen, TWA-a, and TWA-b in 3D positions, velocities, and ages resembles patterns in regions such as Corona Australis, suggesting that TWA is part of a cluster chain from sequential star formation induced by massive stars in Sco-Cen. TWA’s elongation in the opposite direction to that produced by Galactic differential rotation indicates its shape is still influenced by its formation processes and will dissipate in less than 50 Myr due to Galactic forces. These findings unveil the nature of YLAs and low-mass clusters in a new light. We propose that clusters such as ϵ Chamaeleontis, η Chamaeleontis, and TWA were forged by stellar feedback from massive stars in Sco-Cen, while others–such as β Pictoris, Carina, Columba, and Tucana-Horologium–are older and formed differently. Remarkably, all these YLAs and Sco-Cen are part of the α Persei cluster family, a vast kiloparsec-scale star formation event active over the past 60 Myr. This suggests that YLAs are the smallest stellar structures emerging from major star formation episodes and should be common in the Milky Way. Crucially, their formation in regions with intense stellar feedback may have influenced planet formation in these systems.

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

Context. Strongly lensed quasars are fundamental sources for cosmology. TheGaiaspace 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 previousGaiaData 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 allGaiadetections 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 theGaiaarchive. 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 ones. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the uniqueGaiadata for strong lensing studies further.

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, G_BP, or G_RP 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, G_BP, and G_RP 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.

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.

Context. We report the exploitation of a sample of Solar System observations based on data from the third Gaia Data Release (Gaia DR3) of nearly 157 000 asteroids. It extends the epoch astrometric solution over the time coverage planned for the Gaia DR4, which is not expected before the end of 2025. This data set covers more than one full orbital period for the vast majority of these asteroids. The orbital solutions are derived from the Gaia data alone over a relatively short arc compared to the observation history of many of these asteroids. Aims. The work aims to produce orbital elements for a large set of asteroids based on 66 months of accurate astrometry provided by Gaia and to assess the accuracy of these orbital solutions with a comparison to the best available orbits derived from independent observations. A second validation is performed with accurate occultation timings. Methods. We processed the raw astrometric measurements of Gaia to obtain astrometric positions of moving objects with 1D sub-mas accuracy at the bright end. For each asteroid that we matched to the data, an orbit fitting was attempted in the form of the best fit of the initial conditions at the median epoch. The force model included Newtonian and relativistic accelerations to derive the observation equations, which were solved with a linear least-squares fit. Results. Orbits are provided in the form of state vectors in the International Celestial Reference Frame for 156 764 asteroids, including near-Earth objects, main-belt asteroids, and Trojans. For the asteroids with the best observations, the (formal) relative uncertainty σ_a/a is better than 10⁻¹⁰. Results are compared to orbits available from the Jet Propulsion Laboratory and MPC. Their orbits are based on much longer data arcs, but from positions of lower quality. The relative differences in semi-major axes have a mean of 5 × 10⁻¹⁰ and a scatter of 5 × 10⁻⁹.

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.