Localization and mass spectra of various matter fields on Weyl thin brane

TL;DR

This paper investigates how various matter fields behave on a Weyl brane, showing zero modes can be localized and analyzing the mass spectra and gravitational corrections in this five-dimensional model.

Contribution

It provides a detailed analysis of matter field localization and mass spectra on Weyl branes, including gravitational potential corrections, expanding understanding of brane-world models.

Findings

Zero modes of matter fields can be localized on the positive tension brane.

Mass spectra are equidistant for higher states and sparse for lower states.

Newtonian potential correction is proportional to 1/r^3.

Abstract

It has been shown that the thin brane model in a five-dimensional Weyl gravity can deal with the wrong-signed Friedmann-like equation in the Randall--Sundrum-1 (RS1) model. In the Weyl brane model, there are also two branes with opposite brane tensions, but the four-dimensional graviton (the gravity zero mode) is localized near the negative tension brane, while our four-dimensional universe is localized on the positive tension brane. In this paper, we consider the mass spectra of various bulk matter fields (i.e., scalar, vector, and fermion fields) on the Weyl brane. It is shown that the zero modes of those matter fields can be localized on the positive tension brane under some conditions. The mass spectra of the bulk matter fields are equidistant for the higher excited states, and relatively sparse for the lower excited states. The size of the extra dimension determines the gap of the mass spectra. We also consider the correction to the Newtonian potential in this model and it is proportional to $1/r^{3}$.

Tables

Table 1. Table 1: Localization conditions for the zero modes of various bulk matter fields. Here z b subscript 𝑧 𝑏 z_{b} is the size of the extra dimension.

Bulk matter	Lagrangian	$z_b$	condition
Scalar	${ℒ}_0=-\frac{1}{2}{e}^{\lambda \omega }{\partial }_M\Phi {\partial }^M\Phi$	finite	$\lambda >k+\alpha$
		infinite	$\lambda >k+2\alpha$
Vector	${ℒ}_1=-\frac{1}{4}{e}^{\tau \omega }{F}_{MN}{F}^{MN}$	finite	$\tau >\frac{1}{3}\left(k+\alpha \right)$
		infinite	$\tau >\frac{1}{3}\left(k+4\alpha \right)$
Fermion	${ℒ}_{1/2}=\overline{\Psi }{\Gamma }^M\left({\partial }_M+{\omega }_M\right)\Psi -\eta \overline{\Psi }F\left(\omega \right)\Psi$	finite
		infinite

Equations

$$S_5=\frac{M_{\ast }^3}{2}{\int }_{M_5^W}{d}^{5}x\sqrt{-g}{e}^{k\omega }\left[R-\left(4+4k\right){\nabla }_K\omega {\nabla }^K\omega \right].$$

$$R_{LKN}^M$$

$$R_{LN}$$

$$R$$

$${\Gamma }_{MN}^P={\{}_{MN}^P\}-\frac{1}{2}\left({\delta }_N^P{\nabla }_M\omega +{\delta }_M^P{\nabla }_N\omega -g_{MN}{\nabla }_P\omega \right),$$

$$R_{MN}$$

$${\hat{G}}_{MN}$$

$$k\hat{R}$$

$${\hat{G}}_{MN}={\hat{T}}_{MN}.$$

$$d{s}^2$$

$$k{\omega }^{\prime \prime }+(k^2+2k+2){\omega }^{\prime 2}$$

$$+2k\frac{a^{\prime }}{a}{\omega }^{\prime }+3\frac{a^{\prime \prime }}{a}$$

$$(1+k){\omega }^{\prime 2}-2k\frac{a^{\prime }}{a}{\omega }^{\prime }-3\frac{a^{\prime 2}}{a^2}$$

$$2(1+k)({\omega }^{\prime \prime }+3\frac{a^{\prime }}{a}{\omega }^{\prime })+k(1+k){\omega }^{\prime 2}$$

$$-k(\frac{a^{\prime 2}}{a^2}+2\frac{a^{\prime \prime }}{a})$$

$$a(z)$$

$$\omega (z)$$

$${\hat{T}}_{MN}=\frac{2}{\sqrt{-g}}\frac{\delta \hat{S}}{\delta g_{MN}}={\hat{R}}_{MN}-\frac{1}{2}{g}_{MN}\hat{R},$$

$$\rho (z)$$

$$S_{\text{s}}=-\frac{1}{2}\int d^{5}x\sqrt{-g}{e}^{\lambda \omega }{g}^{MN}{\partial }_M\Phi {\partial }_N\Phi ,$$

$$d{s}^2$$

$${\partial }_P\left(\sqrt{-g}{\partial }^P\Phi \right)=0.$$

$$\Phi (x^{\mu },z)=\sum _n{\varphi }_n(x^{\mu }){\chi }_n(z)e^{-\varrho A(z)},$$

$$\left[-{\partial }_z^2+V_{\text{s}}(z)\right]{\chi }_n(z)=m_n^2{\chi }_n(z),$$

$$V_{\text{s}}(z)=\varrho {\partial }_z^{2}A(z)+[\varrho {\partial }_{z}A(z){]}^2.$$

$${\int }_{-z_b}^{+z_b}{\chi }_m(z){\chi }_n(z)dz={\delta }_{mn},$$

$$S_{\text{s}}=-\frac{1}{2}\sum _n\int d^{4}x({\eta }^{\mu \nu }{\partial }_{\mu }{\varphi }_n{\partial }_{\nu }{\varphi }_n+m_n^2{\varphi }_n^2),$$

$$V_{\text{s}}(z)$$

$$V_{\text{s}}(0)$$

$$\lambda >k+\alpha .$$

