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.

The cold dark-matter model successfully explains both the emergence and evolution of cosmic structures on large scales and, when we include a cosmological constant, the properties of the homogeneous and isotropic Universe. However, the cold dark-matter model faces persistent challenges on the scales of galaxies. Indeed, N-body simulations predict some galaxy properties that are at odds with the observations. These discrepancies are primarily related to the dark-matter distribution in the innermost regions of the halos of galaxies and to the dynamical properties of dwarf galaxies. They may have three different origins: (1) the baryonic physics affecting galaxy formation is still poorly understood and it is thus not properly included in the model; (2) the actual properties of dark matter differs from those of the conventional cold dark matter; (3) the theory of gravity departs from General Relativity. Solving these discrepancies is a rapidly evolving research field. We illustrate some of the solutions proposed within the cold dark-matter model, and solutions when including warm dark matter, self-interacting dark matter, axion-like particles, or fuzzy dark matter. We also illustrate some modifications of the theory of gravity: Modified Newtonian Dynamics (MOND), MOdified Gravity (MOG), and $f\left(R\right)$ gravity.

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. 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⁻⁹.

We aim to verify whether refracted gravity (RG) is capable of describing the dynamics of disk galaxies without resorting to the presence of dark matter. RG is a classical theory of gravity in which the standard Poisson equation is modified with the introduction of the gravitational permittivity, which is a universal monotonic function of the local mass density. We used the rotation curves and the radial profiles of the stellar velocity dispersion perpendicular to the galactic disks of 30 disk galaxies from the DiskMass Survey (DMS) to determine the gravitational permittivity. RG describes the rotation curves and the vertical velocity dispersions by requiring galaxy mass-to-light ratios that are in agreement with stellar population synthesis models, and disk thicknesses that are in agreement with observations, once observational biases are taken into account. Our results rely on setting the three free parameters of the gravitational permittivity for each individual galaxy. However, we show that the differences of these parameters from galaxy to galaxy can, in principle, be ascribed to statistical fluctuations. We adopted an approximate procedure to estimate a single set of parameters that may properly describe the kinematics of the entire sample and suggest that the gravitational permittivity is indeed a universal function. Finally, we showed that the RG models of the individual rotation curves can only partly describe the radial acceleration relation (RAR) between the observed centripetal acceleration derived from the rotation curve and the Newtonian gravitational acceleration originating from the baryonic mass distribution. Evidently, the RG models underestimate the observed accelerations by 0.1 to 0.3 dex at low Newtonian accelerations. An additional problem that ought to be considered is the strong correlation, at much more than 5σ, between the residuals of the RAR models and three radially-dependent properties of the galaxies, whereas the DMS data show a considerably less significant correlation, at more than 4σ, for only two of these quantities. These correlations might be the source of the non-null intrinsic scatter of the RG models: this non-null scatter is at odds with the observed intrinsic scatter of other galaxy samples different from DMS, which is consistent with zero. Further investigations are required to assess whether these discrepancies in the RAR originate from the DMS sample, which might not be ideal for deriving the RAR, or whether they are genuine failures of RG.