Massive gauge particles versus Goldstone bosons in non-Hermitian non-Abelian gauge theory

We investigate the Englert-Brout-Higgs-Guralnik-Hagen-Kibble mechanism for non-Hermitian field theories with local non-Abelian gauge symmetry in different regions of their parameter spaces. We demonstrate that the two aspects of the mechanism, that is giving mass to gauge vector bosons and at the same time preventing the existence of massless Goldstone bosons, remain to be synchronized in all regimes characterized by a modified CPT symmetry. In the domain of parameter space where the"would be Goldstone bosons"can be identified the gauge vector bosons become massive and the Goldstone bosons cease to exist. The mechanism is also in tact at the standard exceptional points. However, at the zero exceptional points, that is when the eigenvalues of the mass squared matrix vanish irrespective of the symmetry breaking, the mechanism breaks down as the Goldstone bosons can not be identified and the gauge vector bosons remain massless. This breakdown coincides with the vanishing of the CPT inner product of symmetry breaking vacua defined on the eigenvector space of mass squared matrix. We verify this behaviour for a theory with SU(N) symmetry in which the complex scalar fields are taken in the fundamental as well as in the adjoint representation.


Introduction
Our main objective is to extend the Englert-Brout-Higgs-Guralnik-Hagen-Kibble mechanism [1,2,3,4], hereafter simply referred to as Higgs mechanism, to non-Hermitian field theories with a local non-Abelian gauge symmetry using a pseudo-Hermitian approach. We focus on the two key aspects for which the mechanism was originally developed, that is to give mass to gauge vector bosons and at the same time prevent the existence of massless Goldstone bosons. When keeping the symmetry global one may adopt different starting points for the study of Goldstone phases, such as the field content of local operators, a scattering matrix based on a particle picture or an explicit Lagrangian.
For instance, two dimensional conformal quantum field theories are well understood in terms of their operator content characterized by infinite-dimensional algebras of local conformal transformations [5]. A large class of such theories, minimal models [6], are known to possess a finite operator content and the treatment of unitary and non-unitary theories is formally identical. The simplest massive non-unitary field theory consisting of only one real scalar field describing in its ultraviolet limit the critical point of the Ising model in a purely imaginary magnetic field, the Yang-Lee edge singularity [7,8], is known for a long time to correspond to the non-Hermitian Lagrangian [9,10] (1.1) Exact scattering theories for two-dimensional models have also been identified [11], that can be used to probe the ultraviolet limit most easily by employing the thermodynamic Bethe Ansatz [12]. These techniques have also been employed for hypothetical scattering matrices for massless Goldstone fermions (Goldstinos) [13] and scattering matrices that reduce to them in certain limits [14]. Despite the fact that the Mermin-Wagner theorem prevents the validity of the Goldstone theorem in dimensions d ≤ 2, it was argued in [15] that for certain symmetry groups, e.g. SO(N ) with formally taking N < 2, this restriction can be circumvented so that Goldstone phases maybe be identified in such type of non-Hermitian systems.
Rather than taking an operator content or a scattering matrix as a starting point, one may of course also commence directly with a non-Hermitian Lagrangian. From that perspective it is natural to try to extend techniques and methods developed for the treatment of non-Hermitian quantum mechanics [16,17,18,19] to a quantum field theory setting. Such considerations have been carried out for a scalar field theory with imaginary cubic self-interaction terms [20,21], with a Lagrangian identical to (1.1) but for h = h c without a linear term, deformed harmonic oscillators [22], non-Hermitian versions with a field theoretic Yukawa interaction [23,24,25,26], free fermion theories with a γ 5 -mass term and the massive Thirring model [27], PT -symmetric versions of quantum electrodynamics [28,29,30] and other types of PT -symmetric quantum field theories in higher dimensions [22] than (1.1).
Non-Hermitian quantum field theories are also viewed as possible options to overcome the limitations of the Standard Model of particle physics, which is well known to fail in describing gravity and neutrino oscillations in a consistent manner. In [31] a study was initiated to potentially resolving the latter issue by extending the ordinary two-flavor neutrino oscillation to a non-Hermitian PT-symmetric/pseudo Hermitian setting. A well studies simple Hermitian extension of the Standard Model referred to as the two-Higgs-doubletmodel [32] involves a second Higgs doublet or possibly more [33]. Here we investigate non-Hermitian variants of these models.
Here we are especially interested in complex non-Hermitian scalar field theories and the question of how the aforementioned Goldstone phases manifest in these theories, together with the subsequent extension to the Higgs mechanism [1,2,3,4] in Abelian and non-Abelian gauge theories. These issues have been studied recently by various groups in different approaches, which differ from their very onset: Given a generic action for a complex scalar field theory of the form I = d 4 xL(φ, φ * ), one has two options in a Hermitian theory to derive the equations of motion by means functional variation, either to calculate δI/δφ = 0 or δI/δφ * = 0. Since the standard CPT -theorem [34] applies, the two resulting equations are the same. In contrast, in a non-Hermitian theory one no longer has I = I * , so that the two equations are not only not the same, but in addition one also has the new options δI * /δφ = 0 and δI * /δφ * = 0. In the first approach, we refer to as the "surface term approach", it was suggested [35,36] to take of only the two equations resulting from δI/δφ = 0, δI * /δφ * = 0 and neglect the remaining two. As the resulting equations are in general not compatible, the authors propose to use some non-vanishing surface terms to compensate for the discrepancy. The second approach [37] consists of taking δI/δφ = 0 or δI/δφ * = 0, and when determining the vacuum allowing the real and imaginary parts of the complex scalar field to be also complex. Thus in this approach the field content is doubled or re-defined.
Here we follow an approach, we refer to as the "pseudo-Hermitian approach" [38,39], more aligned to the procedure pursued in non-Hermitian versions of quantum mechanics, in which one employs so-called Dyson maps [40] to transform a non-Hermitian Hamiltonian to a Hermitian Hamiltonian. Since the action I contains a Lagrangian, rather than a Hamiltonian, we need to first Legendre transform the complex Lagrangian L to a non-Hermitian Hamiltonian H, carry out the similarity transformation by means of a Dyson map, while preserving equal time commutation relations, to obtain a Hermitian Hamiltonian h, which we then inverse Legendre transform to a real Lagrangian l A consistent set of equations of motion is then obtained by functionally varying the action s = d 4 xl(ϕ, χ) involving this real Lagrangian l with respect to the real field components ϕ, χ ∈ R of the complex scalar field φ = 1/ √ 2(ϕ + iχ), i.e. δs/δϕ = 0 and δs/δχ = 0. In order to perform the similarity transformation, one needs to canonically quantize the theory first, which is what we will do below.
In fact, this will be the only quantum aspect in our manuscript. Our discussion here is kept classical and we will not carry out a second quantization of the scalar fields, nor BRST quantization for the gauge fields, etc. Supported by our analysis below and in our previous papers [38,39], our view here is that much insight can be gained by a classical treatment, apart from the appeal to the quantum field theory for the canonical commutation relations. One may of course fully analyse the theory in its quantum aspects, but our view is that this should be performed on the transformed, i.e. Hermitian, theory rather than directly on the non-Hermitian theory. In the former scenario one would not expect any issues to arise as one is simply dealing with a Hermitian system, whereas in the latter case one may encounter a number of what we would refer to as "pseudo-problems".
The above mentioned incompatibility of the variational principle is an example for such a problem that already occurs on the level of the classical theory. One may try to fix this problem directly for the non-Hermitian system or simply consider the equivalent Hermitian system in which the problem is entirely absent. Another well known example for a pseudo-problem in the quantum mechanical context is for instance to take the variables x and p, that might appear in the definition of a non-Hermitian Hamiltonian, to be physical observables. However, only the pseudo-Hermitian counterparts or the Dyson mapped quantities ηxη −1 and ηpη −1 can be observed and should be interpreted as being physical. Taking instead the original operators x and p from the non-Hermitian theory one may derive violations of the uncertainly relations and other properties that contradict standard quantum mechanical principles. Again these problems entirely disappear in the equivalent Hermitian system when considering the correct variables.
Previously we have analysed the Goldstone theorem for a non-Hermitian scalar field theory with an Abelian [38] and a non-Abelian symmetry [39]. Here the main purpose is to investigate the extension to the Higgs mechanism. As it is well known, in the standard Higgs mechanism the Goldstone bosons acquires a mass so that our previous findings will inevitably have a bearing on the investigation to be carried out here. Let us therefore recall the key finding from [38,39]: The main object of study has been the eigenvalue spectra of the non-Hermitian squared mass matrix M 2 , obtained by expanding around the symmetry breaking or preserving vacua. The reality of these eigenvalues has been guaranteed by a modified CPT -symmetry of the original action I. Hence, we distinguished in the usual fashion between CPT -symmetry regime characterized by M 2 commuting with this symmetry operator and its eigenstates being simultaneous eigenstates of the symmetry operator. When the latter is not the case, one refers to that regime as the CPT -spontaneously broken regime in which some eigenvalues become complex conjugate pairs. The points in parameter space at which this occurs are commonly referred to as exceptional point. As physical masses are positive and real, we also require the eigenvalues of M 2 to be non-negative. We encountered a special behaviour at the transition points when the eigenvalues become zero, which we referred to [38,39] as zero exceptional points of type I and type II. 1 At the type I points the mass matrix is non-diagonalizable and the continuous symmetry is broken, whereas at the type II points the mass matrix can be diagonalized and the vacuum with broken continuous symmetry re-acquires the symmetry at this point.
Using the above mentioned approaches, the Higgs mechanism was previously studied for Abelian [37,41] as well as non-Abelian gauge theories [36] leading to slightly different findings. In [37] the interesting observation made, that the mass of the gauge vector boson vanishes at the zero exceptional point, was not confirmed in [41]. In addition, for the non-Abelian gauge theories it was found in [36] that the Higgs mechanism even applies in the spontaneously broken CPT -regime. Our aim is here to compare the various observations made using these alternative approaches with a pseudo-Hermitian approach, extend the studies to other models, symmetries and representations within this framework.
Our manuscript is organised as follows: In section 2.1 we introduce first our non-Hermitian model with scalar field taken in the fundamental representation possessing a local SU (2) × U (1)-symmetry, discuss the symmetry breaking vacuum of this theory and subsequently the Higgs mechanism. We indicate how to extend the model from SU (2) to SU (N ). We repeat the same discussion in section 2.2 for anon-Hermitian model with scalar field taken in the adjoint representation. Our conclusions are stated in section 3. The distinction between the two types of exceptional points is crucial and for that reason we include in our appendix a discussion for part of the squared mass matrix that explains the difference.

