Gravitational waves from a strong electroweak phase transition

David Weir [he/him] - davidjamesweir

University of Helsinki

saoghal.net/slides/heidelberg-zoom

Heidelberg online-Teilchentee, 9 July 2020

How can we make this work?

  • View the slides at:
    saoghal.net/slides/heidelberg-zoom
    • They should advance when I advance my slide
    • You can browse, but when I switch slides, you'll be taken to my current slide
    • Also sharing my screen, but movies look worse
  • Open the participants and chat boxes:
    • Use chat to ask questions or seek the floor
    • Feel free to ask me to go slower
    • Maybe yes/no questions to gauge the audience

What happened in the early universe? when the universe was optically opaque? in dark sectors?

After LIGO comes LISA

  • Three laser arms, 2.5 M km separation
  • ESA-NASA mission, launch by 2034
  • Mission adopted 2017 arXiv:1702.00786

LISA: "Astrophysics" signals

LISA: Early Universe Cosmology

Science Investigation 7.2: Measure, or set upper limits on, the spectral shape of the cosmological stochastic GW background.

Operational Requirement 7.2: Probe a broken power-law stochastic background from the early Universe as predicted, for example, by first order phase transitions ...

Electroweak phase transition

  • Process by which the Higgs 'switched on'
  • In the Standard Model it is a crossover
  • Possible in extensions that it would be first order
    ➥ colliding bubbles then make gravitational waves
  • Can probe with experiments (colliders etc.)


Source: arXiv:1206.2942

EW PT in the SM

Work in the 1990s found this phase diagram for the SM:

At $m_H = 125 \, \mathrm{GeV}$, SM is a crossover.
arXiv:hep-ph/9605288 ; arXiv:hep-lat/9704013; arXiv:hep-ph/9809291

How? Dimensional reduction

  • At high $T$, system looks 3D at distances $\Delta x \gg 1/T$
  • Each step involves matching Green's functions in the effective and full theories to the desired order.
  • Handles the infrared problem, light fields can be studied on lattice.

Using the dimensional reduction

  • Simulate DR'ed 3D theory on lattice arXiv:hep-ph/9508379
  • Applied to SM very successfully in the 1990s:
  • NB: if 'new physics' heavy, can integrate out in DR.

Out of equilibrium physics

  1. Bubbles nucleate and grow
  2. Expand in a plasma - create reaction fronts
  3. Bubbles + fronts collide - violent process
  4. Sound waves left behind in plasma
  5. Turbulence; damping

How the wall moves

  • In EWPT: equation of motion is (schematically)
    PRD 46 2668; hep-ph/9503296; arXiv:1407.3132; ... $$ \partial^2 \phi + V_\text{eff}'(\phi,T) + \sum_{i} \frac{d m_i^2}{d \phi} \int \frac{\mathrm{d}^3 k}{(2\pi)^3 \, 2 E_i} \delta f_i(\mathbf{k},\mathbf{x}) = 0$$
    • $V_\text{eff}'(\phi)$: gradient of finite-$T$ effective potential
    • $\delta f_i(\mathbf{k},\mathbf{x})$: deviation from equilibrium phase space density of $i$th species
    • $m_i$: effective mass of $i$th species

Force interpretation

$$ \overbrace{\partial_\mu T^{\mu\nu}}^\text{Force on $\phi$} - \overbrace{\int \frac{d^3 k}{(2\pi)^3} f(\mathbf{k}) F^\nu }^\text{Force on particles}= 0 $$

This equation is the realisation of this idea:

Field-fluid system

Using a flow ansatz for the wall-plasma system:

$$ \overbrace{\partial_\mu T^{\mu\nu}}^\text{Field part} - \overbrace{\int \frac{d^3 k}{(2\pi)^3} f(\mathbf{k}) F^\nu }^\text{Fluid part}= 0 $$

i.e.:

$$ \partial_\mu T^{\mu\nu}_\phi + \partial_\mu T^{\mu\nu}_\text{fluid} = 0 $$

Can simulate as effective model of field $\phi$ + fluid $u^\mu$

astro-ph/9309059

Key parameters for GW production

  • $T_*$, temperature
    • $T_* \sim 100 \, \mathrm{GeV} \longrightarrow \mathrm{mHz}$ today
  • $\alpha_{T_*}$, vacuum energy fraction
    • $\alpha_{T_*} \ll 1$: 'weak'
    • $\alpha_{T_*} \gtrsim 1$: 'strong'
  • $v_\mathrm{w}$, bubble wall speed
  • $\beta/H_*$, 'duration'
    • $\beta$: inverse phase transition duration
    • $H_*$: Hubble rate at transition

