Fluctuation-Induced Friction in Bubble-Wall Dynamics of Cosmological First-Order Phase Transitions

Author(s)

Wei, Dongdong, Guo, Zong-Kuan

Abstract

We study bubble-wall dynamics in cosmological first-order phase transitions in a two-scalar-field model, where the wall is formed by $ϕ$ and an additional real scalar $s$ couples through a portal interaction. We evolve the coupled classical field equations on the lattice and demonstrate that for an initial Bose--Einstein distribution of $s$ fluctuations at the nucleation temperature $T_n$, the resulting patchy background intermittently modulates the local driving pressure on the wall. The wall therefore undergoes alternating episodes of acceleration and deceleration and approaches a quasi-stationary propagation regime with a smaller time-averaged speed than in the decoupled limit. We further identify three familiar propagation profiles -- deflagration, detonation, and hybrid -- distinguished by where the dynamical $s$-sector energy density is concentrated relative to the wall. These effects can impact gravitational wave and baryogenesis predictions.

Figures

Potential slices $V(\phi,s,T_n)$ as functions of $\phi/T_c$ at the nucleation temperature $T_n$, evaluated for five fixed values of $s$ (see legend). Model parameters are the same as in Table~\ref{tab:benchmark_M2}, and $\lambda_{\phi s}=1$.
Caption Potential slices $V(\phi,s,T_n)$ as functions of $\phi/T_c$ at the nucleation temperature $T_n$, evaluated for five fixed values of $s$ (see legend). Model parameters are the same as in Table~\ref{tab:benchmark_M2}, and $\lambda_{\phi s}=1$.
Spacetime evolution of the field $\phi$ along the $z$ axis. We show the line-out $\phi(t,z)$ taken on the central slice ($x=y=0$), i.e., along the $z$ axis through the bubble center. Top: vanishing portal coupling, $\lambda_{\phi s}=0$. Bottom: finite coupling, $\lambda_{\phi s}=7.5$.
Caption Spacetime evolution of the field $\phi$ along the $z$ axis. We show the line-out $\phi(t,z)$ taken on the central slice ($x=y=0$), i.e., along the $z$ axis through the bubble center. Top: vanishing portal coupling, $\lambda_{\phi s}=0$. Bottom: finite coupling, $\lambda_{\phi s}=7.5$.
Spacetime evolution of the field $\phi$ along the $z$ axis. We show the line-out $\phi(t,z)$ taken on the central slice ($x=y=0$), i.e., along the $z$ axis through the bubble center. Top: vanishing portal coupling, $\lambda_{\phi s}=0$. Bottom: finite coupling, $\lambda_{\phi s}=7.5$.
Caption Spacetime evolution of the field $\phi$ along the $z$ axis. We show the line-out $\phi(t,z)$ taken on the central slice ($x=y=0$), i.e., along the $z$ axis through the bubble center. Top: vanishing portal coupling, $\lambda_{\phi s}=0$. Bottom: finite coupling, $\lambda_{\phi s}=7.5$.
Bubble-wall velocity evolution extracted from the isosurface condition $\phi_\ast\equiv(\phi_{\rm false}+\phi_{\rm true})/2$. Black solid: $\lambda_{\phi s}=0$. Blue solid: $\lambda_{\phi s}=7.5$. The horizontal red dashed lines indicate the corresponding terminal velocities estimated at $tT_c=200$.
Caption Bubble-wall velocity evolution extracted from the isosurface condition $\phi_\ast\equiv(\phi_{\rm false}+\phi_{\rm true})/2$. Black solid: $\lambda_{\phi s}=0$. Blue solid: $\lambda_{\phi s}=7.5$. The horizontal red dashed lines indicate the corresponding terminal velocities estimated at $tT_c=200$.
Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Caption Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Caption Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Caption Time evolution of a vacuum bubble in the two-scalar-field system. The simulation volume and discretization are $L T_c=500$, $dx\,T_c=0.1$, and $dt=0.2\,dx$. The parameters entering $V(\phi,s,T)$ are chosen to match the benchmark point used in the $3+1$ dimensions simulations, while the portal coupling is varied as $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ from top to bottom.
Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
Caption Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
Caption Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
Caption Profiles of the $s$ dynamical energy density (left axis, $10^{-1}\rho_s/T_c^4$) and the field $\phi/T_c$ at the fixed time $tT_c=180$, shown as functions of the coordinate $\xi\equiv x/t$. In this representation $\xi$ plays the role of a velocity coordinate (in units of $c$), so the wall position directly visualizes the propagation speed. From top to bottom the portal coupling is $\lambda_{\phi s}=0.65$, $1.1$, and $1.8$ (other potential parameters are kept fixed).
References
  • [1] P. Amaro-Seoane et al. (LISA), Laser Interferometer Space Antenna, (2017), arXiv:1702.00786 [astro-ph.IM].
  • [2] M. Colpi et al. (LISA), LISA Definition Study Report, (2024), arXiv:2402.07571 [astro-ph.CO].
  • [3] W.-R. Hu and Y.-L. Wu, The Taiji Program in Space for gravitational wave physics and the nature of gravity, Natl. Sci. Rev. 4, 685 (2017).
  • [4] Z. Luo, Y. Wang, Y. Wu, W. Hu, and G. Jin, The Taiji program: A concise overview, PTEP 2021, 05A108 (2021).
  • [5] J. Mei et al. (TianQin), The TianQin project: current progress on science and technology, PTEP 2021, 05A107 (2021), arXiv:2008.10332 [gr-qc].
  • [6] J. Luo et al. (TianQin), TianQin: a space-borne gravitational wave detector, Class. Quant. Grav. 33, 035010 (2016), arXiv:1512.02076 [astro-ph.IM].
  • [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, 041301 (2014), arXiv:1304.2433 [hep-ph].
  • [8] J. R. Espinosa, T. Konstandin, J. M. No, and G. Servant, Energy Budget of Cosmological First-order Phase Transitions, JCAP 06, 028, arXiv:1004.4187 [hep-ph].
  • [9] C. Caprini et al., Detecting gravitational waves from cosmological phase transitions with LISA: an update, JCAP 03, 024, arXiv:1910.13125 [astro-ph.CO].
  • [10] D. E. Morrissey and M. J. Ramsey-Musolf, Electroweak baryogenesis, New J. Phys. 14, 125003 (2012), arXiv:1206.2942 [hep-ph].
  • [11] J. M. Cline and K. Kainulainen, Electroweak baryogenesis at high bubble wall velocities, Phys. Rev. D 101, 063525 (2020), arXiv:2001.00568 [hep-ph].
  • [12] X. Wang, F. P. Huang, and X. Zhang, Phase transition dynamics and gravitational wave spectra of strong firstorder phase transition in supercooled universe, JCAP 05, 045, arXiv:2003.08892 [hep-ph].
  • [13] D. Bodeker and G. D. Moore, Can electroweak bubble walls run away?, JCAP 05, 009, arXiv:0903.4099 [hepph].
  • [14] L. Leitao and A. Megevand, Gravitational waves from a very strong electroweak phase transition, JCAP 05, 037, arXiv:1512.08962 [astro-ph.CO].
  • [15] G. D. Moore and T. Prokopec, Bubble wall velocity in a first order electroweak phase transition, Phys. Rev. Lett. 75, 777 (1995), arXiv:hep-ph/9503296.
  • [16] T. Konstandin, G. Nardini, and M. Quiros, Gravitational Backreaction Effects on the Holographic Phase Transition, Phys. Rev. D 82, 083513 (2010), arXiv:1007.1468 [hep-ph].
  • [17] T. Krajewski, M. Lewicki, and M. Zych, Bubble-wall velocity in local thermal equilibrium: hydrodynamical simulations vs analytical treatment, JHEP 05, 011, arXiv:2402.15408 [astro-ph.CO].
  • [18] H. Kurki-Suonio and M. Laine, On bubble growth and droplet decay in cosmological phase transitions, Phys. Rev. D 54, 7163 (1996), arXiv:hep-ph/9512202.
  • [19] J. Ignatius, K. Kajantie, H. Kurki-Suonio, and M. Laine, The growth of bubbles in cosmological phase transitions, Phys. Rev. D 49, 3854 (1994), arXiv:astro-ph/9309059.
  • [20] D. Pı̂rvu, M. C. Johnson, and S. Sibiryakov, Bubble velocities and oscillon precursors in first-order phase transitions, JHEP 11, 064, arXiv:2312.13364 [hep-th].
  • [21] M. Quiros, Finite temperature field theory and phase transitions, in ICTP Summer School in High-Energy Physics and Cosmology (1999) pp. 187–259, arXiv:hepph/9901312.
  • [22] S. Profumo, M. J. Ramsey-Musolf, C. L. Wainwright, and P. Winslow, Singlet-catalyzed electroweak phase transitions and precision Higgs boson studies, Phys. Rev. D 91, 035018 (2015), arXiv:1407.5342 [hep-ph].
  • [23] A. Tofighi, O. N. Ghodsi, and M. Saeedhoseini, Phase transition in multi-scalar-singlet extensions of the Standard Model, Phys. Lett. B 748, 208 (2015), arXiv:1510.00791 [hep-ph].
  • [24] R. Zhou, J. Yang, and L. Bian, Gravitational Waves from first-order phase transition and domain wall, JHEP 04, 071, arXiv:2001.04741 [hep-ph].
  • [25] Y.-J. Li, J. Liu, and Z.-K. Guo, Dynamical backreaction of a mass-acquiring scalar field on first-order phase transitions, Phys. Rev. D 113, 023550 (2026), arXiv:2508.14665 [astro-ph.CO].
  • [26] D. J. Weir, Gravitational waves from a first order electroweak phase transition: a brief review, Phil. Trans. Roy. Soc. Lond. A 376, 20170126 (2018), [Erratum: Phil.Trans.Roy.Soc.Lond.A 381, 20230212 (2023)], arXiv:1705.01783 [hep-ph].
  • [26] D. J. Weir, Gravitational waves from a first order electroweak phase transition: a brief review, Phil. Trans. Roy. Soc. Lond. A 376, 20170126 (2018), [Erratum: Phil.Trans.Roy.Soc.Lond.A 381, 20230212 (2023)], arXiv:1705.01783 [hep-ph].