Massive gravity from a first-quantized perspective

In this work, we investigate the BRST quantization of the massive $\mathcal{N}=4$ supersymmetric spinning particle, with a twofold purpose: exploring different approaches to give mass to spinning particle models and formulating a first-quantized theory for linearized massive gravity on both flat and curved spacetime. Our results suggest that achieving the nilpotency of the BRST charge requires an Einstein spacetime with vanishing cosmological constant as the only viable consistent background. In the course of the analysis, we take the $\mathcal{N}=2$ supersymmetric worldline as an exemplificative model, correctly producing the Proca theory on curved spacetime. Our analysis shows that the associated BRST system uniquely selects the minimal coupling to the background curvature.


Introduction
The theory of General Relativity (GR) has been the subject of modification attempts since the early years following Einstein's formulation.Despite the considerable success of GR in the description of fundamental aspects of our reality, several phenomena persist unexplained, prompting the study of new physics.An intriguing modification of gravitational theories suggests that gravity is propagated by a massive spin 2 particle: a massive graviton.The question of whether the graviton possesses mass remains an unresolved inquiry to date.This prospect, far from being dismissible on the basis of theoretical and experimental arguments [1], would entail profound consequences, which may account for the enduring history of theoretical explorations dedicated to this possibility.The first theory of massive gravity dates back to 1939 when Fierz and Pauli proposed a relativistic action for a massive spin 2 particle on flat spacetime which is now known as Fierz-Pauli (FP) theory [2,3].While appearing deceptively simple in its formulation, the theory already reveals the roots of several issues that have consistently plagued every effort to improve massive gravity beyond the linear level.Specifically, two primary obstacles arise: the vDVZ discontinuity [4,5], namely the different predictions between linear GR and the FP theory in the massless limit, and the presence of a Boulware-Deser (BD) ghost, an extra degree of freedom which inevitably appears at the non-linear level [6].To this date, only a promising non-linear theory of massive gravity free from the DB ghost instability has been formulated, which is known as dRGT theory [7][8][9].The interested reader may refer to excellent reviews for a comprehensive exploration of this subject [10][11][12].This paper represents the first attempt to contribute to the ongoing investigation from an alternative perspective by employing a worldline approach [13].The primary goal is to achieve a first-quantized description of massive gravity.To this end, we exploit the firstquantized models known as O(N ) spinning particles.These are mechanical models with N local supersymmetry on the worldline that have been shown to describe, on Minkowski, a spin s = N 2 particle in first quantization and constitute an alternative to conventional second-quantized field theories in the study of quantum field theory (QFT) [14][15][16][17][18][19][20][21][22][23].In recent years, considerable effort has been invested in exploring the consistent coupling of the spinning particles to more general backgrounds beyond flat spacetime.It was shown in [17] that for N > 2, worldline supersymmetry transformation rules leave the spinning particle action invariant only if the target spacetime is flat.Subsequent progress was made in [24] and in [20], where it was established how to couple the spinning particle to (A)dS spaces and arbitrary conformally flat spaces respectively.A significant advance has been recently achieved realizing that the technique of BRST quantization -originally developed in the context of the path integral quantization of Lagrangian gauge theories [25] -provides a way to investigate different backgrounds rather efficiently.An appealing aspect of BRST quantization relies on the automatic generation of the complete spectrum of spacetime fields essential for the Batalin-Vilkoviski (BV) quantization of the corresponding spacetime quantum field theory [26,27].The BRST quantization of the spinning particle models led to a first-quantized description of Yang-Mills theory from the N = 2 spinning particle [28] and Einstein gravity from the N = 4 spinning particle [29].The latter analysis showed quantum consistency in coupling the model to a more general set of backgrounds, namely to Einstein spaces.Furthermore, this formalism has been extended to include the NS-sector of supergravity in [30].In this work, we intend to apply this procedure to reproduce linearized massive gravity (LMG) on both flat and curved spacetime.As a warm-up, we address the case of the massive N = 2 spinning particle and its connection to the Proca theory.This preliminary investigation serves a dual purpose: introducing the relevant methodologies and exploring the potential impact of the mass as an obstruction to the BRST algebra.One viable approach to incorporate a mass term is the Scherk-Schwarz mechanism [31], which entails a dimensional reduction of the massless model.Recently, another method has been developed consisting of the introduction of auxiliary oscillators describing the internal degrees of freedom of the spinning particle [32].In this work, both methods are employed to confer mass to the graviton on a flat background, allowing for an in-depth exploration of the main features of the massive N = 4 supersymmetric worldline theory.Notably, these two approaches appear to yield different outcomes, with the first method providing the correct BV spectrum and field equations of linearized massive gravity.Conversely, the second approach encounters difficulties already at the initial stage.Subsequently, we attempt to couple the model to a generic curved background, and our findings suggest that the nilpotency of the BRST charge requires the background to be Einstein with cosmological constant set to zero.
The paper is organized as follows.In section 2 we outline the main features of the O(N ) massless spinning particle models, providing the basis for the upcoming analysis.Section 3 provides an overview of the BRST quantization procedure, with emphasis on the first-quantized massless graviton both on flat and curved spacetime.In section 4 we discuss how to confer a mass to a spinning particle model.We begin by taking the N = 2 spinning particle as an illustrative example to elucidate the role of the mass in the BRST algebra, showing that the model correctly reproduces the Proca theory on a general background.We then proceed to give a mass to the graviton and check whether the theory correctly reproduces that of Fierz-Pauli.Finally, in section 5 we attempt to couple the model to a generic curved background, establishing the consistency of the BRST quantization exclusively on a subset of viable spacetimes.Conclusions and possible outlooks are presented in section 6.

