Generalized Lemaítre–Tolman–Bondi spacetime under the influence of electric charge and Palatini f(R) gravity

The objective of this article is to investigate the effects of electromagnetic field on the generalization of Lemaître–Tolman–Bondi (LTB) spacetime by keeping in view the Palatini f(R) gravity and dissipative dust fluid. For performing this analysis, we followed the strategy deployed by Herrera et al. (Phys Rev D 82(2):024021, 2010). We have explored the modified field equations along with kinematical quantities and mass function and constructed the evolution equations to study the dynamics of inhomogeneous universe along with Raychauduary and Ellis equations. We have developed the relation for Palatini f(R) scalar functions by splitting the Riemann curvature tensor orthogonally and associated them with metric coefficients using modified field equations. We have formulated these scalar functions for LTB and its generalized version, i.e., GLTB under the influence of charge. The properties of GLTB spacetime are consistent with those of the LTB geometry and the scalar functions found in both cases are comparable in the presence of charge and Palatini f(R) curvature terms. The symmetric properties of generalized LTB spacetime are also studied using streaming out and diffusion approximations.


Introduction
The recent astrophysical observations indicate the accelerated rate of universe expansion but general relativity (GR) is unable to explain it properly. The existence of an unknown form of matter is suggested by gravitational lensing in the universe named dark matter (DM) while the accelerated expansion is caused due to the existence of dark energy (DE). Einstein's theory is also unable to illustrate the flat rotation curves of galaxies as well as the expansion rate of the universe without invoking DM and DE, which together account for a e-mail: mzaeem.math@pu.edu.pk (corresponding author) b e-mail: zeeshan.math@pu.edu.pk c e-mail: furqanhussainmpf27@gmail.com 96% of the universe's total energy content. Modified gravity theories have been discovered to be one of the most appealing mathematical methods for exploring the universe's rapid expansion. These theories were recently reviewed in the literature to explain various dynamical phases of the universe evolution, particularly the early and late time acceleration eras determined from recent observations. By changing the gravitational part of the standard Einstein-Hilbert action, various dark energy theoretical models are proposed. The f (R) theory is the most fascinating because it offers a viable alternative to DE without taking into consideration DM. The inclusion of the generic function f (R) in the standard Einstein-Hilbert action is the simplest modification of GR knows as f (R) gravity and applying variational principle to the modified Einstein-Hilbert action gives rise to different formalisms of f (R) gravity. The formalism in which we consider the metric and connection as an independent quantity is known as Palatini formalism while the other is known as metric formalism. The Palatini formalism produces a second-order nonlinear differential equation but in the case of metric formalism, we have fourth-order field equations. Some suitable forms of the f (R) model have been developed to describe the accelerated expansion by SnIa data at late times. The Palatini version of f (R) theory generates second-order field equations and present a simple gravitational alternative in the absence of any matter content. Furthermore, the significant advantages of this version are, it provides correct Newtonian limit, it is free of instabilities. Also, the Palatini f (R) vacuum solutions are comparable with GR solutions with a cosmological constant. It is observed that by applying Palatini formalism to the f (R) theories, we can create models with the radiation-dominated, and late-time acceleration phases that are compatible with the recent observations.
The modification of standard Einstein-Hilbert action in order to study the f (R) gravity is manifested as [2] S f (R) = S m + 1 2κ here S M and κ are the matter action and coupling constant and f (R) denotes the non-linear Ricci scalar function, while g represents the determinant of g αζ . The contraction of Ricci tensor R αζ associated with connection symbols is used to construct R in the modified action, meaning that R is the product of geometrical connections. One can get the following equations of motion by applying the variation on the action integral given in (1) with respect to connection symbols μ αζ and metric tensor g αζ as ∇ μ (g αζ √ −gF(R)) = 0, where T αζ represents the standard energy momentum tensor that is independent of connection symbols and F(R) = The equation that compares R (Ricci scalar) with T (Trace of the energy momentum tensor) is obtained by multiplying (2) with conjugate metric tensor g αζ as tells us that the energy-momentum tensor is dependent on R.
For second-order field equations, a standard expression can be obtained by solving (2) and then utilizing the obtained value in (3) of the following form The f (R) field equations can be expressed in alternate form as where G αζ ≡ R αζ − 1 2 g αζ R and is the d Alembert operator and can be defined as ≡ g αζ ∇ α ∇ ζ whereas ∇ ζ indicates the covariant derivative. Moreover is the Palatini f (R) energy-momentum tensor that describes a gravitational contribution. Bhatti et al. [3] examined the stability of charged neutron stars in the presence of modified gravity and constructed the generalized Tolman-Oppenheimer-Volkoff (TOV) equations. They deduced that if uniform entropy and chemical composition are taken under consideration then generalized TOV equations represent stable star configurations that are independent of the generic feature of modified gravity. With the Palatini approach of f (R) gravitation, Kainulainen and co-workers [2] determined the mass equation for the interior geometry and structure of stars. They demonstrated that in general, combining the interior and exterior geometry yield a connection between the density profile of star and gravitational mass that differs from the one found in general relativity. Bamba et al. [4] investigated the thermodynamics of the apparent horizon using equilibrium and non-equilibrium constraints in f (R) gravity. They concluded that understanding the horizon entropy in the equilibrium system is more straightforward than in the non-equilibrium one. They also exhibited that the second law of thermodynamics could be verified from both phantom and non-phantom stages as long as the temperature of the universe outside and inside the apparent horizon stays unchanged. Sharif and Yousaf [5] examined the stability of spherical collapse using the Palatini approach of f (R) gravity. They used Bianchi's identities of energy-momentum tensor to construct the collapse equation and investigated the instability ranges in the Newtonian and post-Newtonian patterns. They deduced that adiabatic index 1 managed the stability of the star.
Bhatti and Yousaf [6] determined the inhomogeneity factors using the dissipative fluid under the influence of electromagnetic field in Palatini f (R) gravity for plane topology. To examine the corresponding inhomogeneity causes, some specific cases are examined with and without dissipation. They investigated anisotropic, dust, and isotropic matter in the presence of charge in a non-radiating scenario. They examined the inhomogeneity factor for a dissipative fluid in the presence of a dust cloud. They inferred that the electromagnetic field raises the matter inhomogeneity, while the additional curvature terms make the system more homogeneous as time passes. Yousaf et al. [7] studied the consequences of charged gravastar models in f (R, T ) gravity using metric formalism. They used the equation of state to investigate the exact models and stable regions of gravastar under the influence of electromagnetic field and taking into account the Reissner-Nordström spacetime to match the exterior geometry with the interior at the hypersurface. They also examined the realistic character of charged gravastar model using a graphical representation. Bhatti et al. [8] examined the impact of density inhomogeneity and local anisotropy on the gravitational collapse in f (R) gravity under the influence of charge. They explored modified Maxwell field equations and established a connection between Weyl tensor and local anisotropy. They acquired an expression for Tolman mass which plays a significant role to study the impact of charge and extra curvature terms on the physical variables, as well as the study of their role in gravitational collapse.
Bhatti and Tariq [9] studied the dynamics of an anisotropic sphere in the presence of charge within the context of f (R) gravity. They used the general pattern provided by Herrera and Barreto to examine the primitive properties of polytropic spheres and constructed physical restrictions for polytropes. They tested the consistency of polytropes through Tolman mass by using conformally flat conditions. They inferred that the electromagnetic field plays the same function as anisotropic pressure. Tremblay and Faraoni [10] investigated the initial value problem of f (R) gravity using the dynamical correspondence with Brans-Dicke gravity for metric and Palatini formalisms. They concluded that the Cauchy problem is well-defined for metric formalism as compared to Palatini formalism under the influence of matter. Sharif and Siddiqa [11] studied the self-gravitating source under the influence of charge in f (R, T ) gravity explaining the situation of expansion and collapse process. They calculated the expansion scalar to examine the process of collapse and expansion and discussed the impact of charge and model parameters on physical variables with the help of graphical representations. They noticed that physical variables changes with time in case of expansion while the charge effects during the collapse process and the model parameters have the same effects in both cases. Garfinkle [12] examined that for type Ia supernova data inhomogeneous spacetimes are used as a cosmological model. They inferred that the model is found to match the data as well as the regular CDM cosmology with certain choices.
Herrera et al. [13] investigated the physical significance of structure scalars for dissipative and dust fluids under the influence of charge. They highlighted the importance of these factors in the scalar functions and addressed their physical outcome. They inferred that the evolution of kinematical variables and modifications in the inhomogeneity factor provided by the above-described factors requires special attention. Sharif and Yousaf [14] examined the dynamical features of cylindrical self-gravitating structures by considering the f (R) model undergoing the process of dissipative collapse and observed the consequences of higher-order curvature terms appeared in the scalar functions. They used the transport equation to study the impact of relaxation time on the evolution of collapsing structure and developed eight modified scalar functions, three Weyl, and two shear scalars that are needed to address the dynamics of cylindrically symmetric spacetime, in contrast to spherical systems. Sharif and Bhatti [15] developed the structure scalars for planer symmetry in the presence of charge and associated them with the fluid variables. They established a connection between Weyl scalar and kinematical quantities with the help of Taub mass formalism and discussed the consequences of the charge on the evolution of matter variables. Yousaf [16] investigated how f (G, T ) terms affect the formation of scalar functions as well as evolution equations. They dealt with a geometry linked with viscous dissipative matter content and established the relationship between Misner-Sharp mass and curvature scalar. They also developed the relation for f (G, T ) scalar functions by decomposing the Riemann tensor and investigating their consequences.
LTB models are the most ancient and fascinating Einstein equation solutions. They define a spherically symmetric distribution of non-dissipative dust that is inhomogeneous and have been used in the study of quantum gravity and gravitational collapse as a cosmological model. LTB spacetimes have the drawback that they do not allow dissipative fluxes. Yousaf et al. [17] examined the evolution of LTB spacetime under the influence of charge in Palatini f (R) gravity. They established the relationship between comoving and tilted observers by considering the effects of charge for dust and imperfect fluid respectively and also studied the f (R) models in LTB geometry in the presence of dissipative fluid. They inferred that dynamics of LTB spacetime have been modified due to the existence of charge and Palatini f (R) correction terms. Bhatti and Yousaf [18] analyzed the impact of f (R) terms for the collapse of LTB spacetime in the presence of the charged field. They examined the apparent horizons to study their consequence during the process of collapse and noticed that the area of horizons reduces by the presence of charge. They also inspected that the impedance in the gravitational collapse is caused by the presence of charge due to which the process slows down and observed that due to a significant role of higher-order terms in the collapse process, the system becomes more consistent. Chakraborty and Saha [19] discussed the dynamics of different horizons in the context of inhomogeneous LTB geometry. They analyzed that the comparison of these horizons reveals that the apparent horizon has a different character from the other horizons and if the apparent horizon accelerates, it will become a trapping horizon. They also described the Kodama vector with a causal character similar to that of the FRW model. Fernandes et al. [20] examined the behavior of LTB spacetime on the massive electromagnetic models and estimated the physical quantities like density and scale factor for Podolsky models and discovered a new type of singularity. Herrera et al. [21] applied the technique of Misner and Sharp to examine the gravitational collapse for dissipative fluids. They linked the transport equation with the dynamical equation in the background of Israel-steward theory together with thermodynamics coefficients and compare the results with previous calculations where these coefficients were ignored. Herrera et al. [22] discussed the properties of dissipative fluids under the influence of shear-free conditions and demonstrated the purpose of dissipative fluxes. They concluded that the magnetic part of the Weyl tensor disappears for geodesic fluid implying that gravitational waves do not emit during the evolution. In this manuscript, we investigate the generalization of LTB spacetime for charged fluid in Palatini f (R) gravity and consider the spherically symmetric spacetime along with dissipative dust fluid by following the strategy deployed by Herrera et al. [1]. In Sect. 2, we formulate the field equations coupled with f (R) gravity and matter variables in the presence of charge. Section 3 includes the structure scalars obtained from the orthogonal decomposition of Riemann tensor and evolution equations. In Sect. 4, we study the properties of LTB spacetime and analyze the behavior of the evolution equation by considering the masses of compact objects with the help of graphical representations. In the same section, we discuss the GLTB spacetime and its symmetric properties using diffusion and streaming-out approximations. Finally, our main results are mentioned in the last section.

