Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes

*Joás Venâncio and Carlos Batista*

#### **Abstract**

The exact electrically charged solutions to the Dirac equation in higherdimensional generalized Nariai spacetimes are obtained. Using these solutions, the boundary conditions leading to quasinormal modes of the Dirac field are analyzed, and their correspondent quasinormal frequencies are analytically calculated.

**Keywords:** quasinormal modes, generalized Nariai spacetimes, Dirac field, boundary conditions

#### **1. Introduction**

Quasinormal modes (QNMs) are eigenmodes of dissipative systems. For instance, if a spacetime with an event or cosmological horizon is perturbed from its equilibrium state, QNMs arise as damped oscillations with a spectrum of complex frequencies that do not depend on the details of the excitation. In fact, these frequencies depend just on the charges of the black hole, such as the mass, electric charge, and angular momentum [1, 2]. QNMs have been studied for a long time, and its interest has been renewed by the recent detection of gravitational waves, inasmuch as these are the modes that survive for a longer time when a background is perturbed and, therefore, these are the configurations that are generally measured by experiments [3–29]. Mathematically, this discrete spectrum of QNMs stems from the fact that certain boundary conditions must be imposed to the physical fields propagating in such background [30]. In this chapter, we consider a higher-dimensional generalization of the charged Nariai spacetime [31], namely, *dS*<sup>2</sup> *<sup>S</sup>*<sup>2</sup> … *<sup>S</sup>*<sup>2</sup> , and investigate the dynamics of perturbations of the electrically charged Dirac field (spin 1/2). In such a geometry, the spinorial formalism [32–34] is used to show that the Dirac equation is separable [35] and can be reduced to a Schrödinger-like equation [36] whose potential is contained in the Rosen-Morse class of integrable potentials, which has the so-called Pöschl-Teller potential as a particular case [37, 38]. Finally, the boundary conditions leading to QNMs are analyzed, and the quasinormal frequencies (QNFs) are analytically obtained [5, 39].

#### **2. Presenting the problem**

In *D* dimensions, the dynamics of general relativity in spacetimes with a cosmological constant Λ is described by the Einstein-Hilbert action<sup>1</sup>

$$S = \frac{1}{16\pi} \int d^D x \sqrt{|\mathbf{g}|} [\mathcal{R} - (D - 2)\Lambda] + \mathcal{S}\_m,\tag{1}$$

electrically charged spinorial field of mass *m* that obeys the Dirac equation minimally coupled to an electromagnetic field in such spacetime. In *D* ¼ 2*d*, this spacetime is formed from the direct product of the de Sitter space *dS*<sup>2</sup> with ð Þ *d* � 1 copies of the unit spheres *S*<sup>2</sup> possessing different radii *Rj*. Thus, the natural line element of the higher-dimensional version of the Nariai spacetime is given by

1

*f r*ð Þ *dr*<sup>2</sup> <sup>þ</sup><sup>X</sup>

*<sup>j</sup>* <sup>þ</sup> sin <sup>2</sup>

1 2 *Q*2 *d*

*<sup>j</sup> ,* (6)

*<sup>j</sup>* is the line element of the *j*th

*<sup>j</sup> :* (7)

*,* (8)

*j*¼2 *R*2 *<sup>j</sup> d*Ω<sup>2</sup>

*θ<sup>j</sup> dϕ*<sup>2</sup>

*<sup>j</sup>* <sup>þ</sup> *<sup>Q</sup>* 2ð Þ *D* � 2

*<sup>j</sup> :* (9)

*<sup>j</sup>* cos *θ<sup>j</sup> dϕj:* (10)

. Usually, the articles consider

� �*,* (11)

*R*1

� ��1*=*<sup>2</sup>

*dx<sup>ν</sup>* ¼ �*f r*ð Þ*dt*<sup>2</sup> <sup>þ</sup>

*, d*Ω<sup>2</sup>

where *Q*<sup>1</sup> is an electric charge and *Qj* are magnetic charges, while *Q* is

*<sup>Q</sup>* <sup>¼</sup> *<sup>Q</sup>*<sup>2</sup>

**<sup>A</sup>** <sup>¼</sup> *<sup>Q</sup>*<sup>1</sup> *rdt* <sup>þ</sup><sup>X</sup>

known as event horizons. The boundary conditions defining QNMs in our

*<sup>j</sup>* <sup>¼</sup> *<sup>d</sup>θ*<sup>2</sup>

*, Rj* ¼ Λ þ

<sup>1</sup> �<sup>X</sup> *d*

This spacetime is a locally static solution of Einstein's equation with a cosmological constant Λ and electromagnetic field *ℱ* ¼ *d***A** whose gauge field A is given by

*d*

*j*¼2 *Qj R*2

spacetime will be posed at these surfaces, as discussed in [39]. For this reason, the dependence of all the components of the field **Ψ** on the coordinates along the

*f r*ð Þ *dr* ) *<sup>x</sup>* <sup>¼</sup> *<sup>R</sup>*<sup>1</sup> arctanh *<sup>r</sup>*

that the coordinate *r* in de Sitter space assume values in the interval *r*∈ð Þ 0*, R*<sup>1</sup> [40–42]. However, this is just justified for de Sitter with *D* >2, but not for *D* ¼ 2; see [39] for more details. By this reason, our domain of interest will be *r*∈ð Þ �*R*1*, R*<sup>1</sup> . In such domain, it is useful to introduce the tortoise coordinate *x* defined by the

The coordinates in the metric are also called static, because they do not depend explicitly on the time coordinate *t*. One may notice that, in this coordinate system, this background has a local Killing vector *<sup>∂</sup><sup>t</sup>* whose norm vanishes at *<sup>r</sup>* ¼ �*R*1. Indeed, *r* ¼ �*R*<sup>1</sup> define closed null surfaces that surround the observer at all times,

*j*¼2

*Q*2

*ds*<sup>2</sup> <sup>¼</sup> *<sup>g</sup>*<sup>0</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

The radii *R*<sup>1</sup> and *Rj* are given by

2 *Q*2

unit sphere *S*<sup>2</sup> as follows

*<sup>R</sup>*<sup>1</sup> <sup>¼</sup> <sup>Λ</sup> � <sup>1</sup>

defined by

equation

**3**

*μνdx<sup>μ</sup>*

where *f r*ð Þ is a function of the coordinate *<sup>r</sup>* and *<sup>d</sup>*Ω<sup>2</sup>

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

*f r*ð Þ¼ <sup>1</sup> � *<sup>r</sup>*<sup>2</sup>

<sup>1</sup> <sup>þ</sup> *<sup>Q</sup>* 2ð Þ *D* � 2

Killing vector *∂<sup>t</sup>* is assumed to be of the form *e*�*iω<sup>t</sup>*

*dx* <sup>¼</sup> <sup>1</sup>

in terms of which the line element Eq. (6) becomes

� ��1*=*<sup>2</sup>

*R*2 1

where R is the Ricci scalar and *Sm* stands for the action of all matter fields f g **Φ***<sup>i</sup>* coupled to gravity appearing in the theory, which can be scalar, spinorial, vectorial, and so on. The least action principle allows to find the equations of motion for the fields *gμν* and **Φ***<sup>i</sup>* which are given, respectively, by

$$\mathcal{R}\_{\mu\nu} - \frac{1}{2}\mathcal{R}\mathbf{g}\_{\mu\nu} + \frac{(D-2)}{2}\Lambda \mathbf{g}\_{\mu\nu} = 8\pi T\_{\mu\nu}, \quad \frac{\delta \mathcal{S}\_m}{\delta \Phi\_i} = 0,\tag{2}$$

where T *μν* is the symmetric stress-energy tensor associated to **Φ***<sup>i</sup>* defined by the equation

$$T^{\mu\nu} = \frac{2}{\sqrt{|\mathbf{g}|}} \frac{\delta \mathbf{S}\_m}{\delta \mathbf{g}\_{\mu\nu}}.\tag{3}$$

Since any symmetry has been imposed, the general solution of the system of Eq. (2) is some metric and fields in the background this metric

$$d\mathbf{s}^2 = \mathbf{g}\_{\mu\nu}(\mathbf{x})d\mathbf{x}^\mu d\mathbf{x}^\nu, \quad \Phi\_i = \Phi\_i(\mathbf{x}).\tag{4}$$

Now, let the pair *g*ð Þ <sup>0</sup> *μν* and **<sup>Φ</sup>**ð Þ <sup>0</sup> *<sup>i</sup>* be a solution for the equations of motion Eq. (2). Then, in order to study the perturbations around this particular solution, we write our fields as a sum of the unperturbed fields *g*ð Þ <sup>0</sup> *μν* and **<sup>Φ</sup>**ð Þ <sup>0</sup> *<sup>i</sup>* and the small perturbations *hμν* and **Ψ***<sup>i</sup>*

$$\mathbf{g}\_{\mu\nu} = \mathbf{g}\_{\mu\nu}^{(0)} + h\_{\mu\nu} \quad \Phi\_i = \Phi\_i^{(0)} + \Psi\_i \tag{5}$$

where by "small" we mean that we neglect the quadratic and higher-order powers of the perturbation fields. Inserting the above equation into Eq. (2), we are left with a set of linear equations satisfied by the perturbed fields *hμν* and **Ψ***<sup>i</sup>* . In general, these equations are coupled, namely, **Ψ***<sup>i</sup>* is a source for *hμν* and vice versa. However, in the special case in which **Φ**ð Þ <sup>0</sup> *<sup>i</sup>* ¼ 0, the equations governing the perturbed fields **Ψ***<sup>i</sup>* can be decoupled from the metric perturbation *hμν* and vice versa. The reason why this happen is that when **Φ**ð Þ <sup>0</sup> *<sup>i</sup>* ¼ 0, the stress-energy tensor T *μν* can be set to zero at first order in the perturbation, since T *μν* is typically quadratic or of higher order in the matter fields and, therefore, can be neglected. Therefore, investigating the linear dynamics of generic small perturbations of the matter fields with T *μν* ¼ 0 is equivalent to studying the test fields **Ψ***<sup>i</sup>* in the background *g*ð Þ <sup>0</sup> *μν* .

In what follows, let us consider a specific matter field **Ψ** propagating in a generalized version of the Nariai spacetime described in Ref. [31]. Here, **Ψ** is an

<sup>1</sup> The coefficient of Λ in *S* can be chosen of several manners. In particular, for any dimension *D*, in order to insure that the pure dS or pure AdS spacetimes are described by *gtt* <sup>¼</sup> <sup>1</sup> � ð Þ <sup>Λ</sup>*=*<sup>3</sup> *<sup>r</sup>*2, as occurs in the case *D* ¼ 4, this coefficient should be ð Þ *D* � 1 ð Þ *D* � 2 .

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes DOI: http://dx.doi.org/10.5772/intechopen.89179*

electrically charged spinorial field of mass *m* that obeys the Dirac equation minimally coupled to an electromagnetic field in such spacetime. In *D* ¼ 2*d*, this spacetime is formed from the direct product of the de Sitter space *dS*<sup>2</sup> with ð Þ *d* � 1 copies of the unit spheres *S*<sup>2</sup> possessing different radii *Rj*. Thus, the natural line element of the higher-dimensional version of the Nariai spacetime is given by

$$ds^2 = \mathbf{g}\_{\mu\nu}^0 d\mathbf{x}^\mu d\mathbf{x}^\nu = -f(r)dt^2 + \frac{1}{f(r)}dr^2 + \sum\_{j=2}^d R\_j^2 d\Omega\_j^2,\tag{6}$$

where *f r*ð Þ is a function of the coordinate *<sup>r</sup>* and *<sup>d</sup>*Ω<sup>2</sup> *<sup>j</sup>* is the line element of the *j*th unit sphere *S*<sup>2</sup> as follows

$$f(r) = 1 - \frac{r^2}{R\_1^2}, \quad d\Omega\_j^2 = d\theta\_j^2 + \sin^2\theta\_j d\phi\_j^2. \tag{7}$$

The radii *R*<sup>1</sup> and *Rj* are given by

**2. Presenting the problem**

*Progress in Relativity*

In *D* dimensions, the dynamics of general relativity in spacetimes with a

ffiffiffiffiffi j j *g* q

