Spinor walls in five-dimensional warped spacetime

We study domain wall solutions of a real spinor field coupling with gravitation in five dimensions. We find that the nonlinear spinor field supports a class of soliton configurations which could be viewed as a wall embedded in five dimensions. We begin with an illuminating solution of the spinor field in the absence of gravitation. In a further investigation, we exhibit three sets of solutions of the spinor field with nonconstant curvature bulk spacetimes and three sets of solutions corresponding to three constant curvature bulk spacetimes. We demonstrate that some of these solutions in specific conditions have the energy density distributions of domain walls for the spinor field, where the scalar curvature is regular everywhere. Therefore, the configurations of these walls can be interpreted as spinor walls which are interesting spinor field realizations of domain walls. In order to investigate the stability of these spinor configurations, the linear perturbations are considered. The localization of the zero mode of tensor perturbation is also discussed.


Introduction
Many nonlinear physical systems support soliton solutions, which represent spatially localized field configurations.Kinks and domain walls are a subclass of co-dimension 1 topological solitons.Kinks are stable configurations interpolating between two of more degenerate vacua in (1 + 1)-dimensional spacetime, while domain walls are the higherdimensional analogues of kinks.They have the following properties: (i) the corresponding solutions are static and depend only on one spatial coordinate; (ii) they are topologically stable and indestructible -once a kink or domain wall is created it cannot spontaneously disappear.They are invaluable for nonperturbative aspects of field theories.A five-dimensional domain wall model in the absence of gravitation was presented in Ref. [1], where fermions can be confined on the wall by the Yukawa coupling between fermions and the background scalar field generating the domain wall.Kinks and domain walls in the presence of gravitation are more remarkable for extra dimension scenarios, especially for brane world frameworks.In the celebrated Randall-Sundrum (RS) II model, gravitation can be localized in the vicinity of a brane, a subspace embedded in five-dimensional anti-de Sitter (AdS) spacetime [2].It was shown that the zero mode of the graviton is localized on the delta-function-like brane and is responsible for the Newtonian potential, while the massive Kaluza-Klein (KK) modes of the graviton make corrections to the Newtonian potential.However, in the RS II model, the bulk scalar curvature is inevitably divergent at the location of the brane, and the extrinsic curvature satisfies the Israel junction condition.Thick branes with smooth warped extra dimensions have no such curvature singularity and could be viewed as domain walls in the presence of gravitation.For previous works on thick branes see Refs.[3][4][5][6][7][8][9][10][11] and recent reviews on thick branes see [12][13][14].
There are many thick brane models constructed by one or more scalar fields, which have achieved some features, such as the localization of matter fields, and the Newtonian potential recovered on the brane.The key point of using scalar fields is that a soliton configuration could be generated.Similarly, there are soliton solutions of spinor fields [15][16][17][18][19].A natural question is whether spinor fields could be used to construct thick branes.It is known that the equations of motion of ordinary spinor fields without selfinteractions are linear and can not produce a soliton configuration.Therefore, nonlinear spinor fields dominate the study of such kinds of solitons.Spinor fields including nonlinear ones are widely used in many fields mainly thanks to their specific behaviour in the presence of the gravitational field.Specifically, spinor fields were explored for particlelike solutions [20], and nonlinear spinor fields were explored to explain dark energy in the quintom scenario [21,22].In addition, as one of the most promising theories for unifying gravitation and other fundamental interactions, string theory at a low energy limit includes the contribution of torsion [23].Nonlinear spinor fields might represent the existence of spacetime torsion (a geometrical entity as spacetime curvature), which is suggested by Einstein-Cartan theory [24].An accepted fact is that torsion does not exist in our current universe.Even though it does, it will be extremely small.To interpret a torsion-free brane, the authors of Ref. [25] proposed an approach to create the illusion of a torsion-free universe from a bulk with torsion.Therefore, nonlinear spinor fields in extra dimension scenarios might provide some ways to shed light on where the torsion of the world is.
Nonlinear spinor fields were early investigated in Refs.[26,27], for many further explorations see the review [28].Heisenberg studied the possibility of nonlinear spinor fields as a theory of elementary particles [29,30].Gürsey proposed a spinor field equation which is similar to Heisenberg's nonlinear generalization of the Dirac equation [31].This equation is the first nonlinear conformal invariant wave equation which is viewed ever as a possible basis for a unitary description of elementary particles.An example of a selfinteracting spinor field, governed by a nonlinear field equation, is the Gürsey equation [31] γ µ ∂ µ ψ + σ( ψψ) which can be derived from the Lagrangian density Here, σ is the coupling strength to measure the Gürsey self-interaction.In order to reserve the conformal invariance of the Lagrangian density, a nonpolynomial form of nonlinear spinor fields was used to write this Lagrangian density.For our interest, the conformal invariance of this Lagrangian density is not required.Nevertheless, we allow the existence of nonlinear spinor fields.The solutions to Gürsey's equation were studied [32].In this paper, we will investigate soliton solutions in the Gürsey-like equation.Elementary fermions mathematically described by spinors are an essential component of our world.In previous studies on brane worlds, spinor fields were introduced for con-structing thick branes in Refs.[33,34], but there is no anticipative warp of the bulk, namely gravitation can not be localized on the brane.We attempt to explore the implementation of gravitational domain walls of spinor fields by extracting insights from domain walls generated by scalar fields.In our exploration, nonlinear spinor fields are considered and the equations of motion are Gürsey-like equations.
This paper is organized as follows.In section 2, we give a brief review of spinor fields in five-dimensional spacetime.In section 3, we introduce a toy model for a spinor field in five-dimensional flat spacetime.In section 4, analytic solutions of spinor fields are studied both in nonconstant curvature and constant curvature spacetimes.In section 5, we investigate the stability of the spinor walls.In section 6, the zero modes of tensor perturbations are discussed.Finally, discussion and conclusions are given in section 7.