What sources GWs?

  • Bubbles nucleate + expand, reaction fronts form...
  • Then:
    1. $h^2 \Omega_\phi$: Bubbles collide - 'envelope phase'
    2. $h^2 \Omega_\text{sw}$: Sound waves set up - 'acoustic phase'
    3. $h^2 \Omega_\text{turb}$: [MHD] turbulence - 'turbulent phase'

  • Sources add together to give observed GW power: $$ h^2 \Omega_\text{GW} + h^2 \Omega_\text{sw} + h^2 \Omega_\text{turb}$$
  • Bubbles nucleate + expand, reaction fronts form...
  • Then:
    1. $h^2 \Omega_\phi$: Bubbles collide - 'envelope phase'
    2. $h^2 \Omega_\text{sw}$: Sound waves set up - 'acoustic phase'
    3. $h^2 \Omega_\text{turb}$: [MHD] turbulence - 'turbulent phase'

  • Sources add together to give observed GW power: $$ h^2 \Omega_\text{GW} + h^2 \Omega_\text{sw} + h^2 \Omega_\text{turb}$$

Why focus on sound waves?

  • "Envelope phase" is short-lived, one-off
    • Relevant only in exotic models where $\alpha_{T_*} \gg 1$ or vacuum transitions...
    • ... but then dynamics aren't so simple arXiv:1802.05712
  • "Turbulent phase" may never occur if Hubble damping happens on a shorter timescale
    • Turbulence: bubble radius/avg. fluid velocity
    • Hubble damping: Hubble time

➤ Focus on the hydrodynamics of the transition
and immediate aftermath.

Weak transitions

$\alpha_{T_*} \ll 1$

arXiv:1704.05871arXiv:1512.06239arXiv:1504.03291arXiv:1304.2433

[with Mark Hindmarsh, Stephan Huber,
Kari Rummukainen + LISA CosWG]

Velocity profile development: detonation vs deflagration

Deflagration $v_\mathrm{w} < c_\mathrm{s}$

Small 3D simulation example

GW power spectra

[NB: curves scaled by $t$: collapse = constant emission]

$v_\mathrm{w} < c_\mathrm{s}$

PTPlot.org

Model ⟶ ($T_*$, $\alpha_{T_*}$, $v_\mathrm{w}$, $\beta$) ⟶ this plot

[Here: $Z_2$-symmetric xSM points from arXiv:1910.13125]

GW reach of SM EFT


Source: arXiv:1903.11604

Strong transitions

$\alpha_{T_*} \sim 1$

arXiv:1906.00480

[with Daniel Cutting, Mark Hindmarsh]

Fluid profile around the wall

Deflagration $v_\mathrm{w} < c_\mathrm{s}$

NB: Deflagration front ends at a shock ($\approx c_\mathrm{s}$)

[Movies: Cosmic Defects Channel]

Deflagration $v_\mathrm{w} < c_\mathrm{s}$

NB: Deflagration front ends at a shock ($\approx c_\mathrm{s}$)

[Movies: Cosmic Defects Channel]

Strong simulation slice 1

[$\alpha_{T_*} = 0.5$, $v_\mathrm{w} = 0.92$ (detonation)], velocity $\mathbf{v}$

Fluid profile around the wall

Deflagration $v_\mathrm{w} < c_\mathrm{s}$

NB: Deflagration front ends at a shock ($\approx c_\mathrm{s}$)

[Movies: Cosmic Defects Channel]

Strong simulation slice 2

[$\alpha_{T_*} = 0.5$, $v_\mathrm{w} = 0.44$ (deflag.)], velocity $\mathbf{v}$

Strong deflagrations: walls slow

At large $\alpha_{T_*}$ reheated droplets form in front of the walls

Strong simulation slice 3

[$\alpha_{T_*} = 0.34$, $v_\mathrm{w} = 0.24$ (deflag.)], vorticity $\nabla \times \mathbf{v}$

Strong deflagrations: vortical modes

GW production can be suppressed

Suppression relative to LISA CosWG ansatz

Conclusions

  • Weak transitions: good estimates of power spectrum
    ptplot.org
  • Strong transitions still an emerging topic:
    ☛ acoustic GW production suppressed
    • Causes: vortical mode production; slower walls
    • Worst for large $\alpha_{T_*}$, small $v_\mathrm{w}$
  • Turbulence still a challenge, work ongoing
  • Anyway: LISA provides a model-independent probe of first-order phase transitions around $100~\mathrm{GeV}$
    ☛ ... but new physics must be dynamical