One of the useful lessons of general relativity is that the Einstein equations are not, by themselves, especially conservative about the kinds of geometries they permit. Smooth Lorentzian metrics can describe black holes, gravitational waves, expanding cosmologies, singularity formation, and even rather exotic global structures. Traversable wormholes are a particularly clean example of this phenomenon. As a matter of differential geometry, it is not difficult to write down a metric containing a throat joining two regions of spacetime. The difficulty appears only when one asks what stress-energy tensor is required to support such a geometry.
In ordinary Einstein gravity, the answer is severe: a traversable throat requires matter which violates the null energy condition. Roughly speaking, the throat wants to collapse inward, and the only way to keep it open is to supply a sufficiently negative radial pressure. This is the point at which the geometry stops being merely a clever construction and becomes a question about the physical matter content of the theory.
The standard Morris-Thorne wormhole metric is
$$ ds^2 = -e^{2\Phi(r)}dt^2 + \frac{dr^2}{1-b(r)/r} + r^2(d\theta^2+\sin^2\theta\,d\phi^2). $$
There are two functions in this ansatz. The function $\Phi(r)$ is the redshift function, and it controls the gravitational redshift and the possible formation of horizons. The function $b(r)$ is the shape function, and it controls the spatial geometry of the wormhole. The throat is located at a radius $r_0$ satisfying
$$
b(r_0)=r_0.
$$
At this radius the coefficient of $dr^2$ becomes singular. This should not be interpreted too quickly as a physical singularity; it is a coordinate effect associated with the choice of the radial coordinate $r$. The real geometric condition is not merely that $b(r_0)=r_0$, but that the surface actually flares outward at the throat instead of pinching off.
To see this, one considers an equatorial spatial slice by setting
$$
t=\mathrm{const},
\qquad
\theta=\frac{\pi}{2}.
$$
The induced two-dimensional metric is then
$$
ds^2
=
\frac{dr^2}{1-b(r)/r}
+
r^2d\phi^2.
$$
We now embed this two-dimensional surface into Euclidean three-space with cylindrical coordinates $(r,\phi,z)$. The Euclidean metric on a surface of revolution $z=z(r)$ is
$$
ds^2
=
\left(1+\left(\frac{dz}{dr}\right)^2\right)dr^2
+
r^2d\phi^2.
$$
Equating this expression with the wormhole slice gives
$$
1+\left(\frac{dz}{dr}\right)^2
=
\frac{1}{1-b(r)/r}.
$$
Hence
$$
\frac{dz}{dr}
=
\pm
\left(\frac{r}{b(r)}-1\right)^{-1/2}.
$$
At the throat, where $b(r_0)=r_0$, this derivative diverges. Geometrically, the embedded surface becomes vertical there. The more important quantity is the second derivative of $r$ with respect to $z$, since this measures whether the surface opens outward. A short calculation gives
$$
\frac{d^2r}{dz^2}
=
\frac{b(r)-rb'(r)}{2b(r)^2}.
$$
Thus the flare-out condition is
$$
b(r)-rb'(r)>0.
$$
At the throat this becomes
$$
b'(r_0)<1.
$$
This is the precise geometric statement that the wormhole opens out rather than closes in. It is a purely spatial condition, but in Einstein gravity it has an immediate dynamical consequence. The field equations convert this geometric inequality into a statement about the stress-energy tensor.
For the Morris-Thorne metric, the energy density and radial pressure in Einstein gravity satisfy, in an orthonormal frame,
$$
\rho(r)=\frac{b'(r)}{8\pi r^2},
$$
$$
p_r(r)
=
\frac{1}{8\pi}
\left[
-\frac{b(r)}{r^3}
+
2\left(1-\frac{b(r)}{r}\right)\frac{\Phi'(r)}{r}
\right].
$$
At the throat, the second term in $p_r$ vanishes because $1-b(r_0)/r_0=0$. Therefore
$$
p_r(r_0)
=
-\frac{1}{8\pi r_0^2}.
$$
Meanwhile
$$
\rho(r_0)
=
\frac{b'(r_0)}{8\pi r_0^2}.
$$
Adding these two expressions gives
$$
\rho(r_0)+p_r(r_0)
=
\frac{b'(r_0)-1}{8\pi r_0^2}.
$$
But the flare-out condition requires $b'(r_0)<1$. Hence
$$
\rho(r_0)+p_r(r_0)<0.
$$
This is exactly a violation of the null energy condition in the radial null direction. So the usual slogan that wormholes require “exotic matter” is not just a qualitative statement. It follows directly from the throat geometry together with the Einstein equations.
Modified gravity enters the story by changing this last step. In Einstein gravity, the relation between geometry and matter is
$$
G_{\mu\nu}=8\pi T_{\mu\nu}.
$$
Thus once the metric is fixed, the stress-energy tensor is fixed. In an $f(R)$ theory, however, the gravitational action is no longer the Einstein-Hilbert action
$$
S
=
\frac{1}{2\kappa^2}
\int d^4x\sqrt{-g}\,R,
$$
but instead
$$
S
=
\frac{1}{2\kappa^2}
\int d^4x\sqrt{-g}\,f(R).
$$
The model considered here is
$$
f(R)=R+\alpha R^n.
$$
This looks like a small deformation of general relativity, but it changes the field equations in an essential way. Varying the action with respect to the metric gives
$$
f_R R_{\mu\nu}
–
\frac{1}{2}f(R)g_{\mu\nu}
=
\nabla_\mu\nabla_\nu f_R
–
g_{\mu\nu}\Box f_R
+
\kappa^2T_{\mu\nu},
$$
where
$$
f_R=\frac{df}{dR}.
$$
For the particular model $f(R)=R+\alpha R^n$, one has
$$
f_R
=
1+\alpha nR^{n-1}.
$$
The new feature is the appearance of derivatives of $f_R$, and hence derivatives of the Ricci scalar $R$. In other words, the curvature itself now contributes additional effective stress-energy terms. It is often useful to rewrite the modified field equations schematically as
$$
G_{\mu\nu}
=
8\pi
\left(
T_{\mu\nu}^{\mathrm{matter}}
+
T_{\mu\nu}^{\mathrm{curv}}
\right).
$$
This formula should not be taken as saying that curvature has literally become matter. Rather, it is a bookkeeping device. The higher-curvature terms can be moved to the right-hand side of the equation and treated as an effective source. This changes the wormhole problem in an important way. The throat may still require null energy condition violation in the effective total source, but the violation need not come entirely from the ordinary matter sector.
The paper studies the shape function
$$
b(r)=re^{-(r-r_0)}.
$$
This is a convenient choice because the basic wormhole conditions can be checked explicitly. First,
$$
b(r_0)=r_0e^{-(r_0-r_0)}=r_0,
$$
so $r=r_0$ is indeed the throat. Next,
$$
\frac{b(r)}{r}=e^{-(r-r_0)}.
$$
Thus
$$
\frac{b(r)}{r}\to 0
\qquad
\text{as}
\qquad
r\to\infty.
$$
This gives the expected asymptotic flatness condition in the radial part of the metric. Differentiating the shape function, one obtains
$$
b'(r)
=
e^{-(r-r_0)}(1-r).
$$
Therefore
$$
b(r)-rb'(r)
=
re^{-(r-r_0)}
–
r e^{-(r-r_0)}(1-r)
=
r^2e^{-(r-r_0)}.
$$
Since $r>0$, this quantity is positive:
$$
b(r)-rb'(r)>0.
$$
So the flare-out condition is satisfied. In particular, at the throat one has
$$
b'(r_0)=1-r_0,
$$
and hence $b'(r_0)<1$ for $r_0>0$, as required.
The authors then consider both constant and variable redshift functions. The variable choice
$$
\Phi(r)
=
\ln\left(1+\frac{r_0}{r}\right)
$$
is especially useful because it is finite for all $r\geq r_0$. Since
$$
e^{2\Phi(r)}
=
\left(1+\frac{r_0}{r}\right)^2,
$$
the metric coefficient $g_{tt}$ never vanishes in the wormhole domain. Thus this choice avoids the formation of an event horizon while still producing a nontrivial tidal structure.
Once this metric is substituted into the $f(R)$ field equations, the resulting expressions for $\rho$, $p_r$, and $p_t$ become rather complicated. This is not surprising. The field equations contain not only $R$ but also derivatives of $f_R$, so the stress-energy components involve higher derivatives of the metric functions. The calculation is still systematic: choose $b(r)$ and $\Phi(r)$, compute the curvature scalar $R$, compute $f(R)$ and $f_R$, insert them into the modified field equations, and then read off the effective density and pressures.
The main point is not that the final formulas are elegant. They are not. The main point is that the modified curvature terms contribute to the effective energy budget of the wormhole. In Einstein gravity, the throat condition directly forces
<
p style=”text-align:center;”>
$$
\rho+p_r<0
$$
for the matter supporting the geometry. In $f(R)$ gravity, the analogous inequality applies to the total effective source:
$$
\rho_{\mathrm{eff}}+p_{r,\mathrm{eff}}<0.
$$
But this effective quantity contains both ordinary matter and curvature contributions:
$$
\rho_{\mathrm{eff}}+p_{r,\mathrm{eff}}
=
(\rho+p_r)_{\mathrm{matter}}
+
(\rho+p_r)_{\mathrm{curv}}.
$$
Thus it becomes possible, at least in principle, for the curvature sector to carry part of the exoticity. The ordinary matter sector may then satisfy the null energy condition in regions where the total effective source violates it.
This is the conceptual mechanism behind the model. The wormhole does not become free of exoticity in an absolute sense. The flare-out condition still demands an effective violation of the null energy condition. What changes is the location of that violation. In general relativity, it must be attributed directly to matter. In modified gravity, some of it can be absorbed into the higher-curvature terms.
This distinction is subtle but important. It means that the statement “wormholes require exotic matter” is not purely a statement about wormhole geometry. It is a statement about wormhole geometry together with a particular set of gravitational field equations. Change the field equations, and the same geometric requirement may be distributed differently between matter and curvature.
Of course, this should not be read as evidence that traversable wormholes are physically realized. An $f(R)$ model must still satisfy many independent tests: stability of the solutions, absence of ghost-like degrees of freedom, compatibility with solar-system and cosmological observations, and sensible behaviour under perturbations. These are serious constraints, and a formal wormhole solution is only the beginning of the story.
Nevertheless, the example is mathematically instructive. It shows that in gravitational physics the boundary between “geometry” and “matter” is more delicate than it first appears. In Einstein gravity the division is sharp: curvature sits on the left-hand side, matter on the right. In modified gravity, higher-curvature terms blur this separation. The same wormhole throat may therefore be interpreted not as being held open entirely by exotic matter, but partly by the effective stress-energy generated by the gravitational action itself.
That is the real lesson. Modified gravity does not merely enlarge the space of possible solutions. It also changes the bookkeeping of what it means for a geometry to be physically supported.