Classical limit for Dirac fermions with modified action in the presence of the black hole

We consider the model of Dirac fermions coupled to gravity as proposed in \cite{VolovikBH}, in which superluminal velocities of particles are admitted. In this model an extra term is added to the conventional Hamiltonian that originates from Planck physics. Due to this term a closed Fermi surface is formed in equilibrium inside the black hole. In this paper we propose the covariant formulation of this model and analyse its classical limit. We consider the dynamics of gravitational collapse. It appears that the Einstein equations admit a solution identical to that of the ordinary general relativity. Next, we consider motion of particles in the presence of the black hole. Numerical solutions of the equations of motion are found which demonstrate that the particles are able to escape from the black hole.


I. INTRODUCTION
The Schwarzschild black hole (BH) solution [2] may be brought to the form, which is especially useful for the consideration of the motion of matter. This is the so -called Painlevé -Gullstrand (PG) black hole [3,4]. In the corresponding reference frame space -time looks like flat space falling down to the center of the BH with velocity that depends on the distance to the center. Such a representation also exists for the Reissner-Nordstrom BH and even for the Kerr BH [5,6]. The structure of the BH solution in the PG reference frame prompts to consider the analogy to the motion of fluid in laboratory. Such an analogy has been considered in the framework of the theory of 3 He superfluid in [7]. It allows to calculate in a demonstrative way the Hawking radiation [8] (see also, for example, [9][10][11][12] and references therein).
The consideration of Dirac fermions in the PG reference frame leads to the unexpected observation, that inside the BH the fermion states with vanishing energy form the surface in momentum space. In equilibrium, at vanishing temperature, it becomes the Fermi surface and separates the region of occupied quantum states from the region of the vacant ones. For the ordinary Dirac fermions such a surface is open and infinite. The analogy with condensed matter physics suggests that the particle Hamiltonian is to be modified in such a way so that the resulting Fermi surface will become finite and closed. It has been proposed in [1], that such a modification occurs due to the Planck physics. The corresponding term has been added to the Dirac Hamiltonian. It leads to several consequences for the BH physics. First of all, the analogy to the BH in the PG reference frame has been found within the class of the recently discovered materials called Weyl semimetals [13][14][15][16][17][18][19][20]. Bloch electrons within those materials behave similarly to the elementary particles. In the special type of such materials called the type II Weyl semimetals (WSII) [22] the dependence of energy of electrons on momenta [21] possesses an analogy to that of the particles in the interior of the BH [23,24].
In [25,26] it was noticed that if a closed finite Fermi surface inside the BH is formed, there should exist particles that escape from the BH without tunneling. In the present paper we take a step back and consider the model proposed in [1] on the classical level. We suppose that a careful consideration of the classical dynamics should precede the more sophisticated discussion of the quantum BH, although the extra term added to the particle Hamiltonian becomes relevant at Planck energies.
This extra term contains the time -like four -vector n µ . In the Painlevé -Gullstrand reference frame it marks the direction of time. In the covariant theory there should be no such preferred direction of time. We assume that it appears as a result of a dynamical symmetry breaking. This symmetry breaking, in turn, results in the appearance of the massless Goldstone modes. Here we shall not discuss the physics of those massless excitations.
First of all we propose the covariant formulation of the discussed model. It contains the vector field n µ , which, after the spontaneous breakdown, acquires its particular value that points in the direction of time in the Painlevé -Gullstrand reference frame. The classical equations of motion for the corresponding point -like particles admit motion with superluminal velocity. Therefore, unsurprisingly, the particles may escape from the black hole already on the classical level. For the discussion of the theories that admit superluminal velocity of particles we refer to [43].
The paper is organized as follows. First of all, in Sect. II we recall the general properties of the Dirac fermions in the PG reference frame in the presence of the charged BH. In Section III we propose the covariant formulation of the model of [1] and derive its classical limit. In Sect. IV we derive the expression for the stress -energy tensor of the noninteracting particles in the PG reference frame. In Sect. V the gravitational collapse in this model is considered.
In Sect. VI we present the results of the numerical solution of the classical equations of motion for the motion of particles in the presence of the existing black hole (at the stage when the gravitation collapse is finished). The physical significance of the results is dicussed in the concluding section VII. In the present paper we will mainly be interested in neutral black holes. However, we start our consideration from the ansatz that corresponds to a charged black hole. The charge is assumed to be small, so that it modifies the BH solution in the small vicinity of the BH center. We will add also another modification in the small vicinity of the center, which will cause the metric to be regular everywhere. Thus we start from the charged black hole metric in the PG coordinates may be considered as a velocity of space falling towards the center of the BH. Q is the charge of the BH, m P is the Plank mass, while M the BH mass. The vielbein is given by and the inverse vielbein E µ a is The metric is equal to The action for the Weyl fermions has the form where symbols R (L) mark the right -handed / left -handed fermions. Their Hamiltonians are The Dirac mass term mixes the right -handed and the left -handed fermions It appears, that the spin connection does not enter this action. Next, following [1], we introduce the term that breaks the Lorentz invariance: Here parameter µ is assumed to be of the order of the Plank mass. For r < r 0 = Q 2 2M the velocity v becomes imaginary. We guess that due to interaction with matter the dependence of v on r is modified within the BH, and v remains real and tends to zero at r = 0 (see also [5]). In the present paper we model this form of v via the following modification Fig. 1 in the system of units with m P = 1 at M = mP 2 and Q = 0.4. The two horizons are placed at and For r > r + there are the ordinary Dirac fermions. At m = 0 the Fermi point appears. Between the two horizons r − < r < r + at m = 0 there is the type II Dirac point, while |v| is larger than light velocity. Neglecting the derivatives of v we come to the following expression for the particle energy Thus assuming that v varies slowly, we come to the conclusion that between the two horizons the particle energy vanishes along the closed surface in momentum space. Its form is represented in Fig. 2 for the particular choice of parameters. It is worth mentioning, that the Hamiltonian of the form of Eq. (11) admits motion with the superluminal velocity on the classical level. Therefore, unsurprisingly, in the considered model the particles are able to escape from the Black hole.