O(N ) massless spinning particles
The so-called O(N ) spinning particles represent prominent examples where gauge symmetries allow for a manifestly Lorentz covariant formulation.This, in turn, gives rise to a constrained Hamiltonian system with first-class constraints [33].These kinds of constraints, denoted as C α (with α spanning the number of constraints), satisfy upon quantization a graded algebra of the form 1 with some structure functions f αβ γ .The phase space action of such models depends on the particle spacetime coordinates x µ together with N real fermionic superpartners Ξ µ I introduced to describe the spin degrees of freedom.In the case of a flat d-dimensional target spacetime, the action is given by2 where a dot indicates a contraction of the internal indices.The theory described by (2.2) should be regarded as a one-dimensional field theory living on the worldline, which provides a first-quantized description of relativistic particles of spin s.A few remarks regarding its field content are in order.
1 With the following notation for the graded commutator: The canonical coordinates (x µ , p µ , Ξ µ I ) upon quantization are subject to the canonical commutation relations (CCR) 3) The worldline supergravity multiplet in one dimension (e, χ I , a IJ ) contains the einbein e which gauges worldline translations, the N gravitinos χ which gauge the worldline supersymmetry, and the gauge field a for the symmetry which rotates by a phase the worldline fermions and gravitinos, the R-symmetry.In (2.2) they act as Lagrange multipliers for the suitable first-class constraints, which are the Hamiltonian H, the supercharges q I , and the R-symmetry algebra generators J IJ , choosing a specific order for the latter to avoid ambiguities Hamiltonian and supercharges together form the following one-dimensional algebra: These constraints have to be introduced to ensure the mass-shell condition and to eliminate negative norm states, ensuring the consistency of the model with unitarity at the quantum level [23].On the other hand, regarding the aforementioned R-symmetry algebra one finds: The gauging of the R-symmetry group is optional and can be used to constrain the model to deliver pure spin s states and have the minimal amount of degrees of freedom: to count them, one can construct the path integral on the one-dimensional torus of the free spinning particles, which has been achieved for all N in [34].Equations (2.5) together with (2.6) form the so-called O(N )-extended worldline supersymmetry algebra: it displays N supercharges q I which close on the Hamiltonian H and which transform in the vector representation of SO(N ), whose Lie algebra is described by the second line of (2.6).In the following, we will refer only to (2.5) as "N SUSY algebra", since it will be appropriate to distinguish the two algebras in light of the future BSRT procedure.It is rather convenient to work with a complex redefinition of the original fermionic variables Ξ I → (ξ i , ξi ), namely for the case of even N , i.e. for particles with integer spin.The respective CCR Accordingly, both the supercharges and the R-symmetry generators split under the complex redefinition: we report the explicit expressions for the particle models of interest, namely the spin one and the spin two.By choosing a specific Fock vacuum, an arbitrary state |Ψ in the Hilbert space is isomorphic to the wavefunction Ψ(x, ξ), on which the conjugated momenta act as derivatives (2.9) Thus, for the case N = 2, the spin one model, one has the splitting q I → (q, q) and J IJ → J as follows where the shift − d 2 in the definition of the R-symmetry generator is a quantum effect stemming from an antisymmetric ordering for the Grassmann variables [35].For the N = 4 case, the spin two, one has q I → (q i , qi ) and the SO(4)-symmetry generators are realized maintaining manifest covariance only under a u(2)-subalgebra J IJ → (J i j , Tr ij , G ij ), explicitly realized as3 .11)Note that the only surviving components of the Tr and G operators are those with i = j.

