Einstein-Cartan-Dirac gravity with $U(1)$ symmetry breaking

Einstein-Cartan theory is an extension of the standard formulation of General Relativity where torsion (the antisymmetric part of the affine connection) is non-vanishing. Just as the space-time metric is sourced by the stress-energy tensor of the matter fields, torsion is sourced via the spin density tensor, whose physical effects become relevant at very high spin densities. In this work we introduce an extension of the Einstein-Cartan-Dirac theory with an electromagnetic (Maxwell) contribution minimally coupled to torsion. This contribution breaks the $U(1)$ gauge symmetry, which is suggested by the possibility of a torsion-induced phase transition in the early Universe, yielding new physics in extreme (spin) density regimes. We obtain the generalized gravitational, electromagnetic and fermionic field equations for this theory, estimate the strength of the corrections, and discuss the corresponding phenomenology. In particular, we briefly address some astrophysical considerations regarding the relevance of the effects which might take place inside ultra-dense neutron stars with strong magnetic fields (magnetars).


I. INTRODUCTION
We have recently witnessed the birth of gravitational wave astronomy, where the LIGO/VIRGO Collaboration reported compelling evidence on the detection of gravitational waves, which is compatible with a scenario of binary black hole mergers predicted by General Relativity (GR) [1,2]. This finding has sparked the interest in probing the strong-field regime of GR via gravitational wave observations of compact objects [3]. At the writing of this paper around a half-dozen events have been detected including a neutron star binary merger [4], with its electromagnetic counterpart compatible with a short gamma ray burst [5]. In the aftermath of these observations, researchers have quickly gone to discuss how well their favourite gravitational models extending GR have fared against them [6][7][8][9][10][11][12][13]. But given the undeniable experimental success of GR [14] why would one consider going beyond it? Modified theories of gravity are indeed motivated by a variety of reasons. On the observational side they are introduced as alternatives to dark matter/dark energy scenarios. On the theoretical side, the need of an ultraviolet completion of GR and the unavoidable existence of space-time singularities deep inside black holes and in the early Universe have been troubling researchers for decades. For some reviews on the motivations of these theories and their phenomenology, see e.g. [15][16][17][18][19][20][21][22][23].
The traditional approach to formulate GR is to consider a symmetric rank-two metric tensor g µν endowed with a pseudo-Riemannian manifold. Parallel transport is mediated by an affine connection Γ α µν , which is assumed * cosmocabral@gmail.com † fslobo@fc.ul.pt ‡ drubiera@ucm.es to be symmetric, Γ α µν = Γ α νµ , and metric-compatible, ∇ Γ α g µν = 0, with a minimal coupling between matter fields and gravity ensuring that the equivalence principle (i.e., the universality of free-fall) is fulfilled. However, it is known that a manifold equipped with zero curvature and non-metricity but non-zero torsion (also known as Weitzenböck space-time) allows to build a theory of gravity that is fully equivalent to GR when the linear action in the torsion scalar is chosen. This is known as the teleparallel equivalent of General Relativity [24][25][26][27]. There is yet another formulation fully equivalent to GR, based on zero curvature and torsion but nonzero non-metricity, called symmetric teleparallel gravity, whose properties have begun to be unravelled very recently [28][29][30][31][32][33].
Given these equivalences in gravitational models under different space-time paradigms, it is natural to ask if there is any guiding principle. Specifically, what spacetime geometry and fundamental degrees of freedom can represent gravity? Indeed, the study of non-Riemann geometries is connected to the topic of gauge symmetries. In the gauge formulation of gravitation, the gauge principle is applied to space-time symmetries leading to non-Riemann geometrical structures such as torsion and non-metricity (depending on the local symmetry group). This principle can provide guidance for the study of the role of non-Riemann geometries in gravitation and in unified field theories, in connection to symmetries in physics and symmetry breaking phase transitions (for detailed reviews on gauge theories of gravity and its applications see [34][35][36] and references therein). The history of the gauge principle is actually very rich, dating back one hundred years to the original proposal by Weyl in his unified field theory of gravity (and electromagnetism). For this purpose, Weyl introduced a space-time (Weyl geometry) with curvature and a non-vanishing trace-vector part of the non-metricity [37]. Later, Kibble gauged the Poincaré global symmetries of Minkowski space-time to arrive at the Einstein-Cartan-Sciama-Kibble theory of gravity (EC for short), within a space-time with curvature and torsion (Riemann-Cartan geometry, or RC). Poincaré gauge theories of gravity (PGTG) have been investigated with special interest on Lagrangians quadratic in the curvature and torsion invariants, and in applications regarding cosmology, gravitational waves, and spherical solutions, see e.g. [38][39][40][41][42].
In this work we shall focus on EC theory [43] (the simplest example of PGTGs), which allows to consistently incorporate the intrinsic angular momentum (spin) of fermionic matter. Thus, in EC theory, in addition to the stress-energy matter sources, the spin energy density is also a source of the gravitational field. Torsion in EC theory becomes important in scenarios where high-spin densities are present. Although Cartan introduced the theory almost one century ago, it continues to trigger interest due to its non-singular solutions (in black holes and cosmology), by the bridge it establishes between fermionic spinors and gravity, and by its elegance and simplicity, as it possesses no free parameters besides Newton's constant. Within its many applications we underline bouncing cosmologies [44][45][46], inflation, cosmological constant and dark energy [47][48][49], perturbations and cosmic microwave background radiation [50,51], phase and signature transitions [52,53], or compact objects [54,55]. For EC gravity coupled to Dirac fields, one also finds applications in particle physics, see e.g. [56][57][58][59][60][61][62].
The need to go beyond EC theory was recognized long ago, mainly due to the fact that the theory is still non-renormalizable [63], but also because quadratic Lagrangians present a natural and theoretically preferable extension [64][65][66] (see also [67,68]). A path recently explored to extend this theory is the analysis of new (nonminimal) couplings between torsion and the matter fields [69][70][71][72][73]. In this sense, the coupling between torsion and electromagnetism has been carefully analysed in the literature with the result that, in general, this can be achieved by either changing the field equations with minimal/nonminimal couplings, or via the constitutive relations between the field strengths F = (E, B) and the excitations H = (D, H) (see [74,75]). Though it is usually assumed that torsion does not minimally couple to the electromagnetic field, since this breaks the U (1) gauge invariance (for details see e.g. [63,74]), the physics of phase transitions in condensed matter systems, superconductivity and early universe is permeated by processes that lead to spontaneous symmetry breaking, and in high density environments torsion can provide a physical mechanism to induce such a symmetry breaking. Since in EC theory torsion vanishes outside the matter sources and is negligible at low densities, the Maxwell equations remain valid for all phenomena that we can presently probe directly.
The main aim of this work is to study the (minimal) coupling of torsion with fermions and electromagnetism (massless spin one bosonic fields). Through the Cartan equations relating the space-time torsion to the mat-ter fields, this coupling induces non-minimal and selfinteractions of the matter fields, and provides a physical mechanism to generate a U (1) symmetry breaking for high densities and fields. In the broken phase, torsion provides an effective mass for the photon, with the electromagnetic potential obeying an extended Proca-like equation.
In our approach we consider first the regime in which torsion is only sourced by fermions, and extend it later to the general case where both fermionic and bosonic fields contribute to torsion via the corresponding spin energy densities. The first case is a simplifying ansatz, where the electromagnetic fields are influenced by the (background) space-time torsion but do not backreact on it. This case serves to illustrate some of the effects in the new dynamics due to the minimal coupling of bosons and torsion. The second case encodes the full dynamics with the bosonic spin contribution to torsion, which induces new non-linearities. Therefore, in this paper we address the most relevant features of EC-Dirac-Maxwell model with U (1) symmetry breaking as well as some of its physical implications.
This paper is organized as follows. In Sec. II we review EC gravity, focusing on the case where fermions are represented by a Dirac field. Sec. III contains the core results of this work, where we extend EC theory by introducing the electromagnetic field minimally coupled to torsion, and find the corresponding dynamics for gravitational, electromagnetic and fermionic sectors in two cases: (i) fermionic background torsion and (ii) the full case, including the bosonic backreaction to torsion via its spin tensor. We conclude in Sec. IV with a broad discussion of the phenomenological implications of these results, including some future perspectives.