Pseudo-Hermitian approach to the Higgs mechanism
In this section we commence by investigating the same model considered in [36] using, however, a pseudo-Hermitian method to compare our results with the findings in [36]. We will observe that the mass spectrum of the fields in the SU (2) fundamental representation coincides with the one found in [36], but the masses for the gauge vector bosons differ and in particular vanish at the zero-exceptional points. We will extend this model to incorporate a SU (N )-symmetry and continue to observe this phenomena also for these more general systems. Finally we will consider a new model for which the fields are taken in a different representation, the adjoint representation of SU (2), making similar observations.

A SU (2)-model in the fundamental representation
We start by applying the pseudo-Hermitian approach to a model with local SU (2) × U (1)symmetry previously studied using the surface term approach in [36]. The model corresponds to the gauged version of the one for which the Goldstone mechanism was studied in [38] Here g, µ ∈ R, m i ∈ R or m i ∈ iR are constants. When compared to [38] we have replaced here as usual the standard derivatives ∂ µ by covariant derivatives D µ := ∂ µ −ieA µ , involving a charge e ∈ R and the Lie algebra valued gauge fields A µ := τ a A a µ . Here the τ a , a = 1, 2, 3, are taken to be Pauli matrices, which when re-defined as i(−1) a+1 τ a are the generators of SU (2). We have also added the standard Yang-Mills term comprised of the Lie algebra valued field strength The two complex scalar fields φ i are taken to be in the representation space of fundamental representation of SU (2). The model described by L 2 admits a global continuous U (1)-symmetry, a local continuous SU (2)-symmetry and two discrete antilinear CPT -symmetries as described in more detail in [38]. Crucially L 2 is not Hermitian, which at this point is simply to be understood as not being invariant under complex conjugation. The Abelian version of L 2 was discussed in [37,39].
As argued in [38], it is useful to decompose the complex fields into their real components Thus simply rewriting the complex scalar fields in equation (2.1) in terms of their real and imaginary components we obtain the following We use here the standard notation * for complex conjugation and † for the simultaneous conjugation with transposition.
The SU (2)-symmetry manifests itself as follows: A change in the complex scalar fields due to this symmetry is δφ k j = iα a T kl a φ l j , where the generators T a of the symmetry transformation are the standard Pauli matrices σ a , a = 1, 2, 3. The infinitessimal changes for the real component fields are then identified as which leave the above Lagrangian invariant. The discrete antilinear CPT ± -symmetries manifest themselves as A noteworthy remark is that it is straightforward to generalise the model from a locally SU (2)-invariant one to a locally SU (N )-invariant one, by extending the sum over k from 2 to N , while keeping the U (1)-symmetry global. In what follow we will focus on N = 2.
A crucial feature of L 2 is that its CPT -invariance translates into pseudo-Hermiticity [42,18], meaning that it can be mapped to a Hermitian Lagrangian l 2 by means of the adjoint action of a Dyson map η as l 2 = ηL 2 η −1 . This may be achieved by the slightly modified version of the Dyson map used in [37,38] We denote here the time-dependence by t ′ to indicate that commutators are understood as equal time commutators for the canonical momenta is not to be viewed as explicitly time-dependent as discussed in much detail for instance in [43]. The adjoint action of η + 2 on the individual fields maps as Thus, we convert the complex Lagrangian into the real Lagrangian We may transform here directly the Lagrangian rather than the Hamiltonian, as suggested in (1.2), since the kinetic energy term is real and the complexity only result from the potential term. Introducing the 2 two-component fields of the form we can re-write the Lagrangians L 2 and l 2 more compactly. Defining the 2 × 2 matrices the real Lagrangian l 2 acquires the form We have simplified here the index notation by implicitly contracting, keeping in mind that we are summing over two separate index sets k ∈ {1, 2} and j ∈ {1, 2}. For instance, we set In this formulation we may think of the real and complex Lagrangians, l 2 and L 2 , as being simply related by a kind of Wick rotation in the field-configuration space One may of course wonder about the negative sign in the kinetic term of l 2 and whether these lead to ghost fields. Thus we finish this subsection with a short discussion that establishes that these signs disappear when the Lagrangian is properly diagonalised. As this argument is the same in the global and local theory, we set the gauge fields to zero without loss of generality for this purpose and consider the corresponding action where S int contains all terms of higher order than Φ 2 . We also assumed that surface terms vanish at infinity when integrating by parts and used I 2 = I, M 2 = IH ± . Next we diagonalise the squared mass matrix as M 2 = T −1 DT and consider next only the integrant of the first term in (2.17) where we introduced the new field Ψ := T ⊤ IΦ and used T −1 = T ⊤ I. The latter relation is derived as follows: We start by defining the right and left eigenvectors of M 2 as respectively. Noting that I † = I and H † ± = H ± the last relation implies that IH ± Iu i = λ i Iu i . Therefore we can express the right eigenvectors in terms of the left eigenvectors as Since the left and right eigenvector form a biorthonormal basis, v i ·u j = δ ij , it follows that T ⊤ IT = I and hence T −1 = T ⊤ I.

