Global evolution of the U(1) Higgs Boson: nonlinear stability and uniform energy bounds

Relying on the hyperboloidal foliation method, we establish the nonlinear stability of the ground state of the $U(1)$ standard model of electroweak interactions. This amounts to establishing a global-in-time theory for the initial value problem for a nonlinear wave-Klein-Gordon system that couples (Dirac, scalar, gauge) massive equations together. In particular, we investigate here the Dirac equation and consider a new energy functional for this field defined with respect to the hyperboloidal foliation of Minkowski spacetime. We provide a novel decay result for the Dirac equation which is uniform in the mass coefficient, and thus allows for the Dirac mass coefficient to be arbitrarily small. Furthermore we obtain energy bounds for the Higgs fields and gauge bosons that are uniform with respect to the hyperboloidal time variable.


Introduction
Main objective. Our primary objective is to study the equations of motion arising from the Higgs mechanism applied to an abelian U (1) gauge theory on a Minkowski background after spontaneous symmetry breaking. We view this model as a stepping-stone towards the full non-abelian Glashow-Weinberg-Salam theory (GSW), also known as the electroweak Standard Model. For background physics information on the models treated in this paper, see for example [1].
In short, we study here a class of nonlinear wave equations which involve the first-order Dirac equation coupled to second-order wave or Klein-Gordon equations. We are interested in the initial value problem for such systems, when the initial data have sufficiently small Sobolev-type norm. We provide here a new application of the hyperboloidal foliation method introduced for such coupled systems by LeFloch and Ma [17], which has been successfully used to establish global-in-time existence results for nonlinear systems of coupled wave and Klein-Gordon equations. This method takes its root in pioneering work by Friedrich [9,10] on the vacuum Einstein equations and by Klainerman [15,16] and Hörmander [12] on the (uncoupled) Klein-Gordon equation, as well as Katayama [13,14]. The issue of nonlinear interaction terms that couple wave and Klein-Gordon equations together, was a challenge tackled in [17] and in the subsequent developments on the Einstein equations [18,19,20] and [8,23].
The model of interest. In the abelian U (1) gauge model, the set of unknowns consists of a Dirac field ψ : R 3+1 → C 4 representing a fermion of mass m g with spin 1/2, a vector field A = (A µ ) representing a massive boson of mass m q with spin 1, and a complex scalar field χ representing the perturbation from the constant minimum φ 0 of the complex Higgs field φ = φ 0 +χ. Note the set of ground states φ 0 = φ 0 (t, x) of the theory are defined by v 2 := φ * 0 φ 0 , however in this work, as in [25], we consider only constant ground states satisfying ∂ µ φ 0 = 0. Given such a constant ground state φ 0 and the three physical parameters m q , m λ , m g , the equations of motion in a modified Lorenz gauge consist of three evolution equations and a constraint equation Here, the quadratic nonlinearities Q A ν , Q χ , and Q ψ are defined later (see (5.2)). Throughout, we use the signature (−, +, +, +) for the wave operator := η µν ∂ µ ∂ ν = −∂ 2 t + ∆, while the matrices γ µ and γ 5 are the standard Dirac matrices (cf. Section 2.1).
The nonlinear stability of the Higgs field. Our main result concerning the system (1.1) is a proof of the global-in-time existence of solutions for sufficiently small perturbations away from the constant vacuum state, defined by the conditions It is convenient to work with the perturbed Higgs field χ = φ − φ 0 as our main unknown. The initial data set are denoted by 4) and these data are said to be Lorenz compatible if (1.5) This elliptic-type system consists of two equations for 11 functions. In particular, it is easily checked that it admits non-trivial solutions, for instance with compact support. Observe also that throughout we use the convention that Greek indices take values in {0, 1, 2, 3} and Latin indices in {1, 2, 3}. In the following statement, we have m g ≥ 0, and the coefficient m −1 g is interpreted as +∞ when m g = 0.
Theorem 1.1 (Nonlinear stability of the ground state for the Higgs boson). Consider the system (1.1) with parameters m q , m λ > 0, m g ∈ [0, min(m q , m λ )] and let N be a sufficiently large integer. There exists ǫ 0 > 0, which is independent of m g , such that for all ǫ ∈ (0, ǫ 0 ) and all compactly supported, Lorenz compatible initial data (in the sense of (1.5)) satisfying the smallness condition (1.6) the initial value problem of (1.1) admits a global-in-time solution (A, χ, ψ) with, moreover, |A| ǫt −3/2 , |χ| ǫt −3/2 , |ψ| ǫ min t −1 , m −1 g t −3/2 . (1.7) Furthermore we obtain the following uniform energy estimates at highest order The assumption m g ∈ [0, min(m q , m λ )] in Theorem 1.1 can be understood both mathematically and physically. Consider the Klein-Gordon equation In order to obtain an estimate of the energy on a constant-time slice t, we use the standard energy estimate Any energy estimate for t ≥ t 0 clearly requires control over the initial term m 2 g w(t 0 ) 2 L 2 (R 3 ) . Since this term depends on the mass parameter m g , we clearly require m g to have an upper bound in order for any energy argument to be made independent of m g . What the upper bound on m g should be is clearly motivated from physics. In physical models, m q and m λ represent the masses of the gauge and Higgs fields respectively, while m g represents any fermion masses. Fermion masses are extremely small compared to both m q , m λ , and so it is particularly relevant to study the system in the entire range m g ∈ [0, min(m q , m λ )].
A simpler model: the Dirac-Proca equations. Keeping (1.1) in mind, it is interesting to first understand a simplified subset of (1.1), called the Dirac-Proca equations. In the Lorenz gauge the equations of motion for this model read where P L = 1 2 (I 4 − γ 5 ). This system describes a spinor field ψ : R 3+1 → C 4 representing a fermion of mass M with spin 1/2 and a vector field A = (A µ ) representing a massive boson of mass m with spin 1. We allow the mass parameters M, m to be positive or zero and, in particular, we will be interested in the case M = 0 (i.e. the standard Dirac-Proca equations) as well as in the case M > 0 (massive field). Observe that the Dirac-Proca system (1.10) contains no Higgs field and so the masses M, m are introduced "artificially" in the model, instead of arising from the Higgs mechanism.
An initial data set for (1.10), say, (1.12) We postpone the statement of our second main result, for this system, to Theorem 4.1 below.
Strategy of proof. The small data global existence problem for (1.10) with m g ≡ 0 and for (1.1) was solved by Tsutsumi in [24,25]. Nonetheless it is very interesting to revisit this system via the hyperboloidal foliation method developed in [17] as it provides a different perspective on the problem and sharper estimates, as detailed below. We thus introduce a hyperboloidal foliation which covers the interior of a light cone in Minkowski spacetime, and we then construct the solutions of interest in the future of an initial hyperboloid. For compactly-supported initial data this is equivalent to solving the initial value problem for a standard t=constant initial hypersurface. Although restricting to compactly supported initial data is a strong assumption, and is not assumed by Tsutsumi in [24,25], we find our results nonetheless very interesting. Furthermore the compact support restriction can be removed, subject to some long computational checks, using the Euclidian-Hyperboloidal Foliation Method developed in [20,21].
A key new insight in our work that arises from using L 2 norms defined on hyperboloids, is that we are able to derive the following energy functional on hyperboloids We can show that this energy is positive definite and controls the norm (s/t)ψ L 2 , which provides us with a notion of energy on hyperboloids and should be compared with the more standard approach based on the standard t = constant hypersurfaces; cf. for example [4]. In the case 0 ≤ M ≪ m we use the Dirac equation (expressed in a first-order form and using this energy functional) combined with Sobolev-type estimates.
Challenges for the global analysis and new insights. There are several difficulties in dealing with the system of coupled wave-Klein-Gordon equations (1.1). The first and most well-known one is that the standard Klein-Gordon equation does not commute with the scaling Killing field S = x µ ∂ µ of Minkowski spacetime, which prevents us from applying the standard vector-field method in a direct manner. Second, the nonlinearities in the Klein-Gordon equation (or the Proca equation) include a bad term describing ψ-ψ interactions, which can be regarded as a wave-wave interaction term and does not have good decay. Tsutsumi [24] was able to overcome this difficulty by defining a new variable whose quadratic nonlinearity only involves the so-called 'strong null forms' Q µj (first introduced in [11]) compatible with the scaling vector field. A similar, but complicated new variable, was also defined in [25]. Due to our use of the hyperboloidal foliation method we do not need to find new variables leading solely to strong null forms, which is an important simplification useful in our future work studying the full GSW theory.
Next, our global-in-time existence result is established under a lower regularity assumption on the initial data. Namely, the boundedness of the initial data in the norm · H N (N ≥ 20) is needed in [24], while we only require N ≥ 6. Most importantly, in the main statement in [24,25], a slow growth of the energy of all high-order derivatives occurs. By contrast, with our method of proof we do not have any growth factor in the L 2 norm (1.8) and the L ∞ norms of the Higgs scalar field and vector boson field. Note that uniform energy bounds were also established in [22]. We refer to Section 7 for a further discussion of these improved estimates.
Furthermore, it is challenging to establish a global stability result uniformly for all m g ≥ 0, while it is relative easy to complete the proof for either m g = 0 or m g a large constant. In the case where the mass is small, say m g = ǫ 2 , if we treat ψ as a Klein-Gordon field then the field decays like However with this decay we cannot arrive at the improved estimate (1/2)C 1 ǫ if we start from the a priori estimate C 1 ǫ, where C 1 is some large constant introduced in the bootstrap method. This is because the "improved" estimates we find are ǫ + C 2 1 instead of ǫ + (C 1 ǫ) 2 . Hence ψ behaves more like a wave component when m g ≪ 1, but since the mass m g may be very small but non-zero, we cannot apply techniques for wave equations. We find it possible to overcome these difficulties by analysing the first-order Dirac equation, which admits the positive energy functional (1.13), and this energy plays a key role in the whole analysis. Furthermore although we obtain logarithmic growth for the L 2 norm of our Dirac field, we obtain interesting L ∞ estimates accessible through our use of the hyperboloidal foliation. We refer to Section 3.3 and Theorem 7.6 for a further discussion of these improved estimates. One of the authors has also recently established further results concerning the zero-mass limit of Klein-Gordon equations [6].
The coupling constants. For the mass parameters (1.2) appearing in the U (1)-model we assume that m q ≃ m λ > 0, m g ≥ 0, (1.14) which depend on the coupling constants q, g and the Higgs constants λ, v. On the other hand, for the mass parameters appearing in the Dirac-Proca model we assume 1 Outline of this paper. In Section 2, we will introduce the hyperboloidal foliation method, the Dirac equation, and the energy for the Dirac component. In Sections 2.1 and 2.2 we give basic definitions for the Dirac equation and hyperboloidal foliation, the latter taken from [17]. In Section 2.3 we define the energy for the Dirac field on hyperboloids and give an energy estimate in Proposition 2.4. Note this will be complemented by other formulations of the energy functional in Section 3. Finally in Section 2.4 we convert the Dirac equation into a second order wave/Klein-Gordon equation and define the appropriate energy functional. Further properties about the hyperboloidal energy functionals are established in Section 3. Therein we provide complementing views on the Dirac energy E H . The first comes from using a Cholesky decomposition for the energy integrand in Section 3.1 and next using a Weyl decomposition for the Dirac spinor in Section 3.2. Then we study the Dirac-Proca model and the abelian model in Section 4 and Section 5, respectively. Next, in Section 6 we discuss the system of equations (1.1) and its nonlinearities are studied. Finally, in Section 7 we rely on a bootstrap argument and prove the desired the stability result.
2 Hyperboloidal energy functionals for the Dirac operator

