Metric-affine bumblebee gravity: classical aspects

We consider the metric-affine formulation of bumblebee gravity, derive the field equations, and show that the connection can be written as Levi-Civita of a disformally related metric in which the bumblebee field determines the disformal part. As a consequence, the bumblebee field gets coupled to all the other matter fields present in the theory, potentially leading to nontrivial phenomenological effects. To explore this issue we compute the weak-field, post-Minkowskian limit and study the resulting effective theory. In this scenario, we couple scalar and spinorial matter to the effective metric, and then we explore the physical properties of the VEV of the bumblebee field, focusing mainly on the dispersion relations and the stability of the resulting effective theory.


Introduction
The consistent inclusion of Lorentz symmetry breaking in a curved space is certainly one of the most important open problems within studies of this phenomenological idea. Traditionally, the Lorentz symmetry breaking has been introduced in two manners, 1) the explicit one, where a constant vector (tensor) is introduced in the theory from the very beginning, and 2) the spontaneous one, where this constant vector (tensor) arises as the vacuum expectation value of some dynamical field (see [1] for a detailed discussion of approaches to Lorentz symmetry breaking in gravity). a e-mail: adria.delhom@uv.es b e-mail: jroberto@fisica.ufpb.br c e-mail: gonzalo.olmo@uv.es d e-mail: petrov@fisica.ufpb.br e e-mail: pporfirio@fisica.ufpb.br It should be emphasized that while the explicit approach is very convenient in the case of linearized gravity, there are essential difficulties with its application in a full-fledged scenario. Indeed, unlike in flat space, in curved space it is highly problematic to define constant vectors (or tensors) because the condition of vanishing covariant derivatives of a given tensor imposes constraints on the background metric which are difficult or even impossible to satisfy. Otherwise, quantum corrections in curved-space extensions of known Lorentz-breaking field theory models will involve an infinite tower of new terms proportional to covariant derivatives of "constant" Lorentz-breaking tensors, which makes the calculations much more complicated (see e.g. [2]). Another difficulty related with this approach is that in a curved spacetime, the group of general covariant transformations plays a double role, being not only the extension of the Lorentz group but also the gauge group. As a result, introducing terms that explicitly break Lorentz symmetry will imply violation of the gauge symmetry as well, which makes the spontaneous symmetry breaking approach much more appropriate (that see discussion in [3]).
In a curved space-time, the most convenient form of introducing spontaneous Lorentz symmetry breaking is the bumblebee model, first proposed in [1], where the breaking mechanism is implemented via a new dynamical vector field with a nontrivial potential characterized by a continuous set of minima. The first studies of modifications of known gravitational solutions within the bumblebee gravity have been carried out in [4]. Various aspects of the bumblebee gravity have been studied in numerous papers, e.g. [5][6][7][8], all of which assume that the underlying geometry is of (pseudo)-Riemannian type. However, in the absence of any empirical evidence supporting that the space-time struc-ture is necessarily Riemannian, or otherwise 1 , in the high energy regime [9][10][11][12][13][14][15][16], it is legitimate to explore other alternatives. In particular, it is well-known that the metricaffine (or Palatini) formulation of gravity theories beyond GR yields field equations which are inequivalent to those obtained in the purely metric approach, and it has been argued that quantum effects could be encoded at low energies by a non-Riemannian counterpart of the space-time geometry [17][18][19][20]. This fact motivates us to initiate the exploration of Lorentz breaking phenomenology of bumblebee gravity in its metric-affine formulation, in which the connection is treated as a geometrical object a priori independent of the metric. In this sense, it is worth mentioning that Lorentzbreaking extensions of both the Standard Model and General Relativity have already been tested experimentally, leading to stringent constraints on Lorentz-breaking parameters, see e.g. [21? -26].
It has been recently shown that some metric-affine theories of gravity beyond GR have a peculiar behavior that makes them depart from their metric counterpart and become potentially testable via elementary particle interactions. In the particular case of minimally coupled Ricci-Based Gravity theories (RBG's), those in which the gravity Lagrangian is a (projective invariant) function of the metric and the Ricci tensor, one finds that the space-time metric picks up two types of contributions, one coming from the integration over the matter sources and which is mainly responsible for the space-time curvature, and another coming from the energymomentum density of the local sources. The latter contribution is responsible for the existence of a nonzero nonmetricity tensor, Q µαβ = ∇ µ g αβ , in regions where the stressenergy of the matter fields is not covariantly constant. The origin of these contributions in the metric can be traced back to the fact that these theories admit an Einstein frame (via a non-conformal transformation) in which the matter sector features new non-linear interactions and the gravity sector is described by GR. Thus, from the original frame of the theory (RBG frame), these non-linear interactions appear encoded in the space-time metric and are responsible for the non-vanishing non-metricity tensor. The existence of the Einstein frame representation for RBG theories was devised in [28,29], and it has been explicitly proven and used for applications for different matter sectors: perfect and anisotropic fluids, scalar fields and non-linear electrodynamics [30][31][32][33][34]. The deviation from the Riemannian condition ∇ µ g αβ = 0 in the RBG frame is intimately related to a departure of the space-time metric from its Minkowskian form whenever matter fields are present even if the effects of curvature (i.e. Newtonian and post-Newtonian corrections) are negligible. As a result, particle physics experiments can be 1 Here Riemannian means that the connection is completely specified by the metric, as opposed to being an a priori independent geometrical entity.
used to place strong constraints on the parameters of those gravity models [35,36].
In this work we will consider a so far unexplored route, in which the bumblebee model plays a key role. In this model the gravity Lagrangian is represented by the Einstein-Palatini GR Lagrangian and an additional (projectively invariant) nonminimal coupling term between the bumblebee field and the affine Ricci tensor. This sort of non-minimal coupling has not yet been considered in detail within metric-affine theories. In addition, keeping the gravitational part exactly as in GR avoids the generation of non-metricity induced by the stress-energy tensor of the matter fields, which allows to disentangle the different contributions that the non-metricity may have. As we will see, the resulting theory admits an exact formal solution for the independent connection which leads to the emergence of a non-metricity tensor generated by the non-minimal bumblebee coupling. The bumblebee field will then become coupled to the rest of the matter fields present in the theory due to its non-minimal coupling with the gravitational sector. Thus, even though the theory is initially formulated following the postulates of metric theories of gravity so as to satisfy the Einstein equivalence principle [37], in the end the resulting coupling of the bumblebee field to the matter sector implies a direct violation of this principle, as one would expect in a Lorentz violating theory. Besides deriving the field equations and discussing their resolution, in this paper we will also discuss the weak field, post-Minkowskian limit with a focus on the modified dispersion relations and stability of scalar and Dirac fields.
The structure of our paper is as follows. In section 2, we present a general family of Ricci-Based Gravity theories (RBGs) with non-minimal couplings between matter and the connection, introducing the metric-affine bumblebee model as a specific implementation. After discussing the field equations of the bumblebee model, in section 3, we study the weak-field, post-Minkowskian limit and the effective dynamics when curvature effects are negligible as compared to the contribution of non-metricity. Section 4 is devoted to the derivation and discussion of the modified dispersion relations and stability corresponding to scalar and spinor fields. We conclude with a summary and discussion of the results. The analysis of metric-affine RBG theories has focused so far on modifications of the gravitational sector but keeping a minimally coupled matter sector, i.e. with the matter fields interacting only with the metric, not with the connection. Such theories can be described by an action of the form where κ 2 = 8πG, Ψ i denotes a collection of minimally coupled matter fields, and R (µν) (Γ )) represents the symmetrized Ricci tensor of the independent connection Γ α µν . We use units h = c = 1. A straightforward way of generalizing the above action is to allow for a non-minimal coupling between some matter fields Φ i and the connection via R (µν) (Γ )), such that (1) turns into where the Φ i have minimal kinetic terms in the sense of [38].
The detail that only the symmetrized part of the Ricci tensor enters the action guarantees the existence of a projective symmetry, which turns out to be relevant for the stability of the theory by ensuring the absence of ghost-like degrees of freedom in RBG theories [39,40] as well as in metric-affine scalar-tensor theories [41].
In the metric-affine framework, metric and connection are treated as independent fundamental fields. Accordingly, their field equations are obtained by extremizing the action without imposing any a priori relation between δ g µν and δΓ α µν . By varying upon the above action, one finds that the metric and connection field equations of the RBG action with non-minimally coupled matter fields are formally identical to those of the minimally coupled version, namely where T α µν = 2Γ α [µν] is the torsion tensor, T M µν is the usual stress-energy tensor of the minimally-coupled matter sector . Following [29], it is possible to show that the torsion terms on the right-hand side of (4) can be eliminated by exploiting the projective symmetry of the theory. After doing this, one finds that, up to a non-physical projective mode, the solution for the connection is given by the Levi-Civita connection of the metric h µν , namely We then see that the inclusion of matter fields that are nonminimally coupled through the Ricci tensor is still compatible with the fact that in RBG theories the connection is the Levi-Civita connection of some metric, which greatly simplifies the process of solving the field equations.
2.2 Building the Einstein frame of generalized RBGs with matter couplings to R µν .
The fact that the connection can be solved as the Levi-Civitta connection of a metric h µν raises the question of whether these generalised RBGs also admit an Einstein frame representation. As we show next, the answer is positive and the Einstein frame can be constructed with the standard procedure (see e.g. [29]) of linearising the action by introducing an auxiliary (symmetric) tensor field Σ µν . Let us sketch how the Einstein-frame is reached for a generalised RBG of the form (2). Consider the actioñ where f = F(g µν , Σ µν , Φ i ) indicates that we have replaced all instances of R (µν) by Σ µν . The field equations of Σ µν are algebraic, and imply that Σ µν = R (µν) provided that δ 2 f δ Σ 2 = 0. Hence, integrating out Σ µν we see that the new action (6) is physically equivalent to the original one (2).
Let us now define the inverse metric h µν by Generally, we can see the above equation as an implicit definition of Σ µν , which should be seen as a function of (g µν , h µν , Φ i ). Thus, by solving algebraically in Σ µν , we would be able to write Σ µν (g µν , h µν , Φ i ). We can then perform a field redefinition in (6) by substituting Σ µν (g µν , h µν , Φ i ) to arrive atS where . Now, we see that the field equations for g µν are also algebraic and, therefore, we can formally solve them to obtain g µν (h µν , Φ i ,Ψ i ). Thus, integrating g µν out by means of this solution we arrive to the Einstein-Hilbert representation of (2), given by We can clearly see here that the above action (9), which by construction is equivalent to 2, formally describes a set of Einstein-like equations for the metric h µν coupled to a matter sector described byS m which will feature new interactions between the Φ i and Ψ i sectors.