$${\chi }_0(z)=\sqrt{\frac{\beta (k+2\alpha -\lambda )}{2\alpha \left[(1+\beta z_b{)}^{\frac{(k+2\alpha -\lambda )}{\alpha }}-1\right]}}{\left(1+\beta |z|\right)}^{\frac{k+\alpha -\lambda }{2\alpha }},$$

$${\chi }_0(z)=\sqrt{\frac{\beta }{2\ln (1+\beta z_b)}}\frac{1}{\sqrt{1+\beta |z|}}(\lambda =k+2\alpha ),$$

$$\overline{|z|}\equiv \beta |z|,{\bar{m}}_n\equiv \frac{m_n}{\beta },{\bar{V}}_s(\bar{z})\equiv \frac{V_s(z)}{{\beta }^2},$$

$$(-{\partial }_{\bar{z}}^2+{\bar{V}}_s)\chi (\bar{z})={\bar{m}}_n^2\chi (\bar{z}),$$

Full text

I Introduction

Over the last decades, there was a great interest to construct a unified theory in higher-dimensional spacetime. The higher-dimensional theory dates back to the work of Kaluza and Klein KK26 in the 1920s. They tried to unify electromagnetism with gravity by assuming that the $g_{5\mu}$ component of the five-dimensional metric represents the electromagnetic potential KK26 . Subsequently, the brane world theories appear as an alternative theory to solve the gauge hierarchy and cosmological constant problems (see e.g., Refs. Rubakov ; RSmodel1 ; Randall2 ; Arkani-Hamed1998 ) with a fundamental idea that the visible universe is localized on a 3-brane embedded in a higher-dimensional bulk.

Most research as regards about brane worlds is based on a Riemann geometry Bajc2000 ; Fuchune ; Hatanaka1999 ; Das2008 ; Volkas2007 ; RandjbarPLB2000 ; Guo1103 ; Liu2009uca ; Shiromizu2000 . The non-Riemannian models can be classified within the language of metric-affine gravity Liu0708 ; Puetzfeld2005 ; Hehl1995 ; Hammond2002 . Generally, there are two types of generalized Riemann geometries, Riemann–Cartan geometry and Weyl geometry Puetzfeld2005 . Cartan introduced the torsion $T$ to represent the antisymmetric piece of the connection, which may be related to spin. When the connection includes a geometric scalar $\omega$, there will be a nonmetricity part called $Q$, which is related to the scalar $\omega$. In Riemann–Cartan geometry, the torsion $T\neq 0$, and the nonmetricity part $Q=0$. When $T=0$, the Riemann–Cartan geometry can degenerate into Riemann geometry. The Weyl geometry is an affine manifold described by a pair ($g_{MN}$, $\omega_{M}$), where $g_{MN}$ is the metric and $\omega_{M}$ is the “gauge” vector as the gradient of a scalar function $\omega$. In the Weyl geometry, $T=0$ and $Q\neq 0$. When $Q=0$, the Weyl geometry can also degenerate into Riemann geometry. As a generalization of Riemann, Weyl geometry has been received much attention in the studies in field theory and cosmology Barbosa2006 ; ChengPRL1998 ; Hochberg1991 ; MoonJCAP2010 . Moreover, in Refs. Barbosa2008 ; LiuJHEP2008 ; Correa1607 ; Arias0202 ; Cendejas0603 ; Liu0708 ; Liu0911 , the authors have studied the Weyl brane world scenario and obtained various solutions of Minkowski thick branes.

In the brane world scenario with all matter fields (scalar, vector, and fermion fields) propagating in the bulk, the zero modes of these bulk fields, i.e., the four-dimensional matter fields in the Standard Model, should be localized on the brane complying with the present observations. So, localization of various bulk matter fields on a brane is an important and interesting issue. References Bajc2000 ; Fuchune ; Furlong1407 showed that the zero mode of a free massless scalar field can be localized on the RS brane or its generalized branes. In order to ensure that the zero mode of vector field can be localized on the brane, the case of generalized RS branes with more extra dimensions was considered Oda2000 ; Sakamura1607 , a dynamical mass term was included in gauge field localization zhao1406 ; alencer1409 , and in Refs. Araujo1406 ; Costa1304 ; Aguilar1401 ; Costa1501 ; sui1703 ; Chumbes1108 ; Cruz1211 ; Zhao1402 ; Alencar1005 ; Alencar1008 the authors introduced the interaction with the background scalar field. Localization of the fermion zero mode on a brane is especially important. It has been demonstrated that localization mechanisms, e.g. Yukawa coupling $\eta\bar{\Psi}F(\phi)\Psi$ Liu0708 ; Volkas2007 ; RandjbarPLB2000 ; Guo1103 ; Fuentevilla1412 ; Cendejas1503 ; Rocha1606 ; Brito1609 ; KoleyCQG2005 ; D. Bazeia0809 ; MelfoPRD2006 ; Almeida0901 ; Almeida2009 ; Davies0705 ; LBCastro1 ; LBCastro2 , derivative scalar Cfermion coupling $\eta\bar{\Psi}\Gamma^{M}\partial_{M}F(\phi)\gamma^{5}\Psi$ Liu1312 , and derivative geometrical–fermion coupling

$\eta\bar{\Psi}\Gamma^{M}\partial_{M}F(R)\gamma^{5}\Psi$ Liu1701 , can ensure that the fermion zero mode can be localized on the brane.

In the RS1 model RSmodel1 , the authors assumed that there are two thin branes located at the boundaries of the extra dimension. One of them is the negative tension brane (called Tev brane, where our universe is resided), and the other is the positive tension brane (called the Planck brane, where the spin-2 gravitons are localized). The warping of the extra dimension in RS1 model can solve the famous gauge hierarchy problem, but the fundamental assumption that the our universe is localized on the negative tension brane can lead to a wrong-signed Friedmann-like equation Csaki99 ; Cline99 ; Shiromizu99 . The authors of Ref. Yang:2011pd presented a generalized RS1 model in the Weyl integrable geometry, which can give a correct Friedmann-like equation since our universe is assumed to be located on the positive tension brane. It was shown that the brane system is stable under the tensor perturbation and the gravity zero mode is localized near the negative tension brane, whereas the matter fields should be confined on the positive tension brane, since our universe is located on this brane.