Dirac spinors and matrices
This section is devoted to analyzing energy functionals for the Dirac equation with respect to a hyperboloid foliation of Minkowski spacetime. Recall the Dirac equation for the unknown with prescribed right-hand side F ∈ C 4 and mass M ∈ R. To make sense of this equation we need to define various complex vectors and matrices. For a complex vector z = (z 0 , z 1 , z 2 , z 3 ) T ∈ C 4 letz denote the conjugate, and z * := (z 0 ,z 1 ,z 2 ,z 3 ) denotes the conjugate transpose. If also w ∈ C 4 then the conjugate inner product is defined by The Hermitian conjugate 2 of a matrix A is denoted by A * , meaning The Dirac matrices γ µ for µ = 0, 1, 2, 3 are 4 × 4 matrices satisfying the identities where η = diag(−1, 1, 1, 1). The Dirac matrices give a matrix representation of the Clifford algebra. We will use the Dirac representation, so that the Dirac matrices take the following form where σ i 's are the standard Pauli matrices: One often uses the following product of the Gamma matrices The γ 5 matrix squares to I 4 and so we also define the following projection operators to extract the 'left-handed' and 'right-handed' parts of a spinor, that is, to extract its chiral parts. We also note the following useful identities 3 For further details on Dirac matrices and the representation considered, see [1].

