Internal symmetry to the rescue: well-posed 1+1 evolution of self-interacting vector fields

Previous studies have identified potential instabilities in self-interacting vector theories associated with the breakdown of the well-posedness of the initial-value problem. However, these conclusions are restricted to Abelian vector fields, leaving room to explore alternative setups, such as non-Abelian vector fields with internal symmetries. Building on this idea, we study the well-posed 1+1 evolution of self-interacting SU(2) vector fields minimally coupled to gravity within the framework of the ’t Hooft-Polyakov magnetic monopole configuration. In this context, we present a counterexample in which self-interacting vector fields retain a well-posed initial value problem formulation. Remarkably, this system exhibits the same characteristic speeds as those found in general relativity (GR) in one spatial dimension. Unlike its Abelian counterpart, we achieve stable numerical evolutions across a wide range of initial conditions within a fully non-linear dynamical background, as evidenced in our time integration algorithm. Although our conclusions are strictly valid for the spherical symmetry case with only magnetic part for the vector field, this study serves as a valuable diagnostic tool for investigating more realistic astrophysical scenarios in three-dimensional settings and under more general background and vector field configurations.

💡 Research Summary

**
The paper addresses a long‑standing concern that self‑interacting massive vector (Proca) theories suffer from a loss of well‑posedness of the initial‑value problem (IVP), a problem that has been demonstrated primarily for Abelian (U(1)) models. The authors ask whether this pathology persists when the vector field carries a non‑Abelian internal symmetry, specifically a global SU(2) symmetry, and whether the additional self‑interactions inherent to non‑Abelian gauge fields can cure the ill‑posedness.

To answer this, they construct a theory described by the action (1), which contains the Einstein–Hilbert term, the kinetic term for an SU(2) vector field (B^a_\mu), a standard mass term (\mu^2 B^a_\mu B^{a\mu}), and two independent quartic self‑interaction terms with couplings (\chi_1) and (\chi_2). By defining an effective self‑interaction parameter (\chi = 2\chi_1 + \chi_2), the model is reduced to two free parameters: the mass (\mu) and the combined self‑coupling (\chi). When (\chi\to0) the theory reduces to the usual non‑Abelian Proca (or massive Yang‑Mills) model.

The analysis is performed under spherical symmetry, using the well‑known ’t Hooft‑Polyakov magnetic monopole ansatz (the “magnetic” sector of the Witten ansatz). This choice sets the electric components of the gauge field to zero and leaves a single dynamical scalar function (w(t,r)) that multiplies the Pauli matrices. Crucially, this configuration automatically satisfies a generalized Lorenz condition (3) without the need for an external constraint, a feature that distinguishes it from the Abelian case where the Lorenz condition must be enforced numerically.

The spacetime metric is written in polar‑areal coordinates (ds^2 = -e^{2A(t,r)}dt^2 + e^{2B(t,r)}dr^2 + r^2 d\Omega^2). Substituting the monopole ansatz into the Einstein and field equations yields a coupled system of nonlinear PDEs: two Einstein constraints (10)–(11) for the metric functions and a wave‑type evolution equation (12) for (w). By introducing first‑order variables (Q = \partial_r w) and (P = e^{B-A}\partial_t w), the authors recast the system into a first‑order form suitable for hyperbolicity analysis.

The principal symbol of the resulting system is computed, and its eigenvalues are found to be (\pm 1) (the speed of light) and 0, exactly as in pure general relativity in 1+1 dimensions. The eigenvectors form a complete basis, confirming that the system is strongly hyperbolic. Consequently, the IVP is well‑posed: solutions exist, are unique, and depend continuously on the initial data.

Numerical experiments are carried out using a fourth‑order Runge‑Kutta time integrator and second‑order finite‑difference spatial discretisation. The authors test a wide range of initial data, including static configurations, Gaussian pulses, and sharp step‑like profiles. In all cases the evolution remains stable for long integration times, the constraints (including the generalized Lorenz condition) stay satisfied to machine precision, and no ghost or tachyonic modes appear. This contrasts sharply with Abelian self‑interacting Proca simulations, where the same numerical scheme quickly develops instabilities and the simulation crashes.

The paper explicitly acknowledges its limitations. The analysis is restricted to spherical symmetry and to the purely magnetic sector; electric components, rotating configurations, or non‑spherical perturbations are not considered. Moreover, the SU(2) symmetry is treated as a global (non‑gauge) symmetry, so the results may not directly carry over to a fully gauge‑invariant Yang‑Mills‑Proca theory. Finally, the work is confined to 1+1 dimensions; extending the hyperbolicity proof and the numerical stability to full 3+1 dimensions remains an open challenge.

Nevertheless, the study provides the first concrete counter‑example to the claim that all self‑interacting massive vector theories are ill‑posed. It demonstrates that internal non‑Abelian symmetries can preserve the hyperbolic structure of the field equations, yielding characteristic speeds identical to those of GR and enabling robust numerical evolutions. The authors suggest several avenues for future research: (i) inclusion of electric components and more general ansätze, (ii) exploration of rotating or axisymmetric backgrounds, (iii) implementation of a fully gauge‑invariant SU(2) Proca–Yang‑Mills system, and (iv) three‑dimensional simulations aimed at astrophysical applications such as black‑hole hair, neutron‑star interiors, or exotic compact objects. This work thus opens a promising path toward viable, well‑posed vector‑field models in both cosmology and high‑energy astrophysics.