2.3
The metric-affine bumblebee model as a particular case.
As a particular case of the class of theories discussed above, we find the metric-affine version of the curved space-time bumblebee model (see [8]), which is defined by an action of the form A key feature of this model is that the bumblebee field B µ has a non-zero vacuum expectation value (VEV) that spontaneously breaks Lorentz symmetry by introducing a privileged space-time direction. This field is coupled non-minimally to the space-time geometry via the B µ B ν R µν term. The parameter ξ characterizes the strength of this non-minimal coupling between the bumblebee field and the affine connection through R µν (Γ ), and we use the minimal coupling prescription for its kinetic term, so that the field strength of the bumblebee field is defined as the exterior derivative of which would allow us to write the Einstein-frame action corresponding to 10 once the metric g µν is solved in terms of the matter fields (B µ ,Ψ i ).
By taking variations of Eq. (10) with respect to the metric, the connection, and the bumblebee field, we obtain the following field equations: where the prime in V ′ denotes derivative of V with respect to its argument. Here T µν is given by 2.4 Solving the connection equation one can follow [29] to show that Eq.(15) is equivalent to Eq.(4) up to an irrelevant projective mode. Thus, as explained above, the connection can be written as the Levi-Civita connection of h µν . In order to find the explicit relation between the space-time metric g µν and the metric h µν , let us rewrite (18) in matrix form as where the hat denotes a matrix, such thatĥ −1 andĥ are the matrix representations of h µν and h µν , respectively. Taking the determinant of (19) we find that h = g det (I + ξ BB) (18), and plugging this back into (18) we arrive at By inverting the above relation, we also have that Using that det (I + ξ BB) = 1 + ξ X, with X ≡ B µ B µ , we can finally write From this last result, one finds that which provides an algebraic relation between g µν with h µν and B µ , though the scalar X ≡ g µν B µ B ν still contains an explicit dependence on g µν . This dependence can be eliminated by noting that Y ≡ h µν B µ B ν = X 1 + ξ X, which completely solves the problem. Hence, we can integrate g µν out of the action by performing a field-redefinition in terms of h µν and B µ . This is useful in order to physically interpret the different elements that contribute to the space-time metric g µν . In fact, using 24, we can write g µν in terms of h µν and B µ in U B and S m , thus finding the Einstein-frame representation of the action (10), given bỹ with where any barred tensor indicates that its indices are raised with h µν .
In this new (Einstein frame) representation it becomes apparent that the bumblebee model can be interpreted as GR coupled to a modified matter sector in which all the matter fields couple to the bumblebee, which also presents new self-interactions encoded in theV term. According to this, the metric h µν satisfies the Einstein equations coupled to a highly non-linear matter sector. This means that h µν will depart from the Minkowski metric only in regions where the Newtonian and post-Newtonian effects are expected to be relevant, i.e. regions with a strong gravitational field. As a result, as it follows from (24), the metric g µν will not only describe the two propagating degrees of freedom of the gravitational field through h µν , but it will also encode information on the local value of the bumblebee field via a conformal factor and a disformal term proportional to B µ B ν .
From the decomposition (24) and the fact that the independent connection satisfies ∇ Γ α h µν = 0, it follows that the non-metricity tensor Q α µν = ∇ Γ α g µν is non-trivial and entirely due to the derivatives of the bumblebee field. Since this field is expected to have a non-trivial VEV that spontaneously breaks Lorentz invariance, this is an example of a gravitationally generated non-metricity tensor that can develop a VEV. In contrast, in RBGs with minimally coupled matter, the non-metricity is associated to derivatives of the stress-energy tensor of the matter fields, which vanish in vacuum. A constant background of non-metricity was assumed in [42], and experimental constraints to all its possible effective couplings to fermions and photons were derived from Lorentz violation searches in Earth laboratories. Since minimally coupled matter fields do not couple explicitly to non-metricity 2 , these constraints do not apply to our model. However, we note that constraints on Lorentzviolating couplings such as those in the Standard Model Extension [43] could translate into constraints on the bumblebee non-minimal coupling ξ . Further work in this direction is currently in progress.