Hyperboloidal foliation of Minkowski spacetime
In order to introduce the energy formula of the Dirac component ψ on hyperboloids, we first need to recall some notation from [17] concerning the hyperboloidal foliation method. We denote the point (t, x) = (x 0 , x 1 , x 2 , x 3 ) in Cartesion coordinates, with its spatial radius r : We write ∂ α := ∂ x α (for α = 0, 1, 2, 3) for partial derivatives and L a := x a ∂ t + t∂ a , a = 1, 2, 3 (2.8) for the Lorentz boosts. Throughout the paper, we consider functions defined in the interior of the future light cone K := {(t, x) : r < t − 1}, with vertex (1, 0, 0, 0). We consider hyperboloidal hypersurfaces H s : is used to denote subsets of K limited by two hyperboloids. 2 In the physics literature, A * is often denoted by A † . 3 In physics, the notationψ := ψ † γ 0 is often used, but we will avoid this here.
The semi-hyperboloidal frame is defined by Observe that the vectors ∂ a generate the tangent space to the hyperboloids. We also introduce the vector field ∂ ⊥ := ∂ t + x a t ∂ a which is orthogonal to the hyperboloids. For this semi-hyperboloidal frame above, the dual frame is given by θ 0 := dt − x a t dx a , θ a := dx a . The (dual) semi-hyperboloidal frame and the (dual) natural Cartesian frame are connected by the relation where the transition matrix (Φ β α ) and its inverse (Ψ β α ) are given by (2.11) Throughout, we use roman font E to denote energies coming from a first-order PDE (see below) and, calligraphic font E to denote energies coming from a second-order PDE. In the Minkowski background, we first introduce the energy E for scalar-valued or vector-valued maps φ defined on a hyperboloid H s : where denotes the rotational vector field. We also write E(s, φ) := E 0 (s, φ) for simplicity. All of our integrals in L 1 , L 2 , etc. are defined from the standard (flat) metric in R 3 , so (2.14)

