arxiv: v1 [math.ap] 7 Feb 2019

arxiv:1902.02685v1 [math.ap] 7 Feb 2019 Global evolution of the U1 Higgs Boson: nonlinear stability and uniform energy bounds Shijie Dong, Philippe G. LeFloch, and Zoe Wyatt February 2019 Abstract Relying on the hyperboloidal foliation method, we establish the nonlinear stability of the ground state of the so-called U1 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 operator and its energy functional defined with respect to the hyperboloidal foliation of Minkowski spacetime. We also provide a decay result for the Dirac equation which is uniform in the mass coefficient, and thus allow for the mass coefficient to be arbitrarily small. Our energy bounds are uniform modulo a logarithm growth with respect to the hyperboloidal time variable. Contents 1 Introduction 2 2 Hyperboloidal energy functionals for the Dirac operator 5 3 Additional properties of Dirac spinors on hyperboloids 12 4 Nonlinear stability of the ground state for the Dirac-Proca model 16 5 Nonlinear stability of the ground state for the U1 model 17 6 Structure and nonlinearity of the models 20 7 Bootstrap argument 22 Laboratoire Jacques-Louis Lions, Centre National de la Recherche Scientifique, Sorbonne Université, 4, Place Jussieu, 75252 Paris, France. Email: dongs@ljll.math.upmc.fr, contact@philippelefloch.org Maxwell Institute for Mathematical Sciences, School of Mathematics, University of Edinburgh, Edinburgh, EH9 3FD, United Kingddom. Email: zoe.wyatt@ed.ac.uk. 1

1 Introduction Main objective. Our primary objective is to study the equations of motion arising from the Higgs mechanism applied to an abelian U1 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 theorygsw, 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 [16], 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 [8, 9] on the vacuum Einstein equations and by Klainerman [14, 15] and Hörmander [11] on the uncoupled Klein-Gordon equation, as well as Katayama [12, 13]. Taking into account nonlinear interaction terms that couple wave and Klein-Gordon equations together, was a challenge tackled in [16] and in the subsequent developments on the Einstein equations [17, 18, 19] and [7, 21]. The model of interest. In the abelian U1 gauge model, the set of unknowns consists of a spinor 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 +χ. Given the ground state φ 0 with norm denoted by v 2 := φ 0 φ 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 m 2 q A ν = Q A ν, χ m 2 q and a constraint equation φ 0 2v 2 φ 0 χ χ φ 0 m 2 λ φ 0 2v 2 φ 0 χ χ φ 0 = Qχ, iγ µ µ ψ m g ψ = Q ψ, 1.1a diva+i m q v φ 0 χ χ φ 0 = 0. 1.1b 2 Here, thequadraticnonlinearitiesq A ν, Q χ,andq ψ aretobedefinedlater, andφ 0 isaconstant complex number with magnitude φ 0 = v. 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. Some important physical parameters are provided together with the mass coefficients above, namely m 2 q := 2q2 v 2 > 0, m 2 λ := 4λv2 > 0, m g := gv 2 0, 1.2 themselves depend on given coupling constants, denoted by λ, g, q, as well as the so-called vacuum expectation value of the Higgs field denoted by v. 2

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 A µ 0, φ φ 0, ψ 0. 1.3 It is convenient to work with the perturbed Higgs field χ = φ φ 0 as our main unknown. The initial data set are denoted by A ν,χ,ψ t 0, = A ν 0,χ 0,ψ 0, t A ν, t χ t 0, = A ν 1,χ 1, 1.4 and these data are said to be Lorenz compatible if a A a 0 = A0 1 iq φ 0 χ 0 χ 0 φ 0, A 0 0 m 2 qa 0 0 = i A i 1 iq φ 0χ 1 χ 1φ 0 +iq χ 0 χ 1 χ 1χ 0 +2q 2 A 0 0 φ 0 χ 0 +χ 0φ 0 +χ 0χ 0 +qψ 0 ψ 0. 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 and m g [0,minm 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 A 0,χ 0,ψ 0 H N R 3 + A 1,χ 1 H N 1 R 3 ǫ, 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 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 A ν m 2 A ν = ψ γ 0 γ ν P L ψ, iγ µ µ ψ +Mψ = γ µ A µ P L ψ, 1.8 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.8 contains no Higgs field and so the masses M,m are introduced artificially in the model, instead of arising from the Higgs mechanism. 3

A set of initial data set for 1.8, say, A ν, t A ν,ψ t 0, = a ν,b ν,ψ 0, 1.9 is called Lorenz compatible if one has b 0 + j a j = 0, a 0 m 2 a 0 = ψ 0P L ψ 0 j b j. 1.10 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.8 with m g 0 and for 1.1 was solved in Tsutsumi in [22, 23]. Yet, it is very interesting to revisit this system via the hyperboloidal foliation method developed in [16] as it provides a different perspective on the problem and somewhat sharper estimates. 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. In the case of the Dirac-Proca model, our theorem refines the results in [22] in several directions. First, the global-in-time existence 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 [22], while we only require N 6. Most importantly, in the main statement in [22], a slow growth of the energy of high-order derivatives may take place. With our method of proof, with the exception of the highest order term, we do not have any growth factor in the or L norms. We refer to Section 7 for a further discussion of these improved estimates. In contrast to [22], for simplicity we assume that the initial data are compactly supported. Uniform energy bounds were also established in [20]. Moreover, inthecase0 M mweanalyzethediracequationexpressedinafirst-order form and derive the following energy functional on hyperboloids E,ψ := ψ ψ x i t ψ γ 0 γ i ψ dx. 1.11 We can show that this energy is positive definite and controls the norm s/tψ, 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]. Challenges for the global analysis. 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 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 [22] 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 [10] compatible with the scaling vector field. A similar, but complicated 4

