On the validity of Strong Cosmic Censorship Conjecture in presence of Dirac fields

A well posed theory of nature is expected to determine the future of an observer uniquely from a given set of appropriate initial data. In the context of general relativity, this is ensured by Penrose's strong cosmic censorship conjecture. But in recent years, several examples are found which suggest breakdown of the deterministic nature of the theory in Reissner-Nordstrom-de Sitter black holes under the influence of different fundamental fields. Nevertheless, the situation has been reassuring for the case of astrophysically meaningful Kerr-de Sitter black hole solutions which seems to respect the conjecture. However, the previous analyses were done considering only the effect of scalar fields. In this paper, we extend the study by considering Dirac fields in Kerr-de Sitter background and show that there exist a parameter space which does not respect the conjecture.


Introduction
In recent years, advancement in technology brought a revolution in observational astrophysics that make it possible to test some of most intriguing predictions of general relativity [1,2]. While, due to the uncertainty in the data, the possibility of a modified theory of gravity is not discarded right away, the stature of Einstein's general relativity as the most successful theory of gravity remains unaltered [3][4][5]. This situation gives confidence to the scientists to propose interesting research work to test the theory in more and more extreme conditions. A possible alternative to this path is to find paradoxes within the theory of general relativity i.e. to check mathematical consistency of the theory. The existence of the Cauchy horizon in Kerr and Reissner-Nordström solutions is one of such paradoxes since the theory loses its predictive power beyond that region. However, soon it was realized that the Cauchy horizons are subject to blue shift instability that can turn them into curvature singularities under the influence of even small perturbations [6][7][8][9]. This phenomena led Penrose to propose the strong cosmic censorship conjecture which can be stated as follows: "for a generic initial data, the maximal Cauchy development (the largest manifold that is uniquely determined by Einstein's field equations from a given set of initial data ) is inextendible as a C 0 metric" [9][10][11][12][13]. This conjecture ensures that the observers who dare to cross the Cauchy horizon are torn apart by the infinite tidal forces. An another way to look at the problem is to consider the effect

Weak solution of Einstein Equation in presence of massless Dirac fields
The fate of strong cosmic censorship conjecture relies on the possibility of finding a solution of Einstein equation at the Cauchy horizon. Even if the metric is not differentiable (but continuous!) at the Cauchy horizon, one can still make sense of Einstein equation there by constituting a weak solution of the equation [14]. This can be understood by considering the effect of linear perturbation on the spacetime. Consider a massless Dirac field which satisfies the equation DΨ = 0, triggers a perturbation in the spacetime. Here D is the Dirac operator that acts on the spinor Ψ. Let the Dirac field act as a first order perturbation which induces a second order perturbation of the metric, denoted by h (2) µν which satisfies the following equation where, O is a second order differential operator and T Ψ µν is the stress-energy tensor for the Dirac field which can be expressed in the following form Even when h (2) µν is not differentiable at Cauchy horizon, we can have a solution of Einstein equation by multiplying Eq. (1) with a smooth, symmetric tensor K µν . By performing integration by parts, we obtain the following equation where, M is the maximal Cauchy development i.e. the region of spacetime uniquely determined by a set of generic initial data [10,15] and L † is the adjoint of the operator L. If this equation is satisfied for any smooth, symmetric function K µν , there exists a weak solution of Einstein equation provided that both sides of the equation remains finite. The requirement is fulfilled when Ψ belongs to H 1 loc . When such solutions exists at the Cauchy horizon, we can extend the metric across it which leads to the breakdown of strong cosmic censorship conjecture.