The symmetry breaking vacuum
The vacuum solutions Φ k 0 , Ψ k 0 by solving δV = 0, with V denoting the potential amounts to solving the two equations Hence in the real component field configuration space the vacuum manifold is a S 3 -sphere with radius R. Consequently, we may consider the equations (2.21) as two eigenvalue equations. Thus, besides the trivial SU (2)-invariant vacuum Φ k 0 = Ψ k 0 = 0, k = 1, 2, we must have zero eigenvalues in both equations, which is equivalent to requiring Since R 2 is positive, this equality imposes restrictions on the parameters g, µ and the possible choices for m 1 ∈ R, m 2 ∈ iR or m 1 ∈ iR, m 2 ∈ R. The corresponding vectors that satisfy equation) (2.21), suitably normalized with regard to the standard inner product, are Imposing now the constraint on R 2 as stated after equation (2.21), a possible solution Hence we recover the symmetry breaking vacuum used in [39].

The Higgs mechanism
Let us now demonstrate how the gauge vector boson acquires a finite mass and how at the same time the emergence of a Goldstone boson is prevented by the Higgs mechanism [1,2,3,4] in the CPT -symmetric regime, at the exceptional points and even in the spontaneously broken CPT -symmetric regime. The mechanism breaks down at the two types of zero exceptional points.
Expanding the potential V in (2.20) around the vacuum specified at the end of subsection 2.1.1 leads to As expected, multiplying the Hessians in (2.25) by I gives back the squared mass matrix we found in [39]. The kinetic term is almost unchanged except for the term involving Ψ 2 The last term corresponds to the mass term of the gauge vector boson that we evaluate to where we used the standard relation τ a † τ b = τ a τ b = δ ab I + iε abc τ c . Therefore we read off the mass of each of the three components of the gauge vector boson as In [39] we identified the physical regions in the parameter space in which the squared mass matrix has non-negative eigenvalues and in which the Goldstone bosons can be identified.
Let us now compare those regions with the values for which the gauge vector boson becomes massive. We immediately see from the expression in (2.29) that the gauge vector boson remains massless when µ 4 = m 4 2 or when R = 0, i.e. µ 4 = −m 2 1 m 2 2 . The two sets of values correspond precisely to the two types of zero exceptional points, type I and II, respectively, at which the squared mass matrix develops zero eigenvalues λ = 0. These points are distinct from standard exceptional points where two eigenvalues coalesce and become complex thereafter, here at λ = µ 4 m 2 2 − m 2 2 . See the appendix for a more detailed explanation about the distinction between these types of exceptional points.
Thus the two aspects of the Higgs-mechanism, i.e. giving mass to the gauge vector boson and at the same time preventing the existence of the Goldstone bosons, remain to go hand in hand. In the CPT -symmetric regime the mechanism applies, but at the zero exceptional points the Higgs-mechanism breaks down as the Goldstone bosons are not identifiable [39] and at the same time the gauge vector boson remains massless. In contrast, at the exceptional point the Goldstone bosons are identifiable [39], although in a different manner, and the gauge vector bosons become massive.
Let us see this in detail by following [39] and replacing m 2 i → c i m 2 i , with c i = ±1 to account for all possibilities in signs. We found that physical regions only exist for the two cases c 1 = − c 2 = 1 and c 1 = − c 2 = −1. For the two cases we may then write noting that m 2 g /m 2 1 only depends on the two parameters m 2 2 /m 2 1 and m 4 2 /m 4 1 similarly as the eigenspectrum of the squared mass matrix [36,39]. We require the right hand side of (2.30) to be positive as depicted in depict in figure 1.
We observe in figure 1 that while the region in which the Goldstone boson can be identified is bounded by exceptional as well as zero exceptional points, the exceptional points lie well inside the region for which the gauge vector boson is massive, i.e. they acquire a mass in the CPT -symmetric regime as well as in the spontaneously broken CPTsymmetric regime. In the CPT -symmetric regime this agrees well with the findings that at these points the "would be Goldstone boson" is prevented from existing as a massless particle. We may think of the sign change in front of the mass terms, c i → −c i , that relates the left to the right panel as a phase transition [44].
Let us now demonstrate this behaviour in detail and expand for this purpose the Lagrangian around the symmetry broken vacuum up to second order in the fields We recall now from [39] that the first two lines of the Lagrangian l 2 can be diagonalized and the Goldstone bosons can be identified in terms of the field content of the model. Furthermore, the Goldstone modes are eigenvectors in the null space of squared mass matrices computed above as Ψ 2 0 and IΨ 2 0 , so that the Goldstone modes are proportional to these two vectors. The explicit forms of the Goldstone fields were found in [38] (equation (3.40) therein), denoted as ψ Gb 5 , ψ Gb 3 and ψ Gb 1 , therein. We express them here as respectively. As expected for the Higgs mechanism the number of "would be Goldstone bosons" equals the amount of massive vector gauge bosons. The fact that the Goldstone modes are inverse proportional to the mass of the gauge bosons explains that they can not be identified for massless gauge bosons. Keeping now only the Goldstone kinetic term from the first two lines of the Lagrangian l 2 and the one involving the gauge fields in equation (2.31), we obtain (2.34) Using the explicit representations of the Pauli matrices, the real and imaginary parts are determined as Finally the Lagrangian in (2.34) can be simplified to where we defined the new vector gauge particle with component fields B a µ := A a µ − 1 mg ∂ µ G a . We may also replace A a µ by B a µ in the field strength F µν so that A µ can be eliminated entirely from the Lagrangian. We see that the Higgs-mechanism applies as long as m g = 0. However, at the zero exceptional points, not only the gauge boson mass vanishes, but the Higgs mechanism no longer applies, in the sense that we can not remove the degrees of freedom of Goldstone bosons.
We summarize the behaviour we found in the different types of regimes in the following Thus we encounter three different types of behaviour: In the CPT -symmetric regime, at the standard exceptional points as well as in the spontaneously broken CPT -symmetric regime the Higgs mechanism applies in the usual way. However, in the latter regime other particles in the theory become non-physical. At the zero exceptional points the vector gauge bosons remain massless and no Goldstone bosons can be identified in the global theory.