In this paper, we will do further research on the localization of various bulk matter fields in this Weyl brane model and check whether the zero modes of these fields can be confined on the positive tension brane. We will consider a dilaton coupling between the bulk matter fields and the background scalar field. It will be found that the coupling is necessary for most of the bulk fields. For the fermion case, we consider the usual Yukawa coupling $\eta\bar{\Psi}F(\omega)\Psi$ to obtain analytical wave functions of the fermion KK modes. This paper is organized as follows. In Sect. II, we give a brief review of the Weyl brane model proposed in Ref. Yang:2011pd . In Sect. III, we investigate localization and mass spectra of the scalar, vector, and fermion fields on the Weyl brane. In Sect. IV, gravitational fluctuations and correction to the Newtonian potential are discussed. Finally, we give a brief conclusion and summary in Sect. V.

II Review of Weyl brane model

The brane model we consider here is based on the five-dimensional Weyl geometric action where the gravity is non-minimally coupled with a background scalar field $\omega$. The action is given by Yang:2011pd

[TABLE]

Here $R$ is the five-dimensional Weyl Ricci scalar, $M^{-3}_{*}\equiv 8\pi G_{5}$ is the fundamental scale of gravity, the Latin letters $M,N,K=0,1,2,3,5$, and $\omega_{,K}$ is a gradient of the Weyl scalar $\omega$. In what follows it is convenient to denote the physical four coordinates by $x^{\mu}$, and the extra dimension by $y$.

In this frame the Weyl Ricci tensors are constructed from $\Gamma_{LN}^{M}$ in a standard way,

[TABLE]

where the Weyl affine connection is expressed as

[TABLE]

with $\{^{P}_{MN}\}$ representing the Riemannian Christoffel symbol. Now with some algebra, the Weyl Ricci tensor $R_{MN}$ can be separated into a Riemannian Ricci tensor part $\hat{R}_{MN}$ plus some other terms with respect to the Weyl scalar function, i.e.,

[TABLE]

where the hatted magnitudes and operators are defined by the Christoffel symbol. Then from the action (1), the equations of motion are

[TABLE]

Further, we can write the field equation (5) in the form of the Einstein equation by moving all Weyl scalar terms to the right-hand side to compose an effective energy-momentum tensor:

[TABLE]

The authors in Refs. Barbosa2006 ; Barbosa2005 ; Barbosa2008 ; Vladimir0904 ; Vladimir0912 ; D. Bazeia1502 ; Marzieh1504 studied thick brane solutions, by considering the metric that satisfies four-dimensional Poincar$\acute{e}$ invariance. The localization and mass spectrum problems of bulk matter fields on the Weyl Minkowski thick branes were discussed in Refs. Barbosa2005 ; Barbosa2008 ; LiuJHEP2008 ; Correa1607 ; Arias0202 ; Cendejas0603 ; Liu0708 ; Liu0911 . We are interested in the Weyl Minkowski thin branes with space $S^{1}/Z_{2}$, for which the metric with four-dimensional Lorentz invariance in a five-dimensional Weyl spacetime is given by

[TABLE]

where $z$ is the conformal coordinate of the extra dimension with $z\in[-z_{b},z_{b}]$ and it relates to the physical coordinate $y$ by a coordinate transformation $dy=a(z)dz$. One assumes that the background scalar field $\omega$ only depends on the extra dimension.

According to the metric (8), the field equations in the bulk can be expressed in terms of the warp factor $a$, and the scalar field $\omega$,as in the following set of equations:

[TABLE]

where the prime denotes the derivative with respect to $z$. The above field equations admit the following brane solution Yang:2011pd :

[TABLE]

where $\alpha=\sqrt{k^{2}+3k+3}$, the parameter $\beta>0$ and $k<-1$.

The five-dimensional effective energy-momentum tensor is given as usual by

[TABLE]

with $\hat{S}$ the effective matter action. In this case the effective energy density is given by the null-null component $\hat{T}_{0}^{0}$ of the energy-momentum tensor,

[TABLE]

This expression directly shows that there are two delta functions appearing at each boundary, so they stand for two thin branes, namely a positive tension brane at the origin and a negative one at the boundary $z_{b}$. These tensions compensate for the effects produced by the bulk component and hence ensure the existence of four-dimensional flat branes. So the brane configuration of this case is similar to that of the RS1 model. However, localization of massless graviton in these two brane scenarios will be quite different, and in Ref. Yang:2011pd one can suppose that the standard model fields are confined on the positive tension brane at $z=0$, so this is crucial to overcome the severe cosmological problem of the RS1 model.

III Localization of various matter fields

In this section, we investigate the localization and mass spectra of various bulk matter fields, including the spin-0 scalar, spin-1 vector and spin-1/2 fermion fields on the Weyl brane. A natural mechanism for localizing spin-0 scalar and spin-1 vector fields on the brane involves the couplings with the background scalar field. For the case of spin-1/2 fermions, we will investigate the usual Yukawa coupling $\eta\bar{\Psi}F(\omega)\Psi$.

III.1 Spin-0 scalar field

Firstly, we explore the localization condition and the mass spectrum for a massless scalar field on the Weyl brane by considering two types of coupling with the background scalar fields, dilaton coupling and Higgs potential coupling Liu2009uca .

III.1.1 The dilaton coupling