III. COVARIANT FORMULATION OF THE THEORY AND ITS CLASSICAL LIMIT
Here we propose the covariant modification of the model of [1] considered in the previous section. Namely, we consider the Dirac field Ψ with action Here we denote D µ = ∂ µ + 1 8 ω abµ [γ a , γ b ] and h µν = g µν − n µ n ν . The 4-vector n in the PG reference frame marks the time direction: Correspondingly n µ = (1, −v, 0, 0) . The spin connection is given by ; g µν = e a µ e b ν η ab . For the details see Ref. [38], and Refs. [39][40][41]. By n we denote the vector field. Its action may be taken in the form: At sufficiently large values of λ this vector field becomes non -propagating. Its vacuum average gives rise to spontaneous symmetry breaking. The 4 -vector n in the Painleve -Gullstrand reference frame is constant and marks the time direction: Correspondingly, n µ = (1, −v, 0, 0). The appearance of this vector breaks the group of general coordinate transformations in four -dimensional space -time to the group of general coordinate transformations in three -dimensional space. As a result three massless Goldstone modes appear corresponding to the broken boosts. We do not discuss here the physics of these massless excitations. However, we do not exclude that such excitations may play a certain role in the formation of dark matter.
The only component of the spin connection (in spherical coordinates) is ω r0r = −ω 0rr = −v ′ (r). This term disappears from the first term in the Painlevé -Gullstrand reference frame. But it is essential for the calculation of the stress -energy tensor.
The first step towards the classical theory is the consideration of the above model with spin neglected. This corresponds to the consideration of the scalar field instead of the Dirac spinor and the corresponding action. Hence the theory is that of a scalar field Φ with the action: There is no precise correspondence between Eqs. (12) and (13). In the transition the spin degrees of freedom and the corresponding terms in the lagrangian were neglected. Variation of this action with respect to Φ gives the wave equation In the semiclassical limit we substitute −i∂ k by momentum p and i∂ 0 by energy E. This gives the following relation for p µ = (E, −p): In the Painlevé -Gullstrand reference frame we get: This equation gives rise to the classical particle Hamiltonian of [1] given by Eq. (11).