where R is the Ricci scalar and *Sm* stands for the action of all matter fields f g **Φ***<sup>i</sup>* coupled to gravity appearing in the theory, which can be scalar, spinorial, vectorial, and so on. The least action principle allows to find the equations of motion for the

where T *μν* is the symmetric stress-energy tensor associated to **Φ***<sup>i</sup>* defined by the

ffiffiffiffiffi j j *<sup>g</sup>* <sup>p</sup>

Since any symmetry has been imposed, the general solution of the system of

*dx<sup>ν</sup>*

Then, in order to study the perturbations around this particular solution, we write

*μν* <sup>þ</sup> *<sup>h</sup>μν,* **<sup>Φ</sup>***<sup>i</sup>* <sup>¼</sup> **<sup>Φ</sup>**ð Þ <sup>0</sup>

be decoupled from the metric perturbation *hμν* and vice versa. The reason why this

order in the perturbation, since T *μν* is typically quadratic or of higher order in the matter fields and, therefore, can be neglected. Therefore, investigating the linear dynamics of generic small perturbations of the matter fields with T *μν* ¼ 0 is equiva-

In what follows, let us consider a specific matter field **Ψ** propagating in a generalized version of the Nariai spacetime described in Ref. [31]. Here, **Ψ** is an

<sup>1</sup> The coefficient of Λ in *S* can be chosen of several manners. In particular, for any dimension *D*, in order to insure that the pure dS or pure AdS spacetimes are described by *gtt* <sup>¼</sup> <sup>1</sup> � ð Þ <sup>Λ</sup>*=*<sup>3</sup> *<sup>r</sup>*2, as occurs in the case

where by "small" we mean that we neglect the quadratic and higher-order powers of the perturbation fields. Inserting the above equation into Eq. (2), we are left with a set of linear equations satisfied by the perturbed fields *hμν* and **Ψ***<sup>i</sup>* . In general, these equations are coupled, namely, **Ψ***<sup>i</sup>* is a source for *hμν* and vice versa. However, in the

*δSm δgμν*

<sup>T</sup> *μν* <sup>¼</sup> <sup>2</sup>

<sup>Λ</sup>*gμν* <sup>¼</sup> <sup>8</sup>*π*<sup>T</sup> *μν, <sup>δ</sup>Sm*

½R � ð Þ *D* � 2 Λ� þ *Sm,* (1)

*δ***Φ***<sup>i</sup>*

*:* (3)

*,* **Φ***<sup>i</sup>* ¼ **Φ***i*ð Þ *x :* (4)

*<sup>i</sup>* and the small perturba-

*<sup>i</sup>* þ **Ψ***i,* (5)

*<sup>i</sup>* be a solution for the equations of motion Eq. (2).

*μν* and **<sup>Φ</sup>**ð Þ <sup>0</sup>

*<sup>i</sup>* ¼ 0, the equations governing the perturbed fields **Ψ***<sup>i</sup>* can

*<sup>i</sup>* ¼ 0, the stress-energy tensor T *μν* can be set to zero at first

*μν* .

¼ 0*,* (2)

cosmological constant Λ is described by the Einstein-Hilbert action<sup>1</sup>

<sup>R</sup>*gμν* <sup>þ</sup> ð Þ *<sup>D</sup>* � <sup>2</sup> 2

Eq. (2) is some metric and fields in the background this metric

*ds*<sup>2</sup> <sup>¼</sup> *<sup>g</sup>μν*ð Þ *<sup>x</sup> dx<sup>μ</sup>*

*μν* and **<sup>Φ</sup>**ð Þ <sup>0</sup>

*<sup>g</sup>μν* <sup>¼</sup> *<sup>g</sup>*ð Þ <sup>0</sup>

lent to studying the test fields **Ψ***<sup>i</sup>* in the background *g*ð Þ <sup>0</sup>

*D* ¼ 4, this coefficient should be ð Þ *D* � 1 ð Þ *D* � 2 .

our fields as a sum of the unperturbed fields *g*ð Þ <sup>0</sup>

*<sup>S</sup>* <sup>¼</sup> <sup>1</sup> 16*π* ð *dDx*

fields *gμν* and **Φ***<sup>i</sup>* which are given, respectively, by

<sup>R</sup>*μν* � <sup>1</sup> 2

Now, let the pair *g*ð Þ <sup>0</sup>

special case in which **Φ**ð Þ <sup>0</sup>

happen is that when **Φ**ð Þ <sup>0</sup>

**2**

tions *hμν* and **Ψ***<sup>i</sup>*

equation

$$R\_1 = \left[\Lambda - \frac{1}{2}Q\_1^2 + \frac{Q}{2(D-2)}\right]^{-1/2},\ R\_j = \left[\Lambda + \frac{1}{2}Q\_j^2 + \frac{Q}{2(D-2)}\right]^{-1/2},\tag{8}$$

where *Q*<sup>1</sup> is an electric charge and *Qj* are magnetic charges, while *Q* is defined by

$$Q = Q\_1^2 - \sum\_{j=2}^d Q\_j^2. \tag{9}$$

This spacetime is a locally static solution of Einstein's equation with a cosmological constant Λ and electromagnetic field *ℱ* ¼ *d***A** whose gauge field A is given by

$$\mathcal{A} = Q\_1 r dt + \sum\_{j=2}^{d} Q\_j R\_j^2 \cos \theta\_j d\phi\_j. \tag{10}$$

The coordinates in the metric are also called static, because they do not depend explicitly on the time coordinate *t*. One may notice that, in this coordinate system, this background has a local Killing vector *<sup>∂</sup><sup>t</sup>* whose norm vanishes at *<sup>r</sup>* ¼ �*R*1. Indeed, *r* ¼ �*R*<sup>1</sup> define closed null surfaces that surround the observer at all times, known as event horizons. The boundary conditions defining QNMs in our spacetime will be posed at these surfaces, as discussed in [39]. For this reason, the dependence of all the components of the field **Ψ** on the coordinates along the Killing vector *∂<sup>t</sup>* is assumed to be of the form *e*�*iω<sup>t</sup>* . Usually, the articles consider that the coordinate *r* in de Sitter space assume values in the interval *r*∈ð Þ 0*, R*<sup>1</sup> [40–42]. However, this is just justified for de Sitter with *D* >2, but not for *D* ¼ 2; see [39] for more details. By this reason, our domain of interest will be *r*∈ð Þ �*R*1*, R*<sup>1</sup> . In such domain, it is useful to introduce the tortoise coordinate *x* defined by the equation

$$d\infty = \frac{1}{f(r)} dr \Rightarrow \mathfrak{x} = R\_1 \text{arctanh}\left(\frac{r}{R\_1}\right),\tag{11}$$

in terms of which the line element Eq. (6) becomes

#### **Figure 1.**

*Illustration of the boundary condition associated to QNMs in our spacetime. The wavy arrows represent the direction of the perturbation field at the boundaries r* ¼ �*R*1*, while the cones are the local light cones. Mathematically, the wavy arrow pointing to the right represents e*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup> , while the wavy arrow pointing to the left represents e*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup> . For more details, see Ref. [39].*

$$ds^2 = \frac{1}{\cosh^2(\mathfrak{x}/R\_1)} \left(-dt^2 + d\mathfrak{x}^2\right) + \sum\_{j=2}^d R\_j^2 d\Omega\_j^2,\tag{12}$$

A field of spin 1*=*2 with electric charge *q* and mass *m* propagating in such spacetime

where A*<sup>α</sup>* stands for the components of the background gauge field. In *D* ¼ 2*d* dimensions, the Dirac matrices <sup>Γ</sup>*<sup>α</sup>* represent faithfully the Clifford algebra by 2*<sup>d</sup>* �

with *<sup>d</sup>* standing for the 2*<sup>d</sup>* � <sup>2</sup>*<sup>d</sup>* identity matrix. The index *<sup>α</sup>, <sup>β</sup>, <sup>γ</sup>* run from 1 to 2*<sup>d</sup>* and label the vector fields of an orthonormal frame f g *e<sup>α</sup>* . In order to solve the Dirac equation, we must introduce a suitable orthonormal frame of vector fields, which in

Γ*<sup>α</sup>* Γ*<sup>β</sup>* þ Γ*<sup>β</sup>* Γ*<sup>α</sup>* ¼ *g eα, e<sup>β</sup>*

*<sup>e</sup>*<sup>1</sup> ¼ �*<sup>i</sup>* cosh ð Þ *<sup>x</sup>=R*<sup>1</sup> *<sup>∂</sup><sup>t</sup> , <sup>e</sup><sup>j</sup>* <sup>¼</sup> <sup>1</sup>

*<sup>e</sup>*<sup>~</sup><sup>1</sup> <sup>¼</sup> cosh ð Þ *<sup>x</sup>=R*<sup>1</sup> *<sup>∂</sup><sup>x</sup> , <sup>e</sup>*<sup>~</sup>*<sup>j</sup>* <sup>¼</sup> <sup>1</sup>

where the index *j* ranges from 2 to *d*. In particular, note that

� � <sup>¼</sup> *δαβ* \$

*g eα, e<sup>β</sup>*

spin connection according to the following relation:

and lowered unpunished. In particular, *ωβγ*

*<sup>σ</sup>*<sup>1</sup> <sup>¼</sup> 0 1

is, then, given by

**5**

<sup>Γ</sup>*<sup>α</sup>*ð Þ <sup>∇</sup>*<sup>α</sup>* � *iq*A*<sup>α</sup>* **<sup>Ψ</sup>** <sup>¼</sup> *<sup>m</sup>***Ψ***,* (14)

*Rj* sin *θ<sup>j</sup>*

*Rj ∂θj ,*

*g ea, e*<sup>~</sup> *b* � � <sup>¼</sup> <sup>0</sup>*,*

8 >>><

>>>:

where *a* and *a*~ are indices that range from 1 to *d*. The index *a* labels the first *d* vector fields of the orthonormal frame f g *e<sup>a</sup>* , while the index *a*~ labels the remaining *d* vectors of the frame f g *e<sup>a</sup>* . The derivatives of the frame vector fields determine the

<sup>∇</sup>*<sup>α</sup> <sup>e</sup><sup>β</sup>* <sup>¼</sup> *<sup>ω</sup> <sup>γ</sup>*

Since the metric *g* is a covariantly constant tensor, it follows that the coefficients of the spin connection with all low indices *ωαβγ* <sup>¼</sup> *<sup>ω</sup> <sup>ε</sup> αβ δεγ* are antisymmetric in their two last indices, *ωαβγ* ¼ �*ωαγβ*. Note that the indices of the spin connection are raised and lowered with *δαβ* and *δαβ*, respectively, so that frame indices can be raised

square brackets are antisymmetrized. The covariant derivative of a spinorial field **Ψ**

Our aim is to separate the Dirac Eq. (14). In order to accomplish this, it is necessary to use a suitable representation for the Dirac matrices. We recall that

*i* 0

� �*, <sup>σ</sup>*<sup>3</sup> <sup>¼</sup> 1 0

<sup>∇</sup>*α***<sup>Ψ</sup>** <sup>¼</sup> *<sup>∂</sup>α***<sup>Ψ</sup>** � <sup>1</sup>

with *∂<sup>α</sup>* denoting the partial derivative along the vector field *eα*.

1 0 � �*, <sup>σ</sup>*<sup>2</sup> <sup>¼</sup> <sup>0</sup> �*<sup>i</sup>*

*<sup>α</sup>* <sup>¼</sup> *<sup>ω</sup>*½ � *βγ*

*g ea*<sup>~</sup>*, e*<sup>~</sup> *b* � � <sup>¼</sup> *<sup>δ</sup>a*~~

*g e*ð Þ¼ *a, e<sup>b</sup> δab,*

*∂ϕj ,*

*b,*

*αβ e<sup>γ</sup> :* (18)

*<sup>α</sup>* , where indices inside the

<sup>4</sup> *<sup>ω</sup> βγ <sup>α</sup>* <sup>Γ</sup>*<sup>β</sup>* <sup>Γ</sup>*<sup>γ</sup>* **<sup>Ψ</sup>***,* (19)

0 �1

� �*,* (20)

(16)

(17)

� �*d,* (15)

is a spinorial field obeying the following version of the Dirac equation:

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

2*<sup>d</sup>* matrices obeying the relation

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

the case of our background is given by

and the gauge field can be rewritten as