BRST quantization
There are at least two methods for dealing with first-class constraints and constructing canonical quantization of gauge systems: the Dirac method and the BRST quantization.
In the former approach, the procedure entails turning the constraints into operators Ĉα and imposing that an element |Ψ of the Hilbert space is physical if it gets annihilated by the constraint operators, i.e.Ĉα |Ψ = 0 .
This procedure has been extensively employed in [20,22] to quantize both massless and massive spinning particle models.It is a relatively straightforward and powerful method, however, its application fails when the constraints algebra ceases to be first-class.These considerations eventually led to the BRST analysis of [28][29][30], which is reviewed in the following to delineate the main steps of the quantization procedure, with the intention of applying them to the massive case.
Take the N = 4 spinning particle with a four-dimensional target space as an instructive example.The internal indices take values i = 1, 2 and the N = 4 worldline SUSY algebra is explicitly realized as The initial step of the BRST procedure entails an enlargement of the Hilbert space to realize ghost-antighost pairs of operators (g α , P α ) associated with each constraint C α , with opposite Grassmann parity of the latter -that is, anticommuting ghosts for bosonic constraints and commuting ghosts for fermionic constraints -and canonical graded commutation relation Hence, we assign the fermionic pair (c, b) to the Hamiltonian, and the bosonic superghost pairs (γ i , β i ) and (γ i , βi ) to the supercharges q i and qi respectively, obeying with ghost number assignments gh(b, The second step consists of constructing the BRST operator Q.It is realized as follows where ǫ(g β ) is the Grassmann parity4 of the ghost g β .Equation (3.7) is exact if the structure functions f αβ γ are constant, and can be derived by demanding the following properties for the BRST operator: • It has to be anticommuting and of ghost number +1.
• It has to act on the operators, corresponding to the original phase space variables prior to quantization, as the gauge transformations with the ghost variables replacing the gauge parameters.This, together with the latter requirement, is enough to constraint the structure to be Q = g α C α + . . .
• Finally, the BRST charge has to be nilpotent.This determines the second structure of (3.7) as can be checked by direct computation of {Q, Q}.
Note that the BRST charge is nilpotent by construction as long as the associated algebra is first-class.In more general cases, higher-order terms may appear and need to be determined by the nilpotency condition.Given its significance for future discussions, we shall show it for the case of the N = 4 spinning particle.The BRST operator associated with the first-class algebra (3.2) is To verify its nilpotency, it is first convenient to define ∇ := γ • q + γ • q, such that for any operator of the form (3.8) one always finds the following potential obstruction-terms to its nilpotency which are vanishing in the present case since Notice the critical role played by the presence of the first-class algebra (3.2) in both the cancellation of (3.9) and of (3.10).The same does not hold anymore in more general cases.
In the next section, we will illustrate how to address this issue.
The careful reader might wonder about the absence of a set of ghosts associated with the so(4) constraints.The key idea of the BRST procedure, as first developed in [28] providing a first-quantized description of the gluon, is to treat the R-symmetry constraints and the SUSY ones on different footings.The formers are imposed as constraints on the BRST Hilbert space, thus defining precisely the general dependence of the "string field" Ψ on the spacetime fields content. 5It is within this restricted Hilbert space that the cohomology of the BRST charge has to be studied.The procedure remains consistent as long as the R-symmetry constraints (2.11) commute with the BRST charge: in order to achieve that, it is necessary to extend (J, Tr, G) → (J , Tr, G) as follows The relevant so(4) generators to be imposed as constraints on the BRST Hilbert space are the two number operators J i := J i i (i not summed), namely the diagonal entries of J i j , the Young antisymmetrizer Y := J 2  1 , and finally the trace Tr, which implement the maximal reduction of the model.We collectively denote these constraints as T α := (J , Y, Tr).The BRST system is defined as follows and its consistency is ensured by [Q, T α ] = 0.This is indeed equivalent to saying that we are studying the cohomology of Q on the restricted Hilbert space of fixed R-charge defined by H red := ker T α .The system above will serve as the starting point for BRST quantization in both the massless and massive cases.

Massless graviton on a flat spacetime
The ghost vacuum |0 is chosen such that it is annihilated by (b, γi , βi ) as in [29], so that a general state |Ψ in the BRST extended Hilbert space is isomorphic to the wavefunction Ψ(x, ξ | c, γ, β), on which the antighosts act as derivatives, i.e.
5 Regarding terminology, we deliberately confuse "string field" and "BRST wavefunction" since Ψ plays the same role in the worldline theory as it would in string field theory: its expansion as a linear combination of first-quantized states display coefficients which correspond to ordinary particle fields [36,37].
In the following, we collectively denote the ghost oscillators as g := (c, γ, β).The BRST operator, making explicit the d'Alembertian = η µν ∂ µ ∂ ν while adjusting the coefficients, acts as while, regarding the so(4) generators where we have explicitly expressed the number operators, counting the number of oscillators with a fixed flavor index, and where we solved ambiguities in the definition of J i by choosing a symmetric ordering for the ghosts.
It is useful for future reference to highlight the intermediate steps of the calculation.The first condition imposed by Ψ ∈ H red , namely The latter expression is essentially a Taylor expansion in powers of c taking the form Ψ = χ + χ * c, where both terms χ and χ * of the wavefunction contain oscillators and spacetime-dependent field components.Requiring the entire wavefunction to have fixed Grassmann parity and fixed ghost number forces the field components in χ * to have opposite parities and ghost numbers decreased by one compared to those in χ.It is important to note that we are not yet in the position to correctly identify (anti)fields: we first need to reduce Ψ to include only the BV spectrum of the theory.Imposing the two remaining conditions we obtain together with a µν = a νµ using Y |Ψ = 0, while from Tr |Ψ = 0 we get Now we can indeed interpret Ψ as a spacetime BV string field, that contains the whole minimal BV spectrum of pure gravity along with auxiliary fields: we assign Grassman parity and ghost number to the component fields such that the entire wavefunction Ψ has total even parity and ghost number zero.It becomes explicit with the following identifications for the fields and for the corresponding antifields Grassman parities and ghost numbers can be read from table 1.The most general state in It contains the graviton h µν and its trace h, an auxiliary vector field v µ , and the diffeomorphism ghost ζ µ , while the remaining components are the corresponding antifields.

BV role Field Grassmann parity Ghost number massless graviton
Table 1.List of fields in the physical sector of the massless N = 4 model with the corresponding ghost number and Grassmann parity.
The final step involves evaluating the field equations of the theory, to verify that the model correctly reproduces a first-quantized representation of linearized gravity.This can be accomplished through the BRST closure equation at ghost number zero, which produces the massless free spin two field equation which implies a vanishing linearized Ricci scalar.To conclude this section, one is left to verify the accurate reproduction of the gauge symmetry: this is achieved from the ghost number zero part of where Λ ∈ ker T α contains the gauge parameters of the associated symmetry while having the same functional form as Ψ, with overall odd parity and ghost number −1, i.e.
where the gauge parameter ε µ has even parity and ghost number zero.As expected, the result is For future reference, the closure equation at ghost number one is reported below: Equations (3.28) represent a set of equations of motion for the diffeomorphism ghost ζ µ (x).

Massless graviton on a curved background
The consistent coupling of the spinning particles to more general backgrounds beyond flat spacetime is a rather delicate matter [38].In previous works, a quantization à la Dirac has been extensively employed.However, when the SUSY algebra fails to be first-class this method loses its validity, potentially leading to the misleading conclusion that quantization is not feasible in more general cases.This limitation was evident in the case of the N = 4 spinning particle, for which only certain restricted backgrounds were found to be viable until it was realized that the technique of BRST quantization offers a way to explore more general backgrounds, as recently discussed in [39,40].In this section, we review the exploration of such possibility.
The coupling of the massless N = 4 spinning particle to a curved background with metric g µν (x) is realized by the covariantization of the derivatives, i.e.
where fermions carry flat Lorentz indices so that ξ µ i := e µ a (x) ξ a i , introducing a background vielbein e a µ (x) and the torsion-free spin connection ω µ ab .The covariant derivative operator ∇µ reproduces the effect of the usual covariant (partial) derivative ∇ µ (∂ µ ) on the tensorial (scalar) components contained in the wavefunction (3.21), e.g.
The presence of a general background manifests itself as deformations of the original BRST system: indeed the supercharges become The main consequence is that the SUSY algebra does not close anymore: where the Laplacian is defined as The BRST operator needs to be deformed as well, requiring the construction of an ansatz due to the fact that the associated algebra is no longer first-class.The general approach involves considering the same general form (3.8), but with the inclusion of more general terms accounting for the obstruction to the first-class character of the associated algebra, namely the curvature in the present case.Therefore, possible non-minimal couplings to the curvature, collectively denoted as ℜ, must be incorporated inside the Hamiltonian.These couplings may act as obstructions to the nilpotency of the BRST operator, and the BRST analysis aims to determine which of these terms persist to ensure a nilpotent Q. Avoiding higher powers of ghost momenta and assuming that derivatives are deformed only through minimal coupling, the ansatz for the BRST charge is where D is the deformed Hamiltonian in its operatorial form with κ a coefficient to be determined, and the ∇ operator has been conveniently redefined as Note that only the so(4) constraints T α (3.15) remain unchanged, therefore the wavefunction Ψ(x, ξ | g) is still expressed as in (3.21).
The BRST analysis starts with the general expression for Q 2 , namely the extension of its flat spacetime counterpart (3.9) with adjusted coefficients: which remains valid regardless of the specific spinning particle model coupled to a general curved background.In general, (3.37) includes two independent obstructions, which in this case are given by Recall that the BRST cohomology is defined on the reduced Hilbert space H red , that is, one needs to evaluate the obstructions on ker T α .Therefore, for the sake of the cohomology, it suffices to ensure nilpotency of the BRST charge when acting on the physical sector of the theory: The effect is for instance that any contribution with at least three barred oscillators is set to zero.Equations (3.38)-(3.39)get then reduced to [D, ∇] The conclusion is that the nilpotency of the BRST charge, in the massless case, is achieved only on Einstein manifolds, i.e.
upon setting κ = 1 2 .Let us clarify this crucial point for the upcoming discussion: while the obstruction given in (3.42) vanishes by itself on Einstein spaces, the one in (3.41) is zero only when acting on the physical sector: The closure equation evaluated using the deformed BRST charge yields the correct field equations for a massless7 graviton on Einstein spaces, namely together with The correct gauge symmetry δh µν = 2∇ (µ ε ν) is obtained as well.

Giving mass to the graviton
In this section, we present two methods to confer a mass to a spinning particle model.The first one, the Scherk-Schwarz mechanism, is analogous to the Kaluza-Klein compactification.In the second method, the auxiliary oscillators approach, the model is treated as a truncation of the RNS open superstring [42].We employ both methods to give mass to the graviton on a flat background and subsequently discuss the results.However, before delving into the graviton case, we begin by examining the N = 2 scenario, in order to introduce the techniques and investigate the role of the mass as a potential obstruction to nilpotency of the BRST charge.

N = 2 massive spinning particle and the Proca Theory
The massless N = 2 spinning particle has been thoroughly examined in previous works, leading to a first-quantized description of massless antisymmetric tensor fields of arbitrary rank [14][15][16][17][18].It has been shown how to provide a first-quantized description of the photon both on a flat and on a curved spacetime [35,40].Additionally, the coupling to a nonabelian background has been explored, reproducing the non-linear Yang-Mills equations [28] upon BRST quantization.
The massive scenario has undergone investigation as well: the procedure à la Kaluza-Klein has been implemented in [15,43] where it has been shown how to embed the model in an arbitrary curved spacetime background.This has allowed for the utilization of this model in the context of worldline quantum field theory formalism, particularly in modeling classical scatterings of compact objects in general relativity (see [44,45] and related literature).More recently, the coupling to both an abelian and a non-abelian vector background field has been investigated in [32], where the authors first introduced the auxiliary oscillators approach.

Mass improvement on flat spacetime
The main features of the N = 2 model can be immediately derived from section 2 by setting I = 1, 2. For a comprehensive analysis of the action, the first-class algebra, and the canonical quantization, we direct the reader to [43].In this work, we will provide a concise overview of the essential elements of the BRST quantization of the massive case on flat and, most notably, on curved background as a toy model for the subsequent analysis, which is not covered in the literature mentioned above.
The procedure consists of starting with the massless theory formulated in one dimension higher, and subsequently reducing the dimensionality of the target space through the introduction of suitable constraints [22].We consider the model to live in a flat spacetime of the form M d × S 1 of D = d + 1 dimensions with coordinates and worldline superpartners given by with the index splitting M = (µ, D).The idea is to gauge the compact direction x D , corresponding to S 1 , by imposing the first-class constraint This results in the phase space action (2.2) having a leftover term ∼ ẋD , which can be regarded as a total derivative and thus dropped.The CCR (2.3) are realized explicitly as with the other (anti)commutators being zero.The N = 2 SUSY constraints get modified by the presence of the mass, taking the form while the R-symmetry generator (2.10) becomes Together they still satisfy the same N = 2 supersymmetry algebra despite the mass improvement which successfully remains first-class.Note the key role played by the surviving fermionic coordinates coming from the extra dimension θ := ξ D and θ := ξD , which are responsible for the introduction of the mass term into the theory.Indeed, in the limit where (θ, θ) → 0 the massless theory must be recovered as a consistency check.
At this point, the BRST quantization on a four-dimensional flat spacetime proceeds smoothly, just as described for the N = 4 case.Upon enlarging the Hilbert space to realize the set of ghost operators (c, γ, γ) with relative momenta (b, β, β) as in (3.4), the BRST charge is constructed as usual and the so(2) constraint has to be extended to include ghost contributions as follows with the usual prescription employed to resolve ambiguities.Nilpotency of Q is still guaranteed, despite the mass improvement, since the associated algebra is first-class and the BRST system is consistent since [Q, J ] = 0.The physical sector is defined as the eigenspace of J with a fixed R-charge − 1 2 ,8 i.e. physical states are eigenstates J |Ψ = − 1 2 |Ψ and are isomorphic to wavefunctions Requiring Ψ to be Grassmann-odd and have ghost number zero, it can be interpreted as a spacetime BV string field displaying the complete minimal BV spectrum of the Proca theory along with an auxiliary field.The Grassmann parities and ghost numbers of the components can be found in table 2. In contrast to the massless case, the physical fields include not only the massive spin one A µ but also the Stückelberg scalar ϕ.The inclusion of the Stückelberg scalar becomes necessary to restore the U (1) gauge symmetry [46], which is broken due to the introduction of the mass.Notably, the spectrum also includes the associated scalar ghost Φ in the spectrum.

BV role Field Grassmann parity Ghost number massive spin one
A µ 0 0 Stückelberg scalar ϕ 0 0 auxiliary scalar φ 0 0 scalar ghost Φ 1 1 Table 2. List of fields in the physical sector of the massive N = 2 model with the corresponding ghost number and Grassmann parity.
The field equations, upon solving for the auxiliary field, are which are indeed the Proca field equations with the scalar ϕ(x) playing the role of Stückelberg field, while the gauge symmetry, from δΨ = QΛ with is where Σ(x) is a local gauge parameter, of even parity and with ghost number zero.The same conclusions hold when employing the auxiliary oscillators approach, although it is not explicitly shown here.The subsequent analysis will elaborate on the limitations of this approach in the context of the graviton case.

Mass obstruction on curved spacetime
We are finally in the position to investigate whether the presence of the mass plays a role when a curved background is considered, generalizing the massless BRST analysis of [40].The coupling of the massive N = 2 spinning particle to gravity is realized by the covariantization of the reduced model just discussed, with the mass already present in the theory.The supercharges are then deformed as follows ∐ := −i ψ a e µ a ∇µ + mθ , ∐ := −i ψa e µ a ∇µ + m θ .where the vanishing of the top-right commutators is guaranteed by the cyclic identity for the Riemann tensor.The corresponding BRST operator is defined as usual Q = c D + ∇ + γγ b, with potential obstructions incorporated into the definition of ∇ with respect to the massless case: The existence of an associated first-class algebra (4.17) ensures the nilpotency of the BRST operator without any further conditions on the background metric, 9 which is expected from QFT considerations.While this has been already established for the massless case, the same holds even in the massive one, as the mass does not affect the algebra, as we now shall work out explicitly.Starting from the general expression one finds the following independent obstructions In the present case, evaluating the BRST cohomology on the reduced Hilbert space, ker J , has the effect of setting to zero any obstruction of the form O AB ZA ZB with O AB arbitrary operators since an arbitrary state in the physical sector has the form This is enough to conclude that indeed without further conditions on the background; in other words, the massive N = 2 spinning particle can be coupled to off-shell gravity.Remarkably, the mass plays no role in the BRST algebra, as it is trivially canceled in (4.20) and commutes with any operator inside (4.21) being constant.It is worth noting that introducing a term proportional to the Ricci curvature inside the Hamiltonian (4.16) D → D + κR, would force κ to be zero (either that or a vanishing Ricci scalar R = 0).This is dictated by the BRST algebra (4.20)-(4.21),as any term of that form is incompatible with the closure of the N = 2 SUSY algebra (4.6).Interestingly, the N = 2 spinning particle appears to select only the minimal coupling to the background, as we shall check from the equations of motion.From the closure equation QΨ = 0, one finds which, upon eliminating the auxiliary field, yield Equations (4.27)-(4.28)represent the field equations of the Proca theory on a general curved spacetime in its Stückelberg formulation.In particular, the N = 2 spinning particle provides a first-quantized version of the minimal extension for a theory of a free massive vector field in curved spacetime, without the inclusion of possible non-minimal couplings to the background in the corresponding spacetime quantum field theory [47,48].Their absence is attributed to the reasons previously discussed.As a final note, the gauge symmetry, using (4.13), reads

N = 4 massive spinning particle and Linearized Massive Gravity
In the remainder of this section, the methods previously discussed are exploited to give a mass to the graviton, considering a flat target spacetime for the time being.

Dimensional reduction approach
The first method proceeds along the lines of what has been shown for the N = 2 case, namely through the reduction of the higher-dimensional massless model.The fermionic coordinates carry a flavor index The same applies to the SUSY constraints (2.4) that satisfy the same first-class algebra despite the mass improvement.Regarding the so(4) symmetry constraints, equations (2.11) become At this point, the BRST quantization proceeds as outlined in section 3, with the BRST system defined by (3.12).The BRST operator (3.14) takes the usual form when acting on wavefunctions Ψ(x, ψ, θ | g) and is still nilpotent, for the same reasons outlined in the previous section.The relevant so( 4) generators (3.15) to be imposed on the extended BRST Hilbert space are extended to commute with Q, including the extra fermionic oscillators The physical subspace has a fixed U (1)×U (1) charge corresponding to 3−d 2 .Let us comment on the fact that only in three spacetime dimensions the wavefunction Ψ ∈ ker J i , which becomes relevant when dealing with the construction of the worldline path integral. 10In four spacetime dimensions, physical states are charged − 1 2 , and we shall consider the shifted operator J = J + 1 2 for the sake of notational simplicity.To emphasize the impact of the introduction of a mass term, we follow the same steps as in the massless case to impose the condition Ψ ∈ H red , with the reduced Hilbert space defined as the kernel of Ji and Y, Tr (4.36), still collectively denoted as T α .From Ji |Ψ = 0 descends the general form of the wavefunction , where in the first line Ψ| m=0 denotes the massless wavefunction (3.16).The condition Y |Ψ = 0 produces, in addition to the "massless" contribution (3.17), while the remaining Tr |Ψ = 0 produces It is now possible to identify the BV spectrum of the theory, which becomes evident after renaming the field components as in (3.19) along with and, for the corresponding antifields, Requiring Ψ to be Grassmann-even and have ghost number zero the field content gets assigned with the corresponding Grassman parities and ghost numbers as reported in table 3. The most general string field Ψ in ker T α reads where, for the sake of simplicity, the antifields content is left implicit.

BV role Field Grassmann parity Ghost number massless graviton
Table 3. List of fields in the physical sector of the massive N = 4 model with the corresponding ghost number and Grassmann parity in the dimensional reduction approach.
A few comments are in order.Firstly, it can be verified that the first two lines of (4.42) correspond to the field content in the massless case (3.21).This is expected since, as previously suggested, up to an obvious redefinition of the scalar field ϕ(x) which can be absorbed into the trace h(x).Taking into account also the other two lines, it is clear that Ψ contains the whole minimal BV spectrum of LMG, with the inclusion, with respect to the massless spectrum, of the two Stückelberg fields A µ and ϕ, a scalar ghost field Φ and an auxiliary scalar field φ [50].At this point, we proceed to analyze the field equations, which are derived from the closure equation at ghost number zero using the BRST operator (4.35).The result is It is possible to solve (4.44e) and (4.44f) for the auxiliary fields v µ and φ respectively, leading to the following equations which can be further simplified combining (4.45b) and (4.45d) to reach the following set of equations Equations (4.46)-(4.48)correspond to the field equations of linearized massive gravity in the Stückelberg formalism, where, as anticipated, ϕ(x) is the Stückelberg field and A µ (x) is its vector counterpart.
A few comments are in order.
• Setting both Stückelberg fields to zero, which in literature is known as the unitary or physical gauge, equations (4.45) can be reduced to the Fierz-Pauli system i.e. the field equations for the theory of a massive spin two field in which the mass term explicitly breaks the gauge invariance.
• At first sight, the massless m → 0 limit is rather peculiar.Instead of resulting in the free and massless spin two field equations (3.23), it produces the following outcome Equation (4.51) represents the field equation for a free-propagating vector field A µ (x), describing a spin one massless particle that becomes decoupled in this limit.On the other hand, (4.52) describes the propagation of a massless scalar field ϕ(x), which, surprisingly enough, is still coupled to the wanna-be free massless graviton field h µν (x) in (4.50).This hints at the vDVZ discontinuity, a peculiarity associated with the Fierz-Pauli formulation of massive gravity.Notably, there exist other LMG formulations in which this discontinuity is absent, as discussed in [51][52][53] and related literature.The origin of the discontinuity is traced back to the fact that the massless limit of a Fierz-Pauli graviton is not a massless graviton, but rather a massless graviton plus a coupled scalar. 11Remarkably, our model captures this characteristic feature of massive gravity.
Let us finalize the analysis by investigating the gauge symmetries of the theory.
Any theory of massive gravity -particularly the Fierz-Pauli theory -is not gauge invariant due to the presence of the mass term, unless redundant degrees of freedom are introduced, as is the case for the Stückelberg trick.To explore this fact from a firstquantized perspective, it is sufficient to consider the closure equation at ghost number one: while in the massless case, the result is a set of dynamical equations for the diffeomorphism ghost (3.28), in the massive scenario things become more intricate due to the presence of another ghost field Φ(x).The outcome manifests the following additional equations: This indicates that both the scalar and the diffeomorphism ghosts are set to zero, which accounts for the breaking of the gauge invariance due to the presence of a mass term.
To display said symmetries, one has to compute the ghost number zero component of δΨ = Q Λ, where Λ is the massive extensions of (3.26), namely where ε(x) and Σ(x) are the two gauge parameters associated with the two gauge symmetries.The outcome is ) as expected from the Stückelberg formulation.Therefore, the massive N = 4 spinning particle produces a gauge invariant formulation of LMG.