Metric and bumblebee equations
Let us now continue with the exploration of the field equations (14) and (16). By taking the trace of (14), a relation between the scalar curvature and the trace of the stress-energy tensor T = g µν T µν can be found in the form which exactly matches the relation in GR. By contracting with one and two B µ fields, we can also find the relations respectively, where X ≡ g µν B µ B ν . Substituting the second of the above equations into the first one, we finally obtain The above results can be plugged back into the metric field equations (14) to obtain an expression for R µν (h) which only involves T µν , B µ and the metric g αβ which, recall, can be explicitly written in terms of h µν and B µ . Therefore, the field equations for h µν can be written in Einstein like form, though their explicit form is cumbersome and does not bring any useful new insight. This confirms that we can interpret the auxiliary metric as we did when (25) was introduced, as the metric that accounts for the cumulative effects of mass and energy.
Regarding the bumblebee field equation, using (29) we are able to get rid of the B ν R µν term in (16), thus arriving at This equation is richer than its equivalent in the metric case, for which the ξ corrections are absent. The new terms induce modifications on the effective potential that depend on the presence of other matter fields, thus implying new phenomenology. In particular, the effective mass of the field may change in regions with high energy densities, potentially leading to large effective masses and chameleon-like screening mechanisms inside and around massive objects.

