1 Introduction

In 2016, with the detection of the gravitational waves, we have witnessed the last and most impressive experimental confirmation of Einstein gravity; in 2012, with the discovery of the Higgs scalar, we have had the final and most important experimental confirmation of the validity of the standard model of particle physics. But despite all the success that the two theories have separately, nevertheless they still suffer the lack of a combined framework.

Even worse is that, theoretically, Einstein gravitation appears to be affected by the Penrose–Hawking theorem, stating that if the dominant energy condition holds then the high density gravitationally-induced formation of singularities is inevitable; in turn this means that the avoidance of singularity formation due to gravity can be done only by invalidating the condition on the energy, and thus what we would have to do would be to modify the structure of the energy tensor by allowing different contributions of matter distribution: this has been attempted in the past several times in [1,2,3,4,5,6,7,8,9,10,11,12] by employing a variety of spinning fields, and it is a general feeling that, if not spin, at least some new type of matter field should enter the game in question. Along this line, it is counter-intuitive to think that the trick might be done by a field that is already known, and which does not even take into account spin, such as the Higgs scalar [13]; furthermore, it seems unlikely that the standard model contain the solution for a problem pertaining to gravity, since the two appear to have little in common. Just the same, Higgs physics can change the energy tensor in its structure, so to make the dominant energy condition invalid, and therefore to have the gravitationally-induced singularity avoided.

2 Higgs sector of standard model with Einstein gravity

As it is known, the standard model (SM) is the model built on the local \(U(1) \times SU(2)\) symmetry group enriched with a mechanism of spontaneous symmetry breaking as induced by the Higgs scalar quartic potential: having the break-down of the symmetry, we can choose the so-called unitary gauge in which, after diagonalization of the mass matrix, we obtain the final Lagrangian, as it is presented in common textbooks; here we will not follow the usual textbook procedure of presenting the Lagrangian, but we will follow the different path of giving what follows from the Lagrangian via the variational approach, not just because it is an alternative way to present the SM, but also because there is more information in the field equations.

The entire set of SM field equations is in [13], and again we have no need to represent it here in full; the sector in which we are interested is the one concerning the Higgs scalar, together with fermions, in Einstein gravity, and so we are going to consider the field equations in the limit in which all gauge fields are neglected, and where massless fields are left out: hence, the Higgs-electron-gravity part is described in terms of the field equations given by

$$\begin{aligned}&\nabla ^{2}H+m_{H}^{2}\left( \frac{H^{2}+3vH+2v^{2}}{2v^{2}}\right) H +\frac{m_{e}}{2v}{\overline{e}}e=0 \end{aligned}$$
(1)

and

$$\begin{aligned}&i\varvec{\gamma }^{\mu }\varvec{\nabla }_{\mu }e- m_{e}\frac{H}{v}e-m_{e}e=0 \end{aligned}$$
(2)

with

$$\begin{aligned} R_{\alpha \mu }- & {} \frac{1}{2}g_{\alpha \mu }R=\frac{i}{8} ({\overline{e}}\varvec{\gamma }_{\alpha }\varvec{\nabla }_{\mu }e -\varvec{\nabla }_{\mu }{\overline{e}}\varvec{\gamma }_{\alpha }e\nonumber \\&+{\overline{e}}\varvec{\gamma }_{\mu }\varvec{\nabla }_{\alpha }e -\varvec{\nabla }_{\alpha }{\overline{e}}\varvec{\gamma }_{\mu }e)\nonumber \\&+\left( \nabla _{\mu }H \nabla _{\alpha }H-\frac{1}{2}\nabla _{\rho }H\nabla ^{\rho }Hg_{\alpha \mu }\right) \\&+\frac{1}{2}m_{H}^{2}H^{2}\left( \frac{H+2v}{2v}\right) ^{2}g_{\alpha \mu } -\frac{1}{8}v^{2}m_{H}^{2}g_{\alpha \mu }\nonumber \end{aligned}$$
(3)

