Majorana Fermions, Supersymmetry Breaking, and Born-Infeld Theory

This review is devoted to highlight some aspects of the relevance of Majorana fermions in rigid supersymmetry breaking in four spacetime dimensions. After introducing some basic facts on spinors, and on their symmetries and reality properties, we consider Goldstino actions describing partial breaking of rigid supersymmetry, then focussing on Born-Infeld non-linear theory, its duality symmetry, and its supersymmetric extensions, also including multi-field generalizations exhibiting doubly self-duality.


Introduction
The present review is aimed to contribute to the Proceedings of the 53rd Course of the International School of Subnuclear Physics (held in Erice in Summer 2015), which also celebrated the Triumph of Ettore Majorana, the great disappeared Sicilian Physicist whose ideas and invention have permeated both theoretical and experimental physics. In Section 3 some basic features of spontaneously broken rigid supersymmetry and partial supersymmetry breaking are described. The Volkov-Akulov and the Born-Infeld Lagrangians naturally emerge in this framework, and they both contribute to the description of partially broken rigid N = 2 supersymmetry.
The electric-magnetic self duality enjoyed by the Born-Infeld theory is then extended in Section 4, in order to describe more general theories, following the work of Gaillard and Zumino [1]. These generalizations involve the Schrödinger formulation of Born-Infeld-type theories, forms of higher degree as well as geometrical mass terms in higher dimensions.

