Scalar-Gradient Bumblebee Field in Lorentz Gravity

• The scalar-gradient bumblebee field is the scalar component emerging from a bumblebee vector with a nonzero vacuum expectation value that breaks Lorentz symmetry.
• It appears in two forms—a Kaluza–Klein-generated dilaton and a massive probe scalar—controlling black-hole dynamics, quasinormal spectra, and cloud formations.
• Lorentz-violating parameters rescale effective radial potentials, altering damping rates, stability conditions, and thermodynamic properties in multiple gravity setups.

Searching arXiv for the cited bumblebee-gravity papers to ground the article in current records. The scalar-gradient bumblebee field, as represented in current Einstein-bumblebee literature, is the scalar sector associated with a bumblebee vector field whose nonzero vacuum expectation value spontaneously breaks Lorentz symmetry. In practice, this sector appears in two distinct but related forms: as a Kaluza–Klein-generated dilaton coupled directly to the bumblebee kinetic and potential terms, and as a massive probe scalar whose propagation, quasinormal spectrum, and stationary clouds are controlled by bumblebee-deformed black-hole geometries. Across these realizations, the common mechanism is that Lorentz-violating parameters rescale the effective radial dynamics, thereby modifying damping rates, cloud existence lines, thermodynamic structure, and horizon-scale phenomenology without always altering the underlying kinematic thresholds (Lessa et al., 2023, Liu et al., 2022, Quan et al., 27 Jan 2025).

1. Conceptual setting in Einstein-bumblebee gravity

Einstein-bumblebee models introduce a vector field $B_\mu$ or $B_M$ with a self-interacting potential that fixes a nonzero vacuum expectation value (VEV). In the four-dimensional rotating black-hole analysis, the theory is defined by

$$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$

and Lorentz violation is encoded in the dimensionless parameter

$\ell=\varrho b_0^2,$

with the rotating background taking the bumblebee field to be purely radial,

$b_\mu=(0,b_r,0,0).$

The dimensionless spin parameter is

$\tilde a=\frac{a}{M},$

assumed small in the slow-rotation treatment (Liu et al., 2022).

In the three-dimensional AdS construction, the rotating BTZ-like background arises from the Einstein-bumblebee action with negative cosmological constant $\Lambda=-1/\ell^2$, and spontaneous Lorentz-symmetry breaking is parametrized instead by $s$. The metric is

$$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$

Here $s>-1$ is required for regularity. A notable structural result is that the horizon radii

$B_M$0

and the horizon angular velocity

$B_M$1

are independent of $B_M$2 (Quan et al., 27 Jan 2025).

This framework suggests that the scalar-gradient sector is not a single field-theoretic object but a class of scalar dynamics governed by the same Lorentz-violating bumblebee background. In that sense, the scalar sector is unified less by field identity than by the way the bumblebee VEV reorganizes effective potentials, radial equations, and black-hole response.

2. Kaluza–Klein origin of the scalar-gradient sector

A higher-dimensional realization makes the scalar component explicit. In the Einstein-Bumblebee-scalar theory obtained from Kaluza–Klein reduction, the starting point is a Kostelecký–Samuel / Einstein-bumblebee model in $B_M$3 dimensions with a VEV condition

$B_M$4

The reduction ansatz

$B_M$5

is constrained by the Einstein-frame condition

$B_M$6

The scalar dilaton $B_M$7 is then defined through

$B_M$8

with KK-induced coupling

$B_M$9

For $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$0,

$$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$1

and the $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$2 case yields the special coupling used later in the thermodynamic analysis (Lessa et al., 2023).

After reduction, the effective $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$3-dimensional theory contains the metric $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$4, the bumblebee vector $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$5, the scalar/dilaton $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$6, and the fluctuation sector

$$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$7

The decomposition is defined by the projectors

$$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$8

so that $$\mathcal{S}=\int d^4x \sqrt{-g}\left[\frac{1}{16\pi G_N}\left( \mathcal{R}+\varrho B^aB^b \mathcal{R}_{ab} \right)-\frac{1}{4}B^{ab}B_{ab}-V \right]+\mathcal{L}_M,$$9 is transverse and $\ell=\varrho b_0^2,$0 is longitudinal. The potential becomes quadratic in the longitudinal mode,

$\ell=\varrho b_0^2,$1

which identifies $\ell=\varrho b_0^2,$2 as the massive mode (Lessa et al., 2023).

For a static, spherically symmetric ansatz with radial VEV, the longitudinal mode does not propagate and becomes a constant $\ell=\varrho b_0^2,$3. The resulting action is

$\ell=\varrho b_0^2,$4