Inhomogeneous LTB geometry coupled with kinematical variables
The diagonal non-static inhomogeneous and irrotational LTB line element is given by the following form where B is a dimensionless metric coefficient and C has similar dimension as r while both depend on t and r . The interior region of the metric is filled by dissipative dust fluid using heat and null radiation. The energy-momentum tensor describes such matter fields as follows where l α , μ, , V α , q α correspond to the null four-vector, energy density, radiation density, four-velocity and the heat flux respectively. The null four vector l α = δ α 0 + B −1 δ α 1 , four velocity V α = δ α 0 , heat flux q α = q B −1 δ α 1 , and unit four-vector χ α = B −1 δ α 1 , in a comoving coordinate system satisfy the following relationships The relativistic matter distribution is fully explained by matter variables like shear tensor, four-acceleration, and expansion scalar. The four-acceleration is defined as a α = V α;ζ V ζ and its non-vanishing components are given by as seen from (11), the occurrence of four acceleration components is completely dependent on Palatini f (R) formalism.
The shear tensor is defined as where h αζ = g αζ + V α V ζ is the projection tensor and the non vanishing elements of shear tensor in f (R) gravity are given by The shear scalar appeared in the above equation is defined as becomes The fractional rate of change in volume at the center of a spherical cloud to proper time is determined by expansion tensor and its trace is known as the expansion scalar and its value in Palatini f (R) gravity takes the form The Modified field equations in the form of Einstein tensor under the influence of electromagnetic field can be written as The electromagnetic tensor for the comoving coordinate system can be used to assess the influence of charge and defined as: where F αβ = −φ α,β + φ β,α is the Maxwell strength tensor that satisfies the Maxwell field equations.
where φ μ , j μ and μ 0 represents the four-potential, fourcurrent and magnetic permeability. For the solution of the Maxwell field equations, we take φ μ = φδ 0 μ , and j μ = σ U μ while charge density and scalar potential are denoted by σ and φ, respectively, which are functions of t and r . We obtain two independent equations while solving the Maxwell field equations of the following form Equation (19) is integrated with respect to r to provide where s = 4π r 0 σ C 2 B F 2 represents the total charge carried by non-interacting particles in a region of radius r . The non-vanishing components of the electromagnetic stress tensor are given by The Palatini f (R) field equations with charged field are given bẏ whereμ e f f , e f f andq e f f each contain the Palatini f (R) extra as well as usual terms whose values are given below , here χ 00 , χ 11 , χ 01 , and χ 22 are the extra curvature terms appeared in the field equations and are addressed in Appendix C, while prime denotes r differentiation and dot denotes t differentiation.
The formula for calculating the mass function m(t, r ) in the presence of charge to determine the amount of matter inside a self-gravitating system is provided by Misner-Sharp [23], which is given by so Eq. (28) can be rewritten as The radial variation of the mass function yields the following equation.
where U indicates the collapsing matter velocity and can be specified as U = D T C where D T = ∂ ∂t , it is assumed to be negative in the case of collapsing spheres. Whileq = q + ,μ = μ + . Integrating partially the Eq. (30), with a Palatini f (R) context reveals the relationship between matter variables and mass function given below The Weyl tensor can be manifested in the context of its electric part for spherically symmetric matter distributions while its magnetic part vanishes and the electric part is defined as in terms of unit four-vector, metric tensor and unit fourvelocity, the electric part of Weyl tensor may be granted as where is the scalar enwraping the products of spacetime curvature. Now, using Eqs. (23), (25), (26), (28) and (34) we have The equation describes the dependence of the Weyl scalar on the matter variables and the mass function under the influence of charge as well as dark source terms.