new variable, was also defined in [23]. Due to our use of the hyperboloidal foliation method we do not need to find new variables leading solely to strong null forms. Furthermore, it is challenging toestablish aglobal stability resultuniformlyforall 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 ψ C 1 ǫ 1 t 3/2. However with this decay we cannot arrive at the improved estimate 1/2C 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 a positive energy functional, and this energy plays a key role in the whole analysis. The coupling constants. For the mass parameters 1.2 appearing in the U1 model we assume that m q m λ > 0, m g 0, 1.12 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 m > 0, M 0. 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 [16]. 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 2.1 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 1 There isnohiggs fieldtogenerate themasses, andwedonotneedtousesubscripts for themass coefficients. 5

ψ C 4 iγ µ µ ψ +Mψ = F, 2.1 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 let z 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 4 z,w := z w = z α w α. α=1 The Hermitian conjugate 2 of a matrix A is denoted by A, meaning A αβ := Ā βα. The Dirac matrices γ µ for µ = 0,1,2,3 are 4 4 matrices satisfying the identities {γ µ,γ ν } := γ µ γ ν +γ µ γ ν = 2η µν I, γ µ = η µν γ ν, 2.2 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 γ 0 I2 0 =, γ i 0 σ i = 0 I 2 σ i, 2.3 0 where σ i s are the standard Pauli matrices: σ 1 0 1 :=, σ 2 := 1 0 0 i i 0, σ 3 := 1 0. 2.4 0 1 One often uses the following product of the Gamma matrices γ 5 := iγ 0 γ 1 γ 2 γ 3 0 I2 =. 2.5 I 2 0 The γ 5 matrix squares to I 4 and so we also define the following projection operators P L := 1 2 I 4 γ 5, P R := 1 2 I 4 +γ 5 2.6 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 γ µ = γ 0 γ µ γ 0, γ 0 γ µ = γ 0 γ µ, {γ 5,γ µ } = 0. 2.7 For further details on Dirac matrices and the representation considered, see [1]. 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. 6

2.2 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 [16] 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 := x = x 1 2 +x 2 2 +x 3 2. 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 := {t,x : t 2 r 2 = s 2 } with s > 1. Also K [s0,s 1 ] := {t,x : s 2 0 t2 r 2 s 2 1 ;r < t 1} is used to denote subsets of K limited by two hyperboloids. The semi-hyperboloidal frame is defined by 0 := t, a := L a t = xa t t + a. 2.9 Observethat thevectors a generate thetangent spacetothehyperboloids. We alsointroduce the vector field := t + xa t a which is orthogonal to the hyperboloids. For this semi-hyperboloidal frame above, the dual frame is given by θ 0 := dt xa t dxa, θ a := dx a. The dual semi-hyperboloidal frame and the dual natural Cartesian frame are connected by the relation α = Φ α α α, α = Ψ α α α, θ α = Ψ α α dxα, dx α = Φ α α θα, 2.10 where the transition matrix Φ β α and its inverse Ψ β α are given by 1 0 0 0 1 0 0 0 Φ β α = x 1 /t 1 0 0 x 2 /t 0 1 0, Ψβ α = x 1 /t 1 0 0 x 2 /t 0 1 0. 2.11 x 3 /t 0 0 1 x 3 /t 0 0 1 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 : E m s,φ := t φ 2 + a φ 2 +x a /t t φ a φ+ t φ a φ +m 2 φ 2 dx a = s/t t φ 2 + a φ 2 +m 2 φ 2 dx 2.12 a = φ 2 + s/t a φ 2 + t 1 Ω ab φ 2 +m 2 φ 2 dx, a a<b where Ω ab := x a b x b a 2.13 denotes the rotational vector field. We also write Es,φ := E 0 s,φ for simplicity. All of our integrals in L 1,, etc. are defined from the standard flat metric in R 3, so φ L 1 f H = s φ ds := φ H s 2 +r 2,x dx. 2.14 s R 3 7

2.3 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 ψ 0 ψ +ψ γ 0 γ j j ψ +imψ γ 0 ψ = iψ γ 0 F. 2.15 The conjugate of 2.1 is Multiplying this equation by ψ gives Adding 2.15 and 2.16 together yields µ ψ γ µ iψ M = if. 0 ψ ψ + j ψ γ 0 γ j ψ imψ γ 0 ψ = if γ 0 ψ. 2.16 0 ψ ψ+ j ψ γ 0 γ j ψ = iψ γ 0 F if γ 0 ψ. 2.17 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 ψ 0 ψ 0 ψ ψ+ψ γ 0 γ j j ψ j ψ γ 0 γ j ψ +2iMψ γ 0 ψ = iψ γ 0 F +if γ 0 ψ. 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. Lemma 2.1. On t = const slices, define the energy E flat t,ψ := ψ ψt,xdx. 2.18 R 3 Then it holds E flat t,ψ = E flat t 0,ψ+ t t 0 R 3 iψ γ 0 F if γ 0 ψ dxdt. 2.19 Such functionals on a constant-time foliation have been considered frequently in the literature, see for instant [6]. The following Lemma gives a new perspective using a hyperboloidal foliation. Lemma 2.2. On hyperboloidal slices define the energy E,ψ := ψ ψ x i t ψ γ 0 γ i ψ dx. 2.20 For solutions to 2.17 this satisfies E,ψ = E 0,ψ+ s s 0 H s s/tiψ γ 0 F if γ 0 ψdxd s. 2.21 8

