Dynamics of evolving cavity in cluster of stars

We examine the evolution of cavities within spherically symmetric cluster of stars in high curvature gravity. For this purpose, we use f(R) gravity through the Starobinsky model to incorporate dark matter effects in the discussion. In particular, we check the physical significance of the f(R) model by associating it with the observational data of stellar object 4U182030. For evolution of the cavity, we consider the purely areal evolutionary phase by assuming that the proper distance (in a radial direction) among neighboring stars remains constant. The analytical solutions are obtained among which a few solutions fulfill the Darmois conditions. It is found that the evolution of the cavity in a cluster of stars is highly controlled by the influence of dark matter.


Introduction
In modern cosmology, dark matter and dark energy have been the most fascinating objectives. Planck's collaborators provided observational data that indicates 5% seen matter, 71% dark energy (DE) and 27% dark matter (DM) [1]. Dark energy is a mysterious form of energy whose presence helps to compress the curvature of space which in turn speeds up the cosmic expansion. Baryonic matter refers to the visible matter in the universe while DM is an invisible hypothetical type of matter which does not interact with electromagnetic radiations. It can only be observed via its gravitational effects on the observable matter.
The theory of f (R) gravity, being an alternative candidate of general relativity (GR), has received a lot of attention due to its simplicity in describing the interaction of DE and DM [2][3][4][5][6][7][8]. This theory is derived by modifying the Ricci scalar in the Einstein Hilbert action to its generic function f (R) that a e-mails: rubab.manzoor@umt.edu.pk; dr.rubab.second@gmail.com b e-mail: saadia.icet@pu.edu.pk (corresponding author) c e-mail: daoudchoudhary1673@gmail.com represents the higher order curvature principle for a gravitational framework. Many researchers used different models of f (R) functions to describe the impact of DM and DE [9][10][11][12][13]. In this context, Starobinsky [14,15] used the second order curvature term, f (R) = R + R 2 , to formulate Einstein's equations with a quantum one loop system. He used this model to describe the exponential cosmic expansion and current time cosmic acceleration for an early time as well as the present time of power-law inflation. This model can also deal with DM issues [16].
Until about a decade ago, star clusters in general, and globular clusters in specific were assumed to be basic formations composed of gravitationally bound stars with a common origin [17]. Many observational surveys like the rotational loops (out to tens of kpc), gravitational lensing (out to 200 kpc), and heat in clusters showed that 95% of a cluster of stars is made up of any unknown segment of DM [18][19][20]. The challenges of galactic rotational curves, stellar clusters and mass differences lead to the presence of DM in cluster evolution [21][22][23]. Thus it exhibits influence of DM in the evolutionary mechanism. This may enable us to unveil some hidden facts that are not yet detectable in the universe's structure. In this context, many useful theoretical works have been done to explore the evolution and development of stellar structures [24][25][26][27][28][29].
At the large scale, about 80% of the universe volume is composed of cosmic voids. These voids are large spaces between clusters of stars, which contain very few or no stellar structures. Cosmic voids are not empty, but they are underdense enough that the phenomenon of gas-stripping galaxies is extremely rare. In this way, these voids provide a unique environment for the description of the evolution and formation of galaxies [30]. It is believed that the actual universe has a sponge-like structure governed by the voids [31]. Many observations indicate that voids of a characteristic scale of 30 h −1 Mpc exist in the present volume of the universe, where h shows the dimensionless Hubble parameter, H 0 = 100 h km s −1 Mpc −1 [32]. However, voids of different sizes, rang-Eur. Phys. J. C (2022) 82:739 ing from minivoids [33] to supervoids [34] can be discovered. According to cold DM ( CDM) simulations, voids should have a density of low-mass halos of about 1/10 the cosmic mean. Many researchers also noticed that the galaxies have a poor tendency to get higher rates of star formation towards void centers. It is important to notice that voids are neither spherical nor empty in general. In simulation results or deep redshift observations, voids are generally referred to spherical vacuum cavities surrounded by a fluid. The issue of cavity evolution followed by a central explosion of a spherically symmetric distribution was firstly discussed by Skripkin [35]. Accordingly, a cavity appears around the center of fluid distribution under a constraint called the Skripkin condition. This topic has been studied in detail, demonstrating that the scalar expansion disappears under the Skripkin conditions (for isotropic fluid with the constant distribution of energy density). A further detailed study on this problem declared that expansion less fluid distribution requires the existence of a cavity within the system [36]. Herrara et al. [37] studied the evolution of expansionfree isotropic fluid and found that the Skripkin model is not consistent with the Darmois junction conditions [38]. Also, expansion-free homogeneous dust fluid fails to attain physical interest because of negative energy density. Herrera et al. [39] also discussed the evolution of the cavity by using the constant proper radial distance between fluid neighboring particles. Sharif and Bhatti [41] described shear-free cavity models with plane symmetry. Yousaf and Bhatti [42] discussed the cavity evolution and instability of relativistic systems. Sharif and Mumtaz [43][44][45] studied the dynamics of collapsing gaseous masses with and without dissipative effects.
In this paper, we discuss evolution of the cavity for a spherically symmetric cluster of stars in the context of f (R) gravity. For this purpose, we consider purely areal evolution in which variation of the proper radial distance between any two infinitesimally small particles of fluid per unit of proper time vanishes [39]. We are particularly interested in analytical models that, despite being relatively simple to assess, incorporate some of the key characteristics of an actual situation. It should be noted that we are only interested in the evolution of cavity, not in the dynamics and conditions of its formation. In this context, the fluid is constrained by two hypersurfaces. The external hypersurface separates the fluid distribution from a Schwarzchild or Vaidya spacetime (depending on whether we assume the evolution to be adiabatic or dissipative) and the internal one is taken as Minkowski spacetime that delimits the cavity.
In general, voids are neither spherical nor empty but for simplicity, we can model vacuum spherical cavities within the cluster of stars as a predecessor of voids. The paper is formatted as follows. Section 2 describes general formalism of f (R) theory. The next section deals with a spherically symmetric cluster of stars in the f (R) model along with junction conditions. In Sect. 4, we apply the Strobinsky model whereas Sect. 5 deals with the purely areal evolution as compared to the radial velocity and describes different models of the cavity in a cluster of stars. Section 6 summarizes our results.