Weak gravitational field, post-Minkowskian limit
Having found the solution of the connection equation and the explicit formal form of the space-time metric, we are now ready to explore the limit of negligible curvature, focusing our attention on how the sources of non-metricity modify the effective theory (post-Minkowskian limit). For this purpose, let us focus on non-gravitational experiments on Earth's surface, so that one can safely neglect all the corrections to the Minkowski metric coming from Newtonian and post-Newtonian corrections. Thus we require 3 h µν ≈ η µν . Let us also consider ξ to be a small coupling and study its perturbative effects. From the above relations between the auxiliary and the space-time metrics (22,23), we can write where B λ B λ = η µλ B µ B λ and ε µν = ξ (B µ B ν − 1 2 B λ B λ η µν ) +O(ξ 2 ). Here we see that even when Newtonian and post-Newtonian corrections can be neglected, the space-time metric is locally departing from its Minkowskian value due to local contributions sourced by the bumblebee field B µ . Since all fields couple to the metric, as a result, in this post-Minkowskian approximation of our theory all the matter fields will couple to the bumblebee due to the unconventional way in which the connection mediates between the geometry and the matter.

Effective dynamics of scalar and fermionic fields.
We will next proceed to study the effective dynamics of scalar and spinor matter fields in this post-Minkowskian scenario. Regarding spinors, it is important to note that they provide a nonzero contribution to the connection equation via torsional terms. Those terms have been omitted in our presentation of the field equations for simplicity. A more careful analysis, along the lines of [28], justifies our choice because in the case of bosonic fields, torsion can be trivialized by a simple choice of projective gauge. For fermions, however, the torsion picks up contributions that cannot be gauged away. However, such terms do not modify the equation satisfied by the symmetric part of the connection, which is still of the form (15) and admits the Christoffel symbols of h µν as solution. Therefore, the resulting effective metric will be the same as (31). The new torsional terms will appear as new fermionic contact interactions on the right-hand side of (14) and will be Planck-scale suppressed (see e.g. [44] or [40] for the RBG case), so that we can neglect them for our purposes.
From (31) it is clear that the inverse effective metric is given by where the indices are raised and lowered by the usual Minkowski metric. Similarly, one has √ −g ≃ 1 − ξ 2 B µ B µ . As a result, the scalar Lagrangian within this approximation is given by Using the definition of ε µν and integrating by parts the former equation, we arrive at where ✷ ≡ η µν ∂ µ ∂ ν . The interaction terms between the bumblebee and the scalar field naturally stem from special features of this metric-affine model, as explained previously. As the bumblebee field develops a non-trivial VEV, we see that the last two terms in (34) carry coefficients for Lorentz violation, which will be discussed in more detail later.
Regarding spinor fields, we start with the Hermitian Lagrangian where e a µ is the tetrad field satisfying e a µ e b ν g µν = η ab . Given that the connection Γ is just the Levi-Civita connection of h µν , in the weak field regime h µν ≈ η µν we can approximate ∇ Γ µ Ψ ≈ ∂ µ Ψ up to Planck scale suppressed and O(ξ 2 ) torsion corrections (see e.g. [36]). Taking this into account and using that e a µ = δ a µ − 1 2 ε µ a = δ a µ + ξ 2 we find that the spinor Lagrangian can be written as (37) Before moving forward, it is worth calling attention to the fact that the combination of all the complex terms present in Eq. (37) result in a Hermitian Lagrangian as a consequence of Eq. (35). Similarly to the scalar sector, the O(ξ ) terms in the above Lagrangian will contribute to Lorentz violation coefficients when the bumblebee acquires a VEV. We will next proceed to analyze their physical implications.