in which H is the Higgs scalar, e is the electronic fermion field, \(g_{\alpha \mu }\) is the metric of the space-time and \(\varvec{\gamma }^{\mu }\) are the Clifford matrices given in tetradic form. Notice that the covariant derivative applied to the scalar \(\nabla _{\mu }H\) is written in a covariant notation but could simply be written as the partial derivative while the covariant derivative applied to the fermion \(\varvec{\nabla }_{\mu }e\) is the spinorial covariant derivative written in terms of the spin connection arising from the non-trivial metric of the space-time itself, and \(R_{\alpha \mu }\) is the Ricci scalar; and finally v is the non-trivial Higgs vacuum with \(m_{H}\) and \(m_{e}\) Higgs and electron masses. Notice that the Higgs scalar vacuum gives indeed rise to the negative cosmological constant term that entails the problematic contribution to the cosmological constant problem.

Once the gravitational field equations are provided, the Penrose–Hawking singularity theorem states that if

$$\begin{aligned}&R^{\rho \sigma }u_{\rho }u_{\sigma }>0 \end{aligned}$$
(4)

holds for some time-like vector \(u^{\alpha }\) then singularities form inevitably [12]. The (4) is the dominant energy condition.

To see if such a dominant energy condition (4) is valid, we need to write the gravitational field equations in the form in which the Ricci tensor is singled out: contracting the gravitational field equations (3) gives

$$\begin{aligned}&-R=\frac{1}{2}m_{e}\frac{H}{v} {\overline{e}}e+\frac{1}{2}m_{e}{\overline{e}}e\nonumber \\&\quad -\nabla _{\rho }H\nabla ^{\rho }H+2m_{H}^{2}H^{2}\left( \frac{H+2v}{2v}\right) ^{2} -\frac{1}{2}v^{2}m_{H}^{2} \end{aligned}$$
(5)

where we have used the fermion field equations and which can be plugged back into (3) to give the final

$$\begin{aligned} R_{\alpha \mu }= & {} \frac{i}{8} ({\overline{e}}\varvec{\gamma }_{\alpha }\varvec{\nabla }_{\mu }e -\varvec{\nabla }_{\mu }{\overline{e}}\varvec{\gamma }_{\alpha }e\nonumber \\&+{\overline{e}}\varvec{\gamma }_{\mu }\varvec{\nabla }_{\alpha }e -\varvec{\nabla }_{\alpha }{\overline{e}}\varvec{\gamma }_{\mu }e) +\nabla _{\mu }H \nabla _{\alpha }H\nonumber \\&-\frac{1}{4}m_{e}\frac{H}{v}{\overline{e}}eg_{\alpha \mu } -\frac{1}{4}m_{e}{\overline{e}}eg_{\alpha \mu }\nonumber \\&-\frac{1}{2}m_{H}^{2}H^{2}\left( \frac{H+2v}{2v}\right) ^{2}g_{\alpha \mu } +\frac{1}{8}v^{2}m_{H}^{2}g_{\alpha \mu } \end{aligned}$$
(6)

which can now be inserted into the dominant energy condition for final evaluation. After this (4) becomes

$$\begin{aligned}&\left[ \frac{i}{4}({\overline{e}}\varvec{\gamma }_{\alpha }\varvec{\nabla }_{\mu }e -\varvec{\nabla }_{\mu }{\overline{e}}\varvec{\gamma }_{\alpha }e) +\nabla _{\mu }H \nabla _{\alpha }H\right. \nonumber \\&\quad -\frac{1}{4}m_{e}\frac{H}{v}{\overline{e}}eg_{\alpha \mu } -\frac{1}{4}m_{e}{\overline{e}}eg_{\alpha \mu }\nonumber \\&\quad -\frac{1}{2}m_{H}^{2}H^{2}\left( \frac{H+2v}{2v}\right) ^{2}g_{\alpha \mu } +\left. \frac{1}{8}v^{2}m_{H}^{2}g_{\alpha \mu }\right] u^{\alpha }u^{\mu }>0\nonumber \\ \end{aligned}$$
(7)

and this is the condition that we are next going to discuss.

3 Effective approximation and comoving frame

In principle the discussion would rely on finding some exact solution and see whether on shell (7) is valid or not.

The actual discussion does not benefit from the leisure of exact solutions, and thus an alternative study must be performed; because the Higgs scalar is massive, the effective approximation can always be done and this allows the Higgs scalar degree of freedom to be integrated away.