Let us consider a massless scalar field $\Phi$ coupled with the dilaton field $\omega$ via the following five-dimensional action:

[TABLE]

where $\lambda$ is a dimensionless coupling constant. For convenience, we rewrite the warp factor so that the metric (8) takes the form

[TABLE]

Following (13) and (14), the equation of motion for the scalar field satisfies the five-dimensional Klein CGordon equation,

[TABLE]

Making the KK decomposition

[TABLE]

with the quantity $\varrho=\frac{3}{2}\frac{k+\alpha-\lambda}{k+\alpha}$, and demanding that $\phi_{n}(x^{\mu})$ satisfies the four-dimensional massive Klein–Gordon equation $(\partial_{\mu}\partial^{\mu}-m^{2}_{n})\phi_{n}=0$, we obtain the equation of the extra-dimensional part $\chi_{n}(z)$ of the scalar KK mode, which can be converted into the following Schrödinger-like equation:

[TABLE]

where $m_{n}$ is the mass of the scalar KK mode $\phi_{n}(x^{\mu})$ and the effective potential $V_{\text{s}}(z)$ reads

[TABLE]

It is clear that the effective potential $V_{s}$ is only dependent on the warp factor. In fact, the Schrödinger-like equation (17) can be written as $\mathcal{H}\chi_{n}=m_{n}^{2}\chi_{n}$, where the Hamiltonian operator is given by $\mathcal{H}=Q^{+}Q$ with $Q=-\partial_{z}+\varrho\partial_{z}A$. Since the operator $\mathcal{H}$ is positive definite, there are no KK modes with negative $m_{n}^{2}$.

By substituting the KK decomposition (16) into the five-dimensional scalar action (13) and introducing the following orthonormality conditions:

[TABLE]

we get the effective low-energy theory in a four-dimensional familiar form

[TABLE]

where we have used Eq. (17).

Due to the explicit forms of $\varrho$ and $A(z)$, the expression of the effective potential $V_{\text{s}}(z)$ corresponding to Eq. (18) is given by

[TABLE]

where $s=\frac{\alpha+k-\lambda}{\alpha}$. The value of $V_{\text{s}}(z)$ at $z=0$ is

[TABLE]

In order to localize the scalar zero mode on the positive tension brane, the effective potential $V_{\text{s}}(z)$ should be negative at $z=0$, which results in the following condition:

[TABLE]

By setting $m=0$ in Eq. (17), we easily get the normalized zero mode of the scalar field for $\lambda\neq k+2\alpha$:

[TABLE]

which is indeed localized on the positive tension brane under the condition (23). When the extra dimension is infinite, the normalization condition $\int^{+\infty}_{-\infty}\chi_{0}(z)^{2}dz=1$ means that the localization condition for the scalar zero mode is much stronger: $\lambda>k+2\alpha$ . For $\lambda=k+2\alpha$ , the normalized scalar zero mode reads as

[TABLE]

and it is also localized on the positive tension brane. But if the extra dimension is infinite, the scalar zero mode is not normalized and hence we cannot get a localized scalar zero mode. When there is no coupling between the scalar and dilaton fields, i.e., the parameter $\lambda=0$, one can verify that the scalar zero mode is also localized on the positive tension brane when $z_{b}$ is finite.

In order to investigate the mass spectrum of the scalar KK modes, we redefine the following dimensionless parameters:

[TABLE]

with which the Schr$\ddot{o}$dinger-like equation (17) is rewritten as

[TABLE]

where the effective potential $\bar{V}_{s}$ reads as

[TABLE]

The general solution of Eq. (27) is given in terms of the combination of Bessel functions

[TABLE]

where $N_{n}$ is the normalization coefficient, $C_{1}$ and $C_{2}$ are the $m$-dependent parameters, $\text{J}_{P}(z)$ and $\text{Y}_{P}(z)$ are the Bessel functions of the first and second kinds.

Then we impose the two kinds of boundary condition (the Neumann boundary condition and the Dirichlet boundary condition) to calculate the mass spectra of scalar KK modes, respectively. For the the Neumann boundary condition $\partial_{\bar{z}}(e^{-\varrho A(\bar{z})}\chi_{n})|_{z=0,z=z_{b}}=0$, we get the mass spectrum which is determined by following condition:

[TABLE]

where $P_{\text{Ns}}$ and $\mathcal{C}_{\text{Ns}}$ are defined as:

[TABLE]

And, for the Dirichlet boundary condition $e^{-\varrho A(\bar{z})}\chi_{n}=0$ at the boundaries $\bar{z}=0$ and $\bar{z}=\bar{z}_{b}$, we get the mass spectrum of the scalar KK modes which satisfies the following relation:

[TABLE]

where $P_{\text{Ds}}$ and $\mathcal{C}_{\text{Ds}}$ are given by

[TABLE]

From Eqs. (30) and (32), we see that the Neumann and Dirichlet boundary conditions give the similar expressions which determine the mass spectrum, whereas the orders of the Bessel functions are different as shown in (30) and (31). By tuning the parameters $k$ and $\lambda$, the same orders of Bessel functions can be obtained, and thus the expression of the mass spectrum derived by the Neumann boundary condition is equivalent to that derived by the Dirichlet boundary condition. Therefore, we only consider the Neumann boundary condition for later discussion.