Lorentz-violating coefficients
Let us now outline the physical effects related to a spontaneous breaking of Lorentz symmetry by the bumblebee VEV B µ = b µ in a weak gravitational field. Generally, observables which couple to b µ will be sensitive to the spontaneous breaking of Lorentz symmetry by the bumblebee field. Since the present model displays non-minimal couplings between the bumblebee field and the matter sources through the non-Riemannian part of the connection, there will arise several Lorentz-violating (LV) coefficients in the effective matter sector once the affine connection has been integrated out.
A straightforward and simple method to examine the effects of the Lorentz violation is to consider the VEV being a fixed-norm vector or, similarly as in Minkowski space, to take all their components to be constant. Such a choice is reasonable because we are dealing with laboratory experiments. Nonetheless, for the sake of completeness, the non-constant b µ case is presented in Appendix A. Therefore, from now on we restrict our analysis to the standard b µ constant case [3]. We note that the allowed values for b µ have been restricted from various experiments, see e.g.
[? ] and references therein. At the same time, let us note that the LV coefficients generated by the non-metric part of the connection in metric-affine bumblebee gravity are analog to some of the parameters of the LV vector sector of the Standard Model Extension. Therefore, data from table D15 in [21] allows us to constrain the norm of b µ to be less than 10 −43 GeV.