Auxiliary oscillators approach
The main idea of the auxiliary oscillators approach consists of enlarging the BRST algebra with additional variables, whose physical interpretation can be seen as to describe the internal degrees of freedom of the N = 4 spinning particle.Thus, the first step involves an extension of the phase space by introducing two canonically conjugated complex bosonic variables, denoted as (α a , ᾱa ), with commutation relation Note that the index a may run over an arbitrary set, including a single value.To preserve the N = 4 worldline supersymmetry, the introduction of four complex fermionic variables (η a i , ηi a ) is necessary.They represent the superpartners of the α's, with anticommutator The next step involves finding a way to deform the SUSY constraints while ensuring the closure of the algebra.A straightforward attempt would be to extend the N = 2 constraints of [32], resulting in (4.59) However, this leads to q i , qj = 2δ j i H, specifically and there is not an obvious redefinition of the Hamiltonian that allows for a closure.Consequently, the resulting algebra is not first-class, and potential issues are anticipated in the quantization process.Indeed, the associated BRST charge Q = c H +γ • q + γ •q −2γ •γ b fails to be nilpotent: The auxiliary oscillators approach may be a viable option, although it is necessary to follow the principles outlined in section 3.2.One should engineer a deformation of the BRST algebra that leads to preventing the obstruction given by (4.61), e.g. through a suitable redefinition of the Hamiltonian and of the so(4) constraints.Further work is needed to explore this avenue, but such investigations are beyond the scope of the current discussion, as a promising method for dealing with a massive graviton is already at hand.