$$\mathcal{A} = Q\_1 R\_1 \tanh\left(\varkappa / R\_1\right) dt + \sum\_{j=2}^d Q\_j R\_j^2 \cos\theta\_j d\phi\_j. \tag{13}$$

In particular, note that the tortoise coordinate maps the domain between two horizons, *r*∈ð Þ �*R*1*, R*<sup>1</sup> , into the interval *x*∈ð Þ �∞*,* ∞ .

The QNMs accounting for an important class of fields are associated to **Ψ** which are solutions to the equations of motion that satisfy specific boundary conditions imposed at the horizons of the spacetime in which the field is propagating; see [5, 6, 43, 44] for more details. In this chapter, we will use the boundary conditions as illustrated in **Figure 1**.

From the mathematical of view, since we are assuming that the time dependence of **Ψ** is *e*�*iω<sup>t</sup>* , this boundary condition means that near the horizons *r* ¼ �*R*1, that is, as *<sup>x</sup>* ! �∞, the radial component of the field **<sup>Ψ</sup>** should behave as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup>* at *<sup>x</sup>* ! <sup>∞</sup>, while it should go as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* at *<sup>x</sup>* ! �∞. The eigenfrequencies of this problem are complex, the reason why they are called QNFs. The real part of the QNFs is associated with the oscillation frequencies of the signal, while the imaginary part is related to its decay in time. This decay in time is closely related to the fact that the event horizon has a dissipative nature.

One interesting feature of this spacetime is that we can compute exactly the QNMs. The exactly solvable systems are usually limits of more realistic systems and allow us to study in detail some properties of a physical process and test some methods which can be used to analyze more complicated systems. Thus they are powerful tools in many research lines. Therefore we expect that the exactly computed QNFs for *D*-dimensional generalized Nariai spacetime may play an important role in future research [27].

#### **3. Dirac equation in** *D***-dimensional generalized Nariai spacetime**

Let us present the construction of a solution to the Dirac equation minimally coupled to the electromagnetic field of *D*-dimensional generalized Nariai spacetime. A field of spin 1*=*2 with electric charge *q* and mass *m* propagating in such spacetime is a spinorial field obeying the following version of the Dirac equation:

$$
\Gamma^a(\nabla\_a - iq\mathcal{A}\_a)\Psi = m\Psi,\tag{14}
$$

where A*<sup>α</sup>* stands for the components of the background gauge field. In *D* ¼ 2*d* dimensions, the Dirac matrices <sup>Γ</sup>*<sup>α</sup>* represent faithfully the Clifford algebra by 2*<sup>d</sup>* � 2*<sup>d</sup>* matrices obeying the relation

$$
\Gamma\_a \Gamma\_\beta + \Gamma\_\beta \Gamma\_a = \mathbf{g} \begin{pmatrix} \mathbf{e}\_a \ \mathbf{e}\_\beta \end{pmatrix} \mathbb{I}\_{d\nu} \tag{15}
$$

with *<sup>d</sup>* standing for the 2*<sup>d</sup>* � <sup>2</sup>*<sup>d</sup>* identity matrix. The index *<sup>α</sup>, <sup>β</sup>, <sup>γ</sup>* run from 1 to 2*<sup>d</sup>* and label the vector fields of an orthonormal frame f g *e<sup>α</sup>* . In order to solve the Dirac equation, we must introduce a suitable orthonormal frame of vector fields, which in the case of our background is given by

$$\begin{aligned} \mathbf{e}\_1 &= -i \cosh \left( \mathbf{x} / R\_1 \right) \boldsymbol{\vartheta}\_l \ , \ \mathbf{e}\_j = \frac{1}{R\_j \sin \theta\_j} \ \boldsymbol{\vartheta}\_{\phi\_j} \end{aligned} \tag{16}$$
 
$$\begin{aligned} \mathbf{e}\_1 &= \cosh \left( \mathbf{x} / R\_1 \right) \boldsymbol{\vartheta}\_{\mathbf{x}} \ , \ \mathbf{e}\_j = \frac{1}{R\_j} \ \boldsymbol{\vartheta}\_{\theta\_j} \end{aligned} \tag{16}$$

where the index *j* ranges from 2 to *d*. In particular, note that

$$\mathbf{g}\left(\mathbf{e}\_{a},\mathbf{e}\_{\beta}\right) = \delta\_{a\beta} \leftrightarrow \begin{cases} \mathbf{g}\left(\mathbf{e}\_{a},\mathbf{e}\_{b}\right) = \delta\_{ab}, \\ \mathbf{g}\left(\mathbf{e}\_{a},\mathbf{e}\_{\dot{b}}\right) = 0, \\ \mathbf{g}\left(\mathbf{e}\_{\dot{a}},\mathbf{e}\_{\dot{b}}\right) = \delta\_{\dot{a}\dot{b}}, \end{cases} \tag{17}$$

where *a* and *a*~ are indices that range from 1 to *d*. The index *a* labels the first *d* vector fields of the orthonormal frame f g *e<sup>a</sup>* , while the index *a*~ labels the remaining *d* vectors of the frame f g *e<sup>a</sup>* . The derivatives of the frame vector fields determine the spin connection according to the following relation:

$$
\nabla\_a \mathbf{e}\_\beta = \boldsymbol{\alpha}\_{a\beta}{}^\gamma \mathbf{e}\_\gamma. \tag{18}
$$

Since the metric *g* is a covariantly constant tensor, it follows that the coefficients of the spin connection with all low indices *ωαβγ* <sup>¼</sup> *<sup>ω</sup> <sup>ε</sup> αβ δεγ* are antisymmetric in their two last indices, *ωαβγ* ¼ �*ωαγβ*. Note that the indices of the spin connection are raised and lowered with *δαβ* and *δαβ*, respectively, so that frame indices can be raised and lowered unpunished. In particular, *ωβγ <sup>α</sup>* <sup>¼</sup> *<sup>ω</sup>*½ � *βγ <sup>α</sup>* , where indices inside the square brackets are antisymmetrized. The covariant derivative of a spinorial field **Ψ** is, then, given by

$$
\nabla\_a \Psi = \partial\_a \Psi - \frac{1}{4} \alpha\_a \, \, ^{\beta \gamma} \Gamma\_\beta \Gamma\_\gamma \Psi,\tag{19}
$$

with *∂<sup>α</sup>* denoting the partial derivative along the vector field *eα*.

Our aim is to separate the Dirac Eq. (14). In order to accomplish this, it is necessary to use a suitable representation for the Dirac matrices. We recall that

$$
\sigma\_1 = \begin{bmatrix} \mathbf{0} & \mathbf{1} \\ \mathbf{1} & \mathbf{0} \end{bmatrix}, \quad \sigma\_2 = \begin{bmatrix} \mathbf{0} & -i \\ i & \mathbf{0} \end{bmatrix}, \quad \sigma\_3 = \begin{bmatrix} \mathbf{1} & \mathbf{0} \\ \mathbf{0} & -\mathbf{1} \end{bmatrix}, \tag{20}
$$

*ds*<sup>2</sup> <sup>¼</sup> <sup>1</sup> cosh <sup>2</sup>

horizons, *r*∈ð Þ �*R*1*, R*<sup>1</sup> , into the interval *x*∈ð Þ �∞*,* ∞ .

and the gauge field can be rewritten as

*left represents e*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup> . For more details, see Ref. [39].*

event horizon has a dissipative nature.

of **Ψ** is *e*�*iω<sup>t</sup>*

**4**

**Figure 1.**

*Progress in Relativity*

ð Þ *x=R*<sup>1</sup>

**<sup>A</sup>** <sup>¼</sup> *<sup>Q</sup>*1*R*<sup>1</sup> tanh ð Þ *<sup>x</sup>=R*<sup>1</sup> *dt* <sup>þ</sup><sup>X</sup>

�*dt*<sup>2</sup> <sup>þ</sup> *dx*<sup>2</sup> � � <sup>þ</sup><sup>X</sup>

*d*

*QjR*<sup>2</sup>

*j*¼2

, this boundary condition means that near the horizons *r* ¼ �*R*1, that is,

In particular, note that the tortoise coordinate maps the domain between two

*Illustration of the boundary condition associated to QNMs in our spacetime. The wavy arrows represent the direction of the perturbation field at the boundaries r* ¼ �*R*1*, while the cones are the local light cones. Mathematically, the wavy arrow pointing to the right represents e*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup> , while the wavy arrow pointing to the*

The QNMs accounting for an important class of fields are associated to **Ψ** which are solutions to the equations of motion that satisfy specific boundary conditions imposed at the horizons of the spacetime in which the field is propagating; see [5, 6, 43, 44] for more details. In this chapter, we will use the boundary conditions as illustrated in **Figure 1**. From the mathematical of view, since we are assuming that the time dependence

as *<sup>x</sup>* ! �∞, the radial component of the field **<sup>Ψ</sup>** should behave as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup>* at *<sup>x</sup>* ! <sup>∞</sup>, while it should go as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* at *<sup>x</sup>* ! �∞. The eigenfrequencies of this problem are complex, the reason why they are called QNFs. The real part of the QNFs is associated with the oscillation frequencies of the signal, while the imaginary part is related to its decay in time. This decay in time is closely related to the fact that the

One interesting feature of this spacetime is that we can compute exactly the QNMs. The exactly solvable systems are usually limits of more realistic systems and allow us to study in detail some properties of a physical process and test some methods which can be used to analyze more complicated systems. Thus they are powerful tools in many research lines. Therefore we expect that the exactly computed QNFs for *D*-dimensional generalized Nariai spacetime may play an important role in future research [27].

**3. Dirac equation in** *D***-dimensional generalized Nariai spacetime**

Let us present the construction of a solution to the Dirac equation minimally coupled to the electromagnetic field of *D*-dimensional generalized Nariai spacetime.

*d*

*j*¼2 *R*2 *<sup>j</sup> d*Ω<sup>2</sup>

*<sup>j</sup> ,* (12)

*<sup>j</sup>* cos *θ<sup>j</sup> dϕj:* (13)

are the Hermitian Pauli matrices and denote the 2 � 2 identity matrix. Using this notation, a convenient representation of the Dirac matrices is the following:

$$\begin{aligned} \Gamma\_{a} &= \underbrace{\sigma\_{3} \otimes \dots \otimes}\_{(a-1)\text{ times}} \sigma\_{1} \underbrace{\otimes \mathbb{I} \otimes \dots \otimes \mathbb{I}}\_{(d-a)\text{ times}} \\ \Gamma\_{\vec{a}} &= \underbrace{\sigma\_{3} \otimes \dots \otimes \sigma\_{2} \otimes \mathbb{I} \otimes \dots \otimes \mathbb{I}}\_{(a-1)\text{ times}} \end{aligned} \tag{21}$$

<sup>Γ</sup>*a***<sup>Ψ</sup>** <sup>¼</sup> <sup>X</sup>

Analogously, we have:

<sup>Γ</sup>*a*<sup>~</sup>**<sup>Ψ</sup>** <sup>¼</sup> <sup>X</sup>

f g*s*

… ⊗ *ξ*

decomposed in the form

*<sup>∂</sup>*~<sup>1</sup> <sup>þ</sup> *<sup>ω</sup>*1~<sup>11</sup>

following relation

**7**

f g*s*

… ⊗ *ξ*

… ⊗ *ξsa*�<sup>1</sup> ⊗ *ξ*

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

*sd* ¼ �*<sup>i</sup>*

… ⊗ *ξsa*�<sup>1</sup> ⊗ *ξ*

below or consult more details in [39]):

*L* ¼

*λ*2 <sup>2</sup> <sup>þ</sup> *<sup>λ</sup>*<sup>2</sup>

8 >>><

>>>:

*ωj*

q

ð Þ *<sup>s</sup>*1*s*2…*sa*�<sup>1</sup> **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sd <sup>ξ</sup>s*<sup>1</sup> <sup>⊗</sup> *<sup>ξ</sup>*

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

*sa* ⊗ *ξ*

ð Þ *<sup>s</sup>*1*s*2…*sa*�<sup>1</sup> ð Þ *isa* **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sd <sup>ξ</sup>s*<sup>1</sup> <sup>⊗</sup> *<sup>ξ</sup>*

*sa*þ<sup>1</sup> ⊗ … ⊗ *ξ*