IV. THE STRESS -ENERGY TENSOR OF THE NON -INTERACTING CLASSICAL PARTICLES
A. General expression for the stress -energy tensor In this section we consider matter that consists of the non -interacting particles. Our aim is to consider the gravitational collapse. Therefore, following [42] we consider the generalization of the Painlevé -Gullstrand spacetime: with The function m(t, r) is to be obtained through the solution of Einstein equations. Although we are going to calculate the stress -energy tensor of the classical system, we prefer to start from the model of the scalar field with action Eq. (13). We will calculate the stress energy tensor of the corresponding quantum system and take the classical limit at the end of the calculation. We are to calculate the stress energy tensor through the variation of action with respect to the variation of metric: In a similar way we calculate the current density as In the semiclassical limit the oscillating factors in the radial wave functions are given by Φ(r, t) ∼ e iS ∼ e −iEt+i r r 0 p(r)dr while the electric current may be identified with the product of the particle density ρ and the velocity of substance. This leads to the following semiclassical relation in the generalized Painlevé -Gullstrand reference frame: Here ρ(x, t) is the particle density. The number of particles in a small volume Ω around the given space -time point is equal to Ω d 3 xρ(x, t). It is supposed that all these particles have equal values of energy E(x, t), velocity V(x, t) =rV (x, t), and momentum p(x, t) =rp(x, t). Those quantities are related to each other via the classical equations of motion The above equations enable us to relate the absolute value of Φ with the physical quantities ρ and V. In an arbitrary reference frame we have the similar relation Here is the density of particles. Therefore, in the Painlevé -Gullstrand coordinates we identify while in the arbitrary reference frame In the same limit the stress -energy tensor is given by Equations of motion give B. The stress energy tensor in the limit µ → ∞ Let us demonstrate, how Eq. (22) gives rise to the conventional stress energy tensor of the non -interacting particles in the limit µ → ∞. We have In the generalized PG reference frame we have This gives Here u j is the four -velocity of the particles/substance. We come to the conventional expression where ǫ = ρm ds dt is the energy density in the rest frame of medium (notice that ρ is the particle density in the given reference frame while ρ ds dt is the particle density in the rest frame).

C. Expression for the stress -energy tensor for finite µ, in the case when the substance is co -moving with the space flow
In the general case the following expression for the stress energy tensor should be obtained: T µν = u µ u ν f u (u, n, µ, ǫ) + n µ n ν f n (u, n, µ, ǫ) + (u µ n ν + u ν n µ )f un (u, n, µ, ǫ) + g µν f g (u, n, µ, ǫ) Scalar functions f u , f n , f un , f g depend on the four -velocity u, the four -vector n, constant µ, and the energy density ǫ = ρmds/dt. In the important particular case, when velocity of substance V everywhere is equal to v the classical equations of motion give p = 0 and we obtain the especially simple result:

V. DESCRIPTION OF THE GRAVITATIONAL COLLAPSE IN THE GENERALIZED PAINLEVÉ -GULLSTRAND COORDINATES
In this section we consider the gravitational collapse of matter that consists of the non -interacting particles. Those particles being placed into the Painlevé -Gullstrand spacetime have the Hamiltonian The generalization of the PG spacetime [42] has the following metric with The function m(t, r) is to be obtained through the solution of Einstein equation. It was shown above that the noninteracting matter with the Hamiltonian of [1] has the stress energy tensor equal to that of the conventional matter in the generalized Painlevé -Gullstrand coordinates in the important particular case, when at each point the velocity of matter is precisely equal to v. The problem for the gravitational collapse of matter with this stress -energy tensor is solved in [42]. We repeat here the solution for completeness.
In the spherical coordinates the Einstein equation 8πT µ ν = R µ ν − (1/2)δ µ ν R receives the form (we use in this section the system of units with m P = 1): Here dots represent the differentiation with respect to time while ′ is the differentiation with respect to the radial coordinate r. From these equations one has the identity For the noninteracting particles (perfect fluid) in the above particular case the energy -momentum tensor T µ ν has the form: Here ǫ is the energy density, and u µ = (1, −v(t, r), 0, 0) the radial four-velocity of the fluid. Correspondingly u µ = (1, 0, 0, 0). For the definiteness let us assume that v is directed along the x axis. Then It follows from Eqs. (34) and (33) that v(t, r) is given by v(t, r) = 2m(t, r) r .
Next, the integration of the above mentioned equations of motion gives and ǫ = 1 6πt 2 , t < 0.
We come to the following pattern of the gravitational collapse. If the velocity of matter at the starting moment t 0 is equal to the function v(r, t 0 ) of the generalized Painlevé -Gullstrand reference frame, and everywhere the threemomentum of the particles vanishes at t = t 0 , then the Einstein equations have the solution given by Eqs. (36), (37), (38). Thus, matter contracts towards the center of the BH together with the "falling" space. This gives m(t, r) = 2r 3 9t 2 m 2 P , v(t, r) = 2r 3|t| (We restore here the expression that contains the Planck mass m P explicitly due to unit considerations.) As a result the position of the horizon (where the velocity equals 1) depends on time: One can see, that at each finite value of t the velocity of space fall vanishes for r = 0 only. The space -time metric remains regular everywhere. The collapse of dust placed within a sphere leads at t → −0 to the formation of the ordinary Painlevé -Gullstrand BH (see [42]). In the next section we will consider the classical motion of particles in the formed BH regularized in the small vicinity of its center as explained in Sect. II. The classical Hamiltonian of the quasiparticles in the presence of the black hole has the form: For the radial motion the corresponding classical equations are: The Fermi surface for the particular value of r crosses the axis of radial momentum at One can see, that the Fermi surface is not formed immediately behind the external horizon. Instead, it is formed at The typical classical trajectories of the particles are calculated and represented in Figs. 3, 4, 5. The particles that fall together with the "vacuum" reach the center of the black hole and stay there. However, if the initial momentum is nonzero, the particles that have fallen down to the black hole receive large values of radial momentum in a small vicinity of the BH center. As a result they escape from the BH. If the initial momentum was directed to the center of the BH, then the particle traverses the BH and escapes it from a diametrically opposite point. If the initial momentum is in the opposite direction, then its velocity reverses the sign close to the center of the BH, the particle is turned back. In the exterior of the black hole the momentum is decreased, and the particle velocity changes the sign again. The particle falls down again thus forming oscillations.
Our numerical data allow to estimate the typical time period of those oscillation for m ≪ µ ∼ m P ≪ M (when the amplitude is of the order of the horizon) T ∼ 20 M/µ 2 (see Fig. 6). This value for the solar mass BH (in seconds) is To conclude, in the present paper we consider the model of noninteracting Dirac fermions with the modified Hamiltonian proposed in [1]. In this model the extra ∼ p 2 term is added to the Dirac Hamiltonian. First of all, we propose the covariant generalization of this model. The resulting field theory is defined in terms of the Dirac spinor field. It depends on the background field n ν that equals n µ = (1, 0, 0, 0) in the Painlevé -Gullstrand reference frame. The field n points into the direction of time in this coordinate system. It appears as a result of the spontaneous symmetry breaking. This symmetry breaking also leads to the appearance of the massless Goldstone modes. The direct interaction term of those modes with matter is suppressed by the factor 1/µ, where µ is of the order of Plank mass. But they are coupled to gravity. The consideration of the physics of those modes is out of the scope of the present paper. But we do not exclude, that, in certain theoretical schemes, they contribute the dark matter. Next, we consider the classical limit of the obtained system. First, neglecting spin we come to the theory of the scalar field with a certain action (depending on n). Next, we consider the semiclassical approximation within this theory, which gives both the classical particle Hamiltonian of [1], and the implicitly defined expression for the stressenergy tensor of medium that consists of the noninteracting particles. It appears, that the Einstein equations admit the gravitational collapse solution identical to that of the system of the conventional noninteracting particles. This solution describes the dust falling together with space -time.
Finally we describe the dynamics of particles in the presence of the existing black holes. It appears, that only those particles remain inside the BH, which fall with velocity v entering the expression for the Painlevé Gullstrand metric. The particles that fall towards the center of the BH with nonzero momentum either pass through the BH and reach infinity or turn back at the small vicinity of the BH center, escape from the BH, and proceed the oscillatory motion. Sure, this means that the particles are able to move with the velocity larger than the speed of light. Although the considered theory is manifestly covariant (as explained in Sect. III), the geodesic lines are already not the solutions of the classical equations of motion of point-like particles. The solutions of those equations may correspond to the space -like pieces of the particle worldlines, which does not contradict the general covariance. According to our estimates for the BH with the solar mass the typical period of the mentioned oscillations (when interactions are neglected) is smaller than one second. This means, that if the effective Hamiltonian of particles indeed receives the considered contribution from Planck physics, then we cannot ignore it in the dynamics: matter that has fallen to the BH escapes from it within the observable period of time.
One of the authors (M.A.Z.) kindly acknowledges useful discussions with G.E. Volovik.