Massless fermions localization on domain walls

TL;DR

This paper studies how massless fermions can be localized on scalar domain walls in a 5D Anti de-Sitter spacetime, revealing a normalizable state influenced by Yukawa interactions and gravitational effects.

Contribution

It demonstrates the existence of a normalizable fermion state on domain walls, independent of $Z_2$ symmetry, due to the interplay of Yukawa coupling and gravity.

Findings

Normalizable fermion state identified on domain walls.

Fermion migration toward the wall edge due to Yukawa-gravity competition.

Localization effect independent of $Z_2$ symmetry.

Abstract

Massless fermions on scalar domain walls are considered. Two walls are established, corresponding to 5-dimensional static spacetime asymptotically Anti de-Sitter, differentiated by the symmetry around the wall, and in each case massless chiral fermions are coupled to the wall by a Yukawa term. We identify a normalizable state associated to the migration of fermions toward the edge of the wall. This effect is generated by the competition between the Yukawa interaction and the gravitational repulsion on the matter fields, and it is independent of the $Z_2$ symmetry of the wall.

Equations

$$S=\int d{x}^5\sqrt{g}\left[\frac{1}{2}R-\frac{1}{2}{g}^{ab}{\nabla }_a\varphi {\nabla }_b\varphi -V(\varphi )\right],$$

$$d{s}^2=e^{2A(y)}{\eta }_{\mu \nu }d{x}^{\mu }d{x}^{\nu }+d{y}^2,$$

$$\rho =-3(2A^{\prime 2}+A^{\prime \prime }).$$

$${\Gamma }^a{\nabla }_a\Psi (x,y)=\lambda \varphi (y)\Psi (x,y).$$

$$\Psi (x,y)={\Psi }_{\text{L}}(x){\psi }_{\text{L}}(y)+{\Psi }_{\text{R}}(x){\psi }_{\text{R}}(y),$$

$$i{\gamma }^{\mu }{\partial }_{\mu }{\Psi }_{\text{R}}^{\text{L}}(x)=0.$$

$$\left({\partial }_y+2A^{\prime }(y)\pm \lambda \varphi (y)\right){\psi }_{\text{R}}^{\text{L}}(y)=0;$$

$${\psi }_{\text{R}}^{\text{L}}(y)=e^{-2A(y)\mp \lambda \int \varphi (y)dy}.$$

$${\psi }_{\text{L}}(y)⟶e^{(2k-\lambda {\varphi }_0)|y|},$$

$$\lambda >\frac{2k}{{\varphi }_0}\equiv {\lambda }_1$$

$${\int }_{-\infty }^{\infty }dy{\psi }_{\text{L}}(y)\simeq {\left[\frac{2\pi }{-2|A^{\prime \prime }(0)|+\lambda {\varphi }^{\prime }(0)}\right]}^{1/2}{\psi }_{\text{L}}(0).$$

$$\lambda >2\frac{|A^{\prime \prime }(0)|}{{\varphi }^{\prime }(0)}\equiv {\lambda }_2$$

$$A(y)=-\delta \ln \cosh \frac{\alpha y}{\delta },$$

$$\varphi (y)={\varphi }_0\arctan \sinh \frac{\alpha y}{\delta },{\varphi }_0=\sqrt{3\delta },$$

$$V(\varphi )=3{\alpha }^2\left[\left(\frac{1}{2\delta }+2\right){\cos }^2\frac{\varphi }{{\varphi }_0}-2\right].$$

$${\psi }_{\text{L}}\sim {\cosh }^{2\delta }\left(\alpha y/\delta \right)\exp \left[\lambda \frac{\delta {\varphi }_0}{\alpha }\left(\frac{\alpha y}{\delta }(2\mathrm{arccot}{e}^{\alpha y/\delta }+\arctan \sinh \alpha y/\delta )\right)\right]$$

$$\times \exp \left[i\lambda \frac{\delta {\varphi }_0}{\alpha }\left(\mathrm{L}{\mathrm{i}}_{2}2(-i{e}^{-\alpha y/\delta })-\mathrm{L}{\mathrm{i}}_{2}2(i{e}^{-\alpha y/\delta })\right)\right]$$

$$\mathrm{L}{\mathrm{i}}_{\mathrm{n}}(\tau )\equiv \sum _{k=1}^{\infty }\frac{{\tau }^k}{k^n};$$

