## Abstract

In this work, we consider a Yukawa modification of the Casimir wormhole. With the help of an Equation of State, we impose Zero Tidal Forces. We will examine two different approaches: in a first approach, we will fix the form of the shape function of the Casimir wormholes modified by a Yukawa term in three different ways and finally a superposition of different profiles. In the second approach, we will consider the original Casimir source modified by a Yukawa term in three different ways and we will deduce the form of the shape function In both the approaches the reference energy density will be that of the Casimir source. Connection with the Absurdly Benign Traversable Wormhole are also discussed.

## Introduction

Yukawa in 1935 [1] proposed to describe nonrelativistic strong interactions between nucleons with the help of a potential whose profile is

This is nothing but the screened version of the Coulomb potential with \(\alpha \) describing the strength of the interaction and \(1/\mu \) its range. This short range interaction has captured the interest of many researchers who have adapted it to the Newtonian potential to understand if it has deviations of the same kind. As a result, the Newtonian gravitational potential between two point masses \(m_{1}\) and \(m_{2}\) (atoms for instance) separated by a distance *r*, acquires a Yukawa correction which formally looks like Eq. (3) . Indeed, one gets

where *G* is the gravitational constant. Potentials of the form (2) have been examined from the astrophysical point of view with a particular attention also on the graviton mass [2,3,4]. It is interesting to note that Yukawa-type forces are also predicted in the context of modified gravity theories [5,6,7,8] and also in bigravity theories [11]. To this purpose, it is important to say that the GINGER experiment offers the opportunity of constraining such theories [9]. On the other hand, in a different framework (MOG), it is possible to obtain black holes and traversable wormholes [10]. This MOG predicts also a variation of the Newton’s constant *G*, in such a way to obtain a Yukawa term which enters the metric. Moreover, a Yukawa term seems to be directly involved in the Galaxy Rotation Curves [30]. Even in the context of Casimir effect, deviations of the Newtonian potential of the form (2) have been considered [12,13,14]. It is interesting to observe that a connection between the Casimir forces and the Yukawa profile has been also introduced in Ref. [23] where Van der Waals himself suggested an interaction potential of the form

with *A* and *B* constants of appropriate dimensions. Since there exists a connection between the Casimir and the Van der Waals forces in the case of relatively large separations when the relativistic effects come into play, one can wonder if Yukawa deformations can play a fundamental role even for Traversable Wormholes. To further proceed we need to recall the Einstein’s Field Equations (EFE)

in an orthonormal reference frame. In such a frame the EFE reduce to the following set of equations

in which \(\rho \left( r\right) \) is the energy density,^{Footnote 1}\(p_{r}\left( r\right) \) is the radial pressure, and \(p_{t}\left( r\right) \) is the lateral pressure. The EFE (5), ( 6) and (7) have been obtained with the help of the line element

representing a spherically symmetric and static wormhole [15, 16]. *b*(*r*) is the shape function, while \(\phi \left( r\right) \) is the redshift function. \(\phi (r)\) and *b*(*r*) are arbitrary functions of the radial coordinate \(r\in \left[ r_{0},+\infty \right) \). A fundamental property of a traversable wormhole is that a flaring out condition of the throat, given by \((b-b^{\prime }r)/b^{2}>0\), must be satisfied [15, 16]. Furthermore, at the throat \(b(r_{0})=r_{0}\) and the condition \(b^{\prime }(r_{0})<1\) is imposed to have wormhole solutions. Another condition that needs to be satisfied is \(1-b(r)/r>0\). For the wormhole to be traversable, one must demand that there are no horizons present, which are identified as the surfaces with \(e^{2\phi }\rightarrow 0\), so that \(\phi (r)\) must be finite everywhere. The last condition is satisfied if we adopt a Zero Tidal Forces model (ZTF) represented by \(\phi ^{\prime }\left( r\right) =0\). Such a condition can be imposed by means of an inhomogeneous Equation of State (EoS) of the form

and

where we have used Eqs. (5) and (6). Equations (9) and (10) lead to

In Ref. [17], we have found that the Casimir wormhole described by

does not satisfy the ZTF condition. In this paper, we will consider the Casimir wormhole shape function deformed by a Yukawa profile satisfying also the ZTF condition: in this way we have the possibility of building a new family of solutions which have a vanishing redshift function. We also assume that the Casimir relationship \(\omega =3\) holds, at least on the throat. There exists another reason to consider a Yukawa deformation to the Casimir wormhole. Indeed, in Ref. [18] we have considered a shape function of the form

obeying Eq. (11) with

and therefore satisfying the ZTF property. However, if we simply assume that

where

then the energy density on the throat becomes

namely the Casimir energy density. Moreover, with the help of the relationship (14) on the throat one gets

and by imposing that \(\omega \left( r_{0}\right) =3\), one finds