The f (R) theory
The f (R) theory is followed by the action [2][3][4][5][6][7][8] Here κ shows the coupling constant, S M represents the matter related segment and f (R) describes a non-linear Ricci scalar part. The metric f (R) field equations are obtained by varying (1) with respect to g ιω as follows where R ιω represents Ricci tensor, f R = d f d R , and T ιω shows ideal energy momentum tensor for matter distribution given by Equation (2), in terms of Einstein tensor, can be composed as with where the superscript D represents dark source parts (may be DE or DM).

Model of cluster of stars
To discuss cavity formation in the cluster of stars, we consider a cluster that is composed of compact stars and exotic material like DM. These types of clusters are found in some dwarf irregular galaxies [40]. Since clustering of stars involves a strong gravitational field thus we model the cluster of stars as a system of self-gravitating fluid having DM and stars as fluid particles. Specifically, we assume a spherically symmetric distribution of star cluster, circumscribed by a spherical surface. The system is considered to be locally anisotropic, experiencing dissipation as a heat flow in diffusion approximation. The general interior metric (in comoving coordi-Page 3 of 15 739 nates) is given as where A, B and C are positive functions of t and r . We take the coordinates x 0 = t, x 1 = r, x 2 = θ, x 3 = φ. The coordinates C and r have the same dimensions whereas A and B are dimensionless. According to Eq. (6), the term Bdr denotes proper radius inside any spherical hypersurface (e) and C is the areal radius. The energy-momentum distribution inside (e) is given by Here μ e f f , p e f f r , q e f f ι , V ι and χ ι represent energy density, radial pressure, dissipation term describing heat flux, four velocity and radial unit four-vector, respectively. The four vectors satisfy the following identities Moreover P e f f r where superscript D represent DM contributions and the terms without superscript indicate baryonic matter terms. The f (R) field equations for the line element (6) are given by The prime and dot denote differentiation with respect to r and t, respectively. The kinematics of an evolving cluster can be described with the help of four acceleration (a α ), expansion scalar (θ ) and shear tensor (σ ιω ) given by where represents the projection tensor. For the metric (6), we have From Eqs. (17) and (20), the only non-zero four acceleration component and its scalar are written as also the expansion scalar is The non-zero components of shear tensor yield and its scalar is given by where The shear tensor can also be written in terms of projection tensor as For the sake of convenience [46,47], the energy momentum tensor (7) can be rewritten in the following form withP Misner and Sharp [48] introduced a mass function m(t, r ) which represents energy within a system [49] given by In order to discuss the dynamical characteristics of the system, we use the proper time and radial derivatives given by [48].
Using Eq. (32), we may define the velocity U of a collapsing or expanding fluid distribution as the differential of areal radius w.r.t the proper time given as Equation (31), through the above equation, can be expressed as Moreover, Eq. (31) along with the field equations (32) and (34) gives [36] whose integration yields Here we consider a regular center for the system and hence m(0) = 0. Using Eqs. (16), (22) and (25), we get which can be rewritten with the help of Eq. (33) as The Weyl tensor describes the effects of tidal forces due to gravitational effects within the system. It is denoted by C αμβν and evaluated in terms of Riemann tensor R ρ αβσ , Ricci tensor R ρ β and Ricci scalar R as The Weyl tensor can further be converted into the electric E αβ and magnetic parts H αβ . For spherically symmetric distribution, the magnetic part vanishes whereas the electric part is defined as or it can be written in terms of four vectors as Here This equation, through the field equations and Eq. (31), can be written in terms of mass function as [36]