Massive graviton on curved spacetimes
In this section, the investigation of a first-quantized massive graviton on a curved spacetime is carried out.Previous works have addressed the challenge of coupling massive higher spin fields -which, in the present formalism means for N > 2 -but the analysis failed in going beyond (A)dS spaces [22].The goal of the present work is to overcome said result with the N = 4 spinning particle.Let us rephrase the objective more clearly: the aim is to provide a worldline formulation of the linear theory of massive gravity, namely of the Fierz-Pauli theory, on a curved spacetime; such a theory describes the propagation of massive spin 2 particle on a non-flat background.
The coupling to a generic curved background is achieved through the covariantization of the reduced model previously discussed, which results only in a mild modification of the massless model of section (3.2).The supercharges are deformed as and together with the Laplacian ∇ 2 form the same obstructed algebra (3.32), with the sole modification arising from the following anticommutator (5. 2) The deformed BRST charge displays still the general form Q = c D + ∇ + γ • γ b, where now The model is considered on Einstein spaces: indeed, it is expected for the massive theory to reproduce the correct results in the m → 0 limit -similarly to the flat spacetime scenario -and the massless BRST charge is nilpotent only when R µν = 1 4 g µν R. It is not possible to obtain a weaker condition for the background in the massive case.Moreover, from a QFT perspective, it is known that an Einstein space is the only space on which a free (massive) graviton can consistently propagate [55][56][57][58].Hence, it is natural to implement this condition right at the beginning of the analysis.
The so(4) constraints T α are given in (4.36), thus the wavefunction Ψ can be read from (4.42).This section aims to investigate if and how the presence of the mass affects the nilpotency of the BRST operator.In contrast to the case of the massive spin one particle on a general background, where nilpotency was not a concern even in the massless case, the situation is different here.The massless N = 4 spinning particle already selects only a specific subset of available backgrounds.The starting point remains and thus, it is necessary to address the independent obstructions Just as in the case of the Proca theory, it is evident that the mass does not enter into the BRST algebra.Specifically, note that the term in (5.6) is canceled both from a symmetrybased argument and by the fact that it annihilates any states Ψ in the physical sector, thereby being set to zero.Additionally, no mass term survives inside the commutator in (5.7).Consequently, the mass does not explicitly obstruct the nilpotency.The expressions remain the same as in the massless case: once evaluated on the reduced Hilbert space H red , taking Einstein manifold simplifications into account, they read

