Hydrodynamics of Relativistic Superheated Bubbles
Abstract
Relativistic, charged, superheated bubbles may play an important role in neutron star mergers if first-order phase transitions are present in the phase diagram of Quantum Chromodynamics. We describe the properties of these bubbles in the hydrodynamic regime. We find two qualitative differences with supercooled bubbles. First, the pressure inside an expanding superheated bubble can be higher or lower than the pressure outside the bubble. Second, some fluid flows develop metastable regions behind the bubble wall. We consider the possible role of a conserved charge akin to baryon number. The fluid flow profiles are unaffected by this charge if the speed of sound is constant in each phase, but they are modified for more general equations of state. We compute the efficiency factor relevant for gravitational wave production.
CPHT-RR050.062024
1 Introduction
A variety of arguments suggest that at least two first-order phase transitions (FOPT) may be present in the phase diagram of Quantum Chromodynamics (QCD) as a function of temperature and baryon chemical potential Stephanov:2004wx ; Alford:2007xm ; Fukushima:2010bq ; Guenther:2022wcr , as sketched in Fig. 1. One is the transition from hadronic matter to quark matter. The other is the transition from a non-superconducting to a color-superconducting phase. While these transitions are well motivated, rigorously establishing their existence is a fundamental open problem in nuclear and particle physics whose solution has resisted both theoretical and experimental attempts for decades.
Gravitational waves produced in neutron star (NS) mergers could provide direct experimental access to these phase transitions Casalderrey-Solana:2022rrn . Numerical simulations of NS mergers based on equations of state (EoS) with a hadronic-quark matter phase transition include Most:2018eaw ; Most:2019onn ; Ecker:2019xrw ; Prakash:2021wpz ; Weih:2019xvw ; Tootle:2022pvd . These studies have shown that the dynamics of the merger results in the formation of regions in which the matter is sufficiently heated and/or compressed that the thermodynamically preferred phase is the quark-matter phase. In other words, the matter in these regions is pushed along the black dotted curve in Fig. 1. For brevity, we will refer to these regions simply as “superheated regions”.
Although no simulation based on an EoS with a color-superconducting phase has been performed, the large baryon densities found in existing simulations make it conceivable that color-superconducting matter may also be formed in NS mergers.
The energy density along the black dotted curve in Fig. 1 displays the multivalued form characteristic of a FOPT, as shown in Fig. 2. The superheated region is pushed from left to right along the lower branch of the phase diagram. Once the superheating is large enough, namely once the region of interest is sufficiently deep into the metastable branch, bubbles of the stable phase begin to nucleate. The point where this happens is labeled “” in Fig. 2, the nucleated phase is labelled “”, and the transition is indicated by a vertical, solid, black arrow. We will refer to these nucleated bubbles as “superheated bubbles”. The dynamics of these bubbles can produce gravitational waves that would provide direct experimental access to the QCD phase transition Casalderrey-Solana:2022rrn . The order-of-magnitude estimate in this reference shows that the frequency of these gravitational waves falls roughly in the MHz range, and that they may be potentially observable with future detectors.
The direction of the “superheated” transition is opposite to that of a cosmological FOPT. In this case the Universe is supercooled, namely pushed from right to left along the upper branch into the upper metastable region, until it transitions down to the lower stable branch, as indicated in Fig. 2 by the dashed, vertical, red arrow. We will refer to the nucleated bubbles in this case as “supercooled bubbles”. This type of transition will also occur in a NS merger if the superheated region cools down again.
The dynamics of relativistic supercooled bubbles has been extensively studied in the literature, motivated by their potential implications for cosmological FOPT — see e.g. Hindmarsh:2020hop for a review. Motivated by their possible role in NS mergers, in this paper we initiate the study of relativistic, charged, superheated bubbles, focusing on their hydrodynamic properties in the simple case of a bag-model, conformal EoS. We find two qualitative differences compared to supercooled bubbles. First, the pressure inside a superheated bubble may be higher or lower than the pressure outside the bubble. In contrast, in supercooled bubbles the inside pressure is always higher than the outside pressure. Second, some fluid flows develop regions behind the bubble wall that are metastable for any choice of parameters in the EoS, which would therefore decay at sufficiently long times. In contrast, in supercooled transitions there always exists a choice of parameters for which the flow behind the wall only explores stable states.
The physics of relativistic and non-relativistic bubbles is expected to be qualitatively similar if the two phases that determine the phase transition have comparable compressibilities in the range of pressures accessed by the bubble dynamics. For example, relativistic supercooled bubbles find a non-relativistic analog in non-relativistic chemical combustions, for this reason often referred to as “relativistic combustions”. Thus, one may expect that these similarities are also present for relativisitic superheated bubbles. However, a qualitative difference appears in a particular case of obvious phenomenological interest: the nucleation and expansion of vapor bubbles in a superheated liquid. In this case the liquid is effectively incompressible, and this leads to a bubble wall velocity that decreases as at asymptotically late times. For completeness, we review this case in Appendix A.
Part of the contents of this paper overlap with those in Ref. Barni:2024lkj , which appeared while this paper was being typewritten.
2 Self-similar profiles
At late times, when the bubbles have become macroscopic in size, the lack of any macroscopic scale implies that these bubbles are expected to become self-similar. This means that the fluid flow profile is a function of only , with the radius of the spherical bubble. Besides the bubble wall and possible shocks in the system, the very diluted gradients mean that a perfect-fluid description should be a valid approximation of the plasma. In this case, the constitutive equations for the stress tensor and the conserved charge take the form
(1) |
with the pressure, the enthalpy, the energy density, the charge density of the conserved charge and the fluid four-velocity. The thermodynamic variables are not all independent but are related to the the temperature and chemical potential of the plasma through the EoS. In this section we will not assume any specific EoS (but later in this section we will assume a constant speed of sound).
The equations of motion are simply the conservation equations of the stress-tensor and the current:
(2) |
As in the neutral case (see e.g. Espinosa:2010hh ), we can project the motion of the flow along the direction parallel and orthogonal to the velocity of the fluid. For this purpose we multiply the conservation equations by and . Using that , we obtain the following equations:
(3) | ||||
The self-similarity assumption implies that
(4) | ||||
where is the number of spatial dimensions. Substituting these expressions into the equations of motion leads to the following set of differential equations:
(5) | |||||
(6) | |||||
(7) |
These can be recast into a more convenient form:
(8a) | ||||
(8b) | ||||
(8c) |
where
(9) |
is the speed of sound and
(10) |
is the Lorentz-boosted velocity.
If each branch of the EoS is conformal, namely, if the pressure in each phase takes the form
(11) |
with an arbitrary function, then the speed of sound is constant . In this case (8) becomes a nested system that can be sequentially solved for , and , similarly to what happens in the neutral case Espinosa:2010hh . The nested structure means that the addition of charge does not change the properties of self-similar flows for conformal EoS. In fact, the last equation simply implies that .
In contrast, if the EoS is non-conformal, the speed of sound will depend on the local properties of the plasma, , making equations (8) coupled to one another. In this case the effect of the charge is already present at the level of the self-similar flow. We leave studies of non-conformal theories for future work.
At the location of the bubble wall or of the possible shocks that will form during the expansion of the bubble, the perfect fluid description will break down. However, we can integrate the conservation laws, in the radial direction, in the rest frame of these discontinuities. Doing so leads to the following matching conditions:
(12) | ||||
where () refers to the plasma in front (behind) the discontinuity, the velocities are measured in the rest frame of the discontinuity, and is the usual Lorentz factor. Rearranging the first two matching conditions we arrive at
(13) |
Eqs. (8) and (13) determine the entire fluid flow as a function of three input parameters: the wall velocity , the nucleation energy density and the nucleation charge density . In the following we will solve these equations in the case of spherical bubbles in three spatial dimensions () and a constant speed of sound.
When is constant we can solve for the velocity in (8) independently of the rest of the variables. By assuming different initial conditions one can obtain all possible flows that are self-similar, shown in Fig. 3. Contrary to the supercooled case, we seek solutions where at (expanding bubbles) the flow has negative speed (energy inflow). In order to produce Fig. 3, we chose . We expect that different values of the speed of sound will change the picture quantitatively but not qualitatively.
For , all flows tend to , doing so from low to high values of . This implies that flows that start at can continuously connect to a state, while flows starting at will inevitably have to develop a discontinuity to do so. Moreover, in the region another discontinuity may occur, for which the matching conditions reduce to . This discontinuity is a shock and its location is displayed in Fig. 3 as a dashed blue curve.
With this information we can already see the type of fluid flows that we can get. For low velocity we expect a so-called “deflagration”, namely a flow in front of the wall that smoothly connects with the asymptotic state at rest. For high velocity, we expect a “detonation”, namely a flow behind the wall that will develop a shock, connecting with the state at rest inside the bubble. Finally, we can also have hybrid solutions, which combine features from both deflagrations and detonations. Hybrid flows consist of a flow behind the wall (as in a detonation) ending in a shock, plus a flow ahead of the wall that smoothly connects with (as in a deflagration).
There is yet another kind of flows, those that are contained in the shaded area of Fig. 3. Given that they would become multivalued, these flows cannot cross the dot-dashed line, and they must follow to higher . As they cannot smoothly connect to a state with , they jump discontinously to it at the shock location. This shock is different to the one present in hybrids and detonations as the flow in motion is behind and not ahead. An entropic analysis shows that the shock front in Fig. 3 increases entropy only if the flow in motion is ahead. Hence, no flow starting in the mentioned area is realizable.
In the next subsection we explore in more detail the three type of solutions just mentioned.
3 Solutions for a conformal equation of state
In order to analyze the fluid flows in more detail, a choice of an EoS is needed. Assuming two branches such that in each of them, we choose
(14) | ||||
where , and reflect the number of degrees of freedom in each phase and “” and “” refer to the high- and the low-energy phases, respectively. The constant is sometimes referred to as the “bag constant” and it measures the energy difference between the two phases in vacuum, that is, at . For the choices and , the phase diagram is shown in Fig. 4: There is a curve of FOPT transition at . We do not choose this EoS because we expect it to provide a realistic description of QCD but because it is the simplest example in which an investigation of superheated bubbles is possible. We will report on more realistic examples elsewhere.
Having specified an EoS, we can now solve Eqs. (8) together with the matching conditions (13) in order to obtain the entire flow.
As mentioned earlier, for low, subsonic wall velocities we have deflagration-type solutions. These are characterized by the fact that the fluid behind the wall is at rest. In contrast, the fluid in front of the wall gets into motion perturbed by the movement of the wall. The fluid in motion continuously connects at large distances with the superheated fluid at rest. An example of a deflagration is shown in the top panels of Figs. 5 and 6. Contrary to what happens in the supercooled case Espinosa:2010hh , the flow velocity is negative everywhere and, notably, there is no shock present in the system. The fluid flow in front of the wall continuously connects to the metastable state at .
In the case of a high, supersonic wall velocity , the solution will be of detonation type. In this case, the wall is moving too fast for the fluid ahead to respond and so it remains at rest. As the wall passes the fluid is set in motion and, as it cannot smoothly connect with a state of , it develops a shock down the stream. The shock then allows to connect with a fluid at rest inside the bubble. An example of a detonation is depicted in the middle panels of Figs. 5 and 6. As in the case of deflagrations, the detonation flow differs qualitatively from the supercooled case. In the latter, a rarefaction wave gets formed behind the wall, continuously connecting with the state at rest inside the bubble at . Instead, in the current case the fluid dragged by the motion of the wall develops a shock that connects with the state inside the bubble.
Finally, for wall velocities close to but below the speed of sound, , a new type of flow can be constructed combining both the flows observed in detonations and deflagrations, giving rise to a hybrid solution. An example of this is displayed in the bottom panels of Figs. 5 and 6. The defining feature is that, in the rest frame of the bubble wall, the flow in front of it moves at the speed of sound . As in any deflagration, the flow ahead of the wall continuously connects with the asymptotic state while the flow behind develops a shock that takes the fluid to rest. Hybrids appear in the current case for subsonic wall speeds, in contrast with the supercooled case Espinosa:2010hh .
The deflagration shown in Fig. 6(top) is an example of a superheated flow in which the pressure inside the bubble (at ) is lower than the pressure outside the bubble (at ), whereas the detonations (middle panel) and the hybrids (bottom panel) shown in Fig. 6 are examples in which the inside pressure is higher than the outside pressure. By constructing a large family of expanding superheated bubbles we have observed that all three types of flows can have either ordering of the pressures — see Fig. 7.
The specific detonation shown in Fig. 6(middle) has been chosen to illustrate that, for some flows, the pressure in a region just behind the wall, in this case in the region , can be negative. This is necessarily lower than the pressure of the -phase, which is always positive regardless of the temperature and of the values of the parameters , and that specify the EoS. This is in stark contrast with supercooled flows, for which there is always a choice of these parameters for which the flow behind the wall only explores stable states. Instead, in the superheated case, at arbitrarily late times there is an arbitrarily large, approximately homogeneous region behind the wall that is thermodynamically disfavoured with respect to the -phase at the same temperature. This region will eventually decay by nucleating bubbles of the -phase if the bubble is allowed to expand for arbitrarily long times. Whether this happens in practice will depend on the nucleation rates in front and behind the bubble wall, on the bubble wall velocities, etc. Nevertheless, from the viewpoint of the dynamics of a single bubble, flows with negative pressure may regarded as thermodynamically unstable. In Fig. 7 these flows correspond to points above the blue curve labelled as . Note that deflagrations cannot suffer from this problem since the phase right behind the wall is directly the stable -phase.
We have also constructed detonations with . In these, the wall connects the outside state with some state below the blue line in Fig. 3, then the flow extends until the blue line and connects with the inside state through a shock. We checked that the entropy production in these solutions can be positive. These solutions could in principle compete with the deflagrations and hybrids. However, they may be unstable and decay into the hybrids if perturbed, in analogy to what happens in the supercooled case with deflagrations and hybrids for Espinosa:2010hh . We leave this stability analysis for future work.
All the flows above can be labeled by three independent parameters, that is, , and or, equivalently, , and . In the current case in which , one can solve for the enthalpy (or energy density) and the fluid velocity first, and then obtain the density for a given choice of . Therefore, for each choice of we can extract a whole family of analogous solutions parameterized by . As in the case of supercooled bubbles, it is useful to characterize this choice in terms of the so-called “strength factor” . Its general definition is111We have inverted the order in the numerator with respect to the supercooled case to ensure that it is positive. (see e.g. Hindmarsh:2020hop ; Ares:2020lbt )
(15) |
In the present case this reduces to
(16) |
so choosing is equivalent to choosing the nucleation energy density in units of . However, for general EoS both the strength factor and the energy density must be simultaneously specified in order to characterize a bubble solution.
The solutions presented in Figs. 5 and 6 depend on the values of and , but are independent on the values of the constants in (14). However, the entropy production depends on these constants in (14), and the solutions will produce positive or negative entropy depending on the values of these constants. Below we comment further on this point.
4 Bounds
We will now see that the presence of a forbidden region in Fig. 3, together with the existence of a lower bound on the energy density in each branch, imposes bounds on the possible values of the parameters that characterize each possible fluid flow.
Let us start with detonations. Given a wall velocity , consider increasing the strength or, equivalently, decreasing the nucleation energy relative to the bag constant . This decreases the energy behind the wall and increases the speed . Clearly, the strongest transition we can consider is that for which the energy density in the fluid behind the wall reaches its vacuum energy, that is, . Precisely at this point we also reach the maximum possible velocity behind the wall, . Solving the matching conditions with these values results in the bound
(17) |
Detonations are therefore confined to the region simultaneously allowed by the conditions (17) and . This region is displayed in pink in Fig. 7.
In the case of deflagrations, a similar upper bound can be obtained for the strength, although it cannot be expressed in terms of a simple analytical formula. Fig. 3 shows that, given a , the speed ahead of the wall can at most reach the boundary of the excluded region. This happens when the speed of the fluid in front of the wall moves at the speed of sound with respect to the wall, . Notice that this is precisely the value at which hybrid solutions begin to exist. The matching conditions at the wall then lead to the following bound for the state inside the bubble:
(18) |
Translating this bound to one on can only be done numerically, and it requires solving for the entire fluid flow for deflagrations. The result can be seen in Fig. 7, where deflagrations are confined to the gray shaded region.
Finally, we can obtain a bound on hybrid solutions. Once again, we see in Fig. 3 that, given a , the flow behind the wall has to fall below the excluded region. This puts the following constraint on the velocity behind the wall:
(19) |
The upper bound, which depends on the speed of sound, means that the wall must move faster than the shock, i.e. the flow must lie below the shaded region in Fig. 3. By solving for the whole flow, these two bounds translate into the two solid black lines that confine the green region in Fig. 7. The lower (upper) bound in (19) corresponds to the upper (lower) bound on in Fig. 7. Notice that the bound on deflagrations happens to precisely match the lower bound on hybrids and therefore these two types of flows do not coexist. The matching of the two bounds could be expected as deflagrations are bounded by the flows that fulfill , which is one of the defining conditions for hybrid flows. Thus deflagrations and hybrids do not compete in the same region of parameters.
We see that, as we increase the strength factor, the first solutions that cease to exist are the hybrid flows. The next ones to be forbidden are detonations. Therefore, for strong transitions only deflagrations survive, but they are confined to regions with very low wall velocities. This picture is opposite to that in the supercooled case (Espinosa:2010hh, ), where deflagrations cease to exist for strong transitions while hybrids and detonations become confined to higher and higher wall velocities.
In Fig. 7 we have added a long-dashed red curve separating the flows for which the pressure inside the bubble is higher (flows above the curve) or lower (flows below the curve) than the pressure outside the bubble. We see that, for a given wall velocity, these two cases occur for phase transitions with high or low strength factor, respectively. Motivated by these findings, we have explored the possible orderings of the pressures for a supercooled phase transition, where the standard picture is that the inside pressure is always higher than the outside pressure. We have constructed a large family of supercooled flows and we have found no counterexample to this picture.
The bounds above are model-independent, in the sense that they follow only from the allowed regions of Fig. 3 and from the existence of a minimum energy density in each phase. As we now discuss, there are two additional bounds that depend on some details of the model, meaning on the numerical coefficients , and that specify the EoS.
The first bound follows from the existence of an upper limit on given by the minimum energy density possible for nucleation, since nucleation can only happen in the stable region of the low-energy phase. This minimum will be attained at some point on the curve of first-order phase transitions, so we have
(20) |
This bound would correspond to a horizontal line in Fig. 7. We do not show it in the figure since its position would depend on all the numerical coefficients in the EoS, but we note that it may compete with the other bounds presented in this section.
The second bound follows from entropic considerations. The flows presented above are described by ideal hydrodynamics, and therefore they do not produce entropy except at the discontinuities. At each discontinuity the entropy production must be non-negative. By Stoke’s theorem this amounts to the requirement that the entropy flow is non-negative across the discontinuity, which results in the condition
(21) |
The entropy densities on either side of the discontinuity, , depend on the numerical coefficients that specify the EoS. The flows allowed by this bound for and are those below the corresponding discontinuous curves in Fig. 7. We have verified that, for , the entropy production is positive for the flows in Figs. 5 and 6. In all the cases that we have examined we have found that the most restrictive bound is the entropic bound. It would be interesting to investigate if this is true in general.
Notice that flows that produce negative entropy are not physically meaningful as expanding, superheated bubbles. However, their time reversal gives rise to physically meaningful flows with positive entropy production that can be interpreted as contracting “drops” that were studied in Rezzolla:1995kv ; Rezzolla:2013dea .
5 Efficiency factor
Ref. Casalderrey-Solana:2022rrn showed that the dynamics of superheated bubbles in NS mergers can give rise to an interesting GW signal in the MHz range. A relevant quantity to estimate the amount of GW production is the so called efficiency factor Hindmarsh:2020hop
(22) |
which is a measure of how much energy has been deposited in the motion of the fluid. We have computed the efficiency factor for several strength factors as a function of the wall velocity in the allowed regions where we can find bubble solutions. The result is summarized in Fig. 8. The empty regions in this figure arise form the disallowed regions, above the solid black curves, in Fig. 7.
The behavior of the efficiency factor for deflagrations and detonations is quite simple, since it is monotonous in both cases. For deflagrations the efficiency factor increases with the wall velocity at fixed strength, whereas for detonations it decreases. In both cases the efficiency factor increases with the strength factor at fixed wall velocity. Hybrids exhibit a richer, non-monotonic behaviour. For low , hybrids connect the deflagrations with the detonations and the efficiency factor has a maximum at some intermediate value of the wall velocity. For moderate strengths, the maximum efficiency factor increases with the strength factor. As the strength is increased, the low-velocity hybrids cease to exist and the efficiency factor becomes a monotonically decreasing function of the wall velocity. At this stage, increasing the strength implies discarding hybrid solutions and the maximum efficiency factor decreases. In summary, the maximum efficiency factor corresponds to a hybrid at some intermediate values of both the wall velocity and the strength factor.
Fig. 8 includes all possible flows allowed by hydrodynamics. However, it must be kept in mind that there are more restrictive, model-dependent bounds that will exclude some of these flows, which in turn may reduce the possible maximum value of the efficiency factor.
6 Large jump in the number of degrees of freedom
We now turn our attention to the particular case in which the jump in the number of degrees of freedom between the two phases is large, that is, . In the supercooled case, this analysis led to an estimate for the wall velocity as a function of the nucleation temperature Sanchez-Garitaonandia:2023zqz . We will obtain an analogous result in the superheated case.
A convenient way to parameterize a large jump in the number of degrees of freedom is to rescale the -phase pressure as , with . This immediately implies that all extensive quantities in the -phase scale as . In turn, this means that the strength factor diverges as
(23) |
In view of Fig. 7, this immediately tells us that the flow must be a deflagration with very small wall velocity, since this is the only type of flow that exists for arbitrarily large strength factors. We have verified that the black solid curve on the left of Fig. 7 at small is given by
(24) |
Comparing the two equations above we conclude that the velocity must scale at small as
(25) |
This can also be seen from the following argument. In the limit the critical temperature and chemical potential approach finite, -independent values . At temperatures and chemical potentials that are not parametrically separated from these the ratios of the following quantities in front and behind the wall obey
(26) |
The pressure jump works slightly differently. By definition, . This means that near the critical curve not only but also are , and therefore so is their difference. It then follows that
(27) |
From this we conclude that fluid velocity behind the wall in the rest frame of the wall is . This is related to the velocity in the lab frame, , through
(28) |
For deflagrations behind the wall and hence
(29) |
in agreement with Eq. (25).
7 Discussion
Motivated by their possible role in NS mergers, in this paper we have studied the hydrodynamics of relativistic superheated bubbles in the presence of a conserved charge akin to baryon number. Equations (8) determine the possible flow profiles. In general, these equations are coupled to one another and one must solve for the velocity and the charge density simultaneously. Here we have observed that this problem simplifies if the speed of sound is constant in each branch of the EoS. In this case one can solve sequentially for and , meaning that the velocity profile is unaffected by the presence of the charge. For simplicity, in this paper we have restricted ourselves to this case by choosing a bag-model EoS for which in both branches. One tantalizing result is that, in the case of a large jump in the number of degrees of freedom between the two phases, only deflagrations are possible and the bubble wall velocity is related to the strength of the transition through (24). Extrapolating to the hadronic-quark matter transition in QCD this would predict a value for the wall velocity around , in agreement with the ballpark value suggested in Casalderrey-Solana:2022rrn . It would be interesting to verify these estimates with a more realistic approximation to the QCD EoS. In this case the velocity and the charge profiles will not decouple from one another and the fluid flows could be qualitatively different.
We have found two qualitative differences between superheated and supercooled flows. The first one is that some superheated detonations and hybrids can develop a necessarily metastable region behind the bubble wall, which at sufficiently long times will decay by nucleating bubbles of the -phase. The second difference is that the pressure inside an expanding superheated bubble can be higher or lower than the pressure outside it. Even in the latter case, the bubble still expands by accreting energy from the outside. This raises an interesting dynamical question about how this steady state is reached from the slightly over-critical bubble that is initially nucleated. The reason is that, for the critical bubble, the inside pressure must be slightly higher than the outside pressure, since it must exactly compensate for the effect of the latter plus the surface tension. It would be interesting to see the time evolution of the inside pressure along the first instants after nucleation in a microscopic model, following the analysis in the supercooled case Bea:2021zsu . An additional insight provided by such a model would be a prediction of the bubble wall velocity from the microscopic dynamics. We will report on these aspects elsewhere.
Acknowledgements
We thank Alexandre Serantes for useful discussions. This work was supported by the “Center of Excellence Maria de Maeztu 2020-2023” award to the ICCUB (CEX2019-000918-M) funded by MCIN/AEI/10.13039/501100011033. YB, DM and JCS acknowledge support from grants PID2019- 105614GB-C21, PID2019-105614GB-C22, PID2022-136224NB-C21, PID2022-136224NB-C22, and 2021-SGR-872. The work of YB is also funded by a “Beatriu de Pinós” postdoctoral program under the Ministry of Research and Universities of the Government of Catalonia (2022 BP 00225) and by a “Maria Zambrano” postdoctoral fellowship from the University of Barcelona. The work of MSG is supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No758759).
Appendix A Non-relativistic, superheated bubbles with an incompressible phase
Superheated bubbles occur in many everyday situations, such as the formation of vapor bubbles in superheated water. In these cases, the matter surrounding the bubble walls is non-relativistic. Additionally, for the relevant range of pressures, the matter outside the superheated bubbles behaves as an incompressible fluid. These two features make the dynamics of those superheated bubbles very different from their relativistic analogues, as described in this work. Specifically, in the non-relativistic case, the bubble size does not grow at a constant rate; instead, the growth rate decreases with time . For completeness, this appendix reviews the dynamics of non-relativistic superheated bubbles, and the discussion here follows closely the review vbblreview
We assume a non-relativistic system in its low-temperature phase that has been superheated past its critical temperature. In this system, characterized by a mass density , the low-temperature phase is incompressible, meaning remains constant. When a bubble of the thermodynamically stable high-temperature phase nucleates, it grows due to entropic preference. This phase has its own mass density, . However, similar to vapor bubbles in superheated water, we do not assume the fluid inside the bubble is incompressible. As with the relativistic case, once the bubble grows sufficiently large, the dynamics inside and outside the bubble can be described by non-relativistic hydrodynamics. The interface between the two phases can then be approximated as a discontinuity in the hydrodynamic fields.
Paralleling the discussion of the relativistic case, we begin by characterizing the fluid velocity in the vicinity of the bubble. This is analogous to the velocity profiles presented in fig. (3). In the low-temperature phase, the mass continuity equation combined with the incompressibility condition implies that the velocity field is divergenceless, i.e., . Under spherical symmetry, the velocity field outside the bubble is radial and can be expressed as:
(30) |
where is the bubble radius and is a function determined by mass conservation, analogous to the charge conservation equation in non-relativistic systems. The velocity field vanishes far from the bubble, ensuring that the fluid at a distance remains at rest.
Since the bubble wall does not accumulate mass, the mass flux on both sides of the wall must be equal. This is equivalent to the last matching condition in Eq. (12).Imposing that the matter inside the bubble is at rest, and considering that as the bubble grows, the change in density between the high and low-temperature phases must be compensated by the influx of mass from the low-temperature phase outside the wall:
(31) |
This expression for the velocity field, combined with the Navier-Stokes (NS) equation for incompressible fluids,
(32) |
where is the kinematic viscosity, determines the form of the pressure field in front of the bubble. By substituting Eqs. (30) and (31) into the NS equation and integrating in the radial direction, we obtain the pressure difference in the low-temperature phase both far away from the bubble, and next to the bubble wall :
(33) |
with . The pressure outside the bubble can be related to the pressure inside the bubble through the non-relativistic analog of the matching conditions for the momentum flux across the wall, given by the first equation in Eq. (12). Taking into account the effect of the bubble’s surface tension, , the fluid stress on both sides of the wall is given by:
(34) |
where the last term comes from the sheer stress contribution of the velocity field, Eq. (31). For very large bubbles, the effect of the surface tension and the viscosity becomes negligible; this is the reason why we have not considered those contribution in the relativistic analysis. Nevertheless, keeping those terms, and together with Eq. (33), we obtain the so called Rayleigh-Plesset (RP) equation
(35) |
which governs the dynamics of non-relativistic bubbles in incompressible fluids.
To complete the analysis, we also must consider the effect of the energy flux across the bubble wall, which in the relativistic case is given by the middle equation in Eq. (12). In general, in the non-relativistic limit, the energy flux is given by
(36) |
with and are the enthalpy and energy per unit mass, is the sheer tensor and the heat conductivity. For spherical bubbles and using Eqs. (30)and (31), the radial component of the flux right outside of the bubble wall is given by
(37) |
This expression shows that non-relativistic bubbles cannot grow at a constant rate, as their relativistic analogues do, at least for normal fluids, such that the density of the fluid phase is larger than the density of the vapour phase, . For a constant and positive , the first to terms of the expression above dominate and the energy flux is positive, which implies that the fluid motion outside the wall takes away energy from the interior of the bubble. This can also inferred from the combination of Eqs. (30)and (31), which show that the velocity of the fluid outside of the bubble is positive, unlike the relativistic case. Therefore, the only way that the latent heat can be compensated is by reducing the internal energy of the fluid outside; when this energy arrives to the critical energy, the pressure difference in the inside and outside of the bubble become equal and the fluid velocity must vanish, as inferred from Eq. (35). Therefore, superheated bubbles cannot expand at constant velocities for arbitrarily long times.
For bubbles to continue growing at asymptotic late times, energy must be transferred toward the interior of the bubble. This is achieved by heat transfer, represented by the last term in Eq. (37), which becomes the dominant mechanism at these late times. In this regime, integrating over the bubble’s surface, the energy balance between the heat flux and the latent heat can be expressed as follows:
(38) |
with the latent heat of the transition. In this limit, the temperature gradient outside the bubble is primarily governed by heat diffusion, yielding with , where is the temperature difference between the wall and the ambient, and is the thermal diffusivity of the fluid outside the wall. Consequently, in this regime, the bubble radius scales as . For a comprehensive derivation of this expression, refer to vbblreview and the references therein.
References
- (1) M.A. Stephanov, QCD Phase Diagram and the Critical Point, Prog. Theor. Phys. Suppl. 153 (2004) 139 [hep-ph/0402115].
- (2) M.G. Alford, A. Schmitt, K. Rajagopal and T. Schäfer, Color superconductivity in dense quark matter, Rev. Mod. Phys. 80 (2008) 1455 [0709.4635].
- (3) K. Fukushima and T. Hatsuda, The phase diagram of dense QCD, Rept. Prog. Phys. 74 (2011) 014001 [1005.4814].
- (4) J.N. Guenther, An overview of the QCD phase diagram at finite and , in 38th International Symposium on Lattice Field Theory, 1, 2022 [2201.02072].
- (5) J. Casalderrey-Solana, D. Mateos and M. Sanchez-Garitaonandia, Mega-Hertz Gravitational Waves from Neutron Star Mergers, 2210.03171.
- (6) E.R. Most, L.J. Papenfort, V. Dexheimer, M. Hanauske, S. Schramm, H. Stöcker and L. Rezzolla, Signatures of quark-hadron phase transitions in general-relativistic neutron-star mergers, Phys. Rev. Lett. 122 (2019) 061101 [1807.03684].
- (7) E.R. Most, L. Jens Papenfort, V. Dexheimer, M. Hanauske, H. Stoecker and L. Rezzolla, On the deconfinement phase transition in neutron-star mergers, Eur. Phys. J. A 56 (2020) 59 [1910.13893].
- (8) C. Ecker, M. Järvinen, G. Nijs and W. van der Schee, Gravitational waves from holographic neutron star mergers, Phys. Rev. D 101 (2020) 103006 [1908.03213].
- (9) A. Prakash, D. Radice, D. Logoteta, A. Perego, V. Nedora, I. Bombaci, R. Kashyap, S. Bernuzzi and A. Endrizzi, Signatures of deconfined quark phases in binary neutron star mergers, Phys. Rev. D 104 (2021) 083029 [2106.07885].
- (10) L.R. Weih, M. Hanauske and L. Rezzolla, Postmerger Gravitational-Wave Signatures of Phase Transitions in Binary Mergers, Phys. Rev. Lett. 124 (2020) 171103 [1912.09340].
- (11) S. Tootle, C. Ecker, K. Topolski, T. Demircik, M. Järvinen and L. Rezzolla, Quark formation and phenomenology in binary neutron-star mergers using V-QCD, 2205.05691.
- (12) M.B. Hindmarsh, M. Lüben, J. Lumma and M. Pauly, Phase transitions in the early universe, SciPost Phys. Lect. Notes 24 (2021) 1 [2008.09136].
- (13) G. Barni, S. Blasi and M. Vanvlasselaer, The hydrodynamics of inverse phase transitions, 2406.01596.
- (14) J.R. Espinosa, T. Konstandin, J.M. No and G. Servant, Energy Budget of Cosmological First-order Phase Transitions, JCAP 06 (2010) 028 [1004.4187].
- (15) F.R. Ares, M. Hindmarsh, C. Hoyos and N. Jokela, Gravitational waves from a holographic phase transition, JHEP 21 (2020) 100 [2011.12878].
- (16) L. Rezzolla, J.C. Miller and O. Pantano, Evaporation of quark drops during the cosmological quark - hadron transition, Phys. Rev. D 52 (1995) 3202 [astro-ph/9502064].
- (17) L. Rezzolla and O. Zanotti, Relativistic Hydrodynamics, Oxford University Press (9, 2013), 10.1093/acprof:oso/9780198528906.001.0001.
- (18) M. Sanchez-Garitaonandia and J. van de Vis, Prediction of the bubble wall velocity for a large jump in degrees of freedom, 2312.09964.
- (19) Y. Bea, J. Casalderrey-Solana, T. Giannakopoulos, D. Mateos, M. Sanchez-Garitaonandia and M. Zilhão, Bubble wall velocity from holography, Phys. Rev. D 104 (2021) L121903 [2104.05708].
- (20) A. Prosperetti, Vapor bubbles, Annual Review of Fluid Mechanics 49 (2017) 221.