<sup>1</sup> ð Þ *t, x* <sup>Ψ</sup>*<sup>s</sup>*<sup>2</sup>

following differential equation (the reader is invited to demonstrate the equation

<sup>2</sup> � *iq*A<sup>1</sup>

by the eigenvalues *λ<sup>j</sup>* of the Dirac operator on unit sphere *S*<sup>2</sup> according to the

*d*

as demonstrated in Appendix A of Ref. [39]. In our frame of vectors, the only

*Rj*

*R*1

cot *θj,*

sinh ð Þ *x=R*<sup>1</sup> *,*

The separation constant *L* in the above equation depends on the angular modes.

where each index *sa* can take the values *sa* ¼ �1, the fields <sup>Ψ</sup>*<sup>s</sup>*<sup>1</sup>

In particular, in the special case of vanishing magnetic charges *Qj*

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

<sup>3</sup> <sup>þ</sup> … <sup>þ</sup> *<sup>λ</sup>*<sup>2</sup>

components of the spin connection that are potentially nonvanishing are

*<sup>ω</sup>*1~<sup>11</sup> ¼ �*ω*11~<sup>1</sup> ¼ � <sup>1</sup>

<sup>~</sup>*jj* ¼ �*ωjj*~*<sup>j</sup>* <sup>¼</sup> <sup>1</sup>

X f g*s*

*sa* ⊗ *ξ*

<sup>Ψ</sup>*<sup>s</sup>*1*s*2…*sd* <sup>¼</sup> <sup>Ψ</sup>*<sup>s</sup>*<sup>1</sup>

<sup>2</sup> � *iq*A<sup>~</sup><sup>1</sup> � *is*<sup>1</sup> *<sup>∂</sup>*<sup>1</sup> <sup>þ</sup> *<sup>ω</sup>*<sup>~</sup>11~<sup>1</sup>

h i � � <sup>Ψ</sup>*<sup>s</sup>*<sup>1</sup>

*sd* <sup>¼</sup> <sup>X</sup> f g*s*

*<sup>s</sup>*<sup>2</sup> ⊗ … ⊗ *ξ*

ð Þ *<sup>s</sup>*1*s*2…*sa* **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sa*�1ð Þ �*sa sa*þ1…*sd <sup>ξ</sup>*

where from the first to the second line we have changed the index *sa* to �*sa*, which does not change the final result, since we are summing over all values of *sa*, which comprise the same list of the values of �*sa*. Moreover, we have used ð Þ *sa* <sup>2</sup> <sup>¼</sup> 1.

*sd ,*

*<sup>s</sup>*<sup>2</sup> ⊗ … ⊗ *ξ*

ð Þ *<sup>s</sup>*1*s*2…*sa sa* **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sa*�1ð Þ �*sa sa*þ1…*sd <sup>ξ</sup><sup>s</sup>*<sup>1</sup> <sup>⊗</sup> *<sup>ξ</sup>*

<sup>2</sup> ð Þ <sup>Φ</sup>2*, <sup>θ</sup>*<sup>2</sup> …Ψ*sd*

*sd :*

All that was seen above are necessary tools to attack our initial problem of separating the general Eq. (14). In order to solve such an equation, we need to separate the degrees of freedom of the field, which can be quite challenging in general. Fortunately, the spacetime considered here is the direct product of twodimensional spaces of constant curvature, which is exactly the class of spaces studied in Ref. [39]. Indeed, in this latter paper, it is shown that the Dirac equation minimally coupled to an electromagnetic field is separable in such backgrounds. In particular, assuming that the components of the spinor field Eq. (24) can be

*sa*þ<sup>1</sup> ⊗ … ⊗ *ξ*

*sa*�<sup>1</sup> ⊗ *ξ*

�*sa* ⊗ *ξ*

*<sup>s</sup>*<sup>1</sup> ⊗ *ξ <sup>s</sup>*<sup>2</sup> ⊗

*sa*�<sup>1</sup> ⊗ *ξ*

�*sa* ⊗ *ξ*

*<sup>s</sup>*<sup>2</sup> ⊗

*<sup>d</sup>* ð Þ Φ*d, θ<sup>d</sup> ,* (29)

<sup>1</sup> <sup>¼</sup> ð Þ *<sup>L</sup>* � *i s*1*<sup>m</sup>* <sup>Ψ</sup>�*s*<sup>1</sup>

*, λ<sup>j</sup>* ¼ �1*,* � 2*,* � 3*,* …*,* (31)

<sup>1</sup> ð Þ *t, x* satisfy the

<sup>1</sup> *:* (30)

, it is determined

(32)

*sa*þ<sup>1</sup> ⊗

*sa*þ<sup>1</sup> ⊗

(27)

(28)

where stands for the 2 � 2 identity matrix. Indeed, we can easily check that the Clifford algebra given in Eq. (15) is properly satisfied by the above matrices.<sup>2</sup> In this case, spinorial fields are represented by the column vectors on which these matrices act. We can introduce a basis of this representation by the direct products of spinors *ξ <sup>s</sup>* given by

$$\mathfrak{F}^+ = \begin{bmatrix} \mathbf{1} & \\ \mathbf{0} & \end{bmatrix}, \quad \mathfrak{F}^- = \begin{bmatrix} \mathbf{0} & \\ \mathbf{1} & \end{bmatrix}, \tag{22}$$

which, under the action of the Pauli matrices, satisfy concisely the relations

$$
\sigma\_1 \mathfrak{F}' = \mathfrak{F}^{-\prime}, \quad \sigma\_2 \mathfrak{F}' = i s \mathfrak{F}^{-\prime}, \quad \sigma\_3 \mathfrak{F}' = s \mathfrak{F}'.\tag{23}
$$

Indeed, in *<sup>D</sup>* <sup>¼</sup> <sup>2</sup>*<sup>d</sup>* dimensions, a general spinor field has 2*<sup>d</sup>* degrees of freedom and can be written as

$$\Psi = \sum\_{\{\epsilon\}} \Psi^{\epsilon\_1 \epsilon\_2 \dots \epsilon\_d} \xi^{\epsilon\_1} \otimes \xi^{\epsilon\_2} \otimes \dots \otimes \xi^{\epsilon\_d},\tag{24}$$

where each of the indices *sa* can take the values "+1" and "�1." Since every *sa* can take just two values, it follows that the sum over f g*<sup>s</sup>* � f g *<sup>s</sup>*1*, s*2*,* …*, sd* comprises 2*<sup>d</sup>* terms, which is exactly the number of components of a spinorial field in *D* ¼ 2*d* dimensions.

In the representation (Eq. (21)), the operator Γ*α*∇*α*, called Dirac operator, is then represented by

$$\Gamma^{d}\nabla\_{a} = \sum\_{a=1}^{d} \left(\Gamma\_{a}\nabla\_{a} + \Gamma\_{\bar{a}}\nabla\_{\bar{a}}\right) = \sum\_{a=1}^{d} \underbrace{\sigma\_{3}\otimes\dots\otimes\mathbb{I}\_{d}}\_{(a-1)\text{ times}}\underbrace{\mathbb{I}\otimes\dots\otimes\mathbb{I}}\_{(d-a)\text{ times}}\tag{25}$$

where

$$D\_a = \sigma\_1 \nabla\_a + \sigma\_2 \nabla\_{\bar{a}},\tag{26}$$

is the Dirac operator on <sup>ℝ</sup><sup>2</sup> with coordinates *xa, ya* f g. The spinorial basis introduced previously is very convenient, since the action of the Dirac matrices on the spinor fields can be easily computed. Indeed, using Eqs. (21), (23), and (24), we eventually arrive at the following equation

 $d$  denotes

<sup>2</sup> In *<sup>D</sup>* <sup>¼</sup> <sup>2</sup>*<sup>d</sup>* <sup>þ</sup> 1, besides the 2*<sup>d</sup>* Dirac matrices <sup>Γ</sup>*<sup>a</sup>* and <sup>Γ</sup>*a*~, we need to add one further matrix, which will be denoted by <sup>Γ</sup>*<sup>d</sup>*þ<sup>1</sup> given by <sup>Γ</sup>*<sup>d</sup>*þ<sup>1</sup> <sup>¼</sup> *<sup>σ</sup>*<sup>3</sup> <sup>⊗</sup> *<sup>σ</sup>*3… <sup>⊗</sup> *<sup>σ</sup>*<sup>3</sup> |fflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflffl} *<sup>d</sup>* times .

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes DOI: http://dx.doi.org/10.5772/intechopen.89179*

$$\begin{split} \Gamma\_{a}\Psi &= \sum\_{\{i\}} (s\_{1\xi 2\ldots s\_{d-1}}) \Psi^{\mu\_{1} \mu\_{2} \ldots \mu\_{d}} \mathfrak{F}^{\epsilon\_{1}} \otimes \mathfrak{F}^{\epsilon\_{1}} \otimes \ldots \otimes \mathfrak{F}^{\epsilon\_{d-1}} \otimes \mathfrak{F}^{-\epsilon\_{d}} \otimes \mathfrak{F}^{\epsilon\_{d+1}} \otimes \\ &\ldots \otimes \mathfrak{F}^{\epsilon\_{d}} = \sum\_{\{s\}} (s\_{1\xi 2\ldots s\_{d}}) \Psi^{\mu\_{1} \mu\_{2} \ldots \mu\_{d-1} (-s\_{d}) s\_{d+1} \ldots s\_{d}} \mathfrak{F}^{\epsilon\_{1}} \otimes \mathfrak{F}^{\epsilon\_{2}} \otimes \\ &\ldots \otimes \mathfrak{F}^{\epsilon\_{d-1}} \otimes \mathfrak{F}^{\epsilon\_{d}} \otimes \mathfrak{F}^{\epsilon\_{d+1}} \otimes \ldots \otimes \mathfrak{F}^{\epsilon\_{d}}, \end{split} \tag{27}$$

where from the first to the second line we have changed the index *sa* to �*sa*, which does not change the final result, since we are summing over all values of *sa*, which comprise the same list of the values of �*sa*. Moreover, we have used ð Þ *sa* <sup>2</sup> <sup>¼</sup> 1. Analogously, we have:

<sup>Γ</sup>*a*<sup>~</sup>**<sup>Ψ</sup>** <sup>¼</sup> <sup>X</sup> f g*s* ð Þ *<sup>s</sup>*1*s*2…*sa*�<sup>1</sup> ð Þ *isa* **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sd <sup>ξ</sup>s*<sup>1</sup> <sup>⊗</sup> *<sup>ξ</sup> <sup>s</sup>*<sup>2</sup> ⊗ … ⊗ *ξ sa*�<sup>1</sup> ⊗ *ξ* �*sa* ⊗ *ξ sa*þ<sup>1</sup> ⊗ … ⊗ *ξ sd* ¼ �*<sup>i</sup>* X f g*s* ð Þ *<sup>s</sup>*1*s*2…*sa sa* **<sup>Ψ</sup>***<sup>s</sup>*1*s*2…*sa*�1ð Þ �*sa sa*þ1…*sd <sup>ξ</sup><sup>s</sup>*<sup>1</sup> <sup>⊗</sup> *<sup>ξ</sup> <sup>s</sup>*<sup>2</sup> ⊗ … ⊗ *ξsa*�<sup>1</sup> ⊗ *ξ sa* ⊗ *ξ sa*þ<sup>1</sup> ⊗ … ⊗ *ξ sd :* (28)

All that was seen above are necessary tools to attack our initial problem of separating the general Eq. (14). In order to solve such an equation, we need to separate the degrees of freedom of the field, which can be quite challenging in general. Fortunately, the spacetime considered here is the direct product of twodimensional spaces of constant curvature, which is exactly the class of spaces studied in Ref. [39]. Indeed, in this latter paper, it is shown that the Dirac equation minimally coupled to an electromagnetic field is separable in such backgrounds. In particular, assuming that the components of the spinor field Eq. (24) can be decomposed in the form

$$
\Psi^{\mathfrak{e}\_1 \mathfrak{e}\_2 \ldots \mathfrak{e}\_d} = \Psi^{\mathfrak{e}\_1}\_1(t, \mathfrak{x}) \Psi^{\mathfrak{e}\_2}\_2(\Phi\_2, \theta\_2) \ldots \Psi^{\mathfrak{e}\_d}\_d(\Phi\_d, \theta\_d), \tag{29}
$$