in agreement with what found in Ref. [36] but with a factor \(\sqrt{3}\) missing. This example suggests that the mixing between the Casimir wormhole and a Yukawa wormhole seems to be promising. The paper is organized as follows: in Sect. 2 we study three different combinations of the Casimir wormhole shape function with a Yukawa term, in Sect. 3 we explore the consequences of a superposition of the profiles considered in Sect. 2, in Sect. 4 we adopt the reverse procedure, namely we fix the form of the energy density and we deduce the form of the shape function, investigating three different profiles. We summarize and conclude in Sect. 5. Units in which \(\hbar =c=k=1\) are used throughout the paper and will be reintroduced whenever it is necessary.

## Casimir–Yukawa wormholes

The Casimir wormhole obtained in Ref. [17] has as a source the original Casimir energy density with a slight but fundamental difference: the plates separation has been promoted to be a variable instead of being a fixed quantity. To satisfy the EFE a non vanishing redshift function has been computed described in Eq. (12). In this section we are interested in examining some modifications of the original Casimir wormhole shape function satisfying the ZTF condition, which can be obtained with the help of the EoS (9). We will take under consideration three shape function profiles. We begin with

### \(b\left( r\right) =r_{0}\exp \left( -\mu \left( r-r_{0}\right) \right) \left( 2+r_{0}/3r\right) \)

The shape function of the Casimir wormhole is defined by

We wonder what are the effects of an additional Yukawa term on the original Casimir shape function whose profile becomes

where \(\mu \) is a positive mass scale to be identified. The original Casimir shape function can be reobtained when \(\mu =0\). The profile (21) satisfies the usual properties, namely the throat condition \(b\left( r_{0}\right) =r_{0}\), the asymptotic flatness and the flare out condition of the throat, written into the form

This is always satisfied together with the property \(1-b(r)/r>0\). Another additional property is

The energy density can be easily computed and we obtain

It is straightforward to see that, for \(\mu \rightarrow 0\), one gets the original Casimir energy density with the plates separation considered as a variable if we make the following identification

which is possible if [17]

However, the identification (26) is inconsistent with the assumption (11) because the relationship (11) leads to a vanishing redshift, while the identification (26) does not, as shown in Ref. [17] Therefore, we consider the following assumption

where *d* is the “*fixed plate distance*”. This identification fixes the scale mass \(\mu \) to the following value

Since \(\mu \ge 0\), one finds that

which is in agreement with what found in Ref. [36] but with a factor \(\sqrt{3}\) missing. Note that

while

Note also that

However, due to the relationship (29), \(\mu \rightarrow \infty \) is equivalent to \(r_{0}\rightarrow \infty \). Therefore this limiting value will be discarded. The second EFE (6) determines the value of the pressure that, differently from the Casimir wormhole, will be computed by imposing the relationship (11). A simple calculation gives

\(\omega \left( r\right) \) has the following properties

which is the original relationship between the energy density and the pressure: in this case the radial pressure. With this assumption, we get

and when

It remains to compute the transverse pressure

which has the following features, for

The SET becomes

where

and

On the throat the SET reduces to

and in the limit \(\mu \rightarrow 0\), one gets

which is verified when the relationship (29) is satisfied. Moreover the SET (46) is in agreement with the SET structure found in Ref. [17]. It is interesting to note that, for \(r\rightarrow \infty \), the SET vanishes reproducing a Minkowski SET. With an abuse of language, one can say that in this limit we find a behavior that looks like a Generalized Absurdly Benign Traversable Wormhole [36]. We say “*it looks like*” because the SET vanishes for a limiting value of the radial coordinate and not for a well determined location in space time. The next profile we are going to examine is

### \(b\left( r\right) =r_{0}\left( 2\exp \left( -\mu \left( r-r_{0}\right) \right) +r_{0}/r\right) /3\)

For the following profile

the Yukawa modification is not distributed over the whole original shape function but only on the constant term. This little displacement has an interesting consequence, because when \(\mu \rightarrow \infty \) we obtain the Ellis–Bronnikov (EB)-like wormhole [28, 29]. Indeed, the EB wormhole is

The shape function (47) satisfies the usual properties, namely the throat condition and so on. For completeness, we write the expression of the flare-out condition, which is

Even in this case, we can easily compute the energy density to obtain

which, on the throat becomes

To fix the value of \(\mu \) we adopt the same procedure of Sect. 2.1 and we find that even in this case the relationship (29) is satisfied. The pressure can be determined by solving the second EFE (6) and by imposing that the relationship (11) be satisfied, namely

This time \(\omega \left( r\right) \) has the following properties

and even in this case the original relationship between the energy density and the pressure is preserved. Thus the radial pressure is

and one finds that

The last quantity to compute is \(p_{t}(r)\), namely

which has the following features, for

To summarize the SET for this particular shape function becomes

where

and

On the throat the SET reduces to

and in the limit \(\mu \rightarrow 0\), one gets

which is in agreement with the SET structure found in Ref. [17] only for \(\mu =0\). Finally, we investigate the following shape function