[D, ∇]
ker Tα = 0 . (5.9) The leftover operator (5.8) is vanishing when acting on the massless physical wavefunction Q 2 Ψ| m=0 = 0.However, this is not the case for the massive graviton, as the wavefunction has changed.When acting on the physical sector, the outcome is Let us highlight that attempts to implement possible couplings involving traces of the Riemann tensor and the newly introduced fermionic variables inside the Hamiltonian (5.3), such as do not solve the issue.Although an appropriate tuning of the coefficients -in particular one finds c 1 + c 2 = −2 -potentially addresses the first obstruction (5.8), these terms would disrupt the second one (5.9),compelling c 1 = c 2 = c 3 = 0.This explains why those terms have been omitted right from the beginning: as already discussed in [29], one cannot have fermions inside ∆R.The only possibility to achieve nilpotency of the BRST operator is to consider Ricci-flat backgrounds, i.e.
R µν (x) = 0 , (5.12) which immediately addresses the remaining obstruction (5.8), yielding a nilpotent BRST charge.This outcome resembles the scenario encountered in string theory, where consistency of the string propagation, worldsheet conformal invariance, implies Ricci flatness to leading order [59].Such an unfortunate result amounts to the difficulties that arise when trying to give a mass to a graviton and trying to generalize it to a curved spacetime.As discussed in the introduction, several obstacles have been encountered on the quantum field theory side, with a promising solution emerging only in recent years.It is plausible that a similar intricate path has to be navigated even within the worldline formalism.The present work may then serve as an initial step in the direction of future developments.
In conclusion, the field equations for a massive graviton on an Einstein background with zero cosmological constant are evaluated using the deformed BRST operator incorporating the appropriate curvature terms.Upon eliminating auxiliary fields, the resulting equations of motion are which can be simplified as Equations (5.14)-(5.16)correspond to the field equations governing linearized massive gravity on a Ricci-flat background within its Stückelberg formulation.The same equations of motion can be derived from a field theory approach: one has to take into account the most general action functional while keeping only a restricted class of non-minimal couplings to the background, namely, the ones that do not excite possible unphysical degrees of freedom [57].The final step involves performing the necessary Stückelberg tricks to restore gauge invariance.As a final remark, it is worth noting that the correct gauge symmetries can be derived from Λ (4.54), yielding