where each index *sa* can take the values *sa* ¼ �1, the fields <sup>Ψ</sup>*<sup>s</sup>*<sup>1</sup> <sup>1</sup> ð Þ *t, x* satisfy the following differential equation (the reader is invited to demonstrate the equation below or consult more details in [39]):

$$\left[\partial\_{\bar{1}} + \frac{a\iota\_{1\bar{1}1}}{2} - iq\mathcal{A}\_{\bar{1}} - i\varepsilon\_{1}\left(\partial\_{1} + \frac{a\iota\_{\bar{1}1\bar{1}}}{2} - iq\mathcal{A}\_{1}\right)\right] \Psi\_{1}^{\varepsilon\_{1}} = \left(L - i\varepsilon\_{1}m\right)\Psi\_{1}^{-\varepsilon\_{1}}.\tag{30}$$

The separation constant *L* in the above equation depends on the angular modes. In particular, in the special case of vanishing magnetic charges *Qj* , it is determined by the eigenvalues *λ<sup>j</sup>* of the Dirac operator on unit sphere *S*<sup>2</sup> according to the following relation

$$L = \sqrt{\lambda\_2^2 + \lambda\_3^2 + \dots + \lambda\_d^2} \quad , \quad \lambda\_j = \pm 1, \pm 2, \pm 3, \dots \tag{31}$$

as demonstrated in Appendix A of Ref. [39]. In our frame of vectors, the only components of the spin connection that are potentially nonvanishing are

$$\begin{cases} \begin{aligned} \boldsymbol{\alpha}\_{1\bar{1}1} &= -\boldsymbol{\alpha}\_{1\bar{1}\bar{1}} = -\frac{1}{R\_1} \sinh\left(\boldsymbol{\alpha}/R\_1\right), \\\\ \boldsymbol{\alpha}\_{\bar{j}\bar{j}} &= -\boldsymbol{\alpha}\_{\bar{j}\bar{j}} = \frac{1}{R\_{\bar{j}}} \cot\theta\_{\bar{j}}, \end{aligned} \tag{32}$$

are the Hermitian Pauli matrices and denote the 2 � 2 identity matrix. Using this notation, a convenient representation of the Dirac matrices is the following:

where stands for the 2 � 2 identity matrix. Indeed, we can easily check that the Clifford algebra given in Eq. (15) is properly satisfied by the above matrices.<sup>2</sup> In this case, spinorial fields are represented by the column vectors on which these matrices act. We can introduce a basis of this representation by the direct products of spinors

*, ξ*

which, under the action of the Pauli matrices, satisfy concisely the relations

Indeed, in *<sup>D</sup>* <sup>¼</sup> <sup>2</sup>*<sup>d</sup>* dimensions, a general spinor field has 2*<sup>d</sup>* degrees of freedom

*<sup>s</sup>*<sup>1</sup> ⊗ *ξ*

where each of the indices *sa* can take the values "+1" and "�1." Since every *sa* can take just two values, it follows that the sum over f g*<sup>s</sup>* � f g *<sup>s</sup>*1*, s*2*,* …*, sd* comprises 2*<sup>d</sup>* terms, which is exactly the number of components of a spinorial field in *D* ¼ 2*d*

In the representation (Eq. (21)), the operator Γ*α*∇*α*, called Dirac operator, is then

*σ*<sup>3</sup> ⊗ … ⊗ |fflfflfflfflfflffl{zfflfflfflfflfflffl} ð Þ *a*�1 times

*d*

*a*¼1

is the Dirac operator on <sup>ℝ</sup><sup>2</sup> with coordinates *xa, ya* f g. The spinorial basis introduced previously is very convenient, since the action of the Dirac matrices on the spinor fields can be easily computed. Indeed, using Eqs. (21), (23), and (24), we

<sup>2</sup> In *<sup>D</sup>* <sup>¼</sup> <sup>2</sup>*<sup>d</sup>* <sup>þ</sup> 1, besides the 2*<sup>d</sup>* Dirac matrices <sup>Γ</sup>*<sup>a</sup>* and <sup>Γ</sup>*a*~, we need to add one further matrix, which will

.

*, <sup>σ</sup>*2*ξ<sup>s</sup>* <sup>¼</sup> *i s <sup>ξ</sup>*�*<sup>s</sup>*

**Ψ***<sup>s</sup>*1*s*2…*sd ξ*

� <sup>¼</sup> <sup>0</sup> 1 � �

*, <sup>σ</sup>*3*ξ<sup>s</sup>* <sup>¼</sup> *<sup>s</sup><sup>ξ</sup>*

*<sup>s</sup>*<sup>2</sup> ⊗ … ⊗ *ξ*

*s*

*Da* ⊗ ⊗ … ⊗ |fflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflffl} ð Þ *d*�*a* times

*Da* ¼ *σ*1∇*<sup>a</sup>* þ *σ*2∇*a*<sup>~</sup>*,* (26)

*σ*<sup>1</sup> ⊗ ⊗ … ⊗ |fflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflffl} ð Þ *d*�*a* times

*σ*<sup>2</sup> ⊗ ⊗ … ⊗ |fflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflffl} ð Þ *d*�*a* times

*,*

*,* (21)

*,* (22)

*:* (23)

*sd ,* (24)

*,* (25)

Γ*<sup>a</sup>* ¼ *σ*<sup>3</sup> ⊗ … ⊗ |fflfflfflfflfflffl{zfflfflfflfflfflffl} ð Þ *a*�1 times

Γ*a*<sup>~</sup> ¼ *σ*<sup>3</sup> ⊗ … ⊗ |fflfflfflfflfflffl{zfflfflfflfflfflffl} ð Þ *a*�1 times

*ξ* <sup>þ</sup> <sup>¼</sup> <sup>1</sup> 0 � �

*<sup>σ</sup>*1*ξ<sup>s</sup>* <sup>¼</sup> *<sup>ξ</sup>*�*<sup>s</sup>*

**<sup>Ψ</sup>** <sup>¼</sup> <sup>X</sup> f g*s*

<sup>ð</sup>Γ*a*∇*<sup>a</sup>* <sup>þ</sup> <sup>Γ</sup>*a*<sup>~</sup>∇*a*<sup>~</sup>Þ ¼ <sup>X</sup>

*ξ*

*<sup>s</sup>* given by

*Progress in Relativity*

and can be written as

dimensions.

represented by

where

**6**

<sup>Γ</sup>*<sup>α</sup>*∇*<sup>α</sup>* <sup>¼</sup> <sup>X</sup>

*d*

*a*¼1

eventually arrive at the following equation

be denoted by <sup>Γ</sup>*<sup>d</sup>*þ<sup>1</sup> given by <sup>Γ</sup>*<sup>d</sup>*þ<sup>1</sup> <sup>¼</sup> *<sup>σ</sup>*<sup>3</sup> <sup>⊗</sup> *<sup>σ</sup>*3… <sup>⊗</sup> *<sup>σ</sup>*<sup>3</sup> |fflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflffl} *<sup>d</sup>* times

and the nonzero components of the gauge field can be written as

$$\mathcal{A}\_1 = -i Q\_1 R\_1 \sinh\left(\mathbf{x}/R\_1\right), \quad \mathcal{A}\_{\dot{\jmath}} = Q\_{\dot{\jmath}} R\_{\dot{\jmath}} \cot \theta\_{\dot{\jmath}}.\tag{33}$$

Now, since the components of the metric are independents of the coordinate *t*, the vector *∂<sup>t</sup>* is a Killing vector for this metric. So, it is useful to assume the following time dependence for the field Ψ*s*<sup>1</sup> <sup>1</sup> ð Þ *t, x*

$$
\Psi\_1^{\epsilon\_1}(t,\varkappa) = e^{-i\alpha t} \psi^{\epsilon\_1}(\varkappa). \tag{34}
$$

which clearly is regular. So, we point out that for this potential both limits (Eqs. (39) and (40)) are finite, and thus there is no reason to demand for a regular

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

Thus, the problem of finding the QNMs is reduced to the searching of the corresponding spectrum of QNFs *ω* of Eq. (36). Most of the problems concerning the QNMs fall into Schrödinger-like equation with real potentials which vanish at both horizons [5], highlighting the fact that the solutions can be taken to be plane waves. However, clearly this is not the case. Although it is possible to make field redefinitions in order to make the potential real, we shall not do this here. For such procedure we refer the reader to [36]. Once an analytical form for the QNFs of Rosen-Morse type potential is not known, we must find an analytical exact solution of Eq. (36) and impose physically appropriate boundary conditions at the horizons,

In order to solve Eq. (36), let us make the following change of variable

with the parameters *α* and *β* being constants conveniently chosen as

*<sup>A</sup>* � *<sup>B</sup>* � *<sup>ω</sup>*<sup>2</sup> <sup>p</sup> *, <sup>β</sup>* ¼ � *<sup>R</sup>*<sup>1</sup>

the functions *H<sup>s</sup>*<sup>1</sup> must be solutions of the following differential equation

This new variable as well as the Ansatz that we have been using are really interesting because in terms of these, it is immediate to see that the functions *H<sup>s</sup>*<sup>1</sup>

satisfy a hypergeometric equation. Indeed, comparing with the standard

*dy*<sup>2</sup> <sup>þ</sup> ½ � *<sup>c</sup>* � ð Þ <sup>1</sup> <sup>þ</sup> *<sup>a</sup>* <sup>þ</sup> *<sup>b</sup> <sup>y</sup>*

<sup>2</sup> <sup>þ</sup> *<sup>α</sup>* <sup>þ</sup> *<sup>β</sup>* <sup>þ</sup>

<sup>2</sup> <sup>þ</sup> *<sup>α</sup>* <sup>þ</sup> *<sup>β</sup>* �

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

*dy*<sup>2</sup> <sup>þ</sup> ½ � <sup>2</sup>*<sup>α</sup>* <sup>þ</sup> <sup>1</sup> � ð Þ <sup>2</sup> <sup>þ</sup> <sup>2</sup>*<sup>α</sup>* <sup>þ</sup> <sup>2</sup>*<sup>β</sup> <sup>y</sup>*

*d*<sup>2</sup>*H<sup>s</sup>*<sup>1</sup>

we find that the constants *a*, *b*, and *c* are given by

8

>>>>>>>><

>>>>>>>>:

*<sup>a</sup>* <sup>¼</sup> <sup>1</sup>

*<sup>b</sup>* <sup>¼</sup> <sup>1</sup>

*c* ¼ 2*α* þ 1*:*

In particular, notice that *y* is defined on the domain *y*∈ð Þ 0*,* 1 with the boundaries *x* ! �∞ being given by *y* ¼ 0 and *y* ¼ 1. In addition to this change of inde-

tanh ð Þ *x=R*<sup>1</sup> *:* (41)

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *A* þ *B* � *ω*<sup>2</sup>

<sup>p</sup> *,* (43)

<sup>1</sup> <sup>þ</sup> ð Þ *<sup>α</sup>* <sup>þ</sup> *<sup>β</sup>* ð Þ <sup>1</sup> <sup>þ</sup> *<sup>α</sup>* <sup>þ</sup> *<sup>β</sup>* � �*H<sup>s</sup>*<sup>1</sup> *:*

*dy* � *abH<sup>s</sup>*<sup>1</sup> <sup>¼</sup> <sup>0</sup>*,* (45)

(44)

(46)

*<sup>ψ</sup><sup>s</sup>*<sup>1</sup> ð Þ¼ *<sup>x</sup> <sup>y</sup><sup>α</sup>*ð Þ <sup>1</sup> � *<sup>y</sup> <sup>β</sup> <sup>H</sup><sup>s</sup>*<sup>1</sup> ð Þ*<sup>y</sup> ,* (42)

2

*dy* � *CR*<sup>2</sup>

*dH<sup>s</sup>*<sup>1</sup>

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 <sup>4</sup> � *CR*<sup>2</sup> 1

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 <sup>4</sup> � *CR*<sup>2</sup> 1

*,*

*,*

r

r

*dH<sup>s</sup>*<sup>1</sup>

*<sup>y</sup>* <sup>¼</sup> <sup>1</sup> 2 þ 1 2

*x* ! �∞, which define the QNFs in a unique way.

pendent variable, if we now set the Ansatz

