Unraveling the mysteries of wormhole formation in Rastall-Rainbow gravity: A comprehensive study using the embedding approach

The present work looks for the possible existence of static and spherically symmetric wormhole geometries in Rastall-Rainbow gravity. Since, the Rastall-Rainbow gravity model has been constructed with the combination of Rastall theory and the gravity's rainbow formalism. Taking advantage of the Karmarkar condition for embedding class one metrics, we solve the modified field equations analytically that describe wormholes for specific choice of redshift function. For specific parameter ranges, the solution represents a traversable wormhole that exhibits the violation of null energy condition and consequently the weak energy condition also. Furthermore, we focus on the wormhole stability via adiabatic sound velocity analysis. This model establishes a strong connection between two model parameters, namely, the Rastall parameters and the Rainbow functions, and how it affects the wormhole solution.


I. INTRODUCTION
Traversable Lorentzian wormholes [1] are of fascinating objects predicted by Einstein's theory of gravity.To be specific, wormholes are handle-like structures that connect two spatially separated locations in our universe or two different universes.Historically, the idea of wormholes was first theorized by Flamm in 1916 [2], shortly afterwards Einstein proposed his general relativity (GR) theory.Almost twenty years later, Albert Einstein and Nathan Rosen [3] proposed the existence of "bridges" between two different points in spacetime.This idea open up the possibility of a shorter route that could reduce travel time and distance also.This simplest form of mathematical model was dubbed "Einstein-Rosen bridges" for what the physicist John Wheeler would later coin a word "wormhole."However, their solution is prohibited by physical rules due to the presence of a physical singularity.But, it was Morris and Thorne [1], who came up with an astonishing idea that could be a two-way passage through them without processing any event horizon.According to them, a traversable wormhole has a throat that connects two asymptotically flat regions.In this setting wormholes throat are identified as a minimal surface area satisfying the flare-out condition [1].On the other hand Morris, Thorne and Yurtsever [4] studied the stability of transforming wormholes into a time machine.Visser systematically systematized all these seminal works in his book "Lorentzian wormholes: from Einstein to Hawking" [5].
The fundamental problem associated with traversable wormholes is the presence of matter that violates the energy conditions to satisfy the flare-out condition [4].Such a matter is referred to as "exotic matter" (the stress-energy tensor of matter that violates the null energy condition).The presence of exotic matter has a repulsive nature that prevents the throat from shrinking.However, Visser et al [6] had pointed out that the matter violating the averaged null energy condition to support traversable wormholes can be made infinitesimally small by means of an appropriate choice of wormhole geometry.Further, the amount of exotic matter can be confined at the wormhole throat by the cut-and-paste technique, yielding a thin-shell wormhole solution can be found in Refs.[7,8].Thin-shell wormholes have been studied by many authors [9][10][11][12].
These objects have been extensively studied within the context of GR.The Morris-Thorne wormholes in the presence of a generic cosmological constant has been studied in [13].Wormhole solutions have also been obtained in support of phantom energy EoS [14][15][16].The key argument is that the energy density of phantom field increases with time and thus provides a notion for the existence of wormholes.Lobo [17] discussed the stability equilibrium configurations for the specific phantom wormhole geometries by applying the linearized spherically symmetric perturbations.In the same line of action, various solutions regarding wormholes were found for generalized Chaplygin gas [18][19][20], varying cosmological constant [21], polytropic phantom energy [22] and ghost scalar fields [23].Subsequently, evolving wormholes have also been studied in support of phantom matter and cosmological constant [24,25], see also Refs.[26,27].
Therefore, it goes without saying that the theoretical construction of wormhole geometries with ordinary matter is a challenging task in gravitation physics.But, we are still far away from the success-line at least in GR context.To overcome this situation stated above, one usually considers modified theories of gravity.Since, the modified theories introduce a new degree of freedom that helps a sustainable wormhole solution while keeping the matter components non-exotic.In alternative theories of gravity, for example, higher order gravity theories [28,29], higher-dimensional cosmological wormholes [30] and in the Einstein-Gauss-Bonnet theory [31][32][33] authors have studied wormhole solutions.In support of this, one can find traversable wormholes in f (R) gravity with ordinary matter [34,35], or sourced by dark matter [36,37].Indeed, wormhole geometries have also been explored in distinct modifications of GR, namely: Extended theories of gravity [38][39][40][41][42][43][44], third-order Lovelock gravity [45,46], hybrid metric-Palatini gravity [47,48] and f (Q) gravity [49][50][51].A similar analysis has been performed in f (R, T ) gravity, see Refs.[52][53][54].Interestingly, authors in [55][56][57] have shown that one can theoretically construct exact wormholes solutions without requiring exotic matter.
In [58], authors have proposed a new modified gravity theory, namely, the Rastall-Rainbow gravity.Since, the Rastall-Rainbow gravity model has been constructed with the combination of Rastall theory [59] and the gravity's rainbow formalism [60].Thus, the closed-form wormhole solution in Rastall-Rainbow gravity has received significant attention from the researchers.The focal point of our discussion is the space-time of embedding class-one.A spacetime is said to be of embedding class N if it can be locally embed in an (4 + N )-dimensional pseudo-Riemannian spacetime.Interestingly, the idea of embedding establishes a geometric relation between the metric functions and their derivatives for general spherically symmetrical line-element [61].This means that one can choose one of the metric functions and generate the other metric function.This idea, the so-called "Karmarkar condition" has been used to address new solutions of relativistic astrophysical solutions, such as compact stars [62][63][64][65] and in wormhole solutions [66,67] also.
The above discussion motivated us to search for wormhole solutions within the Rastall-Rainbow gravity context.The main attraction is the use of Karmarkar conditions on wormhole geometries and investigate the effects of the Rastall-Rainbow parameters on it.Our work plan is: In Section II, we briefly review the Rastall-Rainbow gravity theory and derive the corresponding modified field equations for the configurations under consideration.In the same section we establish a relation between two metric functions via embedding approach.In Section III, we choose a fairly simple redshift function that ensures the space-time is free from horizons and singularities.Using the chosen redshift function, we derive the shape function and studied all the constraints for wormhole solution.The structure equation for standard energy conditions is presented in Section IV.In Section V, we present some features concerning wormhole solution under consideration, specifically, adiabatic sound velocity and energy conditions.Finally, Section VI is devoted to summarize and discuss our findings.

A. Rastall-Rainbow theory
In Einstein's GR, the conservation law of the energy-momentum tensor is stated as T ν µ;ν = 0. Nevertheless, Rastall gravity is a modification of the GR suggested by Rastall [59].In Rastall gravity, the conservation law of the energymomentum tensor is not conserved, and expressed in the following form where λ = 1−λ 16πG .Here, λ stands for the Rastall parameter, which is a constant denoting the deviation from GR and describing the matter field's affinity for coupling with geometry.When λ = 1, the standard law of conservation is recovered.What's more, in the case of a flat space-time where the Ricci scalar R = 0, the standard conservation law also holds.On the other hand, in the presence of Rastall gravity, where λ ̸ = 1 and space-time is non-flat, the above formula captures the modified conservation law effect.The aforementioned formula can be expressed as follows, So, in the case of Rastall gravity, Einstein's equation can be modified to take a slightly different form, expressed as follows, which can be further simplified to the following form, The corresponding expression for the trace of the stress-energy tensor is now obtained in the following way, Accordingly, the modified form of Einstein's equation in Rastall's gravity can be expressed as, In 2002, Magueijo and Smolin [68] introduced the gravitational rainbow as an extension of doubly special relativity to curved spacetimes.The gravitational rainbow involves a distortion of spacetime caused by two arbitrary functions Π(x) and Σ(x), known as rainbow functions.These functions are governed by the following conditions: Here, we are defining x as the ratio between E and E P l i.e., x = E/E P l .In this context, E stands for energy, υ represents momentum, m denotes the mass of a particle under test, and E P l = √ ℏc 5 /G corresponds to the Planck energy.Awad et al [69] and Khodadi et al [70] selected Π(x) = 1 and Σ(x) = √ 1 + x 2 for studying solutions associated with a nonsingular universe.In addition, the exponential form of the rainbow has been employed in [69,71] for investigating gamma-ray bursts.In the absence of test particles, the rainbow functions fulfill the following conditions: In this context, the metric characterizing space-time depends on the energy [60], and can be expressed as, In this case, the energy-dependent Vierbein fields e µ a (x) are connected to the independent Vierbein fields e µ a by the following relationships: Here, the subscript k = (1, 2, 3) stands for the spatial coordinates.In the light of rainbow gravity, Einstein's field equations are modified, based on the assumption that the geometry of space-time is related to the energy of the particle under test.Therefore, all amounts included in the field equations of this theory of gravity are energy-dependent.Contrary to the traditional Einstein field equations, a set of alternative field equations is being incorporated, of which the following is a typical one: where k(x) = 8πG(x).These modifications are of great interest, as they allow us to better elucidate the relationship between gravity, high-energy physics, and the fundamental structure of spacetime.The Rastall-Rainbow gravity model [58] evolves by combining modifications of Rastall gravity and Rainbow gravity.It introduces a unified framework that captures both the non-minimal coupling between gravity and matter, along with the energy-dependent behavior of particles.This model furnishes an expanded description of gravity that goes beyond the classical formulation of GR.In this unified formalism, the field equations are given by, Here, the Rastall parameter λ is supposed to be independent of the particle energy under test.In particular, within the framework of Rastall-Rainbow gravity, one can construct relativistic stars [72][73][74] or exotic objects such as wormholes [75].In the forthcoming sections, we shall concentrate on illustrating the simplest form of wormhole geometry and analyzing its geometric properties.Here, we assume G(x) = 1 for our entire calculation.

B. The wormhole geometry and its corresponding field equations
To construct a shape function satisfying embedding class one and characterizing the geometry of a spherically symmetric and static spacetime, we start with the metric ansatz given by The geometry of spacetime is governed by the functions Φ(r) and Λ(r), which depend on the radial coordinate r.These functions control the gravitational field and are ultimately involved in identifying the non-zero Riemann curvature components ( 13) corresponding to the spacetime configuration under consideration.
The components of Riemann curvature listed above satisfy the well-known Karmarkar condition, In the space-time geometry environment, the embedding class-one concept relates to a very specific form in which the non-zero components of the Riemann curvature tensor can be embedded in a higher-dimensional flat space.This concept enables us to interpret curvature geometrically by means of an embedding function.Following that, we apply the Karmarkar condition by replacing the non-zero components of the Riemann curvature tensor.More precisely, in the case of a static spherically symmetric spacetime, the non-zero components are generally given by, where primes denote derivatives with respect to the radial coordinate r, and e Λ ̸ = 1.The above-mentioned differential equation ( 21) gives the following solution, Here A stands for integration constant.In exploring the properties of wormhole solutions, we start by considering a metric that is spherically symmetric, and time-independent (static).This metric substitutes the standard quantities of GR denoted by e i for spherical symmetry in Eqn. ( 9), leading to the following line element, where A(r) ≡ e 2Φ(r) Π 2 (x) and B(r) ≡ . In this context, Φ(r) and b(r) serve as arbitrary functions of the radial coordinate and are called redshift and shape functions respectively.It should be noted that the metric potentials are influenced by the rainbow functions Π(x) and Σ(x).Nevertheless, it's worth mentioning that the coordinates r, t, θ, and ϕ remain independent of the energy probe particles in gravity's rainbow.
A fundamental feature of wormhole physics concerns the occurrence of a throat connecting two asymptotically flat spacetimes, with a minimum surface radius denoted by b(r 0 ) = r 0 .The radial coordinate covers the range between two patches, i.e. [r 0 , +∞), where r 0 stands for the transition point between the patches.Moreover, the flaring-out condition plays a major role in wormhole's geometry and can be stated as follows [1], In more precise terms, the shape function, indicated by b(r), should comply with the condition b ′ (r 0 ) < 1 to maintain the well-defined nature of the throat.What's more, an auxiliary condition is also necessary for the region lying outside the throat, namely This condition guarantees that the space-time in the region surrounding the throat remains traversable.Importantly, the redshift function called Φ(r), must remain finite throughout spacetime, thereby avoiding the appearance of event horizons in the wormhole geometry.
Our main concern is to develop a model that simulates a wormhole solution.In this regard, we generally employ an anisotropic fluid stress-energy tensor for describing the matter distribution, which is commonly described by, where u µ is the fluid 4-velocity denoted as u µ u µ = −1, χ µ is the unit radial vector satisfying the condition χ µ χ µ = 1 and g µν is the metric tensor.Furthermore, in this context, the energy density is represented by ρ = ρ(r), where r denotes the radial coordinate.Moreover, we have the radial pressure p r = p r (r) and the tangential pressure p t = p t (r), which depend solely on the radial coordinate, r.Now, let us express Eq. ( 4) in its covariant form.On the left-hand side, we construct the geometrical Einstein tensor.On the right-hand side, we introduce an effective stress-energy tensor. where Bearing in mind Eqs. ( 5) and ( 26), the non-zero components of the field equations ( 27) can be reduced to Here, ρ, pr , and pt stand for effective energy density, effective radial pressure, and effective tangential pressure respectively.One define, with When setting λ = 1 and σ = 1, the field equations ( 29)-( 31) can be modified to restore the standard GR definition for an anisotropic fluid sphere.In this modified formulation, there are five unknown quantities: Φ(r), b(r), ρ(r), p r (r), and p t (r), which are governed by three independent differential equations.Nevertheless, this system of equations is undetermined, implying that additional information or suppositions are required to completely capture the wormhole configuration.In our peculiar approach, we have chosen to concentrate on establishing the functions Φ(r) and b(r) as input variables.By doing so, we seek to get a somewhat more specific view of the wormhole.Lastly, it's worth mentioning the conservation equation of the stress-energy tensor T ν µ;ν = λR µ takes the following form

III. TRAVERSABLE WORMHOLE GEOMETRY
In order to construct a wormhole solution, we consider a simple form of redshift function Φ(r) that ensure the space-time is free from horizons and singularities as follows [76], Here ξ represents an arbitrary constant.By comparing equations ( 13) and ( 23), one obtains the following result, The wormhole shape function, denoted b(r), can be obtained by combining Eqs. ( 22), ( 36) and (37), respectively.This enables us to derive the following expression for the wormhole shape function, As stated by Morris and Thorne [1], in order to achieve a traversable wormhole solution, the shape function b(r) must satisfy the following crucial properties: The parameter r 0 corresponds to the wormhole throat radius, and the radial coordinate r satisfies the condition r 0 ≤ r ≤ ∞.To ensure the existence of a traversable wormhole solution, it is crucial that the shape function should satisfies these properties.These properties ensure that geodesics can pass through the throat without encountering singularities or obstructions.However, it's worthwhile to note that fulfilling these fundamental properties is necessary but not sufficient to establish the physical viability of a wormhole solution.Additional factors, such as the energy conditions, stability analysis and the behavior of matter fields, must also be taken into account to construct physically plausible wormholes.By substituting the throat condition b(r 0 ) − r 0 = 0 into Eq.(38), we arrive at the trivial solution r 0 = 0. To address this problem, we introduce a free parameter γ into Eq.(38).Accordingly, Eq. (38) has been modified to take the following form, Based on the information given by condition (2), we can infer the following, Upon replacing the value of A in Eq. ( 38), we can derive an alternative formula for the shape function, b(r).The resulting formula is as follows, where 0 < γ < r 0 .The prescribed range of γ meets both conditions (1) and (2).By using condition (3) with Eq. ( 41), we arrive at the following statement, Consequently, by evaluating Eq. ( 41), we get asymptotically flat traversable wormholes.To demonstrate each of the specified conditions, we have plotted a graph in Fig. 1, which clearly displays the fulfillment of all the necessary criteria by our generated wormhole.In addition, we have employed embedding diagrams to further visualize the wormhole geometry and extract valuable information.In the case of spherical symmetry, we are able to simplify the wormhole analysis by taking an equatorial slice with θ = π 2 .This choice capitalizes on the symmetry of the problem, in effect shrinking it to a two-dimensional slice.By holding the time coordinate at a constant value, t = const, we can investigate a snapshot of the wormhole geometry at a specific moment in time.This simplification allows us to concentrate on the spatial features of the wormhole without taking into account the dynamics related to temporal evolution.With θ = π 2 and t = const, the wormhole metric can be expressed in terms of radial coordinates r and azimuth angle ϕ.Consequently, the resulting two-dimensional geometry can be represented by a line element of the following form, With cylindrical coordinates (r, z, ϕ) in three-dimensional Euclidean space, we can visualize the line element (43).Once visualized, this involves embedding the geometry of the wormhole in the selected coordinate system.The embedding enables us to portray the curved space-time of the wormhole as a surface in three-dimensional Euclidean space.The specific form of the embedding, depending on the line element (43) and the peculiar wormhole model under consideration, can be stated by the following line element, In a three-dimensional Cartesian coordinate system, the embedded surface z ≡ z(r) allows us to visualize the form of the wormhole with axial symmetry as follows, By comparing Eqs. ( 43) and ( 45), we can easily determine that, It's worth mentioning here that Eq. ( 46) reveals an interesting characteristic: at the throat, the embedded surface becomes vertical, indicating that dz dr → ∞.This suggests a rapid change in the height coordinate as the radial coordinate approaches the throat, implying a significant curvature or narrowness in that region.Moreover, if dz dr tends to zero as r approaches infinity, it implies that the height coordinate approaches a constant value away from the throat, resulting in an asymptotically flat behavior.This behavior is indicative of a spacetime that becomes increasingly similar to flat spacetime at large radial distances.The embedded diagram shown in Fig. 1 provides a visual representation of the wormhole geometry.Significantly, this diagram visualizes both the upper universe for z > 0 and the lower universe for z < 0.

IV. ENERGY CONDITIONS
To ensure the physical consistency of solutions in the context of Einstein's equations, energy conditions (ECs) are imposed.These mathematical constraints serve as criteria for determining the validity of a spacetime geometry.The four primary ECs, namely the null energy condition (NEC), weak energy condition (WEC), dominant energy condition (DEC), and strong energy condition (SEC), play a crucial role in ensuring that the solutions satisfy certain physical requirements.By imposing these constraints on the energy-momentum tensor, we can assess the behavior of matter and energy within the framework of GR.According to classical GR, violations of the WEC are closely linked to the existence of exotic matter within wormhole structures.The violation of the WEC is characterized by the following condition where T ξϱ is the energy-momentum tensor, and u ξ represents any timelike vector.Hochberg & Visser [77,78], building onto earlier work by Morris & Thorne [1], extensively investigated the violation of the NEC in the context of wormholes employing exotic matter.They showed that the throat region of a wormhole does indeed violate the NEC.The behavior of the Ω expansion, w ξϱ rotation, and σ ξϱ shear of the congruence associated with the v ξ vector field is governed by the Raychaudhuri equations, which can be expressed as follows, In order for gravity to be attractive, the criterion is that R ξϱ v ξ v ϱ ≥ 0 holds true for all hypersurface orthogonal congruences, where w ξϱ ≡ 0 and σ 2 = σ ξϱ σ ξϱ .Mathematically, this condition can be expressed as, The ECs, including the NEC, WEC, DEC, and SEC, are derived from the Einstein field equations.These conditions impose inequalities on the components of the energy-momentum tensor, which correspond to energy and momentum densities.Mathematically, they can be written as, A violation of the ECs is a strong indication of the presence of exotic matter.More specifically, wormhole configurations necessitate a violation of the NEC.It should be noted that the NEC is the weakest restrictive of the ECs.It should be noted that the NEC is the least restrictive of the energy conditions.It has been shown that the stability of static wormhole configurations relies on the existence of exotic matter.Importantly, however, no experimental evidence has been discovered to prove the existence of exotic matter, raising doubts about the validity and physical achievability of wormhole.

V. EXACT SOLUTIONS FOR THE MATTER CONTENT OF THE WORMHOLE IN RASTALL-RAINBOW FORMALISM
A. Analyzing the properties of shape function Now, we discuss the recently derived shape function, denoted by Eq. ( 41).The zero tidal solution for the wormhole remains completely unaffected by the rainbow function Σ(x), pointing to its robustness independently of the specific form of Σ(x).It's also vital to ensure that the conditions 1 − b(r) r > 0 for r > r 0 , b ′ (r 0 ) < 1, and b(r) − b ′ (r)r > 0 are satisfied, as previously noted.Upon examining Fig. 1, it becomes apparent that for r > r 0 , the expression b(r) − r takes on negative values, indicating that b(r)/r < 1.This implies that with r ≥ r 0 and b ′ (r) < 1, the function b(r) − r exhibits a decreasing trend with respect to r, satisfying the necessary flaring-out condition.Consequently, the derived shape function accurately characterizes the structure of a wormhole.Furthermore, Fig. 1 illustrates the embedding diagram for wormhole, featuring both the upper universe (for z > 0) and the lower universe (for z < 0).This is also verified that the spacetime is asymptotically flat, meaning that dz dr tends to zero as r approaches infinity.

B. Analyzing the adiabatic sound velocity
The stability and physical significance of the solutions can be efficiently evaluated by considering the adiabatic sound velocity, expressed as, where < p > represents the average pressure across the three spatial dimensions, defined as, These fundamental properties hold true while respecting the following constraint, By carefully analyzing Eqs. ( 59), (60), and (61), we can extract the following, where, Thus, we see that sound velocity depends both on the parameters λ and Σ[x].Noted that the parameter λ must be constrained within the regions where both the flaring-out and stability conditions are simultaneously satisfied.Fig. 2 illustrates the parameter space, highlighting the valid region for satisfying these conditions.

C. Generating the exact solutions
Now, by substituting Eq. ( 41) into the Eqs.( 29)-( 31), the energy density and pressure components are readily derived, where, The correlation between energy density (ρ(r)), two pressure components (p r (r) and p t (r)) and the Rainbow function Σ(x) becomes obvious, highlighting the interdependence of these quantities.In Fig. 3, we depict the graphical behaviour of the energy density ρ through specific parameter values: Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r 0 = 1.05, respectively.This graph aims to display the validity of the inequality ρ ≥ 0 for the wormhole solution over a range of values of λ, in particular in the region of λ ∈ [0.4,0.8].However, it's interesting to note that there's an exception in the region λ ∈ [0.4,0.8] where ρ falls within the range [−0.00414, 0], which renders it invalid for the admissible range of λ.This is more precisely stated by saying that what kind of matter contains the Rastall-Rainbow wormhole solution, and thus we search for different ECs in the forthcoming section.This helps us to assess ECs and draw a meaningful conclusion.The derived formulas for NEC, WEC, DEC and SEC are based on the Eqs.( 59)- (61).
• Searching for NEC and WEC along the radial and tangential directions which are indicated by where, • Searching for DEC and SEC along the radial and tangential directions which are indicated by where,   Since, we know that WEC or SEC =⇒ NEC and DEC =⇒ WEC =⇒ NEC.Interestingly, all ECs are automatically violated if NEC is violated.Thus, our main focus lies on the NEC.In Figs. 3 to 8, we depict all ECs i.e., NEC, WEC, DEC, and the SEC along radial and tangential directions.Observing the Figs. 3 and 4, we see that the inequality ρ + p r < 0 holds at the throat of the wormhole and its vicinity only for λ ≤ 0.5.Similarly, the inequality ρ + p t < 0 also holds at the throat and its vicinity of the wormhole model for λ ≤ 0.5 and λ ≥ 0.6.From these observations, we can say that the model violates NEC throughout the spacetime and consequently the violation of WEC also in regions where λ ≤ 0.5 and λ ≥ 0.6.Interestingly, the DEC is also violated at the wormhole's throat and its neighborhood, as evident from Figs. 5, 6 and 7, respectively.On the other hand, the SEC remains positive for λ ≥ 0.5, indicating compliance with the ECs, as shown in Fig. 8.