$${\lambda }_1=\frac{4\alpha }{\pi \sqrt{3\delta }}$$

$${\lambda }_2=\frac{\pi }{2}{\lambda }_1,$$

$${\psi }_{\text{L}}(y)⟶\Theta (y)e^{(2k_{{}_{+}}-{\lambda }_{{}_{+}}|{\varphi }_{{}_{+}}|)y}+\Theta (-y)e^{-(2k_{{}_{-}}-{\lambda }_{{}_{-}}|{\varphi }_{{}_{-}}|)y},$$

$$\lambda >\text{max}\{{\lambda }_{1+},{\lambda }_{1-}\},{\lambda }_{1\pm }=\frac{2k_{\pm }}{|{\varphi }_{\pm }|}.$$

$$A(y)=-\delta \exp \left(-2e^{-\alpha y/\delta }\right)+\delta \text{Ei}\left(-2e^{-\alpha y/\delta }\right)+\beta y,$$

$$\text{Ei}(u)=-{\int }_{-u}^{\infty }{e}^{-t}/tdt.$$

$$\varphi (y)={\varphi }_0\exp \left(-e^{-\alpha y/\delta }\right)-ϵ$$

$$V(\varphi )=18{\left[\alpha \frac{(\varphi +ϵ)}{{\varphi }_0^2}\ln \frac{(\varphi +ϵ{)}^2}{{\varphi }_0^2}\right]}^2-6{\left[\alpha \frac{(\varphi +ϵ{)}^2}{{\varphi }_0^2}\left(1-\ln \frac{(\varphi +ϵ{)}^2}{{\varphi }_0^2}\right)\right]}^2$$

$$+6\beta \left[2\alpha \frac{(\varphi +ϵ{)}^2}{{\varphi }_0^2}\left(1-\ln \frac{(\varphi +ϵ{)}^2}{{\varphi }_0^2}\right)-\beta \right].$$

$${\Lambda }_{-}=-6{\beta }^2,{\Lambda }_{+}=-6(\alpha -\beta {)}^2,0<\beta /\alpha <1.$$

$${\psi }_{\text{L}}\sim \exp \lambda \left[\frac{\delta }{\alpha }{\varphi }_0\text{Ei}\left(-e^{-\alpha y/\delta }\right)+ϵy\right]$$

$$\times \exp \left[2\delta \exp (-2e^{-\alpha y/\delta })-2\delta \text{Ei}\left(-2e^{-\alpha y/\delta }\right)-2\beta y\right].$$

$${\psi }_{\text{L}}(y)\sim \exp \left[-\frac{\alpha }{2e\delta }(\lambda {\varphi }_0-8\alpha ){\left(y-y_0\right)}^2\right]+\mathcal{O}(\alpha y/\delta {)}^3$$

$$y_0\simeq \frac{e\delta \lambda }{2\alpha }\left(\frac{ϵ-{\varphi }_0/e}{\lambda {\varphi }_0-8\alpha }\right)\ll 1.$$

Full text

Keywords: Massless fermions, fermion localization, Yukawa coupling, migratory effect.

1 Introduction

A brane world corresponds to a 4-dimensional hypersurface embedded in a high dimensional bulk, on which gravity is localized [1, 2]. In the most general case, the 4-dimensional sector can be generated by a scalar domain wall interpolating between two asymptotic Anti de-Sitter (AdS) vacuums, with different cosmological constants at each side of the wall, i.e., a solution to the coupled Einstein-Klein Gordon system with a potential with a spontaneously broken discrete symmetry [3, 4, 5, 7, 6, 8, 9, 10, 11]. The gravity fluctuations of this spacetime depend on the wall’s thickness, and are characterized by both a zero mode localized on the wall, and a spectrum of massive states propagating freely. In the thin wall limit [6, 12], the gravitational potential can be calculated, and Newtonian gravity is recovered on the 4-dimensional sector of the scenario.

To localize the matter field, a Yukawa coupling between the bulk fermions and the scalar field of the wall is required; however, in this case, only one chiral massless state (left or right) is localized, while the other is repelled by the gravitational field of the wall [13, 14, 18, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26]. Remarkably, we determined that the minimum bound of the Yukawa constant, obtained from conditions of normalization, is not sufficient to confine the fermion in 4-dimensions with maximum probability on the scalar wall.