Hyperboloidal energy of the Dirac equation
We now derive a hyperboloidal energy for the Dirac equation (2.1). Premultiplying the PDE (2.1) by ψ * γ 0 gives The conjugate of (2.1) is Multiplying this equation by ψ gives Adding (2.15) and (2.16) together yields Note the mass term does not appear in (2.17). Moreover recalling 2Re[z] = z +z for some z ∈ C then we see that (2.17) is the real part of (2.15). It would appear however if we subtracted (2.15) from (2.16), that is the imaginary part of (2.15), then we find However such an expression does not appear to be useful. We return to (2.17) and integrate over regions in spacetime to obtain energy inequalities. These give the following two definitions of first-order energy functionals. Integrating (2.17) over [t 0 , t] × R 3 , and assuming spatially compactly supported initial data, gives the following result.
Then it holds Such functionals on a constant-time foliation have been considered frequently in the literature, see for instant [7]. The following Lemma gives a new perspective using a hyperboloidal foliation.

Lemma 2.2. On hyperboloidal slices H s define the energy
For solutions to (2.17) this satisfies (2.22) Here n and dσ are the unit normal and induced Lebesgue measure on the hyperboloids respectively Using this explicit form of n and dσ gives the result.
which in particular implies the positivity of the energy and re-arranging gives (2.24b).
Using the above positivity property we can now establish the following energy inequality.
Proposition 2.4 (Hyperboloidal energy estimate for Dirac spinor). For the massive Dirac spinor described by (2.1) the following estimate holds Proof. Differentiating (2.21) with respect to s yields

Hyperboloidal energy based on the second-order formulation
Finally in this section we will convert the Dirac equation into a second-order PDE and define associated energy functionals. Apply the first-order Dirac operator −iγ ν ∂ ν to the Dirac equation (2.1) and use the identity (2.2) to obtain Substituting the PDE into the bracketed term gives the following second-order PDE Proof. The conjugate of (2.28a) reads ψ * − M 2 ψ * = G * . Using −∂ t ψ * and −∂ t ψ as multipliers on (2.28a) and its conjugate respectively we obtain Adding these two equations together gives Integrating this equation in the region K [s 0 ,s] , we have Using the explicit form of n and dσ given in (2.23) and noting that 2Re We next use the change of variable formula ∂ i = ∂ i − (x i /t)∂ t to rewrite the energy term: We can now estimate the nonlinearity on the RHS using the change of variables τ = (t 2 −r 2 ) 1/2 and dtdx = (τ /t)dτ dx. In particular we have Thus by differentiating (2.29) and using the above we have Integrating this expression over [s 0 , s] gives the desired result.
We end with a short remark here on positive and negative mass spinors. In the original first-order Dirac equation (2.1) the mass M was defined as a real parameter with no sign restriction. Furthermore the mass did not appear in the hyperboloidal energy E H defined in (2.20). This implies that spinors with equal masses but opposite signs (±M ) would still obey the same energy estimates of Proposition 2.4. This is consistent with the second-order equation (2.27) for the Dirac field. In this equation the mass appears squared, so spinors with equal masses, but of opposite signs, obey the same second order equation. Moreover the mass M 2 appears in the second-order hyperboloidal energy expression E M in (2.28c).
Thus either the mass M should not appear in E H , as we have found, or if the mass M were to appear in E H , it would necessarily need to be invariant under a sign change, so as to agree with the second-order energy estimates involving E M .

