Traversable Morris–Thorne–Buchdahl wormholes in quadratic gravity

The special Buchdahl-inspired metric obtained in a recent paper [Phys. Rev. D 107, 104008 (2023)] describes asymptotically flat spacetimes in pure R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}^{2}$$\end{document} gravity. The metric depends on a new (Buchdahl) parameter k~\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{k}$$\end{document} of higher-derivative characteristic, and recovers the Schwarzschild metric when k~=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{k}=0$$\end{document}. It is shown that the special Buchdahl-inspired metric supports a two-way traversable Morris–Thorne wormhole for k~∈(-1,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{k}\in (-1,0)$$\end{document} in which case the Weak Energy Condition is formally violated, a naked singularity for k~∈(-∞,-1)∪(0,+∞)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{k}\in (-\infty ,-1)\cup (0,+\infty )$$\end{document}, and a non-Schwarzschild structure for k~=-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{k}=-1$$\end{document}.

All we know about exotic matter is that a) it violates our perception of energy, that is, an observer may measure some negative amount of rest energy density and b) it has the ability to sustain wormholes.So far there has been no general theory about exotic matter and most of-if not all-the wormhole solutions derived so far were obtained geometrically upon running the field equations of general relativity (GR) from left to right: The energymomentum tensor (EMT) derived this way was exotic.In this work we rather run the field equations of purely quadratic gravity in the usual way, that is, by solving the field equations.
In generalized and modified theories of GR, explicit EMT are not usually added to the field equations.Instead, extra terms or corrections are introduced to the gravitation sector, which can play the role of exotic matter without truly being exotic matter.An example is the Brans-Dicke (BD) action, d 4 x √ −g φ R − ω φ ∂ µ φ∂ µ φ , which allows the formation of wormholes [25][26][27][28], where the scalar field φ acts as an exotic form of matter that violates the WEC.Another interesting case is the family of f (R) gravity, which in general can be formally cast as a scalar-tensor theory but the scalar field is directly associated with R without its own dynamics.As such, the scalar field does not truly represent an exotic form of matter.In the case of pure R 2 gravity, the scalar field is trivially identical to R.
Violation of the energy conditions, and particularly of the WEC, has direct astrophysical consequences.In [29] it was shown how the non-violation of the WEC constrains the time evolution of both the Hubble parameter and coordinate distance and sets an upper bound for Ω M .Similar conclusions were derived in [30].In the other case, when the WEC is violated, no such constraints emerge.Assuming that the SEC and the DEC hold, a first proof of the Cosmic No-Hair Theorem was given in [31].Another theorem on the future fate of a spacetime, the Lorentzian Splitting Theorem, was proven upon admitting the SEC [32].From a quantum-mechanical point of view, all local energy conditions are violated by quantum fields and also by some classical fields as the non-minimally coupled scalar fields [33,34] and the future-eternal inflating spacetimes [35].However, the scale of violation can be minimized for some cut-andpaste geometric constructions [36], and for type I wormholes without the cut-and-paste construction [5] where the extent of exotic matter has been shown to be inversely proportional to the square of the mass of the wormhole.A detailed compte-rendu of the consequences and plausible violations of the energy conditions are reported in [37], instances include the possibility of formation of cosmological singularities in spatially open or flat spacetimes if the WEC is observed and the possibility of superluminal motion (wrap drive and traversable wormholes) if the WEC is violated.
The pure R 2 theory was first proposed by Buchdahl in the early 1960's [38] and has recently experienced a revival of interest, with active investigations into its theoretical properties as well as its implications for black holes and cosmology [39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54][55][56][57].Among the various extensions of GR, the pure R 2 theory enjoys several distinct advantages.It is a parsimonious theory, having only one term in the action, (2κ where the grav-itational coupling κ is a dimensionless parameter.It is the only theory that is both scale invariant and ghostfree.Regarding the former, it actually possesses a restricted scale invariance under a Weyl transformation, g µν → Ω 2 (x)g µν , where Ω(x) obeys a harmonic condition, Ω = 0, as discovered in [58].Regarding the latter, as a member of f (R), its scalar degree of freedom involves only second-order derivatives when transitioning from the Jordan frame to the Einstein frame, thereby evading the Ostrogradsky instability that often plagues higherderivative gravity [59].Furthermore, pure R 2 gravity has been shown propagate two massless modes: a spin-2 tensor mode and a spin-0 scalar mode, each carrying it own significance [60].On one hand, the massless spin-2 tensor mode indicates the emergence of a long-range potential with the correct Newtonian tail ∼ 1/r [61].On the other hand, the massless spin-0 scalar mode could potentially be responsible for an additional long-range potential, thereby introducing new physics.
One concrete realization of new physics in pure R 2 gravity manifests through the work of Buchdahl in 1962, in which he originated a program aimed at finding vacuum solutions for the theory [38].He was able to make significant progress with his efforts boiling down to solving a non-linear second-order ordinary differential equation (ODE).If an analytical solution to his ODE could be found, then the vacuo solutions he sought would automatically ensue.Unfortunately, Buchdahl deemed the ODE insoluble, prompting him to suspend further pursuit.Consequently, his groundbreaking paper has remained relatively obscure within the gravitational research community for the past sixty years.However, recent advancements made by one of us have revitalized Buchdahl's program and brought it to fruition.Section II in this paper will review its final outcome.
Significantly, the Buchdahl-inspired solutions exhibit non-constant scalar curvature, a distinctive feature resulting from the fourth-derivative nature of the theory.This non-constant scalar curvature is controlled by a new parameter known as the Buchdahl parameter k.Remarkably, these solutions defy the generalized Lichnerowicz "theorem" proposed in [39][40][41], which stipulates that static vacuum solutions of pure R 2 gravity must possess constant scalar curvature exclusively.The Buchdahlinspired solutions evade this "theorem" by circumventing one of the central assumptions [54].The non-constant scalar curvature observed in these solutions is a manifestation of higher-derivative effects, which are encapsulated by the Buchdahl parameter k.
In this paper, we shall use the closed analytical vacuum solution for pure R 2 gravity derived in Ref. [53] to show the formation of a wormhole that connects two asymptotically flat spacetime sheets via a "throat".This wormhole is enabled by the high-derivative nature of the theory without requiring complicated ingredients or true exotic matter.
The paper is structured as follows.In Sec.II we review Buchdahl-inspired metrics obtained in Refs.[52,53,55]; in Sec.III we present two additional representations for Buchdahl-inspired metrics that are asymptotically flat; in Sec.IV we map the special (asymptotically flat) Buchdahl-inspired metrics to the Morris-Thorne ansatz, investigate their properties, and construct a wormhole when the Weak Energy Condition is violated.
In Ref. [52] we advanced a program initiated by Buchdahl in 1962 seeking vacuo configurations for pure R 2 gravity [38].The field equation in vacuo has a static spherisymmetric solution which is expressible in terms of two auxiliary functions p(r) and q(r) per (2) The two functions p and q are coupled via a pair of firstorder evolution-type ODE's: Reflecting the fourth-order nature of quadratic gravity, this solution is specified by four parameters: Λ, k, p 0 and q 0 .When k = 0 the evolution rules recover the Schwarzschild-de Sitter.When k = 0, the Ricci scalar is non-constant and is given by R(r) = 4Λ e −k dr r q(r) (5) It approaches 4Λ at spatial infinity, indicating an asymptotic de Sitter behavior.
In Ref. [53] we further advanced the solutions for the case Λ = 0 and obtained an exact closed analytical form for an asymptotically flat non-Schwarzschild metric, which was called the special Buchdahl-inspired metric, expressible as1 in which k := k rs and ζ := 1 + 3 k2 .It contains two parameters, r s playing the role of a Schwarzschild radius, and k a new (Buchdahl) dimensionless parameter.The solution holds for all value of r ∈ R except at r = 0 and r = r s .The radial direction thus comprises of three sections: 1.The "exterior", r > r s , 2. The "interior", 0 < r < r s , 3. The "repulsive" gravity domain, r < 0. We exclude this unphysical region from our consideration.
Note that the two components g tt and g rr flip their signs at the interior-exterior boundary, r = r s .The Kruskal-Szekeres diagram is analytically constructed in Ref. [53].
Although the special Buchdahl-inspired metric, Eqs. ( 6) and ( 7), is Ricci-scalar flat, viz.R = 0, it is not Ricci flat, hence non-Schwarzschild.Moreover, it can be verified that [55] and (upon taking the trace) That is to say, the solution formally obeys the following equation with the non-vanishing term in the right hand side acting as a "quasi" energy-momentum tensor (EMT) and making the solution non-Schwarzschild.This "quasi" EMT is thus a surrogate of exotic matter which would sustain a wormhole under certain circumstances to be explored in this paper.