On the other hand, it is well known that in absence of $Z_{2}$ symmetry, the gravitational field of the brane shifts the fermion toward the 5-dimensional sector with greater curvature [22]. Here, we show that another shift-like effect exists, regardless of the reflection symmetry of the scenario. For some values of the Yukawa coupling constant, the chiral fermion can be shifted simultaneously at each side of the wall, i.e., the bound state exhibits two relative maximums close to the wall and hence the probability of finding the particle is greater around the wall rather than on it.

To rule out the relationship between the symmetry of the scenario and the fermion’s migration, we consider two models differentiated by their reflection symmetry. The first one corresponds to the member with $Z_{2}$ symmetry of a family of asymmetric domain walls [4], while the second one is an intrinsically asymmetric wall family [9, 21]. In both cases, when the fermion satisfies the normalization condition and it is distributed (centered) around the wall, two peaks at each side of the wall in the probability profile of the fermion are obtained. For the other cases, where only the normalization condition is satisfied, the fermion is always confined outside of the wall.

2 Set-up

Consider the coupled Einstein scalar field system in 5-dimensions (latin index $a,b=0,\cdots,4$; greek index $\mu,\nu=0,\cdots,3$)

[TABLE]

where the metric is given by

[TABLE]

which corresponds to a spacetime with plane parallel symmetry whose energy density is determined by

[TABLE]

Here and in the other sections prime denotes the derivative with respect to $y$.

We are interested in AdS5 domain wall solutions, scenarios for which geometry is determined by a scalar field $\phi$ interpolating between the negative minima of the potential $V(\phi)$. It is well known that the full spectrum of the gravity fluctuations of these scenarios is characterized by a zero mode localized on the wall, $\psi_{\mathrm{g}}\sim e^{2A(y)}$, and a continuum of massive modes propagating freely through the whole 5-dimensional bulk; see [9] for details.

The chiral massless fermions can also be localized on the wall by a Yukawa coupling between the bulk fermions and the scalar field, $\lambda\bar{\Psi}\phi\Psi$ [21]. To evidence this, the Dirac equation for the 5-dimensional spinor in the background (2) must be considered

[TABLE]

To obtain the coordinates representation of the motion equation, it is common to rewrite $\Psi(x,y)$ in terms of its chiral components left (L) and right (R), and factor its spatial degrees of freedom, as shown below

[TABLE]

where $\Psi_{\text{R}}^{\text{L}}(x)\equiv\pm\gamma^{5}\Psi_{\text{R}}^{\text{L}}(x)$, which satisfies the massless 4-dimensional Dirac equation

[TABLE]

Thus, for $\psi_{\text{R}}^{\text{L}}(y)$ we have

[TABLE]

for which

[TABLE]

For a domain wall solution, $e^{A}$ is an integrable and asymptotically vanishing function, such that, $e^{-A}\sim e^{k|y|}$ as $|ky|\gg 1$, with $k=\sqrt{|\Lambda|/6}$, i.e., gravitation produces a repulsive effect on fermions. Hence, for the appropriate values of $\lambda$ one of the $\psi_{\text{R}}^{\text{L}}$ can be normalized. In particular, for $\psi_{\text{L}}$ in the large $y$ limit, it is possible to see that

[TABLE]

where $\phi(y=\pm\infty)=\pm\phi_{0}$ has been considered. Therefore

[TABLE]

in order to cancel the repulsive effect generated by gravity.

The previous restriction for the Yukawa coupling is a necessary condition, but not sufficient to ensure the confinement of $\psi_{\text{L}}$ on the wall. For some values of $\lambda$, we found solutions for $\psi_{\text{L}}$ with two peaks at each sides of the wall, generated by the gravitational competition between the geometry of the spacetime and the Yukawa interaction.

Let $\lambda_{1}$ be the Yukawa critical constant from which (8) is normalizable. Moreover, let $\lambda_{2}$ be the Yukawa threshold from which the maximum of (8) is on the wall ($\lambda_{2}>\lambda_{1}$). Then, for $\lambda<\lambda_{1}$, the fermion is expelled from the thick brane by the gravitational field of the wall; while for $\lambda>\lambda_{2}$ the fermion is localized on the wall. On the other hand, the sector $\lambda_{1}<\lambda<\lambda_{2}$ corresponds to a bounded probability distribution with three critical points: a pair of maximums around the wall and a minimum on it. So, the chiral zero mode of the matter fields is repelled by the thick brane and keeps orbiting around it.