II. EINSTEIN-CARTAN THEORY WITH FERMIONS
A. Einstein-Cartan gravity In general, any affine connection can be decomposed into three independent pieces while the third term is associated to non-metricity Q ρµν ≡ ∇ ρ g µν and is called disformation, Throughout this paper all objects with tilde will refer to expressions computed in a Riemannian space-time. From the general decomposition (1), keeping curvature and torsion but setting non-metricity to zero we obtain a Riemann-Cartan (RC) space-time geometry. In this case the curvature scalars in the Cartan connection and the Levi-Civita connection are related as (4) Now, selecting the linear Lagrangian in the curvature scalar yields the action with the following definitions and conventions: κ 2 = 8πG is Newton's constant, g is the determinant of the spacetime metric g µν , the curvature scalar R = g µν R µν is constructed out of the Ricci tensor R µν (Γ) ≡ R α µαν (Γ), and the matter Lagrangian, L m = L m (g µν , Γ, ψ m ), depends on the metric and the matter fields, collectively denoted by ψ m , and also on the contortion via the covariant derivatives. The above action resembles GR, but the fact that the connection has now an antisymmetric part yields new contributions to the standard Einstein equations and to the dynamical equations for the matter fields, as we shall see below.
To obtain the field equations for EC theory we start by varying the action (5) with respect to the contortion tensor K α βγ , which yields the so-called Cartan equations or is the spin density tensor with dimensions of energy/area, T β ≡ T γ βγ and s β ≡ s γ βγ are the torsion and spin (trace) vectors, respectively 2 . Cartan's equations (6) imply that, 2 It is useful to express the torsion tensor in terms of its irreducible components T λ µν =T λ µν + 2 3 δ λ [ν T µ] + g λσ ǫµνσρT ρ , where the traceless tensor obeysT λ µλ = 0, ǫ λµνρT µνρ = 0, Tµ is the trace vector andT λ ≡ 1 6 ǫ λαβγ T αβγ is the pseudo-trace (axial) vector.
analogously to curvature being sourced by the stressenergy of the matter sources, torsion is sourced by its spin density. These equations are linear and algebraic, which implies that outside regions with spin densities (and, in particular, in vacuum) they vanish identically.
The Lagrangian in (5) can be expressed in the following way with or, explicitly where we used Cartan's equations in the second term and we have neglected the surface term which does not contribute to the field equations. Therefore, variation of the action (5) with respect to the space-time metric g µν yields the generalized Einstein equations, which can be suitably written as whereG µν is the Einstein tensor computed with the Levi-Civita connection in Eq. (1). On the right-hand side, we have the effective stress-energy tensor which we split in the stress-energy tensor of the matter fields T µν = 2 δg µν , with C ≡ − 1 2κ 2 K λ K λ + K αλβ K λβα , which contains the corrections quadratic in torsion U ∼ κ −2 T 2 or in the spin variables U ∼ κ 2 s 2 , via Cartan's equations (6). It is important to point out that, in general, torsion also contributes to the stress-energy tensor T µν , since the covariant derivatives present in the kinetic part of L m introduce new terms depending on torsion via minimal couplings (non-minimal couplings can also be present). Since U ∼ κ 2 s 2 , Eq. (10) defines a typical density ρ C ∼ 10 54 g/cm 3 , known as Cartan's density. This is much higher than the nuclear saturation density, ρ s ∼ 10 14 g/cm 3 , though much lower than Planck's density ρ P ∼ 10 93 g/cm 3 , where quantum corrections to classical gravitation are expected to arise. Therefore, in principle, EC theory can only introduce significant physical effects in environments of very large spin densities, which might arise in the early universe or in the innermost regions of black holes.
Let us also note that, though the contribution in U µν from the spin density to the space-time metric (or curvature) is only significant at or above the Cartan density ρ C , one can show that torsion-induced effects are predicted at smaller densities in the dynamics of fermions (via the so-called Hehl-Datta term, as discussed below). Therefore, the physics of EC-Dirac systems, as well as in those where an additional electromagnetic contribution is considered (as we shall see in this work) should also be studied inside ultra-compact objects such as neutron stars, particularly in magnetars, and even in quark stars. It should also be mentioned that Cartan's energy scale is not always the correct scale to guide the physical reasoning. This is clear for elementary standard model fermions since, as shown in [76], there are important (compensating) physical mechanisms at very small distances (Planck's scale, not Cartan's distance scale which is derived from Cartan's density) which allow the regularization of the self-energy densities around the electroweak scales (in the case of leptons), corresponding to the observed masses (m e , m µ , m τ , ...).
Another aspect of the EC equations (6)-(10) worth emphasizing is that, in the absence of spin density, they boil down to the GR ones. Therefore, EC theory in vacuum does not propagate additional degrees of freedom beyond the two tensor polarizations of the gravitational field that propagate at the speed of light, thus being in agreement with the recent findings resulting from the LIGO-VIRGO Collaboration on the equality of the speed of light and of gravitational waves [5]. Nonetheless, physical mechanisms for generating gravitational waves in the very early universe can also have contributions from spin energy tensor fluctuations around the Cartan density, for example via its effective non-vanishing and time-varying quadrupole moment, or due to phase transitions (including symmetry breaking mechanisms induced by torsion effects)