A. The isotropic coordinates
For the "exterior" section, let us choose a variable r and a function g(r) to fulfill two requirements: With ρ(r) given in Eq. ( 7), solving them or and rendering the metric It is straightforward to see that these expressions are unchanged upon the "image reflection" 4r ζr s ζr s 4r (21) This means that the two separate domains r > ζrs 4 and r < ζrs 4 are reciprocal images, with the value ζr s /4 being a "reflection point".Also, the case k = 0 (viz.ζ = 1) recovers the Schwarzschild metric in Weyl's isotropic coordinates:

Kretschmann invariant
The Kretschmann scalar K := R µνρσ R µνρσ is given in [53] and we shall not reproduce it here.We only report its expression for the isotropic coordinates which, by design, only cover the "exterior" section Remark 1.Interestingly, the "non-analytic" piece in K is isolated in whereas for k = −1 (ζ = 2) In both cases, the Kretschmann scalar carries the same function form.

B. Another representation
Let us define a new radial coordinate r such that Then and The metric given in ( 6)-( 7) can be brought into This representation brings the special Buchdahl-inspired metric under the umbrella of the generalized Campanelli-Lousto solution in Brans-Dicke gravity that we uncover in another report [62].

IV. MORRIS-THORNE-BUCHDAHL WORMHOLES
In terms of x := 1 − rs r ∈ R, the special Buchdahlinspired metric becomes The areal radius in the class is For the exterior: which would have an acceptable root if k ∈ (−1, 0).For the interior: which would have an acceptable root The behavior of R as a function of r is shown in 1. Panel (B) is representative of k ∈ (−1, 0) exhibits a minimum for R(r) in the exterior.We shall bring the metric above to the Morris-Thorne ansatz [2] If we focus on the "exterior" region alone, then (x ∈ (0, 1)) 37) For the exterior, viz.x ∈ (0, 1), let us make a further coordinate transformation In summary, the areal radius is defined as The redshift function and the shape function are given by, respectively, for the region y y * where corresponds to x ext in Eq. ( 33).
Constraint #2.-Minimum value of the R-coordinate, i.e. at the throat of the wormhole, R * being the minimum value of R, given in Eq. (48).space.The equality sign holds only at the throat.This is required in order to ensure the finiteness of the proper radial distance l(R) given in (55) where the ± signs refer to the two asymptotically flat regions which are connected by the wormhole.Note that the condition b(R)/R ≤ 1 assures that the metric component g RR does not change its sign for any R ≥ R * .