According to Eq. (30), we plot the relation between $\mathcal{M}_{N}$ and $\bar{m}_{n}$ in Figs. 1, 2, 3, where the zero points represent the mass of the KK modes. Figure 1 shows that the effect of the parameter $\lambda$ on the mass $\bar{m}_{1}$ of the first massive scalar KK mode. It can be seen that the mass $\bar{m}_{1}$ increases with $\lambda$. On the other hand, there are some singular (vertical) lines in these figures, which come from the zero points of ${\text{Y}_{P_{\text{Ns}}}(\bar{m}_{n})}$ in (31). Since the gaps between two arbitrary adjacent singular lines are almost the same for a set of parameters, we can call the average gap the period $T$. For example, the periods $T$ for the three cases shown in Fig. 1 are 3.13, 3.22 and 3.42, respectively. The period for a special mass spectrum increases with $\lambda$. Figure 2 states the effect of the size of the extra dimension, $\bar{z}_{b}$, on the gap of the massive KK modes. It shows that the number of excited states in a single period increases with $\bar{z}_{b}$, and the gap of massive KK modes decreases with $\bar{z}_{b}$. An interesting phenomenon will appear if we set $\bar{z}_{b}<1$ (and $k=-3$ and $\lambda=5$ at the same time). The massive KK modes will not appear in every period anymore for this case; however, they will appear after severe periods, as shown in Fig. 3. The mass spectra of the scalar KK modes are shown in Fig. 4 for different parameters. The mass gap between the zero mode and the first massive KK mode increases with the dilaton coupling constant $\lambda$ and decreases with the size of the extra dimension $\bar{z}_{b}$. On the other hand, the mass spectrum is almost equidistant for the higher excited states, and relatively sparse for the lower excited states.

III.1.2 The Higgs potential coupling

Next, we consider another coupling, the Higgs potential type coupling between a real scalar $\Phi$ and the background scalar field $\omega$. The five-dimensional action is assumed as Liu2009uca ,

[TABLE]

where the Higgs potential $V(\Phi,\omega)$ is given by

[TABLE]

According to the metric (14) and the action (34), the equation of motion of the five-dimensional real scalar $\Phi$ reads

[TABLE]

Then, making the KK decomposition

$\Phi(x^{\mu},z)=\sum_{n}\phi_{n}(x^{\mu})\chi_{n}(z)e^{-3A/2}$ and demanding that $\phi_{n}(x^{\mu})$ satisfies the four-dimensional massive Klein–Gordon equation $(\partial_{\mu}\partial^{\mu}-m^{2}_{n})\phi_{n}=0$, we obtain the Schr$\ddot{o}$dinger-like equation of the extra–dimensional part $\chi_{n}(z)$

[TABLE]

Here $m_{n}$ is the mass of the $n$-th KK mode, and the effective potential $V_{s}(z)$ is given by

[TABLE]

The fundamental five-dimensional action (34) can be reduced to the effective four-dimensional action of the massive scalars, when $\chi_{n}$ satisfy the following orthonormality conditions:

[TABLE]

The explicit expression of the effective potential $V_{s}(z)$ for the brane solution (10) is

[TABLE]

At $z=0$, the value of $V_{s}(z)$ is

[TABLE]

Considering that $\beta>0$, $k<-1$ and $\sqrt{3}/2<\alpha=\sqrt{k^{2}+3k+3}<-k$, we know that the effective potential $V_{s}(z)$ is negative at $z=0$. Therefore, the scalar zero mode can be localized on the positive tension brane.

For convenience, we redefine some new dimensionless parameters:

[TABLE]

and the Schr$\ddot{o}$dinger-like equation (37) is changed as

[TABLE]

where the expression of the dimensionless effective potential $\bar{V}_{s}$ is

[TABLE]

At the boundary $\bar{z}=\bar{z}_{b}$, the effective potential $\bar{V}_{s}(\bar{z})$ (44) has a positive $\delta$ function (just as an infinite high barrier), and it will yield infinite discrete bound KK modes. Due to the non-analytical solution of the Schr$\ddot{o}$dinger-like equation (43), we can only solve numerically Eq. (43) by fixing the parameter $\bar{\theta}$ and choosing the proper value $\bar{u}$, to ensure the zero mode can be localized on the brane.

The wave functions of the lower bound KK modes are plotted in Fig. 5 and the corresponding mass spectrum is listed as follows:

[TABLE]

where the parameters are set to $k=-2,~{}\bar{\theta}=4.77,~{}\bar{z}_{b}=1$, and $\bar{u}=1.43$. From the numerical solution above, we know that the zero mode can be trapped on the brane by fine tuning the parameters. We can see that the mass spectrum would get denser for the higher excited states while they get sparser for the lower excited states, and the gap of the higher excited states is almost the same for a set of parameters, on account of the dimensionless effective potential $\bar{V}_{s}(\bar{z})$ at the boundary $\bar{z}=\pm\bar{z}_{b}$ having two positive $\delta$ functions, and the effective potential $\bar{V}_{s}(\bar{z})$ just like an infinite high barrier. Therefor, the solutions at the high excited states should be trigonometric functions with the same mass gap. By increasing the value of the parameter $\bar{u}$, the bound state solution with $\bar{m}^{2}_{n}<0$ will appear. Therefore, in order to obtain nonnegative eigenvalues $\bar{m}^{2}$, the mass parameter $\bar{u}$ will have an upper limit when the other parameters are fixed. For $k=-2,~{}\bar{\theta}=4.77,~{}\bar{z}_{b}=1$, as we have chosen here, the upper limit of $\bar{u}$ is $1.43$.

III.2 Spin-1 vector field

Secondly, we turn to the spin-$1$ vector field. A typical five-dimensional action of a vector field coupled with a dilaton field $\omega$ is given as

[TABLE]

where $F_{MN}=\partial_{M}A_{N}-\partial_{N}A_{M}$ is the field strength as usual and $\tau$ is a dimensionless coupling constant. Considering the explicit form of the metric (14), the equations of motion read as

[TABLE]

In order to be consistent with the gauge invariant equation $\oint dzA_{4}=0$, we use the gauge freedom to choose $A_{4}=0$. Under the KK decomposition of the vector field

[TABLE]

