Abstract:We report discovery and spectroscopic follow-up of 21 astrometric binaries containing solar-type stars and dark companions with masses near 1.4 $M_{\odot}$. The simplest interpretation is that the companions are dormant neutron stars (NSs), though ultramassive white dwarfs (WDs) and tight WD+WD binaries cannot be fully excluded. We selected targets from Gaia DR3 astrometric binary solutions in which the luminous star is on the main sequence and the dynamically-implied mass of the unseen companion is (a) more than $1.25\,M_{\odot}$ and (b) too high to be any non-degenerate star or close binary. We obtained multi-epoch radial velocities (RVs) over a period of 700 days, spanning a majority of the orbits' dynamic range in RV. The RVs broadly validate the astrometric solutions and significantly tighten constraints on companion masses. Several systems have companion masses that are unambiguously above the Chandrasekhar limit, while the rest have masses between 1.25 and 1.4 $M_{\odot}$. The orbits are significantly more eccentric at fixed period than those of typical WD + MS binaries, perhaps due to natal kicks. Metal-poor stars are overrepresented in the sample: 3 out of 21 objects (14%) have [Fe/H]$\sim-1.5$ and are on halo orbits, compared to $\sim$0.5% of the parent Gaia binary sample. The metal-poor stars are all strongly enhanced in lithium. The formation history of these objects is puzzling: it is unclear both how the binaries escaped a merger or dramatic orbital shrinkage when the NS progenitors were red supergiants, and how they remained bound when the NSs formed. Gaia has now discovered 3 black holes (BHs) in astrometric binaries with masses above 9 $M_{\odot}$, and 21 NSs with masses near $1.4\,M_{\odot}$. The lack of intermediate-mass objects in this sample is striking, supporting the existence of a BH/NS mass bimodality over 4 orders of magnitude in orbital period.

The Classifications of Double Neutron Stars and their Correlations with the Binary Orbital Parameters

Investigating the Distribution of Double Neutron Stars and Unconventional Component Mass

Be X-Ray Binaries as Progenitors of Double Neutron Stars: Evidence of Similar bi-modal Classifications

Double Neutron Star Populations and Formation Channels

On the Role of Supernova Kicks in the Formation of Galactic Double Neutron Star Systems

Formation of the Double Neutron Star System Psr J1930-1852

On the Formation of Double Neutron Stars in the Milky Way: Influence of Key Parameters

Double neutron star formation via consecutive type II supernova explosions

The Mass Distribution of Galactic Double Neutron Stars: Constraints on the Gravitational-Wave Sources Like GW170817

PSR J2150+3427: A Possible Double Neutron Star System

Novel Model of an Ultra-stripped Supernova Progenitor of a Double Neutron Star

Population Properties of Neutron Stars in the Coalescing Compact Binaries

Mass Ratios of Merging Double Neutron Stars as Implied by the Milky Way Population

VLBA Astrometry of the Galactic Double Neutron Stars PSR J0509+3801 and PSR J1930–1852: A Preliminary Transverse Velocity Distribution of Double Neutron Stars and its Implications

On the Formation of PSR J1640+2224: A Neutron Star Born Massive?

VLBA Astrometry of the Galactic Double Neutron Stars PSR J0509+3801 and PSR J1930-1852: A Preliminary Transverse Velocity Distribution of Double Neutron Stars and Its Implications

Asymmetric mass ratios for bright double neutron-star mergers

Electron-capture supernovae in NS+He star systems and the double neutron star systems

The Simulation of Orbit Decay of Double Neutron Star System PSR J1906+0746 by the Gravitational Wave Radiation

A population of neutron star candidates in wide orbits from Gaia astrometry