Proof. Integrating 2.17 over K [s0,s] gives ψ ψ,ψ γ 0 γ j ψ ndσ ψ ψ,ψ γ 0 γ j ψ ndσ 0 = iψ γ 0 F if γ 0 ψdtdx. K [s0,s] 2.22 Here n and dσ are the unit normal and induced Lebesgue measure on the hyperboloids respectively n = t 2 +r 2 1/2 t, x i, dσ = t 1 t 2 +r 2 1/2 dx. 2.23 Using this explicit form of n and dσ gives the result. To our knowledge the hyperboloidal Dirac energy 2.20 is new. We now show that this energy E,ψ is indeed positive definite. Proposition 2.3 Hyperboloidal energy of the Dirac equation. I. By defining E + s,ψ := ψ x i t γ0 γ i ψ ψ x j t γ0 γ j ψ dx, one has E,ψ = 1 2 E+ s,ψ+ 1 2 which in particular implies the positivity of the energy E,ψ 1 2 Hs s 2 Hs s 2 t 2ψ ψdx, t 2ψ ψdx 0. 2.24a 2.24b 2.24c Proof. Expanding out the bracket in E + s,ψ gives E + s,ψ := = = = ψ ψ 2 x i ψ ψ 2 x i ψ ψ 2 x i and re-arranging gives 2.24b. t ψ γ 0 γ i ψ +ψ xi x k t 2 γ0 γ j γ 0 γ l ψδ ij δ kl dx t ψ γ 0 γ i ψ ψ x ix j t 2 γi γ j ψ dx dx t ψ γ 0 γ i ψ +ψ ψ xi x j t 2 δ ij ψ ψ 2 x i t ψ γ 0 γ i ψ +ψ ψ t2 s 2 t 2 dx = 2E,ψ Hs s 2 t 2ψ ψdx, 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 E,ψ 1/2 E 0,ψ 1/2 + s F f H s d s. 2.25 s 0 9

Proof. Differentiating 2.21 with respect to s yields E H s,ψ 1/2 d d s EH s,ψ 1/2 1 s/t F γ 0 ψ + ψ γ 0 F dx 2 H s s/tψ f H s F f H s. Recalling in Proposition 2.3, we have s/tψ f H s EH s,ψ 1/2. Thus we have d d s EH s,ψ 1/2 F f H s, and the conclusion follows by integrating over [s 0,s]. 2.4 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 η µν µ ν ψ +M iγ ν ν ψ = iγ ν ν F. 2.26 Substituting the PDE into the bracketed term gives the following second-order PDE ψ M 2 ψ = MF iγ ν ν F. 2.27 This provides us with another approach for deriving an energy estimate for the Dirac equation. We now check the hyperboloidal energy coming from 2.27. Lemma 2.5 Second-order hyperboloidal energy estimate for the Dirac equation. For all solution ψ C 4 of ψ M 2 ψ = G, 2.28a one has where E M s,ψ := s E M s,ψ 1/2 E M s 0,ψ 1/2 + G H τdτ, s 0 t ψ 2 + i i ψ 2 +M 2 ψ 2 +2 i x i t Re[ t ψ i ψ ] dx. 2.28b 2.28c 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 t ψ 2 t ψ tψ i 2 i ψ +M2 t ψ ψ = t ψ G, 2 t ψ t ψ i 2 i ψ t ψ +M 2 t ψ ψ = G t ψ. Adding these two equations together gives t t ψ t ψ + i ψ i ψ +M 2 ψ ψ i i t ψ i ψ + t ψ i ψ = t ψ G+G t ψ = 2Re [ G t ψ ]. 10

Integrating this equation in the region K [s0,s], we have t ψ 2 + i 0 = 2 t ψ 2 + i K [s0,s] i ψ 2 +M 2 ψ 2, t ψ i ψ + i ψ t ψ ndσ Re [ G t ψ ] dtdx. i ψ 2 +M 2 ψ 2, t ψ i ψ + i ψ t ψ ndσ Using the explicit form of n and dσ given in 2.23 and noting that 2Re[ i ψ t ψ] = t ψ i ψ+ i ψ t ψ we find s E M s,ψ E M s 0,ψ = 2 Re [ G t ψ ] dtdx. 2.29 H t We next use the change of variable formula i = i x i /t t to rewrite the energy term: = t ψ 2 + i s/t t ψ 2 + i s 0 i ψ 2 +M 2 ψ 2 + xi t ψ i ψ + i ψ t ψ dx t i ψ 2 +M 2 ψ 2 dx. Wecannowestimatethenonlinearity ontherhsusingthechangeofvariablesτ = t 2 r 2 1/2 and dtdx = τ/tdτdx. In particular we have 2 Re [ G t ψ ] τ/tdx 2 τ/t t ψ H τ G H τ H τ 2E M τ,ψ 1/2 G H τ. Thus by differentiating 2.29 and using the above we have d dτ E Mτ,ψ 1/2 G H τ. 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 equation2.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. 11