Scalar field
Let us first study the Lagrangian for a scalar field propagating on top of the bumblebee background, which takes the form where s µν = ξ b µ b ν . The O(ξ ) terms will typically induce LV coefficients through the VEV of the bumblebee field. The s µν term constitutes a modification of the standard kinetic term which can be encoded in an effective metric for the scalar field of the form g µν eff = η µν − ξ s µν . Hence, a "wrong" signature of the LV coefficient s µν could trigger ghost-like instabilities around strong enough bumblebee backgrounds. Note however that in that case, the perturbative expansion would break down since ξ b 2 would be O(1), and a full non-perturbative analysis would be required. The correction to the mass term in (34) can also be encoded in an effective mass of the form m 2 eff = m 2 (1 − (ξ /2)b 2 ) which could also trigger tachyonic-like instabilities for a space-like bumblebee VEV (again non-perturbative effects could play a non-negligible role).
In order to explore potential instabilities in more detail, let us analyze the particular cases of a time-and space-like constant b µ by working out their respective dispersion relations. Starting with a time-like VEV b µ = [b, 0, 0, 0] we find the dispersion relation This dispersion relation is healthy for positive values of ξ but could potentially develop ghost-and tachyonic-like instabilities for negative values of ξ , although in this case the higher-order terms might be relevant to the discussion. Considering a space-like choice for b µ = [0, b] we end up with Again, this dispersion relation is healthy for positive values of ξ . For negative values of ξ a tachyonic-like instability as well as a ghost-like instability (in directions which are non-orthogonal to b) could potentially arise, and again the higher-order terms might be relevant to the discussion if they appear.

Dirac field
Let us now turn our attention to the spin 1/2 fields. To explore the physics of our interest in a more convenient way, we will work with the decomposition of LV coefficients that is more commonly used in the literature [45,46]. To that end, let us rewrite the weak-field spinor action 37 as follows where Γ µ and M are elements of the 16-dimensional Clifford algebra defined by the Dirac gamma matrices. We can thus expand them in the usual basis of this algebra as Γ µ = e µ I + (δ µ α + c µ α )γ α + d µ α γ 5 γ α + i f µ γ 5 + where c µ α , d µ α , e µ , f µ , g µ λ α , a µ , k µ and l µν are LV coefficients. Comparing 37 to 41 and 42, we find the non-zero LV coefficients We see that within the metric-affine bumblebee model, the LV coefficients that appear provide a modification of the fermionic mass through m eff and a modification of the standard kinetic term through c µ α . In general, these will introduce modifications in the dispersion relation of spin 1/2 fields. To that end, we first notice that the modified Dirac equation is and multiplying on the left by iΓ µ ∂ µ + M we arrive at By using now the relations we find the following dispersion relation As we did for the scalar field, let us particularize for a constant bumblebee background of both time-and space-like types. Starting with a time-like VEV b µ = [b, 0, 0, 0] we get the dispersion relation While for positive values of ξ this dispersion relation is perfectly well-behaved, for negative values of ξ there could be instabilities if b 2 ξ ≥ 1, but higher order terms would be non-negligible in this case and 51 would not be trustable anymore. For a (constant) space-like vector b µ = [0, b] we obtain For negative values of ξ this is well behaved, though if it is positive then ghost-like instabilities could arise beyond the perturbative level (see [47] for a discussion on the typical energy scales at which these instabilities become relevant).