Majorana Fermions and Supersymmetry
Supersymmetry deals with basic spinorial entities called Majorana spinors. This notion originates from a fundamental observation of the Italian (Sicilian) physicist Ettore Majorana, who noticed that the four-dimensional Dirac equation admits "real" solutions if the γ-matrices are suitably chosen [2]. This is the so-called Majorana representation, in which the γ-matrices are real. Therefore the Dirac equation admits manifestly real solutions.
Spinors which are eigenstates of γ 5 are called Weyl spinors. From the previous properties, it follows that Weyl spinors have two components instead of four, and they are complex in D = 4.
Since we have been mentioning them a few lines above, it is here worth recalling also the anticommutator sectors of supersymmetry algebras (without central extensions) in D-dimensional Minkowski spacetime: where P µ is the momentum operator. The Latin capital indices run over the number N of independent spinor charges : A, B = 1, ..., N ; for D even, the automorphism group of the supersymmetry algebra (named R-symmetry group) reads In D = 11, if the massless spectrum is bound to have 2 as the highest spin, then N max = 1, and the corresponding M-theory superalgebra with central extensions reads as follows : D even : each ψ ± has 2 D− 2 2 components (which double for complex or pseudoreal ones); (20) D odd : each ψ has 2 D− 1 2 components (which double for pseudoreal ones).
Consequently, (pseudo)real spinor components of dimension 8, 16 and 32 are possible for D 6, D 10 and D 11, respectively.
The dimension D mod 8 determines the symmetry of the tensor products of spinor representations (denoted by ψ), which are defined by the morphism map where Λ k denotes the space of k-forms. The resulting symmetry of A, depending on D and k, is presented in Table 1. For a comprehensive treatment of spinors and Clifford algebras for arbitrary (s, t)-signatures of spacetime, see e.g. [3].
It should be recalled that the massless Dirac equation ((1) with m = 0) has an extra symmetry, which in the Majorana representation takes the form this is indeed a U(1) rotation, since γ 5 is anti-involutive : (γ 5 ) 2 = −1.
However, this symmetry is broken by the mass term for a Majorana fermion (called Majorana mass). Defining a complex Dirac spinor ψ D as where ψ 1 and ψ 2 are Majorana spinors, one can them have a U(1) symmetry which rotates ψ 1 into Then, one can define Therefore, ψ L and χ L (and ψ R and χ R ) have opposite U(1) phases, whereas ψ L and χ R (and ψ R and Thus, the Dirac mass term, in this notation, is mψ L χ L + h.c., which is of course U(1)-invariant.
In principle, if one gives up the U(1) symmetry, a Dirac fermion can have three types of mass terms, namely : where m is named Dirac mass, and M and N Majorana masses. As a result, the matrix has two eigenvalues : In particular, for N = 0 and M ≫ m, the eigenvalues become This result is at the basis of the see-saw mechanism [4,5].
Concerning the superpartner of the graviton, this is a spin-3 2 Rarita-Schwinger field, called gravitino. Under the Lorentz group SO(1, D −1) it is a vector-spinor ψ µα . If it is massless (namely, in the case of unbroken supersymmetry), it lies in an irreducible, γ-traceless representation of the little group SO(D − 2); so, assuming real (semi)spinors, the gravitino has a number of (on-shell ) independent components given by which for example yields D = 11 : # (ψ µα ) = 8 · 2 4 = 128; The simplest Lagrangian density of supergravity couples Einstein gravity with the Rarita-Schwinger density, yielding Assuming 2 as the highest spin in the massless spectrum (pertaining to the graviton), after [6] it is known that Poincaré supersymmetry coupled to gravity is possible for D 11, while conformal supersymmetry is possible for D 6; this result also stems from the classification of superalgebras done by Kac [7].
3 Spontaneous Symmetry Breaking of Rigid Supersymmetry : Goldstino Actions As holding for any action invariant under continuous symmetries, in theories invariant under one or more supersymmetry transformations there are conserved Majorana vector-spinor Noether currents [8] (in van der Waerden notation): The corresponding charges generate the (N -extended) supersymmetry algebra, whose non-trivial part in D = 4 (which we consider throughout this Section) is given by (14).
In N = 1 spontaneous breaking, the order parameter µ enters a term linear in the supercurrent [9,10] where G is a Majorana field, the Goldstino; note that µ has dimension 2 in natural units.
The supersymmetry current algebra implies that so that By generalizing this for several supercharges, one can conclude that (in absence of central extensions, as well as in the rigid case -no gravity coupled -, as we assume throughout this Section), either all supersymmetries are unbroken or they are all broken at the same scale f : One can evade this result by modifying the current algebra via terms not proportional to δ B A , namely by adding a contribution proportional to a constant matrix C B A in the adjoint of SU(N ) (i.e., traceless and Hermitian): As it holds in general for Goldstone particles, the Goldstino G is a massless fermion, with lowenergy self-interactions described by a Lagrangian invariant under a non-linear realization of supersymmetry, and depending on both N and C B A , because it can be shown that If k supersymmetries are unbroken, the Hermitian matrix V 1 + C has rank N − k in the vacuum, For N = 2, there are various examples realizing a spontaneous supersymmetry breaking [11,12,13,14,15], and the previous analysis in terms of SU(2) invariants was performed in [16]. The rigid limit of N = 2 supergravity theory reproducing partial supersymmetry breaking was considered in [17]. The unique Casimir of SU(2) is the squared norm of a 3-vector ξ x constructed via electric and magnetic Fayet-Iliopoulos (FI) terms where Q x = (m Λ x , e x Λ ), and Λ = 1, . . . , n, with n denoting the number of N = 2 vector multiplets. The wedge product ∧ is with respect to the symplectic structure inherited from the rigid special Kähler geometry endowing the vector multiplets' scalar manifold, and the matrix C B A acquires the explicit form with eigenvalues Thus, when V = √ ξ x ξ x , λ + = 0 and N = 2 rigid supersymmetry is broken down to N = 1, with order parameter The last term in the r.h.s. of Eq. (44) corresponds to a "vector" central charge Z A µ|B [18,19]; the only field not being inert is the gauge field of N = 2 vector multiplets.
For a general N -extended rigid supersymmetry broken down to N = 0 (corresponding, from the previous reasoning, to C B A = 0), the low-energy effective Lagrangian for N Goldstino fields is an obvious extension of the Volkov-Akulov (VA) action [20] derived in the N = 1 case, and in general reads It should be stressed that this Lagrangian contains a finite number of terms, because of the nilpotency of the Goldstino fields G A . Furthermore, L V A has a U(N ) R-symmetry, being invariant under the following non-linear supersymmetry transformations : On the other hand, the Goldstino action for N = 2 rigid supersymmetry partially broken down to N = 1 is described by the supersymmetric Born-Infeld (SBI) theory : this is given by a nonlinear Lagrangian L SBI (G α , F µν ) for an N = 1 vector multiplet whose chiral superfield strength W α = D 2 D α V contains the goldstino G at θ = 0 and the self-dual Maxwell field strength at the next order in θ, with the following two properties : where L BI is given by (94) further below.
L SBI exhibits two types of super-invariances: a manifest N = 1 supersymmetry, linearly realized in N = 1 superspace, and a second supersymmetry which is non-manifest and nonlinearly realized. In terms of microscopic parameters, the scale of the broken supersymmetry is the SU(2)-and symplectic-invariant quantity (50); see Sec. 4.4 further below.
It is here worth pointing out that the Goldstino action for a partial breaking N → N − k of rigid supersymmetry is only known for N = 2, corresponding to the non-linear limit of a quadratic action of N = 2 vector multiplets endowed with FI terms. The simplest case is given by a single N = 2 vector multiplet, which under the spontaneous breaking N = 2 → 1 reduces to an N = 1 vector multiplet W α and to an N = 1 chiral multiplet X. The presence of a non-vanishing matrix C B A allows the chiral multiplet X to acquire a non-zero mass m X : For m X large enough, X can be integrated out, giving rise to a non-linear theory for the fields G α and F µν of the N = 1 vector multiplet, described by the aforementioned L SBI (G α , F µν ) [21,23].

Electric-Magnetic Duality, Born-Infeld and Generalizations
Electric-magnetic (e.m.) duality is one of the most fascinating symmetries of non-linear (and thus interacting) theories of gauge fields.

Maxwell Theory
Free Maxwell theory is the prototype of e.m. duality invariant theory. In the vacuum, the Eqs. of and they are invariant under SO(2) rotations: This latter transformation, in terms of electric and magnetic fields Also in this notation it is immediate to see that Maxwell Equations in absence of sources (F oi = E i , Note that the Hamiltonian is e.m. duality invariant, while the Lagrangian is not. It should also be pointed out that the e.m. duality symmetry is not an internal symmetry, since it rotates tensors with pseudo-tensors (as such, it can be regarded as a sort of "bosonic chiral transformation").
Electric-magnetic dualities are transformations among two-form field strengths and their duals. In order to generalize the e.m. duality to p-forms, one demands that D/2 = p + 2, so that one has (p + 1)-form gauge fields coupled to sources, which are p-extended objects: p = 0 (pointlike), p = 1 (string), p = 2 (membrane), p = 3 (three-brane), etc.: When D/2 = p + 2 we have that D − p − 2 = p + 2, and p-dimensional extended objects can source both electric and magnetic fields. We thus have dyons, extended objects that carry both an electric and a magnetic charge. In all other cases, electric and magnetic objects extends in different space dimensions, (p, D − p − 4) respectively, so dyons cannot exist.

Non-Linear Electric-Magnetic Duality
The very idea of duality is to generalize the duality of free Maxwell theory to the case of interactions.
Interactions introduce non-linearities, which can be due either to the Maxwell field itself, or to other matter fields (as in electrodynamics) which interact with the e.m. field.
Non-linear (pure) electromagnetism can be seen as an effective theory of Maxwell theory in a medium that has a nonlinear response to the electric and magnetic fields. These nonlinearities are captured by the relations expressing the electric displacement − → D and the magnetic field − → H in terms of the magnetic induction − → B and the electric field − → E . The field equations (in absence of sources) are thus divided in a set of linear ones that are independent from the medium and a set describing the nonlinearities, i.e. the properties of the medium, hence these equations are called constitutive relations, Maxwell theory is recovered when In a theory based on a Lagrangian density L = L − → E , − → B , one has the Eqs.
Note that Eqs. (68) are invariant under the linear transformations  It is convenient to rewrite the nonlinear equations and the duality rotations using a relativistic formalism. We set as usual F µν = − → E , − → B and G µν = − → H , − − → D , so that in the vacuum: The dynamical equations then read dF = 0, dG = 0 i.e., while the constitutive relations become and, if we have a Lagrangian, In order to have compatibility with the equations of motion, this demands the following integrability conditions : The duality rotations correspond to transform which equivalently reads In absence of scalars, or if scalars are inert under duality rotations, the e.m. duality group is at most U(n) (in particular, for n = 1, one retrieves Maxwell e.m. duality group: U(1) ∼ SO (2)).
The constraints on the n field strength F Λ and G Λ in this case read For n = 1 (namely, U(1)), one gets For general e.m. theories where also matter fields are present the constitutive relations will also depend on these fields, that we generically denote ζ i , (and on their derivatives that we omit writing) so that we have G µν = G µν (F, ζ i ). Duality rotations now will include also a (nonlinear) transformation on the matter fields ζ i .
As shown by Gaillard and Zumino [1] (and reviewed in [24]) in this general case, under an infinitesimal duality rotation parametrized by the Lagrangian varies by where δL ≡ L(F + δF, ζ i + δζ i ) − L(F, ζ i ) is the total infinitesimal variation of the Lagrangian under the duality rotations of the field strengths and of the matter fields. If we expand δF in terms of F and G using the infinitesimal version of the duality rotation (82) then we obtain If there are no matter fields the left hand side is zero, the duality rotation group reduces then to We notice that in the special case of Maxwell theory G = ⋆F and therefore the infinitesimal Maxwell theory (that acts non-locally on the gauge potential A) was studied in [26]; however it does not correspond to a finite off shell symmetry. As we see from the second term in (90) in order for F ⋆ G to be a total derivative we have to use also the equations of motion dG = 0, that in Maxwell theory simply read d ⋆ F = 0.
Two important results are obtained from the previous formulae for general nonlinear theories.
1. Any Lagrangian density can be written as 1 4 F ⋆ G plus a duality invariant Lagrangian L inv , The term L inv provides information about the constitutive relations (i.e., the medium), in particular L inv = 0 vanishes in linear Maxwell theory without coupling to matter. On the other hand, for (quadratic) theories coupled to matter L inv provides information about some matter couplings (such as Pauli-like terms). In supergravity it is in fact possible to construct pairs of antisymmetric Lorentz tensors (H µν (ζ) , I µν (ζ)) (typically bilinear in fermions) which transforms as the couple (F, G). In this case, one has 2. The energy-momentum tensor is duality invariant if duality is a symmetry of the theory. In the absence of matter fields the trace of the energy momentum tensor is proportional to L inv : 1 4 Θ µ µ = L inv .

Born-Infeld Theory
The Born-Infeld (BI) model [27] is the simplest non-linear e.m. duality invariant theory. It was introduced to remove the divergence of the electrons self-energy in classical electrodynamics. This is obtained by having an upper bound on the electric field; indeed, for In fact, L BI has unique properties against instabilities created by the medium, and not many generalizations are known in the multi-field case (i.e., in presence of more than one Maxwell field).
Furthermore, by expanding L BI in power series of the coupling constant 1 µ one obtains and hence for µ → ∞, i.e., for weak fields, one recovers Maxwell theory.
The equations of motion of BI theory are a particular realization of the constitutive relations (here µ is a dimensionful coupling constant typically present in nonlinear theories, which was previously omitted for sake of brevity). Moreover, they satisfy the e.m. duality conditions (87).
By recalling (94) and the definition (80), the explicit dependence of G on F can be computed to The Hamiltonian depends on the variables − → D , − → B , and is obtained by a Legendre transform with respect to − → E , namely by defining It is worth stressing that the Hamiltonian is perfectly regular for all values of these fields variables: In particular in the limits of weakly and of strongly valued fields: As yielded by the property (54), the BI Lagrangian L BI is the bosonic part of the N = 1 non-linear Goldstino Lagrangian for the aforementioned partial breaking N = 2 → 1 of rigid supersymmetry; in this context, the Maxwell field is the N = 1 superpartner of the spin-1 2 Goldstino field [28,29,11]. For various recent developments on BI theory, supersymmetry, and e.m. duality, see e.g. [30,31,32,33,34].

Supersymmetry, Special Geometry, and Multi-Field Extension
The supersymmetric BI Lagrangian follows from an algebraic constraint in Superspace among various multiplets [35,36].
This constraint enforces a non-linear relation between the two N = 1 supermultiplets X and W α that build the N = 2 vector multiplet [11,37,38]: where m 1 corresponds to a particular choice of a magnetic FI term. The non-linear relation (105) determines X as a non-linear function of D 2 W 2 and D 2 W 2 and implies the nilpotency constraints The supersymmetric BI Lagrangian then acquires the form where e 1 and e 2 are electric FI terms whose labels correspond to the first two directions in an SU(2) triplet. For canonically normalized vectors, one can then recognize that the scale µ and the theta-angle ϑ can be expressed in terms of m 1 , e 1 and e 2 as µ = √ m 1 e 2 , ϑ = e 1 e 2 . (108) In the bosonic limit, the quadratic constraint becomes where F + is the self-dual curvature of the Maxwell field : and F is an auxiliary complex scalar field such that In the multi-field case, the generalization of L SBI rests on the constraints (A = 1, ..., n) [21,22] and the complete Lagrangian reads The θ 2 component of eq. (112), where F A + are self-dual field strength combinations and F A is the auxiliary field, is the multi-field generalization of the BI constraint (109) that is induced by the Special Geometry. Moreover, it holds that d ABC is a rank-3 completely symmetric tensor, encoding the coefficients of the cubic term of the holomorphic prepotential of the aforementioned rigid special Kähler geometry characterizing the N = 2 vector multiplets' scalar manifold, whose classification is based on the singularity structure of cubic varieties [22,23].
The real parts of Eqs. (114) are n quadratic equations that are generally coupled. Letting they take the form On the other hand, the imaginary parts of Eqs. (114) are n linear equations for ℑF A : The bosonic part of the multi-field Lagrangian (113) can be expressed in terms of the F A and of the real magnetic changes m A . It can also be written in terms of additional complex charges The matrix C AB is needed whenever the matrix d AB := d ABC q C is not positive definite.
Moreover, by a change of symplectic basis one could also eliminate the real parts e 1A of the electric charges, which multiply combinations of the fields that are total derivatives. These Lagrangians combine, in general, a quadratic Maxwell-like term with additional higher-order contributions.
For n = 1, or whenever the matrix d AB is positive definite, one is not compelled to introduce the C AB and the Lagrangian takes the simpler form (111). In all cases, however, the difficult step in the construction of the Lagrangians is the solution of the quadratic constraints, and in particular of the non-linear ones given in (117).

Schrödinger Formulation
The general study and classification of constitutive relations admitting e.m. duality symmetry can be attacked by adopting complex and chiral variables, indeed these variables simply transform with a phase under duality symmetry rotations. We therefore define 2 : so that the e.m. duality condition (87) equivalently reads The corresponding chiral combinations read where F ± = 1 2 (F ± i ⋆ F ), G ± = 1 2 (G ± i ⋆ G). The equations of motion of linear electromagnetism in these variables read The constitutive relations in these new variables, studied in [39], and independently and more generally in [40,41], read as follows: These are six relations if the arbitrary function I (the "Lagrangian" of these new variables) is real.
Moreover, duality invariance is obtained if I is invariant under duality rotations. This is the case We note that the function I is U(1)-duality-invariant because the Lorentz scalar (T − ) 2 (T − ) 2 is duality invariant, indeed under a U(1) duality rotation of angle α, the fields T and T transform with opposite phases e −iα and e iα .
By use of a Legendre transform [39,42,25]  Furthermore, it was shown by Schrödinger [43] (see also [44]) that the constitutive relations for the BI theory are In the general case the following constitutive relationsà la Schrödinger define a nonlinear and self-dual theory [42] ⋆ where we recall that u = 2|(T − ) 2 | = |T 2 | + |T ⋆ T | , and where we require lim u→0 f (u) = 0, so to recover Maxwell theory for weak fields. For example the duality conditions (121) follow contracting (127) with ⋆T µν . Contracting (127) with T µν and with T µν it can also be shown that where In the first Ref. of [41] the function t(u) corresponding to the BI Lagrangian density was found through an iterative procedure order by order in u; the first coefficients of its power series expansion were observed to match those of a particular generalized hypergeometric function, yielding to the conclusion that The relation between the constitutive relationsà la Schrödinger and those in terms of the I(u) function follows from the relations between the T variables and the chiral T ± variables. This allows to determine f (u) from I(u) and vice versa. Explicitly, given f (u), one obtains t = t(u) from (131); the function I(u) is then given by

Generalizations
Other non-linear (square-root) BI-like type Lagrangian densities in D dimensions have been suggested in [45,46], starting from alternative D = 4 Goldstino Lagrangians studied in [12,37,47]. Goldstino multiplet in rigid supersymmetry [47,46]. They read and they have the property of being doubly self-dual under Moreover, after a single duality the Lagrangian ends up containing two forms of the same degree, and it can be recast in the following manifestly U(1) e.m. -invariant form [46]: Actions with the full U(n) (where n is the number of Maxwell vector fields) were proposed in [48], but they are currently not available in closed form, even for n = 2.
On the other hand, we can present three relatively simple extensions of the BI Lagrangian to a pair of curvatures F 1 and F 2 whose field equations possess U(1), U(1) × U(1) and SU(2) duality symmetries. All these models solve subsets of the non-linear Gaillard-Zumino constraints (85) and (86) in the case Λ, Σ = 1, 2, i.e. of : The Lagrangian with U(1) invariance reads and was obtained as an application of double duality between two forms of the same degree. This U(1) corresponds to the constraint (85) for Λ = 1 and Σ = 2, namely (138).
The other two examples are simply different complexifications of the two-field system, whose manifest U(1) invariance corresponds to the single constraint of Eq. (86). In particular, letting the case with U(1)×U(1) invariance, where the second U(1) constraint corresponds to the equation reads This is actually a particular case of the U(n, n) dualities of [48] for n = 1, where in the absence of scalars U(1, 1) reduces to its maximal compact subgroup U(1) × U(1). Finally, the case with SU(2) invariance reads where the other two constraints correspond to the SU(2) generators From these examples it is manifest that the standard BI system admits several types of inequivalent complexifications that differ in their quartic couplings.
A fourth option, with maximal U(2) duality, is naturally captured by the construction of [48], but its Lagrangian, as have been mentioned before, is not known in closed form.
All the above non-linear actions can be massive by introducing Green-Schwarz terms [45], i.e.
couplings to another gauge field of the form A Stückelberg mechanism is then generated, and the non-linear mass terms turn out to take the same functional form as the original non-linear curvature terms; in other words, the (D − p − 1)form gauge field has been eaten to give mass to the original p-form gauge field (whose field strength is H p+1 ). The simplest example of this class of Lagrangians is the D = 4 BI action used to make an antisymmetric field B µν massive : starting point is thus the master action where we have introduced a dimensionless parameter g, the counterpart of the parameter k that accompanies the two-form kinetic term. The parameter µ is the BI scale factor, with mass-squared dimension, which sizes the non-linear corrections.
As above, a massive variant of the BI action principle would be obtained by eliminating H after moving to a first-order form where it is unconstrained. However, as we have seen, the additional field is just a standard Stückelberg mode, so that for brevity we can just display the gauge-fixed Proca-like Lagrangian for the massive BI vector, where L BI (g, µ, F ) = µ 2 8 g 2 1 − 1 + 4 In the massless case, the self-duality of the BI Lagrangian would translate into the condition that L BI g , µ , F µν (A) = L BI g ′ = 1 g , µ ′ = µ g 2 , F µν (C) , where C is the dual gauge field. On the other hand, in the presence of the Green-Schwarz term (m = 0), one can eliminate the vector altogether and work, in the dual formulation, solely in terms of the two-form B µν . The self-duality of the massless BI theory then implies that the dual action involves a BI mass term and reads The massless limit can be recovered reintroducing the gauge invariant combination m B + dC before letting m → 0. In this fashion, the limiting Lagrangian describes a massless two-form, dual to a scalar, and a dual massless vector C.
This procedure could be relevant when the system is coupled to N = 2 supergravity, in which case one of the two gravitino would belong to a massive spin-3 2 multiplet.