3 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 A = P P, 3.1 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. Proposition 3.1 Hyperboloidal energy for the Dirac equation. II. There exists a lower triangular matrix P with real and positive diagonal entries such that E,ψ = Pψ Pψdx, 3.2 and specifically P = which can also be expressed as s/t 0 0 0 0 s/t 0 0 x 3 /t x 1 /t ix 2 /t 1 0, 3.3 x 1 /t+ix 2 /t x 3 /t 0 1 P = s/t+1 2 I 4 + s/t 1 γ 0 x i +δ ij 2 t γ0 γ j. 3.4 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 norm, namely ψ ψ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,ψ 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 A = I 4 +N j γ 0 γ j = A. Also N j γ 0 γ j N k γ 0 γ k = N j N k γ 0 γ 0 γ j γ k = N j N k γ j γ k = N j N j I 4. Then for all z C 4 we have 0 Az Az = 1+N j N j z I 4 z +2z N j γ 0 γ j z 2 z I 4 z +z N j γ 0 γ j z = 2z Az. 12

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 I2 0 A = I 4 +N 2 02 σ j j 0 2 I 2 σ j 0 2 I2 0 = 2 02 σ +N j 0 2 I j 2 σ j. 0 2 Here I 2 and 0 2 represent the 2 2 identity and zero matrices respectively. Calculate the second term above using the Pauli matrices: 0 σ j 0 1 0 i 1 0 N3 N N j σ j = N 0 1 +N 1 0 2 +N i 0 3 = 1 in 2. 0 1 N 1 +in 2 N 3 Define ω := N 1 +in 2 and recall N i R. Thus we have N I 3 ω 2 A = ω N 3 N 3 ω. I ω N 2 3 Consider now 2 2 complex matrices B,C,D such that N B 0 B 0 I 3 ω 2 = ω N 3 C D C D N 3 ω. I ω N 2 3 This implies the following identities D D = I 2, C D = D N3 ω C =, ω N 3 B B +C C = I 2. N3 ω If we let D = I 2 and C = then we must solve ω N 3 B N3 ω N3 ω 1 0 B = I 2 = λ, ω N 3 ω N 3 0 1 where λ := 1 N 2 3 + ωω = 1 N2 1 +N2 2 +N2 3. Indeed λ = 1 r/t2 = s/t 2 0 so we can take B = λi 2 = s/ti 2. 3.2 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. 13

This provides another convenient way to study Dirac equations. Decompose the spinor ψ and source term F as u+v F+ +F ψ =, F =, 3.5 u v F + F 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 v +imu = if +, + u+imv = if. 3.6 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: E± σ s,u := u u± x j t u σ j u dx. 3.7 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: E σ ± s,w 1 2 Hs s 2 Furthermore for solutions u,v to 3.6 we have E σ + s,u+e σ s,v 1/2 E σ + s 0,u+E σ s 0,v 1/2 + s t 2w wdx 0. 3.8 s 0 F + H s+ F H sd s. 3.9 Proof. Step 1. Using the Dirac representation 2.3 and the decomposition 3.5, the PDE 2.1 becomes I2 0 0 I 2 u+v 0 + u v 0 σ j σ j 0 Defining ± := 0 ±σ i i this becomes + u+ v +im v + u u+v j +im u v u+v = i u v Adding and subtracting the two rows above gives the following v +imu = if +, + u+imv = if. u+v = i u v F+ +F F + F. 3.10 F+ +F. 3.11 F + F Following a similar approach to deriving 2.17, we multiply the first and second equation by v and u respectively. u 0 u+u σ j j u+imu v = iu F, v 0 v v σ j j v +imu v = iv F +. 14