Embedding
With y * defined in Eq. ( 46), the shape function in Eq. ( 43) is In the embedding diagram, the MT ansatz Obviously dz/dR diverges at y = y * , meaning that in the embedding diagram (the lower panel of Fig. 2), z(R) is vertical at y = y * ,the "throat" of the wormhole.The function z(R) is a combination of Appell hypergeometric functions that is not particularly illuminating and hence will not be produced here. 3 Nevertheless, the proper radial distance is simpler to obtain: As an example, the upper panel of Fig. 2 plots the proper radial distance for k = −0.5, viz.B ≈ −1.134, y * ≈ 0.0628, r * ≈ 1.14 r s , R * ≈ 1.7 r s .A "throat" is manifest at R = R * .
In the range of k ∈ (−∞, −1) ∪ (0, +∞), The Kretschmann invariant diverges at y = 0 (i.e., r = r s ), indicating a physical singularity on the interior-exterior boundary, y = 0.As a result, the spacetime is not geodesically complete, and the geodesics terminate at the physical singularity.The ζ−Kruskal-Szekeres (KS) diagram previously constructed in Ref. [53] is reproduced in Fig. 3 here for the reader's convenience.In the left panel of Fig. 3, we also show the radial infalling motion of a massive particle along the trajectory A → B, with the particle eventually hitting the physical singularity, represented by point B on the interior-exterior boundary.
In the range of k ∈ (−1, 0), the Kretschmann invariant likewise diverges at y = 0 (i.e., r = r s ), indicating a physical singularity on the interior-exterior boundary, y = 0.However, since the areal radius possesses a minimum value at y * = B+1 B−1 ∈ (0, 1), we can generate a wormhole solution by "gluing" the region y * ≤ y < 1 (corresponding to r * ≤ r < +∞ as shown in Panel (B) of  diagram.The right panel of Fig. 3 shows how this is done.The wormhole "throat" is represented by the two red lines, r = r * , which further split Region (I) into (Ia) and (Ib), and Region (III) into (IIIa) and (IIIb).That is to say, we are connecting an asymptotically flat exterior sheet (i.e., Region (Ia)) with another asymptotically flat exterior sheet (i.e., Region (IIIa)), which are mirror images of each another with respect to a sign flip (T, X) ↔ (−T, −X) in their ζ−KS coordinates.
In this construction, both the upper and lower sheets of the wormhole correspond to r * ≤ r < +∞ and they approach asymptotic flatness at spatial infinity.The two sheets are smoothly connected at the "throat" at y = y * , where z(R) becomes vertical, hence ensuring a smooth connection.The upper and lower sheets are distinguished by the ± sign in the proper length parameter l, per Eq. ( 55), which runs continuously from −∞ to ∞ as a traveler moves from the lower sheet of the wormhole to the upper one.
In the right panel of Fig. 3, the radial infalling motion of a massive particle is depicted by the trajectory A → B in Region (Ia), with B lying on the "throat".The particle then emerges at point C (which also lies on the "throat" and is opposite to point B on the ζ−KS diagram) then continue on the path C → D in Region (IIIa).Note that points B and C represent the same spacetime event.

Violation of the Weak Energy Condition
Formally, the geometric form for the Weak Energy Condition requires that G µν t µ t ν ≥ 0 for every futurepointing timelike vector t µ ; e.g., see Ref. [63].In particular, G 00 ≥ 0.
The special Buchdahl-inspired metric has the 00−component of the Einstein tensor For k ∈ (−1, 0), the exterior region exhibits a wormhole "throat" as can be seen in Panel (B) of Fig. 1.At the same time, G 00 < 0 ∀r for k ∈ (−1, 0), thus violating the Weak Energy Condition.
As Morris and Thorne envisioned [1], in a traversable wormhole, light rays that enter it at one mouth then reemerge at its other mouth have a cross-sectional area initially decreasing and then increasing.In order for this phenomenon to occur, there necessarily be some "gravitational repulsion" near the "throat", exerting influence on the light rays.In Eq. ( 56), the magnitude of G 00 monotonically increases as one approaches the interiorexterior boundary, r = r s .For k ∈ (−1, 0), the negative and dominant G 00 component thus acts like gravitational repulsion.We expect that it could leave signatures, distinguishing a wormhole from a black hole [64].
Simultaneously, the interior-exterior boundary is a set of naked singularities.Consequently, the spacetime for k ∈ (−1, 0) accommodates both a wormhole connecting two asymptotically flat exterior sheets (viz.Regions (Ia) and (IIIa) in Fig. 3) and a set of naked singularities (belonging to Regions (Ib) and (IIIb)), thereby representing a non-trivial geometrical structure in this situation.
As mentioned in the concluding remark of Section II, the right hand side of Eq. ( 10) corresponds to a "quasi" energy-momentum tensor, defined as (modulo a multiplicative constant) It was shown in Ref. [55] that despite the vanishing Ricci scalar throughout spacetime, this "quasi" energymomentum tensor remains well-defined and is identical to the G 00 component presented in Eq. ( 56).It effectively serves as a surrogate to exotic matter required to sustain a wormhole.

V. CONCLUSION
In a previous work [53], we derived a special Buchdahlinspired metric that describes asymptotically flat spacetimes in pure R 2 gravity.This metric, expressed in a closed analytical form, enabled us to construct a Kruskal-Szekeres diagram representing the maximal analytic extension for the metric.The Buchdahl parameter k in the metric is a new parameter that reflects the higherderivative nature of the pure R 2 action.
In this paper, we present several additional advancements.Firstly, we describe two additional representations of the metric.Secondly, we examine the metric within the framework of the Morris-Thorne ansatz.For values of k falling within the ranges (−∞, −1) and (0, +∞), the interior-exterior boundary constitutes a naked singularity.However, for the range k ∈ (−1, 0), the areal radius has a minimum value in exterior region.Despite the geodesic incompleteness of the solution in this situation, where geodesics terminate on the singularity, it is possible to shield the singularity by removing the region of space neighboring the singularity and gluing the Kruskal-Szekeres symmetric copy of the remaining space region to that same region.The Morris-Thorne-Buchdahl wormhole constructed this way consists of a pair of asymptotically flat spacetime sheets connected at a "throat" that allows two-way passage.
Thirdly, we find that when a wormhole is formed, that is, when k ∈ (−1, 0), the Weak Energy Condition is formally violated, even though no exotic matter is in presence.Therefore, pure R 2 theory can support a wormhole without the need for truly exotic matter in the energy-momentum tensor or complicated ingredients in the gravitation sector such as torsion, non-metricity, or non-locality.Our work opens up new avenues for exploring the fascinating properties of wormholes and naked singularities in higher-derivative gravity theories.For wormholes and naked singularities of quadratic relativity, all known potential astrophysical observations, including light deflection, precession, shadow, and quasi-periodic oscillations are the subjects of a future investigation [65] and cannot be carried out in this work.
Another important question, even more important than all that has been said, is the stability of the wormhole discussed in this work.Stability analysis is a more involved issue [66][67][68][69][70] as this necessitates to perform a perturbation analysis of the metric, which will make the subject of another subsequent paper.However, based on the generic analysis made in [71] the wormhole is likely to be stable.

Figure 3 :
Figure 3: ζ−Kruskal-Szekeres diagrams for the asymptotically flat Buchdahl-inspired spacetimes for k ∈ (−∞, −1) ∪ (0, +∞) (the case of naked singularity, shown in left panel) and for k ∈ (−1, 0) (the case of wormhole, shown in right panel).Left panel: On the infalling radial timelike trajectory (blue line), a particle in Region (I) eventually hits the naked singularity, r = rs.Right panel: The wormhole "throat" is depicted by the red lines, which further split Region (I) into (Ia) and (Ib), and Region (III) into (IIIa) and (IIIb).On the radial trajectory A → B → C → D, a particle in Region (Ia) first enters the wormhole mouth at point B (hence, an infalling motion) then escape into Region (IIIa) by emerging at the other mouth at point C (hence, on outgoing motion).As the two red lines are "glued" together to form a wormhole that connects Region (Ia) and Region (IIIa), the two opposite points B and C represent the same spacetime event.