where $\epsilon=\frac{k+\alpha-3\tau}{2(k+\alpha)}$, and the orthonormality conditions

[TABLE]

the action (46) can be reduced to an effective one including the four-dimensional massless vector field (the zero mode) and a set of massive vector fields (the massive KK modes):

[TABLE]

where $f^{(n)}_{\mu\nu}=\partial_{\mu}a^{(n)}_{\nu}-\partial_{\nu}a^{(n)}_{\mu}$ is the field strength of the four-dimensional vector field. In addition, the extra–dimensional part $\rho_{n}(z)$ satisfies the following Schrödinger-like equation:

[TABLE]

where $m_{n}$ is the mass of the vector KK mode and the effective potential $V_{\text{v}}(z)$ is

[TABLE]

The explicit expression of the effective potential $V_{\text{v}}(z)$ reads

[TABLE]

where b=$\left(k-5\alpha-3\tau\right)$. The value of the effective potential $V_{\text{v}}(z)$ at $z=0$ is

[TABLE]

For localizing the vector zero mode on the positive tension brane, the effective potential $V_{\text{v}}(z)$ should be negative at $z=0$, which is equivalent to the following condition:

[TABLE]

By setting $m=0$ in Eq. (51), we can get the normalized zero mode of the vector field for $~{}\tau\neq\frac{k+4\alpha}{3}$:

[TABLE]

where $d=k+4\alpha-3\tau$. It is localized on the positive tension brane when the extra dimension is finite under the condition (55). It is easy to see that the vector zero mode can also be localized on the brane when $\tau=0$ if the extra dimension is finite. If the parameter $\tau=\frac{k+4\alpha}{3}$, the normalized vector zero mode is

[TABLE]

The vector zero mode (57) can be localized on the positive tension brane only for the case of a finite extra dimension. However, we note here that if the extra dimension is infinite, the zero mode can also be localized on this brane for $\tau>\frac{k+4\alpha}{3}$.

By defining the dimensionless potential

[TABLE]

Eq. (51) can be rewritten in a dimensionless form,

[TABLE]

where the dimensionless parameters $\bar{z}$ and $\bar{m}_{n}$ are defined as in Eq. (26). For simplicity, we only require that the KK modes satisfy the Neumann boundary condition $\partial_{\bar{z}}(e^{-\epsilon A(\bar{z})}\rho_{n})=0$ at the boundaries $\bar{z}=0$ and $\bar{z}=\bar{z}_{b}$. Using the boundary condition at $\bar{z}=0$, we get the general solution of Eq. (51),

[TABLE]

where $N$ is the normalization coefficient and

[TABLE]

With another boundary condition at $\bar{z}=\bar{z}_{b}$, we obtain the spectrum of the vector KK modes determined by the following condition:

[TABLE]

We plot the relation between $\mathcal{M}$ and $\bar{m}$ In Figs. 6 and 7 for large and small $\bar{z}_{b}$, respectively. Figure 6(a) and Fig. 6(b) imply that the mass of the first excited state increases with the coupling constant of dialton, $\tau$. And Fig. 6(b) and Fig. 6(c) show that in a single cycle the number for the excited states increases with the size of the extra dimension, $\bar{z}_{b}$. In Fig. 7, we adjust the size of the extra dimension $\bar{z}_{b}=0.1$ and $\bar{z}_{b}=0.06$ and let $k=-3,\tau=20$ at the same time; the result shows that the excited states do not exist in each cycle and they only emerge after multiple periods. In addition, Fig. 8 shows that the spectrum interval approaches a constant for the higher excited states, while for the lower excited states it is relatively sparse.

III.3 Spin-1/2 fermion field

Furthermore, we investigate the localization and mass spectrum of spin-1/2 fermions on the Weyl brane. For studying the localization of fermions on thick branes, Refs. Liu0708 ; Volkas2007 ; Fuentevilla1412 ; Cendejas1503 ; Rocha1606 ; Brito1609 ; KoleyCQG2005 ; D. Bazeia0809 ; MelfoPRD2006 ; Almeida0901 ; Almeida2009 ; Davies0705 ; LBCastro1 ; LBCastro2 have introduced the Yukawa coupling. Here, we consider the following five-dimensional Dirac action for a fermion coupled with the background scalar field $\omega$:

[TABLE]

where $\omega_{M}$ is the spin connection defined as

$\omega_{M}=\frac{1}{4}\omega_{M}^{\bar{M}\bar{N}}\Gamma_{\bar{M}}\Gamma_{\bar{N}}$ with $\omega_{M}^{\bar{M}\bar{N}}$ given by

[TABLE]

The non-vanishing components of the spin connection $\omega_{M}$ for the background metric (14) are

[TABLE]

where $\hat{\omega}_{\mu}=\frac{1}{4}\bar{\omega}_{\mu}^{\bar{\mu}\bar{\nu}}\Gamma_{\bar{\mu}}\Gamma_{\bar{\nu}}$ is the spin connection derived from the metric $\hat{g}_{\mu\nu}(x)=\hat{e}_{\mu}^{~{}\bar{\mu}}(x)\hat{e}_{\nu}^{~{}\bar{\nu}}(x)\eta_{\bar{\mu}\bar{\nu}}$. Then we can obtain the equation of motion

[TABLE]

where $\gamma^{\mu}(\partial_{\mu}+\hat{\omega}_{\mu})$ is the Dirac operator on the brane. For the current case of flat brane, the spin connection on the brane vanishes, i.e., $\hat{\omega}_{\mu}=0$.

We make a general chiral decomposition,

[TABLE]