### \(b\left( r\right) =r_{0}\left( 2+r_{0}\exp \left( -\mu \left( r-r_{0}\right) \right) /r\right) /3\)

For the following profile

the Yukawa modification is now put only on the variable term. Even in this modification, we have an interesting consequence, because when \(\mu \rightarrow \infty \) we obtain a constant term smaller than the throat. The shape function (47) satisfies the usual properties, namely the throat condition and so on. For completeness, we verify if the flare out condition is satisfied. We find that

and on the throat one gets

The energy density is straightforward to obtain since

and for \(\mu \rightarrow \infty \), one finds

On the throat we obtain

The second Einstein’s field equation (6) determines the value of the pressure that and, even in this case, we impose that the relationship (11) be satisfied. This implies that the redshift function vanishes and that

This time \(\omega \left( r\right) \) has the following properties

As we can see, from the relationship (76) , one finds that \(\omega \left( r\right) \) is divergent: this is a consequence of the EoS. Indeed for \(r\rightarrow \infty \), the energy density (70) vanishes because of the presence of the damping exponential overall, while into the pressure the damping exponential appears only in the constant term. For this reason, this profile will be discarded. In the next section, we explore a profile which is a superposition of the previous profiles with the aim of generalizing as much as possible the features of a Yukawa–Casimir wormhole.

## Superposing traversable wormholes shape functions

In this section we will consider a linear combination of the previous profiles described by the following shape function

with \(\mu ,\nu >0,\alpha \ge 0\) and \(c\in \mathbb {R} \). Note that when \(\alpha =2/3\), \(c=1\) and \(\mu =\nu =0\), it is immediate to see that the Casimir wormhole shape function is obtained. When \(\alpha =1\), we find a pure Yukawa wormhole discussed in Ref. [18] as well as for \(c=0\) and \(\mu =\nu \). For \(\alpha =0\) and \(c=1\), one finds the Yukawa modification to the EB wormhole. Finally, note that for \(\mu =\nu =0\) and \(c<-1\), the wormhole is no more traversable. As a first step we examine under what conditions the flare-out property is satisfied. From

one finds

From the Eq. (80) we can easily compute the energy density

and on the throat we obtain

As we can see, \(\rho \left( r_{0}\right) \) can be considered as a function of the throat. In order to fix the wormhole throat, we find the stationary point of \(\rho \left( r_{0}\right) \), assuming

Plugging \(\bar{r}_{0}\) into Eq. (83), one finds

where we have imposed that, even in this case, the source is described by Eq. (16). A solution of the previous equation is given by

Plugging \(\bar{\mu }\) into \(\bar{r}_{0}\) of Eq. (84), one finds

The value of \(\alpha \) can be determined with the help of the relationship (11) which, in this case, becomes

and on the throat one finds

\(\omega \left( r_{0}\right) \) can be further reduced to the following simple expression

where we have used Eqs. (86) and (87) and where we have imposed the Casimir relationship between pressure and energy density on the throat. Plugging (90) into (87), one gets

which is the same result of Eq. (29). With the help of the Eqs. (86), (87) and (90), it is possible to show that \(\bar{r}_{0}\) represents the minimum of \(\rho \left( r_{0}\right) \). This means that if we want to have a wormhole throat with a radius smaller than \(\bar{r}_{0}\), we need to have an increasing of negative energy. To complete the analysis, we compute the transverse pressure

which, on the throat, becomes

With the help of the relationships (91), (86) and (90), one finds

The analytic form of the SET is quite complicated. However, it becomes very simple on the throat

and in particular in correspondence of the minimum \(\bar{r}_{0}\), it simplifies to

**Remark** It is important to observe that thanks to the EoS (9), the form of the SET is

By construction the SET (97) is divergenceless, but it is not traceless. However, it is always possible to rearrange the previous SET (97) in such a way to extract the traceless part. Indeed

where \(T_{\mu \nu }^{T}\) is the traceless part of the SET (97). It is interesting to observe that by imposing the following condition

one finds that

independently on the form of \(\omega \left( r\right) \). This means that either by decomposing the SET like in Eq. (98) and fixing \(\omega \left( r_{0}\right) =1\) or by fixing \(\omega \left( r_{0}\right) =3\) without the decomposition (98), it is always possible to preserve the fundamental relationship between pressure and energy density. Note also that, from the point of view of the wormholes throat size, the choice (99) or the choice \(\omega \left( r_{0}\right) =3\), do not change the size of the throat size, as it should be. In the next section, we are going to examine the reverse procedure, namely we fix the form of the energy density and we will deduce the form of the shape function.

## Traversable wormholes with a Yukawa energy density profile

In this section we change the strategy and we fix our attention on some energy density profiles modified by a Yukawa term and with the help of Eq. (5), we will deduce the form of the shape function. Three different forms will be examined. We begin with the following profile

### \(\rho \left( r\right) =-r_{0}\rho _{C}\frac{e^{-\mu \left( r-r_{0}\right) }}{r}\)