Spinor fields in five-dimensional spacetime
We firstly review the theory of spinor fields in five-dimensional spacetime.There are many different conventions in use for spinor fields.Our conventions about spinor fields are closest to those in Ref. [35].Introduction to spinor fields refers to Refs.[36,37].
In general, gamma matrices are introduced to construct generators of Lorentz group in spinor representation and can be chosen in any representation.Considering the chiral nature of spinor fields, the gamma matrices which we adopt in five-dimensional flat spacetime are where the indices with ¯are the ones of flat spacetime coordinates, I is a 2 × 2 unity matrix, and σ ī are Pauli matrices with ī = 1, 2, 3.These gamma matrices satisfy the Clifford algebra {γ M , γ . The action of a spinor field in flat spacetime reads where ψ = iψ † γ 0 is the conjugation of the spinor ψ which generally has four complex components in five-dimensional spacetime.This potential term may include self-interactions of the spinor field.The equations of motion for the spinor field are 3) The energy-momentum tensor (on-shell), i.e., the Belinfante tensor, for the spinor field in the action (2.2) is arrived at [38] T It is well known that there are no spinor representations of the general linear group.A local tangent space should introduce local inertial frames, vielbeins, which are necessary to introduce spinor fields in curved spacetime.The relation between the spacetime metric and the tangent space metric is where the indices with and without ¯are the local Lorentz frame indices and the general coordinate indices, respectively; e M M = ∂x M ∂x M are the vielbeins.They also satisfy the orthogonality conditions The gamma matrices in curved spacetime are related to the ones in flat spacetime through the vielbeins as They satisfy the Clifford algebra {γ M , γ N } = 2g M N .The covariant derivative of a spinor field is given by Here ψ = iψ † γ 0 is the conjugation of the spinor ψ with four complex components in fivedimensional curved spacetime; ω

Ā B
M can be explicitly written in terms of the vielbeins as [39] ω and S Ā B , the generator of spinor representation of Lorentz group, is where Σ Ā B is the Pauli operator.
The action for a spinor field minimally coupling with gravitation reads where e = √ −g is the determinant of the vielbeins e M M .A potential term with selfinteractions of the spinor field could be included in the action (2.11).The variation of the action (2.11) with respect to the spinor field yields While varying the action with respect to the vielbeins e M M gives the Belinfante tensor for the spinor field The energy-momentum tensor (2.13) is not only conserved but also symmetric.Therefore, it can act as a source of the gravitational field.Generally, solutions of a spinor field will be rather complicated in the presence of gravitation.Therefore, we will not exactly restrict the potential term V ( ψ, ψ) for the spinor field ψ.There are many possibilities to construct Lorentz invariant terms responsible for the potential of the spinor field.The simplest Lorentz invariant term is ψψ.Without loss of generality, we adopt the following potential term for the spinor field in the actions (2.2) and (2.11) which is a function of the Lorentz invariant ψψ and generally involves a nontraditional potential.It is worth mentioning that the mass term m ψψ will be vanished for a massless spinor field.Before exploring solutions of the spinor field coupling with gravitation, we will begin with an illuminating solution for the spinor field in the absence of gravitation.

