Wormhole formation in massive gravity: an analytic description

The present study analyses the wormhole solution both in the dRGT-f(R, T) massive gravity and Einstein massive gravity. In both the models, the anisotropic pressure solution in ultrastatic wormhole geometry gives rise to the shape function that involves massive gravity parameters γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \gamma $$\end{document} and Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Lambda $$\end{document}. However, the terms consisting of γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \gamma $$\end{document} and Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Lambda $$\end{document} acts in such a way that the spacetime loses asymptotic flatness. Similar to the black hole solution in massive gravity, this inconsistency arises due to the repulsive effect of gravity which can be represented by the photon deflection angle that goes negative after a certain radial distance. It is investigated that the repulsive effect induced in the massive gravitons push the spacetime geometry so strongly that the asymptotic flatness is effected. On the other hand, in this model, one can have a wormhole with ordinary matter at the throat that satisfies all the energy conditions while the negative energy density is sourced by massive gravitons. Finally, using the TOV equation, it is found that the model is stable under the hydrostatic equilibrium condition.


I. INTRODUCTION
Wormholes are the smooth bridges between two different universes, or sometimes between two distant parts of the same universe.The concept was first put forward by Einstein and Rosen in their famous Einstein-Rosen bridge in 1935 [1], where Misner and Wheeler first coined the term 'Wormhole' later in 1957 [2].The exact solution of the Einstein Field Equation for traversable wormhole were first successfully examined much later by Morris and Thorne in the year 1988 [3,4].They investigated that the energy-momentum component for such a wormhole always violates the null energy condition [3,5] which is the weakest classical energy condition that in turn violates all other energy conditions.So, for the construction of traversable wormhole one needs a matter with negative energy density called the exotic matter.For example, one may consider the well-known Ellis wormhole solution of General Relativity for instance [6][7][8][9][10][11][12].However, if traversable wormholes are constructed in modified theories of gravity, the requirement of exotic matter is highly reduced.There are numerous studies that deals with the wormhole geometries and their corresponding energy conditions in various modified gravities.For example, one may go through the wormhole solutions in f (R) gravity [13,14], in f (R, T ) gravity [15][16][17][18], and in other theories [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33].
Ever since the traversable wormhole model was presented, it was a point of interest if wormholes can be constructed with ordinary matter.Recently, a study has been introduced to discuss the wormhole solutions both in Einstein gravity and in modified gravity theories, and it is established that we may have a wormhole with ordi-nary matter in modified gravity, which may satisfy all the energy conditions [34].Although, the usual matter can be ordinary, but the effective geometric matter which is the matter source of modified gravity violates the usual null energy condition.There are various studies conducted on the ground that obtained non-exotic matter wormholes [35][36][37][38][39][40][41][42][43][44][45].
There are models of modified gravity theory as an extension of Einstein's GR to address cosmological phenomenon.The massive gravity theory is an important candidate in this context.The discovery of gravitational waves due to merger of two Black holes and massive stars as detected by LIGO and VIRGO [46,47], there is a severe constraint in the mass of the graviton, and thus, massive gravity theory has gained significant interest in theoretical regime.The de Rham-Gabadadze-Tolley (dRGT) massive gravity [48,49] is an interesting nonlinear generalization to overcome the so-called van Dam-Veltman-Zakharov discontinuity [50,51] of Fierz-Pauli linear massive gravity [52].The dRGT theory is utterly promising as it avoids the long standing issue of Boulware-Deser ghost instability [53] that usually arises on the introduction of nonlinear generalization of massive gravity theory.In the simplest dRGT model of static spherically symmetric solution, there arises two additional characteristics scales (i.e., γ and Λ) as compared to the Schwarzschild solution of General Relativity [54].In the theory, the massive graviton generates these parameters (one of which is the effective cosmological constant) [55,56] which can be used to address the dark matter and dark energy problems of the galactic and extragalactic scenario in a single theory.The nonlinear graviton interaction generating density and pressure which can behave as dark energy, and as a result, the galaxy rotation curve is influenced due to generation of dark matter halo.Thus, the massive gravitons may act as the dark matter halo resulting in asymptotically flat rotation curves.In [57], the dRGT model is quite successfully fitted with observational data for Milky Way rotation curves, and LSB galaxies without considering any additional dark matter.The consistency of dRGT model with the NFW profile has also been noted in this context.
Further, the massive gravitons play the role of an anisotropic fluid, which in turn is a kind of dark energy.As a consequence, one may able to explain the late time cosmic acceleration, although, it can not justify the past inflation.Subsequently, in a rectified framework, it is found that the massive gravity theory is consistent [58] with the Plank 2018 data [59] and its combination with BK18 and BAO [60].
The phenomenology in dRGT massive gravity for compact objects is interestingly promising.Various studies have already been conducted on Black hole, Black string, rotating Black string solutions, stability and greybody factor on charged Black hole and Black strings in dRGT model which are modified significantly from the conventional studies due to the presence of massive gravitons [61][62][63][64][65][66][67][68].
On the other hand, the existence of wormhole demands the violation of null energy condition.At the same time, a natural outcome in the massive graviton density is that it can be negative in a certain region where the radial pressure is positive.So, in the massive gravity theory, the massive graviton energy-momentum tensor naturally fulfil the feature of energy condition violation [69].Thus, the existence of wormhole in the massive gravity can be a natural outcome, and the discussion of this particular object has an additional importance on its own.There are only couple of studies present on wormhole solutions in dRGT massive gravity with its own limitations.Tangphati et.al. investigated the traversable wormhole solution in f (R) massive gravity [70], where they considered an exponential shape function in order to discuss the asymptotic geometry.But, due to their restricted choices of shape function, they were unable to obtain a geometry which includes the effect of massive gravitons.Alongside, Kamma et.al. presented the wormhole solution in conventional Einstein gravity on the background of massive gravitons [71], and obtained both the redshift function and shape function which includes massive gravity parameters.Although, they discussed about negative pressure tension in wormholes, but they didn't shed any light on the asymptotic geometry triggered by the redshift and shape function.
So, in this article, we shall investigate the wormhole solution in f (R, T ) massive gravity (the field equations are discussed in section II), and discuss the corresponding anisotropic pressure solution in the ultrastatic wormhole geometry in section III.It is observed that a repulsive gravity effect is present and the asymptotic structure is violated in the model which is not so well-known in wormhole geometries.Hence, to represent the relationship between asymptotic structure and the repulsive gravity in a more generalized scenario, we examined wormhole solution in dRGT extension of Einstein gravity in section IV.On the other hand, the stability of the model is analyzed in section VI followed by the discussions of energy con-ditions and matter content in section V. Finally, section VII is dedicated for the conclusions and discussions.
We consider the natural units throughout the study, i.e., G = c = 1.