where \(\rho _{C}\) has dimensions of an energy density and \(\mu \) is a positive mass scale parameter. We can identify \(\rho _{C}\) with the value expressed by (16). Equation (101) can be easily integrated to obtain

where we have used the condition \(b\left( r_{0}\right) =r_{0}\). It is immediate to verify that the shape function (102) satisfies the asymptotic flatness and the flare-out condition. To have ZTF, Eq. (11) must be imposed, that it means

As one can see, for \(r\rightarrow \infty \), \(\omega \left( r\right) \rightarrow \infty \). However, one can adopt another strategy to have a finite \(\omega \left( r\right) \). Indeed from the shape function (102), we can find that there exists \(\bar{r}>r_{0}\), such that \(b\left( \bar{r}\right) =0\), where

where \(\mathrm {W}\left( x\right) \) is the Lambert function defined mathematically as the multivalued inverse of the function \(x\exp \left( x\right) \),

If \(-1/e<x<0\), there are two real solutions, and thus two real branches of \(\mathrm {W}\) [33,34,35]. Inspired by the Absurdly Benign Traversable Wormhole (ABTW) and its generalization, the Generalized Absurdly Benign Traversable Wormhole (GABTW) [36], we define the shape function (102) in such a way that

where \(\bar{r}\) has been defined in (104). As a consequence also \(\omega \left( r\right) \) behaves in the same way and therefore also the radial pressure. Nevertheless the energy density does not vanish because \(\bar{r}\) does not set to zero its value. However \(\rho \left( \bar{r}\right) \) can be very small and therefore even the transverse pressure. Therefore outside the region defined by \(r\ge \bar{r}\), one obtains a quasi-Minkowski spacetime. To complete the analysis, we compute the value of \(\omega \left( r\right) \) on the throat. We find

in agreement with what found in Ref. [36] and with Eq. (19).

### \(\rho \left( r\right) =-\frac{\rho _{C}}{2}\left( \alpha +\beta r_{0}\frac{e^{-\mu \left( r-r_{0}\right) }}{r}\right) \)

The second energy density profile we are going to consider is obtained with a small modification of the profile (101)

As we can see, this is a linear combination between the original Casimir profile and the Yukawa profile (101). Note that for \(\mu =0\), \(\alpha =\beta =1\) and \(r=r_{0}\), we obtain the pure Casimir energy density. Note also that this profile is a generalization of the potential (2). Differently from the profile (101), here we can choose \(\alpha \) and \(\beta \) in such a way to have

The motivation for this choice will be clarified below. The shape function can be easily computed and we find

It is easy to see that for \(r\gg r_{0}\)

This behavior is due to the constant term in (108) which is dominant and produces a divergent shape function. However, since \(\beta \) is not fixed, we can impose that \(b\left( \bar{r}\right) =0\), where \(\bar{r}\) is the same of Eq. (109). Plugging the value of \(\alpha \) found in Eq. (109) into Eq. (110), and by imposing that \(b\left( \bar{r}\right) =0\), one finds

and

Thus if we assume that for \(r>{\bar{r}}\), \(b\left( \bar{r}\right) =0\), we get a feature similar to the ABTW. Moreover, to have ZTF, \(\omega \left( r\right) \) must be

and, with the help of Eqs. (109) and (112), we get

For \(r\rightarrow {\bar{r}}\), \(\omega \left( r\right) \) is an indeterminated form of the kind 0/0. However close to the point \(r=\bar{r}\), the shape function can be approximated by

while \(b^{\prime }\left( r\right) \) can be approximated by

Thus, even in this case, we can assume that

On the throat the analytic form of \(\omega \left( r\right) \) is far to be simple. Indeed, we find

However, one finds that

while

However when \(\bar{r}\gg r_{0}\), we get

which is finite and positive. Therefore we can conclude that from an energy density of the form (108), it is possible to extract another shape function which generalizes an ABTW. It is important to observe that such a generalization is realized because of the presence of a repulsive Yukawa–Casimir profile, otherwise for a choice of the form

the energy density outside the region \(r>\bar{r}\) is not Minkowski. Although interesting, the profile (108) has the defect of having a way to compare the throat radius with the original Casimir source, like in Eq. (107).

### \(\rho \left( r\right) =\frac{r_{0}\rho _{C}}{r}\left( \alpha e^{-\mu \left( r-r_{0}\right) }-\left( 1-\alpha \right) e^{-\nu \left( r-r_{0}\right) }\right) \)

To this purpose, we fix our attention on an energy density profile which can reproduce both a Yukawa behavior and an ABTW. This is represented by

When \(\alpha =0\), we obtain the profile (101), while when \(\alpha =1\), we obtain its repulsive version. \(\forall \alpha \ne 0,1\), we have a linear superposition of the Yukawa–Casimir profile. The combination of an attractive and repulsive potential is suggested also by the potential (108) together with the option (109) . Note that the option (109) is relevant only if one wishes to reproduce an ABTW. If such an option is not adopted the existence of a traversable wormhole is guaranteed the same. For the profile (124), it is immediate to calculate the form of the shape function

