Cavity Evolution and Instability Constraints of Relativistic Interiors

In this manuscript, we have identified the dynamical instability constraints of a self-gravitating cylindrical object within the framework of $f(R,T)$ theory of gravity. We have explored the modified field equations and corresponding dynamical equations for the systematic constructions of our analysis. We have imposed the linear perturbations on metric and material variables with some known static profile up to first order in the perturbation parameter. The role of expansion scalar is also examined in this scenario. The instability regimes have been discussed in the background of Newtonian and post-Newtonian limits. We found that the dark source terms due to the influence of modification in the gravity model is responsible for the instability of the system.


Introduction
The global properties of our universe have been studied extensively by the researchers within the background of general theory of relativity (GR). Many new insights of astrophysics and cosmology have revealed the unexpected picture of the universe. The latest data of some reliable sources such as supernovae surveys and cosmic microwave background radiation have put forth the dark side of the universe. Therefore, one has to accept that current matter-energy contents and the evolutionary picture of the universe is astonishing and requires some explanation. During the last couple of decades, the discussions in the standard GR due to the observational evidence is unable to describe the key features of the present low energy universe without imposing certain assumptions. Particularly, it requires the inclusion of some exotic contribution in the matter contents of the cosmos in order to study the dynamics of stars and their clusters as well as the current accelerated expansion of the universe.
The modified gravity approach indicates that the accelerated expansion is due to the influence of modification in gravity for late/early time universe. Some generalizations have been proposed for GR since its inception, most of them have not passed the test of time. It is worth mentioning that any reasonable gravity theory should reduce to Newtonian (N) gravity for a slowly moving weak source. An additional degree of freedom exist generically in any modification of GR. There exist different modifications in the Einstein-Hilbert (EH) action describing the dark components (i.e., dark energy and dark matter) of the current accelerated expanding universe. Dark energy is used to explain the cosmic speed up and dark matter is used to explain the emergence of large scale structures in the universe. Thus, a wide class of gravity theories exist to study the dark side side of the universe via the enhancement of the gravitational force. The f (R, T ) gravity theory is one of such generalization to Einstein's theory of gravity which constitutes on the matter and geometry coupling. In this theory, the Lagrangian for EH action includes the extra degrees of freedom along with trace of stress energy tensor.
The exploration of instability regimes for collapse process strengthened the study of astronomical and astrophysical theories. The gravitational test through pulsar-timing experiments have motivated to study the stability issue. Harada [1] presented the stability analysis in scalar tensor theory for spherically symmetric star configuration and extracted the range of instability from the first order derivative of coupling function. In the study of gravitational instability theory, the amplifications in the density perturbations are responsible for the generation of cosmic structures in the early universe. For a time dependent mass density, the instability for a gravitating system has been discussed in the framework of f (R) theory [2]. Bamba et al. [3] have studied the matter instability describing the curvature inside the sphere in f (R) gravity theory, which is one of the most important criteria to check the validity of any modified gravity theory. It has been the matter of interest for relativists that dynamical equations are developed due to the Bianchi identities and Einstein's equations. Also, an exact evolution equation for Lagrangian can be obtained through Bianchi identities representing the gravitational tidal field. Sharif and Yousaf [4] have explored instability conditions of stellar systems at both N and post-Newtonian (pN) eras with different backgrounds in modified gravity theories.
Nojiri and Odintsov [5] studied the behavior of modified gravity on some solar systems test for Newton law corrections and matter instabilities with the effective cosmological constant regime in the late and early universe. Tiret and Combes [6] presented the dynamical evolution of spiral galaxies in modified dynamics, which is compared to the gravity with dark matter with numerical simulations. Bogdanos and Saridakis [7] explored the perturbative instabilities by imploding scalar and tensor perturbations within a flat background in Hořava gravity. Some astrophysical test for modified gravity theories have been presented by Jain et al. [8] using low-redshift distance indicators to carry out tests and mainly focused on particular stages on the evolution of supergiant and giants to observe distinct observational signatures.
The motion of the matter can be characterized by the fluid parameters like four-acceleration, shear tensor, expansion scalar and vorticity tensor (which is zero for spherical stars). The significance of shear scalar and its vanishing have been brought forward by many researchers for self-gravitating stars. The expansion scalar measures the change in the volume element of fluid configuration during the evolution and its absence leads to the formation of a cavity within the system. This is based on the reason that during the evolution the system is expanding leading to an increase in the volume element due to the increase in the external boundary. The increase in the volume is compensated by the formation of a cavity inside system by imposing the the expansion free condition and the innermost shell would be away from the center.
Initially, Skripkin [9] analyzed the appearance such kinds of vacuum cavity during the evolution of spherically symmetric models which has significance in the modeling of voids. Later, it was found [10] that the Skripkin model has no compatibility with the Darmois matching conditions [11]. In the same paper, they also examine that spherical stars evolving under the expansion free condition must have inhomogeneous energy density. The formation of cavity using the kinematical quantities different from zero expansion condition has also been investigated in the literature [12]. Herrera et al. [13] also explored the instability eras for spherically symmetric collapsing model due to zero expansion in the fluid configuration. Sharif and Bhatti [14] explored instability conditions of cylindrically symmetric self-gravitating systems coupled with charged expansion-free anisotropic matter distribution.
In the study of stellar structures, it is common to model the star interior with perfect fluid which implies the same pressure in the interior of the compact object. Some theoretical advances in recent decades indicate the deviation from isotropic pressure particularly in high density regimes to study their properties. Weber [15] proposed that strong magnetic fields play role for generating pressure anisotropy inside a compact star. It is also observed that anisotropy of pressures is present in wormholes [16] or gravastars [17], so-called exotic solutions of the field equations. General relativistic stellar models gain crucial significance due to the existence of pressure anisotropy in the matter distribution [18]. Sharif with his research fellows [19] investigated different effects of physical parameters on the dynamical instability of selfgravitating collapsing stars. The existence of various compact objects in the realm of f (R) gravity have been investigated [20]. Recently, Yousaf et al. [21] explored the importance of energy density inhomogeneities in the study of stellar collapse.
The current paper presents a full analytical approach have been initiated to understand the instability regimes of cylindrical object within the physical background of f (R, T ) dark sources. The format of this paper is as follows. The next section explores some basics of f (R, T ) theory of gravity including the modified field equations and kinematical quantities. In section 3, we have provided the perturbation scheme up to first order to analyze the stability of our gravitating source. In section 4, we have found the collapse equation by exploring the dynamical equation c.g.s unit systems. Section 5 investigates the instability ranges under N limit with the zero expansion condition. The last section concludes our main findings.