Exact solutions in flat spacetime
Solutions of a free spinor field in four-dimensional flat spacetime have been obtained in the literature.We intend to study solutions of a spinor field in five-dimensional flat spacetime, which may provide us insights into the subject.We use coordinates x µ and y to label the four-dimensional hypersurface and an extra spatial dimension, respectively.The five-dimensional spacetime is characterized by the background metric where η µν is the four-dimensional Minkowski metric, and y = x 5 is the extra-dimensional coordinate.The spinor field is set to depend only on the extra spatial dimension y.
Normally, a spinor field has four complex components in five-dimensional spacetime.For the sake of simplicity, however, we assume that the solution of a real spinor field in five-dimensional bulk spacetime has the following form [33] ψ(y) = (p(y), 0, q(y), 0) T , ( which has two nonvanishing real components p(y) and q(y).Here, "T" denotes the transposition of a matrix.The Lorentz invariant ψψ for the real spinor field assumption (3.2) is given by ψψ = 2p(y)q(y).
Substituting the real spinor field assumption (3.2) and the potential term (2.14) into Eq.(2.3), the equations of motion for the spinor field are reduced to where the prime denotes the derivative with respect to the extra-dimensional coordinate y.For the ansatz (3.1) and the assumption (2.14), we find the following solution Here, c 1 and c 2 are integration constants and f (2c 1 ) ≡ ∂V ∂( ψψ) .It is obvious that the Lorentz invariant ψψ = 2p(y)q(y) = const.For the case of a free massive spinor field with the solutions of p(y) and q(y) can be written as where m is the mass parameter of the spinor field ψ.
One can derive a first integral from Eq. (3.3) as where C 1 is an integration constant.It is clear that the solution (3.4) and the simpler one (3.6)satisfy the first integral (3.7).The energy density of the spinor field with respect to a static observer U M = (−1, 0, 0, 0, 0) is generally given by Further, we will investigate exact solutions of a spinor field in the presence of gravitation.

Spinor walls in warped spacetime
We explore solutions of a spinor field minimally coupling with gravitation in five-dimensional asymptotically AdS spacetime in which the extra spatial dimension is curved ("warped").Such system differs from the system composed of gravitation and a scalar field.Generally speaking, a coupled system consisting of different components of the spinor field is complicated, especially in the case of the spinor field coupled to gravitation.For the sake of simplicity, a real spinor field will be taken into account.Considering the Z 2 symmetry of the hypersurface and the energy density distribution of the spinor field, we speculate the warp form of the hypersurface and the abstract form of the spinor field solution.Further, we will employ a reconstruction technique that the action of the spinor field is obtained according to the above speculations regarding the hypersurface and the spinor field solution.As in section 3, we suppose that the real spinor field possesses an abstract form as Eq.(3.2).

The setup
Within the context of general relativity, the bulk Einstein equation is In this paper, we are interested in static flat walls with the four-dimensional Poincaré symmetry, for which the general form of the metric (the general ansatz for a domain wall) is given by ds 2 = a 2 (y)η µν dx µ dx ν + dy 2 , ( where a(y) is the warp factor, η µν is the four-dimensional Minkowski metric, and y = x 5 is the extra-dimensional coordinate.Note that the warp factor a(y) and the background spinor field ψ(y) are merely functions of y for static flat walls.Using the metric ansatz (4.2), the spin connection ω Ā B M (2.9) and the generator of spinor representation S Ā B (2.10) can be worked out.The nonvanishing components of With the above metric (4.2) and spin connection (4.3), Eqs.(4.1) and (2.12) are reduced to and respectively.From Eq. (4.4), we have From Eq. (4.5), the following first integral can be found where C 2 is an integration constant.The first integral (4.7) is a generalization of Eq. (3.7) in the presence of gravitation.The Lorentz invariant does not depend on the form of the potential V ( ψψ).This implies that the dynamics of the spinor field is different from that of a scalar field, especially in the presence of gravitation.In warped spacetime, the first integral (4.7) will be a significant constraint on the solutions of the spinor field.
It is worth pointing out that, however, only three equations are independent in this system in the light of the contracted Bianchi identity.However, in Eqs.(4.4) and (4.5), arbitrary three equations are not completely independent.Thus, we need to consider the rest equation of them.We take Eqs.(4.4a), (4.5a), and (4.5b) as the beginning of solving this system.We will see that the solutions of Eqs.(4.4a), (4.5a), and (4.5b) have similar characteristics as the solutions in flat spacetime.Actually, one can also obtain a set of solutions by Eqs.(4.5) and (4.6).However, these solutions are equivalent to the ones which we will present below.
Since this system satisfies the equations of motion, the energy density could be of the general form For a static observer U M = (−1/a(y), 0, 0, 0, 0) lying on the wall, the distribution of the observed energy density along extra dimension y would be It is worth noting that in this energy density, the contribution of the cosmological constant need to be subtracted, so that the minimum energy density (if it exists) vanishes.
There are two significant parameters for domain walls with Z 2 symmetry [40], which will be employed in the following discussions.
• The first one is the tension of the wall, which measures the energy per unit volume on the wall.The tension of the wall is defined as the integral of the energy density along the coordinate y where g is the determinant of the metric and Σ is the integration interval of the fifth dimension y.For an infinite extra dimension, the integration interval will be (−∞, ∞).Generally speaking, if the tension of the wall is finite, the solution of the wall is regular.
• The second one is the thickness of the wall, which evaluates the typical length scale of the scale variation of the wall.The thickness of the wall in the domain wall case is defined as the reciprocal of the factor k appearing in front of the y.We denote the thickness of the wall by σ = 1/k and the above definition leads to where c ranges from 0% to 100%.For the following Z 2 domain wall solution in flat spacetime c is approximately 64%.Here, k is a parameter and the length scale 1/k is related to the thickness of the wall.
By means of the curvature of the bulk, the bulk spacetime could be classified into two kinds: the constant and nonconstant curvature spacetimes.We will discuss solutions in these two kinds of spacetimes in the following two subsections.

Nonconstant curvature solutions
We will investigate exact solutions of a spinor field in warped spacetime.Firstly, we begin with a given bulk spacetime which is characterized by Eq. (4.2) with the warp factor a(y) = sech n (ky) . (4.13) An obvious property is that the warp factor (4.13) is smooth rather than singular at the location of y = 0.Such property indicates that the curvature singularity at the location of the codimension one hypersurface will be eliminated.The profile of the warp factor (4.13) for positive n is shown in Fig. 1a.The bulk spacetime with the warp factor (4.13) with nonzero n is asymptotically AdS at the boundary of the extra dimension.Generally speaking, the n in this warp factor (4.13) could be negative, zero, and positive.As we will demonstrate in section 6, if n is positive, a localized graviton zero mode can be obtained.If n = 0, the bulk spacetime is flat.Nevertheless, we also discuss the case of negative n, because it has similar behaviour to the Janus solution [41], one of the early examples on the study of defects and interfaces in field theories.By solving the equations of motion, we find a solution of spinor domain wall generated by the spinor field (3.2) with the following components p(y) and q(y) ) and the potential where we have defined a new function D(y) as and c 1 2c 2 is identified as the mass parameter (we can also see it from Eq. (3.6)).Here, sgn is the sign function and B is the incomplete beta function.According to the assumption (3.2), the function D(y) (4.16) should be real.Due to this incomplete beta function B cosh 2 (ky) −2n − 1 2 , 1 2 is a complex function, we need to find the conditions that the function D(y) (4.16) is real.If we take n to be zero, positive integers, and positive half-integers, the function D(y) is real.In addition to these values, the function D(y) is complex.In this case, the imaginary component should be discarded because it is introduced in the process of solving the equations, where the variable y has been analytically continued to the complex domain.From the perspective of the equations of motion, the resulting solution of a real differential equation will be real.In essence, the new defined function D(y) can always be a real function, owing to the real component of the incomplete beta function B can be canceled by using the integration constant c 3 .For the potential term (4.15), this spinor field has a nonlinear potential as a correction to the linear term ψψ.The additional constant in the potential (4.15) could be accounted for the effective cosmological constant.In fact, the equations of motion (4.4) and (4.5) require that the integration constant c 1 is zero.For comparison with the solutions in flat spacetime, we have retained the integration constant c 1 so that c 1 = 0 is a special case.
To demonstrate the behaviour of the warp factors and the corresponding energy densities for the following spinor field solutions, we will henceforth introduce the dimensionless quantities ỹ = ky and ρ (ỹ) = ρ(ỹ)/k 2 in the figures below to display their profiles.
For the warp factor (4.13), the energy density is calculated as ρ(y) = 3k 2 n(2n + 1)sech 2 (ky) .(4.17) If n takes 0 or − 1 2 , the energy density will be zero.The warp factor (4.13) and the corresponding energy density of the spinor field are depicted in Fig. 1.Note that, we have deducted the contribution of the effective cosmological constant so that the minimum (c) Warp factor a(ỹ) Figure 1: Plots of a(ỹ) and ρ (ỹ), where a(ỹ) is the warp factor in Eq. (4.13) and ρ(ỹ) is the energy density in Eq. (4.9).The parameter n is set to 1, − 1 4 , and −1.The warp factors (left panels) with a blue solid line, orange dashed line, and green dashing line correspond to the energy densities (right panels) with a blue solid line, orange dashed line, and green dashing line, respectively.value of the energy density is zero.As shown in Fig. 1, for − 1 2 < n < 0, the energy density increases with |y| and approaches to a constant as ỹ → ±∞.There is a trap of the energy density at ỹ = 0.For n < − 1 2 and n > 0, the energy density decreases with |y| down and approaches to zero as ỹ → ±∞.The energy density of the spinor field centers on the vicinity of the hypersurface.The more concentrated the energy density of the spinor field is, the more warped the extra dimension is.There is a wall between two sides of the hypersurface.The wall is similar to a domain wall, a codimension one soliton generated by scalar fields.This solution of the spinor field can be viewed as a thick domain wall interpolating between two asymptotically AdS vacua.
In the case of the warp factor (4.13), according to (4.10), we obtain the tension of the wall in the domain (−∞, ∞) for n > 0 where Γ is the Euler gamma function.The tension of the trap in the domain where the special function 2 F 1 is the hypergeometric function.The tension of the wall in the domain The integral (4.10) is divergent for n < − 1 2 and − 1 2 < n < 0, and is convergent for n > 0. Consequently, the solution (4.14) and (4.15) is globally regular for n > 0 and is locally regular for − 1 2 < n < 0 and n < − 1 2 .For the case where the bulk field is a scalar field, the profile functions of the Jauns solution [42] typically are Here, the parameter γ ∈ 0, 1 √ 2 .In the Janus solution, the scalar field solution only varying spatially resembles a domain wall which is a typical soliton configuration.There is a possibility of a wall-like configuration for the spinor field.The warp factor (4.21) is also smooth at the location of y = 0 and the curvature singularity at the location of the codimension one hypersurface will be eliminated.A common property of the warp factors (4.13) and (4.21) is that the bulk spacetime is asymptotically AdS at the boundary ỹ → ∞.
As a deformation of the warp factor (4.13), the Janus solution (4.21) is rather difficult to solve the equations of motion (4.4) and (4.5).For the convenience of solving the equations of motion, we take β ≡ 1 − 2γ 2 , β ∈ [0, 1] and the cosh(2ky) is replaced by the cosh(ky) in Eq. (4.21).Therefore, we begin with the warp factor of the form One can obtain the p(y) and q(y) of Eq. (3.2) as follows and the potential of the spinor field where the function E(y) is given by Here, c 1 2c 2 could also be viewed as the mass of the spinor field.In the potential (4.25), there are nonlinear terms as a correction and an extra constant term could be accounted for the effective cosmological constant.It is worth stressing that according to the requirement of the equations of motion, the integration constant satisfies c 1 = 3 4 k 2 β 2 (4 + β 2 ).Nevertheless, we still retain c 1 as an arbitrary constant to compare with the solutions in flat spacetime.We depict the warp factor (4.23) and the corresponding energy density in Fig. 2. For the warp factor (4.23) which is similar to the warp factor (4.13) with n = −1, the corresponding solution is not globally regular.From the aforementioned results, one can find that there is no singularity in the geometry of the bulk.This resembles a thick domain wall.A remarkable property is that this solution has a holographic interpretation in terms of the theory of living on the boundary [41].
The second analytic solution of the spinor field in warped spacetime will be exhibited in the following context.We choose the following warp factor [4] a(y) = cos n (ky) . (4.27) The warp factor is smooth at the position of the hypersurface as well.The profile of the warp factor is shown in Fig. 3.This indicates that the curvature singularity is eliminated at the location of the codimension one hypersurface.For the warp factor (4.27) with negative n, the bulk spacetime will be singular at the boundary ỹ = ± π 2 .Therefore, the range of ỹ is from − π 2 to π 2 .The case n = 0 corresponds to a flat bulk spacetime.We find a solution of the spinor field (3.2) with the following p(y) and q(y) and the potential is determined by where the function F (y) is given by and c 1 2c 2 is viewed as the mass.In this case, the incomplete beta function is real.Therefore, the function F (y) (4.30) is a real function which satisfies the assumption (3.2).Similar to the first solution, the equations of motion (4.4) and (4.5) require that the integration constant c 1 is zero.For comparison with the solutions in flat spacetime, we have retained the integration constant c 1 so that c 1 = 0 is a special case.In the case of the warp factor (4.27), the energy density is computed as If n takes 0 or 1 2 , the energy density will be a constant.The warp factor and the energy density of the spinor field are plotted in Fig. 3.The energy density at the ỹ direction boundaries − π 2 and π 2 is positive infinity for 0 < n < 1 2 and is negative infinity for n < 0 and n > 1  2 .Nevertheless, from the local perspective, these configurations can be seen as a local wall for n < 0 and n > 1  2 or a local trap for 0 < n < 1 2 .In line with the definition (4.10), the tension of the wall in the domain (−σ/2, σ/2) for n < 0 and n > 1  2 is given by The tension of the trap in the domain This result implies that the wall or trap is locally regular.If the thickness of the wall is extended to the boundary of the extra dimension σ = π, the wall will be singular.
The integral (4.10) is divergent for any n, because the energy density at the ỹ direction boundaries − π 2 and π 2 is positive infinity for 0 < n < 1 2 and is negative infinity for n < 0 and n > 1  2 .Nevertheless, from the local perspective, the integral (4.10) is locally convergent and regular.Therefore, these configurations can be seen as a local wall for n < 0 and n > 1  2 or a local trap for 0 < n < 1 2 .We will show the third solution which displays an analogous behaviour of a warped bulk compared with the first solution.We choose the warp factor as follows a(y) = exp −(ky) 2 . (4.34) The warp factor possesses a simplified form of the normal distribution and is smooth at the position of the hypersurface.Therefore, there do not exist singularities in the bulk.In this warped spacetime, a spinor domain wall is generated by the spinor field (3.2) with the following p(y) and q(y) and the potential is determined by where the function G(y) is given by with erf the error function.Obviously, the function (4.37) is a real function in the real domain.Similar to previous solutions, when c 1 is zero, this solution satisfies the equations of motion (4.4) and (4.5).
Substituting the warp factor into Eq.(4.9), one obtains the energy density as follows The warp factor and the energy density of the spinor field are depicted in Fig. 4. As shown in Fig. 4, the energy density approaches to the negative infinity as ỹ → ±∞.From the local perspective, there is a local wall at the position ỹ = 0, because the boundary of the bulk along ỹ can not be reached physically.It is worth mentioning that we always set the minimum value of the energy density of the local wall as zero so that the vacuum energy density is zero and an effective cosmological constant is introduced.
In accordance with (4.10), the tension of the wall in the domain Such integral (4.10) of the energy density does not converge in the domain (−∞, ∞).It implies that the wall is locally regular.If the thickness of the wall is extended to the boundary of the extra dimension σ = ∞, the wall will be singular.We have investigated spinor wall solutions in the bulk with nonconstant curvature.Among these solutions, constant curvature cases have been included.In the following subsection, we will specifically discuss them.All these solutions definitely satisfy the first integral (4.7).As mentioned above, we retain the integration constant c 1 so that c 1 = 0 (c 1 = 3 4 k 2 β 2 (4 + β 2 ) for the Janus solution above) is a special case which determinately satisfies the equations of motion (4.4) and (4.5).The solutions in curved spacetime can be regarded as a correction to the ones in flat spacetime.Yet from the perspective of the solutions with c 1 = 0 (c 1 = 3 4 k 2 β 2 (4 + β 2 ) for the Janus solution above), the solutions in curved spacetime are rather different from the ones in flat spacetime since the mass term in these potentials of the spinor field actually is vanishing.Therefore, the solutions in curved spacetime could also not be considered as a correction on the ones in flat spacetime.
It is worth noting that the three sets of solutions of the spinor field are static.All the solutions of the spinor field can lead to a wall-like configuration.The wall solutions with the warp factors (4.13) (including the deformation (4.23)) and (4.27) are smooth at the location of the hypersurface and have a well-defined thin wall limit.The wall solution with the warp factor (4.34) is also smooth at the location of the hypersurface but has a definite thickness.Moreover, in the case of the warp factor (4.13), its energy density mainly localizes in the vicinity of the hypersurface even though the extra dimension is infinite.In the case of the warp factors (4.27) and (4.34), their energy densities are negative infinity at boundaries and there are local walls from the local perspective.
Whether these soliton-like walls are stable or not is still crucial.Localization of gravitation is also worth investigating since it relates to effective gravitational theories on these walls.These issues will be addressed in the following sections.

Constant curvature solutions
Having exhibited three sets of solutions for the spinor field, we intend to discuss special solutions which correspond to three types of bulk spacetimes with constant curvature, namely, zero scalar curvature, de Sitter (dS), and AdS bulk spacetimes.One can verify that the three types of five-dimensional bulk spacetimes have three kinds of metrics characterized by Eq. (4.2) with analytic warp factors [43], which are listed in Tab. 1.

Warp factors a(y) Scalar curvature R(y)
Bulk spacetimes  We plot their profiles of the warp factors in Fig. 5a.The warp factors cos has an inflection point at ỹ = 0 and there is a scalar curvature singularity at ỹ = 0.The bulk spacetime is cut into two identical constant curvature parts by a hypersurface.The sharp point of this warp factor leads to a δ-function in the second-order equations of motion for the general relativity case, which resembles the RS II brane world model.This case can be interpreted as the appearance of an exotic wall.
We will present the three solutions for the zero scalar curvature, dS, and AdS bulk spacetimes, respectively.As in the case of nonconstant curvature solutions, below we still retain the integration constant c 1 so that c 1 = 0 is the case satisfying Eqs.(4.4) and (4.5).For the bulk spacetime with zero scalar curvature, there is a spinor field solution with the following p(y) and q(y) and An obvious property for this solution is that there is an additional nonlinear term in the potential.As shown in Fig. 5b, one can find that the energy density of the spinor field for this solution is infinite at the position ỹ = 0. Therefore, a scalar curvature singularity occurs at the position of ỹ = 0.The energy density approaches to zero as ỹ → ±∞.We could believe that an exotic spinor wall exists at ỹ = 0.According to the definition (4.10), the tension of the exotic wall in the domain The tension of the exotic wall is not globally regular, but is locally regular.
For the dS bulk spacetime, a spinor field solution is given by and with the redefined function where 2 F 1 is the hypergeometric function and is real.Therefore, H(y) is a real function.
As shown in Fig. 5c, the energy density for this solution is valid for ỹ ranging from − π 5 to π 5 , and mainly distributes in the domain − π 5 , π 5 with two divergence points of the energy density at the boundaries.This configuration can be viewed as a local trap.
The tension of the trap in the domain ; sin 2 5kσ 4 .
The tension of the trap is not globally regular, but is locally regular.For the AdS bulk spacetime, an analytic solution is obtained as and with the redefined function ; cosh 2 5ky 2 , ( where the hypergeometric function 2 F 1 is complex.Therefore, the redefined function I(y) is complex.The imaginary component of this complex function should be discarded because it is introduced in the process of solving the equations, where the variable y has been extended to a complex domain.As shown in Fig. 5d, the energy density increases and approaches to a constant as ỹ → ±∞.There is a trap of the energy density at ỹ = 0.
According to (4.10), the tension of the trap in the domain The tension of the trap is locally regular.So far we have investigated spinor wall solutions with constant curvature.These solutions satisfy the first integral (4.7).However, from the perspective of the energy density, the solution with zero scalar curvature leads to an exotic wall with an infinite large energy density at the position of the hypersurface.The energy density for the solution with zero scalar curvature has a singularity since the corresponding bulk spacetime has a scalar curvature singularity at ỹ = 0.The energy density for the dS solution is not bounded from above and this solution is not globally regular.The energy density for the AdS solution converges to a constant as ỹ tends to the boundary and this solution is also not globally regular.Nonetheless, from the local perspective, the energy densities of the dS and AdS solutions have local traps.Moreover, the stability of the above walls and the localization of gravitation need to be considered.

Stability of spinor walls
Having found spinor wall solutions, we are obliged to examine perturbations about this background with the purpose of studying their stability.In general, according to the four-dimensional Poincaré transformations, the linear perturbations of the metric (4.2) and the energy-momentum tensor can be mathematically decomposed into transversetraceless (TT) tensor mode, transverse vector modes, and scalar modes.In our case, the perturbation of the metric will couple to the perturbation of the bulk spinor field.The spinor field here behaves like a scalar field in spacetime.Thus, the perturbation of the bulk spinor field is similar to the perturbation of a scalar field.As a consequence, the stability of the spinor walls can be analyzed with the help of methods in [8][9][10][11].
For the sake of convenience, we introduce the transformation dz = dy/a(y) for the metric (4.2).The linear perturbation to the metric (4.2) in the conformal coordinate z can be parameterized as δ (1) where δ (1) signifies linear perturbation.Its inverse will be The TT tensor mode h µν satisfies The vector modes ξ µ and ζ µ are both transverse The four components θ, ω, φ, and χ in Eq. (5.1) are scalar modes of the metric perturbation.
The vielbeins can be separated into their background and perturbation parts: The unperturbed vielbeins and its inverse are respectively According to Eq. (5.1), the corresponding perturbation part of the vielbeins1 can be parameterized as δ (1) e N M = a(z) and the inverse is where ǫ αβγ is the permutation tensor, φ 3 is a pseudoscalar, ζ ι 3 is a pseudovector, and (5.9) The spinor field can be separated into its background and perturbation parts: 1) ψ . (5.10) The evolution equations of the linear perturbations are Rδ (1) g M N = δ (1) T M N , (5.11) 1) ψ = 0 , (5.12) .13)and the perturbation of the energy-momentum tensor is given by δ (1) T M N = 1 2 δ (1) ψγ (M ∇ N ) ψ + ψδ (1)  1) ψ . (5.14) Here, the perturbations of the gamma matrices and the covariant derivative of the spinor field ψ are respectively.Using Eq. (2.9) the spin connection in the linear perturbation can be computed as: (5.17) The nonvanishing components of the perturbed spin connection are δ (1) δ (1) ω 5 β µ =δ (1) For Eqs. (5.12) and (5.13), the equation of motion of the tensor perturbation leads to the following constrains δ (1) e µ ν γ ν ∇ µ ψ + e µ ν γ ν δ (1) ω µ ψ = 0 , (5.20) ∇ µ ψδ (1) e µ ν γ ν − ψδ (1) ω µ e µ ν γ ν = 0 . (5.21) In fact, the perturbation δ (1) ψ of the spinor field ψ will contribute to the scalar modes of the perturbation.A complete investigation of the stability of this system includes tensor mode, vector modes, and scalar modes.For the sake of simplicity, we only consider the tensor perturbation here.As a consequence, the perturbation of the energy-momentum tensor is vanishing for the tensor perturbation.We eliminate the gauge degrees of freedom not by choosing a gauge but by working only with gauge-invariant quantities.The tensor mode h µν is gauge-invariant [8].We obtain the following tensor perturbation equation (4) where (4) = η µν ∂ µ ∂ ν is the four-dimensional d'Alembert operator.
After performing the separation of variables, the equation of motion of the tensor mode is separated into the Klein-Gordon equation for the four-dimensional part and a secondorder differential equation for the extra-dimensional part.For the extra-dimensional part of the tensor mode, we can transform it into a canonically normal mode and obtain a Schrödinger-like equation for the extra-dimensional part [10] −∂ 2 z + W (z) Ψ(z) = m 2 Ψ(z) . (5.23) Here, the effective potential W (z) has the following form with Ω(z) = 1 2 ∂ z ln a 3 (z) . (5.25) Consequently, the Schrödinger-like equation (5.23) can be factorized as a supersymmetric quantum mechanics form where the "Hamiltonian" Q(z) Q † (z) is a positive definite Hermitian operator and there are no normalizable negative energy graviton modes.The supersymmetric form of the Schrödinger-like equation demonstrates that any tachyonic modes with m 2 < 0 are absent for the tensor mode.Therefore, these spinor walls can exist stably under the tensor perturbation in such spacetime with the metric (4.2).It is worth stressing that this result is valid only for a regular bulk spacetime without singularities in curvature invariance.From this perspective, the stability of the spinor wall solution in the spacetime with zero scalar curvature is still unknown.The spinor wall solution in the spacetime with zero scalar curvature represents an exotic wall which has a singularity at the location of the wall.