From SU (2) to SU (N )
We will now follow the same line of reasoning as in the previous subsection and generalize our model from possessing a SU (2)-symmetry to one with a SU (N )-symmetry. For this purpose we simply replace the Pauli matrices in all our expressions by the traceless and skew-Hermitian N × N -matrices corresponding to the SU (N )-generators T a with a = 1, . . . , (N 2 − 1). The vacua are still determined by the solutions of the eigenvalue problem (2.21) with zero eigenvalue condition The zero eigenvalue condition implies that the vacuum manifold is a S 2N −1 -sphere with radius R. This follows from the fact that SU (N ) acts on the 2N dimensional space spanned by (ϕ 0 1 ) i , (χ 0 1 ) i , i = 1, . . . , N , with norm equal to R 2 . On this space SU (N − 1) simply permutes the fields amongst themselves, hence acting as a stabilizer or isotropy subgroup. Thus the vacuum manifold corresponds to the coset SU (N )/SU (N − 1) ∼ = S 2N −1 .
As we discussed in detail in [39], we may utilize the symmetry of the Lagrangian to transform the vacua into convenient forms without changing the eigenvalue spectrum of the mass matrix. Thus using the elements T ∈ SU (N )/SU (N − 1) ⊂ SU (N ) we may transform the vacuum into the form satisfying the constraint (2.38). Let us now use this SU (N )-symmetry breaking vacuum to calculate the mass of the gauge vector boson. Taking the proper SU (N )-algebra rather than the physicist's version, as in the last subsection for SU (2), we also change e → ie.
Dropping here the kinetic term reported in (2.26) and considering only the relevant term in the Lagrangian we obtain We evaluate the last factor using the identity where the g abc and f abc are completely symmetric and anti-symmetric tensors, respectively. We note that from which we read off the masses m (a) g of the N 2 − 1 gauge vector bosons. We note that once again they vanish at the zero exceptional points, but now for all SU (N )-models.

