Quantum fluctuations of W and Z bosons do not become unstable in the same parts of the model studied. A numerical map finds that transverse Z modes are stable across most of the scanned parameter space, with unstable exceptions appearing as narrow lines. Longitudinal modes occupy a much broader, band- or horn-like unstable region, while the W-boson calculation adds further structures at very low momentum.
The result is a map of possible behavior in a mathematical field model, not a measurement of particles emerging from an experiment. The preprint, posted as arXiv:2608.19855v1 on 20 August 2026, does not calculate particle numbers or decay yields.
A map of an oscillating Higgs background
The calculation follows time-dependent W- and Z-boson modes coupled to a nonlinear classical Higgs wave. It varies the initial field value and the square of the spatial momentum, then identifies which combinations produce stable oscillations and which lead to unstable solutions. The analysis treats transverse and longitudinal polarizations separately, allowing the two types of motion to be compared directly.
The background belongs to a family of dn-type nonlinear Higgs solutions. In the setup reported here, the Higgs field does not pass through zero, so the W and Z fields remain massive throughout each oscillation. That restriction defines the part of parameter space the maps are intended to describe.
To make the calculation, the authors adapt massive-vector quantization to the changing Higgs background. They derive separate equations for the two polarizations and impose the Proca constraint, which links the components of the massive vector field.
The boundary between stable and unstable motion
The researchers reduce the mode equations to Hill-type equations and use Floquet theory to classify the resulting motion. The calculation treats a Hill discriminant d greater than 2 as unstable and d less than 2 as stable. A Floquet index then characterizes the solution and its behavior in the unstable regime.
This boundary can be extremely sharp. In one selected Z transverse comparison, the same field parameter of 1.4 was paired with squared momenta of 0.38 and 0.37. The first point was labelled stable, with d = 1.998; the second was labelled unstable, with d = 2.001 and a reported Floquet index of 0.004463i.
For the Z transverse scan, the stated coupling was g_A = g_Z = 0.1374 and the Higgs self-coupling parameter was λ = 0.13. The numerical maps used grids ranging from 141 to 400 divisions along the field axis and 400 divisions along the squared-momentum axis, depending on the panel.
Z modes split sharply by polarization
The transverse Z result is the narrower of the two pictures. Most scanned combinations were stable, and the unstable exceptions formed line-shaped structures rather than a broad continuous zone. The reported transverse instability was classified as parametric amplification, meaning that the changing background can strengthen an oscillation; the analysis found no spinodal component in this sector.
The longitudinal Z map is more expansive and more varied. Its unstable region was described as band- or horn-shaped. Selected cases were classified as spinodal when the longitudinal coefficient a_L became negative, and as parametric when a_L remained positive. The model therefore contains both classifications in the longitudinal sector.
The wider longitudinal region is a feature of the modelled longitudinal mode equation, not a laboratory observation. It also does not by itself provide a particle abundance, because the preprint does not calculate cosmological particle creation through Bogoliubov coefficients.
The W map adds low-momentum features
The W-boson scan preserves the main line or band pattern seen in the Z analysis but adds extra unstable structures near zero squared momentum. With g_A = g_W = 0.1066, the additional features appear for field parameter values below 0.4 and around 1.4, both at squared momentum near 0.
A selected W example shows how the numerical classification appears in a mode solution. At field parameter 0.2 and squared momentum 0.01, the reported initial absolute mode-function value was 0.31 and the Floquet index was 0.1207. The example represents an unstable mode in the calculation, not a direct count of W particles produced.
The low-momentum structures indicate a different instability pattern for W modes within the stated model inputs. The study does not establish that they produce more W particles, more decay products or a measurable effect during an electroweak phase transition.
A map is not yet a particle-production prediction
The central gap is the conversion from unstable mode functions to particle abundances. The paper does not calculate cosmological particle creation with Bogoliubov coefficients, and it leaves the longitudinal-mode particle-creation problem specifically unaddressed. The reported Floquet growth therefore identifies candidate unstable behavior without supplying a number of particles.
The calculation also omits expansion of the Universe while searching for parameter sets that trigger instability. Its conclusions are based on analytical equations and numerical Hill/Floquet calculations for the specified time-dependent classical Higgs background, rather than on an evolving cosmological calculation.
The longitudinal scan excludes a numerically singular region near field parameter 0 when squared momentum lies between 0 and 1.0. The study also restricts its initial conditions to avoid zeros of the classical Higgs field. These choices mean that the reported map covers a defined slice of the model rather than every possible background or initial state.
Further work would need to construct a suitable longitudinal Bogoliubov framework and include cosmological expansion and damping before the instability regions could be translated into particle-creation predictions. The analysis also identifies fermionic and top-quark fluctuations, and the possibility of a new asymptotic equilibrium state, as questions beyond the present calculation.
Paper data and sources
Original title: Instability diagram of the massive gauge quantum fields around the nonlinear massive classical wave solution
Authors: Yoshio Kitadono, Tomohiro Inagaki
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text