The matching conditions
We consider the exterior spacetime [24] given by the following equation in order to derive the junction conditions for exterior and interior spacetime.
where the symbols r and ν referred to the null and retarded time coordinates while M is a function of retarded time coordinate. Darmois junction conditions are used to match smoothly the interior and exterior geometry at the surface which produces where = signify that both quantities are calculated on the boundary surface.

Dynamical equations and structure scalars
Here, we define the tensors Y αζ and X αζ and dynamical variables in the context of f (R) gravity. We also manifested a connection between curvature scalar and fluid parameters in the presence of electromagnetic field and Palatini f (R) terms by formulating the modified Ellis equations. Bel [25] studied the orthogonal decomposition of Riemann curvature tensor R αγ ζ δ which has been applied to describe the evolution and structure formation of self-gravitating objects in the background of different Astrophysical scenarios. We introduce the tensors Y αζ and X αζ by adopting the pioneer work of Bel, which are the elements of the orthogonal splitting to obtain structure scalars in Palatini f (R) gravity as where the right and double dual of the Riemann tensor can be expressed as R * αζ γ δ = 1 2 η γρδ R ρ αζ and * R * αγ ζ δ = 1 2 η ρ αγ R * ρζ δ and the Levi-Civita tensor is denoted by η ργ δ . The tensors Y αζ and X αζ can be written as The tensors Y αζ and X αζ can be decomposed in their trace and trace-free parts. Each of these terms in its turn (in the spherically symmetric case) may be expressed through a single scalar function (named as the structure scalars). These variables were firstly introduced, defined and discussed in detail for the first time by Herrera et al. [26]. Following the same scheme as developed in [26], we can find Y T , X T , Y T F and X T F as where using Eqs. (31), (35) and (44), we obtain where the values of M 3 , are mentioned in Appendix A. The null radiation and dissipative variables as well as the correction terms can be used to express the scalar Y T F . In [26], it was also stated that there is a relation between the structure scalar Y T F and Tolman mass. In literature, it has been shown that for geodesic case, the scalar Y T F supervise the stability of the shear-free condition and the trace part of Riemann second dual tensor is found to be dependent on the energy density with extra curvature terms due to f (R) Palatini gravity, while the remaining scalar functions are dependent on the anisotropic stress tensor.
We can construct a couple of differential equations to continue our investigation and to examine the expansion and motion of particles by using the method addressed by Ellis [27]. These equations describe a relationship between fluid variables and the Weyl tensor along with Palatini f (R) terms and can be found aṡ (μ e f f + e f f )Ḃ B