with

$\ell=\varrho b_0^2,$5

The transverse mode behaves as an effective Maxwell field, whereas the longitudinal mode acts as a cosmological-constant-like source. The dilaton multiplies both sectors through $\ell=\varrho b_0^2,$6, so the scalar gradient modulates both gauge-like and Lorentz-violating contributions (Lessa et al., 2023).

3. Static black-hole solutions and thermodynamic status

The four-dimensional reduced theory admits static, spherically symmetric solutions with

$\ell=\varrho b_0^2,$7

The field equations are supplemented by the electric-type field strength

$\ell=\varrho b_0^2,$8

In the decoupling limit $\ell=\varrho b_0^2,$9, the metric becomes

$b_\mu=(0,b_r,0,0).$0

namely a charged de Sitter–Reissner–Nordström solution (Lessa et al., 2023).

For nonzero $b_\mu=(0,b_r,0,0).$1, the dilatonic branch is specified by

$b_\mu=(0,b_r,0,0).$2

and

$b_\mu=(0,b_r,0,0).$3

The solution has several distinguished limits: for $b_\mu=(0,b_r,0,0).$4 it reduces to Schwarzschild; for $b_\mu=(0,b_r,0,0).$5 it has one horizon and strong dilaton effects; and for $b_\mu=(0,b_r,0,0).$6 it has no horizon and is a naked singularity. The positivity of $b_\mu=(0,b_r,0,0).$7 and the requirement of spontaneous Lorentz breaking impose

$b_\mu=(0,b_r,0,0).$8

and

$b_\mu=(0,b_r,0,0).$9

These inequalities control both horizon existence and the physicality of the Lorentz-violating potential (Lessa et al., 2023).

The thermodynamic analysis is performed mainly for the KK value $\tilde a=\frac{a}{M},$0. The local stability criteria are

$\tilde a=\frac{a}{M},$1

With $\tilde a=\frac{a}{M},$2, the temperature is

$\tilde a=\frac{a}{M},$3

the entropy is

$\tilde a=\frac{a}{M},$4

and the first law takes the form

$\tilde a=\frac{a}{M},$5

The heat capacity is positive when

$\tilde a=\frac{a}{M},$6

so the solution is thermodynamically stable in that regime, whereas the isothermal charge susceptibility is negative,

$\tilde a=\frac{a}{M},$7

indicating electrical instability (Lessa et al., 2023).

4. Scalar perturbations of slowly rotating Einstein-bumblebee black holes

A distinct realization of the scalar sector is the massive Klein–Gordon field on a slowly rotating black hole in Einstein-bumblebee gravity. The background is a Kerr-like bumblebee metric expanded to $\tilde a=\frac{a}{M},$8, with outer and inner horizons

$\tilde a=\frac{a}{M},$9

The construction is explicitly approximate: the field equations are violated only at order $\Lambda=-1/\ell^2$0, so for sufficiently small $\Lambda=-1/\ell^2$1 and $\Lambda=-1/\ell^2$2 it is acceptable as a slowly rotating solution (Liu et al., 2022).

The scalar is decomposed into spherical harmonics, and because the background is axisymmetric the $\Lambda=-1/\ell^2$3-modes decouple. At first order in $\Lambda=-1/\ell^2$4, the radial problem becomes Schrödinger-like,

$\Lambda=-1/\ell^2$5

with effective potential

$\Lambda=-1/\ell^2$6

At second order, the scalar equation develops explicit $\Lambda=-1/\ell^2$7 mixing. The coupling is organized using angular identities involving

$\Lambda=-1/\ell^2$8

and a single master field is obtained through

$\Lambda=-1/\ell^2$9

The final scalar master equation is

$s$0

This second-order construction is the preferred one for the scalar problem because it captures rotational corrections more faithfully (Liu et al., 2022).

Quasinormal modes are computed with a matrix method and Leaver’s continued fraction method, under ingoing boundary conditions at the horizon and outgoing conditions at infinity. For the second-order treatment,

$s$1

with

$s$2

The matrix method discretizes the radial equation on $s$3, while the continued fraction method uses a Frobenius expansion in $s$4. For the scalar case, the standard three-term Leaver recursion applies. The comparison at $s$5 shows that the second-order approximation agrees with exact Kerr frequencies to about the $s$6 level up to $s$7 for $s$8, and the two numerical methods agree very well, with differences below $s$9 (Liu et al., 2022).

The main physical result is spectral rather than geometric: increasing $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$0 decreases $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$1, so the scalar perturbation decays more slowly, whereas $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$2 changes only mildly. The authors interpret this as a direct imprint of the bumblebee field on the effective potential, since $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$3 and $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$4 rescale both the centrifugal barrier and the frame-dragging correction (Liu et al., 2022).