The critical bound $\lambda_{2}$ can be determined by a saddle point evaluation of (8). By assuming that $\psi_{\text{L}}$ is peaked near the brane, $y=0$, and that $\psi_{\text{L}}^{\prime}(0)=0$ and $\psi_{\text{L}}^{\prime\prime}(0)<0$, the following is obtained

[TABLE]

Therefore, to preserve it as a real amount, the following is required

[TABLE]

Now, to determinate the hierarchy between $\lambda_{1}$ and $\lambda_{2}$ we will consider two specific solutions to the Einstein-Klein Gordon system (1), differentiated by the reflection symmetry of each scenario.

3 Symmetric walls

A solution to the coupled system (1) is given by

[TABLE]

[TABLE]

and

[TABLE]

This solution was obtained in [4] and corresponds to a two-parameter domain wall interpolating asymptotically between two AdS5 spacetime, with a cosmological constant given by $\Lambda=-6\alpha^{2}$. Also, it is well known that the zero mode of gravity fluctuations is normalizable, and hence standard gravitation can be recovered in 4-dimensions.

With regard to massless fermions, from (8)

[TABLE]

where $\operatorname{Li_{n}}$ is the polylogarithm function of order $n$ given by

[TABLE]

and from (10) and (12) we find

[TABLE]

and

[TABLE]

respectively.

Fig.1 shows the fermions profile for some values of $\lambda$. We observed that for $\lambda>\lambda_{2}$, $\psi_{\text{L}}$ is normalizable, with maximum probability on the wall. On the other hand, if $\lambda_{1}<\lambda<\lambda_{2}$, a change occurs in the behavior of $\psi_{\text{L}}$; in this case, the competition between gravitational interaction and Yukawa coupling divides the maximum of profile into two maximums around the wall; i.e., the fermions prefer to orbit around the structure.

4 Asymmetric walls

Before considering a thick domain wall with different cosmological constants at each side of the wall, $\Lambda_{+}$ and $\Lambda_{-}$, it is convenient to estimate the asymptotic behavior of $\psi_{\text{L}}$ in this scenario, i.e.,

[TABLE]

where $\phi(y=\pm\infty)=\phi_{\pm}$ and $k_{\pm}=\sqrt{|\Lambda_{\pm}|/6}$ have been considered.

In order to obtain a normalizable solution for $\psi_{\text{L}}$, from (20) reads

[TABLE]

Thus, two critical Yukawa constants associated to $\Lambda_{\pm}$ are obtained, which are the result of the lack of $Z_{2}$ symmetry.

Now, suppose that $\lambda_{1-}>\lambda_{1+}$ and that $\lambda_{1-}>\lambda_{2}>\lambda_{1+}$. In this case, although $\lambda_{2}$ is greater than one of the asymptotic constants, the migratory effect may not be present because $\lambda_{2}<\lambda_{1-}$. On the other hand, there is the possibility that $\lambda_{2}>\lambda_{1-}$ and even though the fermion is normalizable, it is more likely to orbit around the wall. To analyze this, let us consider one specific solutions for the domain wall.

An asymmetric solution is determined by

[TABLE]

where $\alpha$ and $\delta$ are positive constants and Ei is the exponential integral given by

[TABLE]

The scalar field and potential for this scenario are

[TABLE]

and

[TABLE]

This solution was reported in [9] and it corresponds to a two-parameter family of plane symmetric static domain walls without $Z_{2}$ symmetry. The scalar field interpolates between two asymptotic vacua, $\phi_{-}=-\epsilon$ for $y<0$ and $\phi_{+}=\phi_{0}-\epsilon$ for $y>0$, with different cosmological constants

[TABLE]

Now, equation (8) provides that

[TABLE]

Notice that, constants $\epsilon$ and $\beta$ fix the asymptotic values of the scalar field and the cosmological constants, in such a way that, neither for the fermion nor for the graviton it is possible to find normalizable solutions as $\epsilon=\beta=0$, because $\phi_{-}=\Lambda_{-}=0$.