and for \(r\rightarrow \infty \), we find

which can be set to zero if

Plugging Eq. (127) into (125), one finds

which is useful to compute the function of the EoS \(\omega \left( r\right) \)

This time the function \(\omega \left( r\right) \) goes to zero for large values of *r*, while on the throat one gets

Even in this case, if we desire to extract information on the throat size, we need to compare \(\omega \left( r_{0}\right) \) with a physical source like the Casimir source. To do this, we assume that

To do calculations in practice, it is useful the following setting

and Eq. (131) becomes

whose solution is

To constraint *r* to be very small, we observe that the r.h.s. of Eq. (134) vanishes when

\(m_{-}\) will be discarded because is the negative root. Note that it is not necessary to have a vanishing *x*, rather we need an *x* with a value of the order of \(10^{-10}\) or greater. This is due to the rescaling in (132) setting the size of the wormhole throat to be of the order of

Therefore we can conclude that with a linear combination of two Yukawa–Casimir profiles, actually a difference of them, one finds a traversable wormhole with a throat that can be fine tuned with respect to the original Casimir source.

## Conclusions

In this paper we have taken under examination the modification of the Casimir wormhole examined in Ref. [17] which uses, as a source, the negative energy density of the Casimir device. Differently from Ref. [17], this time we have imposed the ZTF condition to obtain different solutions. We have found that the ZTF condition can be imposed only if we modify the form of the energy density or the form of the shape function. To this purpose we have considered Yukawa type modifications of the original profile. The motivation for this choice stands in the attempt to detect signals of variations of the ordinary gravitational field even for TW. In this context, it appears interesting to note that, recently, investigations on how it is possible to see signals of a TW have been considered, especially in the framework of extended theories of gravitation [37,38,39,40]. To further proceed, we have divided the paper in two parts: the first part is devoted to the analysis of the modification of the shape function with a Yukawa term and the second part is devoted to the modification of the original Casimir energy density with an appropriate Yukawa term. In the first part we have examined three different profiles having in common the Casimir wormhole shape function, namely

and we have included a Yukawa modification of the type \(\exp \left( -\mu \left( r-r_{0}\right) \right) \) acting on every single term and globally. Two of the three profiles have shown features compatible with the throat size estimated in Ref. [36], while one of them has developed a divergent behavior for large values of the radial variable *r*. I recall to the reader that, in this paper, we have examined the Casimir energy density with the plates separation considered as a parameter and not as a variable. This choice has led to have a huge throat size instead of a Planckian one like in Ref. [17]. As a further analysis, we have also considered a superposition of different categories of TW. Even in this case the resulting size of the wormholes throat is huge and compatible with the size of a GABTW described in Ref. [36]. I recall again to the reader that the huge size of the wormhole throat has been found by imposing that the inhomogeneous function of the EoS (11) at the throat has a constant value compatible with the ordinary Casimir relationship \(p=3\rho \). In the second part, we have fixed the form of the energy density and we have deduced the form of the shape function with the help of the first of the EFE (5). Even in this case, we have analyzed three different profiles. Since every of these profiles has produced a correction to the size of the throat at infinity, we have considered the possibility of taken another generalization for the ABTW. In particular we have found a value of the radial variable, located at \(r=\bar{r}\) where the shape function vanishes and we have truncated the region outside \(r=\bar{r}\). In this way the pressure outside the region \(r=\bar{r}\) vanishes. However, for the profile of Sect. 4.1, the energy density does not vanish for \(r\ge \bar{r}\), it is small because of the exponential but not nought. Therefore the structure of an ABTW or GABTW cannot be reproduced. On the contrary, for the profiles discussed in Sects. 4.2 and 4.3, it is possible to reproduce an ABTW in a generalized form different by the GABTW at the price of introducing a repulsive potential, namely we have the difference of two Yukawa profiles: one attractive and one repulsive. I recall the reader that an ABTW is defined by the following shape and redshift functions

Therefore outside the location \(r=r_{0}+a\), the spacetime is Minkowski. The same happens for a GABTW, where one finds