Localization of gravitation
In what follows, we will discuss localization of the zero mode of the tensor perturbation for different warp factors investigated in section 4. In the process of obtaining a four-dimensional effective gravitational theory from high dimensions, a four-dimensional effective coupling constant could be obtained.To quadratic order in the tensor perturbation, the gravitational action would be (5) 2 d 5 x δ (2) −g (5) R + −g (5) δ (2) R + δ −g (5) δR , ( where δ (2) denotes second-order quantities in the tensor perturbation.Taking the TT condition (5.3) and the separation of variables in the coordinates (x µ , z) into account, one obtain a four-dimensional effective theory with the Einstein-Hilbert action and additional terms: The four-dimensional part of such canonical action is anticipated to describe gravitation observed experimentally.The four-dimensional effective Planck mass M Pl is expressed as follows Here, Σ is the integration interval of the fifth dimension z.For an infinite extra dimension, the integration interval will be (−∞, ∞).An arbitrary KK mode is normalizable if and only if the following normalization condition is satisfied In order to obtain the four-dimensional Newtonian potential, the tensor zero mode is required to be localized.From Eq. (5.23), we obtain a zero eigenvalue solution [45,46], namely the zero mode of the graviton Here, C 1 and C 2 are integration constants.There exists the possibility of the normaliza-tion of the zero mode with this general solution.Furthermore, the graviton zero mode propagating in the bulk is subjected to boundary conditions [47].For our case, the vanishing variation of the action of the zero mode at the boundary ∂Σ of the fifth dimension z leads to the boundary condition Notice that h µν is assumed to vanish at the four-dimensional boundary x ρ → ±∞, and corresponding boundary terms are automatically zero.Accordingly, the boundary condition (6.6) can again be satisfied by imposing either Dirichlet condition: or Neumann condition: The zero eigenvalue solution of the Schrödinger-like equation (5.23) is required to satisfy these boundary conditions.The Dirichlet condition allows the zero mode with the second particular solution The Neumann condition leads to the zero mode with the general solution (6.5).Especially, the first particular solution identically meets the Neumann condition, which is independent of the concrete form of the warp factor.From this result, the solutions of the warp factor that we have investigated above coincide with the requirement of the Neumann condition.Essentially, one does not rule out the possibility that the second particular solution (6.9) satisfies the boundary condition (6.6) for some warp factors.In the following, we discuss the localization of the graviton zero mode only for the first particular solution (6.10).
For an AdS bulk or an asymptotically AdS one, the warp factors may be exponentially divergent at the fifth-dimensional boundary z → ±∞ and the first particular solution (6.10) presumably is not normalizable.Thus, only the (asymptotically) properties of the bulk geometry are not sufficient.The normalization of the first particular solution (6.10) will be discussed in the following.Inserting (6.10) into the normalization condition (6.4), one obtains Here, Σ ′ is the integration interval in the y coordinate.
For the warp factor (4.13), we have ∞ −∞ sech 2n (y) dy = 4 n kn 2 F 1 (n, 2n; n + 1; −1) .(6.12) As a result, if n > 0 the integral (6.11) is convergent and the zero mode of the graviton (6.10) is normalized.The warp factor in the Janus solution (4.21) or (4.23) leads to a divergent result of the integral (6.11) and the graviton zero mode is not localized.We also conclude that for the warp factor cosh 2 ky corresponding to the AdS constant curvature spacetime, the graviton zero mode is not localized.For the warp factor (4.27), the result will be Consequently, if n > − 1 2 the integral (6.11) is convergent and there is a normalized graviton zero mode.Therefore, the warp factor for the spinor wall solution with the dS bulk spacetime achieves a normalizable graviton zero mode.We conclude that for the warp factor cos 2 5 5 2 ky corresponding to the dS constant curvature spacetime, the graviton zero mode is localized.For the warp factor (4.34), the result is In this case, the zero mode of the graviton is localized.From these results, spinor walls with warp factors (4.13), (4.27), and (4.34) can localize the graviton zero mode under specific conditions.To summarize, under our assumptions about the background, it is shown that the gauge-invariant perturbation corresponding to the tensor mode of the geometry is localized on the walls.The localized zero mode of the tensor perturbation results in the four-dimensional general relativity and hence the four-dimensional Newtonian potential.The massive KK modes of the tensor perturbation will make a correction ∆U(r) ∼ 1/r 2 concerning the Newtonian potential [2,6].

Discussion and conclusions
In the above solutions, we found the solution corresponding to the warp factor (4.13) with n > 0 satisfies the conditions that the energy density is regular; the solution is stable under the tensor perturbation; and the graviton zero mode can be localized on the wall.For the warp factor (4.13), there exists another solution, for which the energy density is regular; the solution is stable under the tensor perturbation; while the graviton zero mode cannot be localized on the wall.The solution corresponding to the warp factor (4.23) also belongs to this case.Moreover, we found an interesting result that a global trap can stably exist, but its energy density is not regular and the graviton zero mode cannot be localized on the trap.For the solutions with the warp factors (4.27) and (4.34), we believe that local walls or traps can exist stably, whether or not the graviton zero mode is localized.The constant curvature solutions in the dS and AdS bulk spacetimes are special cases of the solutions with the warp factors (4.34) and (4.13) respectively.The constant curvature solution in the bulk spacetime with zero scalar curvature corresponds to an exotic wall.
We comment on the divergences in the above cases.For the constant curvature solution in the bulk spacetime with zero scalar curvature, the wall which has a curvature singularity at the origin can be considered as a thin wall.For naked singularities at the boundaries of the bulk, the walls/traps can be physically meaningful.One reason is that in renormalization group flows to nonconformal theories, the AdS horizon is replaced by a naked singularity.This singularity is physical in the sense that the singular behavior corresponds to strong coupling effects such as confinement or screening in the boundary theory.The other reason comes from analyzing the spectrum of gravity from a fourdimensional point of view.The KK states couple to matter on the thin domain wall and cause small violations of the four-dimensional conservation of energy and momentum.Spaces with naked singularities are physically acceptable only if one imposes boundary conditions that guarantee the four-dimensional conservation of energy and momentum.
In summary, we have investigated analytic solutions of a spinor field in five-dimensional spacetime.We mainly focused on analytic solutions of the spinor field for flat and warped spacetimes.First, we considered a toy model in which a real spinor field provides an illuminating solution.We found a set of solutions in the absence of gravitation and three sets of solutions of the spinor field in the presence of gravitation.These solutions with gravitation are discussed for three kinds of warp factors which correspond to different warped geometries.Further, we discussed domain wall solutions of the spinor field in the five-dimensional zero scalar curvature, dS, and AdS bulk spacetimes.In order to identify stability of these walls, we studied perturbations of these wall systems.The proof demonstrate that these spinor walls are stable under the tensor perturbation.Moreover, it was shown that the zero mode of the graviton can be localized on the walls in the fivedimensional bulk spacetime with nonconstant curvature.It implies that the Newtonian potential is recovered.
In many works in which analytical solutions of domain walls or thick branes were found, scalar fields were widely considered because not only it is easy to be solved but also matter fields can be localized by employing a natural mechanism.In this work, we attempted to explore the possibility of domain walls generated by a real spinor field.We found that, if such field exists, appropriate choices of the spinor potential would generate a rich variety of behaviours, quite different from their widely studied scalar field counterparts.The energy density of the spinor field concentrates in the vicinity of walls except for traps.Moreover, nonlinear spinor fields may shed light on torsion of the bulk spacetime, and present a torsionless four-dimensional hypersurface which may stands for a world with the illusion of torsion free [25].
Whether these spinor wall solutions can be seen as branes or not is yet important to verify.Branes require that matter fields are confined on them or bounded in the vicinity of them.Therefore, the dynamics of matter fields in the presence of a wall configuration is also a topic worth studying.In this work, we studied solutions for a real spinor field.However, from the perspective of fundamental particles, fermions in the standard model are described by complex spinors with four components.It is interesting to investigate solutions for a complex spinor field in warped five-dimensional geometries.So far, we only considered a simple Lorentz invariant ψψ for a spinor field in our solutions.However, other Lorentz invariants in the potential term V ( ψ, ψ) of the Lagrangian density, such as ψγ 5 ψ, current-current interaction density ψγ µ ψ ψγ µ ψ , and ψγ 5 γ µ ψ ψγ 5 γ µ ψ , are also rational.We will leave this in future works.In this work, we only considered a wall with the four-dimensional Poincaré symmetry.Naturally, it is worth exploring other walls which possess other symmetries, such as the four-dimensional dS or AdS symmetry.

Figure 3 : 4 , and − 1 4 .
Figure3: Plots of a(ỹ) and ρ (ỹ), where a(ỹ) is the warp factor in Eq. (4.27) and ρ(ỹ) is energy density in Eq. (4.9).The parameter n is set to 1,1  4 , and − 1 4 .The warp factors (left panels) with a blue solid line, orange dashed line, and green dashing line correspond to the energy densities (right panels) with a blue solid line, orange dashed line, and green dashing line, respectively.

5 2 5 5 2 , π 5 . The warp factor 5 2 ky 2 5
ky are smooth in their domain of definition.The warp factor cos 2 ky is only well-defined in the interval ỹ ∈ − π 5 ρ (ỹ) for the AdS bulk spacetime

Figure 5 :
Figure 5: Plots of the warp factor a(ỹ) and the corresponding dimensionless energy density ρ (ỹ), where a(ỹ) is the warp factors in Tab. 1 and ρ(ỹ) is the energy density in Eq. (4.9).The warp factors (the upper left panel) with a blue solid line, orange dot-dashed line, and green dashed line correspond to the energy densities (the rest panels) with a blue solid line, orange dot-dashed line, and green dashed line, respectively.

Table 1 :
Table for three background spacetimes and their corresponding warp factors. 2AdS