Additional properties of Dirac spinors on hyperboloids 3.1 Hyperboloidal energy based on a Cholesky decomposition
Our first task is to obtain a hyperboloidal energy for the Dirac field ψ expressed in terms of a product of complex vectors z(ψ) * z(ψ). Such an expression is then easily seen to be positive semi-definite which, clearly, is in contrast to the form given in (2.20) and Proposition 2.3.
Recall the standard Cholesky decomposition: any Hermitian, positive-definite matrix A can be decomposed in a unique way as where P is a lower triangular matrix with real and positive diagonal entries. In particular if A is positive semi-definite then the decomposition exists however one loses uniqueness and the diagonal entries of P may be zero. We now prove the following result. and specifically

3)
which can also be expressed as P = (s/t) + 1 2 The above expression is quite natural and resembles what is known for the wave equation: the factor x i /t comes from Stokes' theorem applied to hyperboloids and we cannot expect to fully control the standard L 2 norm, namely Hs ψ * ψ dx. Proof.
Step 1. Existence of the decomposition. Before we proceed with the derivation of the identity, we present an argument showing that such a decomposition exists by proving positive semi-definiteness. For simplicity of notation, let N i := x i /t. The integrand of E H (s, ψ) can be written as ψ * Aψ where A := I 4 + N j γ 0 γ j . Here the spatial indices are contracted with δ ij , so that N j γ j = δ ij N i γ j . Note A is hermitian since Then for all z ∈ C 4 we have We used that N j N j = (r/t) 2 ≤ 1 which holds in the light-cone K. Thus A is positive semi-definite.
Step 2. Computing the matrix P . With respect to the Dirac representation (2.3) we have Here I 2 and 0 2 represent the 2 × 2 identity and zero matrices respectively. Calculate the second term above using the Pauli matrices: Define ω := N 1 + iN 2 and recall N i ∈ R. Thus we have Consider now 2 × 2 complex matrices B, C, D such that This implies the following identities If we let D = I 2 and C = N 3ω ω −N 3 then we must solve