Junction conditions
We assume Vaidya spacetime as an exterior spacetime (i.e., we consider massless outgoing radiation case) defined by where M(υ) and υ denote the total mass and retarded time, respectively. In the absence of thin shell geometry, Darmois junction conditions are used for the the matching of Vaidya spacetime to the non-adiabatic sphere, along the surface r = constant = r (e) . These junction conditions include consistency of the 1st and 2nd fundamental forms, based on matching hypersurfaces [38]. The results of Darmois conditions yield and Here L represents the total luminosity of cluster given by with as the entire luminosity determined by an observer in the state of rest at infinity. Equation (47) along with the field equations provides In the case of cavity formation, we also apply matching conditions between the solution and the Minkowsky spacetime to delimit this cavity on the boundary surface. If (i) denotes the boundary surface between the cluster fluid and the cavity, then matching conditions associated with Minkowski spacetime within the cavity and the cluster distribution give q e f f (e) = P e f f r .
Let us consider an empty cavity, then we have L

The Starobinsky model
The f (R) theory of gravity is consistent with DM and DE at a large scale (stellar object scales like clusters). Starobinsky generalized Einstein Hilbert's action with higher order curvature term to construct field equations related to quantum one loop distribution. This derived model can describe cosmic exponential expansion as well as current time cosmic acceleration for both early and current times of power-law inflation. It is defined as with ε as a positive real number. It can solve various issues like inflation, cosmic acceleration and DE [14,15]. Moreover, it has been shown that the respective model can be used for DM problems according to WMAP data for ε = 1 6M 2 with M = 2.7 × 10 −12 GeV, (ε = 2.3 × 10 22 Ge/V 2 ). The square-order curvature term behaves like an extra degree of freedom (scalar graviton) whose interaction with standard model particles yields an abundance of the thermal energy which provides DM regions. The scalar graviton of mass M < 10 −12 generates Yukawa force of attraction between DM particles having dissimilar masses [16]. For f (R) → R, GR solutions are recovered. The field equations (13)-(16) for this model become

