From sine-Gordon to vacuumless systems in flat and curved spacetimes

TL;DR

This paper introduces a family of models interpolating between sine-Gordon and vacuumless systems using a deformation method, and explores their implications in thick brane scenarios, including stability and gravitational geometry transitions.

Contribution

It presents a novel deformation approach to generate models transitioning from well-minimized to non-minimized systems, with applications to thick brane and gravitational geometry.

Findings

Reproduces sine-Gordon and vacuumless models via deformation.

Analyzes stability and renormalization group flow of brane scenarios.

Shows gravity transitions from warped to Minkowski geometry.

Abstract

In this work we start from the Higgs prototype model to introduce a new model, which makes a smooth transition between systems with well located minima and systems that support no minima at all. We implement this possibility using the deformation procedure, which allows the obtention of a sine-Gordon-like model, controlled by a real parameter that gives rise to a family of models, reproducing the sine-Gordon and the so-called vacuumless models. We also study the thick brane scenarios associated with these models and investigate their stability and renormalization group flow. In particular, one shows how gravity can change from the 5-dimensional warped geometry with a single extra dimension of infinite extent to the conventional 5-dimensional Minkowski geometry.

Equations

$$\mathcal{L}(\varphi ,{\partial }_{\mu }\varphi )=\frac{1}{2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi -V(\varphi ).$$

$${\partial }_{\mu }{\partial }^{\mu }\varphi +\frac{dV}{d\varphi }=0.$$

$$\frac{d^2\varphi }{d{x}^2}=\frac{dV}{d\varphi }.$$

$$T^{\mu \nu }={\partial }^{\mu }\varphi {\partial }^{\nu }\varphi -{\eta }^{\mu \nu }\mathcal{L}.$$

$$\rho (x)$$

$$=$$

$$V(\varphi )=\frac{1}{2}{W}_{\varphi }^2,$$

$$\rho (x)=\frac{1}{2}{\left({\varphi }^{\prime }\mp W_{\varphi }\right)}^2\pm \frac{dW}{dx}$$

$$\frac{d\varphi }{dx}=\pm W_{\varphi }.$$

$$E_{BPS}$$

$$V(\varphi )=\frac{1}{2}(1-{\varphi }^2{)}^2,$$

$$j^{\mu }={ϵ}^{\mu \nu }{\partial }_{\nu }W(\varphi (x)),$$

$$Q={\int }_{-\infty }^{\infty }{j}^{0}dx=W(\varphi (\infty ))-W(\varphi (-\infty )).$$

$$\varphi \to f(\chi ).$$

$$\mathcal{L}(\varphi ,{\partial }_{\mu }\varphi )=f_{\chi }^2\mathcal{L}(\chi ,{\partial }_{\mu }\chi )$$

$$\mathcal{L}(\chi ,{\partial }_{\mu }\chi )=\frac{1}{2}{\partial }_{\mu }\chi {\partial }^{\mu }\chi -U(\chi ),$$

$$U(\chi )=\frac{V\left(\varphi \to f(\chi )\right)}{f_{\chi }^2}.$$

$$\frac{d\chi }{dx}={\bar{W}}_{\chi },$$

$$E=|\bar{W}\left(\chi (\infty )\right)-\bar{W}\left(\chi (-\infty )\right)|,$$

$$\chi (x)=f^{-1}(\varphi (x)),$$

$$f_{\lambda }(\chi )=\tanh \left(\frac{1}{\theta \sqrt{2-{\theta }^2}}{\text{tanh}}^{-1}\left(\theta \frac{\text{sc}\left(\chi ,\lambda \right)}{\sqrt{2-{\theta }^2}}\right)\right)$$

$$\text{sc}(\chi ,\lambda )=\frac{\text{sn}(\chi ,\lambda )}{\text{cn}(\chi ,\lambda )},$$

$${\text{cn}}^2(\chi ,\lambda )+{\text{sn}}^2(\chi ,\lambda )$$

$${\text{dn}}^2(\chi ,\lambda )+\lambda {\text{sn}}^2(\chi ,\lambda )$$

$$\text{U}(\chi ,\lambda )=\frac{1}{2}\frac{{\left(2{\text{cn}}^2(\chi ,\lambda )-{\theta }^2\right)}^2}{{\text{dn}}^2(\chi ,\lambda )}.$$

$$\text{W}(\chi ,\lambda )$$

$$U(\chi ,0)$$

$$U(\chi ,1)$$

$${\chi }^{\prime }=\frac{2\text{cn}(\chi ,\lambda {)}^2-{\theta }^2}{\text{dn}(\chi ,\lambda )}.$$

$${\chi }_{L,\lambda }(x)={\text{sc}}^{-1}\left(\frac{\sqrt{2-{\theta }^2}}{\theta }\tanh \left(\theta \sqrt{2-{\theta }^2}x\right),\lambda \right)+2n{\text{K}}_{\lambda }.$$

$${\chi }_{L,1}(x)$$

$${\chi }_{L,0}(x)$$

Full text

I Introduction

Topological solutions in Field Theory are related to many phenomena in Physics vachaspati ; v ; ms . In particular, when we have models in 1+1 dimensions involving a potential written in terms of a scalar field and with a set of degenerate minima, such solutions represent transitions between consecutive minima and are called kinks. The minima define the possible vacua states at the quantum level, and are distributed over the many values ​​for the field. Such distribution can appear in the most diverse ways, for instance, in the $\phi^{4}$-model one has only one topological sector, defined between two consecutive minima, while in the sine-Gordon model caudrey one has an infinite copy of the same sector, always between two consecutive and well localized minima. A system that is quite different from this perspective is the vacuumless model vilenkin1 ; vilenkin2 . In this model we still have a topological sector, which can be interpreted as connecting two minima of the scalar potential, but now they are located at infinity. In this case, the field solution is asymptotically divergent and has infinite amplitude, but keeps its topological character well behaved bvacless . It is worth mentioning that vacuumless systems appear in a diversity of contexts in high energy physics vilenkin1 ; vilenkin2 ; bvacless ; bbn ; bbl ; bbg ; dutra .

The aim of this paper is to construct a model that, with the proper variation of a given parameter, makes the transition between systems with well-located minima and systems that support no minima at all. The model presented here has as limit cases the sine-Gordon and the vacuumless models. Using the deformation procedure developed in Ref. blm , we find a field transformation that takes us to a new system with a double sine-Gordon-like behavior, which has two infinitely degenerate sets of solutions. One of these sets transits between the sine-Gordon kink and the vacuumless solution. The second set coincides with the previous one in one of its limits and, except from a phase, in the other limit case it is destroyed, remaining only the zero energy solutions. Thus, in one of the limits of the model the two phases coincide, generating infinite copies of the same topological sector, and in the other limit only one topological sector survives.

Models in Field Theory motivate generalizations of the Randall-Sundrum model rs1 ; rs2 in the presence of scalar fields gw ; fre ; csaki ; gremm ; d ; bm ; bmm , which are known as thick branes. For this reason, in this work we also study the thick brane scenario generated by the scalar field model which we first introduce and study in the flat spacetime. In particular, we find a brane that transits between the sine-Gordon brane gremm and a flat 5-dimensional spacetime with two zero energy solutions for the scalar field. The idea here is similar to the case investigated in bc , in which the authors propose a braneworld scenario where the brane changes from a thick to a thin behavior. In the current study, however, we describe a mechanism in which a single parameter can be used to control the brane profile, contributing to change the 5-dimensional warped geometry into a flat geometry.

Although the braneworld model that we explore below is more involved, it also supports analytical solutions. Thus, in the 5-dimensional spacetime with a single extra dimension of infinite extent, we also analyze the stability of the braneworld scenario against tensorial perturbations. In the sense of Gauge/Gravity Duality maldacena ; skenderis1 , where the extra dimension can be identified with the energy scale of the holographic dual field theory, we also study its implications for the renormalization group flow (the RG flow) in the dual Field Theory, since stability of the gravitational sector has relevant information about the dual model skenderis .

The subject to be explored in the current work is organized as follows. In Sec. II we review several aspects of the first order formalism for a single real scalar field and its relationship with the Bogomol’nyi-Prasad-Somerfield (BPS) solutions bps . Then, we discuss the deformation procedure and show how it acts to generate the new model in the 2-dimensional spacetime. We go on and study in Sec. III the new model and its topological solutions, including how the energy and energy density behave. In addition, we also show that the solutions obtained are stable. Moreover, in Sec. IV we analyze the properties of the thick branes that can be constructed from the model. We also investigate stability of the brane against metric fluctuations and study implications of the RG flow for the dual Field Theory. We end the work in Sec. V, adding some comments and conclusions.

II Generalities

II.1 First order formalism

The behavior of a scalar field is usually encoded in a Lagrangian density having the general form

[TABLE]

Here $\phi$ is the scalar field and $V(\phi)$ is the potential of the model, which determines how the field behaves. In the flat spacetime with $(1,1)$ spacetime dimensions, the metric tensor becomes $\eta_{\mu\nu}=\text{diag}(1,-1)$, so the scalar field only depends on the two coordinates $x$ and $t$, i.e., $\phi=\phi(x,t)$. For simplicity, we also work with dimensionless fields and coordinates. The equation of motion for the scalar field derived from the Lagrangian (1) is given by

[TABLE]

As the Lagrangean (1) is Lorentz invariant, we can focus on static solutions, since traveling waves can be obtained from a Lorentz boost. For static configurations the equation (2) becomes a second order differential equation given by

[TABLE]

Another important quantity we are interested in is the energy-momentum tensor

[TABLE]

In particular, its $T^{00}$ component provides the energy density of the solution we are looking for. We represent it by $\rho(x)$, which is explicitly given by

[TABLE]

Note that for positive-definite energy, the potential must be non-negative, i.e., $V(\phi)\geq 0$. A powerful tool in the treatment of these models is the use of an auxiliary function, denoted by $W(\phi)$, which is introduced as follows

[TABLE]

where $W_{\phi}={dW}/{d\phi}$. In this case the expression for the energy density (5) becomes

[TABLE]

and the equation (3) can now be given as the first order differential equations

[TABLE]

It implies that the quadratic term in (7) disappears and, as a consequence, the energy density is only related to the $x$-derivative of $W$. Thus, the energy of the model is determined only by the asymptotic behavior of the function $W$ in the coordinate space; that is, one can write

[TABLE]

In this case the energy (9) is called BPS energy bps .

Relevant phenomena occurs when the system under analysis presents a set of degenerate minima. In these situations each pair of consecutive minima form distinct topological sectors that, in turn, have different solutions. These solutions are called kinks. The simplest case of models having such properties is the well-known Higgs prototype or $\phi^{4}$-model, defined by the potential

[TABLE]

which has two degenerate minima at $\phi=\pm 1$ and a topological sector having a kink solution given explicitly by $\phi(x)=\tanh(x)$.

An interesting way to characterize a kink is the existence of a topological current. Here, we define it by bvacless

[TABLE]

where $\epsilon^{\mu\nu}$ is the antisymmetric symbol in two dimensions with $\epsilon^{01}=1$. Associated with this current (11) we have a topological charge given by

[TABLE]

Despite the similarities in the values of the charge (12) and the BPS energy (9), they have fundamental differences. While the BPS energy is associated with a continuous symmetry and can be identified from the Noether Theorem, the charge (12) is associated with the topology of each solution, and results from the transition of the topological solution in between two minima of the potential.

II.2 The deformation procedure

Once we know the behavior of a given model, with its characteristics and general behavior, it is interesting to look for new well-behaved models. In this sense, the Deformation Method blm is a powerful tool to find new models in Field Theory. The method consists of choosing a theory with the Lagrangian ${\cal L}(\phi,\partial_{\mu}\phi)$ having the form (1) and then perform a transformation of the type

[TABLE]

We then get

[TABLE]

where

[TABLE]

with the potential $U(\chi)$ given by

[TABLE]

This is a field redefinition, but if we consider the model (15) described by the potential (16) as a new model, and call it the deformed model, in this case there is a $\bar{W}$ function such that the first-order differential equations

[TABLE]

with $\bar{W}_{\chi}=W_{\phi}(\phi\to f(\chi))/f_{\chi}$, provide solutions for the field $\chi(x)$. In particular, among the characteristics of the new model, we can highlight that the total energy of the solution is

[TABLE]

since the deformed model also has a first order structure. However, the deformed model engenders another important property: the transformation (13) allows us to find the solution for the field $\chi(x)$ without the need to deal directly with equation (17). What we have is that the solution to the new field is obtained by inverting the transformation (13). Thus, one finds

[TABLE]

where $\phi(x)$ represents a solution for the previous model. In this sense, the equation (19) represents the main link between the $\phi$-model and the deformed $\chi$-model. For more details, see blm .

III Model

In this work we follow the deformation procedure and introduce a new model generated by the deformation function

[TABLE]

which is applied to the $\phi^{4}$ model (10). Here $\theta$ is a real parameter $(\neq\pm\sqrt{2})$ and the function $\text{sc}^{-1}(\chi,\lambda)$ is one of the Jacobi Elliptic Functions, defined by the ratio

[TABLE]

where $\text{sn}(\chi,\lambda)$ and $\text{cn}(\chi,\lambda)$ are the Jacobi elliptic sine and cosine, respectively. Here we have

[TABLE]

where $\lambda$ is a parameter in the interval $[0,1]$. In particular for $\lambda=0$ we have $\text{sn}(\chi,0)=\sin(\chi)$, $\text{cn}(\chi,0)=\cos(\chi)$ and $\text{dn}(\chi,0)=1$, where we retrieve the basic trigonometric relations, and for $\lambda=1$ we have $\text{sn}(\chi,1)=\text{tanh}(\chi)$ and $\text{cn}(\chi,1)=\text{dn}(\chi,1)=\text{sech}(\chi)$, which lead us to the hyperbolic functions. The new model found from the deformation function (20) and the potential (10) is

[TABLE]

Its behavior is depicted in Fig. 1. The potential (23) has $\mathbb{Z}_{2}-$symmetry and is invariant under transformation $\phi\to\phi+2K_{\lambda}$, where $K_{\lambda}=\text{cn}^{-1}(0,\lambda)$. In this case we find $W_{\chi}=\pm(2\text{cn}^{2}(\chi,\lambda)-\theta^{2})/\text{dn}(\chi,\lambda)$, which implies that the $W$ function of the model is

[TABLE]

where $\text{am}(\chi,\lambda)=d\left(\text{dn}(\chi,\lambda)\right)/d\chi$ is the Jacobi amplitude. Although we initially present the model with two parameters, we are interested in situations where only one parameter is necessary to describe the changes in the model. So we assume that $\theta=\theta(\lambda)$. We also impose on $\theta(\lambda)$ the conditions $\theta(0)=1$ and $\theta(1)=0$ in order to find the vacuumless model bvacless for $\lambda=1$ and the sine-Gordon model caudrey for $\lambda=0$. As a consequence we obtain the particular cases $W(\chi,1)=4\tan^{-1}(e^{\chi})-\pi$ and $W(\chi,0)=\sin(2\chi)/2$, so the sine-Gordon and vacumless models appear in the system as

[TABLE]

The parameter $\lambda$ has an interesting behavior: when it increases from [math] to unity, the model changes from the sine-Gordon to the vacuumless model. Physically, it transforms the periodic sine-Gordon potential (25a) to a non-periodic one, the hyperbolic potential (25b) which defines the vacuumless model. The distinct solutions appear in Fig. 2 and in Fig. 3 one illustrates how the energy density becomes more and more diffuse, as $\lambda$ increases from zero to unity.

Once the model has been presented, we must solve the first order equation

[TABLE]

Equation (26) has two infinite set of solutions. The first one we call large kinks and are given by

[TABLE]

The behavior of the solution (27) is depicted in Fig. 2 for $n=0$. It asymptotically approaches $\chi_{L,\lambda}(\pm\infty)=\pm\text{sc}^{-1}(\sqrt{2-\theta^{2}}/\theta,\lambda)$ and in the vicinity of the origin behaves like $\chi_{L,\lambda}(x\simeq 0)\simeq(2-\theta^{2})x+\mathcal{O}(x^{2})$. We note that as $\theta\rightarrow 0$, $\chi_{L,\lambda}(\pm\infty)$ tends to diverge and $\chi_{L,\lambda}(x\simeq 0)$ remains well-behaved. For $\theta\to 1$ we do not have pathologies in the solution. Thus, the large kink describes, as a function of the $\lambda$ parameter, systems that transit between the sine-Gordon kink and the vacuumless solution. Particularly, we have

[TABLE]

The second set of solutions we find from equation (26) is

[TABLE]

and we call it small kink. Its shape is depicted in Fig. 3 for $n=0$, except for the phase $K_{\lambda}$. The solution (29) asymptotically approaches $\chi_{S,\lambda}(\pm\infty)=\pm\text{sc}^{-1}(\theta/\sqrt{(2-\theta^{2})(1-\lambda)},\lambda)+K_{\lambda}$ and in the neighborhoods of $x=0$ behaves as $\chi_{S,\lambda}(x\simeq 0)\simeq K_{\lambda}+\theta^{2}x/\sqrt{1-\lambda}+\mathcal{O}(x^{2})$. The small kink has a phase $K_{\lambda}$, which causes its topological sector to move away from the center of the potential (23) to infinity as $\lambda\to 1$. With a suitable choice for $\theta(\lambda)$ we can destroy the topological sectors associated with small kinks when $\lambda=1$. Thus, only the topological sector of the vacuumless solution remains at that point. If we drops out the phase in the small kink (29), it is possible to note that the zero energy solutions also remains, but at infinity. A simple choice in this direction that also obeys the conditions $\left(\theta(0),\theta(1)\right)=(1,0)$ is $\theta(\lambda)=\sqrt{1-\lambda}$. With this choice for $\theta(\lambda)$ we can explicitly rewrite the solutions of large and small kinks as

[TABLE]

So now we have a double sine-Gordon-like model with two field solutions coming from the two manifest topological sectors. Both kinks retrieve the sine-Gordon solution when $\lambda=0$, but as $\lambda$ grows, such solutions have distinct properties. Large kink becomes more diffuse until reaching the solution of the vacuumless model, which has divergent amplitude. It implies that at this point all topological sectors of the model are sent to infinity, except for the sector that is at the center of the potential (23). Small kink also becomes more diffuse as $\lambda\to 1$, but is destroyed when $\lambda=1$. As we shall see later, it happens because as the associated topological sector moves away from the center of the potential its energy approaches zero. Moreover, one can show that the mass of the meson in the minima of the potential is given by $m_{\lambda}^{2}=U^{\prime\prime}(\chi=\chi_{min})=4(1-\lambda^{2})$, where $\chi_{min}$ is a minimum of (23). Here, we choose to write the model in terms of that quantity, whenever possible.

We can now perform the analysis of the energy densities of the model. The energy density of the large kink (27) is given by

[TABLE]

and its shape is depicted in Fig. 4. It makes the transition between the curves $\rho_{L,1}(x)=4/(1+4x^{2})$ and $\rho_{L,0}(x)=\text{sech}^{2}(2x)$. Asymptotically equation (31) decay as $\rho_{L,\lambda}(x\to\infty)\simeq 8(1-\lambda)(1+\lambda)^{2}e^{-2m_{\lambda}x}/(2+\lambda)+\mathcal{O}(e^{-4m_{\lambda}x})$, where we found the particular case $\rho_{L,0}(x\to\infty)\simeq 4e^{-4x}$. We don’t have information for the case $\lambda=1$ in the general asymptotic expansion for the energy density, but a direct approach in $\rho_{L,1}(x)$ lead us to the expression $\rho_{L,1}(x\to\infty)\simeq 1/x^{2}+\mathcal{O}(x^{-4})$, which decays much slower than the exponential. This change in the asymptotic behavior is due to the mass scale of the quantum meson which, for $\lambda=1$, is zero. On the other hand the behavior of (31) in the vicinity of $x=0$ is given by $\rho_{L,\lambda}(x\simeq 0)\simeq(1+\lambda)^{2}-(1+\lambda)^{3}(4-\lambda^{2}-\lambda)x^{2}+\mathcal{O}(x^{3})$, which shows that the central portions of the energy density increases and becomes a bit more concentrated as $\lambda\to 1$.

For the small kink we have the energy density given by

[TABLE]

and its behavior is depicted in Fig. 5. Equation (32) transits between the curves $\rho_{S,1}(x)=0$ and $\rho_{L,0}(x)=\text{sech}^{2}(2x)$. Asymptotically it decay as $\rho_{S,\lambda}(x\to\infty)\simeq 8(1-\lambda)(1+\lambda)^{2}e^{-2m_{\lambda}x}/(2+\lambda)+\mathcal{O}(e^{-4m_{\lambda}x})$, revealing that we can not distinguish the energy densities of the fields (27) and (29) when $x\to\infty$. For $x\simeq 0$ we have $\rho_{S,\lambda}(x\simeq 0)\simeq(1-\lambda)-(1-\lambda)^{2}(\lambda+4)x^{2}+\mathcal{O}(x^{3})$. It shows that, despite the similar asymptotic behavior (for $\lambda\neq 1$), the evolution of (31) and (32) in terms of the parameter $\lambda$ is very different. The height and width of the energy density (31) grow with $\lambda$, becoming more diffuse but still with a localized profile. Meanwhile, as $\lambda$ grows, the width of the energy density (32) increases, but its central portion decreases and becomes less concentrated, which implies that as the solution evolves it becomes more diffuse, but also delocalized. As a consequence, the area under the energy density (32) become smaller until it disappears, at $\lambda=1$.

Integrating (31) and (32) we find the energy of large and small kinks, which are given by

[TABLE]

and

[TABLE]

respectively. The corresponding behaviors are depicted in Fig. 6. Here one observes that the energy of the solutions (27) and (29) are bounded. The expression (33) is a monotonically increasing function of the $\lambda$ parameter confined in the interval $[1,2\pi]$, where $E_{L}(0)=1$ and $E_{L}(1)=2\pi$, and the expression (34) is a monotonically decreasing function of the $\lambda$ parameter confined in the interval $[0,1]$, with $E_{S}(0)=1$ and $E_{S}(1)=0$. Moreover, it is easy to show that (33) and (34) are related by $E_{L}(\lambda)=E_{S}(\lambda)+\frac{\pi}{\lambda}\left(2-(2+\lambda)\sqrt{1-\lambda}\right)$.

III.1 Linear stability

In this section we analyze the stability of the solutions of the models presented so far. The usual procedure is to take a time-dependent perturbation around the static solution written in the form $\chi(x,t)=\chi(x)+\sum_{n}\eta_{n}(x)\cos(\omega_{n}t)$, for small $\eta_{n}(x)$, and then insert $\chi(x,t)$ into (2). The procedure gives

[TABLE]

which is a Schrödinger-like equation with a stability potential given by

[TABLE]

Inserting the large kink solution (27) into (36) we find

[TABLE]

which is depicted in Fig. 7. For $\lambda=1$ we have $v_{L,1}(x)=4\left(8x^{2}-1\right)/\left(4x^{2}+1\right)^{2}$, which is a volcano potential, and for $\lambda=0$ we have $v_{L,0}(x)=4-8\text{sech}^{2}(2x)$, which is a reflectionless potential. Equation (37) has a global minimum at $v_{L,\lambda}(0)=(1+\lambda)\left(\lambda^{2}+\lambda-4\right);$ this minimum increases or decreases, depending on $\lambda$ being above or below the point at $\tilde{\lambda}=\frac{1}{3}\left(\sqrt{13}-2\right)$. It shows that the potential (37), as $\lambda$ grows, becomes deeper and after $\lambda=\tilde{\lambda}$ it comes back to its initial deepness.

The change in the shape of the stability potential, in this case, is due to the behavior of the meson mass in the minima of (23). As $\lambda\to 1$ we have $m_{\lambda}^{2}\to 0$. Note that $v_{L,\lambda}(0)<0$ for all $\lambda$ and translational invariance requires the existence of at least one bound state, so the transition between the reflectionless and the volcano shapes in this case describes the transition from systems with massive meson to systems having a massless meson. This is another relevant physical behavior induced by the parameter $\lambda$, which will lead to distinct possibilities when used to describe braneworld scenarios, as we discuss in Sec. IV.

Now, inserting the small kink solution (29) in the general expression for the stability potential (36), we find

[TABLE]

which is depicted in Fig. 8. For $\lambda=1$ we have $v_{S,1}=0$ and for $\lambda=0$ we have $v_{S,0}=4-8\text{sech}^{2}(2x)$, as expected. At the center we have $v_{S,\lambda}(0)=(\lambda-1)(\lambda+4)$, showing that the potential (38) becomes shallower as $\lambda$ grows, finally desappearing when $\lambda=1$. Now, as $\lambda$ grows, the stability potential keep its reflectionless shape. As a consequence, the bound states of (38) becomes less expressive as $\lambda\to 1$ and disappears when $\lambda=1$.

Asymptotically, both (37) and (38) approaches $v_{L,\lambda}(\pm\infty)=v_{S,\lambda}(\pm\infty)=m_{\lambda}^{2}$. Moreover, one notes that the hamiltonian ${H}=-d^{2}/dx^{2}+v(x)$ can be rewritten as ${H}=S^{{\dagger}}S$ with $S^{{\dagger}}=-d/dx-W_{\phi\phi}$ by using the equation (12). Thus we have $\omega_{n}^{2}\geq 0$ for all values of $n$, since $H$ is non-negative. As equation (35) does not admit negative energy modes, the stability of the solution is then ensured.

The translational invariance of the solutions we presented so far implies the existence of at least one bound state for each topological sector, which is given by the zero mode of the equation (35). In formula, if we represent the zero-mode by $\eta_{0}(x)$, it is the derivative of the field solution,

[TABLE]

For large and small kinks the shapes are depicted in Fig. 9, where one can observe how (37) holds the nice behavior of its zero-mode, and how the zero mode of (38) gets smaller until it disappears when $\lambda=1$.

IV Braneworld

Models described by scalar fields have direct applications in Gravitation, providing braneworld scenarios for thick branes. In this context, the scalar field acts as a source of gravity around the brane, and thus describes how gravity behaves throughout the bulk. The system we are interested consists of a 3-brane embedded in a (4+1) spacetime with an extra dimension of infinite extent. The background geometry can be written in terms of a static warped metric given by

[TABLE]

Here $a,b=0,...,4$, $\mu,\nu=0,...,3$, $ds^{2}_{4}=\eta_{\mu\nu}dx^{\mu}dx^{\nu}$, and the $y$-coordinate describes the extra spatial dimension. The functions $A(y)$ and $e^{A(y)}$ are called warp function and warp factor, respectively, and are assumed to depend only on the extra dimension.

In this braneworld scenario, we are interested in models that can be described by the action

[TABLE]

where $\mathcal{L}(\phi,\partial_{a}\phi)=\frac{1}{2}g_{ab}\partial^{a}\phi\partial^{b}\phi-U(\phi)$ is the Lagrangian for the scalar field and, for simplicity, we assume $4\pi G_{5}=1$. We also assume that $\phi=\phi(y)$, i.e., the scalar field only depends of the extra dimension.

The Einstein equations that follows from the action (41) are

[TABLE]

where $G_{ab}$ is the Einstein tensor and $T_{ab}$ is the energy-momentum tensor, similar to the expression (4). The $00$ and the $44$ components of (42) are given by, respectively

[TABLE]

with the prime denoting derivation in respect to the coordinate $y$. With equations (43) at hands we can subtract the first equation from the second to obtain

[TABLE]

Equation (44) provides a way to rewrite the system in terms of first order equations. For this, we introduce a function $W(\phi(y))$ in the equations through the relation

[TABLE]

As a consequence, the equation providing the solution for the scalar field is now

[TABLE]

In order to solve the equations of motion, these two first order equations requires that the potential obeys

[TABLE]

The above equations (45) and (46) constitute the first order framework and can be used to analyze possible scenarios of thick brane that can be generated by the models presented in the previous section.

The analysis of the thick branes scenario generated by the large kink (27) is similar to the case of small kinks (29), so we concentrate on the brane generated by the small kink. To simplify the investigation, let us perform a shift in the field, $\phi\to\phi-K_{\lambda}$. In this case, one gets

[TABLE]

Due to the shift performed over $\phi$, the field solution of (46) becomes

[TABLE]

Its shape can be observed in Fig. 3, when we make the changes $x\to y$ and $\phi_{\lambda}(y)=\chi_{S,\lambda}(y)-K_{\lambda}$. Note that for $\lambda=1$ the field approaches a constant, $\phi_{\lambda=1}(\pm\infty)=\pm\text{csch}^{-1}\left(\sqrt{2}\right)$.

With these ingredients we analytically solve the equation (45), together with the boundary conditions $A(0)=A^{\prime}(0)=0$, to find

[TABLE]

where

[TABLE]

Here $\text{Li}_{2}(y)$ is the polylogarithmic function. The shape of $e^{2A(y)}$ is depicted in Fig. 10. We can observe that in the vicinity of the center the warp factor behaves like $A(|y|\simeq 0)\simeq(1-\lambda)y^{2}+\mathcal{O}(y^{4})$, which implies that as $\lambda\to 1$, the brane becomes less and less localized, becoming effectively delocalized for $\lambda=1$. Meanwhile, in the asymptotic regime $|y|\to\infty$ the thick brane approaches AdS vacua with cosmological constant $\Lambda_{5}\sim-W(\phi(\pm\infty))^{2}$, and we have $A(|y|\to\infty)\simeq-\frac{1}{3}E_{S}(\lambda)|y|$, where we used the relations (18) and (34). As $E_{S}(\lambda)$ decreases as $\lambda$ grows (starting at $\lambda=0$, where we have the sine-Gordon brane) the thick branes generated by this model approaches AdS vacua with smaller cosmological constant as $\lambda\to 1$. In particular, for $\lambda=1$ the brane seems to fill the entire space, so we have no graviton localization anymore at this point, with all the modes dispersed through the extra dimension. The braneworld scenario is destroyed in the limit $\lambda\to 1$.

From equations (43) we find that the potential has to obey (47) and one notes that it keeps the correct form for vacuum stability in gravitational theories boucher ; townsend . In particular, for $\lambda=1$ the potential $U(\phi)$ vanish, and the spacetime solution reaches a flat 5-dimensional geometry. At this point, only the dynamical term $\sim\partial^{a}\phi\partial_{a}\phi$ survives in the action (41), having as solution the constants $\phi_{\lambda=1}(\pm\infty)$. The presence of two values possible for the scalar field in the flat background is due the original $\mathbb{Z}_{2}$-symmetry, which act as a memory related to the kind of thick brane system from which it is derived.

The energy density of the model is given by

[TABLE]

and its shape is depicted in Fig. 11. It is known that models derived from the first order equations (45) and (46) have zero energy. It happens because we can rewrite the energy density as a total derivative, $\rho(y)=\frac{d}{dy}\left(W\,e^{2A}\right)$. As $W(\phi(\pm\infty))$ is finite and asymptotically the warp function falls off as $e^{-constant^{2}|y|}$, the integral of $\rho(y)$ over all space must vanish. Here we see that when $\lambda=0$, the energy density is well-concentrated around the origin and, as $\lambda$ increases, it becomes more and more diffuse until finally disappearing at $\lambda=1$.

The above results describe an interesting scenario, in which the parameter $\lambda$ may be used to control the physical properties of the 5-dimensional spacetime. If $\lambda$ increases from [math] to unity, it may change the spacetime from a braneworld model with a single extra spatial dimension of infinite extent to a 5-dimensional Minkowski spacetime with no graviton localization.

IV.1 Metric fluctuations

In this section we analyze the stability of the gravitational sector. For this, we perform a redefinition of variable $dy^{2}\rightarrow e^{2A(z)}dz^{2}$ in (40), which allows us to rewrite the metric in a conformally flat scenario $\tilde{g}_{ab}=e^{2A(z)}\eta_{ab}$. With a linear perturbation the metric becomes

[TABLE]

In the transverse-traceless gauge ($\partial_{\mu}h^{\mu\nu}=0$ and $h_{\mu}^{\mu}=0$) the conformal Einstein tensor is ${\bar{G}}_{ab}=-\frac{1}{2}\partial_{c}\partial^{c}h_{ab}$ and the linearized Einstein tensor is given by

[TABLE]

In this way the $\mu\nu-$components of $G_{ab}^{(1)}$ are

[TABLE]

and the linearized energy-momentum tensor becomes

[TABLE]

where the prime denotes the derivative in relation to variable $z$. When using linearized Einstein equations, $G_{\mu\nu}^{(1)}=2T_{\mu\nu}^{(1)}$, we obtain the equation for $h_{\mu\nu}$, which is $-\partial_{c}\partial^{c}h_{\mu\nu}+3A^{\prime}h^{\prime}_{\mu\nu}=0$. At last, the redefinition $H_{\mu\nu}=e^{-ipx}e^{3A/2}h_{\mu\nu}$ allows us to rewrite the equation for $h_{\mu\nu}$ as

[TABLE]

Note that equation (57) has the form of a supersymmetric Schrödinger equation, where the stability potential is

[TABLE]

wich is depicted in (12). We can observe that (58) keep its volcano shape, as usual, but as $\lambda$ grows it becomes less expressive, and we have no stability potential for $\lambda=1$. Equation (57) has the factorized form $S^{+}S^{-}\psi=p^{2}\psi$, where $S^{\pm}=\left(\pm\partial_{z}+3A^{\prime}/2\right)$. In particular, the zero-energy solution of (57) is given by $\psi_{0}(y)=e^{3A(y)/2}$, where the asymptotic behavior of $A(y)$ ensures gravity localization around the brane for $0\leq\lambda<1$. Since the Hermitian operator $S^{+}S^{-}$ is non-negative, we have no normalizable negative gravitons modes, and it ensures system stability. At $\lambda=1$, we find a flat space, so the stability analysis as presented here fails. In this case, what guarantees the stability of this space are the positive mass theorems on asymptotically flat spacetimes, which are valid for dimensions $\leq 7$ schoen .

IV.2 RG Flow

As pointed out in skenderis ; skenderis1 , by rewriting the metric (40) in the form

[TABLE]

where $u=e^{A(y)}$, allows us to interpret the function $u(y)$ as the renormalization scale of some Quantum Field Theory in the sense of Gauge/Gravity Duality maldacena . Note that, due to the boundary conditions on the deformation factor, $u$ is limited to the range $[0,1]=[e^{A(\pm\infty)},e^{A(0)}]$. Thus, domain wall solutions such as those presented here naturally lead us to confining regimes with a UV-cutoff located at $u=1$ in the dual Field Theory. In this context, the solution $\phi(y(u))$ is identified with the running coupling of the system (see bbl ; kiritsis and references therein), so that the $\beta$-function is defined by

[TABLE]

to describe the RG flow in the dual theory.

Denoting the critical points of (47) by $\phi_{\lambda,\infty}$, we can identify them with the zeros of $W_{\phi}$, which lead us to AdS (or flat) vaccua. Expanding the $\beta$-function around its critical point we have $\beta(\phi)\simeq\beta(\phi_{\lambda,\infty})+\beta^{\prime}(\phi_{\lambda,\infty})(\phi-\phi_{\lambda,\infty})+\mathcal{O}\left((\phi-\phi_{\lambda,\infty})^{2}\right)$. Since in the critical points we have $\beta(\phi=\phi_{\lambda,\infty})=0$, we can find the following expression for the running coupling

[TABLE]

Note that if $\beta^{\prime}(\phi_{\lambda,\infty})<0$, $\phi_{\lambda,\infty}$ is a UV fixed point when $u\to\infty$, and if $\beta^{\prime}(\phi_{\lambda,\infty})>0$, $\phi_{\lambda,\infty}$ is a IR fixed point for $u\to 0$. For $\beta^{\prime}(\phi_{\lambda,\infty})=0$ we have a conformal theory.

For the model presented here, we have

[TABLE]

where $\phi_{\lambda,\infty}=\text{sc}^{-1}\left(\frac{1}{\sqrt{1+\lambda}},\lambda\right)$. Note that (62) is a monotonically increasing function of $\lambda$ and we have no divergences in the running coupling since $u\in[0,1]$, but for any $\lambda$ the IR regime $(u\to 0)$ at $\phi=\phi_{\lambda,\infty}$ is well defined in the dual field theory, even for the 5d Minkowski setup. As a consequence, none of the solutions has a conformal dual field theory.

V Comments

In this work we studied a sine-Gordon-like model, which is controlled by a real parameter that continuously connects the sine-Gordon and the vacuumless models. The model appears as a deformation of the $\phi^{4}$ model, and the real parameter is $\lambda$: for $\lambda=0$ one gets the standard sine-Gordon model and for $\lambda=1$ it reproduces the so-called vacuumless model. However, for $\lambda$ in the interval $(0,1)$ one gets a double sine-Gordon model, which contains two distinct topological sector, the large and the small sectors, which give rise to the large and small kinks, respectively.

As it was shown, in the 2-dimensional spacetime, the energy of the large kink in the large sector varies from $E_{L}(\lambda=0)=1$ to $E_{L}(\lambda=1)=2\pi$, and in the case of the small sector one gets $E_{S}(\lambda=0)=1$ and $E_{S}(\lambda=1)=0$. We then see that the small sector, which is degenerate to the large sector at $\lambda=0$, disappears as $\lambda=1$, with the large sector becoming the topological sector of the vacuumless model.

In the 5-dimensional case, we considered a warped geometry with a single extra dimension of infinite extent and studied the new braneworld scenario described in the small sector. In this scenario, the brane energy density is such that the brane energy vanishes, independently of the value of $\lambda$. If we see the model with $\lambda$ increasing from zero to unity, it then nicely describes a way to change a 5-dimensional warped geometry which is asymptotically AdS into a 5-dimensional Minkowski geometry. However, if $\lambda$ is supposed to run in the reverse sense, decreasing from unity to zero, the model could do the reverse, changing the 5-dimensional Minkowski geometry into a braneworld scenario with a warped geometry which is asymptotically AdS5.

As we have shown, the model is stable under tensorial fluctuations in the metric and of current interest, so one should now investigate how it modifies Newton’s law, and how fermion and gauge fields can be entrapped into the brane as $\lambda$ varies in the interval $[0,1]$. Another issue of current interest concerns the variation of $\lambda$: the present investigation cannot tell the value of $\lambda$, so one should search for this considering other arguments. An interesting possibility could be to investigate the conformational entropy associated with the current braneworld model, to see how it behaves as $\lambda$ varies in the interval $[0,1]$. This has been recently investigated in other contexts in e1 ; e2 ; e3 ; e4 and in references therein, and may provide important guide towards the physical realization of gravity localization in the present model, since the entropy could perhaps suggest the better way $\lambda$ should vary, increasing or decreasing in the interval $[0,1]$.

Acknowledgements.

We thank the Brazilian agencies CAPES and CNPq for financial support.