Of course \(1/\mu \) plays the rôle of *a* and viceversa. Note that it is the exponent in Eq. (138) and in Eq. (139) that plays a key rôle to determine the Minkowski structure for \(r\ge r_{0}+a\) or \(r\ge r_{0}+1/\mu \). Such a property is completely absent for the profiles discussed in Sects. 4.1, 4.2 and 4.3 and one must build a profile that potentially can develop such a property. This is the reason why a repulsive Yukawa profile is necessary to have a vanishing value outside a certain region. Of course this is related to the attempt to reproduce some of the features of an ABTW. If one abandons this request, the repulsive potential is not fundamental. However, our insistence to reproduce the features of an ABTW is justified by the fact the negative energy density is concentrated in a very small region of the space and there is no redshift. Coming back to the profile (124) in Sect. 4.3, it represents again a difference of two Yukawa profiles, and behaves “*like*” a ABTW, because for a sufficiently large values of *r*, \(\rho \left( r\right) \), \(p_{r}\left( r\right) \) and \(p_{t}\left( r\right) \) vanish. Nevertheless, because of the exponentials, it is not necessary that the radial value *r* needs to be really large. Another interesting feature of the profile (124) is that, this time, we can fine tune the throat size down to acceptable values, which is exactly what one needs. Always on the side of phantom energy I proposed the idea of *Self-Sustained Traversable Wormholes*, namely TW sustained by their own quantum fluctuations [41,42,43,44,45]. Even in this case, because the quantum fluctuation carried by the graviton behaves like the ordinary Casimir effect, we found that no need for phantom contribution is necessary. On this context, in a next paper we will explore how behaves a system formed by the *Casimir TW*, here analyzed, and the corresponding self-sustained TW version. To conclude, we have also to point out that in the context of *Self-Sustained Traversable Wormholes*, namely TW sustained by their own quantum fluctuations [41,42,43,44,45], could be interesting to consider how the Casimir source and the quantum fluctuation carried by the graviton combine to see if the GABTW can be self-sustained in this context. In this picture, the Casimir source could be interpreted as the switch on to power the traversability of the wormhole.

## Data Availability Statement

This manuscript has no associated data or the data will not be deposited. [Authors’ comment: This is a theoretical study and no experimental data has been listed.]

## Notes

- 1.
However, if \(\rho \left( r\right) \) represents the mass density, then we have to replace \(\rho \left( r\right) \) with \(\rho \left( r\right) c^{2}.\)

## References

- 1.
H. Yukawa, Proc. Phys. Math. Soc. Jpn.

**17**, 48 (1935) - 2.
D. Borka, P. Jovanović, V. Borka Jovanović, Constraining the range of Yukawa gravity interaction from S2 star orbits. JCAP

**11**, 050 (2013). arXiv:1311.1404 [astro-ph.GA] - 3.
A.F. Zakharov, P. Jovanovic, D. Borka, V. Borka Jovanovic, Constraining the range of Yukawa gravity interaction from S2 star orbits II: bounds on graviton mass. JCAP

**05**, 045 (2016). arXiv:1605.00913 [gr-qc] - 4.
A.F. Zakharov, P. Jovanovic, D. Borka, V. Borka Jovanovic, Constraining the range of Yukawa gravity interaction from S2 star orbits III: improvement expectations for graviton mass bounds. JCAP

**04**, 050 (2018). arXiv:1801.04679 [gr-qc] - 5.
S. Capozziello, V. Borka Jovanović, D. Borka , P. Jovanović, Phys. Dark Univ. 100573 arXiv:2004.11557 [gr-qc]

- 6.
I. De Martino, R. Lazkoz, M. De Laurentis, Analysis of the Yukawa gravitational potential in f(R) gravity I: semiclassical periastron advance. Phys. Rev. D

**97**(10), 104067 (2018). arXiv:1801.08135 [gr-qc] - 7.
M. De Laurentis, I. De Martino, R. Lazkoz, Analysis of the Yukawa gravitational potential in f(R) gravity II: relativistic periastron advance. Phys. Rev. D

**97**(10), 104068 (2018). arXiv:1801.08136 [gr-qc] - 8.
S. Capozziello, M. De Laurentis, Extended theories of gravity. Phys. Rep.

**509**, 167 (2011). arXiv:1108.6266 [gr-qc] - 9.
S. Capozziello, C. Altucci, F. Bajardi, A. Basti, N. Beverini, G. Carelli, D. Ciampini, A. D. V. Di Virgilio, F. Fuso, U. Giacomelli, E. Maccioni, P. Marsili, A. Ortolan, A. Porzio, A. Simonelli, G. Terreni , R. Velotta, Constraining theories of gravity by GINGER experiment. Eur. Phys. J. Plus

**136**(4), 394 (2021) [Erratum: Eur. Phys. J. Plus**136**(5), 563 (2021)] arXiv:2103.15135 [gr-qc] - 10.
J.W. Moffat, Black Holes in Modified Gravity (MOG). Eur. Phys. J. C

**75**, 175 (2015). arXiv:1412.5424 [gr-qc] - 11.
Z. Berezhiani, F. Nesti, L. Pilo, N. Rossi, Gravity modification with Yukawa-type potential: dark matter and mirror gravity. JHEP

**0907**, 083 (2009). arXiv:0902.0144 [hep-th] - 12.
M. Bordag, G.T. Gillies, V.M. Mostepanenko, New constraints on the Yukawa-type hypothetical interaction from the recent Casimir force measurement. Phys. Rev. D