A SU (2)-symmetric model in adjoint representation
As we have demonstrated, the gauge vector boson becomes massive for the SU (N )-symmetric model in the CPT -symmetric regime and at the exceptional point when the fields are taken to be in the representation space of the fundamental representation. On the other hand the Higgs-mechanism breaks down at the zero exceptional points. Remarkably it still applies when the CPT -symmetry is broken, although in that regime other particles acquire complex masses so that the region is non-physical. Recall that the regions in which the CPT -symmetry holds are identified in figure 1.
Let us now see whether we encounter a similar behaviour when the fields are taken in adjoint representation. We consider here a slightly different non-Hermitian SU (2)-invariant Lagrangian where as in equation (2.2) we take g, µ ∈ R, m i ∈ R or m i ∈ iR, to be constants. The two complex scalar fields are expressed as φ i = φ a i T a , i = 1, 2 and a = 1, 2, 3, where the T a are the three SU (2)-generators in the adjoint representation that, up to a factor of 2, satisfy the same algebra as the Pauli spin matrices, that is [T a , T b ] = iε abc T c . Hence, the adjoint representation is (T a ) bc = −iε abc , i.e. to be explicit such that Tr(T a T b ) = 2δ ab and therefore T r(φ 2 ) = 2 3 a=1 φ a φ a . The SU (2)-symmetry in the adjoint representation for each generator T a is therefore so that the infinitesimal changes to the fields φ a i result to We will utilize this expression below. In more a compact form the Lagrangian in (2.42) can be expressed equivalently as where repeated indices are summed over the appropriate index sets i, j, µ, ν ∈ {1, 2} and a, b ∈ {1, 2, 3}. The matrix M 2 is defined as (2.47) and E as in (2.10). The covariant derivative in the adjoint representation acting on a complex field takes on the form Pursuing here a pseudo-Hermitian approach we perform a similarity transformation on the Lagrangian in (2.46) with Dyson map that maps the complex Lagrangian L ad 2 to a real Lagrangian where the matrix H is defined as (2.51) and I as in (2.10).

