HydroGrav: Precise hydrodynamics and gravitational waves for cosmological phase transitions
Pith reviewed 2026-06-29 04:17 UTC · model grok-4.3
The pith
HydroGrav computes gravitational wave spectra from first-order phase transitions using the exact equation of state from the effective potential rather than bag or μν approximations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
HydroGrav builds self-similar hydrodynamic fluid profiles from the exact equation of state taken directly from the effective potential and inserts them into the sound shell model to obtain gravitational wave spectra. For a Z2-symmetric extension of the Standard Model the resulting profiles and spectra differ from those generated by the bag and μν equations of state both in peak amplitude and in spectral shape. Parameter-space scans identify the regions where the discrepancies are largest, and the consequent change in LISA signal-to-noise ratio after four years is estimated.
What carries the argument
Self-similar fluid profiles generated from the exact equation of state of the effective potential, which replace bag and μν approximations inside the sound shell model for gravitational wave spectrum computation.
If this is right
- Fluid profiles differ between the exact equation of state and the bag or μν approximations in the Z2 model.
- Gravitational wave spectra obtained via the sound shell model show shifts in peak amplitude and spectral shape.
- Certain regions of the model's parameter space exhibit the largest discrepancies between exact and approximate equations of state.
- LISA signal-to-noise ratios after four years change when the exact equation of state is used instead of simplified approximations.
Where Pith is reading between the lines
- The HydroGrav method can be applied to additional particle-physics models to test whether simplified equations of state remain adequate.
- If the reported differences persist across models, existing gravitational-wave forecasts for phase transitions may require systematic revision.
- Direct comparison of HydroGrav outputs against full 3D simulations would provide an independent check on the self-similar assumption.
Load-bearing premise
The self-similar hydrodynamic profiles obtained from the exact equation of state remain valid representations of the bubble-wall dynamics throughout the transition, and the sound-shell model correctly maps those profiles onto the gravitational-wave spectrum without additional model-dependent corrections.
What would settle it
A three-dimensional numerical hydrodynamic simulation of bubble expansion in the same Z2-symmetric model that produces fluid velocity and enthalpy profiles inconsistent with HydroGrav's self-similar exact-EOS solutions would falsify the claimed mapping to gravitational wave spectra.
read the original abstract
We present HydroGrav, a C++ code used to construct self-similar fluid profiles, using the exact equation of state determined directly from the effective potential, for any particle physics model capable of producing a first-order electroweak phase transition. HydroGrav also supports the bag and $\mu\nu$ (or improved bag) equations of state and includes an implementation of the sound shell model for computing the corresponding gravitational wave spectra. Using this framework, we compare the fluid profiles and gravitational wave spectra for the simplified (bag and $\mu\nu$) and exact equations of state for a $\mathbb{Z}_2$-symmetric extension of the Standard Model. Furthermore, we perform a scan across the parameter space of this model to identify regions where the simplified and exact equations of state differ in peak amplitude and spectral shape. Finally, we estimate the effect of using the exact equation of state on the signal-to-noise ratio across the parameter space, as measured by LISA after a 4-year mission.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents HydroGrav, a C++ code that constructs self-similar fluid profiles from the exact equation of state derived directly from the effective potential for any model with a first-order electroweak phase transition. The code also implements the bag and μν equations of state and the sound-shell model for gravitational-wave spectra. For a Z2-symmetric extension of the Standard Model the authors compare fluid profiles and GW spectra between the exact and simplified EOS, scan the model parameter space to locate regions of difference in peak amplitude and spectral shape, and estimate the resulting change in LISA signal-to-noise ratio after a 4-year mission.
Significance. If the central comparisons are robust, the work supplies a practical tool for quantifying the systematic error introduced by simplified EOS approximations in GW forecasts from cosmological phase transitions, directly informing the interpretation of LISA data and the design of future detectors.
major comments (1)
- [Sound-shell model implementation and results section] The central claim that exact-EOS profiles produce identifiable differences in GW peak amplitude and spectral shape rests on feeding self-similar solutions directly into the sound-shell model. The manuscript provides no quantitative validation (e.g., comparison against 3-D hydrodynamic simulations) that quantifies the modeling error of the sound-shell decomposition when the EOS is non-polytropic and temperature-dependent in a model-specific way; this validation is load-bearing for the reported differences and the LISA SNR estimates.
minor comments (2)
- The abstract states that a parameter-space scan was performed but does not report the scanned ranges, the number of points, or the fraction of parameter space where differences exceed a stated threshold; adding these numbers would improve reproducibility.
- Notation for the exact EOS (derived from the effective potential) versus the bag and μν cases should be introduced with a single equation or table early in the text to avoid repeated re-definition.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting an important point about the sound-shell model. We respond to the major comment below.
read point-by-point responses
-
Referee: The central claim that exact-EOS profiles produce identifiable differences in GW peak amplitude and spectral shape rests on feeding self-similar solutions directly into the sound-shell model. The manuscript provides no quantitative validation (e.g., comparison against 3-D hydrodynamic simulations) that quantifies the modeling error of the sound-shell decomposition when the EOS is non-polytropic and temperature-dependent in a model-specific way; this validation is load-bearing for the reported differences and the LISA SNR estimates.
Authors: We agree that the manuscript does not contain a direct quantitative validation of the sound-shell model via comparison to 3D hydrodynamic simulations for the case of non-polytropic, temperature-dependent EOS derived from a specific effective potential. The sound-shell model is a standard approximation in the literature for computing GW spectra from self-similar fluid profiles (as used in prior works with bag and μ u EOS), and the central contribution of this work is to implement and compare the exact EOS within that established framework. Performing dedicated 3D simulations to quantify the additional modeling error for model-specific EOS lies outside the scope of the present manuscript. We will revise the results section to include an expanded discussion of the assumptions underlying the sound-shell model and its applicability to exact EOS, thereby clarifying the context and limitations of the reported amplitude and shape differences as well as the LISA SNR estimates. revision: yes
Circularity Check
No significant circularity; comparisons between independently solved EOS cases
full rationale
The manuscript implements HydroGrav to solve self-similar fluid profiles directly from the exact equation of state extracted from the effective potential, and separately from the bag and μν EOS. These profiles are then fed into the sound-shell model to generate GW spectra for comparison across the Z2 model parameter space. No equation or result reduces by construction to a fitted quantity defined from the same data, no self-citation chain supplies a load-bearing uniqueness theorem, and the central comparisons rest on distinct numerical implementations of the EOS rather than renaming or self-referential prediction. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The sound shell model maps fluid velocity and enthalpy profiles to the gravitational-wave spectrum without additional model-dependent corrections.
- domain assumption Self-similar solutions accurately capture the late-time hydrodynamic evolution of expanding bubbles.
Reference graph
Works this paper leans on
-
[1]
Witten,Cosmic Separation of Phases,Phys
E. Witten,Cosmic Separation of Phases,Phys. Rev. D30(1984) 272
1984
-
[2]
M. Kamionkowski, A. Kosowsky and M.S. Turner,Gravitational radiation from first order phase transitions,Phys. Rev. D49(1994) 2837 [astro-ph/9310044]. [3]LISAcollaboration,Laser Interferometer Space Antenna,1702.00786. [4]LISAcollaboration,LISA Definition Study Report,2402.07571
Pith/arXiv arXiv 1994
-
[3]
Caprini et al.,Science with the space-based interferometer eLISA
C. Caprini et al.,Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions,JCAP04(2016) 001 [1512.06239]. [6]LISA Cosmology Working Groupcollaboration,Cosmology with the Laser Interferometer Space Antenna,Living Rev. Rel.26(2023) 5 [2204.05434]
Pith/arXiv arXiv 2016
-
[4]
A. Mazumdar and G. White,Review of cosmic phase transitions: their significance and experimental signatures,Rept. Prog. Phys.82(2019) 076901 [1811.01948]
arXiv 2019
-
[5]
C. Caprini et al.,Detecting gravitational waves from cosmological phase transitions with LISA: an update,JCAP03(2020) 024 [1910.13125]
Pith/arXiv arXiv 2020
- [6]
-
[7]
M. Hindmarsh, S.J. Huber, K. Rummukainen and D.J. Weir,Gravitational waves from the sound of a first order phase transition,Phys. Rev. Lett.112(2014) 041301 [1304.2433]
Pith/arXiv arXiv 2014
-
[8]
M. Hindmarsh, S.J. Huber, K. Rummukainen and D.J. Weir,Numerical simulations of acoustically generated gravitational waves at a first order phase transition,Phys. Rev. D92 (2015) 123009 [1504.03291]
Pith/arXiv arXiv 2015
-
[9]
M. Hindmarsh, S.J. Huber, K. Rummukainen and D.J. Weir,Shape of the acoustic gravitational wave power spectrum from a first order phase transition,Phys. Rev. D96(2017) 103520 [Erratum: Phys.Rev.D 101, 089902 (2020)] [1704.05871]
arXiv 2017
-
[10]
Steinhardt,Relativistic detonation waves and bubble growth in false vacuum decay,Physical Review D25(1982) 2074
P.J. Steinhardt,Relativistic detonation waves and bubble growth in false vacuum decay,Physical Review D25(1982) 2074
1982
-
[11]
J. Ignatius, K. Kajantie, H. Kurki-Suonio and M. Laine,Growth of bubbles in cosmological phase transitions,Phys. Rev. D49(1994) 3854 [astro-ph/9309059]
Pith/arXiv arXiv 1994
-
[12]
J.R. Espinosa, T. Konstandin, J.M. No and G. Servant,Energy Budget of Cosmological First-order Phase Transitions,JCAP06(2010) 028 [1004.4187]
Pith/arXiv arXiv 2010
-
[13]
L. Leitao and A. Megevand,Hydrodynamics of phase transition fronts and the speed of sound in the plasma,Nucl. Phys. B891(2015) 159 [1410.3875]
Pith/arXiv arXiv 2015
- [14]
- [15]
-
[16]
C. Tian, X. Wang and C. Bal´ azs,Gravitational waves from cosmological first-order phase transitions with precise hydrodynamics,Eur. Phys. J. C85(2025) 1091 [2409.14505]
arXiv 2025
- [17]
-
[18]
M. Hindmarsh,Sound shell model for acoustic gravitational wave production at a first-order phase transition in the early Universe,Phys. Rev. Lett.120(2018) 071301 [1608.04735]. – 36 –
Pith/arXiv arXiv 2018
-
[19]
M. Hindmarsh and M. Hijazi,Gravitational waves from first order cosmological phase transitions in the Sound Shell Model,JCAP12(2019) 062 [1909.10040]
arXiv 2019
-
[20]
A. Roper Pol, S. Procacci and C. Caprini,Characterization of the gravitational wave spectrum from sound waves within the sound shell model,Phys. Rev. D109(2024) 063531 [2308.12943]
arXiv 2024
-
[21]
C.L. Wainwright,CosmoTransitions: Computing Cosmological Phase Transition Temperatures and Bubble Profiles with Multiple Fields,Comput. Phys. Commun.183(2012) 2006 [ 1109.4189]
Pith/arXiv arXiv 2012
-
[22]
J.E. Camargo-Molina, B. O’Leary, W. Porod and F. Staub, Vevacious: A Tool For Finding The Global Minima Of One-Loop Effective Potentials With Many Scalars,Eur. Phys. J. C73(2013) 2588 [1307.1477]
Pith/arXiv arXiv 2013
-
[23]
A. Masoumi, K.D. Olum and J.M. Wachter,Approximating tunneling rates in multi-dimensional field spaces,JCAP10(2017) 022 [1702.00356]
arXiv 2017
-
[24]
Fowlie,A fast C++ implementation of thermal functions,Comput
A. Fowlie,A fast C++ implementation of thermal functions,Comput. Phys. Commun.228 (2018) 264 [1802.02720]
Pith/arXiv arXiv 2018
- [25]
-
[26]
Sato,SimpleBounce : a simple package for the false vacuum decay,Comput
R. Sato,SimpleBounce : a simple package for the false vacuum decay,Comput. Phys. Commun. 258(2021) 107566 [1908.10868]
arXiv 2021
- [27]
- [28]
- [29]
-
[30]
TransitionSolver
P. Athron, C. Balazs, A. Fowlie, L. Morris, C. Wang and Y. Zhang, “TransitionSolver.” https://github.com/TransitionSolver/TransitionSolver, 2026
2026
-
[31]
A. Ekstedt, P. Schicho and T.V.I. Tenkanen,DRalgo: A package for effective field theory approach for thermal phase transitions,Comput. Phys. Commun.288(2023) 108725 [2205.08815]
arXiv 2023
-
[32]
A. Ekstedt, O. Gould and J. Hirvonen,BubbleDet: a Python package to compute functional determinants for bubble nucleation,JHEP12(2023) 056 [2308.15652]
arXiv 2023
-
[33]
A. Ekstedt, O. Gould, J. Hirvonen, B. Laurent, L. Niemi, P. Schicho et al.,How fast does the WallGo? A package for computing wall velocities in first-order phase transitions,JHEP04 (2025) 101 [2411.04970]
arXiv 2025
- [34]
-
[35]
C. Caprini, A.S. Midiri, S. Procacci and A. Roper Pol,Fluid perturbations from expanding bubbles in first-order phase transitions,2604.02240
-
[36]
J. Matuszak and C. Tasillo,TransitionListener v2.0 – Robust gravitational wave predictions for cosmological phase transitions,2605.15259
-
[37]
Wang,LeWRON: Agentic Analysis of Electroweak Phase Transitions,2606.19425
I.R. Wang,LeWRON: Agentic Analysis of Electroweak Phase Transitions,2606.19425
-
[38]
M. Laine and A. Vuorinen,Basics of Thermal Field Theory, vol. 925, Springer (2016), 10.1007/978-3-319-31933-9, [1701.01554]. – 37 –
-
[39]
J. Choi and R.R. Volkas,Real Higgs singlet and the electroweak phase transition in the Standard Model,Phys. Lett. B317(1993) 385 [hep-ph/9308234]
Pith/arXiv arXiv 1993
-
[40]
S.W. Ham, Y.S. Jeong and S.K. Oh,Electroweak phase transition in an extension of the standard model with a real Higgs singlet,J. Phys. G31(2005) 857 [hep-ph/0411352]
Pith/arXiv arXiv 2005
-
[41]
Fonseca,GroupMath: A Mathematica package for group theory calculations,Comput
R.M. Fonseca,GroupMath: A Mathematica package for group theory calculations,Comput. Phys. Commun.267(2021) 108085 [2011.01764]
arXiv 2021
-
[42]
X. Wang, C. Tian and F.P. Huang,Model-dependent analysis method for energy budget of the cosmological first-order phase transition,JCAP07(2023) 006 [2301.12328]
arXiv 2023
-
[43]
Laine,Bubble growth as a detonation,Phys
M. Laine,Bubble growth as a detonation,Phys. Rev. D49(1994) 3847 [hep-ph/9309242]
Pith/arXiv arXiv 1994
-
[44]
Kosowsky, M.S
A. Kosowsky, M.S. Turner and R. Watkins,Gravitational Waves from First Order Cosmological Phase Transitions,Phys. Rev. Lett.69(1992) 2026
1992
-
[45]
A. Megevand and A.D. Sanchez,Analytic approach to the motion of cosmological phase transition fronts,Nucl. Phys. B865(2012) 217 [1206.2339]
Pith/arXiv arXiv 2012
-
[46]
A. Megevand and A.D. Sanchez,Detonations and deflagrations in cosmological phase transitions, Nucl. Phys. B820(2009) 47 [0904.1753]
Pith/arXiv arXiv 2009
-
[47]
L. Leitao and A. Megevand,Spherical and non-spherical bubbles in cosmological phase transitions,Nucl. Phys. B844(2011) 450 [1010.2134]
Pith/arXiv arXiv 2011
-
[48]
L. Leitao and A. Megevand,Gravitational waves from a very strong electroweak phase transition, JCAP05(2016) 037 [1512.08962]
Pith/arXiv arXiv 2016
-
[49]
L. Leitao and A. Megevand,Hydrodynamics of ultra-relativistic bubble walls,Nucl. Phys. B905 (2016) 45 [1510.07747]
Pith/arXiv arXiv 2016
-
[50]
T. Konstandin and J.M. No,Hydrodynamic obstruction to bubble expansion,JCAP02(2011) 008 [1011.3735]
Pith/arXiv arXiv 2011
-
[51]
Chodos, R.L
A. Chodos, R.L. Jaffe, K. Johnson, C.B. Thorn and V.F. Weisskopf,New extended model of hadrons,Physical Review D9(1974)
1974
-
[52]
T.V.I. Tenkanen and J. van de Vis,Speed of sound in cosmological phase transitions and effect on gravitational waves,JHEP08(2022) 302 [2206.01130]
arXiv 2022
- [53]
-
[54]
W.-Y. Ai, B. Laurent and J. van de Vis,Model-independent bubble wall velocities in local thermal equilibrium,JCAP07(2023) 002 [2303.10171]
arXiv 2023
-
[55]
X. Wang, C. Tian and C. Bal´ azs,Self-consistent prediction of gravitational waves from cosmological phase transitions,2409.06599
-
[56]
S.-J. Wang and Z.-Y. Yuwen,The energy budget of cosmological first-order phase transitions beyond the bag equation of state,JCAP10(2022) 047 [2206.01148]
arXiv 2022
-
[57]
Landau and E
L. Landau and E. Lifshitz,Fluid Mechanics(1987)
1987
- [58]
- [59]
-
[60]
Gyulassy, K
M. Gyulassy, K. Kajantie, H. Kurki-Suonio and L. McLerran,Deflagrations and detonations as a mechanism of hadron bubble growth in supercooled quark-gluon plasmas,Nuclear Physics B237 (1984) 477. – 38 –
1984
-
[61]
Enqvist, J
K. Enqvist, J. Ignatius, K. Kajantie and K. Rummukainen,Nucleation and bubble growth in a first-order cosmological electroweak phase transition,Phys. Rev. D45(1992) 3415
1992
-
[62]
H. Kurki-Suonio and M. Laine,Supersonic deflagrations in cosmological phase transitions,Phys. Rev. D51(1995) 5431 [hep-ph/9501216]
Pith/arXiv arXiv 1995
-
[63]
H.-K. Guo, K. Sinha, D. Vagie and G. White,Phase Transitions in an Expanding Universe: Stochastic Gravitational Waves in Standard and Non-Standard Histories,JCAP01(2021) 001 [2007.08537]
arXiv 2021
- [64]
-
[65]
C. Tian, X. Wang and C. Bal´ azs,DeepSSM: an emulator of gravitational wave spectra from sound waves during cosmological first-order phase transitions,JCAP08(2025) 060 [2501.10244]
arXiv 2025
-
[66]
A. Megevand and S. Ramirez,Bubble nucleation and growth in very strong cosmological phase transitions,Nucl. Phys. B919(2017) 74 [1611.05853]
Pith/arXiv arXiv 2017
-
[67]
D. Cutting, M. Hindmarsh and D.J. Weir,Vorticity, kinetic energy, and suppressed gravitational wave production in strong first order phase transitions,Phys. Rev. Lett.125(2020) 021302 [1906.00480]
arXiv 2020
- [68]
-
[69]
G.D. Moore and T. Prokopec,How fast can the wall move? A Study of the electroweak phase transition dynamics,Phys. Rev. D52(1995) 7182 [hep-ph/9506475]
Pith/arXiv arXiv 1995
-
[70]
S.J. Huber and M. Sopena,The bubble wall velocity in the minimal supersymmetric light stop scenario,Phys. Rev. D85(2012) 103507 [1112.1888]
Pith/arXiv arXiv 2012
-
[71]
D. Bodeker and G.D. Moore,Can electroweak bubble walls run away?,JCAP05(2009) 009 [0903.4099]
Pith/arXiv arXiv 2009
-
[72]
D. Bodeker and G.D. Moore,Electroweak Bubble Wall Speed Limit,JCAP05(2017) 025 [1703.08215]
Pith/arXiv arXiv 2017
-
[73]
K. Schmitz,New Sensitivity Curves for Gravitational-Wave Signals from Cosmological Phase Transitions,JHEP01(2021) 097 [2002.04615]
Pith/arXiv arXiv 2021
-
[74]
E. Thrane and J.D. Romano,Sensitivity curves for searches for gravitational-wave backgrounds, Phys. Rev. D88(2013) 124032 [1310.5300]
Pith/arXiv arXiv 2013
-
[75]
Cutler,Angular resolution of the LISA gravitational wave detector,Phys
C. Cutler,Angular resolution of the LISA gravitational wave detector,Phys. Rev. D57(1998) 7089 [gr-qc/9703068]
Pith/arXiv arXiv 1998
-
[76]
LISA Science Study Team,Lisa science requirements document, Tech. Rep. ESA-L3-EST-SCI-RS-001, European Space Agency, https://www.cosmos.esa.int/documents/678316/1700384/SciRD.pdf(2018). – 39 –
2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.