**56**, 6 (1997) - 13.
M. Bordag, B. Geyer, G.L. Klimchitskaya, V.M. Mostepanenko, Stronger constraints for nanometer scale Yukawa-type hypothetical interactions from the new measurement of the Casimir force. Phys. Rev. D

**60**, 055004 (1999). arXiv:9902456 [hep-ph] - 14.
B.W. Ninham, M. Boström, C. Persson, I. Brevik, S.Y. Buhmann, Bo. E. Sernelius, Casimir forces in a plasma: possible connections to Yukawa potentials. Eur. Phys. J.

**D 68**, 328 (2014). arXiv:1409.1032 [quant-ph] - 15.
M.S. Morris, K.S. Thorne, Am. J. Phys.

**56**, 395 (1988) - 16.
M. Visser,

*Lorentzian Wormholes: From Einstein to Hawking*(American Institute of Physics, New York, 1995) - 17.
R. Garattini, Casimir Wormholes. Eur. Phys. J.

**C 79**(11), 951 (2019). arXiv:1907.03623 [gr-qc] - 18.
R. Garattini, Traversable Wormholes and Yukawa Potentials. arXiv:1907.03622 [gr-qc]

- 19.
H. Casimir, Proc. Kon. Ned. Akad. Wetenschap

**51**, 793 (1948) - 20.
M. Sparnaay, Attractive forces between flat plates. Nature

**180**, 334 (1957) - 21.
M. Sparnaay, Measurement of attractive forces between flat plates. Physica

**24**, 751 (1958) - 22.
S. Lamoreaux, Demonstration of the Casimir force in the 0.6 to 6 mm range. Phys. Rev. Lett.

**78**, 5 (1997) - 23.
P.W. Milonni,

*The Quantum Vacuum: An Introduction to Quantum Electrodynamics*(Academic Press Inc., New York, p. 99) (1994) - 24.
S.W. Kim, arXiv:1302.3337v1 [gr-qc]

- 25.
P. Martin-Moruno, M. Visser, JHEP

**1309**, 050 (2013). arXiv:1306.2076 [gr-qc] - 26.
M. Bouhmadi-Lopez, F.S.N. Lobo, P. Martin-Moruno, JCAP

**1411**, 007 (2014). arXiv:1407.7758 [gr-qc] - 27.
M. Visser, S. Kar, N. Dadhich, Phys. Rev. Lett.

**90**, 201102 (2003) - 28.
H. Ellis, Ether flow through a drainhole—a particle model in general relativity. J. Math. Phys.

**14**, 104 (1973) - 29.
K.A. Bronnikov, Acta Phys. Pol. B

**4**, 251 (1973) - 30.
P. Mishra, T.P. Singh, Phys. Rev. D

**88**, 104036 (2013). arXiv:1108.2375 [astro-ph.GA] - 31.
K.K. Nandi, Y.Z. Zhang, R.G. Cai, A. Panchenko, Phys. Rev. D

**79**, 024011 (2009). arXiv:0809.4143 [gr-qc] - 32.
O. Genç, Some aspects of Morris–Thorne wormhole in scalar tensor theory. arXiv:1709.06431 [gr-qc]

- 33.
J.H. Lambert, Acta Helv.

**3**, 128 (1758) - 34.
J.H. Lambert, in Nouveaux M émoires de l’Acad émie Royale des Sciences et Belles-Lettres de Berlin (German Academy of Sciences Berlin, 1772)

- 35.
D. Veberic, Comput. Phys. Commun.

**183**, 2622 (2012) - 36.
R.Garattini, Generalized absurdly benign traversable wormholes powered by Casimir energy. Eur. Phys. J. C

**80**(12), 1172 (2020). arXiv:2008.05901 [gr-qc] - 37.
V. De Falco, M. De Laurentis, S. Capozziello, Epicyclic frequencies in static and spherically symmetric wormhole geometries. Phys. Rev. D

**104**(2), 024053 (2021). arXiv: 2106.12564 [gr-qc] - 38.
V. De Falco, E. Battista, S. Capozziello, M. De Laurentis, Reconstructing wormhole solutions in curvature based Extended Theories of Gravity. Eur. Phys. J. C

**81**(2), 157 (2021). arXiv: 2102.01123 [gr-qc] - 39.
V. De Falco, E. Battista, S. Capozziello, M. De Laurentis, Testing wormhole solutions in extended gravity through the Poynting–Robertson effect. Phys. Rev. D

**103**(4), 044007 (2021). arXiv:2101.04960 [gr-qc] - 40.
V. De Falco, E. Battista, S. Capozziello, M. De Laurentis, General relativistic Poynting-Robertson effect to diagnose wormholes existence: static and spherically symmetric case. Phys. Rev. D

**101**(10), 104037 (2020). arXiv:2004.14849 [gr-qc] - 41.
R. Garattini, Class. Quantum Gravity

**22**, 1105 (2005). arXiv:gr-qc/0501105 - 42.
R. Garattini, Class. Quantum Gravity