Hyperboloidal energy based on the Weyl spinor representation
Yet one more approach in deriving energy estimates is obtained by expressing the Dirac spinors in terms of Weyl spinors and then studying the energy of Weyl spinors (3.7) instead. This provides another convenient way to study Dirac equations. Decompose the spinor ψ and source term F as where u, v : R 1+3 → C 2 are Weyl spinors and F ± ∈ C 2 . Defining ∂ ± := ∂ 0 ± σ i ∂ i the PDE (2.1) can be shown to be equivalent to A Dirac-Klein-Gordon system with respect to such a Weyl spinor decomposition has been studied, albeit in the low-regularity setting, by Bournaveas [4]. Following a similar approach to Section 2.3 we find an analogous hyperboloidal Weyl spinor energy: Similar to Propositions 2.3 and 2.4 we can prove positivity and an energy estimate for E σ ± .
Proposition 3.2. For a C 2 -valued function w the following holds: Proof.
Step 1. Using the Dirac representation (2.3) and the decomposition (3.5), the PDE (2.1) becomes Defining ∂ ± := ∂ 0 ± σ i ∂ i this becomes Adding and subtracting the two rows above gives the following Following a similar approach to deriving (2.17), we multiply the first and second equation by v * and u * respectively.
One then adds these equations to their conjugate to obtain the following: Note the mass terms appear above. However if add these equations together we obtain which does not contain a term involving M . This equation is the analogous Weyl spinor version of (2.17). Clearly integrating (3.12) over K [s 0 ,s] gives the energy functional E σ ± (s, u) defined in (3.7).
Step 2. Next we establish that for a C 2 -valued function w the following holds: The idea is in the spirit of Proposition 2.3. Observe that the sigma matrices are Hermitian and satisfy the following anti-commutator relation: {σ i , σ j } = 2δ ij I 2 . Then we have Thus we have Step 3. Next, let us show that the following hyperboloidal energy inequality holds for the Weyl spinor equation (3.6) Namely, integrating (3.12) over K [s 0 ,s] we obtain Differentiating in s and noting that (s/t)u, (s/t)v L 2 (Hs) ≤ u, v L 2 (Hs) ≤ (E σ

Sobolev-type estimates for a Dirac spinor
In this section we will obtain novel decay estimates for the Dirac spinor. First we recall the following standard Sobolev estimate on hyperboloids. where the summation is over Lorentz boosts L.
This estimate will be stated again in Section 7 and proved for general functions, see Proposition 7.4. In this section, we will adapt this Sobolev-type inequality (3.13) to include boosts which commute with the Dirac operator iγ ν ∂ ν . Following Bachelot [3] we introduce modified boosts L a that differ from L a by a constant matrix (3.14) It then holds that [ L a , iγ ν ∂ ν ] = 0. The following result is a simple extension of Proposition 3.3. where L denotes a modified Lorentz boost.
We can combine this Sobolev estimate (3.15) with an appropriate energy bound in a standard way to obtain the following 'weak' decay estimate for the Dirac field.
Proof. It is not hard to show Then together with the Sobolev inequality on hyperboloids (3.13), it gives the desired pointwise decay estimate (3.16).
Since t < s 2 within the cone K, the decay rate ( We notice that there is no information about the mass M in the Dirac energy E H . Nonetheless by closely studying each component of the functional E H (s, ψ), we obtain the following estimate.
Proposition 3.6. Consider a solution to the Dirac equation
Step 1. We first express the Dirac equation (2.1) in the semi-hyperboloidal frame to get It is easy to see that the first and the third terms are bounded by Then recalling the Sobolev inequality on hyperboloids (7.10), it suffices to show Step 2. We note that Proposition 2.3 implies which suggests us to study the commutator L a , γ 0 − γ a (x a /t) .
A computation gives which only contains good terms. An induction shows L J , γ 0 − γ a (x a /t) with |J| ≤ 2 only contains good terms. Hence the proof is complete.
4 Nonlinear stability of the ground state for the Dirac-Proca model 4

.1 The Dirac-Proca model as a PDE system
Using the tools of Section 2, we will now discuss the Dirac-Proca model, the gauge condition and choice of initial data. This leads us to the second main stability Theorem 4.1 below, which can be proved using the methods of Sections 6 and 7. Without fixing a gauge, the field equations for the Dirac-Proca model with unknowns A µ and ψ read Recall P L = 1 2 (I 4 − γ 5 ) was defined in (2.6). Here ψ : R 3+1 → C 4 represents a massless Dirac fermion with spin 1/2 and A µ : R 3+1 → R represents a massive boson (the Proca field) of mass m 2 with spin 1. As discussed in [24], (4.1) is equivalent to the following system It can easily be shown that ∂ µ A µ satisfies a homogeneous Klein-Gordon equation. Thus we may treat ∂ µ A µ = 0 as a constraint provided we specify the initial data set A ν 0 , A ν 1 , ψ 0 at some time t 0 > 0: to satisfy the following two 'Lorenz compatibility' conditions For more generality, see also [3], we consider (4.1) with an artificial mass M added to the spinor: The mass parameters M ≥ 0 and m > 0 are constants, and we will study both cases M = 0 and M > 0. Again the initial data will be taken to satisfy (4.4). This elliptic-type system of two equations for nine scalar functions admits non-trivial compactly supported solutions. For example one may choose: A j 0 ∈ C c (R 3+1 ) for j = 1, 2, 3, ψ 0 ∈ C c (R 3+1 ) such that each component of ψ 0 is the same, and all remaining initial data is trivial.

Main result for the Dirac-Proca model
We now state our result for the Dirac-Proca model. The proof of the theorem below will clearly follow from the the proof we will develop for our main result (Theorem 7.1), and so we omit it.
For the model under consideration in this section, the ground state of the theory is simply given by Theorem 4.1. Consider the Dirac-Proca system (1.10) with M ≥ 0, m > 0, and let N be a sufficiently large integer. There exists ǫ 0 > 0, which is independent of M , such that for all ǫ ∈ (0, ǫ 0 ) and all compactly supported, Lorenz compatible initial data in the sense of (1.12), satisfying the smallness condition the initial value problem of (1.10) admits a unique global-in-time solution (ψ, A ν ). Furthermore, the following decay result holds 5 Nonlinear stability of the ground state for the U (1) model

The U(1) model as a PDE system
We now treat the Higgs mechanism applied to a U (1) abelian gauge field. This gives some exposure to the problems coming from the Higgs field that we will meet when dealing with the full GSW model in future work. From the Lagrangian of this theory; see (5.6) below, and in a suitably modified Lorenz gauge; see (5.5) below, the field equations for this model read with quadratic or higher order terms given by Here χ := φ − φ 0 : R 1+3 → C is the perturbation from a constant vacuum state φ 0 satisfying φ * 0 φ 0 = v 2 and ∂ µ φ 0 = 0. The field ψ : R 3+1 → C 4 represents a fermion of mass m g with spin 1/2 and A µ : R 1+3 → R represents a massive boson of mass m 2 q with spin 1. Furthermore the mass coefficients m 2 q = 2q 2 v 2 , m 2 λ = 4λv 2 , m g = gv 2 (5.3) depend themselves on given coupling constants denoted by q, g, λ, as well as the so-called "vacuum expectation value" of the Higgs field, denoted by v.
The initial data set are denoted by which are said to be 'Lorenz compatible' if they satisfy The derivation of (5.5) will follow from results of the following Section 5.2, in particular the gauge choice condition given in Lemma 5.3. Similar to the constraint equations (4.4) for the Dirac-Proca model, we also obtain in (5.5) an elliptic-type system of only two equations for eleven functions. Clearly non-trivial solutions with compact support can be constructed. Our main result was stated in the introduction (Theorem 1.1) and the proof will be provided in Section 7.

The abelian action and U(1) invariance
The Lagrangian we consider is where we use the following definitions for the Higgs potential, gauge curvature and gauge covariant derivatives: Furthermore λ, v, g, q are constants. A calculation shows that the Euler-Lagrange equations for (5.6) are the following together with (5.7). These PDEs can also be expressed as Proof. The U (1) gauge symmetry induces the following transformations where α = α(t, x) is some arbitrary function of space and time. This implies By the commutation of partial derivatives we immediately see that F µν is invariant under U (1) gauge transformations. Inserting these transformations into the Lagrangian one obtains Thus L is invariant under the transformation (5.10).

Propagation of an inhomogeneous Lorenz gauge
Next, with a similar aim to that of Section 4.1, we will turn the PDE (5.9) into one of definite type by specifying a particular gauge for the vector field A µ .
Lemma 5.2 (The inhomogeneous Lorenz gauge). Let X be a suitably regular scalar field. Consider the modified system Suppose the initial data for (5.13) satisfy Then as long as the solution to (5.13) exists, it will satisfy div A = X.
Proof. To propagate the gauge choice div A = X imposed on the initial data we will take the divergence of the first equation in (5.13) and use the second evolution equation (5.13): where in the last line we used the specific Higgs potential (5.7).

Gauge choice for the abelian model
The Higgs potential has a non-zero minimum at φ 0 := ve iθ/v where and θ : R 1+3 → R is arbitrary. There is ambiguity in this minimum state due to the U (1) symmetry represented by θ. A particular choice of θ will break this U (1) symmetry, and such a scenario is termed 'spontaneous symmetry breaking'. We consider perturbations of the form where φ 0 is constant in space and time and has magnitude |φ 0 | = v. The following result is a consequence of Lemma 5.2 by choosing X = −iq(φ * 0 χ − χ * φ 0 ). Lemma 5.3. Suppose the initial data satisfy the following gauge condition Then the Euler-Lagrange equations for (5.6) are equivalent to those written in (5.1).
6 Structure and nonlinearity of the models 6

.1 Aim of this section
In this section we will transform the variables in (5.1) in order to treat difficult terms in their nonlinearities. The nonlinearities, defined in (5.2), which require additional work to control are of the following two types ψ * γ 0 γ ν ψ, ψ * γ 0 ψ, (6.1) appearing in Q A ν and Q χ respectively. If the mass of the Dirac spinor is small or zero, that is 0 ≤ m g ≪ min(m q , m λ ), then the nonlinearities in (6.1) will have insufficiently fast decay for the bootstrap argument to close1. To address these issues, we employ a transformation introduced by [24], see (6.3) and (6.9) below. Finally although the equation satisfied by χ is of ambiguous type, we find that χ can be decomposed into two components, with each component satisfying a Klein-Gordon equation.

Hidden null structure from Tsutsumi
Following Tsutsumi [24] we can uncover a null structure using a particular transformation. Significantly however, we will not need to reduce to the 'strong null forms' via any complicated transformations, as in [25, eq (2.6)]. For complex-valued functions Φ(t, x), Ψ(t, x) : R 3+1 → C n , recall the null form

Define a new variable
This satisfies the non-linear Klein-Gordon equation in which the nonlinearities are where (2.27) and (5.2) were used to compute G ψ . Note the nonlinearity ψ * γ 0 γ ν ψ in (6.4) now appears with a good factor m 2 g . We will control this term using Propositon 3.6.

Decomposition of χ
In order to study the behaviour of χ, it is more convenient to consider equations for the following two variables Since the following identity holds it is equivalent to estimate either χ or χ ± . These new variables satisfy the Klein-Gordon equations with

Following Tsutsumi again [24], we define the new variable
which satisfies the following Klein-Gordon equation with the nonlinearity The final term 2m g ψ * γ 0 ψ here cancels the problematic term in Q χ + and any other nonlinearities of the form ψ * γ 0 ψ now appear with a good factor of m 3 g in front. Note that in the case 0 ≤ m g ≪ min(m q , m λ ), the second order formulation (2.27) of the Dirac equation is more like a nonlinear wave equation. In this case we do not have good bounds for either the L 2 f or L ∞ norm of ψ and this is why the null structure of (6.4) and the factor m 2 g in front of the ψ-ψ interaction are needed. In the case m g ∼ min(m q , m λ ), better Klein-Gordon bounds are available and we do not require such a particular structure anymore.

Overview
This section is devoted to using a bootstrap argument to prove Theorem 1.1. After the treatment of the equations in Section 6, we remind one that in the following analysis we will deal with the unknowns We also provide here the following theorem, which details the L 2 and L ∞ estimates of the unknowns.
2) and the following L ∞ estimates sup (t,x)∈Hs 3) where C 1 is some constant to be determined later.
Strategy of the proofs of Theorems 1.1 and 7.1. The proofs are based on a bootstrap argument. In Section 7.2 we recall standard estimates for null terms, various commutators and Sobolev-type estimates on hyperboloids. The bootstrap assumptions will be made in (7.14). These bootstraps, combined with some standard commutator estimates and Sobolev-type inequalities, will lead to certain weak estimates in (7.16) and (7.17). In Section 7.4 we will use our first-order hyperboloidal energy to upgrade our estimates for the Dirac component, namely in Theorem 7.6 and Corollary 7.7. In Section 7.5 we obtain estimates for the transformed variablesÃ ν andχ + defined above. Putting all of this together we finally are able to close our bootstrap assumptions.

