
Probing Formation of Double Neutron Star Binaries around 1mHz with LISA
Lucy O. McNeill∗and Naoki Seto
Department of Physics, Kyoto University, Kyoto 606-8502, Japan
(Dated: October 11, 2022)
We propose a novel method to examine whether Galactic double neutron star binaries are formed
in the LISA band. In our method, we assign an effective time fraction τto each double neutron
star binary detected by LISA. This fraction is given as a function of the observed orbital period
and eccentricity and should be uniformly distributed in the absence of in-band binary formation.
Applying statistical techniques such as the Kolmogorov-Smirnov test to the actual list of τ, we can
inspect the signature of the in-band binary formation. We discuss the prospects of this method,
paying close attention to the available sample number of Galactic double neutron star binaries
around 1mHz.
I. INTRODUCTION
Double neutron star binaries (hereafter DNSBs) are
bountiful astrophysical targets. They have been detected
as radio pulsars [1, 2], and the current sample has orbital
periods Pfrom 1.9 hours [3] to six weeks [4]. This sam-
ple contains around 20 DNSBs, and searches might be in-
complete at the faint end of the luminosity function (also
limited by the beaming fraction). Shorter period (P <
1 hour) DNSBs ought to exist in the Galaxy as well.
However, due to Doppler smearing and shorter merger
timescales, their detection would be more difficult than
the longer period ones [5].
The Laser Interferometer Space Antenna (LISA) is
planned in the 2030’s [6] and is sensitive to gravita-
tional waves (GWs) around 0.1-100mHz. It will detect
all Galactic DNSBs in the frequency range f&1.5mHz
(corresponding to the orbital period P= 2/f .20 min),
unlike the longer-period radio sample. Observationally
motivated estimates [7] suggest that at least dozens of
DNSBs exist in the Galaxy at f&1.5mHz. Numerical
galactic modelling [8, 9] predicts that LISA will detect
altogether from a few to upto hundreds of DNSBs.
DNSB formation depends on the complex interplay
between many astrophysical processes. In the general
picture, the binary must survive two supernova explo-
sions. Preceding these supernovae, various mechanisms
have been theorized with respect to mass loss/exchange
in the binary [10–13] after hydrogen burning. These pro-
cesses all play a role in determining the separation at
formation, if the binary survives.
However, the related efficiencies and rates in popula-
tions are not well established (see [14] for a review). In
particular, it’s unclear whether there is a critical mini-
mum orbital period (or maximum orbital frequency) for
isolated DNSB formation.
Dynamical encounters in star clusters is an alternative
pathway for short period DNSB formation, though their
contribution to the LISA sample is estimated to be small
[15, 16]. Also, for the dynamical scenario it will be diffi-
∗email: mcneill@tap.scphys.kyoto-u.ac.jp
cult to solidly estimate the distribution function for the
orbital periods of the generated DNSBs.
In this paper, we are interested in the possibility of
DNSB formation specifically at f&1mHz. We hereafter
call this channel as “in-band” (mHz) DNSB formation, or
simply “injections”. Considering the aforementioned the-
oretical uncertainties, it will be fruitful to observationally
examine the in-band formation in a model independent
manner.
We thus develop a statistical method to examine the
in-band formation with LISA (see also [17] for forma-
tion between the LISA band and the lower frequencies
already probed by the radio sample). Recently, several
studies have proposed to statistically deal with multiple
LISA sources in the Galaxy. Among others, a large num-
ber (∼104) of white dwarf binaries (WDBs) will be a
powerful data set for various astronomical analyses (e.g.,
probing the Galactic structure [18, 19]). In this context,
one of the authors suggested to measure the fluxes of the
Galactic WDBs in frequency space [20]. He pointed out
that the measurement will enable us to follow the col-
lective evolution of the WDBs, resulting in mergers or
stable mass transfers.
One might imagine that we can get some information
about the in-band formation by similarly measuring the
DNSB flux with LISA at various frequencies. Unfortu-
nately, LISA will detect much fewer DNSBs than WDBs,
and the small number statistics will severely limit the
flux approach for DNSBs. On the other hand, unlike
WDBs, DNSBs can be well regarded as point particle
systems in the LISA band, and their long-term orbital
evolution from GW emission can be predicted quite ac-
curately [1]. Considering these pros and cons of DNSBs,
we newly introduce the effective time fraction τ, corre-
sponding to the fraction of total time that each binary
has spent in the mHz band. Without in-band forma-
tion, the fraction τshould be uniformly distributed. By
analysing the observed list of τ, we can examine potential
in-band formation through its deviation from a uniform
distribution.
The basic assumption of our study is that the mHz
DNSB population is in the “steady state” [21]. This is
a reasonable assumption, since a DNSB passes through
the mHz band in the timescale of Myr but the Galactic
arXiv:2210.04407v1 [astro-ph.HE] 10 Oct 2022