where $\psi_{Ln}(x)=-\gamma^{5}\psi_{Ln}(x)$ and $\psi_{Rn}(x)=\gamma^{5}\psi_{Rn}(x)$ are the left- and right-chiral components of a four-dimensional Dirac field, respectively. Then we can show that $\psi_{L}(x)$ and $\psi_{R}(x)$ satisfy the four-dimensional massive Dirac equations: $\gamma^{\mu}(\partial_{\mu}+\hat{\omega}_{\mu})\psi_{Ln}(x)=m_{n}\psi_{R_{n}}(x)$ and $\gamma^{\mu}(\partial_{\mu}+\hat{\omega}_{\mu})\psi_{Rn}(x)=m_{n}\psi_{L_{n}}(x)$, and the KK modes $f_{Ln}(z)$ and $f_{Rn}(z)$ satisfy the following coupled equations

[TABLE]

From the above coupled equations, we can obtain the Schrödinger-like equations for the left- and right-chiral KK modes of the fermion

[TABLE]

where the effective potentials are given by

[TABLE]

In order to get the effective four-dimensional action for a massless fermion and a set of massive Dirac fermions

[TABLE]

we need to introduce the following orthonormality conditions:

[TABLE]

We set $F(\omega)=\partial_{z}(e^{\upsilon\omega})$ as a simple example. The explicit expressions of the effective potentials read

[TABLE]

One can easily calculate the values of $V_{L}(z)$ and $V_{R}(z)$ at $z=0$:

[TABLE]

It can be shown that the left- and right-chiral fermion zero modes cannot be localized on the positive tension brane at the same time. For localizing the zero mode of the left-chiral fermion on the positive tension brane, the effective potential $V_{L}(z)$ should be negative at $z=0$, which requires that $\eta v<0$, since $\beta$ and $\alpha$ are positive.

In the following, we suppose the left-chiral fermion has the zero mode localized on the positive tension brane. So we only consider the left-chiral fermion (the right-chiral fermion will have the same mass spectrum with $m_{n}>0$). The solution of the left-chiral fermion zero mode is

[TABLE]

For the case of infinite extra dimension, in order to localize the zero mode $f_{L0}(z)$ on the brane, the following normalization condition should be satisfied:

[TABLE]

which turns out to be $v<\frac{1}{3}\left(k+\alpha\right)~{}(<0)$ and $\eta>0$. Next, we consider the case of finite extra dimension, for which we only need $\eta v<0$ in order for the zero mode $f_{L0}$ to be localized on the positive tension brane. We will investigate the localization and mass spectrum for the special case $v=\frac{1}{3}(k+\alpha)<0$. The effective potential for the left-chiral fermion KK modes is reduced to

[TABLE]

Now, in order to transform Eq. (68) into a dimensionless one, we redefine the following dimensionless potential:

[TABLE]

We assume the KK modes satisfy the Neumann boundary condition $\partial_{\bar{z}}(e^{-2A(\bar{z})}f_{Ln})=0$ at the boundaries $\bar{z}=0$ and $\bar{z}=\bar{z}_{b}$. By using the boundary condition at $\bar{z}=0$, we get the general solution of Eq. (68),

[TABLE]

where

[TABLE]

Another boundary condition will give the mass spectrum of the fermion KK modes:

[TABLE]

where $B_{P_{\text{f}}}$ and $B_{P_{\text{f-1}}}$ are defined as follows:

[TABLE]

We can obtain the mass spectrum of the left-chiral fermions by numerical calculation, and show some results in Figs. 9, 10, and 11. Figure 9 shows that the mass of the first massive mode increases with the value of the Yukawa coupling constant $\eta$, and the number of the excited states in a single period increases with the size of the extra dimension $\bar{z}_{b}$. Figure 10 indicates that, for the case of $\bar{z}_{b}<1$, the excited states will not appear in every period and they only emerge after multiple periods. We plot the mass spectrum in Fig. 11 for different parameters. It can be seen that the spectrum interval approaches a constant for the higher excited states due to the plane wave behavior of the wave function, while it is relatively sparse for the lower excited states.

IV Correction to Newtonian potential

Finally, we consider the tensor fluctuations in this model. Now we add a small perturbation $a^{2}(z){h}_{\mu\nu}(x,z)$ to the background metric (8):

[TABLE]

Following Ref. dewolfe , we suppose the tensor perturbation ${h}_{\mu\nu}$ satisfies the transverse-traceless (TT) condition: ${h}^{\mu}_{~{}\mu}=\partial^{\nu}{h}_{\mu\nu}=0$. The equation of motion for ${h}_{\mu\nu}$ is given by Yang:2011pd

[TABLE]

Furthermore, ${h}_{\mu\nu}$ can be decomposed in the form

[TABLE]

where $H(z)=(1+\beta|z|)^{1/3}$. The four-dimensional mass $m_{n}$ of a KK excitation is defined by the Klein–Gordon equation $\Box^{(4)}{\varepsilon^{(n)}_{\mu\nu}(x)}=m_{n}^{2}\varepsilon^{(n)}_{\mu\nu}(x)$. We obtain a Schr$\ddot{\text{o}}$dinger-like equation from Eq. (87),

[TABLE]

where the effective potential $V(z)$ is given by

[TABLE]

The general solution of Eq. (89) is a linear combination of the Bessel functions,

[TABLE]

Imposing the Neumann boundary condition $\partial_{z}{h}_{\mu\nu}(x,z)=0$ at the boundaries $z=0$ and $z=z_{b}$, we get

$\alpha_{n}=-J_{1}(\frac{m_{n}}{\beta})/Y_{1}(\frac{m_{n}}{\beta})$ and

[TABLE]

The solution of the graviton zero mode $\Psi_{0}(z)$ is

[TABLE]

which is localized on the negative tension brane for finite $z_{b}$. Note that the graviton zero mode cannot be normalized anymore when the extra dimension is infinite, which is very different from the RS1 model. The massless zero mode sector of the action (1) leads to us the effective four-dimensional gravitational theory,