Standard Estimates: null forms, commutators and pointwise estimates
We first illustrate estimates for the quadratic null terms, which, roughly speaking, reveal the following bounds ∂ α u∂ α v s 2 /t 2 ∂u ∂v .
Lemma 7.2. We have the following estimate for the quadratic null term T αβ ∂ α u∂ β v with constants T αβ and u, v sufficiently regular.
One refers to [17] for the proof. We next recall the estimates for commutators, also proved in [17]. Heuristically speaking, the following lemma provides the relation in which Γ takes values from the set Lemma 7.3. Assume a function u defined in the region K is regular enough, then we have [∂ I L J , ∂ α ]u ≤ C(|I|, |J|) [ Furthermore we can obtain L ∞ estimates by recalling the following Sobolev-type inequality on hyperboloids [17]. where the summation is over Lorentz boosts L. Note the implied constant is uniform in s ≥ 2, and one recalls that t = s 2 + |x| 2 on H s .
Lemma 7.5 (Gronwall-type inequality). Let u(t) be a non-negative function that satisfies the integral inequality where b(t) is non-negative function for t ≥ t 0 . Then it holds (7.13)

Bootstrap assumptions and basic estimates
By the local well-posedness of semilinear PDEs, there exists an s 1 > s 0 in which the following bootstrap assumptions hold for all s ∈ [s 0 , s 1 ] E mg (s, ∂ I L J ψ) 1/2 ≤ C 1 ǫ log s, |I| + |J| = N. (7.14) If we can prove the refined estimates then we are able to assert that s 1 cannot be finite, which in turn implies a global existence result for (1.1).
Combining the bootstrap assumptions (7.14) with the estimates for commutators in Lemma 7.3, the following sets of estimates are obtained: Combining these estimates with Proposition 7.4 and Proposition 3.6 the following hold: (7.17)

First-order energy estimate for the Dirac field
To obtain decay estimates for the Dirac component ψ, a standard method is to analyse the second-order form of the Dirac equation (2.27). This is then a semilinear Klein-Gordon equation with mass m 2 g and so there are now standard techniques to estimate the nonlinearity; see for example [2] and [18]. However, the right-hand side term appearing in our wave equation (2.27) does not decay sufficiently fast for this argument to close, which is due to the possibly vanishing mass m 2 g ≥ 0. Thus at this point we recall Proposition 2.3 and the lower bound (2.24c) for the energy E H . This motivates us to analyse the first-order form of the Dirac equation in the following Theorem to obtain certain improved L 2 and L ∞ estimates for ψ that are uniform in m g . As a consequence, one has the following sup-norm estimate for ψ: sup Hs t∂ I L J ψ ǫ + (C 1 ǫ) 2 , |I| + |J| ≤ N − 2. (7.20) Proof.
As a consequence, we have the following sup-norm estimates for ψ.