Physical significance of the f (R) model
Here we examine physical importance of the function f (R) = R + R 2 for stellar bodies. For this purpose, we analyze the viability of parameter associated with a compact object. Thus we consider a simple model that relies on spatial astral density (a theory same as the de Vaucoulour's statement in the exterior zone that avoids a fixed center). Hernquist and Jaffe [50,51] firstly proposed such kinds of models with central astral densities proportional to r −1 and r −2 . These models lead to a variety of energy density distributions with various central slopes described by Hereα represents a scaling radius and M shows the total mass which is in proportion to r τ at the center. The value of τ is fixed as [0, 3) whereas τ = 1 and τ = 2 restricts the model to Hernquist and Jaffe models, respectively. The Krori-Barua anstaz method is used to investigate the metric functions related to interior spacetime [52]. The interior spacetime is given by whereã =Br 2 +C and b =Ãr 2 . To represent a realistic model for compact stellar object, we employ the Schwarzschild spacetime as exterior geometry having an asymptotically flat and static outside area. The Schwarzschild Eur. Phys. J. C (2022) 82:739 metric is given by We are now applying boundary conditions to both spacetimes associated with exterior and interior geometries. At the condition r = R, the continuity of exterior and interior line elements produces where the superscripts (−) and (+) represent the interior and exterior surfaces of the stellar object, respectively. Thus, Eq.
(63) implies Here we examine the physical significance of the function f (R) = R + R 2 for stellar systems. For this, we analyze the viability of parameter associated with a compact object. So, we consider a simple model that relies on spatial astral density (a theory same as the de Vaucoulour's statement in the exterior zone that avoids a fixed center). Hernquist and Jaffe [50,51] firstly proposed such kinds of models with central astral densities proportional to r −1 and r −2 . These models lead to a variety of energy density distributions with various central slopes described bỹ Substituting the values ofÃ,B andC in Eqs. (56)-(59), we can explore the behaviors of energy density, radial, and tangential pressures associated to the matter distribution.

The purely areal evolution condition and radial velocity
Here we discuss a new definition of collapsing velocity U , previously, it is defined as the change of areal radius (C) per unit proper time. The velocity can also be defined as a variation of the infinitesimal proper radial distance between two adjacent points (δl) per unit of proper time, i.e., D T (δl).
It has been shown that this infinitesimal rate of change is related to the shear and expansion effects, given by [36] which, through Eqs. (22) and (25), becomes Using Eqs. (22), (25), (34) and (70), we obtain and This shows that the areal velocity U , being the variation of C (areal radius) of a layer of fluid particle, is generally different from D T (δl) which is the relative velocity of neighboring layers of fluid particles present within the clusters. Equation whose integration with respect to r gives with ζ(t) being an integration function. It can be noticed that in the situation when all the spherical surfaces of a cluster of stars including the center (r = 0) is filled with baryonic as well as non-baryonic fluids, we should imply the regularity condition (ζ = 0). As we are interested in the formation of the cavity around the center, such condition is not needed. In other case, Eq. (74) with (71) leads to Hence, for non-dissipative case (q e f f = 0 ⇒ q M = q D = 0), the purely areal evolution condition suggests U = U (t). This result is certainly unsuitable with a regular symmetry center except (U = 0). To have a purely areal evolution condition to be compatible with a time dependent stage (U = 0), this yields • The cluster has no symmetry center. Or • The center is encircled by a compact spherical region of other spacetime, appropriately matched to the rest of the fluid inside the cluster.
Here we are going to reject the first case because this indicates the unusual geometry of a star cluster without a center. Furthermore, for the second case, we select Minkowski's inner vacuum spherical vacuole. It can be noticed that dissipation due to DM (q D ) might be some sort of squandering of non-baryonic particles. This type of dissipation is due to the gravitational energy of DM. Now a question arises here, can we consider a non-dissipation case in the presence of DM?
The evidence of DM is based upon its gravitational effects so the dissipation due to DM gravitational effects can not be neglected easily. Hence during the evolution of the cluster of stars having DM, the non-dissipative case is very special or unusual.
Let us now assume an alternative case that entirely belongs to areal dissipative evolution (D T (δl) = 0). If the gravitating fluid (for baryonic as well as non-baryonic parts) fills the entire spherical cluster, we get a symmetry center and there is no cavity surrounding the center. In this situation, we need to put ζ = 0 such that Eq. (75) becomes which is congruent with a regular symmetric center. For this situation, we shall consider a cavity enclosing the center in an ad hoc manner. The following subjective argument can propose this assumption. Consider an outwardly dissipative condition (q e f f > 0) which implies that all the terms inside the integral are positive such that Eqs. (72) and (76) give > 0, U > 0. Inversely, for q e f f < 0, we get < 0, U < 0. Since we are dealing with dissipation due to matter and DM so one of the sources for this dissipation is the gravitational energy. Moreover, according to Kelvin-Helmholtz evolution phase [55], when all the outgoing dissipative flux originates from the gravitational energy, we should anticipate contraction instead of expansion. In contrast, inwardly directed flux of dissipation leads to expansion. The above discussion shows that dissipation due to heat flux behaves differently from dissipation due to gravitational energy. So the outgoing dissipation due to DM (q D > 0) contributes to the contraction phase of the evolution and inversely, q D < 0 helps in expansion of the evolving cluster. Hence, we observe that the condition D T (δl) = 0 (purely areal evolution phase) shows up to be particularly appropriate for describing evolution of a cluster of stars with a cavity enclosing the center [39]. Finally, using Eqs. (71) and (72) in (26), the purely areal evolution phase can be described as Next, we shall explore some models for cluster of stars having purely areal evolution.

Models of cavities in clusters
In this section, we explore the general characteristics of models of clusters that satisfy the purely areal evolution case.
Here we use the general concept similar to that suggested by Skripkin [35]. According to his proposal, an explosion at the center causes an overall expansion throughout the fluid, resulting in the development of a cavity around the center. But, here is the distinction that we are considering purely areal evolution D T (δl) = 0 rather = 0. In this context, Eq. (70) impliesḂ = 0 (Ċ = 0) which gives B = B(r ) and we can take B = 1 without the loss of generality. As previously indicated, the physical importance of this sort of model comes from the fact that the conditionḂ = 0 requires the existence of a center-based cavity. For this case, the field equations become Also, the non-zero parts of the Bianchi identities T (e f f )−αβ ;β = 0 are given by In specific geodesic condition, we have A = 0 which implies A = 1. At this point, the field equations (79)-(81) become It follows from Eq. (31) that Equation (44) implies This along with Eq. (51) implies that the respective model satisfies Dramois junction conditions only for its absorbing dissipative energy behavior, i.e., q e f f < 0 otherwise, we get M = 0, m < 0. Furthermore, the condition q e f f < 0 implies two possibilities either q < −q D or q D < −q.

Model consistent with Darmois conditions
Now we consider some simplified analytical models that are not based upon thin shells either (e) or (i) , but are consistent with Darmois conditions. In this context, the nondissipative models are found to be the simplest ones.
Taking q e f f = 0 along with Eq. (79), integration gives where h 1 (t) is an integration function of t. Without loss of generality, after re-parameterizing t, we can assume with It can be noticed from Eqs. (34) and (91) that the velocity U becomes 8π P e f f r = 1 8π P From Eqs. (96) and (97), it can be observed that Computing function of mass (31) with (91) and B = 1, we have We shall now explore some specific cases.

Conformally flat case
If the spacetime among r = r (e) and r = r (i) is conformally flat, i.e., = 0, then Eq. (101) gives where W 1 is an arbitrary function of r . After integration, we get Here h 2 (t) is an integration function of t. A comparison of Eqs. (99) and (103) gives such that h 2 can be obtained from the junction condition (52) given by Consequently, this sort of models are described through a single function W 1 (r ) whose selection depends upon the satisfaction of the remaining Darmois conditions. Furthermore, the evolution is shear-free for isotropic non-dissipative spherically symmetric conformally flat spacetimes. But this is not applicable in anisotropic case. Thus, the models examined here are necessarily considered anisotropic.

Tangential pressureless case
Considering P e f f ⊥ = 0, integration of Eq. (98) yields where W 2 (r ) is a function of r . Using Eq. (106) into (96), we have Also, Eq. (99) implies In order to get viable models, we need to consider a particular type of energy density or mass function. For instance, let us consider In this way, Eqs. (37) and (106) imply Using Eqs. (100), (106) and (110), we get Substituting Eq. (110) into (99), we have where whose integration gives Eur. Phys. J. C (2022) 82:739 Here r 0 (t) is an integration function of t. Assessing Eq. (116) on (i) , we get This shows the 1st order nonlinear equation for C (i) that can be determined for every r 0 (t) function. The result of this integration, combined with Eq. (116), provides the absolute information required to get r and t dependency for all metric and physical variables. It can be found that energy density within the fluid distribution remains positive and regular. Additionally, the condition r 0 (t) with 0 <Ċ (i) C −1 ≤ 1 guarantees the presence of positive pressure which is smaller than the energy density.

Physical conditions for stellar object
The densities and gravitational fields of the whole star cluster are different than its component star. Sometimes it becomes complex to explore the entire density and gravitational effects of the cluster of stars in strong field regimes. So for the sake of convenience and to check the physical significance of the obtained cavity model, we choose a binary star 4U 182030 as a test star from the component of compact stars cluster and check its physical features under the influences of DM. In this way, we may suggest the mutual relation between a compact star and DM in cluster of stars or galaxies. The physical viability of a compact star depends on the following conditions to be satisfied throughout the configuration: • The metric potentials A(r ), B(r ) and the matter components μ, P r , P ⊥ must be well defined at the center as well as regular within star. • The positivity of energy density, i.e., μ ≥ 0 is required throughout the stellar object interior. Being positive at the center, it should be monotonically decreasing to the boundary inside the star, mathematically dμ dr ≤ 0. • The radial and the tangential pressures should be positive, i.e., P r ≥ 0, P ⊥ ≥ 0, whereas their gradient must be negative inside the stellar object configuration, i.e., d P r dr ≤ 0 and d P ⊥ dr ≤ 0. At the boundary, the radial pressure must vanish while the tangential pressure may not be zero.
• The causality condition must be satisfied for physical viability of the model, i.e., the velocity of sound should be smaller than 1 in the stellar object interior. Mathematically, 0 ≤ V r ≤ 1 and 0 ≤ V ⊥ ≤ 1, where V r = d P r dμ and V ⊥ = d P ⊥ dμ represent the radial and transverse velocity of sound, respectively.
• The smooth matching of the interior and exterior metric functions at the boundary is also required. • The adiabatic index must be greater than 4 3 for the stability of stellar object configuration.
6.1 Physical analysis with observational data of 4U 1820 − 30 Here, we discuss physical interpretation of the solutions obtained for cavity model. For physical compatibility of stellar model, we combine our calculations with the data obtained for the pulsar 4U 1820 − 30 [53,54].

Metric potentials
Firstly, we examine the physical viability of the model by using Eqs. (61), (66)-(68) and explore the physical behavior of gravitational potentials A(r ) and B(r ) graphically as shown in Fig. 1. It is observed that this model has finite values of A(r ) = eB r 2 +C , B(r ) = eÃ r 2 at the center (r = 0) of the star. One can also find that A (r = 0) = 0, B (r = 0) = 0 which depict the regularity of metric potentials at the center and also show that the metric is well-behaved throughout the interior of star. This implies fulfillment of the requirement for realistic star.

Energy density, pressures and their gradient
For the sake of physical analysis in f (R) gravity, let us begin with the main components of matter distribution, i.e., energy density and pressures (radial and tangential). To incorporate DM effects, we take the Starobinsky model f (R) = R + R 2 . We have plotted graphical results in 3D with effects of DM for anisotropic sphere in f (R) gravity by varying radial component r as well as parameter (DM parameter). In this way, we find their role in the emergence of various physical factors. To provide a comparison between the stellar object model in f (R) and GR, we use data of binary pulsar 4U 1820 − 30. Moreover, the 2D illustration provides the corresponding results in GR for = 0.
First of all, we analyze the behavior of energy density graphically corresponding to DM ( = 0) as well as baryonic matter (in GR limits = 0). The behavior of matter density μ with DM effects as a function of r and is shown in Page 11 of 15 739  Fig. 2 (upper panel). It can be seen that the energy density remains positive throughout the star's interior and decreases monotonically towards the boundary inside the stellar object interior but does not vanish at the boundary as increases.
The pressures under the influence of DM in the radial and tangential directions are also shown in Fig. 2 (middle and lower panels, respectively). It is noticed that the radial and the tangential pressures remain positive, i.e., P r ≥ 0, P ⊥ ≥ 0 throughout the cluster. For increasing values of , the radial pressure associated with DM rise to the maximum value at the center but vanishes at the boundary for all choices of (middle panel). Furthermore, the behavior of matter tangential pressure is also shown in Fig. 2 (lower panel). The tangential pressure decreases throughout the star but at the center, its value reaches a maximum as the value of increases. Thus, the pressures for both cases remain decreasing in the radial direction while the radial pressure disappears at the boundary as it should be.
We also present our results with = 0 for energy density and pressures in GR limits (Fig. 3). The resultant illustration in 2D represents the positivity of all matter components throughout the cluster. All the three matter components remain positive throughout the interior and tend to decrease towards the boundary, whereas the tangential pressure vanishes at the boundary. Also, they attain maximum value at the center but not more than the values with DM effects = 0. It is noticed from the above analysis that for the increasing values of , the energy density and pressure components of DM become dominant in contrast to the matter-energy density and pressures. Thus, the results in presence of DM are more consistent with observations on a cluster of stars like  galaxy rotational curves, according to which, at the center of the cluster the density and pressure become enormous [56]. For physical viability of our cavity model, the energy density and pressure gradients should also be negative. Here we show the results for gradients in tabular form with DM effects. Table 1 shows that dμ dr ≤ 0, d P r dr ≤ 0 and d P ⊥ dr ≤ 0 within the cavity by taking = 100 and different choices of radii.
The result for energy density and pressure gradients with = 0 is given in Fig. 4 which shows consistency of our results with GR, i.e., negativity of gradients within the cluster.

Energy conditions
We also analyze the viability of the stellar object model by energy conditions to be satisfied by the matter components of DM. It is mentioned here that for WEC, the energy density is always positive, i.e., μ ≥ 0 throughout the interior of a cluster as already discussed, thus we need to find the behavior of remaining energy conditions for DM. The values of energy  Table 2.
For DM, it is observed that WEC, NEC and SEC are satisfied in the presence of DM. Moreover, DEC is violated at the center, whereas μ − P r ≥ 0 and μ − P ⊥ ≥ 0 for all other choices of r in the interior of the star. We have plotted the behavior of energy conditions for the matter case ( = 0) graphically as shown in Fig. 5. We find that all the energy conditions are satisfied throughout the star's interior in GR limits.

Velocity of sound
For the physical acceptability of the respective model for star, we examine causality conditions to be satisfied within the interior. The radial and transverse velocities are defined by We provide the behavior of radial and transverse velocities versus r and = 100 in Table 3. It can be observed, from the table, that in the presence of DM we have 0 ≤ V r ≤ 1 and 0 ≤ V ⊥ ≤ 1 at the boundary. This shows that the DM affects the propagation of the speed of sound. For matter case, we also perform a graphical analysis and observe the fulfillment of causality condition, i.e., 0 ≤ V r ≤ 1 and 0 ≤ V ⊥ ≤ 1 (Fig. 6). It is mentioned  here that the requirement of smooth matching of the interior and exterior metric functions at the boundary is already performed in Sect. 5.2.

Adiabatic index
The adiabatic index is related to the stability of a relativistic anisotropic stellar object configuration. Mathematically, it is defined as For physical viability of stellar object model, must be greater than 4 3 leading to stable configurations. We analyze  the behavior of with DM effects as given in Table 4 which appears consistent with the required condition. We plot the behavior of adiabatic index for ordinary matter as shown in Fig. 7 which also depicts > 4 3 .

Conclusion
The study of evolving cluster of stars has been of great interest for the description of phenomena like DM, voids and galactic structures. It is believed that the cavity model described by purely evolution conditions might be used as a void predecessor. Voids have made a major contribution to the scientific understanding of the cosmos. In this paper, we have investigated the evolution of the cavity within a spherically symmetric compact stars clusters. Generally, voids are neither spherical nor empty (according to redshift and simultaneously observations), however for sake of simplicity, they are shown as a vacuum spherical cavity surrounded by a fluid. We have resembled the cluster by a self-gravitating fluid where the collection of fluid particles represents a collection of stars (baryonic matter) and DM particles (nonbaryonic matter). To accommodate DM's contribution in the discussion, we have adopted the idea of high curvature gravity ( f (R) gravity). We have investigated the impact of nonbaryonic matter on the evolution of cavities within the cluster by using Starobinsky model, f (R) = R + R 2 . Specifically, we have analyzed the significance of the Starobinsky model for the evolving cluster. For this purpose, we have explored the behaviors of matter and DM by using data associated with star 4U 1820 − 30. It is found that the role of DM is dominant in contrast to the baryonic matter.
We have applied purely areal evolution condition to discuss the presence of cavity around the center. All dynamic equations under this condition have been derived and various models have been discussed. The obtained models are associated with various features like expansion-free, geodesic and conformally flat conditions. These models satisfy Darmois junction conditions for both hypersurfecs (i) and (e) in the presence of DM and normal matter. We have found that the evolution of cavity in cluster of stars is deeply controlled by the effects of DM. Finally, we perform the physical analysis of the obtained cavity models by relating them with the data of star 4U 1820 − 30 and investigated the physical viability of the star configuration. In this context, we have analyzed behavior of various physical factors within the cavity and compared our results with that in GR. It is found that the contribution of density and pressure becomes maximum at the center which gradually tends to decrease towards the boundary whereas their gradients remain negative throughout the stellar object interior. It is worth mentioning that our results with DM are well consistent with GR except the case for tangential pressure. For DM, the DEC is violated at the center of stellar object cluster as compared to the conditions in GR. Moreover, our stellar object cluster model is stable with > 4 3 in the presence of DM. We conclude that the inclusion of DM significantly controls the evolutionary mechanism of cavity in cluster of stars.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There is no observational data related to this article.] 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 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.