One then adds these equations to their conjugate to obtain the following: 0 u u+ j u σ j u+imu v v u = iu F if u, 0 v v j v σ j v+imv u u v = iv F + if +v. Note the mass terms appear above. However if add these equations together we obtain 0 u u+v v+ j u σ j u v σ j v = iu F if u+iv F + if +v, 3.12 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 [s0,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: E σ ± s,w 1 2 Hs s 2 t 2w wdx 0. 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 w ± x j t σj w w ± x j t σj w dx = w w ±2 x j t w σ j w + x jx k t 2 w σ j σ k w dx = w w1+r/t 2 ±2 x j t w σ j w dx Thus we have E σ ± s,w = 1 2 = 2E σ ±s,w Hs s 2 t 2w wdx. w ± x j t σj w w± x j t σj w dx+ 1 s 2 Hs 2 t 2w wdx 0. Step 3. Next, let us show that the following hyperboloidal energy inequality holds for the Weyl spinor equation 3.6 E σ + s,u+e σ s,v 1/2 E σ + s 0,u+E σ s 0,v 1/2 + s Namely, integrating 3.12 over K [s0,s] we obtain E σ +s,u+e σ s,v = E σ +s 0,u+E σ s 0,v+ s s 0 F + H s + F H sd s. d s s/tiu F if u+iv F + if+vdx. s 0 H s Differentiating in s and noting that s/tu,s/tv u,v E σ + s,u + E σ s,v 1/2 gives d 1/2 E σ ds +s,u+e s,v σ F+ H τ + F H τ. 15

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 A ν m 2 A ν + ν µ A µ = ψ γ 0 γ ν P L ψ, iγ µ µ ψ = γ µ A µ P L ψ. 4.1 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 [22], 4.1 is equivalent to the following system A ν m 2 A ν = ψ γ 0 γ ν P L ψ, iγ µ µ ψ = γ µ A µ P L ψ, µ A µ = 0. 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: 4.2 A ν t 0, = a ν, t A ν t 0, = b ν, ψt 0, = ψ 0. 4.3 to satisfy the following two Lorenz compatibility conditions b 0 + j a j = 0, a 0 m 2 a 0 = ψ 0 P Lψ 0 j b j. For more generality, see also [3], we consider 4.1 with an artificial mass M added to the spinor: A ν m 2 A ν = ψ γ 0 γ ν P L ψ, iγ µ µ ψ +Mψ = γ µ 4.5 A µ P L ψ. 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 cr 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. 4.2 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 A 0, ψ 0. 4.6 4.4 16

Theorem 4.1. Consider the Dirac-Proca system 1.8 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.10, satisfying the smallness condition a ν,ψ 0 H N + b ν H N 1 ǫ, 4.7 the initial value problem of 1.8 admits a unique global-in-time solution ψ,a ν. Furthermore, the following decay result holds A ǫt 3/2, ψ ǫmin t 1,M 1 t 3/2. 4.8 5 Nonlinear stability of the ground state for the U1 model 5.1 The U1 model as a PDE system We now treat the Higgs mechanism applied to a U1 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 m 2 q Aν = Q A ν, χ m 2 q φ 0 2v 2 φ 0 χ χ φ 0 m 2 λ φ 0 2v 2 φ 0 χ+χ φ 0 = Qχ, iγ µ µ ψ m g ψ = Q ψ, 5.1 with quadratic or higher order terms given by Q A ν := iq χ ν χ ν χ χ +2q 2 A ν χ φ 0 +φ 0χ+χ χ +qψ γ 0 γ ν ψ, Q χ := 2iqA µ µ χ+q 2 χ φ 0 χ χ φ 0 +q 2 A µ A µ φ0 +χ +2λχ χφ 0 +2λχ φ 0χ+χ φ 0 +χ χ g φ 0 +χ ψ γ 0 ψ, Q ψ := g φ 0 χ+χ φ 0 +χ χ ψ qγ µ A µ ψ. 5.2 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λv2, 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 A ν,χ,ψ t 0, = A ν 0,χ 0,ψ 0, t A ν, t χ t 0, = A ν 1,χ 1, 5.4 which are said to be Lorenz compatible if they satisfy a A a 0 = A 0 1 iq φ 0χ 0 χ 0φ 0, A 0 0 m2 q A0 0 = ia i 1 iq φ 0 χ 1 χ 1 φ 0 +iq χ 0 χ 1 χ 1 χ 0 +2q 2 A 0 0 φ 0 χ 0 +χ 0 φ 0 +χ 0 χ 0 +qψ 0 ψ 0. 5.5 17

The derivation of 5.5 will follow from results of the following Section 5.2, in particular 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. 5.2 The abelian action and U1 invariance The Lagrangian we consider is L = 1 4 F µνf µν D µ φ D µ φ Vφ φ iψ γ 0 γ µ D µ ψ +gφ φψ γ 0 ψ, 5.6 where we use the following definitions for the Higgs potential, gauge curvature and gauge covariant derivatives: Vφ φ := λφ φ v 2 2, F µν := µ A ν ν A µ, D µ φ := µ iqa µ φ, D µ ψ := µ iqa µ ψ. 5.7 Furthermore λ, v, g, q are constants. A calculation shows that the Euler-Lagrange equations for 5.6 are the following µ F µν = iqφ D ν φ iqd ν φ φ+qψ γ 0 γ ν ψ, D µ D µ φ = V φ φφ gφψ γ 0 ψ, iγ µ D µ ψ = gφ φψ, together with 5.7. These PDEs can also be expressed as A ν ν diva 2q 2 A ν φ φ = iq φ ν φ ν φ φ +qψ γ 0 γ ν ψ, φ V φ φφ iqφ µ A µ = 2iqA µ µ φ+q 2 A µ A µ φ gφψ γ 0 ψ, iγ µ µ ψ gφ φψ = qγ µ A µ ψ. 5.8 5.9 Lemma 5.1. L has a U1-gauge symmetry. Proof. The U1 gauge symmetry induces the following transformations A µ A µ := A µ + µ α, φ φ := e iqα φ, ψ ψ := e iqα ψ, 5.10 where α = αt,x is some arbitrary function of space and time. This implies D µ φ D µ φ = e iqα D µ φ, D µ ψ D µ ψ = µ e iqα ψ iqa µ + µ αe iqα ψ = e iqα D µ ψ. 5.11 By the commutation of partial derivatives we immediately see that F µν is invariant under 18

U1 gauge transformations. Inserting these transformations into the Lagrangian one obtains L = iψ γ 0 γ µ D µψ 1 4 F µνf µν D µφ D µ φ Vφ φ +φ φ gψ γ 0 ψ = ie iqα ψ γ 0 γ µ µ e iqα ψ i e iqα ψ, iqγ µ A µ + µ αe iqα ψ 1 4 F µνf µν D µ φ e iqα e iqα D µ φ Vφ e iqα e iqα φ+gφ φe iqα ψ γ 0 e iqα ψ 5.12 = iψ γ 0 γ µ D µ ψ 1 4 F µνf µν D µ φ D µ φ Vφ φ+gφ φψ γ 0 ψ = L. Thus L is invariant under the transformation 5.10. 5.3 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 A ν 2q 2 A ν φ φ = ν X +iq φ ν φ ν φ φ +qψ γ 0 γ ν ψ, φ V φ φφ iqφx = 2iqA µ µ φ+q 2 A µ A µ φ gφψ γ 0 ψ, iγ µ µ ψ gφ φψ = qγ µ A µ ψ. Suppose the initial data for 5.13 satisfy diva X t0, = 0, t diva X t0, = 0. 5.13 5.14 Then as long as the solution to 5.13 exists, it will satisfy diva = X. Proof. To propagate the gauge choice diva = 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: diva X = iqφ φ φ φ+2q 2 ν A ν φ φ+q ν ψ γ 0 γ ν ψ = iqφ 2iqA µ µ φ+iqφ µ A µ +q 2 A µ A µ φ+v φ gφψ γ 0 ψ iq 2iqA µ µ φ iqφ µ A µ +q 2 A µ A µ φ +V φ gφ ψ γ 0 ψ φ +2q 2 ν A ν φ φ+qψ γ 0 γ µ µ ψ +q µ ψ γ 0 γ µ ψ = iq φ V φ V φ φ = 0, where in the last line we used the specific Higgs potential 5.7. 5.4 Gauge choice for the abelian model The Higgs potential has a non-zero minimum at φ 0 := ve iθ/v where Vφ 0 φ 0 = 0, 5.15 19

and θ : R 1+3 R is arbitrary. There is ambiguity in this minimum state due to the U1 symmetry represented by θ. A particular choice of θ will break this U1 symmetry, and such a scenario is termed spontaneous symmetry breaking. We consider perturbations of the form χ := φ φ 0, 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 diva+iqφ 0 χ χ φ 0 t 0, = 0, t diva+iqφ 0 χ χ φ 0 t 0, = 0. 5.16 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 use various techniques to treat the variables in 5.1 and 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 zero, m g = 0, then the nonlinearities in 6.1 will have insufficiently fast decay. In this case it is sufficient to employ a transformation introduced by [22], see 6.3 and 6.9 below, to overcome such problems. 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. 6.2 Hidden null structure from Tsutsumi Following Tsutsumi [22] we can now uncover a null structure using a particular transformation. However we will not need to reduce to the strong null forms via any complicated transformations [23, eq 2.6]. For complex-valued functions Φt,x,Ψt,x : R 3+1 C n, recall the null form Q 0 Φ,Ψ := 0 Φ 0 Ψ Φ Ψ. 6.2 Define a new variable à ν := A ν + q m 2 qψ γ 0 γ ν ψ. 6.3 This satisfies the non-linear Klein-Gordon equation m 2 q Ãν = QÃν, 6.4 20

in which the nonlinearities are QÃν := 2q m 2 0 ψ,γ qq 0 γ ν ψ+ q m 2 qg ψ γ0 γ ν ψ + q m 2 ψ γ 0 γ ν G ψ + 2q q m 2 m 2 g ψ γ 0 γ ν ψ, q G ψ := m g Q ψ iγ ν ν Q ψ = m g gφ 0 χ+χ φ 0 +χ χψ qγ µ A µ ψ igγ ν ν φ 0χ+χ φ 0 +χ χψ +iqγ ν γ µ ν A µ ψ gφ 0χ+χ φ 0 +χ χm g ψ +Q ψ 2igA µ µ ψ ga µ γ µ m g ψ +Q ψ, 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 to compensate for the lack of derivatives hitting ψ. In the case 0 m g minm 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 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. However, in the case m g minm q,m λ, better bounds are available and we do not require such a null structure any more. 6.3 Decomposition of χ In order to study the behaviour of χ, we find it more convenient to consider equations for the following two variables χ ± := φ 0 χ±χ φ 0. 6.5 Since the following identity holds χ + 2 + χ 2 = 2v 2 χ 2, 6.6 it is equivalent to estimate either χ or χ ±. These new variables satisfy the Klein-Gordon equations with Q χ+ := 2iqA µ µ χ +q 2 χ 2 +q2 A µ A µ 2v 2 +χ + χ + m 2 λ χ + = Q χ+, 6.7 χ m 2 qχ = Q χ, 6.8 +4λv 2 χ χ+2λχ 2 + +2λχ + χ χ g 2v 2 +χ + ψ γ 0 ψ, Q χ := 2iqA µ µ χ + +q 2 χ χ + +q 2 A µ A µ χ +2λχ χ + +2λχ χ χ gχ ψ γ 0 ψ. Similar to the previous section, we introduce the new variable which satisfies the following Klein-Gordon equation with the nonlinearity Q χ+ := 4 m g m 2 λ Q 0 ψ,γ 0 ψ+2 m g m 2 λ χ + := χ + 2m g m 2 ψ γ 0 ψ, 6.9 λ χ + m 2 λ χ + = Q χ+ 6.10 G ψ γ0 ψ +2 m g m 2 λ ψ γ 0 G ψ 4 m3 g m 2 ψ γ 0 ψ +Q χ+ +2m g ψ γ 0 ψ. λ 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. 21

7 Bootstrap argument 7.1 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 à ν = A ν + q m 2 q ψ γ 0 γ ν ψ, χ + = χ + 2m g m 2 ψ γ 0 ψ, χ, ψ, 7.1 λ which satisfy equations Ãν m 2 qãν = QÃν, χ + m 2 λ χ + = Q χ+ χ m 2 qχ = Q χ, iγ µ µ ψ m g ψ = Q ψ. We also provide here the following theorem, which details the and L estimates of the unknowns. Theorem 7.1. Under the smallness assumptions in Theorem 1.1 the solution satisfies the following energy estimates s/t I L J µ A ν,χ f + s/t µ I L J A ν,χ f C 1ǫ, I + J N, I L J A ν,χ f C 1ǫ, I + J N, s/t I L J ψ f C 1ǫ, I + J N, s/t α I L J ψ f C 1ǫlogs, I + J = N, and the following L estimates sup t,x t 1/2 s α I L J A ν,χ +t 1/2 s I L J α A ν,χ C 1 ǫ, I + J N 2, sup t 3/2 I L J A ν, I L J χ C 1 ǫ, I + J N 2, t,x sup t 1/2 s I L J ψ C 1 ǫ, I + J N 2, t,x sup t,x where C 1 is some constant to be determined later. 7.2 t 1/2 s α I L J ψ C 1 ǫlogs, I + J = N 2, 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 the estimates for null terms and various commutators. The bootstrap assumptions will be made in 7.10. These bootstraps, combined with some standard commutator estimates and Sobolev-type inequalities, will lead to certain weak estimates in 7.12 and 7.14. In Section 7.4 we will use our first-order hyperboloidal energy to 7.3 22

upgrade our estimates for the Dirac component, namely in Theorem 7.7 and Corollary 7.8. 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. 7.2 Estimates for null forms and commutators 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. I L J T µν µ u ν v I 1 L J 1 a u I 2 L J 2 β v + I 1 L J 1 β u I 2 L J 2 a v I 1 + I 2 I, J 1 + J 2 J, a,β +s/t 2 I 1 + I 2 I, J 1 + J 2 J I 1 L J 1 t u I2 L J 2 t v. 7.4 One refers to [16] for the proof. We next recall the estimates for commutators, also proved in [16]. Heuristically speaking, the following lemma provides the relation Γ I L J u I L J Γu, in which Γ takes values from the set { α, α, α β,s/t α }. 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 β I L J u, 7.5 [ I L J, a ]u C I, J [ I L J, α ]u C I, J [ I L J, α β ]u C I, J J < J,β I < I, J < J,b I < I, J < J,β I I, J < J,γ,γ I L J s/t α u s/t α I L J u +C I, J b I L J u +t 1 β I L J u +t 1 I I, J J I I, J J,β I L J u, 7.6 β I L J u, 7.7 γ γ I L J u, 7.8 I I, J J,β s/t β I L J u. 7.9 23

7.3 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 mq s, I L J A ν 1/2 +Emλ s, I L J χ 1/2 C1 ǫ, I + J N, If we can prove the refined estimates E mg s, I L J ψ 1/2 C 1 ǫ, I + J N 1, E mg s, I L J ψ 1/2 C 1 ǫlogs, I + J = N. E mq s, I L J A ν 1/2 +Emλ s, I L J χ 1/2 1 2 C 1ǫ, I + J N, E mg s, I L J ψ 1/2 1 2 C 1ǫ, I + J N 1, 7.10 7.11 E mg s, I L J ψ 1/2 1 2 C 1ǫlogs, I + J = N, then we are able to assert that s 1 cannot be finite, which in turn implies a global existence result for 1.1. Combining bootstrap assumptions 7.10 and the estimates for commutators in Lemma 7.3, the following sets of estimates are immediately obtained: s/t I L J µ A ν,χ f + s/t µ I L J A ν,χ f C 1ǫ, I + J N. I L J A ν,χ f C 1ǫ, I + J N, m g I L J ψ f + s/t µ I L J ψ f C 1ǫ, I + J N 1, m g I L J ψ f + s/t µ I L J ψ f C 1ǫlogs, I + J = N. 7.12 FurthermorewecanobtainL estimatesbyrecallingthefollowingsobolev-typeinequality on hyperboloids [16]. Proposition 7.4 Sobolev-type inequality on hyperboloids. For all sufficiently smooth functions u = ut,x supported in the region {t,x : x < t 1}, then for s 2 one has sup t 3/2 ut,x L J u H f H, 7.13 s s J 2 where the summation is over Lorentz boosts L. Combined with the estimates 7.12 the following also hold: t 1/2 s α I L J A ν,χ +t 1/2 s I L J α A ν,χ C 1 ǫ, I + J N 2, sup t,x sup t 3/2 I L J A ν, I L J χ C 1 ǫ, I + J N 2, t,x sup m g t 3/2 I L J ψ +t 1/2 s α I L J ψ, I L J α ψ C 1 ǫ, I + J = N 3, t,x sup t,x m g t 3/2 I L J ψ +t 1/2 s I L J ψ, I L J α ψ C 1 ǫlogs, I + J N 2. 24 7.14

7.4 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 standardtechniques to estimate the nonlinearity; see for example[2] and[17]. 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. We first adapt our Sobolev-type inequality on hyperboloids 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 L a := L a + 1 2 γ 0γ a. 7.15 It then holds that [ L a,iγ ν ν ] = 0. The following result is a simple extension of Proposition 7.4. Corollary 7.5. Suppose ψ = ψt, x is a sufficiently smooth spinor field supported in the region {t,x : x < t 1}, then it holds for s 2 sup t 3/2 ψt,x LJ ψ H f H, 7.16 s s J 2 where L denotes a modified Lorentz boost. We will also make use of the following Gronwall inequality. Lemma 7.6 Gronwall-type inequality. Let ut be a non-negative function that satisfies the integral inequality t ut C + bsus 1/2 ds, C 0, 7.17 t 0 where bt is non-negative function for t t 0. Then it holds ψ. t 2 ut C + bsds. 7.18 t 0 We are now in a position to obtain certain improved energy and sup norm estimates for Theorem 7.7. Under the same assumptions as Theorem 7.1, the Dirac field ψ satisfies s/t LJ ψ f ǫ+c 1ǫ 2, J N, 7.19 sup t 1/2 s L J ψ ǫ+c1 ǫ 2, J N 2. 7.20 As a consequence, one has the following sup-norm estimates for ψ: sup t I LJ ψ ǫ+c 1 ǫ 2, I + J N 2. 7.21 25

Proof. Step 1. Recall the equation in 5.1 for the Dirac field iγ µ µ ψ m g ψ = Hψ, H := gφ 0χ+χ φ 0 +χ χ gγ µ A µ. Since iψ γ 0 H ih γ 0 ψ = 0 we have the following conserved energy E,ψ = E 0,ψ, 7.22 then according to the inequality 2.24c in Proposition 2.3, we are able to initialise the induction argument by s/tψ f ǫ. Step 2. For induction purposes, assume s/t L J ψ f ǫ+cǫ2 holds for 0 J k 1 N 3, and now consider the case 1 J = k N 2. Act L J on the Dirac equation above to obtain with γ µ µ L J ψ+im g LJ ψ = ih L J ψ ir, R := L J Hψ H L J ψ, H = gφ 0 χ+χ φ 0 +χ χ gγ µ A µ. Observe that R contains only terms, up to some constant matrices, of type L J 1 H L J 2 ψ, J 1 + J 2 J, J 2 k 1. 7.23 Using 2.25 and consequently by Lemma 7.4 and the induction assumption s E, L J ψ 1/2 E 0, L J ψ 1/2 + R f H s d s s 0 s ǫ+ t/s L J 1 H L s/t L J 2 ψ H s d s f H s s 0 J 1 + J 2 J J 2 J 1 ǫ+c 1 ǫ 2 s s 0 s 3/2 d s, 7.24 which gives s/t LJ ψ ǫ, J k. f Hs Step 3. The above analysis shows for J N 4 LJ s/t LJ ψ f H ǫ+c s 1ǫ 2. J 2 Thus by the Sobolev inequality 7.13, we deduce sup t 1/2 s L J ψ ǫ+c 1 ǫ 2, J N 4. 26

Step 4. We now consider the case J = N 1. An energy estimate yields E, L J ψ 1/2 ǫ+ + s s 0 J 1 + J 2 N 1, J 1 > J 2, J 2 N 4 J 1 + J 2 N 1, J 1 J 2, J 2 N 2 ǫ+c 1 ǫ 2 s s 0 s 3/2 d s, LJ 1 H LJ 2 ψ f H s L H s t/s LJ 1 H s/t LJ 2 L ψ H s f H s d s 7.25 which implies s/t L J ψ f H s ǫ+c 1ǫ 2. 7.26 The same analysis also applies to the case J = N. And repeating Step 3 gives 7.20 for J N 2. As a consequence, we have the following sup-norm estimates for ψ. Corollary 7.8. It holds that 7.5 Refined estimates sup t 1/2 s I L J ψ ǫ+c1 ǫ 2, I + J N 2. 7.27 In this final subsection we close our bootstrap argument. For this to work we move to the transformed vector field Ãν defined in 6.3 and the transformed scalar field χ + defined in 6.9, which are heuristically of the form A ν = Ã ν +O ψ 2, χ + = χ + +O ψ 2. Thus using the estimates for A ν and χ + coming from 7.12 and 7.14, together with the previous energy and sup-norm estimates for ψ, the following estimates for Ãν and χ + hold s/t I L J µ Ãν, χ + f H + s s/t µ I L J Ãν, χ + f H C s 1ǫ, I + J N, I L J Ãν, χ + f H C s 1ǫ, I + J N, sup t 1/2 s α I L J Ãν, χ + +t 1/2 s I L J α Ãν, χ + C 1 ǫ, I + J N 2, sup t 3/2 I L J Ãν, χ + C 1 ǫ, I + J N 2. We first look at the energy for ψ in the case I + J = N E mg s, I L J ψ 1/2 ǫ+c ν,µ s s 0 I L J µ Hψ f H s + I L J H µ ψ f H s 7.28 d s. 27

By noting I L J µ Hψ f H s µ I 1 +I 2 =I,J 1 +J 2 =J I 1 + J 1 N 2,µ + I 1 +I 2 =I,J 1 +J 2 =J I 1 + J 1 N 1,µ t/s I 1 L J 1 µ H L s/t I 2 L J 2 ψ H s f H s s/t I 1 L J 1 µ H f H s t/s I 2 L J 2 ψ L H s C 1 ǫ 2 s 3/2 log s+c 1 ǫ 2 s 1 C 1 ǫ 2 s 1, 7.29 we obtain E mg s, I L J ψ 1/2 ǫ+cc 1 ǫ 2 logs. 7.30 Similarly for I + J N 1 we obtain E mg s, I L J ψ 1/2 ǫ+cc 1 ǫ 2. 7.31 In order to obtain estimates for A ν, we first bound the energy for Ãν s E m s, I L J Ã ν 1/2 ǫ+ I L J γ 0 γ ν µ ψ, µ ψ s f H s +m2 g I L J ψ γ 0 γ ν ψ f H s 0 + I L J ψ 2 A+ψA ψ +ψ 3 ψ f H s d s ǫ+cc 1 ǫ 2. 7.32 Next, recalling definition 6.3 we use Young s inequality to obtain for all I + J N E mq s, I L J A ν 1/2 3/2E mq s, I L J Ã ν 1/2 +3/2E mq s, I L J ψ 2 1/2 3/2ǫ+C I 1 L J 1 ψ L H Es, I 2 s L J 2 ψ 1/2 ǫ+c 1 ǫ 2. I 1 +I 2 =I,J 1 +J 2 =J I 1 +J 1 N 2 7.33 A similar procedure gives the refined estimates for χ + while the refined estimates for χ E mλ s, I L J χ + 1/2 3/2ǫ+CC 1 ǫ 2, 7.34 E mq s, I L J χ 1/2 ǫ+cc 1 ǫ 2 7.35 can be obtained directly. A combination of 7.34 and 7.35 gives the refined estimates for χ E mλ s, I L J χ 1/2 ǫ+c 1 ǫ 2. 7.36 By choosing C 1 sufficiently large and ǫ sufficiently small, we arrive at the refined bounds 7.11. This shows global existence and thus completes the proof of Theorem 7.1. 28