Dirac equation in Kerr-de Sitter spacetime and the criteria of violation of Strong Cosmic Censorship Conjecture
In the previous section, we have found that the presence of Dirac fields in Kerr-de Sitter spacetimes can lead to the violation of strong cosmic censorship, if the spinor field Ψ belongs to H 1 loc at Cauchy horizon. In this section, we closely inspect this condition. We solve the Dirac equation near the Cauchy horizon of Kerr-de Sitter black holes and then rewrite the condition in terms of black hole parameters. We start with a Kerr-de Sitter black hole spacetime in Boyer-Lindquist coordinate (t, r, θ, φ) whose line element can be expressed as follows [28] where, Here, M is the mass of the black hole, a is the black hole rotation parameter and Λ > 0 is the cosmological constant. Throughout the paper, superscript ' * ' denotes the complex conjugate of a quantity. Since we want to check the validity of strong cosmic censorship conjecture in presence of positive cosmological constant, we choose the values of the black hole parameters M , a and Λ in such a way that the spacetime possesses three distinct horizons. The position of Cauchy, event and cosmological horizon which we denote by r − , r + and r c respectively, can be found by solving the equation ∆ r (r) = 0. The properties of Dirac fields in this spacetime can be best understood in the framework of Newman-Penrose formalism [29,30].
Here, we choose our null tetrad to be [28,31,32] in the (t, r, θ, φ) coordinate. We can easily verify that the only non-vanishing inner product combination is given by the normalization condition, l · n = −1 and m ·m = 1. The advantage of this choice is that the tetrad vectors are regular across the Cauchy horizon. To see that we write the tetrad in outgoing Eddington-Finkelstein coordinate (u, r, θ, ϕ) by using the following transformation [33] Under this transformation the tetrad vectors take the following form Note that, the tetrad vectors are regular at the Cauchy horizon. In Newman-Penrose formalism, the equation for a massless Dirac field Ψ can be written as four coupled differential equations as follows [29] (D + − ρ) where, D = l · ∇, ∆ = n · ∇, δ = m · ∇,δ =m · ∇ are the directional covariant derivative along the tetrad vectors and α, β, γ, , π, ρ and τ are the spin coefficients (for details see [29,30]). Here, F 1 , F 2 , G 1 , G 2 denote the spinor components such that Ψ = (F 1 , F 2 , −G 2 , G 1 ) T . For our choice of tetrad (given by Eq. (4)), the non-vanishing spin coefficients are given by [31,32] where, ∆ θ = ∆ θ sin 2 θ. Due to presence of timelike and angular Killing vectors, the Dirac field can be decomposed as Moreover, if we take the following transformation [29,32] Eq. (6) can be decomposed into radial and angular parts which can be written as follows [31,32] where, Here, λ is a constant of separation, V (r) = (1 + α)[ω(r 2 + a 2 ) − am] and H(θ) = (aω sin θ − m csc θ).
Since we are interested in the behavior of the Dirac field near the Cauchy horizon, it is convenient to adopt the outgoing Eddington-Finkelstein coordinate (u, r, θ, ϕ) [33]. Moreover, if we set R + (r) = (∆ r ) 1 2 R + (r) and R − (r) = R − (r), Eq. (7) can be written in a more symmetric form D+ − R+ − = λ R− + . Near the Cauchy horizon, the radial equations can be written as follows and Ω − corresponds to the angular velocity of the black hole at the Cauchy horizon. Here, dr * = (1 + α)(r 2 + a 2 )dr/∆ r and dr φ = (1 + α)a dr/∆ r . The above equations have two independent solutions such that the spinor field can be written as follows Here, R+ − (r) represent some smooth functions of r which are non-vanishing at the Cauchy horizon and p = i(ω − mΩ − )/κ − where κ − corresponds to the surface gravity of the black hole at the Cauchy horizon. Given the solution of Dirac equation near the Cauchy horizon, we need to check whether the Dirac field belongs to H 1 loc or not in order to investigate the possibility of having weak solutions of Einstein equations. In other words, we have to check whether ∂ µ Ψ is locally square integrable which boils down to investigate the finiteness of the integral of quantity ∼ (r−r − ) 2(p−1) . Thus, the condition for Ψ to remain in H 1 loc then reduces to the following inequality The existence of weak solutions of Einstein's equation at Cauchy horizon is guaranteed by this condition which leads to a possible violation of strong cosmic censorship conjecture. In our study, we need to focus on the dominant mode contributions only which are the least damped modes of the quasinormal spectrum.