*<sup>α</sup>* <sup>¼</sup> *<sup>R</sup>*<sup>1</sup> 2

hypergeometric differential equation

*y*ð Þ 1 � *y*

*y*ð Þ 1 � *y*

**9**

*d*2 *H<sup>s</sup>*<sup>1</sup>

solution in this point.

Inserting this field along with the gauge field Eq. (33), and taking into account the first relation of the Eq. (32) into the Eq. (30), we end up with the following coupled system of differential equations:

$$
\left[\frac{d}{d\mathbf{x}} + \dot{\mathbf{s}}\_1 \boldsymbol{\omega} + \left(\dot{\mathbf{s}} \, q \, \mathbf{Q}\_1 \mathbf{R}\_1 - \frac{\mathbf{1}}{2\mathbf{R}\_1}\right) \tanh\left(\mathbf{x}/\mathbf{R}\_1\right)\right] \boldsymbol{\upmu^{\prime 1}} = \frac{\left(L - \dot{\mathbf{s}} \mathbf{1}\_1 \mathbf{m}\right)}{\cosh\left(\mathbf{x}/\mathbf{R}\_1\right)} \boldsymbol{\upmu^{-\prime 1}}.\tag{35}
$$

In order to solve these equations, we should first decouple the fields *ψ<sup>s</sup>*<sup>1</sup> and *ψ*�*s*<sup>1</sup> . Eliminating *ψ*�*s*<sup>1</sup> we obtain a second-order equation for *ψ<sup>s</sup>*<sup>1</sup> . Indeed, we can prove that the fields *ψ<sup>s</sup>*<sup>1</sup> satisfy the following second-order ordinary differential equation

$$
\left[\frac{d^2}{d\mathbf{x}^2} + \alpha^2 - V(\mathbf{x})\right] \boldsymbol{\psi}^{\boldsymbol{\epsilon}\_1} = \mathbf{0},\tag{36}
$$

which is a Schrödinger-like equation with *V* being a potential of the form

$$V(\mathbf{x}) = A + B \tanh\left(\mathbf{x}/R\_1\right) + \frac{C}{\cosh^2\left(\mathbf{x}/R\_1\right)},\tag{37}$$

where the parameters *A*, *B*, and *C* are given by

$$\begin{cases} A = \frac{1}{4R\_1^2} - q \, Q\_1(is\_1 + q \, Q\_1R\_1^2), \\\\ B = -\frac{a}{R\_1} \left( is\_1 + 2q \, Q\_1R\_1^2 \right), \\\\ C = m^2 + L^2 + \frac{1}{4R\_1^2} + q^2 \, Q\_1^2R\_1^2. \end{cases} \tag{38}$$

These are known as potentials of Rosen-Morse type, which are generalizations of the Pöschl-Teller potential [37, 38]. It is straightforward to see that this potential satisfies the following properties:

$$V \to \begin{cases} A + B & \text{at} \quad x \to +\infty, \\ A - B & \text{at} \quad x \to -\infty. \end{cases} \tag{39}$$

In many cases, the potential function *V* is regular at *r* ¼ 0ð Þ *x* ¼ 0 , in particular *V* can be equal to a constant different from zero. In fact, in our case, we find that

$$V \to A + C \quad \text{at} \quad x \to 0,\tag{40}$$

#### *Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes DOI: http://dx.doi.org/10.5772/intechopen.89179*

which clearly is regular. So, we point out that for this potential both limits (Eqs. (39) and (40)) are finite, and thus there is no reason to demand for a regular solution in this point.

Thus, the problem of finding the QNMs is reduced to the searching of the corresponding spectrum of QNFs *ω* of Eq. (36). Most of the problems concerning the QNMs fall into Schrödinger-like equation with real potentials which vanish at both horizons [5], highlighting the fact that the solutions can be taken to be plane waves. However, clearly this is not the case. Although it is possible to make field redefinitions in order to make the potential real, we shall not do this here. For such procedure we refer the reader to [36]. Once an analytical form for the QNFs of Rosen-Morse type potential is not known, we must find an analytical exact solution of Eq. (36) and impose physically appropriate boundary conditions at the horizons, *x* ! �∞, which define the QNFs in a unique way.

In order to solve Eq. (36), let us make the following change of variable

$$y = \frac{1}{2} + \frac{1}{2} \tanh\left(\mathfrak{x}/\mathfrak{R}\_1\right). \tag{41}$$

In particular, notice that *y* is defined on the domain *y*∈ð Þ 0*,* 1 with the boundaries *x* ! �∞ being given by *y* ¼ 0 and *y* ¼ 1. In addition to this change of independent variable, if we now set the Ansatz

$$
\psi^{s\_1}(\mathbf{x}) = \mathcal{Y}^a(\mathbf{1} - \mathbf{y})^\beta H^{s\_1}(\mathbf{y}),
\tag{42}
$$

with the parameters *α* and *β* being constants conveniently chosen as

$$a = \frac{R\_1}{2}\sqrt{A - B - a^2} \quad , \quad \beta = -\frac{R\_1}{2}\sqrt{A + B - a^2} \tag{43}$$

the functions *H<sup>s</sup>*<sup>1</sup> must be solutions of the following differential equation

$$y(1-y)\frac{d^2H^{\epsilon\_1}}{dy^2} + \left[2a+1-(2+2a+2\beta)y\right]\frac{dH^{\epsilon\_1}}{dy} - \left[\mathcal{CR}\_1^2 + (a+\beta)(1+a+\beta)\right]H^{\epsilon\_1}.\tag{44}$$

This new variable as well as the Ansatz that we have been using are really interesting because in terms of these, it is immediate to see that the functions *H<sup>s</sup>*<sup>1</sup> satisfy a hypergeometric equation. Indeed, comparing with the standard hypergeometric differential equation

$$y(\mathbf{1} - \mathbf{y}) \frac{d^2 H^{\mathbf{e}\_1}}{d \mathbf{y}^2} + \left[ c - (\mathbf{1} + a + b)y \right] \frac{d H^{\mathbf{e}\_1}}{d \mathbf{y}} - ab H^{\mathbf{e}\_1} = \mathbf{0},\tag{45}$$

we find that the constants *a*, *b*, and *c* are given by

$$\begin{cases} a = \frac{1}{2} + a + \beta + \sqrt{\frac{1}{4} - \text{CR}\_1^2} \\\\ b = \frac{1}{2} + a + \beta - \sqrt{\frac{1}{4} - \text{CR}\_1^2} \\\\ c = 2a + 1. \end{cases} \tag{46}$$

and the nonzero components of the gauge field can be written as

the vector *∂<sup>t</sup>* is a Killing vector for this metric. So, it is useful to assume the

<sup>1</sup> ð Þ¼ *t, x e*

Ψ*s*1

2*R*<sup>1</sup>

*d*2

where the parameters *A*, *B*, and *C* are given by

8

>>>>>>>>>><

>>>>>>>>>>:

*V* !