3F
F Bq 2F B e f f + Structure scalars (43)-(46) for line element (8), in terms of metric coefficients are given by These structure scalars are interlinked with the simulation of compact objects in f (R) gravity and also in that scenario, we use the Raychaudhuri equation for investigating the consequences of R 2 terms. In the absence of Palatini f (R) extra curvature and charge terms, we obtain the dynamical variables as mentioned in [1] from the Eqs. σ The Eq. (57) explains the expansion evolution with the help of the f (R) structure scalar Y T with the addition of dark source terms while Eq. (58) describes the connection between the structure scalar Y T F and shear scalar σ e f f which illustrate the impact of electromagnetic field and Palatini f (R) gravity on the shearing motion. Moreover, the relationship between Weyl tensor and fluid variables in terms of dynamical variables can be found as Now, the Eq. (59) is reduced to the following form for nondissipative dust fluid ifμ e f f = 0 then Eq. (60) becomes integrating Eq. (61), which produces The energy density is being regulated by the structure scalar X T F for non-dissipative dust fluid as observed from the Eq. (62). This is also true for an isotropic fluid as described in [26].

Lemaître-Tolman-Bondi spacetime and its generalization
In this section, we address the properties of LTB spacetime and analyze the behavior of the evolution equation under the influence of charge by considering the masses of compact objects. We generalize our LTB spacetime by incorporating the heat flux and acquire structure scalars for that generalized spacetime. To find similarities between LTB and generalized spacetime, we match the structure scalars corresponding to LTB and GLTB. While on the other hand, we discuss symmetric properties of generalized spacetime based on the assumption that they correspond to LTB properties using two cases: diffusion and streaming out limit. We consider the nondissipative dust fluid to obtain the general mechanics for LTB spacetime. By substituting q = 0 in Eq. (24) and then integrating we obtain where κ serves as a constant for integration. We get the Lemaitre-Tolman-Bondi-metric in the following form by substituting the Eq. (63) in (8).
The LTB spacetime is generally coupled with an inhomogeneous dust source but it is important to mention that the general source corresponding to the metric appeared in (64) is an anisotropic fluid [29]. Now the element of Bianchi identities take the simplified form for non-dissipative dust fluid  Fig. 1 The plot is for κ < 0, m = 1.4 and s = 0.5 that corresponds to closed universe which explains that after expanding to maximum size, the universe falls into itself known as big cruch integrating Eq. (65), we get where D 5 = exp and h(r ) is an integration constant. Now using the Eqs. (28) and (63) we obtaiṅ where The structure scalars Y T F and Y T do not incorporate the term κ in their formulation as shown in Eqs. (69) and (70). Now again integrating the Eq. (24), by including the heat effects we obtain (74) In the absence of dissipation Eq. (74) reduces to we acquired the GLTB spacetime by using Eqs. (73) and (8) in the following form Structure Scalars Y T F , X T F , Y T and X T for the line element (76) are given by To continue our analysis and obtain some specific solutions, we need to impose some constraints. The criteria for choosing such constraints will be based on how similar the obtained solutions to LTB geometry so we shall regard the extensions of LTB geometry based on symmetric properties and structure scalars to analyze the impact of charged fluid and modified gravity. We shall study the transport equation to obtain the temperature profile for the diffusion case.