Strong Cosmic Censorship Conjecture for Kerr-de Sitter black holes in presence of Dirac field
In this section, we compute the quasi-normal modes of Kerr-de Sitter black hole in presence of Dirac field using numerical methods. But before that let us put the radial and angular equations given by Eq. (7) and Eq. (8) in standard Teukolsky form which read as follows [28] ∆ −s r ∂ r ∆ s+1 r ∂ r R + (r) + 4i(α + 1)rsω + 2(1 − α)s − 2αr 2 (s + 1)(2s + 1) where, x = cos θ, s = 1/2 and ξ = aω. Note that, in the non-rotating limit (a → 0), the angular Teukolsky equation Eq. (13) gives the separation constant as λ → l(l − 1) − s 2 + s. This can be used to define the angular eigen mode number l which satisfies the following relation, l ≥ max(|m|, |s|) [34]. These transformed equations allow us to use the method developed by [28], who showed that Eq. (13) can be transformed into Heun's equation which, in turn, give us a three-term recurrence relation for the angular  equation. This three term recurrence relation can be rewritten in terms of a continued fraction equation which we denote by P 1 (λ, ω) = 0. The radial Teukolsky equation Eq. (12) has five regular singularities at r + , r − , r c , −(r + + r − + r c ) and spatial infinity. The quasi-normal modes are defined as the eigen values of ω with Ψ satisfying the following boundary condition : there are only outgoing waves at the cosmological horizon r c and only ingoing waves at the event horizon r + . In order to satisfy the boundary condition, we write R + (r) as the multiplication of a function y(z) which is regular at the boundary and a factor which is divergent at r + and r c , i.e., where, dB/dr = V (r)/∆ r and z = (r − r + )(r c − r − )/(r − r − )(r c − r + ). By inserting Eq. (14) into Eq. (12) and expressing y(z) as a Frobenious series of the form ∞ n=0 a n z n , we get a seven-term recurrence relation. Due to convergence problem of seven term recurrence relation [35], it is better to reduce the seven-term recurrence relation into a three term recurrence relation using Gaussian elimination method [36,37]. Similar to the angular equation, this three term recurrence relation can also be written as an infinite continued fraction equation which we denote by P 2 (λ, ω) = 0. By solving the angular and radial continued fraction equations simultaneously, one gets the desired quasi-normal modes. However, instead of using the Gaussian elimination procedure to find the three-term recurrence relation for the radial equation, we have employed the Mathematica package developed by Jansen [38] to find the quasinormal modes. Here, we have solved the angular equation P 1 (λ, ω) = 0 and radial equation for y(z) iteratively, taking the approximate value of λ in Ref. [28] as the initial guess value. In order to understand the effect of Dirac particles on strong cosmic censorship conjecture, we need to look for modes for which the parameter β ≡ − Im(ω)/κ − becomes greater than 1/2. As stated earlier, we are interested in the dominant modes (least damped modes) of the quasinormal spectrum for a given value of angular eigen mode number l for the calculation of the parameter β. From Fig. 1, we can see that the dominant mode is always corresponds to modes with m = l for any values of mass scaled cosmological constant ΛM 2 [14,39]. Moreover, the imaginary part of the quasinormal frequency decreases with the increase of both the rotational parameter a and the cosmological constant Λ. In this regard, our result is fully consistent with the results presented in [39]. Moreover, near to the extremity, a rapid decrease of the value of Im(ω) is observed. The variation of the quantity Im(ω)/κ − with respect to the rotational parameter a for different values of ΛM 2 is presented in Fig. 2. Here, we consider only the least damped modes for a given value of angular eigen mode number i.e. the modes with m = l. It is interesting to see that for smaller values of mass scaled cosmological constant (ΛM 2 ≈ O(10 −3 )), modes corresponding to l = 1/2 dominate the quasinormal spectrum for certain values of a. However, for larger values of ΛM 2 , the eikonal modes (corresponding to the modes with large l value) become dominant. As evident from Fig. 2, strong cosmic censorship conjecture is respected for larger values of ΛM 2 . But for smaller values of mass scaled cosmological constant (ΛM 2 0.02), we are able to find a parameter space for which strong cosmic censorship gets violated. Moreover, in this parameter space, the violation is more severe for smaller values of ΛM 2 since β becomes greater than 1/2 for smaller values of rotational parameter a. So the assurance that strong cosmic censorship is always respected in astrophysical black holes [14,24] is no longer valid in the presence of Dirac fields.  Figure 2: The variation of Im(ω)/κ − as a function of (a/a max ) has been presented for different values of ΛM 2 . The parameter β corresponds to the least damped modes of the quasi-normal spectrum. Here, we have considered only the modes with m = l. In each of these plots, red, green and blue curves corresponds to value of Im(ω)/κ − for l = 0.5, l = 1.5 and l = 10.5 respectively. The vertical lines in the plots presented in upper-left and upper-right corner corresponds to the value of rotation parameter a for which strong cosmic censorship conjecture gets violated. As evident from the plots, for smaller values of ΛM 2 , we can find a parameter space which does not respect the the conjecture.