The effective approximation is implemented by requiring that within the Higgs field equation (1) the dynamical term be irrelevant compared to the mass term as

$$\begin{aligned}&m_{H}^{2}\left( \frac{H^{2}+3vH+2v^{2}}{2v^{2}}\right) H +\frac{m_{e}}{2v}{\overline{e}}e\approx 0 \end{aligned}$$
(8)

which in the limiting case of high densities further approximates down to the easier

$$\begin{aligned}&H^{3}+\frac{vm_{e}}{m_{H}^{2}}{\overline{e}}e\approx 0 \end{aligned}$$
(9)

allowing to integrate the Higgs scalar degree of freedom.

When equation (9) is plugged into (7) we obtain

$$\begin{aligned}&\left[ \frac{i}{2}({\overline{e}}\varvec{\gamma }_{\alpha }\varvec{\nabla }_{\mu }e -\varvec{\nabla }_{\mu }{\overline{e}}\varvec{\gamma }_{\alpha }e) +\frac{1}{4}\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}} g_{\alpha \mu }\right] u^{\alpha }u^{\mu }>0\nonumber \\ \end{aligned}$$
(10)

in which we have kept only the most relevant terms.

When instead equation (9) is plugged into (2) and we multiply the result by the conjugate spinor we get

$$\begin{aligned}&i{\overline{e}}\varvec{\gamma }^{\mu }\varvec{\nabla }_{\mu }e +\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}}=0 \end{aligned}$$
(11)

where again we have kept only the most relevant terms and which can be split in space and time parts as

$$\begin{aligned}&i{\overline{e}}\varvec{\gamma }^{t}\varvec{\nabla }_{t}e= -i{\overline{e}}\vec {\varvec{\gamma }}\cdot \vec {\varvec{\nabla }}e -\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}} \end{aligned}$$
(12)

in a form that will be needed to evaluate the time derivative of the fermion field in the dominant energy condition.

In fact, time-like vectors can always be boosted into a comoving frame, and there (10) takes the form

$$\begin{aligned}&\frac{i}{2}({\overline{e}}\varvec{\gamma }_{t}\varvec{\nabla }_{t}e -\varvec{\nabla }_{t}{\overline{e}}\varvec{\gamma }_{t}e) +\frac{1}{4}\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}}>0 \end{aligned}$$
(13)

which with (12) becomes

$$\begin{aligned}&-\frac{i}{2}({\overline{e}}\mathbf { \varvec{\gamma }}\cdot \vec {\varvec{\nabla }}e -\vec {\varvec{\nabla }}{\overline{e}}\cdot \vec {\varvec{\gamma }}e) -\frac{3}{4}\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}}>0 \end{aligned}$$
(14)

as the dominant energy condition that we have to discuss.

4 The scaling behaviour

To discuss the dominant energy condition (14) we may start from considering that in absence of the Higgs sector the only term that remains is the dynamical term

$$\begin{aligned}&K=-\frac{i}{2}({\overline{e}}\vec {\varvec{\gamma }}\cdot \vec {\varvec{\nabla }}e -\vec {\varvec{\nabla }}{\overline{e}}\cdot \vec {\varvec{\gamma }}e) \end{aligned}$$
(15)

with K kinetic energy, always positive, hence recovering the old Penrose–Hawking result. But from a more general perspective, when this Penrose–Hawking analysis is done accounting for the Higgs sector, we get

$$\begin{aligned}&K-\frac{3}{4}\left( \frac{m_{e}{\overline{e}}e}{\sqrt{vm_{H}}}\right) ^{\frac{4}{3}}>0 \end{aligned}$$
(16)

displaying a remarkable difference: indeed, as an analysis of the scaling behaviour indicates, the kinetic energy term cannot suppress the Higgs-induced fermionic term.

In the evaluation of the dominant energy condition in this general case, both contributions are to be kept and the two have opposite effects, so that (16) no longer holds as a necessity: in case the kinetic energy dominates, some singularities could still form; but when the Higgs-induced fermionic term dominates, singularities are avoided.