Transport equation
To calculate the temperature of the self-gravitating system under the diffusion approximation (q e f f = q e f f , e f f = 0), we will apply a transport equation for dissipative fluids derived from the Muller-Israel-Steward second-order phenomenological theory. The heat flux transport equation can be written as where the heat flow, thermal conductivity, temperature, and relaxation time are represented by q e f f , K , T, and ζ , respectively. Set ζ = 0 to restore the Eckart-Landau equation. In the case of a spherically symmetric distribution, the transport equation has only one independent element which is given by in the shortened edition of the theory, the last term in Eq. (81) is absent. As a result, Eq. (82) implies if e f f = 0 then we obtain from Eq. (52) where where g(r ) is an arbitrary function, and the temperature of the matter distribution can be determined in terms of metric coefficients using Eqs. (82) and (83).

Extension through scalar functions
We have drawn the following conclusions from our analysis of structure scalars. The structure scalars Y T F and Y T entirely determine the evolution of expansion and shear scalar in the case of dissipative dust fluid and the difference between their terms is the same. In the case of LTB spacetimes, the term κ is absent from the expressions of Y T and Y T F . Next we examine the structure scalars Y T and Y T F by keeping in view the above conclusions and maximal similarity condition occur between LTB and its generalized spacetime(GLTB). Therefore by comparing Eqs. (69) and (78), we havė which exactly match with Eq. (103) of [1] but in expression (86), the value of K (t, r ) contains the extra f (R) terms and its value is mentioned in (74). integrating Eq. (86) produces where C 1 (r ) is a constant of integration. Using Eqs. (74) and (75), we have integration of Eq. (87) yields where C 2 (r ) is another arbitrary constant of integration. Taking t Derivative of Eq. (88), we obtaiṅ comparing the Eqs. (74) and (90), we have The relation between integration constants C(r ) and C 2 (r ) can be found by comparing the Eqs. (88) and (91) Based on the suppositions put on Y T and Y T F , we can obtain a GLTB. Taking the temporal and radial dependence of R from (91), the dissipative flux of GLTB may be derived up to two functions of r . The regularity condition C 1 (0) = 0 must be satisfied by the function C 1 (r ), and its disappearance returns to the original LTB solution. The rest of the physical variables are presented by the field equations. To explain the procedure described earlier, we consider an example from the parabolic subclass of LTB which corresponds with Eq.
(110) of [1] because we are following the strategy mention in [1]. Therefore, we take: The parameters for parabolic LTB solution are f (r ) and T (r ). Using Eq. (93), we manifest the integral dt (C ) 2 as plugging τ 3 = −t + T (r ) in Eq. (94) and integration yields: where b = 2 3 T (r ) f (r ) and a = f (r ). Putting Eq. (95) in (91) and then utilizing (93) we obtain This expression exactly corresponds with Eq. (113) of [1] but the termq e f f contains the extra curvature terms and its value is mentioned in Eq. (27). When we desired LTB perturbations then the amplitude of dissipation is being regulated by the function C 1 (r ), which can be adjusted as little as desired.
The results have shown that these GLTB can be transformed into LTB when we rule out dissipative fluxes. For a smaller value of C 1 (r ) associated with LTB geometry, the physical features of GLTB have logical consequences.

Symmetric properties
In this subsection, we present another way to obtain a GLTB spacetime under the impact of Palatini f (R) gravity and electromagnetic field. This study is based on the assumption that the symmetric properties of GLTB spacetime correspond to LTB. The Lie derivative of stress tensor w.r.t. ξ is defined as Taking the Lie derivative of the effective energy-momentum tensor in the presence of charge using (97) we get the following equations: 1

Diffusion approximation
Following the same strategy as deployed by Herrera et al. [1], we use Eqs. (24) and (73) using Eqs. (73), (84) and (101), we get the following The Eq. (104) can be manifested in the following form multiplying Eq. (105) by q e f f B, we get now, from Eq. (106) we may write as using Eqs. (101) and (107), we get: Now using the Eqs. (84), (107) and (108), we have which exactly corresponds with Eq. (102). Integration of Eq. (102) yields where C (r ) denotes the integration constant. For a specific situation, using the value of C(t, r ) by using Eq. (93), we obtain inserting Eq. (111) in (84), we get the following relation The GLTB spacetime derived using diffusion approximation is defined by four functions of r , although one function vanishes due to radial coordinate invariance. Two functions may be discovered using the LTB seed model, and the remaining one can be determined using the Eq. (106) under the assumption that the LTB solution has the same radial dependence as ξ . Moreover, Eq. (112) is used to calculate the temperature profile using the transport equation for diffusion case. These results are compatible with those of [1].

Streaming out approximation
The Bianchi identities take the following form in the case of streaming out approximations (q e f f = 0, e f f = 0) as obtained by Herrea et al. in the case of GR [1].
(98)-(100) turn into the following form by keeping in view the assumption ( e f f = 0, q e f f = 0) [1], Solving Eqs. (117), (118) and (119), we have To indulge the property of maximal similarity between LTB geometry and GLTB, we assume that ξ 0 = F(t), and in this scenario, the Eqs. (117)-(119) take the form [1] F(t)(μ e f f + e f f ) . + ξ 1 (μ e f f + e f f ) The Eq. (123) may be written as or Eq. (123) may also be rewritten as under these assumptions the Eq. (46) becomes Or Eq. (126) may be rewritten as The following approach can be used to obtain GLTB based on LTB (in the streaming out approximation) in the presence of f (R) gravity and electric charge: Assume that the radial dependence of the LTB solution is comparable to the vector field ξ 1 . Take the value of ( e f f B 2 ) from Eq. (127) and insert it into Eq. (125).
Using a definite LTB model with a specified vector field ξ , assuming F(t) = 1 and integrating Eq. (125), we obtain the effective radiation density with respect to r . These results are compatible with those of GR under usual limit [1].

Discussion
The observations of type Ia supernovae have motivated the researchers to study the inhomogeneous LTB spacetime for determining the accelerated rate of the universe expansion. In this article, we have studied the inhomogeneous LTB geometry and its properties in the framework of f (R) gravity. The f (R) gravity enables us to a better understanding of the inflationary era and late-time accelerated rate of the universe expansion. We have considered the dissipative dust fluid along with spherically symmetric spacetime to study the consequences of charge and f (R) terms. To continue further, we have evaluated the Maxwell-field equations for compact objects and have seen that the field equations contain the higher-order terms due to the Palatini formulation. We have computed the matter variables in the context of f (R) gravity which has a substantial role in the field of relativistic astrophysics and found that the four acceleration components appear only due to the Palatini version. To observe the amount of matter inside the geometrical object, we utilized the interpretation for mass function presented by Misner and Sharp and devised the Darmois junction conditions for smooth fitting of the solutions. We have established an expression for Weyl scalar in terms of fluid variables along with dark source terms that demonstrated the consequences of the charge. The method used by Ellis has been adopted to investigate the evolution equations to examine the inhomogeneities that arise due to the presence of the Weyl tensor.
Herrera et al. [26] computed a set of scalar functions using the orthogonal decomposition method to express the physical behavior of the matter variables and these scalars are known as structure scalars. We set up the expressions for tensors Y αζ and X αζ in f (R) gravity and further break down into their trace and trace-free parts. The energy density of the fluid is addressed by the scalar X T with extra curvature terms that identify the outcomes of an unknown form of energy and matter. The scalar Y T carried the information about the pressure anisotropy with higher-order terms due to modified gravity that signifies the existence of DM and DE. These scalars are also served to express the solutions of field equations using a static spacetime and we have associated them with metric coefficients using field equations. The previously attained differential and evolution equations are expressed in the form of scalar functions through modified field equations.
We generalized LTB spacetime by following the strategy as mentioned in [1] in the background of Palatini f (R) gravity and electric charge. We studied the properties of LTB geometry and analyzed the behavior of evolution equation by considering the masses of compact objects through graphical representation and found the structure scalars for LTB geometry using (63). The LTB spacetime has been extended for dissipative fluid which is known as generalized LTB spacetime (GLTB) and we have calculated the scalar functions for GLTB with the inclusion of charge and curvature terms. We compared the structure scalars obtained from LTB and GLTB to achieve the maximum similarity between two spacetimes and studied the symmetric properties of LTB and GLTB geometry using two techniques. we have estimated a transport equation to obtain the temperature profile for the diffusion case. In this article, the generalization of LTB spacetime has been modified and our calculation has generalized the previously obtained results under the influence of charge in f (R) gravity. Our results are compatible with those of GR [1] under the usual limit.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A
Appendix B