5. Stationary scalar clouds in rotating BTZ-like backgrounds

In the rotating BTZ-like black hole of Einstein-bumblebee gravity, the massive scalar field $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$5 of mass $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$6 is separated as

$$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$7

leading to a radial equation in which the bumblebee parameter enters as an overall factor $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$8. For the non-extremal case $$ds^2 = -f(r)\, dt^2 +\frac{1+s}{f(r)}\, dr^2 + r^2 \left(d\varphi - \frac{j}{2r^2}dt\right)^2, \qquad f(r)=\frac{r^2}{\ell^2} - M +\frac{j^2}{4r^2}.$$9, the coordinate

$s>-1$0

maps the horizon to $s>-1$1 and the AdS boundary to $s>-1$2. The radial equation then reduces to hypergeometric form with parameters

$s>-1$3

together with corresponding $s>-1$4, $s>-1$5, and $s>-1$6 (Quan et al., 27 Jan 2025).

The horizon condition selects the ingoing branch by setting $s>-1$7 in the general solution. At the AdS boundary, the vanishing energy-flux condition yields a Robin boundary condition,

$s>-1$8

with Dirichlet at $s>-1$9 and Neumann at $B_M$00. The matching relation is

$B_M$01

Stationary clouds occur at the superradiant threshold

$B_M$02

for which $B_M$03 and $B_M$04. The resulting cloud quantization condition determines the existence lines in black-hole parameter space (Quan et al., 27 Jan 2025).

The dependence on $B_M$05 and $B_M$06 is highly structured. Increasing $B_M$07 makes clouds exist for smaller background mass at fixed $B_M$08, whereas increasing $B_M$09 makes clouds exist for larger background mass at fixed $B_M$10. Because both parameters enter through the combination $B_M$11, different pairs $B_M$12 can generate the same existence line. The explicit examples

$B_M$13

produce the same cloud-supporting curve. This degeneracy is only macroscopic: the radial profiles remain different (Quan et al., 27 Jan 2025).

Only the fundamental mode $B_M$14 supports stationary clouds. The associated QNM analysis shows that the imaginary part crosses zero only for the fundamental left-moving mode, identifying the cloud as a marginally bound state at the onset of superradiant instability. For Dirichlet and Neumann boundary conditions, the quasinormal frequencies are analytic, and throughout this construction the superradiance condition remains the standard one,

$B_M$15

because $B_M$16 is unchanged by $B_M$17 (Quan et al., 27 Jan 2025).

6. Interpretation, limits, and recurrent misconceptions

Several conclusions recur across these works. First, Lorentz-symmetry breaking does not uniformly manifest through horizon kinematics. In the BTZ-like problem, the horizon radii and $B_M$18 are unchanged by $B_M$19, even though the scalar spectrum, Robin threshold, and cloud existence lines are shifted. A common misconception is therefore that the bumblebee parameter must modify the superradiance inequality itself; in this system it does not (Quan et al., 27 Jan 2025).

Second, the scalar sector is not always a propagating new degree of freedom. In the KK-reduced static spherical background, the longitudinal bumblebee mode becomes a constant $B_M$20 and acts as an effective cosmological-constant-like source, while the propagating transverse mode behaves like a Maxwell field. The massive mode is therefore frozen into the background rather than behaving as an additional dynamical scalar in that setup (Lessa et al., 2023).

Third, the slowly rotating four-dimensional black-hole metric used for scalar QNMs is not an exact rotating solution. Its validity rests on the statement that the field equations are violated only at order $B_M$21, so the formalism is controlled only for sufficiently small Lorentz violation and spin. Within that regime, however, the second-order slow-rotation treatment is a quantitatively useful benchmark because it reproduces the Kerr limit more accurately as $B_M$22 (Liu et al., 2022).

Taken together, these results suggest a coherent picture: the scalar-gradient sector in bumblebee gravity is primarily a mechanism for reshaping effective radial dynamics. In one branch it generates dilaton-dressed black holes with a Maxwell-like transverse mode and a cosmological-constant-like longitudinal mode; in another it controls scalar damping, mode mixing, and cloud formation around Lorentz-violating black holes. The most robust signature across the cited systems is not a universal deformation of horizon kinematics, but a systematic reweighting of effective potentials and spectral thresholds by the Lorentz-violating parameters $B_M$23, $B_M$24, and $B_M$25 (Lessa et al., 2023, Liu et al., 2022, Quan et al., 27 Jan 2025).