B. Einstein-Cartan-Dirac theory
Let us consider, as the matter sector in the action (5), a free Dirac fermionic field with mass m minimally coupled to torsion. The corresponding Lagrangian density can be expressed as [77] for spinors ψ and their adjointsψ = ψ + γ 0 , and where the covariant derivatives are defined as where D µ andD µ are the (Fock-Ivanenko) covariant derivatives built with the Cartan connection and the Levi-Civita connection, respectively, and γ µ are the in-duced Dirac-Pauli matrices 3 obeying {γ µ , γ ν } = 2g µν I, where I is the 4 × 4 unit matrix and g µν is the space-time metric 4 . The Lagrangian can be expressed as Given the definition of the spin tensor in (7), we have This expression (as the previous one) is valid for any Dirac field minimally coupled to the RC spacetime geometry, it does not depend on any particular theory of gravity.
For this matter source the spin tensor is totally antisymmetric, i.e., 5 and is expressed in terms of the Dirac (axial) spin vector as (see e.g. [36,44,[56][57][58]) This Dirac pseudo-vector field will play a crucial role later. Accordingly, the Cartan equations simplify since torsion is completely antisymmetric and Eq. (6) becomes therefore, the above Lagrangian introduces an effective spin-spin interaction induced by torsion.
Using the Cartan equations we then have The dynamical stress-energy tensor of Dirac fermions in EC theory is then given by 3 The usual constant Pauli-Dirac matrices γ c , which obey γ a , γ b = 2η ab I, are related to the γ µ matrices via γ µ e a µ = γ a , where e a µ are the tetrads satisfying gµν = η ab e a µ e b ν and e a µ e µ b = δ a b , e c ν e µ c = δ µ ν and η ab is the Minkowski metric. 4 The Fock-Ivanenko covariant derivatives of spinors, in Riemann geometry, are given byDµψ = ∂µψ +(1/2)w ab µ σ ab ψ andDµψ = ∂µψ − (1/2)w ab µψ σ ab , wherew ab µ = −w ba µ are the spin connection components related to the Levi-Civita connection (also known as Ricci rotation coefficients), w ab µ =w ab µ + K ab µ is the RC spin connection and σ ab ≡ (1/2)γ [a γ b] are the generators of the Lorentz group in the spinor representation (GL(2,C)). 5 Here we used the identities . Moreover, we can also compute the form of the torsion-induced corrections on the right-hand side of the Einstein equations (10) where in this case C becomes simplified due to the fact that for Dirac spinors contortion is completely antisymmetric. Using the expressions above, the EC field equations (6)-(10) becomẽ and respectively. Cosmological solutions in EC-Dirac theory have been investigated in detail in the literature, see e.g. [34,36,44,55,58].

III. EINSTEIN-CARTAN-DIRAC-MAXWELL THEORY
Let us now generalize the action (5) to incorporate a minimal coupling between torsion and the electromagnetic field (for recent works on torsion-matter couplings see e.g. [69][70][71][72][73]). This can be directly implemented at the level of the matter fields by assuming the matter Lagrangian density where the Dirac Lagrangian is the same as in Eq. (12) therefore including a minimal coupling to the RC geometry, while A µ is the electromagnetic four-potential and j ν the electric charge current density of fermions. On the other hand, we now have the generalized Maxwell Lagrangian in a RC space-time (satisfying local Poincaré invariance) written as where λ is a coupling parameter setting the system of units, and the generalized field strength tensor is defined as (note that K λ where ∇ is the covariant derivative in RC space-time constructed with the independent connection Γ λ µν in Eq. (1), whileF µν = ∂ µ A ν − ∂ ν A µ is the standard field strength tensor when torsion is neglected. This expression of the electromagnetic field puts forward that, due to the presence of torsion in the minimal coupling, the second term in Eq. (28) breaks the U (1) local symmetry. More explicitly, the Lagrangian density in Eq. (27) becomes We will next proceed with the derivation of the gravitational, electromagnetic and fermionic field equations corresponding to the action (5) with the U (1)-breaking term just introduced.
A. Gravitational sector

Fermionic background torsion
Let us assume a background torsion resulting from the spin density of fermionic fields. Variation of the action (5) with respect to the metric for the above matter sources yields the gravitational equations in this case, which can be suitably written as withT and the new term (compare to Eqs. (24)) induced by the minimal and non-minimal interactions between fermions and electromagnetic fields takes the form where It is the (minimal) coupling of torsion with fermions and bosons that give rise to the non-minimal interactions between the matter fields, once the Cartan equations are used to replace the torsion components by the matter field variables. Under the ansatz that torsion is exclusively resulting from matter fields with half-integer intrinsic spin (fermions), it becomes completely antisymmetric and the Cartan equations in this case are still given by Eq. (19). This assumption is already implicit in the second term of Einstein's equations above. This choice corresponds to keeping the (pseudo-trace) axial vector part of torsion as the only non-vanishing components. Under such an ansatz, the interaction part (32) in the effective stress-energy tensor introduces terms both linear and quadratic in the spin density, all of which depend on the electromagnetic quantities. Accordingly, the value of the coupling constant λ determines the scale at which the electromagnetic contribution becomes nonnegligible. Thus, in principle, one could test torsion effects at spin densities smaller than the Cartan one, for sufficiently high electromagnetic fields, which suggests that new gravitational (metric) effects could be present in the core of magnetars [78] and (hypothetical) quark stars [79]. Significant effects are expected for polarized matter because the linear terms will not average to zero and can introduce stronger torsion (spin) contributions to the Einstein equations at lower densities given a sufficiently high electromagnetic potential. In summary, torsion contributions to the metric field equations scale with κ 4 s 2 for the pure EC spin density correction, and κ 4 λsF A and κ 6 λs 2 A 2 for the linear and quadratic U (1) symmetry breaking terms, respectively.
In general, the macroscopic description of a physical system is achieved through an averaging procedure. For simplicity, we will consider physical systems where the spin density obeys an approximate random distribution. This simplification (which is not valid in the presence of sufficiently intense magnetic fields that tend to align the spins) allows us to neglect the terms linear in the spin density and consider only the quadratic ones. For Dirac fermions, if we consider only the terms quadratic in torsion, we obtain after some algebraic manipulations Here, as usuals 2 ≡s λs λ , A 2 ≡ A λ A λ ands · A ≡s λ A λ .

Full approach including the spin contribution from the electromagnetic sector
If we also consider, besides the fermionic spin, the contribution from the generalized electromagnetic Lagrangian (29) to the (total) spin tensor, s λαβ = s M λαβ + s D λαβ , we obtain where s D α βγ is Dirac's spin tensor and represents the electromagnetic contribution to the spin tensor, i.e, s M λµν = δL U(1)break M /δK µνλ , which also depends on torsion due to the minimal coupling previously introduced. The new Cartan equations can be written as since Dirac's (completely antisymmetric) spin tensor has zero trace vector. With a bit of algebra we get where τ α βγ ≡ T α βγ + T γ δ α β − T β δ α γ is the modified torsion tensor and the first term on the right-hand side corresponds to Dirac's spin tensor previously introduced. These expressions show that it is not trivial to separate the purely geometric torsion functions from the matter fields.
Let us contract the indices α and γ to obtain which we shall use to find an expression for the torsion trace vector. From the equation above one gets the result T β A β = 0 and this can be used after contracting Eq. (38) with A α A γ to arrive at Therefore, from this system of two equations we easily solve for T β Proceeding in a similar manner, by contracting Eq. (38) with A α , after some algebra it is finally possible to transform the Cartan equations into a form in which the geometric torsion is separated from the matter fields, that is (41) which can be further simplified down to and we denotes  (45) This expression for torsion as a function of the matter fields can then be replaced in the matter Lagrangian (26), i.e, in the bosonic sector, and in the Dirac Lagrangian, which after some algebra can be written as meaning that Dirac fields only interact with the axial vector part of torsion,T λ ≡ 1 6 ǫ λαβγ T αβγ . This expression is valid for any Dirac field minimaly coupled to a RC spacetime geometry, regardless of the gravitational theory.
Since bosons are also contributing to the torsion with the corresponding spin tensor, the axial torsion vectoȓ T λ has now a new contribution (besides that of the Dirac axial spin vector). From Eq. (42), we obtain where we have dropped the D symbol in the Dirac axial spin vector. Substituting in the Dirac Lagrangian above we get, after some algebra The first term is Dirac's Lagrangian on a (pseudo) Riemann space-time, while the other terms come from the corrections of a RC geometry where torsion (given by the Cartan equations) is due to the spin tensors of fermionic spinors and electromagnetic fields. The first term inside the parenthesis corresponds to the well known spin-spin (axial-axial) contact interaction. Due to the presence of new fermionic-electromagnetic interactions induced by torsion, we see that the spin-spin contact interaction is now modulated at each point by the strength of the electromagnetic 4-potential (squared). The spin-spin effect is therefore affected locally by the electromagnetic potential at very high densities and fields, due to the κ 4 factor. The other two terms represent further (fermionic) spin-electromagnetic interactions. In the first of these, significant at very high densities and fields, the relative orientation (alignment) between the spin vector and the electromagnetic potential is relevant, which might suggest that this interaction could involve precession effects and possibly generate anisotropies in the spin distribution, for example via a macroscopic (averaged) alignment of the fermionic spin.
To proceed with the analysis of the gravitational sector, we need to compute the effective stress-energy tensor T eff µν = T µν + U µν . The dynamical stress-energy tensor δg µν can be computed from the matter Lagrangian (26), using Eqs. (46) and (50), with the torsion corrections sourced by the new Cartan equations (42). These torsion-induced corrections to the stressenergy tensor T µν correspond to non-minimal interactions and also self-interactions of the matter fields. We can now write whereT µν =T Dirac µν +T Max µν is the standard stress-energy tensor for the matter-fields in curved space-time (that is, the one of GR), Π Mint µν comes from the bosonic Lagrangian (29) and includes non-minimal boson-fermion interactions (induced by torsion) and also bosonic selfinteractions. Its explicit expression can be derived from Eq. (46) using the Cartan equations (42), and gives rise to a quite long expression containing many terms (non-minimal and self-interactions). If we consider the regime of approximate random distributions of fermionic spin, corresponding to unpolarized matter (in the sense of zero macroscopic intrinsic magnetic moment), we end up with a simplified expression, since all terms that scale linearly with the fermionic spin are neglected. On the other hand, Ξ Dint it corresponds to non-minimal fermion-boson interactions (induced by torsion) and also spin-spin fermionic selfinteractions. This tensor can be computed exactly as assuming random fermionic spin distributions, therefore neglecting the last term in (50).
Now we can compute the correction to the bosonic La-grangian as The last term here depends on the spinors via Dirac axial vectors λ and represents non-minimal boson-fermion interactions, while every other term in that expression corresponds to self-interactions 6 6 The third and fourth terms in the Lagrangian above can be rewritten as where the functions h(A) and t(A) are given by respectively. In the expression (54) only the last four terms depend on the spinors via the Dirac spin axial vector. The term L M self g µν corresponds to the purely electromagnetic terms of Eq.(53), i.e, the self-interactions. In absence of electromagnetic potentials we recover the second term of Eq.(21) of the EC-Dirac model. Let us illustrate these considerations with an example. Assuming homogeneity and isotropy (FRWL metric) and A µ = (φ(t), 0, 0, 0), which meansF µν = 0 ands αβγ = 0, we obtain for the combination above whose components give the energy density and the pressure terms respectively. Let us focus now our attention upon the geometrical term U µν = 2 With a bit of algebra we obtain Note that torsion appears in the first three terms, where it has to be replaced by the corresponding matter (spin) quantities via Cartan's Eq. (42). Formally, the tensor U µν obeys the following expression and the resulting expression, as a function of the matter fields, upon substitution of torsion via Cartan's equations in (45) is rather long and complicated. To simplify things we will assume again spatial homogeneity and isotropy, which allows us to write the torsion tensor therefore we find and Using these expressions, after some algebra we obtain and the corresponding stress-energy tensor is This completes the calculations for the effective stressenergy tensor T eff µν = T µν + U µν , that enters the righthand side of the Einstein equations As one can see from the expressions obtained in the right-hand side of this equation, if quarks can form a condensate in vacuum, corresponding to non-zero vacuum expectation values for the quadratic terms ∼s 2 in the above equations, then the model predicts an effective cosmological constant modulated by the bosonic dynamics, providing a dynamical dark energy from vacuum condensates. This is possible thanks to the non-minimal and self-interactions involving the square of the Dirac (axial) spin.
The total matter Lagrangian will give rise to extended Dirac and electromagnetic equations. To compute this we will now analyze the bosonic and fermionic field dynamics.

B. Electromagnetic sector
Let us now obtain the electromagnetic field equations. Consider the action constructed from Eq. (27) together with the usual source term j λ A λ . Varying it with respect to the vector potential A µ yields which can be conveniently rewritten as where we have defined the torsion-induced four-current (67) with K λ ≡ K α λα . As can be seen in the expression for the Lagrangian in Eq. (29), or in the field equations (66), the terms quadratic in the contortion or, equivalently, in the spin density, resemble Proca-like terms. From this analogy, the coupling between the electromagnetic fourpotential and the space-time torsion provides an effective mass for the photon m 2 γ ∼ λT 2 in physical environments where the U (1)-breaking phase transition takes place. The terms linear in torsion, on the other hand, reveal new physical effects due to the coupling between electromagnetism and torsion, which in this framework become significant for spin densities much lower than the Cartan density. That is, way before the manifestation of new metric effects that are implicit in Eqs. (30) and (32), the torsion (spin) of fermions start interacting significantly with electromagnetism, affecting Maxwell's dynamics. This is another motivation to consider physical effects of the full dynamics in astrophysical and cosmological environments with spin densities below the Cartan threshold, as in the core of neutron stars and in the early universe.
To complete the electromagnetic sector of the dynamics, we include the generalized conservation equation, which can be written as Alternatively, one can write where is the commutator of covariant derivation of an antisymmetric (0, 2)-tensor in RC space-time. We see that Dirac's current is not conserved, therefore, from a probabilistic semi-quantum description point of view the particle number can change due to intense gravitational fields.

Fermionic background torsion
Assuming the ansatz of a completely antisymmetric background torsion, as in the case where torsion comes from the background Dirac fermionic fields, we get the same form of the field equations but the torsion-induced current gets simplified (71) According to the minimal coupling between torsion and electromagnetic fields, as it is apparent from Eq. (28), only the antisymmetric part of the contortion tensor enters the electromagnetic sector in a RC space-time (at the Lagrangian level). However, for fermions both torsion and contortion are totally antisymmetric, so we have dropped out the brackets for antisymmetrization. In that case it is useful to express the Maxwell Lagrangian with torsion contributions, Eq. (29), as Under the assumption of the random spin distribution, from Eq. (71) and using the Cartan equations (19), we obtain the following (spin) torsion-induced four-current which arises from the interaction between the fermionic axial vector field and the electromagnetic 4-potential.

Full approach including the contribution of the generalized electromagnetism to the spin tensor
In this case, the Cartan equations are given by Eq. (42). We will consider for convenience the generalized current in the following form where we have used the fact that contortion is antisymmetric in the first two indices and also that K ν [λµ] = T ν λµ and K λ = 2T λ . Now, given the fact that the total matter Lagrangian can be written as where L D includes bosonic-fermionic interactions and is given by Eq. (50) and L M is given in Eq. (29), upon applying the variational principle with respect to the electromagnetic potential, we get a new generalized Maxwell equation in Eq. (66) given bỹ where comes from the (effective) Dirac Lagrangian (50) as .
Using now Eq. (42) we obtain a long expression for the torsion-induced current J ν with non-linear terms. One can also use the effective Maxwell Lagrangian in Eq. (53) to obtain where the last term is computed as and we have introduced the definitions These highly involved expressions can be interpreted as non-linear electrodynamics with non-minimal couplings between fermionic matter (spinors) and electromagnetic fields induced by the RC space-time geometry. These equations are simplified in two cases: (i) matter with a random distribution of fermionic spins, where we neglect all quantities linear in the Dirac spin, leaving only the quadratic ones which do not vanish after macroscopic averaging and (ii) the case of homogeneity and isotropy, with A = (φ, 0, 0, 0), ands = 0 =F . In the first case we obtain and in the second case, the simplified J ν is simply where the non-linearities (in the electromagnetic quantities) disappear and the equation above corresponds exactly to what we had in the first approach in Eq. (74).

Fermionic background torsion
Let us consider first the case in which the matter fields are fermionic spinors. The variation of the Dirac action in a RC space-time (given by the Lagrangian density in Eq. (12)) with respect to fermionic fields yields the Fock-Ivanenko-Heisenberg-Hehl-Datta equation [80] i γ µD µ ψ − mψ = where torsion was substituted by its source, the spin density of Dirac fermions, via the Cartan equations. Now we introduce electromagnetic fields minimally coupled to torsion, but without backreacting on it. In this case, the variation of the action (5) The Hehl-Datta term, ∼ κ 2 2 (ψγ ν γ 5 ψ)γ ν γ 5 ψ, which is cubic in the spinors, is already present in the usual EC-Dirac theory. This term represents a spin-spin contact interaction inside fermionic matter. For charged antifermions, after performing the charge conjugation oper- All cubic terms, similarly to the term having the fermionic charge, have flipped sign after the Ctransformation relative to the mass term. This behaviour is connected to the fact that the corresponding effective Lagrangian terms behave in an opposite manner under a C-transformation in relation to the rest of the terms in the Lagrangian [60].
It has been shown that the Hehl-Datta term, which corresponds to an effective axial-axial spinor interaction of repulsive nature, can provide important physical effects in the particle domain [56,57,[59][60][61][62]76], including a valid mechanism for generating a residual matter/antimatter asymmetry in the context of baryogenesis in cosmology, and has been shown to posses other applications, such as an effective cosmological constant [61] and non-singular configurations [62]. Such a term can be derived from an effective interaction Lagrangian of the form L int Hehl−Datta ∼ κ 2sµs µ . Analogously, the new cubic terms we have derived also come from similar effective Lagrangian terms quadratic in Dirac's axial (spin) vector L eff ∼ κ 4 λs 2 A 2 , and are induced from the coupling between torsion and the electromagnetic potential. These terms correspond to the quadratic ones appearing in Eq. (29). Therefore, the axial-axial or spin-spin contact interaction effect is potentially enhanced (at very high densities) by the presence of the electromagnetic four-potential. Moreover, in general the four-potential propagates, therefore a richer dynamics is induced in the effective spin-spin interaction. This scenario is of course compatible with the fact that we have broken the local (gauge) U (1) invariance under a phase transition above a certain critical value of the spin density. Accordingly, the vector potential that appears explicitly in the dynamical equations can be thought as representing physical degrees of freedom 7 .

Full approach, including the contribution of bosonic fields to the spin tensor
Previously, using Eq. (42) we arrived at the fermionic Lagrangian given by If we consider the total matter Lagrangian including the contribution from the bosonic (electromagnetic) side, we obtain the following extended Dirac (cubic) equation where and we have while T γµν M is the purely bosonic part of the torsion tensor. This equation can be considered in the approximation of space-time flatness and also in the non-relativistic limit. One can then solve the energy levels problem which is expected to reveal a kind of hyperfine structure that could be used to probe for the existence of torsion with high resolution spectrography. In fact, the correction terms in (50) can be interpreted as effective interaction potentials with ϕ ≡s 2 , χ ≡ A 2 , ζ ≡ A ·s and we neglected the term linear ins, for simplicity. Such analysis is currently under study and will be developed in a future work. To close this section, let us mention that for anti-particles we have: which is not exactly the same dynamics, suggesting possible applications for asymmetries and baryogenesis.

IV. CONCLUSION AND DISCUSSION
In this work we have studied the Einstein-Cartan-Dirac-Maxwell model with U (1) symmetry breaking and discussed its physical relevance. We considered a Dirac field and an electromagnetic field minimally coupled to torsion, which induces rich gravitational dynamics and non-linear fermionic and bosonic dynamical equations, including non-minimal and self-interactions. We considered two regimes: i) one in which torsion is sourced by fermions and ii) the full case with the contribution from both fermions and bosons to the total spin tensor entering in Cartan's equations.
In general, the effects for the space-time metric only become important at very high (spin) densities, as in the usual EC theory. For example, in the first approach with torsion generated by fermionic spin, torsion (or spin) contributions to the metric field equations scale with κ 4s2 for the pure EC correction, while the model with the U (1) symmetry breaking studied here introduces terms both linear and quadratic with torsion, that scale as κ 4 λsF A and κ 6 λs 2 A 2 , respectively. This has to be compared with the κ 2T µν contribution from the usual stressenergy tensor in GR. Thus, for very strong electromagnetic fields/potential, the term linear in the spin density could become important (in polarized matter) at densities slightly (but not significantly) below Cartan's typical density. On the other hand, the effects of torsion in the electromagnetic and fermionic sectors require a more careful analysis.
Let us discuss the electromagnetic dynamics. The generalized Maxwell theory include terms linear in torsion (also in the spin density) that become significant at densities much lower than Cartan's density, which should be taken into account in strong gravity regimes such as in the interior of astrophysical compact objects (neutron stars, magnetars, quark stars) and in the early Universe. These terms are non-negligible for polarized matter, i.e., for non-random spin distributions and, consequently, the presence of strong magnetic fields provide the adequate physical conditions for the study of the phenomenology associated with these corrections. For approximately random spin distributions, i.e., for unpolarized matter, only the quadratic terms (in torsion or in the spin density) are non-vanishing with its phenomenology being related to much higher densities. In any case, the presence of strong electromagnetic fields (potential) tend to enhance such effects.
When the U (1) symmetry is broken the corresponding (Noether) charge current is not conserved. Although the fermionic charge density and number density of the fermions is not conserved locally in this model, the equations suggest interpreting the terms of geometric origin as effective charge currents that compensate and balance the non-conservation of the usual charge current. In other words, by following this interpretation the space-time geometrodynamics would gain physical features, such as effective mass, spin or charge currents, when it couples to matter. When these terms are considered, then a new conserved quantity is clear. Another way to see this is to deduce the phenomenology associated to such an interpretation and search for possible observational tests of the predictions. In this context, this type of models where the stress-energy tensor or the charge current is not conserved in the usual sense, predict the creation of particles from the energy available in the space-time geometrodynamics, in strong gravity environments.
When the contribution from the bosonic sector to the spin tensor is taken into account, then the bosonic field propagates on a RC spacetime and backreact on its geometry. Since torsion in EC theory is given by an algebraic expression of the matter fields, one then gets non-minimal couplings between these but also selfinteractions. Therefore, we obtain effectively a non-linear dynamical equations for the bosonic fields. In fact, just as in the case of fermions where a linear Dirac field in RC space-time of the EC theory is equivalent to a nonlinear spinor in GR, also here the linear electromagnetic Lagrangian in the RC space-time leads to an effective non-linear electrodynamics in GR. Non-linear dynamics in the matter fields can emerge naturally from the (minimal) couplings of these fields with the extended spacetime geometries of gauge theories of gravity.
In the case of fermionic fields in EC theory, torsion effects can also become significant in environments where the density is lower than Cartan's density. This is not so commonly mentioned in the literature, on the contrary, much emphasis is put on the fact that in EC theory the effects of torsion in Einstein's equations, i.e., for the metric, are only significant at extremely high densities such as those found in the very early Universe or inside black holes. Since the Cartan equations imply K ∼ κ 2s , after its substitution in the Dirac equation i γ λ D λ ψ−mψ = 0, one obtains the (cubic) Hehl-Datta equation where the torsion-induced term will become significant at (spin) densities comparable to any strong-gravity regime where GR effects become important.
Let us stress that the Hehl-Datta term, which is related to an effective axial-axial (spin-spin) repulsive interaction, has been studied in connection to different physical mechanisms important for particle physics and cosmology, such as non-singular black holes, matter/anti-matter asymmetry and energy-levels, etc. Analogously, in our U (1) symmetry-breaking model similar cubic terms are present that are quadratic in the electromagnetic fourpotential. In this case, these torsion-induced corrections scale with κ 4 , which means that the corresponding physical effects (on the energy levels, generalized effective Feynmann diagrams, etc) will only become relevant at extremely high densities (Cartan's density or above, but still lower than Planck density). In this model, the mini-mal coupling between the electromagnetic potential and torsion induce, at the dynamical equation level, a nonminimal coupling between fermions and electromagnetic potential/fields, in the generalized Dirac equation. Formally, this follows after the substitution of torsion by its corresponding spin density source via Cartan's equations. The new terms are both linear and cubic in the spinors. The former introduces effects that will become relevant around the same densities as for the original Hehl-Datta term. These considerations motivate further study on the full EC-Dirac-Maxwell dynamics inside astrophysical compact objects.
Finally, let us mention several cosmological, astrophysical and particle physics applications that can be worked out from the theory considered in this work. In Cosmology one expects the possibility of non-singular models as in the usual EC model, and new physics during the torsion-dominated era. One should also expect the production of gravitational waves from the transitions between primordial phases: from the U (1)-broken phase to the U (1)-restored phase, and from the usual torsiondominated phase of EC to the radiation phase. These transitions can contribute to a stochastic gravitational wave background of cosmological origin, with possible imprints from the physics beyond the standard model.
On the other hand, the standard EC theory can prevent black hole singularities and, therefore, the research on whether one can have equilibrium configurations in compact objects denser than neutron stars, before the appearance of an horizon, is of utmost relevance. In our model we have physical mechanisms induced by torsion that act as an effective repulsive interaction, which could possibly provide the required pressure to balance the selfgravity of a newly born (unstable) neutron star. After the coalescence of two neutron stars in models of GW emission, it is usually assumed that the resulting object stabilizes to a neutron star or decays into a black hole (directly or after some relaxation time), due to GR instabilities, but in modified gravity, torsion/spin effects should allow for other equilibrium configurations, i.e, stable compact objects denser than neutron stars.
To conclude, it is necessary to investigate whether the astrophysical data about the final object might be reinterpreted using models with torsion. In our view, there are good motivations to consider gravitational models where non-Riemannian geometries, fermionic spin densities, and phase transitions become important, which can be tested with astrophysical, cosmological and gravitational wave observations. Work along these lines is currently underway.