The SU (2)-symmetry preserving and breaking vacua
To find the different types of vacua φ 0 , we need to solve again δV = 0. The corresponding functional variation of the Lagrangian in (2.50) leads to the three sets of equations H + 2gR 2 E φ 0 a = 0, a = 1, 2, 3, where the N a are normalization constants. Given the solution in (2.53), the relation for R 2 imposes the additional constraint R 2 = N 2 1 + N 2 2 + N 2 3 on these constants. Expressing the Lie algebra valued vacuum field φ 0 i = φ 0 i a T a in the matrix form of the adjoint representation (2.43) we obtain We can now apply the SU (2)-symmetry to the vacuum state in the form as solutions for φ vac . Evidently, these states are linearly dependent as According to Goldstone's theorem the states v 0 i should be eigenvectors of the squared mass matrix with eigenvalue 0. As only two of them are linearly independent we expect to find two massless Goldstone bosons, which in our gauged model correspond to "would be Goldstone bosons". Hence the SU (2)-symmetry has been broken down to a U (1)symmetry, so that the group theoretical argument predicts two Goldstone bosons equal to the dimension of the coset SU (2)/U (1).

The squared mass matrix
Expanding the Lagrangian in equation (2.46) about the vacuum solution gives 60) where the last two terms originate from expanding the covariant kinetic term. The Hessian matrix is then computed by differentiating (2.52) once morê from which we obtain the non-Hermitian squared mass matrix as The entries in the rows and columns of M 2 are labeled as (φ 1 1 , φ 1 2 , φ 2 1 , φ 2 2 , φ 3 1 , φ 3 2 ) =: Ψ. The six eigenvalues λ of M 2 are then computed to with κ := 3µ 4 /2m 2 2 − m 2 2 /2 − m 2 1 . We can now verify that the three vectors v 0 i in (2.56)-(2.58), corresponding to the infinitesimal changes of the vacuum (2.53) under the action of the SU (2)-symmetry, are indeed vectors in the null space of M 2 . Due to their linear dependence we may choose two of them to be associated with the two massless "would be Goldstone bosons".
We notice that the eigenvalues in (2.63) do not depend on the choice of the three normalization constants N a , since all of these vacua are equivalent as they are related by SU (2)-symmetry transformations. The physical regions of the model are determined by the requirement that the eigenvalues are real and positive. Taking now account of the possibility that m i ∈ R or m i ∈ iR, by allowing for different signs in front of the m 2 i terms in setting m 2 i → c i m 2 i , we find that the model does not possess any physical region when c 1 = c 2 = ±1 and physical regions when c 1 = −c 2 = ±1 as argued also in the previous section.