Now, to simplify we will consider the case $\beta=3\alpha/e^{2}$, where the wall is centered in $y=0$, without loss of generality. However, the location of the fermion is shifted with respect to the brane. To demonstrate this, let us calculate $\psi_{\text{L}}$ near $y=0$,

[TABLE]

where

[TABLE]

Thus, the wave function of the fermion is centered at $y_{0}$ with a width given by $\sqrt{e\delta}/\sqrt{(\lambda\phi_{0}-8\alpha)\alpha}$. The fermion’s shift was verified numerically beyond the approximation $\alpha y/\delta\ll 1$, and the results are shown in Fig.2.

Unlike in this article, the shift has previously been used to locate fermions on the wall: in [27] where the spacetime is flat with an orbifold geometry $\mathbb{R}^{4}\times\text{S}^{1}$; and in [22] where the 5-dimensional spacetime is generated by a static double domain wall, in order to have the graviton and the fermion in different sub-walls.

For a scenario generated by a single wall with one non-compact extra dimension, shifting the fermion does not favor the creation of a brane world, since the maximum of probability is not on the wall. Therefore, in agreement with (29), $\epsilon=\phi_{0}/e$ is required to center the fermion on the wall and under this circumstance it is possible to evaluate the migratory effect.

From (21) and (12), the following is obtained

[TABLE]

and

[TABLE]

in such a way that $\lambda_{2}>\lambda_{1-}>\lambda_{1+}$.

In Fig. 3, we plot the zero mode fermion for some values of $\lambda$. For the migratory sector, $\lambda_{1+}<\lambda<\lambda_{2}$, we can see two interfaces associated to the bounded state of $\psi_{\text{L}}$. The difference in the amplitude of the peaks is generated by the gravitational repulsion applied on the fermions, which expels them towards the region of greater curvature. For $\lambda>\lambda_{2}$ the probability’s maximum is in $y=0$, and therefore, $\psi_{\text{L}}$ is confined on the wall.

5 Discussion

The Yukawa coupling counteracts the repulsion exerted on the fermions by the gravitation of the wall, such that, utilizing a critical value of the Yukawa constant $\lambda>\lambda_{1}$ (obtained by asymptotic analysis of the wave function), ensures the normalization of one of the chiral states. However, we found a range of values for the Yukawa constant, $\lambda_{1}<\lambda<\lambda_{2}$ (by the Laplace method), where the fermion is in a normalizable state orbiting around the wall.

We calculated the Yukawa bounds, $\lambda_{1}$ and $\lambda_{2}$, in two scenarios generated by a thick domain wall with or without reflection symmetry.

On the domain wall solution with reflection symmetry [4], $A(y)$ and $\phi(y)$, the chiral massless fermion $\psi_{\text{L}}$ is normalizable around the wall. On the other hand, the symmetric wall is part of a family of asymmetric solutions where $A(y)\rightarrow A(y)+\beta y$ and $\phi(y)\rightarrow\phi(y)-\epsilon$ [22], and where the zero mode of the fermion is determined by

[TABLE]

the fermion is shifted from the wall by the term $e^{-2\beta y}$. Therefore, it is not possible to evaluate the migratory effect on another member of the family of solutions, different than the $Z_{2}$-limit case: $(\beta,\epsilon)\rightarrow(0,0)$.

To evaluate the migratory effect on a scenario without $Z_{2}$ symmetry, an intrinsically asymmetric domain wall [9], $A(y)$ and $\phi(y)$, is required. This scenario is also a member of a two-parameter family, $A(y)+\beta y$ and $\phi(y)-\epsilon$, but without the $Z_{2}$-limit, and to ensure that $\psi_{\text{L}}(y)$ is centered in the maximum of $\rho(y)$, $\beta=3\alpha/e^{2}$ and $\epsilon=\phi_{0}/e$ must be utilized.

In both cases the constraints imposed on the fermion in order to get it on the wall, lead to a hierarchy between $\lambda_{2}$ and $\lambda_{1}$, such that, regardless of the scenario’s symmetry, the migratory effect could be present.

6 Acknowledgements

This work was supported by IDI-ESPOCH under the project titled Domain walls and their gravitational effects.