(

satisfies the following properties:

**8**

*<sup>A</sup>* <sup>¼</sup> <sup>1</sup> 4*R*<sup>2</sup> 1

*<sup>B</sup>* ¼ � *<sup>ω</sup> R*1

*<sup>C</sup>* <sup>¼</sup> *<sup>m</sup>*<sup>2</sup> <sup>þ</sup> *<sup>L</sup>*<sup>2</sup> <sup>þ</sup>

*dx*<sup>2</sup> <sup>þ</sup> *<sup>ω</sup>*<sup>2</sup> � *V x*ð Þ " #

*V x*ð Þ¼ *A* þ *B* tanh ð Þþ *x=R*<sup>1</sup>

which is a Schrödinger-like equation with *V* being a potential of the form

� *qQ*<sup>1</sup> *is*<sup>1</sup> <sup>þ</sup> *qQ*1*R*<sup>2</sup>

*is*<sup>1</sup> <sup>þ</sup> <sup>2</sup>*qQ*1*R*<sup>2</sup>

� �*,*

1 4*R*<sup>2</sup> 1

These are known as potentials of Rosen-Morse type, which are generalizations of the Pöschl-Teller potential [37, 38]. It is straightforward to see that this potential

> *A* þ *B* at *x* ! þ∞*, A* � *B* at *x* ! �∞*:*

In many cases, the potential function *V* is regular at *r* ¼ 0ð Þ *x* ¼ 0 , in particular *V* can be equal to a constant different from zero. In fact, in our case, we find that

� �*,*

1

<sup>þ</sup> *<sup>q</sup>*<sup>2</sup> *<sup>Q</sup>*<sup>2</sup> 1*R*2 1*:*

� �

� �

following time dependence for the field Ψ*s*<sup>1</sup>

coupled system of differential equations:

*dx* <sup>þ</sup> *is*1*<sup>ω</sup>* <sup>þ</sup> *isqQ*1*R*<sup>1</sup> � <sup>1</sup>

*d*

*Progress in Relativity*

equation

A<sup>1</sup> ¼ �*iQ*1*R*<sup>1</sup> sinh ð Þ *x=R*<sup>1</sup> *,* A*<sup>j</sup>* ¼ *QjRj* cot *θj:* (33)

*<sup>ψ</sup>s*<sup>1</sup> ð Þ *<sup>x</sup> :* (34)

*<sup>ψ</sup><sup>s</sup>*<sup>1</sup> <sup>¼</sup> <sup>0</sup>*,* (36)

*ψ*�*s*<sup>1</sup> *:* (35)

*,* (37)

(38)

(39)

*<sup>ψ</sup>s*<sup>1</sup> <sup>¼</sup> ð Þ *<sup>L</sup>* � *is*1*<sup>m</sup>* cosh ð Þ *x=R*<sup>1</sup>

*C* cosh <sup>2</sup>

1

*V* ! *A* þ *C* at *x* ! 0*,* (40)

ð Þ *x=R*<sup>1</sup>

Now, since the components of the metric are independents of the coordinate *t*,

<sup>1</sup> ð Þ *t, x*

�*iωt*

Inserting this field along with the gauge field Eq. (33), and taking into account the first relation of the Eq. (32) into the Eq. (30), we end up with the following

tanh ð Þ *x=R*<sup>1</sup>

In order to solve these equations, we should first decouple the fields *ψ<sup>s</sup>*<sup>1</sup> and *ψ*�*s*<sup>1</sup> . Eliminating *ψ*�*s*<sup>1</sup> we obtain a second-order equation for *ψ<sup>s</sup>*<sup>1</sup> . Indeed, we can prove that the fields *ψ<sup>s</sup>*<sup>1</sup> satisfy the following second-order ordinary differential

Such an equation admits two linearly independent solutions whose linear combination furnishes the following general solution:

$$H^{t\_1}(\mathbf{y}) = D\_2 F\_1(a, b, c; \mathbf{y}) + E \mathbf{y}^{(1-c)} \,\_2F\_1(\mathbf{1} + a + c, \mathbf{1} + b + c, \mathbf{2} - c; \mathbf{y}),\tag{47}$$

where <sup>2</sup>*F*<sup>1</sup> is the hypergeometric function and *D* and *E* are arbitrary integration constants. Given the hypergeometric solution for *Hs*<sup>1</sup> is known, one can immediately find the general solution for *ψs*<sup>1</sup> . Indeed, from Eqs. (42), (46), and (47), we conclude that the solution of Eq. (36), which is regular at the origin, can be written as

$$\begin{split} y^{\varepsilon\_1} &= (\mathbf{1} - \mathbf{y})^{\frac{1}{2}(a+b-c)} [Dy^{\frac{1}{2}(c-1)} \,\_2F\_1(a,b,c;\mathbf{y}) \\ &+ Ey^{-\frac{1}{2}(c-1)} \,\_2F\_1(\mathbf{1}+a-c, \mathbf{1}+b-c, \mathbf{2}-c; \mathbf{y})]. \end{split} \tag{48}$$

In order to fix the integration constants *D* and *E*, we need to apply the appropriate boundary conditions. Inverting the Eq. (41) we find that, near the boundaries *x* ! �∞, the relation between the coordinates *x* and *y* assumes the simpler form

$$\begin{array}{ll} \mathsf{y} \simeq e^{+2\mathbf{x}/R\_1} & \text{at} & \mathsf{x} \to -\infty, \\\\ \mathbf{1} - \mathbf{y} \simeq e^{-2\mathbf{x}/R\_1} & \text{at} & \mathbf{x} \to +\infty. \end{array} \tag{49}$$

*a* ¼ *iR*<sup>1</sup>

*b* ¼ �*iR*<sup>1</sup>

*<sup>c</sup>* <sup>¼</sup> <sup>1</sup>

Ψ<sup>þ</sup> <sup>1</sup> ð Þ *t, x* � �

, this means that Ψ<sup>þ</sup>

Ψ<sup>þ</sup> 1 � �

such a case, from Eq. (52), the field Ψ<sup>þ</sup>

we are led to the following constraint:

*ω* ¼ �

q

*e*�*iω<sup>t</sup>*

**11**

q

q

<sup>2</sup> <sup>þ</sup> *is*<sup>1</sup> *qQ*1*R*<sup>2</sup>

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *<sup>μ</sup>*<sup>2</sup> <sup>þ</sup> *<sup>q</sup>*2*Q*<sup>2</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

1*R*2 <sup>1</sup> <sup>þ</sup> *<sup>L</sup>*<sup>2</sup>

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *<sup>μ</sup>*<sup>2</sup> <sup>þ</sup> *<sup>q</sup>*2*Q*<sup>2</sup>

<sup>1</sup> � *ωR*<sup>1</sup> � �*:*

In particular, the following relations hold

at the following behavior of the solution at *x* ! ∞:

*<sup>x</sup>*!þ<sup>∞</sup> <sup>≃</sup> *<sup>E</sup>* <sup>Γ</sup>ð Þ *<sup>c</sup>* � *<sup>a</sup>* � *<sup>b</sup>* <sup>Γ</sup>ð Þ <sup>2</sup> � *<sup>c</sup>*

*<sup>x</sup>*!�<sup>∞</sup> <sup>¼</sup> *D e*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup> <sup>e</sup>*

1*R*2 <sup>1</sup> <sup>þ</sup> *<sup>L</sup>*<sup>2</sup>

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

þ ð Þ 1 þ *s*<sup>1</sup>

þ ð Þ 1 þ *s*<sup>1</sup>

ð Þ *<sup>c</sup>* � <sup>1</sup> *<sup>=</sup>R*<sup>1</sup> ¼ �*is*1*<sup>ω</sup>* <sup>þ</sup> *is*1*qQ*1*R*<sup>1</sup> � <sup>1</sup>

ð Þ *a* þ *b* � *c =R*<sup>1</sup> ¼ �*iω* � *iqQ*1*R*<sup>1</sup> þ *s*<sup>1</sup>

*iqQ*1*R*1� <sup>1</sup> 2*R*1 � �*<sup>x</sup>*

lower values of *x*. Then, since the time dependence of the field Ψ<sup>þ</sup>

Γð Þ 1 � *a* Γð Þ 1 � *b* � �*<sup>e</sup>*

Using the Eq. (53), we find that these constraints translate to

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *<sup>m</sup>*<sup>2</sup> <sup>þ</sup> *<sup>q</sup>*<sup>2</sup>*Q*<sup>2</sup>

1*R*2 <sup>1</sup> <sup>þ</sup> *<sup>L</sup>*<sup>2</sup>

<sup>þ</sup>*<sup>E</sup>* <sup>Γ</sup>ð Þ *<sup>a</sup>* <sup>þ</sup> *<sup>b</sup>* � *<sup>c</sup>* <sup>Γ</sup>ð Þ <sup>2</sup> � *<sup>c</sup>* Γð Þ *a* � *c* þ 1 Γð Þ *b* � *c* þ 1 � �*<sup>e</sup>*

Now we are ready to impose the boundary conditions. Obviously, without loss of generality, we can consider that the spin *s*<sup>1</sup> is already chosen and fixed at *s*<sup>1</sup> ¼ þ or *s*<sup>1</sup> ¼ � since the QNFs should not depend on the choice of *s*<sup>1</sup> ¼ �. Let us impose, for instance, the boundary conditions for the component *s*<sup>1</sup> ¼ þ of the spinorial field. In this case, using the identity Eq. (54) along with the Eq. (34), we eventually arrive

Now, **Figure 1** tells us that the field is assumed to move toward higher values of *x* at the boundary *x* ! �∞, while at the boundary *x* ! �∞ it should move toward

*<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup>* at *<sup>x</sup>* ! þ∞. Thus, from Eq. (55), we conclude that we must set *<sup>D</sup>* <sup>¼</sup> 0. In

<sup>1</sup> becomes

Finally, to satisfy the QNM boundary condition near the boundary at *x* ! ∞, we must eliminate the term *e*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* of the above equation. Since *E* cannot be zero (as otherwise the field would vanish identically), we need the combination of the gamma functions to be zero. Now, once the gamma function has no zeros, the way to achieve this is to let the gamma functions in the denominator diverge, Γð Þ¼ 1 � *a* ∞ or Γð Þ¼ 1 � *b* ∞. Since the gamma functions diverge only at nonpositive integers,

1 � *a* ¼ �*n* or 1 � *b* ¼ �*n,* where *n*∈ f g 0*,* 1*,* 2*,* … *:* (58)

þ *i R*1 *n* þ 1 2

�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup> e*

1 <sup>4</sup> � *<sup>i</sup>ωR*<sup>1</sup> 2

> 1 <sup>4</sup> � *<sup>i</sup>ωR*<sup>1</sup> 2

� �

� �

� *i*ð Þ 1 � *s*<sup>1</sup>

2*R*<sup>1</sup>

<sup>þ</sup> *E e*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup> <sup>e</sup>*

<sup>1</sup> should behave as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* at *<sup>x</sup>* ! �∞, while it should go as

*iqQ*1*R*1� <sup>1</sup> 2*R*1 � �*<sup>x</sup>*

�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup> <sup>e</sup>* �*iqQ*1*R*1<sup>þ</sup> <sup>1</sup>

2*R*1 � �*<sup>x</sup>*

� �*,* (59)

*:*

1 2*R*<sup>1</sup>

� *iqQ*1*R*1� <sup>1</sup>

2*R*1 � �*<sup>x</sup>*

*:* (56)

(57)

<sup>1</sup> is of the type

� *i*ð Þ 1 � *s*<sup>1</sup>

*qQ*1*R*<sup>2</sup> 1 <sup>2</sup> *,*

> *qQ*1*R*<sup>2</sup> 1 <sup>2</sup> *,*

*,* (54)

*:* (55)

(53)

Thus, taking into account the latter relation and using the fact that at *y* ¼ 0ð Þ *x* ! �∞ the hypergeometric function <sup>2</sup>*F*1ð Þ¼ *a, b, c;*0 1, one eventually obtains that near the boundary *<sup>x</sup>* ! �<sup>∞</sup> the field *<sup>ψ</sup><sup>s</sup>*<sup>1</sup> behaves as

$$\left. \psi^{\epsilon\_1} \right|\_{\mathbf{x} \to -\infty} \simeq D\epsilon^{(c-1)\mathbf{x}/\mathcal{R}\_1} + E\epsilon^{-(c-1)\mathbf{x}/\mathcal{R}\_1}.\tag{50}$$

On the other hand, in order to apply the boundary conditions at *y* ¼ 1 ð Þ *x* ! ∞ , it is useful to write the hypergeometric functions as functions of 1ð Þ � *y* , so that they become united at the boundary. This can be done by rewriting the hypergeometric functions appearing in Eq. (48) by means of the following identity [45]:

$$\begin{aligned} 2F(a,b,c;y) &= \frac{\Gamma(c)\Gamma(c-a-b)}{\Gamma(c-a)\Gamma(c-b)} F(a,b,a+b-c+1;1-y) \\ &+ \frac{\Gamma(c)\Gamma(a+b-c)}{\Gamma(a)\Gamma(b)} (\mathbf{1}-y)^{(c-a-b)} \,\_2F(c-a,c-b,c-a-b+1;1-y), \end{aligned} \tag{51}$$

where Γ stands for the gamma function. Doing so, and using Eq. (49), we eventually arrive at the following behavior of the solution at *x* ! þ∞:

$$\begin{split} \left. \Psi^{\epsilon \mathbf{1}} \right|\_{\mathbf{x} \to +\infty} & \approx \left[ D \frac{\Gamma(c-a-b)\Gamma(c)}{\Gamma(c-a)\Gamma(c-b)} + E \frac{\Gamma(c-a-b)\Gamma(2-c)}{\Gamma(1-a)\Gamma(1-b)} \right] e^{-(a+b-c)\mathbf{x}/R\_{1}} \\ & + \left[ D \frac{\Gamma(a+b-c)\Gamma(c)}{\Gamma(a)\Gamma(b)} + E \frac{\Gamma(a+b-c)\Gamma(2-c)}{\Gamma(a-c+1)\Gamma(b-c+1)} \right] e^{(a+b-c)\mathbf{x}/R\_{1}}. \end{split} \tag{52}$$

Now, from parameters Eqs. (38) and (43), we find that the constants appearing in the hypergeometric equation can be written as

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes DOI: http://dx.doi.org/10.5772/intechopen.89179*

$$\begin{split} a &= iR\_1\sqrt{\mu^2 + q^2 Q\_1^2 R\_1^2 + L^2} + (1 + s\_1)\left(\frac{1}{4} - i\nu \frac{R\_1}{2}\right) - i(1 - s\_1)\frac{qQ\_1 R\_1^2}{2}, \\ b &= -iR\_1\sqrt{\mu^2 + q^2 Q\_1^2 R\_1^2 + L^2} + (1 + s\_1)\left(\frac{1}{4} - i\nu \frac{R\_1}{2}\right) - i(1 - s\_1)\frac{qQ\_1 R\_1^2}{2}, \\ c &= \frac{1}{2} + i\nu\_1 (qQ\_1 R\_1^2 - \nu R\_1). \end{split} \tag{53}$$

In particular, the following relations hold

Such an equation admits two linearly independent solutions whose linear com-

where <sup>2</sup>*F*<sup>1</sup> is the hypergeometric function and *D* and *E* are arbitrary integration constants. Given the hypergeometric solution for *Hs*<sup>1</sup> is known, one can immediately find the general solution for *ψs*<sup>1</sup> . Indeed, from Eqs. (42), (46), and (47), we conclude

<sup>2</sup>ð Þ *c*�1

In order to fix the integration constants *D* and *E*, we need to apply the appropriate boundary conditions. Inverting the Eq. (41) we find that, near the boundaries *x* ! �∞, the relation between the coordinates *x* and *y* assumes the simpler form

<sup>þ</sup>2*x=R*<sup>1</sup> at *<sup>x</sup>* ! �∞*,*

�2*x=R*<sup>1</sup> at *<sup>x</sup>* ! þ∞*:*

Thus, taking into account the latter relation and using the fact that at *y* ¼ 0ð Þ *x* ! �∞ the hypergeometric function <sup>2</sup>*F*1ð Þ¼ *a, b, c;*0 1, one eventually obtains

On the other hand, in order to apply the boundary conditions at *y* ¼ 1 ð Þ *x* ! ∞ ,

*F a, b, a* ð Þ þ *b* � *c* þ 1*;*1 � *y*

Γð Þ *c* � *a* � *b* Γð Þ 2 � *c* Γð Þ 1 � *a* Γð Þ 1 � *b*

Γð Þ *a* � *c* þ 1 Γð Þ *b* � *c* þ 1

it is useful to write the hypergeometric functions as functions of 1ð Þ � *y* , so that they become united at the boundary. This can be done by rewriting the hypergeometric functions appearing in Eq. (48) by means of the following

<sup>2</sup>*F*1ð Þ *a; b;c*; *y*

*<sup>x</sup>*!�<sup>∞</sup> <sup>≃</sup> *D e*ð Þ *<sup>c</sup>*�<sup>1</sup> *<sup>x</sup>=R*<sup>1</sup> <sup>þ</sup> *E e*�ð Þ *<sup>c</sup>*�<sup>1</sup> *<sup>x</sup>=R*<sup>1</sup> *:* (50)

<sup>2</sup>*F c*ð Þ � *a, c* � *b, c* � *a* � *b* þ 1*;*1 � *y ,*

*e*

*e*

�ð Þ *a*þ*b*�*c x=R*<sup>1</sup>

ð Þ *<sup>a</sup>*þ*b*�*<sup>c</sup> <sup>x</sup>=R*<sup>1</sup> *:*

<sup>2</sup>*F*1ð1 þ *a* � *c;* 1 þ *b* � *c;* 2 � *c*; *y*Þ�*:*

that the solution of Eq. (36), which is regular at the origin, can be written as

<sup>2</sup>ð Þ *<sup>a</sup>*þ*b*�*<sup>c</sup>* ½*Dy*<sup>1</sup>

<sup>2</sup>*F*1ð Þ 1 þ *a* þ *c,* 1 þ *b* þ *c,* 2 � *c;y ,* (47)

(48)

(49)

(51)

(52)

bination furnishes the following general solution:

*<sup>ψ</sup>s*<sup>1</sup> <sup>¼</sup> ð Þ <sup>1</sup> � *<sup>y</sup>* <sup>1</sup>

<sup>þ</sup> *Ey*�<sup>1</sup>

<sup>2</sup>ð Þ *c*�1

*y*≃*e*

1 � *y* ≃*e*

that near the boundary *<sup>x</sup>* ! �<sup>∞</sup> the field *<sup>ψ</sup><sup>s</sup>*<sup>1</sup> behaves as

*ψs*1 j

Γð Þ *c* � *a* Γð Þ *c* � *b* <sup>2</sup>

Γð Þ*c* Γð Þ *a* þ *b* � *c*

Γð Þ *c* � *a* � *b* Γð Þ*c* <sup>Γ</sup>ð Þ *<sup>c</sup>* � *<sup>a</sup>* <sup>Γ</sup>ð Þ *<sup>c</sup>* � *<sup>b</sup>* <sup>þ</sup> *<sup>E</sup>*

<sup>þ</sup> *<sup>D</sup>* <sup>Γ</sup>ð Þ *<sup>a</sup>* <sup>þ</sup> *<sup>b</sup>* � *<sup>c</sup>* <sup>Γ</sup>ð Þ*<sup>c</sup>*

in the hypergeometric equation can be written as

<sup>Γ</sup>ð Þ *<sup>a</sup>* <sup>Γ</sup>ð Þ *<sup>b</sup>* ð Þ <sup>1</sup> � *<sup>y</sup>* ð Þ *<sup>c</sup>*�*a*�*<sup>b</sup>*

eventually arrive at the following behavior of the solution at *x* ! þ∞:

where Γ stands for the gamma function. Doing so, and using Eq. (49), we

<sup>Γ</sup>ð Þ *<sup>a</sup>* <sup>Γ</sup>ð Þ *<sup>b</sup>* <sup>þ</sup> *<sup>E</sup>* <sup>Γ</sup>ð Þ *<sup>a</sup>* <sup>þ</sup> *<sup>b</sup>* � *<sup>c</sup>* <sup>Γ</sup>ð Þ <sup>2</sup> � *<sup>c</sup>*

Now, from parameters Eqs. (38) and (43), we find that the constants appearing

identity [45]:

*ψs*1 j

**10**

<sup>2</sup>*F a, b, c;y* ð Þ¼ <sup>Γ</sup>ð Þ*<sup>c</sup>* <sup>Γ</sup>ð Þ *<sup>c</sup>* � *<sup>a</sup>* � *<sup>b</sup>*

þ

*<sup>x</sup>*!þ<sup>∞</sup> ≃ *D*

*<sup>H</sup>s*<sup>1</sup> ð Þ¼ *<sup>y</sup> <sup>D</sup>*2*F*1ð Þþ *a, b, c;y Ey*ð Þ <sup>1</sup>�*<sup>c</sup>*

*Progress in Relativity*

$$(c - 1)/R\_1 = -is\_1o + is\_1qQ\_1R\_1 - \frac{1}{2R\_1},\tag{54}$$

$$(a+b-c)/R\_1 = -i\nu - iqQ\_1R\_1 + s\_1\frac{1}{2R\_1}.\tag{55}$$

Now we are ready to impose the boundary conditions. Obviously, without loss of generality, we can consider that the spin *s*<sup>1</sup> is already chosen and fixed at *s*<sup>1</sup> ¼ þ or *s*<sup>1</sup> ¼ � since the QNFs should not depend on the choice of *s*<sup>1</sup> ¼ �. Let us impose, for instance, the boundary conditions for the component *s*<sup>1</sup> ¼ þ of the spinorial field. In this case, using the identity Eq. (54) along with the Eq. (34), we eventually arrive at the following behavior of the solution at *x* ! ∞:

$$\left. \Psi\_{1}^{+}(t, \mathbf{x}) \right|\_{\mathbf{x} \to -\infty} = D e^{-i\alpha(t+\mathbf{x})} e^{\left(iqQ\_{1}R\_{1} - \frac{1}{2R\_{1}}\right)\mathbf{x}} + E e^{-i\alpha(t-\mathbf{x})} e^{-\left(iqQ\_{1}R\_{1} - \frac{1}{2R\_{1}}\right)\mathbf{x}}.\tag{56}$$

Now, **Figure 1** tells us that the field is assumed to move toward higher values of *x* at the boundary *x* ! �∞, while at the boundary *x* ! �∞ it should move toward lower values of *x*. Then, since the time dependence of the field Ψ<sup>þ</sup> <sup>1</sup> is of the type *e*�*iω<sup>t</sup>* , this means that Ψ<sup>þ</sup> <sup>1</sup> should behave as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* at *<sup>x</sup>* ! �∞, while it should go as *<sup>e</sup>*�*iω*ð Þ *<sup>t</sup>*þ*<sup>x</sup>* at *<sup>x</sup>* ! þ∞. Thus, from Eq. (55), we conclude that we must set *<sup>D</sup>* <sup>¼</sup> 0. In such a case, from Eq. (52), the field Ψ<sup>þ</sup> <sup>1</sup> becomes

$$\begin{split} \left. \Psi\_{\mathbf{1}}^{+} \right|\_{\mathbf{x} \to +\infty} & \simeq E \left[ \frac{\Gamma(c - a - b) \Gamma(2 - c)}{\Gamma(1 - a) \Gamma(1 - b)} \right] e^{-i a (t - \mathbf{x})} e^{\left( i q Q\_{1} R\_{1} - \frac{1}{2 \mathbf{f}\_{1}} \right) \mathbf{x}} \\ & + E \left[ \frac{\Gamma(a + b - c) \Gamma(2 - c)}{\Gamma(a - c + 1) \Gamma(b - c + 1)} \right] e^{-i a (t + \mathbf{x})} e^{\left( -i q Q\_{1} R\_{1} + \frac{1}{2 \mathbf{f}\_{1}} \right) \mathbf{x}} . \end{split} \tag{57}$$

Finally, to satisfy the QNM boundary condition near the boundary at *x* ! ∞, we must eliminate the term *e*�*iω*ð Þ *<sup>t</sup>*�*<sup>x</sup>* of the above equation. Since *E* cannot be zero (as otherwise the field would vanish identically), we need the combination of the gamma functions to be zero. Now, once the gamma function has no zeros, the way to achieve this is to let the gamma functions in the denominator diverge, Γð Þ¼ 1 � *a* ∞ or Γð Þ¼ 1 � *b* ∞. Since the gamma functions diverge only at nonpositive integers, we are led to the following constraint:

$$\mathbf{1} - \mathfrak{a} = -n \quad \text{or} \quad \mathbf{1} - \mathfrak{b} = -n, \quad \text{where} \quad n \in \{0, 1, 2, \ldots\}. \tag{58}$$

Using the Eq. (53), we find that these constraints translate to

$$
\omega = \pm \sqrt{m^2 + q^2 Q\_1^2 R\_1^2 + L^2} + \frac{i}{R\_1} \left( n + \frac{1}{2} \right), \tag{59}
$$

which are the QNFs of the Dirac field propagating in *D*-dimensional generalized Nariai spacetimes. The real part of a QNF is associated with the oscillation frequency, while the imaginary part is related to its decay rate. At this point, it is worth recalling that *L* is a separation constant of the Dirac equation that is related to the angular mode of the field.

field. Although bosonic fields like scalar, electromagnetic, and gravitational fields can exhibit superradiant behavior in four-dimensional Kerr spacetime [46],

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes*

*DOI: http://dx.doi.org/10.5772/intechopen.89179*

curiously, this is not the case for the Dirac field [36]. Thus, it would be interesting to investigate whether an analogous thing happens in the background considered

here [47].

**Author details**

Joás Venâncio<sup>1</sup>

**13**

\* and Carlos Batista<sup>2</sup>

provided the original work is properly cited.

1 Physics Program, Universidade Federal de Pernambuco, Brazil

\*Address all correspondence to: joasvenancio@df.ufpe.br

2 Department of Physics, Universidade Federal de Pernambuco, Brazil

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

Likewise, imposing the boundary condition to the component *s*<sup>1</sup> ¼ � of the spinorial field, we find that we must set *E* ¼ 0 at Eq. (50) and then *c* � *a* ¼ �*n* or *c* � *b* ¼ �*n*, with *n* being a nonnegative integer. This, in its turn, leads to the same spectrum obtained for the component *s*<sup>1</sup> ¼ þ as expected, namely, Eq. (59).

#### **4. Conclusions**

In this chapter we have investigated the perturbations on a spinorial field propagating in a generalized version of the charged Nariai spacetime. Besides the separability of the degrees of freedom of these perturbations, one interesting feature of this background is that the perturbations can be analytically integrated. They all obey a Schrödinger-like equation with an integrable potential that is contained in the Rosen-Morse class of integrable potentials. Such an equation admits two linearly independent solutions given in terms of standard hypergeometric functions. This is a valuable property, since even the perturbation potential associated to the humble Schwarzschild background is nonintegrable, despite the fact that it is separable. We have also investigated the QNMs associated to this spinorial field. Analyzing the Eq. (59), namely,

$$
\rho\_{\rm D} = \pm \sqrt{m^2 + q^2 Q\_1^2 R\_1^2 + L^2} + \frac{i}{R\_1} \left( n + \frac{1}{2} \right), \tag{60}
$$

it is interesting to note that the imaginary parts of the QNFs, which represent the decay rates, do not depend on any details of the perturbation; rather, they only depend on the charges of the gravitational background through the dependence on *R*1. On the other hand, the real parts of the QNFs depend on the mass of the field and on the angular mode of the perturbations. Another fact worth pointing out is that the fermionic field always has a real part in its QNFs spectrum, meaning that it always oscillates. This is not reasonable. Indeed, for Klein-Gordon and Maxwell perturbations in the *D*-dimensional Nariai spacetime, their QNFs are equal to [39].

$$\begin{split} \rho\_{\rm KG} &= \pm \sqrt{m^2 + \sum\_{j=2}^{d} \frac{\ell\_j \left(\ell\_j + 1\right)}{R\_j^2} - \frac{1}{4R\_1^2}} - \frac{i}{R\_1} \left(n + \frac{1}{2}\right), \\\rho\_{\rm MF} &= \pm \sqrt{\sum\_{j=2}^{d} \frac{\ell\_j \left(\ell\_j + 1\right)}{R\_j^2} - \frac{1}{4R\_1^2}} - \frac{i}{R\_1} \left(n + \frac{1}{2}\right), \end{split} \tag{61}$$

where ℓ*<sup>j</sup>* and *mj* are integers, ∣*mj*∣ ≤ ℓ*j*, and ℓ≥0. Due to the negative factor �1*<sup>=</sup>* <sup>4</sup>*R*<sup>2</sup> 1 � � inside the square root appearing in the bosonic spectrum, it follows that for small enough *R*1, along with small enough mass and angular momentum, the argument of the square root can be negative, so that this term becomes imaginary.

To finish, we believe that a good exercise is to calculate the QNFs of the gravitational field in *D*-dimensional generalized charged Nariai spacetime. Research on the latter problem is still ongoing and, due to the great number of degrees of freedom in the gravitational field, shall be considered in a future work. The next interesting step is the investigation of superradiance phenomena for the spin 1*=*2

*Quasinormal Modes of Dirac Field in Generalized Nariai Spacetimes DOI: http://dx.doi.org/10.5772/intechopen.89179*

field. Although bosonic fields like scalar, electromagnetic, and gravitational fields can exhibit superradiant behavior in four-dimensional Kerr spacetime [46], curiously, this is not the case for the Dirac field [36]. Thus, it would be interesting to investigate whether an analogous thing happens in the background considered here [47].