Conclusions
In this work, we have investigated the issue of consistently coupling massive higher-spin fields to a general curved background within the wordline formalism, with our primary focus being that of a first-quantized formulation of linearized massive gravity.During this investigation, we have explicitly shown the BRST quantization of the massive N = 2 spinning particle coupled to off-shell gravity, which, to the best of our knowledge, has never been investigated before in the literature with these tools.Our results show that the model reproduces the Proca theory on curved spacetime with the specific selection of the minimal coupling to the background, improving knowledge built on previous results and offering potential utility in the context of worldline computations employing this model.Subsequently, we have addressed directly the issue of a first-quantized massive graviton, discussing the correct reproduction of the Fierz-Pauli theory on a flat spacetime through the dimensional reduction of the higher-dimensional massless N = 4 spinning particle.Additionally, we addressed the challenges associated with the auxiliary oscillators approach.Finally, we tried to couple the massive particle to an Einstein background.Our findings suggest that a Ricci-flat spacetime emerges as the only available background for consistency at the quantum level.Let us stress that the theory constructed here is the worldline counterpart of the linear theory of massive gravity, namely the Fierz-Pauli theory from a QFT perspective: section 5 has been devoted to the correct reproduction, starting from the BRST system of the spinning particle model (3.12), of the field equations and the gauge symmetries of LMG in the Stückelberg formulation.In particular, from the quantum field theory side, the same results can be derived from the FP action put on a curved spacetime, which is the theory for a massive spin 2 propagating on non-flat geometry and which can be obtained as the first approximation of a general massive gravity theory around a given fixed background.Even from the quantum field theory side, the construction of interacting theories of massive gravity is much harder, and a well-defined non-linear completion of the theory has been only recently constructed: achieving such a result in the first quantization formalism -if feasible -would require a considerable effort and is beyond the scope of the present work.
With these results at hand, this work may be regarded as a first step towards the realization of a first-quantized massive graviton on an Einstein spacetime with non-zero cosmological constant, either by identifying consistent improvements in the dimensional-reduced BRST system, or fully pursuing the auxiliary oscillators procedure, contingent on resolving the seemingly simpler yet already intricate flat spacetime scenario first.As a preliminary step in both instances, an intriguing prospect involves investigating whether modifying the model could result in a first-quantized partially massless graviton [60,61].There are possible extensions worth exploring, such as relaxing the constraint on the BRST Hilbert space to subalgebras, which should in principle lead to the inclusion inside the BV spectrum of the N = 0 supergravity, i.e. the particle theory that has in its spectrum the graviton, the dilaton, and the antisymmetric Kalb-Ramond tensor field.Such extensions might be pursued along the lines of [30], investigating possible couplings to the background fields through deformations of the BRST charge.Furthermore, the massive N = 4 worldline model developed here could be exploited in future computations.One potential application involves investigating the one-loop divergences of massive gravity, reproducing, and potentially extending, the results recently obtained from the quantum field theory side [62,63].To pursue this application, it is necessary to establish a suitable one-loop quantization of the massive model and study the corresponding partition function on the one-dimensional torus.The initial step would be to determine the correct measure on the moduli space of the supergravity multiple, implementing the correct projection on the massive gravity contribution, following the analysis of [64].The computation of the heat kernel/effective action coefficients should proceed along the lines of [65], where the counterterms necessary for the renormalization of the one-loop effective action of pure gravity were computed using the massless N = 4 spinning particle.This project is currently underway [49].
ers and the students of the 29th W.E. Heraeus Summer School "Saalburg" for Graduate Students on "Foundations and New Methods in Theoretical Physics" for generating an inspiring scientific environment, which led to meaningful discussions.I also acknowledge hospitality during the final stage of the work by the LMU Munich.