The would be Goldstone bosons
Let us now identify the two massless Goldstone bosons ψ Gb 1,2 in the different PT -regimes by the same procedure as previously explained in [38,39], with the difference that they will be made to vanish due to the presence of the gauge bosons. In terms of the original scalar fields in the model we identify the Goldstone bosons by evaluating where the matrix U diagonalizes the squared mass matrix by U −1 M 2 U = D with diag D = (λ 1 , λ 2 , λ 3 , λ 4 , λ − , λ + ) and diagÎ = {I, I, I}. In the PT -symmetric regime the similarity transformation U is well defined by where the v i are the eigenvectors of M 2 . Up to normalizations constants for each eigenvector, we obtain in our example the concrete expressions Here we should stress that M 2 is only the analogue to a Hamiltonian and we are not attempting to construct an associated Hilbert space. For convenience we take now N 1 = N 2 = 0, N 3 = R and compute , and ψ Gb 2 := We note that det U = λ 3 λ 4 (λ − − λ + )µ 6 R 4 , indicating the breakdown of these expressions at the exceptional points when λ − = λ + , the zero exceptional point when λ 3 = λ 4 = 0 and at the trivial vacuum when R = 0, as previously observed in [38,39]. However, at the exceptional point we may still calculate the expressions for the Goldstone boson when taking into account that in this case the two eigenvectors v − and v + become identical. In order to obtain two linearly independent eigenvectors when the squared mass matrix is converted into its Jordan normal form we multiply two entries of the vector v + by some arbitrary constants α = β as With this change the matrix U becomes invertible as det U = λ 3 λ 4 (β − α)(m 2 2 + κ)N 2 1 µ 6 R 2 . We may now evaluate the expression in (2.65) obtaining the same formulae for the Goldstone bosons as in (2.68). At the zero exceptional point it is not possible to identify the Goldstone in terms of the original fields in the model.

The mass of the vector gauge boson
Finally we calculate the mass of the gauge vector bosons by expanding the minimal coupling term in equation (2.50) around the symmetry breaking vacuum (2.55) where we used the standard identity ε abc ε ade = δ bd δ ce − δ be δ cd . A convenient choice for the normalization constants N i that is compatible with (2.53) and diagonalizes (2.69) is to set two constants to zero and the remaining one to R. For instance, taking N 1 = N 2 = 0, N 3 = R the only nonvanishing terms in (2.69) are Thus for µ 4 = m 4 2 and R = 0 we obtain two massive vector gauge bosons m (1) g and m (2) g , that is one for each "would be Goldstone boson". When µ 4 = m 4 2 , that is then model is at the zero exceptional point of type I, the gauge mass vector bosons remain massless. This feature is compatible with our previous observations in [38,39] and above, that at these points the Goldstone bosons can not be identified.
We notice here that the two massive vector gauge bosons are proportional to the CPT -inner product of the symmetry broken vacuum solution (2.72) Hence the vanishing of the mass for the vector gauge bosons at the two types of zero exceptional points can be associated to the vanishing of the CPT -inner product at these points. This is reminiscent of the vanishing of the CPT -inner product at the standard exceptional points, which is responsible for interesting phenomena such as the stopping of light at these locations in the parameter space [45,46]. We note, however, a key difference between the two scenarios: While the CPT -inner product in (2.72) is devised on the eigenvector space of squared mass matrix, the latter is a CPT -inner product on the Hilbert space.

