Instability diagram of the massive gauge quantum fields around the nonlinear massive classical wave solution

Abstract: The stability and instability in the time dynamics of quantum fluctuation of the massive gauge field coupling to the nonlinear massive wave solution is studied. In particular, the instability in the transverse polarization and longitudinal polarization modes are obtained as the two dimensional plot of the initial field value parameter and the spatial momentum of the quantum massive gauge fields with the help of the theory of the Hill's equation. We present the formalism to analyze the massive gauge field by taking into account Proca constraint and we found that our formalism can predict the time dynamics of the unstable quantum mode by the Floquet index. The resulting polarization-resolved Floquet maps show that the transverse modes possess only narrow parametric-resonance bands, whereas the longitudinal modes exhibit substantially broader regions generated by both parametric and spinodal instabilities. We also find additional low-momentum instability regions for the $W$ boson that are absent or strongly suppressed in the $Z$ sector.

• The paper develops polarization-resolved Floquet instability diagrams for electroweak gauge bosons fluctuating around exact Jacobi-elliptic Higgs backgrounds, showing narrow transverse resonance bands and broader longitudinal instability regions.
• The analysis finds that longitudinal modes undergo both parametric resonance and spinodal growth because their Hill-equation coefficient can become negative, while transverse instabilities arise only through parametric resonance.
• The paper identifies additional low-momentum instability pockets for W bosons, including a representative Floquet index of 0.1207i, suggesting stronger W production than Z production in nonlinear electroweak phase-transition backgrounds.

Overview and motivation

Kitadono and Inagaki analyze the stability of quantum fluctuations in the massive electroweak gauge fields $W^\pm$ and $Z$ propagating on a classical, time-dependent Higgs background. The background is not the usual small oscillation about the vacuum but a fully nonlinear solution of the scalar equation of motion expressed through Jacobi elliptic functions, so the nonlinearity of the Higgs potential is retained exactly at the classical level (2608.19855). The central result is a set of polarization-resolved instability diagrams: two-dimensional maps over the initial field amplitude $\tilde\phi_0$ and squared dimensionless momentum $\tilde k^2$ that classify each mode as stable or unstable via Floquet theory for Hill's equations.

The work extends the authors' earlier scalar analysis to vector fields, which requires handling the Proca constraint and the subtleties of longitudinal polarization quantization with a time-dependent mass. It is also positioned relative to recent work questioning whether the static effective potential adequately captures time dynamics when particle creation is included; the instability maps here are intended as a diagnostic tool for identifying which parameter regions dominate particle production during an electroweak phase transition.

Classical background: Jacobi elliptic solutions

The classical Higgs field obeys $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$, whose general bounded solutions are parameterized by the initial value $\phi_0 = \tilde\phi_0 v$:

• $1/dn$ branch ($0 < \tilde\phi_0 < 1$): $\phi_{\rm cl} = \phi_0/dn(p_1\cdot x, k_1)$
• $dn$ branch ($Z$0): $Z$1
• $Z$2 branch ($Z$3): $Z$4

The paper restricts attention to the two $Z$5-type branches, a choice with physical consequences: since $Z$6, the Higgs field never vanishes and the gauge boson masses $Z$7 remain nonzero throughout each oscillation period. This distinguishes the setup from earlier studies using $Z$8 backgrounds, where repeated zeros of the mass occur. In the limit $Z$9 the solution reduces smoothly to the vacuum expectation value, so it describes a quasi-excited state of the classical Higgs condensate. The analysis is performed in the rest frame of the classical four-momentum, and cosmological expansion is neglected — an assumption the authors state explicitly, deferring scale-factor effects to future work.

Quantization with the Proca constraint

Expanding around the classical solution, $\tilde\phi_0$0, the quadratic Lagrangian gives the vector EOM

$\tilde\phi_0$1

with $\tilde\phi_0$2 and $\tilde\phi_0$3. Taking the divergence yields the Proca constraint $\tilde\phi_0$4, which constrains solutions rather than gauges because the effective mass is time dependent.

Transverse modes satisfy $\tilde\phi_0$5 and $\tilde\phi_0$6, so the constraint is automatically satisfied in the rest frame of the background. The mode function $\tilde\phi_0$7 obeys a Klein–Gordon-type Hill equation with coefficient $\tilde\phi_0$8, structurally identical to the scalar fluctuation equation analyzed previously by the same authors.

Longitudinal modes are more involved. Following Marañón-González and Navarro-Salas, the polarization vector itself carries time dependence through an auxiliary function $\tilde\phi_0$9, and the temporal component $\tilde k^2$0 is nondynamical — algebraically determined by the spatial mode via $\tilde k^2$1. Imposing both the Proca constraint and the Wronskian normalization fixes $\tilde k^2$2 in terms of $\tilde k^2$3, and after eliminating the first-derivative term via the field redefinition $\tilde k^2$4, one obtains a Hill equation for $\tilde k^2$5 with coefficient containing, among other terms,