**24**, 1189 (2007). arXiv:gr-qc/0701019 - 43.
R. Garattini, F.S.N. Lobo, Class. Quantum Gravity

**24**, 2401 (2007). arXiv:0701020 [gr-qc] - 44.
R. Garattini, F.S.N. Lobo, Phys. Lett. B

**671**, 146 (2009). arXiv:0811.0919 [gr-qc] - 45.
R. Garattini, F.S.N. Lobo, Phys. Rev. D

**85**, 024043 (2012). arXiv:1111.5729 [gr-qc]

## Acknowledgements

This work has been supported by the project “ Traversable Wormholes: A Road To Interstellar Exploration” of Texas A & M Engineering Experiment Station (TEES), The Texas A & M University System.

## Author information

### Affiliations

### Corresponding author

## Appendix A: Features of the superposition of traversable wormholes

### Appendix A: Features of the superposition of traversable wormholes

In this section we are going to explore some of the features of the profile ( 79) which is quite general to include the other profiles discussed in Sect. 2 of this paper. We begin to examine the proper length which is defined as

where the “±” depends on the wormhole side we are. In the case of the shape function (79), we know that it vanishes exponentially for \(r\rightarrow +\infty \). This is true for \(\mu \ne 0\) and \(\nu \ne 0\). For instance for \(\mu =\nu =0\) and \(c=1\), one finds that the shape function is represented by the EB wormhole (48) , whose proper length is

For this reason, for the other cases it is sufficient to consider what happens close to the throat. In general, we can write

Thus the approximated proper length becomes

A further approximation leads to

and, in the case of the shape function (79) , (A6) assumes the form

The argument of the denominator is positive if and only if the flare-out condition is satisfied. However, we can use the constraint (86) and (90) to have a better estimate. We find

or by means of the constraint (87), we can write the proper length only in terms of the plates separation

In a similar way, to compute the embedded surface, we need to evaluate

which, close to the throat, becomes

To further investigate the properties of the Traversable Wormholes described by the shape function (79), we consider the computation of the total gravitational energy for a wormhole, defined as [31]

where

is the total mass and

is the proper mass. In particular we find for the total mass

where we have used the relationship (91) and we have momentarily reintroduced the speed of light. For \(M_{\pm }^{P}\), we can estimate the value of the integral close to the throat, following what has been done for the proper length and the embedded surface. We can write

where the “±” depends one the wormhole side we are. Thus the total gravitational energy (A12) becomes

Even for the total energy, this is true for \(\alpha \ne 0\), \(\nu \ne 0\) and \(c\ne 1\). Indeed for \(\alpha =\nu =0\) and \(c=1\) we can write the total gravitational energy of the EB wormhole (48) which reduces to

It is interesting to note that the total energy is concentrated completely on the throat and at infinity vanishes showing a screening mechanism: in other words, the “*imprint at infinity*” disappears [16]. Another important traversability condition is that the acceleration felt by the traveller should not exceed Earth’s gravity \(g_{\oplus }\simeq 980\) \(cm/s^{2}\). In an orthonormal basis of the traveller’s proper reference frame, we can find

and in this case, because \(\Phi (r)=0\), the traveller has no acceleration, which is in agreement with Ref. [15]. As regards the lateral tidal forces, we find

where we have used the relationship (11). This is a constraint about the velocity with which observers traverse the wormhole. \(\eta \) represents the size of the traveller which can be fixed approximately equal, at the symbolic value of 2 *m* [15]. If we assume a constant speed *v* and \(\gamma \simeq 1\), close to the throat, the lateral tidal constraint becomes

If the observer has a vanishing *v*, then the tidal forces are null. Note that the total time defined by

and the total proper time given by

are the same for the profile (79) because the redshift function is nought. Assuming that *v* is approximately constant, we can use the estimate (A21) to complete the evaluation of the crossing time which approximately is

However, we can use the constraint (84), (86), (87) and (90) to have a better estimate of the crossing time which becomes

where we have considered a possible time trip in going from one station located in the lower universe, say at \(l=-l_{1}\), and ending up in the upper universe station, say at \(l=l_{2}\). Following Ref. [15], we have located \(l_{1}\) and \(l_{2}\) at a value of the radius such that

that it means \(1-b\left( r\right) /r\simeq 1\) which is in agreement with the estimates found in Ref. [15]. The last property we are going to discuss is the “ total amount” of ANEC violating matter in the spacetime [27] which is described by Eq. ( A27). For the metric (79), one obtains

where we have approximated the expression close to the throat and where we have defined

After the integration, we find

and the result is finite. Therefore we can conclude that, in proximity of the throat the ANEC can be arbitrarily small as it should be.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Funded by SCOAP^{3}

## About this article

### Cite this article

Garattini, R. Yukawa–Casimir wormholes.
*Eur. Phys. J. C* **81, **824 (2021). https://doi.org/10.1140/epjc/s10052-021-09634-3

Received:

Accepted:

Published: