Hyperbolically symmetric sources in Palatini f(R) gravity

A thorough examination of static hyperbolically symmetric matter configuration in the context of Palatini f(R) gravitational theory has been carried out in this manuscript. Following the work of Herrera et al. (Phys. Rev. D 103: 024037, 2021) we worked out the modified gravitational equations and matching conditions using the Palatini technique of variation in Einstein–Hilbert action. It is found from the evaluations that the energy density along with the contribution of dark source terms is inevitably negative which is quite useful in explaining several quantum field effects, because negative energies are closely linked with the quantum field theory. Such negative energies may also assist in time-travel to the past and formation of artificial wormholes. Furthermore, we evaluated the algebraic expressions for the mass of interior hyperbolical geometry and total energy budget, i.e., the Tolman mass of the considered source. Also, the structure scalars are evaluated to analyze the properties of matter configuration. Few analytical techniques are also presented by considering several cases to exhibit the exact analytical static solutions of the modified gravitational equations.


Introduction
The dilemma of present-day accelerated expansion of the universe, i.e., the dark energy problem is an elementary issue in theoretical physics of modern history. Certain cosmological observations influenced by the survey of type-Ia Supernovae demonstrated that the cosmological models which require self-deceleration of the universe under its gravitational pull are un-natural. Few astrophysical scenarios have convincingly proven that at some time in the past, our universe went through a changeover to an expansion stage of de Sitter type. a e-mail: mzaeem.math@pu.edu.pk (corresponding author) b e-mail: zeeshan.math@pu.edu.pk c e-mail: zohatariq24@yahoo.com The elucidation of all such unanticipated observational outcomes needs a detailed knowledge about the cosmological evolution which is not possible without Einstein's invincible theory of gravitation, i.e., general theory of relativity (GR). Despite the fact that GR produces precise results at every inch in the space, there exist some intractable issues like inflation, speedy expansion of universe, dark matter and dark energy. These have been a matter of interest of many scientists and urged them to curiously look for the solutions. This laid the foundation for the discovery of extended gravitational theories (also termed as modified gravitational theories (MGT)) that were supposed to satisfactorily elucidate the dark energy origin, the interaction between dark matter and dark energy and the late-time expansionary phase of our universe, thus paving a path for the forthcoming research. Unfortunately, no extended theory has appropriately elucidated the abovementioned cosmological mysteries, however, they have provided new insights into the constitution and structure formation of our universe from its birth to current evolutionary stage and other intricate issues involving the dynamics of evolution. The f (R) theory of gravity, being less sophisticated, straightforwardly generalizes GR by using f (R) (a generic function of Ricci scalar R) instead of R itself in the action principle. The Palatini version is taken into account, i.e., variation w.r.t. the metric tensor and the affinity (connection symbol) is deployed by considering them independent of each other.
The matter configuration inside a spacetime interrelates with its geometric curvature via field equations. To model a specific gravitational phenomena, the solutions to field equations are needed that can be acquired by imposing certain constraints e.g., symmetry assumptions, conformal flatness and other post-Newtonian approximations because of their high non-linearity. One out of such solutions is the renowned Schwarzschild solution for empty space and it is believed to be a unique static asymptotically flat solution to the field equations possessing spherical symmetry. From the word spherical symmetry, it is obvious that such a solution owns three killing vectors (mentioned in Eq. (2) in [1]), whereby invariance under rotations is depicted. In the form of polar coordinates, for R > 2M i.e. outside the horizon, one may express Schwarzschild metric as having the signature (−, +, +, +). For R < 2M, the metric tensors g tt and g rr interchange their signs causing the signature to become (+, −, +, +) which may depict a mathematical solution but physically it is insignificant. For R = 2M, an identifiable singularity emerges that can easily be eliminated by making use of coordinate transformation. A noticeable fact here is that the time-like killing vector outside the horizon turns into a space-like killing vector inside the horizon (implied by the fact that g tt and g rr interchange their signs). So, the coordinate transformation results in the loss of staticity inside the horizon. The staticity of the Schwarzschild solution over the whole spacetime is vital as it stems from the vacuum and sphericity conditions of the metric. So, in order to extend the solution from outside the horizon to its interior while maintaining the staticity everywhere, an absolutely different interpretation of horizon is requisite. Such an effort leads to a static solution that characterizes the spacetime as a combination of two absolute 4-manifolds: one representing the region exterior to the horizon and the other indicating the region inside the horizon. Both of these manifolds fabricate a single continuous topological manifold [2]. These θ − φ submanifolds possess spherical symmetry in the region exterior to the horizon and hyperbolic symmetry inside the horizon and meet each other i.e. R = 2M at a single curve θ = 0. Thus, for R < 2M, we acquire having the killing vectors as given by Eq. (4) in [1]. The motion of particles from exterior region of the horizon to its interior would be possible only when both of the manifolds intersect at θ = 0 which leads to the assumption of radial motion only. The plane as a whole is treated as a single mathematical manifold having a phase change at the point R = 2M because of signature change. Herrera et al. [1] propounded a comprehensive analysis on the motion of particles in the inner horizon of a blackhole with hyperbolic symmetry by investigating its geodesic equations in detail. The obtained outcomes are then compared with the previously obtained results for the case of spherically symmetrywhich revealed that a repulsive force in the horizon prevents the test particles from reaching the center. Herrera and Witten [2] put forward an alternative explanation of Schwarzschild black hole with a restriction that its solution must be static outside as well as inside the horizon. As an outcome of this restriction, they had to change the signature by implementing a transformation of one of the angle variables. The authors showed that a change of symmetry is producedwith a phase transition on the surface. Harrison [4] traced out the exact solutions of gravitational equations ofmotion by following themethod of separation of variables. Using these solutions along with the interchanging of variable technique, few other solutions are also obtained. All the solutions are then analyzed via geometrical and physical interpretations. Gaudin et al. [5] took into account a massless scalar field and investigated static solutions of Einstein's gravitational equations. Using some type of duality transformations, they developed a link of such solutions with the Kantowski-Sachs solutions. They also worked on the physical attributes of limiting case of a vacuum hyperbolic spacetime. The Schwarzschild solution indicating the exterior region of an astronomical configuration is widely recognized as an appropriate gadget for analyzing the gravitational field outside a spherically symmetric mass. The region interior to the horizon had no appropriate mathematical formulation to explicate it [3]. In fact, it was propounded at first by Harrison [4] and later became the subject of interest of many researchers [5][6][7][8][9][10].
Description of gravitational phenomena in competition to GR has been the most discussed topic for the past many years. Nojiri and Odintsov [11] overviewed the gravitational alternatives to GR i.e. the modified theories by considering various kinds of models in f (R), f (G) and f (R, G) gravity. They showed that few such theories went successfully through the Solar system tests and have excellent cosmological structures i.e. these are capable of describing the latetime epoch naturally. They also mentioned the possibility to elucidate the coincidence issue that manifests the expansion of the cosmos. Meng and Wang [12] investigated the Palatini formalism to modify the gravity by utilizing an ln R term first introduced by Nojiri and Odintsov [13] in the action integral. They demonstrated that in Palatini formalism, the ln R term can explain the current exponential expansion and reduces to Friedmann evolution for the region having high red shift. They also discussed the issue associated with the extension of ln R gravity by the addition of R m terms. Bertolami et al. [14] derived equations of motion for highly massive particles in the context of f (R) gravity theory. They demonstrated that an extra force is generated when an explicit coupling between lagrangian matter density and arbitrary function f (R) is achieved. Such an extra force is found to be orthogonal to the four-velocity. In weak field limit, they also obtained the corresponding law of acceleration. Paul et al. [15] considered a polynomial function of the curvature R known as the lagrangian density f (R) in the Einstein-Hilbert action for vacuum. By taking the Robertson-Walker metric, highly non-linear field equations are obtained that are not adequate for analytic solution, so the authors adopted numerical technique to analyze the evolution of FRW universe. The models worked out by the authors may help in reproducing the current expansion phase of the universe. Lobo and Oliveira [16] constructed traversable wormholes in the background of f (R) gravity theory. They imposed the condition that the matter in the wormhole fulfills the energy conditions which reveals that the higher order derivatives of the curvature terms violate the null energy condition depicting that the higher order terms support such wormhole geometries. Olmo et al. [17] provided an overview of the primary results such as the Tolman-Oppenheimer-Volkoff equation, mass-radius relations and Chandrasekhar's limit etc. in the background of modified theories of gravity. They combined the constraints and outcomes of both relativistic and non-relativistic stars and built a bridge between the observations and efforts of researchers working on variety of theories. Clifton et al. [18] presented a complete survey on the modified gravity theories along with their cosmological influences. They also reviewed the efforts to construct a Parameterized Post-Friedmannian technique useful to constrain the deviations from GR.
Solutions to Einstein's gravitational equations possessing hyperbolic symmetry have been an important topic among the researchers. Herrera et al. [19] considered dissipative fluid with hyperbolic symmetry and studied their general properties and found that the energy density comes out to be negative. The region in the center contains no element of the fluid so this region is described with the help of a vacuum cavity. By taking into account the quasi-homologous regime, few exact solutions are procured analytically.
Riemann tensor can be divided into certain parts by using its dual along with four velocity vector to attain few tensor quantities which in turn produce few scalar functions called the structure scalars. Such scalars appear to be connected with the fluid variables and play a crucial role in understanding the dynamics of astronomical object. Lanczos [20] found that the antidual part of Riemann tensor can be reduced to a contracted tensor quantity. From their analysis, the procured third order tensor having sixteen components and found that such a tensor can assist in integrating the field equations of a quadratic action by reducing them from fourth to second order. Herrera et al. [21] inspected static axially symmetric matter sources by evaluating certain dynamical equations e.g the field equations and the conservation equations. The inhomogeneity factors with the help of structure scalars are identified and few exact solutions are procured analytically, one out of which elucidates an incompressible spheroid with pressure isotropy. Bhatti et al. [22,23] constructed the modified structure scalars for spherically symmetric dissipative system by taking into account the Palatini f (R) gravity theory from the point of view of a tilted observer with and without the presence of electromagnetic force. They evaluated certain dynamical equations including the Raychaudhuri equation and the Bianchi identities. With the help of theses dynamical equations, they analyzed the structure formation as well as the evolutionary stages of astronomical object. Bhatti and Tariq [24] used the technique of splitting the Riemann tensor into three parts and obtained few structure scalars in the context of GR in the presence of an electromagnetic force for a spherically symmetric dissipative fluid. The equations exhibit that the presence of electric charge results in an increase in the mass of spherical body.
A locally-defined set of vector fields that are linearly independent and orthonormal to each other are much more beneficial as compared to the coordinate basis. Such a set is termed as orthonormal tetrad and has been extensively employed in the literature. Recently, Bhatti et al. [25] pursued a comprehensive study on static hyperbolically symmetric sources by extending the work of Herrera et al. [26] to the case where the electromagnetic force is non-zero. Writing down the gravitational equations, they worked out few analytical solutions along with their corresponding generating functions. Tariq et al. [27] explored the pressure isotropy condition for a spherically symmetric dissipative astronomical system in the background of Palatini f (R) gravity by working out a differential equation for Weyl scalar. They demonstrated that the factors like energy density, shear stress and dissipative flux are responsible for inducing the pressure anisotropy in the system. Herrera et al. [28] elucidated the changeover of an axially symmetric, radiating and dissipative fluid to its nonradiating state and evaluated the circumstances under which such a system reaches its equilibrium state. Herrera et al. [29] carried out a detailed investigation on the collapse of axially symmetric sources in GR by deploying all the basic equations. The equations represent a general dissipative fluid with pressure anisotropy. The techniques used by the authors have both analytical as well as numerical applications. Structure scalars are also identified and their role in the fluid dynamics is demonstrated. Herrera et al. [30] analyzed the exterior region of an axially symmetric source occupied by a null fluid produced as a result of dissipation process. The part played by the vorticity and its connection with the radiations is also determined.
This manuscript is assembled in the following pattern: Sect. 2 provides a bird's-eye view of the Palatini technique to arrive at the gravitational field equations. Further, it states the line element for hyperbolical symmetry along with the source, the orthonormal tetrad and the matching conditions. Section 3 incorporates the Palatini field equations, the hydrostatic equilibrium condition and the interior mass of the hyperbolic source. Section 4 comprises an expression for the intrinsic curvature tensor i.e. the Riemann tensor and an important relation between the Weyl scalar and the fluid variables. Section 5 includes an expression for the Tolman mass which is a measure of total energy budget of the source. Uti-lizing the breakdown technique of the Riemann tensor, Sect. 6 expresses certain modified scalar functions called the structure scalars. Section 7 demonstrates certain static solutions for the considered hyperbolic symmetry by taking different cases along with their sub-cases. Section 8 winds up the findings of the manuscript with a short summary and the values for Greek letters used in the dynamical equations representing dark source elements are provided in the Appendix.

Modified action principle and hyperbolically symmetry source
Modified gravitational theories come in handy when the gravitational field at large scales does not behave well in GR which is derived form the variation of standard Hilbert-Einstein action. In this scenario, it is found from the observations that a more general action is requisite to fully comprehend the gravitational interaction. One of the possible approach to attain such a generalized action is associated with the modification of geometric portion of the Hilbert-Einstein action. This modification aids in unraveling various cosmological mysteries including the dark matter problem. The most convenient and unsophisticated modification of the geometric portion is the substitution of Ricci scalar R with f (R) (a generalized function of R that contains higher order curvature invariants). Based on this generalized function, this theory is named as f (R) gravitational theory and its principal constituent is that extra degrees of freedom can be accumulated into an effective form of stress-energy tensor which could generate dark source effects. Following is the expression for such a modified Hilbert-Einstein action with S M nd S G signifying the matter lagrangian and gravitational part of action respectively To extract the modified field equations out of it, three different procedures exist in literature i.e. the metric approach, the Palatini approach and the metric-affine approach [31]. The Palatini approach relies on the fundamental idea that the connection α πν is independent of the spacetime metric g πν i.e the geodesic structure of spacetime is dissociated from its metric structure. The implementation of such a variation leads us to the following modified field equations where ∂g πν . Variation of Eq. (2) for independent Levi-Civita connection yields A connection between the trace of stress-energy tensor T and Ricci invariant R is as follows For vacuum, the Ricci invariant holds a constant value that can be defined using f (R). So, we may use another metric tensor notation i.e. h πν = f R g πν from where cosmic attributes of vacuum space in the context of f (R) theory can be observed [32]. Mathematically, the following relation may be established between this new form of metric tensor and the connection Making use of Eqs. (3) and (4), we acquire with ∇ σ and signifying the covariant derivative and the d'Alembert operator. Concisely, we can write this equation in the form of a single set of modified field equations as The efforts to determine a matter configuration sourcing the metric defined in Eq. (1) revealed that this metric basically indicates the interior region of the horizon which, in accordance with the observations, can possess hyperbolic symmetry only. We analyze such a hyperbolic distribution of matter that may or may not be bounded by the surface e from the exterior but needfully bounded from the interior by the surface i . The line element is as follows with spacetime variables ν and λ dependent on the radial coordinate only. The configuration of matter and energy indicated by the two-rank symmetric tensor T δζ is given as with μ and δζ indicating the energy density and the anisotropic pressure respectively and V δ is the fluid fourvelocity defined in covariant notation as (e ν/2 , 0, 0, 0). If orthonormal bases exist at some point, it provides an aid in specifying a locally inertial frame to express the components of tensors as observed by an observer at rest in that frame. The observer then possess a set of space axes (orthonormal) provided by spatial vectors e η . The combination of spatial axes with the observer's 4-velocity yields orthonormal tetrad. This observer's tetrad assists in splitting various tensors into their space-like and time-like parts and in gravitational theories, it helps in drawing out physical quantities from geometric objects. Following [26], the canonical orthonormal tetrad for our case is defined as which (after some manipulations provided in [26]) transforms the stress-energy tensor as Here, P x x , P xy , P yx , P yy and P zz illustrate the value of stresses in the directions indicated by their subscripts. Furthermore, the anisotropic tensor takes the below-mentioned form where P x x − P zz = . The isotropic stress has the form P = P xx 3 + 2P zz 3 . Since the interior region becomes a hyperboloid of two sheets [2], it is evident that the matter particles have no access to the region at the center of symmetry. Thus, a vacuum cavity emerges around the central region which demonstrates that the interior region becomes bounded too. The surface that bounds the matter from outside the horizon is represented as e having equation r = r e = constant whereas the surface that bounds the matter from inside is represented as i having equation r = r i = constant. The junction conditions for both the surfaces imply

Modified gravitational equations
Palatini f (R) equations acquired by substituting the nonzero components of tensors in Eq. (5) take the following 8π 8π The conservation equation law i.e. T δ β;δ = 0 generates the hydrostatic equilibrium equation as follows The expression for Misner-sharp mass function m = m(r ) for the interior region of the horizon is Using this value in the field equation associated with the energy density i.e. Eq. (9), we can relate the energy density with the geometric mass within the horizon as This expression depicts that geometric mass is negative which is unphysical because it is absolute i.e. it can never admit negative values. To avoid such a conflict, we have to assume negative energy density which is physically possible (dark energy is an example). So, we replace −μ by |μ| which yields the following outcome The radial derivative of the spacetime variable ν turns out to be equal to Equation (12) on substituting the value of dν/dr becomes Physically, this expression depicts that in the context of Palatini f (R) gravity, the pressure gradient, the passive gravitational mass density (P r − |μ|) along with the active gravitational mass, the pressure anisotropy and the dark source constituents must cancel the effect of each other to make the conservation of mass and energy possible.

The intrinsic curvature tensor
A symmetric trace-free tensor quantity of rank two defined by E δζ = C δγ ζ π V γ V π is termed as an electric part of Weyl tensor. It aids in the comprehension of tidal forces acting on an astronomical object in the manifold. Using the orthonormal tetrad, we may express it as E δζ = K δ K ζ + h δζ /3 with being equal to From the modified gravitational equational equations (9)-(11) and Eqs. (13) and (18), one may acquire 3m Evaluating the radial derivative of this expression along with the usage of Eq. (14), we attain Evidently, the dark source elements have a direct relation with which depicts that such elements give rise to the gravitational pull of the astronomical object which, in turn, results in the greater tidal force generated by it. For the mass function, we find When compared with Eqs. (43) and (44) of [26], it is evident that the curvature amendment causes an increase in the tidal force as well as the geometric mass of the object under consideration. For a specific f (R) model, one may be able to describe the extent to which the increment is produced e.g. for the Starobinsky model f (R) = R + α R 2 where the curvature amendment is quadratic in R, the dark source elements would involve curvature invariants up to second order for the case of Palatini f (R) and fourth order for the case of metric f (R) gravity.

Total energy budget of the source
The mass function describing the total mass and energy content of the source is called as the Tolman mass i.e. the active gravitational mass. Its mathematical expression mentioned in Eq. (45) of [26] produces the following result for the case of hyperbolic symmetry in the background of Palatini f (R) gravitational theory.
With the aid of Palatini field equations, we find which can be integrated to acquire an interesting relation connecting the Tolman mass with the spacetime variable ν in the presence of dark source elements as follows Using Eq. (16), with From Eq. (14), it is evident that we avoided the negativity of Misner-sharp mass by replacing −μ with |μ|. In the same manner, 4πr 3 P R − ξ 2 2 must be less than m f R to obtain the negative value for Tolman mass function which would indicate the repulsive nature of gravitational field around the considered astronomical object.
To comprehend the kinematics of fluid particles, let us evaluate the magnitude of 4-acceleration of the fluid configuration defined by a δ = V δ;β V β . In the context of Palatini f (R) gravity, it turns out to be equal to Relating it with the Tolman mass, we can write which points out that if 4πr 3 P R − ξ 2 2 < m f R , the acceleration of the fluid particles act radially inward. This provides additional support to the fact that gravitational force in our case has repulsive nature. Employing the r -derivative of Eq. (23) and using Eq. (25), we find Completing the derivative on the L.H.S by multiplying by r 4 and integrating w.r.t r , we get To relate the Tolman mass with the fluid variables i.e. the energy density inhomogeneity and pressure anisotropy, we substitute the value of in the above expression and find In comparison with Eq. (54) of [26], the above expression depicts that the dark source elements also affect the total mass and energy distribution of the system. Just like the geometric mass, the Tolman mass is also affected by including higher order curvature invariants in the theory.

Modified structure scalars
The intrinsic curvature tensor i.e., the Riemann tensor has four indices. When acted upon by the Levi-Civita tensor, it yields its dual which can further be contracted with a unit time-like vector (4-velocity in our case) to obtain a new tensor quantity having two indices. Such formalism results in the resolution of Riemann tensor into its electric and magnetic parts which have great significance as they lead to certain scalar functions called the Structure Scalars. The mathematical procedure is well-defined by Herrera et al. in [33] and utilized by Bhatti et al. in [22] for the case of Palatini f (R) gravity. So, we have The tensors in Eqs. (27) and (28) can further be decomposed into trace and trace-free parts as follows and where It is clearly visible that each scalar function is associated with some characteristics of fluid configuration, for instance, the structure scalar X T is associated with the description of energy density of the fluid in addition to the dark source elements. Using Eq. (20), the trace-free scalar Y T F , comes out to be directly related with the energy density inhomogeneities and dark source elements as follows Similarly, for Y T F , we find Furthermore, the combination of trace-free scalar functions X T F and Y T F is found to be associated with the local stress anisotropy with the addition of dark source elements i.e.
On comparing Eqs. (30)- (36) with Eqs. (69)-(74) in [26], it becomes evident that the curvature amendments affect the structure scalars too. To explore the scalars Y T and Y T F , we get back to the expressions for Tolman mass in Eqs. (22) and (26). Using Eqs. (32) and (36), we find Evidently, the total mass inside the fluid (bounded between the radius r and the surface r e ) can be expressed by using the structure scalar Y T F . This scalar explains all the basic characteristics of the fluid configuration i.e the energy density inhomogeneity and the anisotropic stress that generate complexity in a system. For this reason, it can be termed as the complexity factor. Y T is directly linked with the Tolman mass density.

Few static solutions
This section copes with the issue of determining all the possible static hyperbolically symmetric solutions i.e. the values for the spacetime variables in the context of Palatini f (R) gravitational approach. Using the format of [26], we first develop a general formalism and then implement certain constraints to acquire the possible solutions to the modified gravitational equations. From Eqs. (10) and (11), which on using new variables e −λ = y, z = ν /2 + 1/r becomes Solving this first order differential equation and putting back the value of y, we find e λ(r) = z 2 e 2z+ 4 zr 2 dr r 6 2 z r 8 8π where A denotes the integration constant. This formalism reveals that all the possible solutions can be well-described with the aid of two generating functions z and . So, for the fluid variables, we acquire Now, we proceed towards the determination of certain explicit solutions to modified gravitational equations along with their generating functions.

Conformally flat constraint
The simplest constraint that may allow a static solution to Eqs. (9)-(11) is the conformally flat constraint. It is based on the hypothesis that the Weyl scalar vanishes i.e. = 0. This constraint when used in Eq. (18) allows the completion of derivatives which in turn produces the following result ∂ ∂r which on using e −λ = y and ν = 2w /w can be rewritten as It is a first order differential equation in y. Its solution can be determined as where B emerges as the integration constant and the values forf 1 (r ),f 2 (r ) are given as In terms of spacetime variables, we can rewrite the obtained solution as with γ symbolizing the constant of integration. Utilization of junction conditions assists in the evaluation of γ as Usage of Eq. (39) generates the following outcome where σ appears as integration constant. Upon using junction conditions, σ is evaluated to be equal to So, the conformally flat condition resulted in the evaluation of only one generating function. To determine the complete solution, assumption of another constraint is requisite. For this reason, we consider vanishing radial pressure constraint i.e. P r = 0. Then Eq. (10) implies Employing both the constraints (P r = 0 = ) in Eq. (18), we find which on assuming 2g = e −λ + 1 generates the following result Integrating w.r.t radial coordinate, we determine the following relation The fluid variables turn out to be equal to For e ν to be positive, g > 2/3 which coincides with the result of GR (Eq. (101) in [26]). This leads to a possibility that a thin shell emerges at the boundary r = r min in Palatini f (R) gravity too. The generating functions , z are determined as follows

Solution obtained from zero complexity factor constraint
Just like GR, the scalar Y T F in Palatini f (R) determines the complexity generated in the fluid configuration. So, it would be beneficial to look for a model for which the complexity factor vanishes i.e. Y T F = 0. Specifically, we need to assume an additional constraint besides this to obtain a realistic mathematical model. For this purpose, we take P r = 0. Thus .
The implementation of vanishing complexity factor leads to the following outcome Using together the above two equations along with Eq. (24), the spacetime variable e ν is found to be equal to The contribution of curvature amendments is clearly visible. For f R = 0 and f (R) = R, the result reduces to the findings of GR given by Eq. (129) in [26]. Employing Y T F = 0, we determine the following expression in terms of g Completing the derivative and integrating, we obtain withC 2 appearing as integration constant. The fluid variables are determined as which clarifies that the fluid configuration in Palatini f (R) gravitational theory is also restricted by the minimum value of the radial coordinate fulfilling the constraint g(r min ) < 1 (because otherwise the fluid variable would become undefined). Such a minimum value i.e. r min can be determined from Eq. (44). The existence of vacuum cavity, a thin shell and the negativity of Tolman mass follow in the same way as defined for GR [26] with a significant difference that r min determined from Eq. (44) involves the contribution of dark source elements too. This model possesses the generating functions defined as follows It is evident that the curvature amendments affect the generating functions of the model too.

Solution obtained from stiff equation of state
To determine a mathematical model representing a realistic astronomical object, a thermodynamical equation called the equation of state (EoS) is used. One of its types i.e. the stiff EoS is used for the fluid that is harder to compress and requires more support against the gravitational force e.g. Neutron stars. We utilized it here because it may provide of the possible solutions to Palatini field equations as it did for GR. The constraint that is followed by such EoS is the equivalence of energy density and radial pressure i.e. |μ| = P r . Hydrostatic equilibrium condition under this constraint takes the following form To determine the generating functions, we need another constraint except |μ| = P r because it is insufficient to obtain a solution to highly non-linear modified field equations. This calls for the assumption of two additional constraints mentioned below.
• Vanishing tangential pressure i.e. P ⊥ = 0 Considering P ⊥ = 0, the value of P r is determined by integrating Eq. (45) and producing Pr dr , withK expressing the integration constant. This expression demonstrates that the energy density as well as the radial pressure are influenced by the higher order curvature invariants of the theory. The mass function and the spacetime variables under these circumstances turn out to be as follows In GR, the spacetime variable ν for such a model came out to be constant which represented the vanishing of the active gravitational mass. Here, we found that ν is not constant rather it is a function dependent on the dark source elements. This depicts that for the same mathematical model (obtained by using the constraints |μ| = P r and P ⊥ = 0) under the influence of Palatini f (R) gravity, the strength of the gravitational field generated by the astronomical body is non-zero and depends on the force exerted by the dark source constituents only. So, for Palatini f (R) gravity, the active gravitational mass is non-zero and the passive gravitational mass density i.e. P r −|μ| vanishes identically since we are working with the constraint P r = |μ|. • Vanishing complexity factor i.e. Y T F = 0 When stiff equation of state and the vanishing complexity factor constraint are used simultaneously, the hydrostatic equilibrium condition yields with ω (D) symbolizing the dark source elements. The solution to this second order differential equation is Clearly, the dark source constituents add to the stress acting in the radial direction.ã andb correspond to positive constants of integration. The mass function for such a model turns out to be Using this equation and Eq. (13), one can easily determine the spacetime variable λ. The value for other spacetime variable ν can be found by substituting these results in Eq. (15). Further, we assume that the fluid configuration is bounded from the exterior region by the surface e to obtain It is evident that the active gravitational mass and the passive gravitational mass density are negative and zero respectively for this kind of model in Palatini f (R) gravity.

Discussion
The Palatini f (R) gravity is quite straightforward and uncomplicated generalization of GR, because it renders GR without requiring any kind of link between the metric and the connection. When tested via solar system tests, its results were found to be consistent with the observations. The fascinating mysteries of the universe like blackholes, gravitational waves and wormholes also exist in this theory [34][35][36]. It incorporate no additional degrees of freedom instead, the dynamical equations are changed by the addition of new terms indicating the curvature amendments of the theory. Using a specific type of model, such curvature amendments provide an aid to comprehend the dark matter or dark energy epoch of the universe. The physical and mathematical models are modified by the curvature invariants involved in the theory to describe the expansionary phase of universe. One of the significant advantages of using Palatini f (R) gravitational theory is that it works well also for the models that show certain inconsistencies when utilized with the metric version of f (R) theory. An example is the cosmic speed-up problem modeled by f (R) = R − μ 4 /R which when used with metric f (R) gravity exhibits certain deviations from the observations, but when used with Palatini version, generates coherent results.
We have conducted a comprehensive study on the characteristics of hyperbolically symmetric fluid configurations with anisotropic stresses. The findings are parallel to the case of GR [26] where the manifold represents a 2-surface sphere from the exterior and hyperboloid of two sheets from the interior 2-surface. It is found that the fluid particles cannot reach the central region which results in the formation of a vacuole in that region. Negative energy density is requisite for the mass function to be positive in Palatini f (R) gravity too. The contribution of dark source constituents can easily be seen in the expressions of the Weyl scalar and the mass function which demonstrates that the presence of dark source influence the tidal forces as well as the mass of hyperbolically symmetric astronomical object. The Tolman mass exhibits that in addition to the baryonic matter, the fluid is composed of dark source elements also. Few feasible solutions to such a type of symmetry in the manifold are then worked out, each one of which is defined by a different set of constraints to deal with the non-linearity of the modified gravitational equations. The obtained models resemble GR models with a considerable difference that their dynamical equations exhibit a physical behavior influenced by the higher order curva-ture invariants of the theory. All the outcomes are consistent with the findings of GR under the restriction f (R) = R and f R = 0.