[TABLE]

where $\hat{R}^{(4)}$ is the four-dimensional Riemannian Ricci scalar made out of $g^{(4)}_{\mu\nu}=\eta_{\mu\nu}+h^{(0)}_{\mu\nu}(x)$. Thus the effective four-dimensional Planck scale $M_{Pl}$ is Yang:2011pd

[TABLE]

For the light modes ($m_{n}/\beta\ll 1$), the solution (91) in the long-range limit $\alpha_{n}\approx\frac{\pi m_{n}^{2}}{4\beta^{2}}\ll 1$ becomes

[TABLE]

Correspondingly, the spectrum of KK gravitons is Yang:2011pd

[TABLE]

where $x_{n}$ satisfies $J_{1}{(x_{n})}=0$, and $x_{1}$=3.83, $x_{2}$=7.02, $x_{3}$=10.17, $\cdots$. The normalization constant $N_{n}$ is determined by guobin Add

[TABLE]

Imposing the approximation $J_{0}(x)\approx\sqrt{\frac{2}{\pi x}}\cos(x-\frac{1}{4}\pi)$ and the approximate formula of the zero point of $J_{1}(x_{n})$, $x_{n}\approx(n+\frac{1}{4})\pi$, we get $N_{n}\approx\sqrt{\frac{\pi x_{n}}{2z_{b}^{2}}}$. Then the normalized KK modes are

[TABLE]

For simplicity, the effective gravitational potential between two point-like sources with mass $M_{1}$ and $M_{2}$ is obtained by the exchange of the graviton zero mode and massive KK modes, and it can be expressed as Randall2 ; Y.shirman

[TABLE]

where the four-dimensional Newtonian potential is caused by the zero mode, and the correction term is produced by the massive KK modes. Combining Eqs. (95), (99), and (100), we can get the correction term to the Newtonian potential:

[TABLE]

For $r\gg z_{b}$ , the summation term tends to $e^{-x_{n}r/z_{b}}\approx 0$, and the correction can be ignored. For the case of $r\ll z_{b}$ , the correction term cannot be ignored and it can be calculated as follows:

[TABLE]

It shows that, for $r\ll z_{b}$, the Newtonian potential is corrected by a $1/r^{3}$ term, which is the leading term. For $r\gg z_{b}$, the correction can be ignored, and we recover the Newtonian potential $U(r)\propto 1/r$. According to the gravitational experiments of the correction to the Newtonian potential Eingorn11 ; Iorio11 ; Linares13 ; Luo2016 , we know that the size of the extra dimension, $z_{b}$, should not be larger than the micron scale.

V Conclusion and summary

To summarize, we have investigated the localization and mass spectra of various bulk matter fields on the Weyl brane. We first gave a brief review of the thin brane arising from a five-dimensional Weyl integrable spacetime. Then, we investigated localization of the zero modes for various bulk matter fields (i.e., scalar, vector, and fermion fields) on the brane, and we got the mass spectra of the fields by redefining some dimensionless parameters. We also considered the correction to the four-dimensional Newtonian potential from the massive KK gravitons.

It was found that the zero modes of various bulk matter fields can be localized on the positive tension brane under some conditions, which are collected in Table. 1. It can be seen that the localization conditions for the case of infinite extra dimension are stronger than the case of a finite extra dimension. When the extra dimension is finite, the scalar and vector zero modes can be localized on the positive tension brane even if there is no interaction with the background scalar field (i.e., $\lambda=0$ and $\tau=0$).

For the scalar field, we considered two types of couplings with the background scalar field. For the case of the scalar–dilaton coupling, we found that the mass of the first massive scalar KK mode increases and decreases with the dimensionless coupling constant $\lambda$ and the size of the extra dimension $\bar{z}_{b}$, respectively, and the gap of the mass spectrum decreases with $\bar{z}_{b}$. When the dimensionless size $\bar{z}_{b}>1$, the number of the excited states in a single period, shown in Fig. 2, increases with $\bar{z}_{b}$. When $\bar{z}_{b}<1$, the excited states do not appear in each period, they will emerge after several periods. For the case of the Higgs potential coupling, we fixed the parameter $\bar{\theta}$ and chose the proper value $\bar{u}$, to ensure the zero mode can be localized on the brane. In order to avoid negative $\bar{m}^{2}_{n}$, the mass parameter $\bar{u}$ in the Higgs potential should have an upper limit when the other parameters are fixed.

For the vector field, it was shown that the mass of the first massive vector KK mode increases with the coupling constant $\tau$, and the gap of the mass spectrum decreases with the size of the extra dimension.

For the fermion field, we introduced the usual Yukawa coupling with $F(\omega)=\partial_{z}e^{\upsilon\omega}$, and we found that the left-chiral fermion zero mode can be localized on the positive tension brane at some conditions. We calculated the mass spectrum of the fermion KK modes for $v=\frac{1}{3}(k+\alpha)$ as an example. The mass of the first excited state increases with the Yukawa coupling constant $\eta$. The size of the extra dimension can also affect the mass gap in the same way as the cases of the scalar and vector fields. The spectra of various bulk matter fields show the same phenomenon: that the spectrum interval approaches a constant for the higher excited states, while it is relatively sparse for the lower excited states.

For the TT tensor perturbation of the gravitational field, the zero mode is localized on the negative tension brane for the finite extra dimension, and cannot be normalized anymore when the extra dimension is infinite. For the case of the finite extra dimension, the gravity zero mode gives the Newtonian potential, while the gravity massive KK modes will give a correction to the Newtonian potential by a $1/r^{3}$ term when $r\ll z_{b}$.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grants Nos. 11522541 and 11375075) and the Fundamental Research Funds for the Central Universities ((Grants No. lzujbky-2016-k04)).