$\tilde k^2$6

which arises directly from the classical Higgs EOM. This term can drive the periodic coefficient negative, opening the possibility of spinodal (tachyonic-type) instability absent in the transverse sector. The reference mass entering the longitudinal polarization vector is defined from the zero mode of the Fourier expansion of $\tilde k^2$7, giving $\tilde k^2$8, which correctly reproduces $\tilde k^2$9 as $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$0.

Instability diagrams via Floquet theory

Both mode equations are cast as Hill equations $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$1 with periodic coefficients of period $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$2, where $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$3 is the complete elliptic integral. Stability is determined by the discriminant $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$4 built from canonical fundamental solutions: $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$5 implies stability, $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$6 exponential growth with Floquet index $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$7 fixed by $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$8. Numerical integration uses SM couplings $\partial^2\phi - \mu^2\phi + \lambda\phi^3 = 0$9, $\phi_0 = \tilde\phi_0 v$0, and $\phi_0 = \tilde\phi_0 v$1, chosen to reproduce the measured $\phi_0 = \tilde\phi_0 v$2 and Higgs masses; Wronskian violations are controlled at the $\phi_0 = \tilde\phi_0 v$3–$\phi_0 = \tilde\phi_0 v$4 level.

The main findings are:

Feature	Transverse modes	Longitudinal modes
Unstable region shape	Narrow line-shaped bands	Broad band/horn-shaped domains
Instability mechanism	Parametric resonance only ($\phi_0 = \tilde\phi_0 v$5 always)	Both parametric and spinodal ($\phi_0 = \tilde\phi_0 v$6 goes negative)
$\phi_0 = \tilde\phi_0 v$7 vs $\phi_0 = \tilde\phi_0 v$8 difference	Similar structure	Additional low-$\phi_0 = \tilde\phi_0 v$9 instabilities for $1/dn$0

Explicit solutions confirm the classification: near-boundary parameters such as $1/dn$1 yield exponentially growing modes with small imaginary Floquet indices (e.g., $1/dn$2), while adjacent stable points show bounded quasi-periodic oscillations whose long-period components lengthen as the boundary is approached.

Physical interpretation

The paper's most substantive claim is that the polarization hierarchy traces back to the structure of the longitudinal coefficient function. Because $1/dn$3 is manifestly positive definite, transverse instabilities are exclusively parametric-resonance bands. By contrast, the term $1/dn$4 inherited from the Higgs EOM renders $1/dn$5 negative during part of each cycle for certain parameters, producing spinodal instability alongside parametric amplification. The authors interpret this asymmetry as reflecting the fact that gauge boson masses originate entirely from the Higgs condensate, and that the longitudinal polarization is precisely the mode tied to mass generation — the expanded longitudinal instability region makes this connection visible at the level of quantum fluctuations.

Notably, these results demonstrate that gauge-field instability does not require zeros of the gauge boson mass: even on $1/dn$6-type backgrounds where masses never vanish, broad longitudinal instability domains emerge. For the $1/dn$7 boson, additional low-momentum instability pockets appear near $1/dn$8 and $1/dn$9 at $0 < \tilde\phi_0 < 1$0, with a representative point $0 < \tilde\phi_0 < 1$1 showing a Floquet index of $0 < \tilde\phi_0 < 1$2 — substantially stronger growth than any corresponding $0 < \tilde\phi_0 < 1$3-mode example. The authors suggest this implies $0 < \tilde\phi_0 < 1$4-boson decay products would dominate over $0 < \tilde\phi_0 < 1$5 ones if such a nonlinear background were realized during the electroweak phase transition.

Limitations and open questions

Several caveats are stated plainly in the paper. Cosmological expansion is neglected throughout, justified only heuristically by redshift suppression of the nonlinear background at late times; incorporating the scale factor is deferred. Particle production is not computed: while Bogoliubov coefficients for the transverse mode follow from the scalar case, no suitable adiabatic basis or Hamiltonian density has been constructed for the longitudinal mode, whose complicated mode equation obstructs the standard formalism. Fermionic fluctuations — notably the top quark, with its large Yukawa coupling — are entirely omitted. Finally, connecting the instability maps to the conjectured new equilibrium state of Herring, Cao, and Boyanovsky requires damping effects from expansion, which remain unaddressed; the authors propose treating the scalar case first as a tractable stepping stone.

Conclusion

This paper constructs the first systematic polarization-resolved Floquet analysis of electroweak gauge boson fluctuations on non-vanishing nonlinear Higgs backgrounds. Its principal results are the narrowness of transverse parametric-resonance bands versus the broad, spinodal-inclusive longitudinal instability domains, and the appearance of coupling-dependent low-momentum instabilities unique to the $0 < \tilde\phi_0 < 1$6 sector. The diagrams provide a concrete screening tool for identifying dominant particle-production channels in scenarios involving an oscillating, nonlinear Higgs condensate, though quantitative predictions await a completed treatment of longitudinal-mode particle creation and cosmological damping.