The f (R, T ) Gravity and Cylindrical Systems
The notion of f (R, T ) gravity as a possible modifications in the gravitational framework of GR received much attention of researchers. This theory provides numerous interesting results in the field of physics and cosmology like plausible explanation to the accelerating cosmic expansion [22,23,24]. The main theme of this theory is to use an algebraic general function of Ricci as well trace of energy momentum tensor in the standard EH action. It can be written as [25] S where g, T are the traces of metric as well as standard GR energy-momentum tensors, respectively while R is the Ricci scalar. There exists variety L M in literature which corresponds to particular configurations of relativistic matter distributions. Choosing L M = µ (where µ is the system's energy density) and varying the above action with respect to g αβ , the corresponding f (R, T ) field equations are given as follows where is the effective energy-momentum tensor representing modified version of gravitational contribution coming from f (R, T ) extra degrees of freedom while G αβ is an Einstein tensor. Further, ∇ α represents covariant derivation dT , ∇ α ∇ α and df (R,T ) dR operators, respectively.
The system under consideration is modeled as a cylindrical stellar object whose relativistic motion is characterized by three dimensional timelike surface represented by Σ (e) . This boundary demarcated our manifold into two different interior and exterior portions. These regions are denoted by V − and V + , respectively. The V − region can be described with the help of the following non-rotating diagonal spacetime [26] while spacetime for V + is [28] ds 2 where M is a cylindrical gravitating mass, ν is the retarded time and ζ indicates arbitrary constant. The mathematical formula describing fluid distribution within the cylindrical relativistic interior is [29] T where P φ , P r and P z are stresses corresponding to φ, r and z directions, respectively. Here V β and K β , S β are four velocity and four-vectors, respectively which under the following comoving coordinate system obey some relations. These are given as follows The scalar variable controlling expansion and contraction of matter distribution is known as expansion scalar. This can be obtained through Θ = V α ;α mathematical expression. The expansion scalar associated with cylindrically symmetric relativsitic interior is where over dot symbolizes temporal partial differenration. The corresponding Ricci scalar is where prime means radial partial differentiation. The f (R, T ) field equations (2) for the metric (3) give the following set of equations where Now, we are interested to formulate two equations describing the dynamical evolution of cylindrical relativistic interiors framed within f (R, T ) background. In f (R, T ) gravitational theory, the divergence of stress-energy tensor is non-vanishing and is obtained as The divergence of f (R, T ) energy-momentum tensor yields the following cou-ple of continuity equationṡ These are the required couple of dynamical equations obtained through contracted Bianchi identities of the f (R, T ) effective stress energy tensor. It is well-known that these dynamical equations assist enough to help to analyze the dynamical evolution of stellar system collapse with the passage of time. This also help to explore total energy variation within the collapsing celestial self-gravitating systems in regard with time and adjacent boundaries. Mac-Callum et al. [27] gave a nice way to discuss perturbed boundary conditions in joining matter filled interior with the asymptotically flat vacuum exterior solution. In the above equations, the quantities D 0 and D 1 are functions of t and r and representing extra curvature dark source terms emerging from f (R, T ) gravitational field. The quantities D 0 and D 1 describe f (R) corrections in the energy variations of the collapsing cylindrical relativistic interior associated with time and adjacent surfaces, respectively. These corrections are given in Appendix A.
The matter quantity of cylindrical collapsing stellar geometry can be defined through gravitational C-energy which was proposed by Thorne [30]. This is obtained as where l indicates specific cylindrical length. Before calculating its variations among adjacent surfaces of cylindrical anisotropic fluid distribution, we shall define some operators. The proper and radial derivative operators are defined as follows The relativistic velocity of the collapsing stellar interior can be obtained with the help of proper derivative operator as Using Eqs. (20) and (22), we obtaiñ Next, from Eqs. (20), (21) and (23), it follows that This equation tells us provides the total energy variation between adjacent surfaces within the matter configuration. From the above equation it is evident that the first part of the right hand side is due to energy density and f (R) higher curvature quantities. The last term in the above equation is ( T 01 /BA)(U/Ẽ) < 0. In this term, the quantity U represents collapsing matter velocity. It is well-known that U is less than zero for the collapsing fluid models. Thus, the last term (containing U and non-attractive f (R) corrections) lessens fluid energy influences in the evolutionary system phases. Equation (16) upon integration yields It is known that zero expansion leads to the emergence of boundary surfaces in which one (outer) demarcates the relativistic interior fluid from the exterior spacetime while the second (inner) separates Minkowskian core from the matter distribution. Under null expansion scalar framework, the relativistic fluid evolves without being compressed. For example, during expansion self-gravitating system, the increase in volume of matter configuration leads to expansion of the external boundary surface which can be counterbalanced by a similar expansion of the internal surface to make Θ zero. Thus, zero expansion scalar triggers evolution of relativistic system in such a way that the inner most shell moves away from the central point thereby causing the emergence of vacuum core. Due to this zero expansion, matter sources could be effective for the voids explanation. Voids are, roughly speaking, underdense areas incorporating substantial amount of information on the cosmological environment [31]. Voids provide a reliable guide to study the large scale cosmic structure formation. They are more rich in modified gravity [32] as compared to GR as this extended gravity theory is more likely to host large structures with smaller radii.
The continuity of Eqs. (3) and (4) over Σ (e) can be obtained by using Darmois matching conditions [33] where = means that measurements are performed over outer hypersurface. The matching conditions over Σ (i) lead to m(t, r) To present f (R, T ) gravity as a cosmologically and theoretically consistent theory, the selection of its models is very important. Here, we consider particular class of models as follows This model involves the explicit non-minimal curvature matter coupling. We now consider f (R, T ) power law type model given by Such functional f (R, T ) configurations match to the Lagrangian form mentioned in Eq. (29). Here, we take All GR solutions can be obtained by taking limit λ → 0.

Perturbation Scheme
Perturbation theory gives us a mathematical technique that assists enough to find an approximate solution of a differential equation. After applying perturbation scheme, one can break corresponding equations into "solvable/static" and "perturbed" parts. The impact of perturbed terms in the equation keeps on decreasing, controlled by the perturbation parameter. Here, we consider perturbation scheme that was proposed by Herrera et al. [34]. In this scheme, we take α to be perturbation parameter with α ∈ [0, 1] and consider effects upto first order. We assume that the system initially is in the state of hydrostatic equilibrium. Due to this cylindrical scale factors as well as fluid variables are independent of temporal coordinate. We shall represent static configurations of corresponding variables by zero subscript. Upon perturbations, all these structural variables depend upon the same time dependence η(t), which eventually gives same time dependence to Ricci scalar. The perturbation scheme is where R 0 is the static form of Ricci invariant whose value is The f (R, T ) field equations (8)- (11) under static background with C 0 = r turn out to be where (S) ψ ii indicate static configurations of the corresponding dark source components and are given by Equations (8)-(10) under non-hydrostatic state take the following form where χ ′ i s are dark source terms coming due to f (R, T ) gravity. These terms contain static combinations of metric variables and are mentioned in Appendix A. After applying our radial perturbation approach, we have seen that first dynamical equation (18) is trivially obeyed while the second dynamical equation (19) turns out to be where D 1S is the static form of D 1 and is found to be .
The static as well as non-static portions of cylindrical relativistic C-energy function turn out to be the expansion scalar in its perturbed formulation is The equation relating b and B 0 is found after perturbing Eq.(8) as The first and second conservation laws (18) and (19) after using perturbation scheme provide the following non-static perturbed distributions aṡ where χ 1 constitutes gravitational effects coming from effective energy density of the cylindrical stellar objects, while D 3 contains f (R, T ) higher curvature corrections. These terms are mentioned in Appendix A. Equations (53) and (54), known as dynamical equations, would be very useful in the discussion of collapsing behavior of relativistic stellar interiors. The matching condition (27) in account of perturbation provides P r0 Using (00) field equation in the second of the above equation, we get where Equations (48) and (57), after some manipulation, gives the following second order partial differential equation where Equation (58) has two solutions with two different behaviors and these behaviors are totally independent of each other. In this paper, we are interested to find the unstable constraints of evolving cylindrical compact object in modified gravity. Further, it is mentioned earlier that our system was initially in complete hydrostatic equilibrium. Then, it enters into the collapsing phase by reducing its areal radius. Therefore, we now restrict our perturbations in such a way that all radial perturbed functions, i.e., a, b, c and e, are positive definite, which eventually making ω 2 Σ (e) > 0. The solution associated with Eq.(58) is obtained as follows

N & pN Terms and Collapse Equation
In this section, we shall express second dynamical equation into centimetergram-second (cgs) units and then indicate terms relating to N, pN and parameterized post Newtonian (ppN) epochs. This would be done by expanding cgs second dynamical equation upto O( 1 C 4 ), where C indicates light speed. For N and pN epochs, we shall consider the following approximations where G is the gravitational constant. Equation (45) gives us following peculiar form of double derivative of A 0 with respect to radial coordinate (61) Equations (44) and (51) provide the following first radial derivatives of A 0 and B 0 as We use Eqs.(60), (61) and (63) in static form of second dynamical equation (50) which after converting into cgs system turns out to be It is worthy to mention that expansion of the above equation provide terms relating to some specific eras with details as follows terms of O(C 0 ) indicates contribution at N era, The complete expansion of above equation The f (R, T ) field equations are filled with complicated derivatives of radial and temporal coordinates, the exploration of their generic solutions is a painstaking task that yet now has not been accomplished. However, certain restrictions with some physical background could assists enough to find their solution. In this perspective, we consider expansion-free evolution of cylindrical compact objects against linear perturbation. The expansion-free condition can be obtained from Eq.(52) as follows We are now interested to calculate f (R, T ) expansion-free cylindrical collapse equation. This would be furnished by considering the second perturbed dynamical equation along with the supposition that cylindrical evolution is supported by null expansion scalar background. Thus, using Eq.(54), junction conditions (55), (56), (59) and above expansion-free constraint, we obtain the following modified collapse equation over the exterior boundary surface where quantities Φ and Ω e are described in Appendix A. It is worthy to stress that in the above equation, the term Φ describes the gravitational contribution in expansion-free systems. This contribution is same as in expansion cylindrical systems evolution. However, the quantity Ω e is responsible for the emergence of cylindrical central Minkowskian core during system evolution. The gravitational effects coming from Ω e causes the naked singularity appearance during cylindrical collapse. This comes from the fact that the quantity Ω e contains all those terms that have been evaluated by performing expansion-free condition. When the system fluid moves inward, during collapsing phenomenon, with null expansion rate, there will be a blowup of the shearing scalar at the center. It is well-known from the work of Joshi et al. [35] is that strong shear could produce hindrances in the appearance of apparent horizon, thus producing a platform for naked singularity formation. Thus, Ω e gave a way to discuss naked singularity emergence in a simple way. On making b B 0 = − c r , one can remove all of the expansion-free effects in the above collapse equation.

Instability Constraints at Both N & pN Epochs
In this section, we explore unstable regions of collapsing cylindrical stellar model at both N and pN eras. Under pN limits, Eq.(66) in relativistic units yields For the onset of instability, one needs to satisfy above relation. Since majority of the above terms are positive, therefore instabilities will appear because of negative terms in the above equation. For that reason we consider the following constraints to be obeyed. .
These are the required instability constraints that the system must satisfy in order to enter in the unstable window during evolution. In other words, the cylindrical relativistic anisotropic interior will be unstable as long as it obeys above relations. It is worthy to stress that in order to consider above equation (67) as an instable constraint, we need to consider that all the terms on both sides of the equation are definite positive. Thus, we take Φ > 0 and Ω e > 0. It is seen that this constraint depends merely on the static configurations of fluid as well as f (R, T ) variables. Now, we consider N order effects, then Eq.(67) reduces to Now, we use constant curvature condition firstly in above equation and then in Eq. (24). After using simultaneously these equations, we get where tilde shows that terms are evaluated under constant curvature condition. Now, we use ansatz µ 0 = ξr n in which n ∈ (−∞, ∞), while ξ is any positive real number constant. This type of ansatz is well justified because any function of a single variable can be expanded in a series form of that variable. Considering n = −3 and n = −4, we get from the above equation This is the required instability constraint at N era for the cylindrically symmetric collapsing systems framed within f (R, T ) gravity. This indicates that our relativistic instability constraint depends on radial dependant fluid and f (R, T ) model variables. It is worthy to stress that this instability constraint is independent of stiffness parameter, which generally has utmost relevance in the discussion dynamical instability of any stellar object. Here, we also need to suppose that all quantities coming on the both sides of the above equation are non-zero and non-negative.

Summary and Discussion
In this paper, we have explored some dynamical constraints which are essential for a cylindrical object to be in a physically stable state. We have constructed our analysis quiet systematically by forming the modified field equations within the background of f (R, T ) gravity theory. The cylindrical system is chosen to be filled with anisotropic matter in the interior while the exterior region is considered vacuum. The dynamical equations are obtained from contracted Bianchi identities and some useful kinematical variables are also explored including the expansion scalar. We have presented the linear perturbation technique for metric and matter variables with some known static profile of cylindrical object. Initially, the system is assumed to be at rest and then gradually enters into the non-static phase with the same time dependence on the metric coefficients.
We have perturbed all the relevant equations to construct the collapse equation. The general collapse equation is obtained by using the conservation laws and field equations up to first order in perturbation parameter. Also, we have constructed a real static solution describing a collapsing state at large past time. For which we have considered that all the metric functions in the static background are positive indicating the cylindrical line element to be Lorentz invariant. Moreover, we have categorized N, pN and ppN eras by expanding the collapse equation up to order of C −4 .
For this purpose have converted our system to the c.g.s unit systems because we know that every term related with some power of speed of light have some physical interpretation and describes some useful regimes. Since the terms associated with the zeroth order of speed of light, i.e., C 0 corresponds to N era and the terms linked with C −2 provides the information of pN regime. Similarly, the terms appearing with the order of C −4 present the era of ppN. Further, we have imposed the expansion free condition on the collapse equation due to physical significance of this constraint. Such systems are consistent with those astronomical objects which have an inner cavity after the central explosion. A physical application of our study is possible in those astrophysical objects which have cavity in the interior region. This is due to the fact that cavity formation in the the expansion free case is compensated by the increment in the boundary surface during the overall expansion.
Generally, the adiabatic index (Γ) which describes the rigidity in the fluid distribution indicate the instability regimes for a gravitating source in the presence of expansion scalar. Particular values of Γ (i.e., Γ < 4 3 for spherical systems and Γ < 1 for cylindrical objects with perfect matter in the interior) exists in literature indicating the unstable phase of the relativistic body. However, in the absence of expansion scalar, this factor Γ have no such importance to evaluate the unstable regions of the relativistic system. To examine the unstable regions of a self-gravitating cylindrical object with zero expansion, it should satisfy the requirements (67) and (70). The violation of these constraints describes the stable configuration of the collapsing system. We found that the instability range depends upon the dark source terms originating due the f (R, T ) theory of gravity as well as on the length of the cylinder. Moreover, the material profile like anisotropic pressure and energy density also control the stability of the object during the evolution.
The exploration of the zero expansion condition could be closely linked with the study of voids which are sponge like structures and can be explained with the Minkowskian cavity inside it. Thus, the zero expansion condition asserts the existence of Minkowskian cavity at the center. The potential applications of our dynamical analysis are present in those astronomical objects which are carrying a central Minkowskian cavity. Finally, we would like mention here that all our results correspond to the instability constraints obtained in GR under the particular limit, i.e., f (R, T ) = R [14].

Appendix A
The dark source terms D 0 and D 1 of Eqs. (18) and (19) are given as follows The components of extra curvature terms mentioned in Eqs.(47)-(49) are .