VI. DISCUSSION AND IMPLICATIONS
The physical importance of Rastall-Rainbow gravity is profound regardless of any modified gravity theories.Since, the Rastall-Rainbow gravity model has been constructed with the combination of Rastall theory and the gravity's rainbow formalism.The present article aims to explore the possible existence of static and spherically symmetric wormholes in Rastall-Rainbow gravity.The modeling has been carried out using the well known embedding class I techniques to embed a four-dimensional space-time into a pseudo-Riemannian space-time.This procedure helps us to establish a connection between two metric potentials, and thus we are able to solve the modified field equations analytically.
More specifically, we obtained exact wormhole solutions by considering a specific choice of redshift function and solve the field equations analytically.We apply stringent constraints to our model parameters (λ, Σ[x]) from the flaring-out conditions and investigate the effects of both parameters on it.Notably, the resulting shape function obtained from the embedding approach is depending on both model parameters λ and Σ[x].Taking into account these constraints, we further investigate the energy conditions and showed that wormhole space-time violates the NEC at the throat and its vicinity, and consequently the violation of WEC also.Additionally, we investigate the overall stability of the configuration by determining the adiabatic sound velocity and explicitly computed the stable region, as shown in Fig. 2. Finally, we can say that the present approach may constitute a stable wormhole solution in Rastall-Rainbow gravity theory.

DATA AVAILABILITY STATEMENT
This manuscript has no associated data, or the data will not be deposited.(There is no observational data related to this article.The necessary calculations and graphic discussion can be made available on request.)

2 2 .FIG. 2 .
FIG. 2. The graph consists of two panels: the left panel displays a contour plot, while the right panel presents a threedimensional plot.It illustrates the behavior of the parameter space (λ, r), and highlights the region where the adiabatic sound speed condition is satisfied.The wormhole solution is characterized by specific parameter values, namely Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r0 = 1.05.

FIG. 3 .
FIG. 3. The graph consists of two panels: the left panel displays a contour plot, while the right panel presents a threedimensional plot.It illustrates the behavior of the energy density (ρ) for our wormhole solution with different values of λ in the range of λ ∈ [0.4,0.8].The wormhole is characterized by specific parameter values, namely Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r0 = 1.05.

FIG. 4 .
FIG. 4. The graph consists of two panels: the left panel displays a contour plot, while the right panel presents a threedimensional plot.It illustrates the behavior of the NEC (WEC) i.e., (ρ + pr) for the wormhole solution with different values of λ in the range of λ ∈ [0.4,0.8].The wormhole is characterized by specific parameter values, namely Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r0 = 1.05, respectively.

FIG. 7 .
FIG. 7. The graph consists of two panels: the left panel displays a contour plot, while the right panel presents a threedimensional plot.It illustrates the behavior of the DEC i.e., (ρ − pt) for our wormhole solution with different values of λ in the range of λ ∈ [0.4,0.8].The wormhole was characterized by specific parameter values, namely Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r0 = 1.05, respectively.

FIG. 8 .
FIG. 8.The graph consists of two panels: the left panel displays a contour plot, while the right panel presents a threedimensional plot.It illustrates the behavior of the SEC i.e., (ρ + pr + 2pt) for the wormhole solution with different values of λ in the range of λ ∈ [0.4,0.8].The wormhole was characterized by specific parameter values, namely Σ(x) = 0.95, γ = 0.95, ξ = 0.001 and r0 = 1.05, respectively.