Summary and conclusions
In this work we have formulated the bumblebee model within the metric-affine formalism. By solving the equations of motion of the connection, we have been able to express the metric as a function of this field and of an auxiliary metric which accounts for the standard effects of the gravitational interaction. To our knowledge, this is the first time that a solution for the connection in a curvature-based metric-affine gravity theory with spontaneously broken Lorentz symmetry has been found. The methods used to solve the connection are analog to those commonly employed in Ricci-based gravity theories (see e.g. [28,29,48]), but its qualitative properties are rather different from those theories. In particular, while in RBGs the metric picks up local corrections that depend on the stress-energy tensor of the matter fields, here those corrections are entirely determined by the bumblebee field itself. This occurs due to the absence of higher-curvature terms in the action and to the non-minimal coupling of the bumblebee to gravity via the Ricci tensor.
From our analysis it follows that, unlike in RBGs, where the non-metricity is given by gradients of the energy-momentum density, in our model it is the gradient of the bumblebee field that generates a non-vanishing non-metricity. The fact that the bumblebee has a VEV which breaks Lorentz symmetry allows for a background non-metricity that could fit with the proposal in [42]. This is the first gravitationally-induced non-metricity model with a VEV that we are aware of.
An immediate effect of the metric dependence on the bumblebee field is that all the matter fields couple to it. Since this coupling is not a gauge one it can generate, for instance, the coupling of a vector to a neutral scalar. In the weak-field limit, we have seen that the resulting theory looks like a bumblebee coupled to matter with non-linear interactions in Minkowski space. Therefore, this theory is naturally treated as an effective theory where the role of the energy scale is played by ξ −1/2 (see e.g. [49]). Typically, for the VEV solution we have checked that the theory can present ghost and tachyonic-like instabilities depending on the range of the non-minimal parameter ξ and the explicit form of the VEV. A natural continuation of this study consists in calculating quantum corrections in this effective theory. The search for astrophysically relevant solutions beyond the weak-field limit presented here is also another research avenue currently under consideration.
Appendix A: Implications on taking b µ non-constant Let us start this appendix by considering the scalar Lagrangian, which takes the form . Note the presence of the additional coefficient t µ in relation to the standard case (it vanishes for constant bumblebee VEVs). This coefficient introduces an imaginary term in the scalar dispersion relation. The modified scalar dispersion relations now looks Note the existence of imaginary terms in the above dispersion relation, which only vanish if one considers a frame where t µ p µ = 0. This undesired property leads to complex eigenvalues of the Hamiltonian operator, which turns out to be non-Hermitian. As a consequence, we conclude that a non-constant b µ produces an effective Minkowskian theory with serious problems in its dispersion relation, as displayed in the former equation. As pointed out in [1], this justifies our choice of taking the Lorentz coefficient b µ as a constant in order to avoid undesired effects such as a non-Hermitian Hamiltonian and violation of energy conservation due to the fact that a non-constant b µ would play the role of an external field whose presence would break space-time homogeneity. For spinor fields, the dispersion relation is modified in relation to the constant b µ as − Γ µ Γ ν ∂ µ ∂ ν − Γ µ (∂ µ Γ ν )∂ ν − iΓ µ ∂ where σ µν = i 2 [γ µ , γ ν ] we find the following dispersion relation ξ mγ µ ∂ µ b 2 + iγ µ γ α ∂ µ a α + iγ µ γ β p i ∂ µ c iβ + + i 2 a α σ αi p i − 2ma µ γ µ . (A.7) Notice that the derivative terms dependent on a µ are accompanied by gamma matrices and vanish if the bumblebee background varies slowly enough. In whole analogy with the scalar field situation, spinor fields suffer from the same instabilities.