Moreover, cases with \(K=0\) but \({\overline{e}}e\ne 0\) entail unavoidable failure of (16), and singularities become impossible in general; these situations are for instance those in which the fermionic fields can form condensed states.

We shall discuss one further example next.

5 Some comments

We have seen that in presence of the Higgs sector the dominant energy condition is given by (16) with K being the kinetic energy, always positive, while the other term is the Higgs contribution proportional to the fermionic density, always negative: hence (16) does not hold in all cases where the former term is irrelevant compared to the latter term, and singularities are avoided; and (16) can not hold in all cases in which \(K=0\) but \({\overline{e}}e\ne 0\) and there singularities are impossible. Cases of zero kinetic energy but non-zero fermion density are those given by cold and dense distributions of matter, like condensates; however, this is even more so for extreme situations such as those that are encountered in a Black Hole. As a consequence we should expect Black Holes to display no singularity.

Notice that it is rather reasonable that the Higgs scalar produce an interaction between fermions with an attractive character: indeed, this can be seen as the remnant of the attractive interaction between left-handed and right-handed chiral parts that the Higgs needs to exert in order to give mass to the fermion. In fact, it is the interplay between chiral parts what gives mass to the fermion field.

Notice also that it may look paradoxical that this attractive interaction be responsible for a counter-balance of the gravitational force, but there is quite a simple explanation for this phenomenon in the Einsteinian theory of gravitation: in Einstein theory of gravity the gravitational force is thought to be encoded within the curvature of the space-time, which in general gives rise to gravitational pull because matter distributions in general have a positive energy. However, an attractive interaction, entailing a potential well, dispenses negative contributions inside the total energy, and in situations in which these are the most relevant of the contributions the energy reverts to negative, and the space-time curvature relaxes.

As we have said, this may look rather counter-intuitive, but it is a natural consequence of the Einsteinian gravity, where the gravitational field is the space-time curvature and as such it is susceptible of being positive as well as being negative; and similarly it is the natural consequence of the presence of the Higgs sector, and the fact that it does give rise to a generally attractive interaction.

It is also important to notice that Einstein gravity with the Higgs scalar but with no fermion field would still be insufficient to have the singularity avoidance.

6 Conclusion

In this paper, we have considered the SM after symmetry breaking, and isolating only the Higgs-electron sector, we have exhibited the field equations, and we have considered the Einstein gravitational field equations sourced by the energy tensor of this scalar-fermion sector; to discuss the influence of the scalar-fermion sector in the dominant energy condition, we have assumed the effective approximation for the Higgs field and we have boosted into the frame comoving with the electronic field: we found that the dominant energy condition in this case accounted for the kinetic term, always positive, plus the fermion-Higgs term, always negative. This makes the dominant energy condition fail to be valid as a necessity; we have recalled cases where the dominant energy condition does not hold, and cases in which the dominant energy condition cannot hold at all. Therefore, singularity formation is sometimes avoided, and in some situations it is totally impossible.

We have stated that such circumstances are those that we have for example in Black Holes: because Black Holes are gravitational structures constituted by fermion fields in interaction through the Higgs scalar, being very cold and highly dense they are precisely the case of object for which the dominant energy condition cannot be valid, so that the gravitational singularity formation is impossible.

As a consequence of this fact, we have that Black Holes cannot be gravitationally singular as a mathematical consequence of the theorem of Penrose and Hawking.

However, they may still be gravitationally singular as a consequence of other mechanisms: as a start, all these results may change considerably without the effective approximation on the Higgs field; additionally, there may be differences of dramatic importance when employing field equations before symmetry breaking in SM. Finally, Black Holes may still form gravitational singularities in Einstein gravity’s extended versions or its alternatives.

For example, the lack of an effective approximation, or even worse the thorough lack of the symmetry breaking, makes it unreasonable to have these results stretched as to include the gravitational singularities we could find in Big Bang cosmologies, and therefore we cannot have any confidence in arguing that the Big Bang itself could turn out to be free of any gravitational singularity.

There is a similar impossibility in having these results enlarged to Einstein gravity’s extensions or alternatives, where the different structure of the background requires a complete revision of the singularity theorems.

In these cases, Black Holes or cosmologies free of any singularity may still be possible [14].