
in the cosmological case the metastable phase is supercooled. In this first paper we will take
advantage of the physical picture that emerges from the large body of work devoted to the
cosmological case to give an order-of-magnitude estimate of the frequency and amplitude
of the GW signal in a neutron star merger. A more refined analysis will be presented
elsewhere [14].
In the cosmological case the separation of scales is provided by the fact that the
expansion rate of the Universe, H1,withHthe Hubble rate, is much longer than the
microscopic time scale given by the inverse of the local temperature, T1. As the Universe
expands and cools down it eventually enters the metastable phase. At some point bubbles
of the stable phase begin to nucleate. These bubbles then grow and collide and eventually
leave behind a superposition of long-lived sound waves propagating on the stable phase.
At sufficiently late times turbulence may also develop. Although each of these processes
contributes to the GW spectrum, it is believed that the dominant contribution comes from
collisions of sound waves with one another. As a first approximation, we will assume that
the same may be true in a NS merger. We will first determine the peak frequency of the
GWs and then we will estimate their characteristic strain.
The peak wavelength of the produced GWs is given by the characteristic wavelength
of the sound waves, which in turn is determined by the mean bubble separation, R.To
estimate the latter we recall that, once the system is in the metastable phase, the bubble
nucleation rate per unit volume takes the form
(t)
V=µ(t)4eS[µ(t)] ,(2.1)
where Sis the action of the critical bubble at the chemical potential µ
. This scale, which is the dominant one in NS, would be replaced by the temperature in
the cosmological case. The phase transition starts when the nucleation rate in the available
volume of metastable phase becomes comparable to the characteristic evolution rate of the
system. NS merger simulations show that the typical size of a HoCS is L⇠1 km. In the
cosmological case this would be the Hubble radius, H1. The evolution rate has units of
inverse time and measures the rate at which the physical properties of the system change,
for example 1
⌧=1
µ
dµ
dt .(2.2)
NS simulations show that this is of order ⌧⇠1 ms. In the cosmological case this would be
the Hubble expansion rate, H1. Using these values, together with µ⇠1 GeV, we see that
the phase transition starts when the action of the critical bubble reaches the numerical
value
S⇠
for the HoCS with a characteristic time scale set by the merger ⌧⇠1ms. Its relation
to the heating/compression rate is,
Recent NS merger simulations show that a typical size LHoCS ⇠1km for the HoCS
with a characteristic time scale set by the merger ⌧⇠1ms. Its relation to the heat-
ing/compression rate is,
–3–
A
C
Non-superconducting
Superconducting
in the cosmological case the metastable phase is supercooled. In this first paper we will take
advantage of the physical picture that emerges from the large body of work devoted to the
cosmological case to give an order-of-magnitude estimate of the frequency and amplitude
of the GW signal in a neutron star merger. A more refined analysis will be presented
elsewhere [14].
In the cosmological case the separation of scales is provided by the fact that the
expansion rate of the Universe, H1,withHthe Hubble rate, is much longer than the
microscopic time scale given by the inverse of the local temperature, T1. As the Universe
expands and cools down it eventually enters the metastable phase. At some point bubbles
of the stable phase begin to nucleate. These bubbles then grow and collide and eventually
leave behind a superposition of long-lived sound waves propagating on the stable phase.
At sufficiently late times turbulence may also develop. Although each of these processes
contributes to the GW spectrum, it is believed that the dominant contribution comes from
collisions of sound waves with one another. As a first approximation, we will assume that
the same may be true in a NS merger. We will first determine the peak frequency of the
GWs and then we will estimate their characteristic strain.
The peak wavelength of the produced GWs is given by the characteristic wavelength
of the sound waves, which in turn is determined by the mean bubble separation, R.To
estimate the latter we recall that, once the system is in the metastable phase, the bubble
nucleation rate per unit volume takes the form
(t)
V=µ(t)4eS[µ(t)] ,(2.1)
where Sis the action of the critical bubble at the chemical potential µ
. This scale, which is the dominant one in NS, would be replaced by the temperature in
the cosmological case. The phase transition starts when the nucleation rate in the available
volume of metastable phase becomes comparable to the characteristic evolution rate of the
system. NS merger simulations show that the typical size of a HoCS is L⇠1 km. In the
cosmological case this would be the Hubble radius, H1. The evolution rate has units of
inverse time and measures the rate at which the physical properties of the system change,
for example 1
⌧=1
µ
dµ
dt .(2.2)
NS simulations show that this is of order ⌧⇠1 ms. In the cosmological case this would be
the Hubble expansion rate, H1. Using these values, together with µ⇠1 GeV, we see that
the phase transition starts when the action of the critical bubble reaches the numerical
value
S⇠
for the HoCS with a characteristic time scale set by the merger ⌧⇠1ms. Its relation
to the heating/compression rate is,
Recent NS merger simulations show that a typical size LHoCS ⇠1km for the HoCS
with a characteristic time scale set by the merger ⌧⇠1ms. Its relation to the heat-
ing/compression rate is,
–3–
Figure 1. Two possible phase transitions in QCD, indicated by the solid red curves. Tand µare
the temperature and the baryon chemical potential, respectively. The dotted black curve on the
left panel shows a possible evolution of a region of a NS merger as this region is heated and/or
compressed. The points dubbed A,A0and Ccorrespond to the states shown in Fig. 3.
refer to these regions as “Hot or Compressed Spots”, or HoCS for short — see Fig. 2. The
characteristic size of these HoCS is [12]
L'5 km .(1.1)
To our knowledge, no simulation based on an EoS that includes a transition to a
color-superconducting phase has been performed. Nevertheless, existing simulations have
shown that the merger leads to the formation of cold quark-matter regions in which the
baryon density can be ten times larger than the nuclear saturation density. This makes it
conceivable that color-superconducting matter may be formed in NS mergers.
Numerical simulations have shown that the presence of HoCS affects the overall merger
dynamics, mainly because of a softening of the EoS, and leads to detectable modifications
of the GW spectrum in the usual kilo-Hertz (kHz) frequency range. One may think of
this as a “macroscopic” consequence of the phase transition. However, to the best of our
knowledge, the phase transition dynamics in HoCS has not been previously considered.
The purpose of this paper is to show that this dynamics leads to a characteristic signal
in the Mega-Hertz (MHz) range. One may think of this as a “microscopic” consequence
of the phase transition. We will see that this conclusion is extremely robust since it only
depends on generic properties of a FOPT. Therefore we do not need to commit ourselves
to a particular type of FOPT. The only assumptions we need to make are that some FOPT
is present in QCD and that this is accessed by the merger dynamics.
The key idea is that the millisecond characteristic time scale for the evolution of the
merger is much longer than the characteristic nuclear time scale, 1 fm '10−23 s. This
huge separation of scales means that, from the viewpoint of QCD processes, the merger
evolution is adiabatic to an extremely good approximation. In turn, this implies that the
HoCS are initially born as carefully prepared metastable regions of superheated and/or
– 2 –