II. THE FIELD EQUATIONS
To investigate the f (R, T ) gravity model in the context of de Rham-Gabadadze-Tolley (dRGT) massive gravity, we consider the action in the following form where f (R, T ) is function of R and T , U is the selfinteracting potential of the graviton with graviton mass m g , L m is the matter Lagrangian, and g is the determinant of the metric tensor g µν .Here, the potential U is defined as ( where U 2 , U 3 and U 4 are given by Here, [K] represents the trace of K µ ν , where (K i ) µ ν = K µ ρ1 K ρ1 ρ2 ...K ρi ν , ϕ a is the Stückelberg field, and the reference fiducial metric F ab has the explicit form given by where the unitary gauge is fixed as, ϕ a = x µ δ a µ , so that Now, by varying the action with respect to the metric g µν , we arrive to the field equation for modified dRGTf (R, T ) massive gravity theory given by where f R (R, T ) and f T (R, T ) are the differentiation of f (R, T ) with respect to R and T respectively, and □f R = g µν ∇ µ ∇ ν f R .We also have the variation of trace of energy-momentum tensor of the matter field, T = g µν T µν as where Θ µν and T µν are given by Here, T µν can also be written in the principal pressure terms as where u µ is the timelike unit vector, χ µ is a spacelike unit vector orthogonal to timelike unit vector, such that u µ u µ = −1 and χ µ χ µ = 1.Assuming the universal choice of Lagrangian matter density L m = ρ, we have Θ µν = −2T µν + ρg µν .Alongside, χ µν is the massive graviton tensor, given by where the parameters α and β are defined by According to the definition, which is nothing but the energy momentum tensor of the massive gravity sector.Using Eq. ( 11), one can readily calculate the density and pressure components ρ (g) (r) and p (g) r,⊥ (r) which are given by [57,70,72,73] where effective cosmological constant Λ and a new parameter γ are introduced and written in the linear combination of the parameters α and β, given by 17) Here at this point, after applying some simple mathematical calculations on Eq. ( 6), we arrive to the final field equation. where Note that, when f (R, T ) ≡ f (R), such that f T (R, T ) = 0, we get back the usual f (R) massive gravity solution [70].In this field equation, the energy-momentum tensor component (T µν ) of f (R, T ) gravity represents the interaction between matter and curvature, and one may interpret this as the curvature-matter coupling occurs due to the exchange of energy and momentum between the both.Otherwise, the total energy-momentum tensor which is the sum of f (R, T ) and massive gravity sectors, adds the interaction of massive gravitons with the curvature-matter coupling.It can be expressed as (21) However, f (R, T ) gravity proposed by Harko et al. [74] has a critical issue that violates the energy conservation law, and it leads to non-geodesic motion of test particles.But, an alternative approach of f (R, T ) gravity proposed by Chakraborty [75] showed that the form of the field equation remains similar if one takes into account of the conservation of EM tensor.As a result, test particles move in geodesic orbits, and the choice of Lagrangian is not completely arbitrary.For homogeneous and isotropic model of the universe, the approach refers to the field equation which is equivalent to the Einstein gravity with non-interacting 2-fluid system, one of which is the usual perfect fluid in modified gravity and the other one shows exotic nature.
Recently, in [76], it is showed that in f (R, T ) gravity, one may have a wormhole with non-exotic matter where the ordinary EM tensor satisfies the NEC, and the extra curvature caused by the modified gravity manipulate violation of NEC, so that the total EM tensor violates NEC.In a general perspective of modified gravity, the existence of ordinary matter wormhole is briefly discussed in [34].For another example of non-exotic matter wormhole in f (R, T ) gravity, readers are referred to [77].Now, it is mathematically evident that the massive gravity EM tensor holds the property of anisotropic dark energy, as p (g) r = −ρ (g) , which may successfully manipulate the phenomenological aspects of various compact objects and astrophysical phenomenons [57,[61][62][63][64][65][66][67][68].Thus, in f (R, T )-massive gravity, the massive gravity and curvature-matter coupled system has the massive gravitons and effective geometric matter (due to the curvature) that can exhibit exotic nature.Hence, if the EM sector of usual matter satisfies the energy conditions and is dominated by the sum of gravitons and geometric matter, we can have a usual matter wormhole.Now, to construct traversable wormhole solution, we consider the Morris-Thorne line element given by [3] where dΩ 2 = dθ 2 + sin 2 θdϕ 2 , Φ(r) and b(r) are the redshift function and shape function respectively.The minimal requirement of traversability of this particular wormhole geometry demands the following conditions: 1.The wormhole is constructed by connecting two asymptotic flat regions at the throat.The throat radius is defined by a global minimum r = r 0 , so the radial coordinate runs in the interval r ∈ [ r 0 , ∞) .
2. The redshift function Φ(r) must be finite everywhere in order to avoid the presence of horizons and singularities.So, e Φ(r) > 0 everywhere for r > r 0 .
In this context, the ultrastatic wormhole is a particular point of interest which defines the zero-tidalforce wormhole.Here, Φ(r) = 0, so that e Φ(r) = 1 i.e., in a frame free from gravitational acceleration, a particle dropped from rest remains at rest [3,78].

The flaring-out condition
must hold at or near the throat r = r 0 .
4. The above mentioned conditions imply that b(r 0 ) = r 0 and b ′ (r 0 ) ≤ 1 for all r ≥ r 0 , where the equality of b ′ (r 0 ) only holds at the throat.Further, for r > r 0 ⇒ b(r) < r.
The redshift and shape function must obey these conditions (readers are referred to [3,70,[79][80][81], for details).Now, the Einstein tensor components for the spacetime metric are given by [3] Therefore, using these components, the calculation of field equation (from Eq. ( 18)) of f (R, T ) massive gravity is rather straightforward to obtain, Now, to construct wormhole solution, one can consider restricted choices of Φ(r), b(r) and f (R, T ) among others.The other way considers specific pressures (isotropic/anisotropic) or equations of state for p r , p t .Although, fixing the shape function prior to the solution may loose the effect of massive gravity on the shape.For an example, we can consider the wormhole solution in f (R) massive gravity by Tangphati et.al. [70].In this present study, we will consider the anisotropic pressure solution in f (R, T ) massive gravity and exhibit the effect of massive gravitons on the wormhole shape function.

III. ANISOTROPIC WORMHOLE SOLUTION
For the anisotropic pressure fluid wormhole solution in dRGT-f (R, T ) massive gravity, we consider a simple choice of f (R, T ), e.g., f (R, T ) ≡ (αR + βT ) [98], where α = 1 and β = 2λ brings forth the well known linear choice of f (R, T ) cosmology, i.e., f (R, T ) ≡ (R + 2λT ) [74].We have f R = α, f T = β, and f ′ R = f ′′ R = 0 for the specific choice.On the other hand, ultrastatic wormhole is chosen for the particular solution which results Φ ′ (r) = Φ ′′ (r) = 0, and we can directly compute the reduced EM components to get, The anisotropic pressure fluid equation is given by p t = σp r , where σ = 1 defines the isotropic pressure condition (i.e., p t = p r = p).At the same time, σ must be less than zero for anisotropic pressure fluid to ensure the asymptotic flatness of the spacetime.Hence, for our particular choices, we can directly obtain the shape function by using equations ( 30) and (31).
where C is the constant of integration.It is expected to note that the shape function effectively depends on the massive gravity parameters (γ and Λ), and hence on the mass of the graviton.
It is worth noting here that for σ = −1 and γ = Λ = 0, the equation reduces to the so called Ellis wormhole [6] shape function.
However, one can obtain the position of wormhole throat from the root of 1 − b(r) r = 0, which (for our area of interest) eventually gives us where For Ellis wormhole, Eq. ( 33) becomes It is clear that the equation has four roots, and one of the real root determines the throat radius.One can easily impose condition 4 (from the 5 minimum requirements for wormhole geometry, section II) i.e., b(r = r 0 ) = r 0 to determine C, so that the throat radius r 0 can be specified.Hence, the final shape function is given by b The result is much expected where massive gravity free (i.e., γ = Λ = 0) expression turns out the usual solution of anisotropic pressure fluid in f (R, T ) gravity.However, due to the presence of r and r 2 respectively with γ and Λ in the second term of the shape function changes the asymptotic structure of the spacetime.The effect is very much similar to the black hole solution in dRGT massive gravity.We refer [82] for the discussion in black hole asymptotic structure.
To explore the shape function characteristics discussed in section II, we plot b(r), b ′ (r), b(r)/r, and the flaringout condition (i.e., −rb ′ (r) + b(r)) in Fig. 1.It is observed that, though all the properties are satisfied near the throat, the conditions given by b ′ (r) < 1, b(r)/r < 1 for all r > r 0 , and b(r)/r → 0 for r → ∞ are significantly violated after certain radial distance.So, it is straightforward to conclude that the asymptotic flatness of the spacetime vanishes after certain distance due to the presence of massive gravitons.We can interpret this as the result of repulsive effect of gravity, which is another unique feature of the dRGT massive gravity theory.A similar effect arises in the black hole solution in massive gravity where the asymptotic structure changes due to the said effect [82].It should be noted here that repulsive behaviour of gravity is also found in the BHT massive gravity [83], and in other modified gravities with exotic matters and energy [84][85][86][87][88].
Coming back to the discussion of non-asymptotic geometry, for r → ∞, b(r)/r did not approach ′ 0 ′ , rather it goes to ′ ∞ ′ due to the presence of linear r and r 2 terms with γ and Λ respectively in Eq. (36), resulting the loss of asymptotic flatness.It is evident and can be exhibited mathematically that other choices of f (R, T ) in ultrastatic wormhole geometry also involve such r and r 2 terms with γ and Λ in the shape function.Similar circumstances can be observed in other gravities such as Einstein's gravity in the background of massive gravitons (which will be discussed later in this literature), and in f (R) massive gravity wormholes [70].Hence, we can conclude that it is a generic feature of wormholes in dRGT massive gravity.Now, we need to check the presence of repulsive gravity effect in the wormhole solution to verify if it has any role on the asymptotic structure.To easily verify the phenomenon, one can introduce the photon deflection angle on the wormhole, which eventually goes negative in a spacetime if repulsive gravity acts on the photons [82].

A. Repulsive Behaviour of Gravity
To check the deflection angle of photons from null geodesics, we first introduce a general spherically symmetric and static line element given by [89,90] The geodesic equation, which relates the momenta oneforms of a freely falling body and background geometry is given by [90] where λ is the affine parameter.One can immediately tell that if the components of g αν are independent of x β for a fixed index β, then p β is a constant of motion.Hence, if we consider only the equatorial slice by setting θ = π/2, then all the g αβ are independent of t, θ, ϕ in Eq. ( 38), i.e., one can obtain the respective killing vector fields δ µ α ∂ ν with α as a cyclic coordinate.Now, we can set the constants of motion p t and p ϕ as where E and L are the energy and angular momentum of the photon respectively.Thus, we have where the overdot represents the differentiation w.r.t.affine parameter λ.Again, the radial null geodesic can be obtained easily as However, one can write the equation for the photon trajectory in terms of impact parameter µ = L/E, as Now, one can obtain the deflection angle of photon by considering a source of photon radius r s causing the geometry, then the photons can hit the surface only when an existing solution r 0 obeys the condition r 0 > r s , and ṙ2 = 0. Here, r 0 is the distance of closest approach or turning point.In that case, the impact parameter becomes and it is obvious that for weak gravity limit, µ ≈ C(r 0 ).Thus, if a photon is coming from polar coordinate limit given by lim r→∞ r, − π 2 − α 2 , passes through the turning point at (r 0 , 0) and approaches lim r→∞ r, π 2 + α 2 , then the deflection angle of photon is defined as this α, which is a function of r 0 [91].We can compute it from Eq. (42) as For the choices of metric coefficients in wormhole geometry, the deflection angle becomes One can now readily exhibit the deflection angle of photons in massive gravity by numerically integrating the above equation after imposing the shape function given by Eq. (36), which is shown in Fig 2.
It is interesting to observe that the deflection angle becomes negative after a certain value of r 0 .We can interpret this as the effect of repulsive gravity, where neglecting massive gravity parameters γ and Λ (i.e.neglecting massive gravity), one can comfortably verify the The plot exhibit the negative deflection angle representing the repulsive nature of gravity.It also verifies the reason of non-asymptotic behaviour.
non-existence of negative deflection angle [92].At the same time, referring Fig. 1, it is verified that the inconsistencies in the spacetime structure like non-asymptotic flatness occurs beyond the radial distance where the repulsive effect of gravity dominates.So, we can conclude that the loosening of asymptotic structure results from the repulsive behaviour of gravity.However, the repulsive anisotropic dark energy nature is inherent in massive gravitons.So, their presence may physically cause the repulsive gravity effect in dRGT massive gravity.

IV. WORMHOLE IN EINSTEIN-MASSIVE GRAVITY
To start the discussion of Einstein gravity with the background of dRGT massive gravity, the action is considered as We can now vary the action with respect to g µν after considering the massive gravity and other parameters from section II and get However, one can easily derive the above results just by considering f (R, T ) = R, i.e., f R (R, T ) = 1, in f (R, T ) massive gravity.
In Eq. ( 47), the EM tensor of the Einstein gravity follows the energy conservation law, and the total EM tensor defined by the coupled perfect fluid-massive graviton system is given by where 'E ′ in the superscript represents the EM components for Einstein gravity.So, one can directly write ∇ µ T tot µν = 0, ∇ µ T µν = 0 and ∇ µ T (g) µν = 0.In Einsteinmassive gravity, the construction of wormhole can be possible with non-exotic matter source, where the violation of energy conditions may manipulated by the exotic nature of massive gravitons.Now, to derive the calculations of field equation for Morris-Thorne wormhole, we may use Eq. ( 23), (24), and (25) in Eq. ( 47) and obtain, Here, we are about to construct an anisotropic wormhole solution with restricted choice of redshift function, e.g.Φ = constant (i.e.Φ ′ = Φ ′′ = 0), which represents the ultrastatic wormhole geometry.Hence, the equation of state for the anisotropic pressure solution, p t = σp r provides the shape function of the wormhole.By using the boundary condition at the throat junction, i.e. b(r 0 ) = r 0 , the final shape function is given by b The shape function is very much similar to that obtained in f (R, T ) gravity, and the visualization of the Morris-Thorne type shape function properties, i.e. b(r), b ′ (r), b(r)/r, and (−rb ′ (r) + b(r)) are plotted in Fig. 3.The plots violate the asymptotic flatness after a certain radial distance, as it also contains linear r and r 2 terms with γ and Λ respectively.However, if the expression is freed from the massive gravity parameters, i.e. γ = Λ = 0, asymptotically well-behaved shape function is recovered.The discussion of repulsive effect of gravity in this ground is indeed necessary, and we directly use Eq.(45) to plot the numerical solution of deflection angle of photons in dRGT extension of Einstein gravity.The existence of negative deflection angle indicate the presence of repulsive gravity, where all the other comprehensive discussions are similar to the f (R, T )-massive gravity.The plot is exhibited in Fig. 4, and comparing it with Fig. 3, it is satisfied that the Morris-Thorne wormhole shape function properties are violated after the radial distance where the deflection angle goes negative.So, also in Einstein-massive gravity, the asymptotic flatness of the traversable wormhole is lost due to the repulsive effect of gravity, and it is evident that it arises from the anisotropic dark energy nature (repulsive nature) of massive gravitons.

V. ENERGY CONDITIONS
Discussing energy conditions in traversable wormhole configuration is highly important to have an idea of the matter content on the wormhole.As discussed, the throat of usual wormholes must be threaded by a matter of negative energy density that violates the Null energy condition (NEC) and is termed as 'exotic matter'.The matter is used to keep the wormhole throat open, thus making it traversable.
Here, in this section, we will discuss the energy conditions for wormhole configuration consecutively in f (R, T )-massive gravity and Einstein-massive gravity.
From the definitions and mathematical fundamentals, one can summarize the energy conditions in principal pressure forms which is given by (i) NEC : Now, one can use standard mathematical tools to calculate ρ, ρ + p r , ρ + p t , ρ − |p r |, ρ − |p t |, and ρ + p r + 2p t from Eq. ( 29), (30), and (31) by imposing the shape function of Eq. (36).Using these six terms, all the four energy conditions for f (R, T )-massive gravity can be investigated.The mathematical results of the terms are plotted in Fig. 5 for α = 1, β = 2, σ = −1, and r 0 = 1.Alongside, the summaries of the results from the plots are listed in Table I.
From Fig. 5 and Table I, we can conclude that all the energy conditions including the energy density are satisfied throughout the wormhole from r = 0.87 and r = 0.83, respectively for γ = 0.01, Λ = 0.7 and γ = −0.01,Λ = 0.7.On the other hand, one can also conclude that the energy condition components does not differ much for the signs of γ.However, they are highly dependent on the effective cosmological constant (Λ).For γ = 0.01, Λ = −0.7,all the components are completely negative, resulting violation of energy conditions.Alongside, for smaller γ and Λ terms, i.e., for less effects of massive gravity, the energy conditions are violated at the wormhole throat, but gets satisfied after a certain radial distance.
Therefore, it is interesting to note that, for the first two choices of γ and Λ, i.e., γ = 0.01, Λ = 0.7 and γ = −0.01,Λ = 0.7, the wormhole can be constructed with ordinary matter, which satisfies all the energy conditions.The situation is investigated in [34] with great details, where it is established that we can have wormholes with ordinary matter in modified theories of gravity.However, the violation of NEC is customary at the wormhole throat to keep it open for traversability.They investigated that the presence of effective geometric matter (e.g. the curvature EM component in f (R, T ) gravity) can play the role of energy condition violation.Few of the recent studies in f (R, T ) gravity [76,77] deal with the wormhole models with non-exotic matters, where the curvature matter source acts as the exotic matter and performs in the violation of energy conditions.Here, in this current scenario, the matter source of massive gravitons (which acts as the anisotropic dark energy) added with the curvature matter source, violate the energy condition components.So, the usual matter source can be ordinary matter.However, the repulsive dark energy nature in massive gravitons is so strong that it brings forth the repulsive gravity effect in the spacetime and pushes the geometry so strongly that the asymptotic flatness is effected.
On the other hand, the Einstein-massive gravity does not differ much from our linear f (R, T )-massive gravity model.By imposing α = 1, and β = 0 in Eq. ( 29), (30), (31), and (36), we can readily recover EM components  and the shape function for Einstein gravity.Mathematically, the EM components are very loosely dependent on the β values, and graphically, the energy conditions possess adequately similar behaviour for α = 1, and β = 0.For γ = 0.01, Λ = 0.7 and γ = −0.01,Λ = 0.7, it also contains ordinary matter wormhole.But here, only the massive gravity sector adds the energy condition violations.So, it is evident that the discussions for Einsteinmassive gravity model are similar to that of the f (R, T ) gravity.

VI. EQUILIBRIUM ANALYSIS
The equilibrium condition for our present work is provided by the generalized Tolman-Oppenheimer-Volkov (TOV) equation which is given by The above equation is an important and elegant method and can be used to examine the stability condition of the astrophysical solutions including wormhole/compact objects.Readers may refer to [77,[93][94][95][96][97] for detailed study.
The Eq. ( 53) can also be written as which provides the equilibrium condition for the worm-hole.Here, where, F a denotes the force due to the anisotropic matter of the wormhole, F g is the gravitational force, and F h is the hydrostatic force.F a arises due to the modification of the gravitational Lagrangian of the Einstein-Hilbert action.However, it is clear from Eq. ( 54) that for the system to be in equilibrium, the sum of the three different forces must be equal to zero.In Fig. 6, the behaviour of F a , F g , and F h are shown for a particular choice of the parameter values where the energy conditions are satisfied.However, the value of F g is zero due to the assumption of Φ to be constant (i.e. the gravitational force has no effect on our model).From the figure, one can visualize that the other two forces are exactly same and opposite to each other.Thus, it is evident that the equilibrium of forces is achieved due to the combined effect of the three force terms, and hence this supports the stability of the system.

VII. DISCUSSIONS
In this work, we have presented the wormhole solution in the dRGT massive gravity extension of f (R, T ) Here F g = 0 for the ultrastatic wormhole choice.
gravity and Einstein gravity.For the anisotropic pressure wormhole solution, we considered a linear choice of f (R, T ), i.e., f (R, T ) ≡ αR + βT , which is a straightforward choice of the linear f (R, T ) cosmology model, i.e., f (R, T ) ≡ R + 2λT .One of the advantages in this particular choice is that, one can easily reduce this to Einstein gravity by considering α = 1 and β = 0. On the other hand, from this model, f (R)-massive gravity can also be recovered by suitable choices of parameters.Although, in f (R, T ) gravity, the energy conservation law is violated [74], but still wormhole solutions (also with non-exotic matter) can be constructed in this theory.For example, readers can go through the following literature [76,77].However, one may also consider conservation of energy density, where interestingly the field equations remains similar [75], and the possibility of non-exotic matter wormhole can be investigated.The shape function obtained from the anisotropic pressure solution in the model involves product terms of r and r 2 respectively with the massive gravity parameters γ and Λ, which causes non-asymptotic flat geometry in the spacetime after a certain radial distance.Although, we can extend the region of asymptotically well-behaved spacetime by considering smaller values of γ and Λ.However, a globally complete asymptotic flat wormhole solution can be achieved by imposing restricted choices of shape function.The non-asymptotic flatness and the reason is similar to the black hole solution in massive gravity where γr and Λr 2 terms arise in the metric coefficients [82].
The repulsive effect of gravity in a spacetime can be examined by the negative deflection angle of photons.It is investigated that, in an extended spacetime description with massive gravity, the deflection angle goes negative due to the terms consisting of γ and Λ.Hence, the nonasymptotic flatness and the presence of repulsive gravity is the generic feature of dRGT massive gravity, and the massive gravitons source the effects both in the black holes and wormholes.
On the other hand, if we investigate the equation of state (p r = ωρ) to obtain the shape function, we get One can readily verify that the expression is very similar to the anisotropic pressure solution, so it also replicates all the properties of the anisotropic solution.
The dRGT-Einstein gravity is the generalization of f (R, T ) model with α = 1, and β = 0, so it reproduces all the effects and analysis present in the f (R, T ) model.The shape function violates the Morris-Thorne properties, and the asymptotic flatness is lost just beyond the radial distance where the repulsive gravity effect arises due to the massive gravitons.
The discussion of wormhole solution in massive gravity has a particular point of interest as the energymomentum tensor arising from the massive gravitons naturally fulfil the violation of null energy conditions.Further, the introduction of f (R, T ) modification in massive gravity is more viable to construct wormholes with wellbehaved matter source.In wormholes of modified gravity, there are two kinds of matter sources, one of which is the usual matter and another is the effective geometric matter which arises from the geometry of modified gravity [34].In the current study, it is established that for different values of γ and Λ, there are large number of possibilities of ordinary matter wormhole that satisfies all the energy conditions.Here, the usual matter source of wormhole is ordinary matter, where the matter source of massive gravitons coupled with the curvature term of f (R, T ) gravity act as the geometric matter that sources the negative energy density necessary to keep the wormhole throat open.It is due to the exotic anisotropic dark energy nature of coupled curvature and massive gravitons.But for Einstein-massive gravity, it is only the massive gravitons that sources the exotic component in non-exotic matter wormhole.In this context, it is customary to note that the strong repulsive nature of massive gravitons give rise to the repulsive effect of gravity and pushes the geometry for non-asymptotic flatness, both in f (R, T ) and Einstein massive gravity.
Finally, the Tolman-Oppenheimer-Volkov (TOV) equation is investigated for the equilibrium condition analysis.Although, the gravitational force is zero due to the constant redshift function, but the wormhole model is stable by the interaction of anisotropic and hydrostatic force terms.

TABLE I :
Numerical results for the possible regions where respective energy conditions are satisfied.