Conclusion
In recent years, several examples are found which suggest a breakdown of determinism in Reissner-Nordström-de Sitter black holes under the influence of several fundamental fields [18,19,[21][22][23][24][25][26]. Anyway, astrophysically meaningful Kerr-de Sitter black hole solutions seems to respect the conjecture [14,24]. However, the previous analyses were done considering only the effect of scalar fields. In this paper, we extended the study by considering Dirac fields in Kerr-de Sitter background. By considering the effect of linear perturbations on the spacetime metric of our interest, we found that there exist weak solutions of Einstein equation at the Cauchy horizon, if the parameter β defined by Eq. (11) become greater than 1/2 which leads to the violation of strong cosmic censorship conjecture. Comparing our result with Ref. [14], where the effect of scalar fields are considered, we see that the criteria for strong cosmic censorship violation remains the same in the presence of Dirac fields also. Once this criteria is obtained, we performed detailed numerical computation to find lowest lying quasinormal modes to determine the value of β as a function of rotational parameter a for different values of mass scaled cosmological constant ΛM 2 . In Fig. 1 and Fig. 2, we presented our main results. From these figures, it is clear that the lowest lying quasinormal modes always correspond to the m = l mode. Moreover, an increased value of rotational parameter results in modes with smaller decay rate. Near to the extremity, a rapid decrease of decay rate is observed. The decay rate also decreases with the increase of ΛM 2 . From Fig. 2, it is clear that the value of the parameter β always remains smaller than 1/2 for larger values of ΛM 2 . Hence, we can conclude that the strong cosmic censorship conjecture is always respected in "large" rotating black holes (black holes with ΛM 2 0.02). However, for smaller values of ΛM 2 , we found a parameter space where β becomes larger than 1/2. This, in turn, implies that in "smaller" black holes, Dirac fields can be smoothly extended beyond Cauchy horizon. Hence, in presence of Dirac fields, even the rotational parameter can not save the strong cosmic censorship conjecture for a certain parameter range. Moreover, for a fixed value of cosmological constant, black holes with smaller masses are more prone to violate the conjecture. Note that, in our work, we have considered the effect of linear perturbation only. It may be possible that the violation of strong cosmic censorship can be prevented by considering the non-linear or quantum gravitational effects. However, Cardoso et al. have shown that for Einstein-Maxwell-scalar field system, the non-linear perturbations can not save the conjecture [25]. Moreover, a recent study showed that the violation of this conjecture in 2 + 1 dimensional BTZ black holes even when quantum correction terms are added [40]. However, the studies in 3 + 1 dimensional black holes remain inconclusive. In particular, it has shown that the quantum correction can either amplify or suppress the blue shift instability in Reissner-Nordström black holes [41]. It would be interesting to see the those effects on Einstein-Dirac field system, which we leave for future work. Recently, Dafermos and Shlapentokh-Rothman [15] suggested an interesting proposal that the strong cosmic censorship can still be saved if one starts with rough initial data. This idea is further supported by Ref. [42], where the authors have studied a coupled gravitational and electro-magnetic perturbation and showed that in order to save the strong cosmic censorship, one must consider physically reasonable but slightly less smooth initial data. It will be interesting to see whether this proposal can save the conjecture in presence of Dirac field also. However, this is beyond the scope of this paper.