Conclusions
Employing a pseudo-Hermitian approach we found that the Higgs mechanism applies in the usual way in the CPT -symmetric regime by giving a mass to the vector gauge bosons and preventing Goldstone bosons to exist, which was also found in [37,41] using different approaches. As in [41] we also observed that in the spontaneously broken CPT -symmetric regime the vector gauge bosons become massive and the Higgs mechanism is in tact. However, as in this regime other particles acquire complex masses it has to be discarded as non-physical for that reason. Even though technically one needs to treat the standard exceptional point differently from the other regimes, the main principle of the Higgs mechanism still holds up. In contrast to the finding in [41], we observed that the Higgs mechanism breaks down at the zero exceptional points, which was also observed in [37]. We find the same characteristic behaviour, i.e. the matching of the amounts of massive vector gauge bosons and "would be Goldstone bosons", for the complex scalar fields taken in the fundamental as well as in the adjoint representation. The vanishing of the mass for the vector gauge bosons coincides with the vanishing of the CPT -inner product on the eigenvector space of squared mass matrix.
Obviously there are many interesting extensions to these investigations, such as for instance the treatment of models with a more involved field content or different types of continuous symmetries. is preserved on both sides. The zero exceptional point occurs when two eigenvalues coalesce at zero, hence the name.
We consider here a (3 × 3)-matrix of a very generic form that occurs for instance as a building block of the squared mass matrix in the model discussed in [39], see equation (3.48) therein, Here we carry out the discussion for a Hamiltonian H, having in mind the analogy to the squared mass matrix. The determinant is easily computed to det H = Aκ − CW 2 , κ := V 2 −BC. In order to obtain a zero eigenvalue, λ 0 = 0, we enforce now the determinant to vanish by setting A = W 2 C/κ. The other two eigenvalues then become respectively. We abbreviatedκ := ± κ (κ + V 2 − C 2 ) ± 2κ 3/2 V andκ := √ V 2 − C 2 . The first set of eigenvalues in (A.4) correspond to the standard exceptional point and the second set to the zero exceptional point.
Next we calculate the bi-orthonormal basis from the normalized left and right eigenvectors u i , v i , i = 0, ±, respectively, for with U = diag(1, −1, 1) and normalization constants N 0 = κ 2 + W 2κ2 , N ± = V 2 κ 2 + W 2 [V 2 − C(C + λ ± )] 2 − (C + λ ± ) 2 κ 2 . By construction these vectors satisfy the orthonormality relation u i · v j = δ ij . We observe now that at the standard exceptional point the two eigenvectors for the non-normalised (N ± become zero at the exceptional points) eigenvalues λ e ± coalesce, which distinguishes exceptional points from standard degeneracy. The left and right eigenvectors become in this case v e,r ± = κ, V √ κ − κ, C √ κ , v e,r 0 = (−Cκ, Cκ, Vκ) , (A.7) v e,l ± = κ, κ − V √ κ, C √ κ , v e,l 0 = (−Cκ, −Cκ, Vκ) , (A.8) with v e,l ± · v e,r ± = 0, and v e,l 0 · v e,r 0 = C 2 κ 2 −κ 2 + V 2κ2 . (A.9) Similarly, at the zero exceptional point the eigenvectors for the eigenvalues λ 0 and λ 0 − coalesce, which qualifies this point also to be called "exceptional" in the standard terminology. In this case the left and right eigenvectors become with v e,l + · v e,r + = C 3 + BV 2 − 2CV 2 2 , and v 0,l · v 0,r = 0. (A. 12) In order to understand the key difference between these two types of exceptional points we consider at first the eigenvalues (A.2) near the critical values in (A.5). Concerning the standard exceptional points we note that the two eigenvalues become identical when τ → 0. Thus writing τ /C 2 = (W 2 − (W e ) 2 )(W 2 −W ), withW being the second root of the polynomial in W 2 under the square root, it is now clear that if we consider the eigenvalues as functions of W 2 the argument of the square root has different signs for W 2 = (W e ) 2 + ǫ and W 2 = (W e ) 2 − ǫ. Hence the eigenvalues are real on one side of the exceptional point in the W 2 -parameter space and complex on the other. In contrast none of the eigenvalues becomes complex in the neighbourhood of the critical value W 0 .
For completion we also report the Dyson map and hence the metric operator for which the same behaviour may be observed. Using the operator that diagonalises the non-Hermitian Hamiltonian H we verify the pseudo and quasi Hermiticity relations We observe that the map break down at both exceptional points, i.e. det η = 0 for the critical values W e and W 0 , and on one side of the standard exceptional point. In all other regions of the parameter space it holds. Thus we find the same behaviour as already observed for the analysis of the eigenvalues.