Implications for vanishing complexity in dynamical spherically symmetric dissipative self-gravitating fluids

The complexity factor, originally based on a probabilistic description of a physical system, was re-defined by Herrera et al. for relativistic systems. This involves an assessment of the energy density inhomogeneity, anisotropic and shear stresses, and in the case of radiating collapse, the effects of heat flux. Already well integrated into the modelling of static configurations, the complexity factor is now being studied with respect to dynamical, self-gravitating systems. For static systems, the constraint of vanishing complexity is typically used however for the non-static case, the physical viability of the vanishing condition is less clear. To this end, we consider the ideal case of vanishing complexity in order to solve for the time-dependent gravitational potentials and generate models. We find that vanishing complexity constrains the metric to be of a form similar to that of Maiti’s conformally flat metric.


Introduction
The process of gravitational collapse of massive objects is a challenging topic in relativistic astrophysics, having been pioneered by Oppenheimer and Snyder [1]. As the end state of the collapse process is approached, it is likely that the physical quantities calculated via modelling processes become inaccurate and less representative of the actual physics at work. Certainly, in the case where the remnant formed is a black hole, an accurate description of the final stage of collapse might be elusive. This is not due to shortcomings in the theory of general relativity, but rather a result of the modelling process which often involves setting up initial static configurations [2,3] with associated boundary conditions and equations of state [4]. Initial conditions can also include an a e-mail: bogadi.robert@gmail.com (corresponding author) b e-mail: megandhreng@dut.ac.za c e-mail: moyos@dut.ac.za inhomogeneous energy density [5] which is relevant to our study involving complexity. Constraints based on a quasistatic system in the distant past might lose relevance as the collapse process proceeds. The progression of the collapse is also assumed to proceed in a well-behaved manner, mainly governed by the heat flux boundary condition of Santos [6] in the case of non-adiabatic collapse. So far, this has been successful in generating models and may well be accurate for most of the collapse process. However, near the time of remnant formation where energy densities and pressures escalate to extreme measures, additional non-linear effects might arise, rendering the initial conditions and assumptions made in setting up the model of little relevance. In an effort to generate models, other physical constructs such as gravitational decoupling [7] and complexity factor analysis are now being considered. A definition of complexity from the statistical point of view was developed by López-Ruiz et al. [8] and involves the interplay between entropy and information content. In the case of relativistic objects, it is more suitable to seek a measure of complexity in terms of the field equations of relativity. This was done by Herrera [9] for static systems and also for the dynamical case [10] which can by applied to the problem of gravitational collapse. With respect to relativistic objects, complexity was considered to assist in discerning equations of state of static compact objects [11]. The complexity factor focuses on the degree of energy inhomogeneity and can be coupled with pressure anisotropy, shear stress and in the case of radiating collapse, the heat flux. Much work has already been done in the case of static systems in which the constraint of vanishing complexity has been implemented. In the case of dynamical systems the situation is less clear and the complexity factor might likely have a specific, non-vanishing behaviour at certain stages during collapse. Initial investigations have already been conducted in this respect [12]. In this work, we assume the contrary that the dynamical complexity factor vanishes, in order to establish the restriction on the models obtainable under this regime. The framework has then been set up for future work whereby the complexity factor could be linked to a physically desirable non-vanishing function of spacetime. This could then lead to more suitable models in which improved stability and more realistic physical features are obtainable.

Dissipative collapse
For the modelling of a radiating star undergoing shear-free gravitational collapse, the interior spacetime is described by the spherically symmetric line element in which the metric functions A and B, describing the gravitational potentials, are yet to be determined. This shear-free metric was studied in detail by Bonnor et al. [13]. The energy momentum tensor for a fluid which dissipates heat to the exterior is given by where ρ and p are the energy density and isotropic pressure of the fluid. The heat flow vector q a is orthogonal to the velocity vector so that q a u a = 0. The Einstein field equations describing the interior of the stellar fluid are where the dot and prime are derivatives with respect to time and radial distance respectively. We obtain the condition of pressure isotropy by equating (4) and (5) and we note that this is time-independent. Time dependence can then be introduced by promoting the constants of integration to functions of time. This was recognised by Ivanov [14].
Since the star is radiating energy, the exterior spacetime is described by the Vaidya metric [15] where v is the retarded time and m is the total mass inside the comoving surface Σ forming the boundary of the star. Matching of the interior line element (1) to the exterior spacetime (8) is achieved via the junction conditions, originally proposed by Santos [6]. These are where m Σ is the total mass within a sphere of radius r Σ and (10) represents the conservation of the momentum flux across the boundary Σ.
Lastly, we note that the shear-free approximation is prone to instability and in some respects, pressure anisotropy assists in overcoming this problem [16]. In our application to modelling, we evaluate the adiabatic index in order to highlight this aspect.

Complexity
The complexity factor for radiating, self-gravitating systems, as defined by Herrera [10] is given by (12) which in terms of our shear-free metric and associated field equations (3)- (6) gives The integration in (13) can be completed to yield 3.1 Isotropic systems In the case of isotropic systems (δ = 0), we find that vanishing complexity is ensured if where C(t) is a promoted constant of integration. A general solution for B is then given by where R(t) and k(t) arise from integration and transformation.
The pressure isotropy condition, obtained via the Einstein field equations, gives which may be integrated after substituting the form for B(r, t) in (16). A solution is then given by where ξ(t) and ζ(t) are further constants of integration which have been promoted to functions of time. These are to be set according to boundary conditions, equations of state, or additional physical constraints.

Anisotropic systems
For a more general, anisotropic system (δ = 0), the complexity factor (14) is given by and for vanishing complexity, a first integration yields which reduces the problem of finding solutions to a singlegenerating function.

Application to modelling
The application of the vanishing complexity condition has become an integral part of modelling relativistic objects within the static regime. For dynamical systems, the complexity factor as given according to the above definition is less likely to vanish especially in the case of gravitational collapse [12]. This is expected due to the causal nature of dynamical processes whereby the establishment of an equilibrium with respect to pressure anisotropy, energy density inhomogeneity and heat flux cannot occur instantaneously. Nevertheless, the vanishing complexity condition might still be valid for quasi-static systems or in our case, at early stages in the collapse process. We see that for the isotropic case, vanishing complexity results in metrics which are similar to those of Maiti and Bergmann [17][18][19]. On application to compact, radiating fluid spheres, the field equations give homogeneous energy density profiles as is the case with the simple Schwarzschild interior solution [13,20].
Maiti's solution can represent a fluid with heat flux in a conformally flat spacetime and the general form is given by Banerjee et al. [21] also considered this metric and explored its generality in terms of conformal flatness. We make use of the results obtained from vanishing complexity for the isotropic case and set the gravitational potentials as follows, where f (t) is some function of time, and the static part is given by We thus chose (k(t) → 1 and R(t) → 2R f (t)) in (16) and (ξ(t) → f (t) and ζ(t) → ζ f (t)) in (18). This is almost that of the Maiti form, apart from the modified time-dependence on metric function A. Metric functions with the above separable form in radial distance and time have been utilised by Tewari [22] and is a more general form of that used by Bonnor et al. [13].
The associated field equations for the above system (21)- where ρ s and p s denote the energy density and pressure respectively of the initial static configuration. These are given by 8π p s = 1 Using the boundary condition ( p s | Σ = 0) for the initial static configuration, we obtain Considering the time-dependent, collapsing phase, the appropriate boundary condition is ( p Σ = (q B) Σ ) which provides an equation for the temporal dependence, A first integration of (31) yieldṡ and further integration, We now consider the collapse of a 5M shear-free fluid of initial radius r Σ B 0 (r Σ ) = 2.159×10 3 km and initial density ρ 0 = 2.363×10 8 g/cm 3 . The initial density then determines the value of R to be R = 2.608 × 10 4 km. These parameters may be compared with data from the extensive and in-depth study done on Supernova 1987A [23]. The associated model parameters and calculated quantities are given in Table 1.

Discussion
As noted in application to modelling, imposing a vanishing complexity constraint results in metric line elements of the Maiti and Bergmann form. The Maiti form was selected with the gravitational potentials set in variable separable form, providing a static configuration consistent with the Schwarzschild interior solution. This method has been used by Pinheiro and Chan [24] with metric functions chosen such that only the spatial part incorporated time dependence. We followed the approach of Tewari and Charan [25] in which both metric functions A and B are spatially and temporally separable. In particular, the time dependence was not removed from gravitational potential A(r, t) = A 0 (r ) f (t) as it is common to set A(r, t) → A 0 (r ) [13,24]. By including time dependence, we incorporate horizon formation and comparison may be made with systems which follow the perturbative approach [26]. Gravitational collapse models with heat flow and without horizon formation have also been studied [27]. From Table 1, we see that horizon formation occurs at f ≈ 0 when most of the mass (99.3%) has been radiated as seen in Fig. 1. In Figs. 2 and 3 we see that the energy density and pressure become exceedingly large within the last 100 ms. The luminosity peaks at about 57 ms before horizon formation as seen in Fig. 4, which is about two orders of magnitude slower than that predicted for the gravitational collapse of unstable neutron stars [4]. It could be that a non-vanishing complexity might provide for additional processes or pathways which are conducive to more rapid gravitational collapse. The inclusion of a string field has already shown that inhomogeneity promotes earlier horizon formation [26]. In addition, nonvanishing complexity might assist with stability. In Fig. 5, it is evident that our radiating model suffers from instability towards the centre and instability appears to grow as the collapse proceeds. This is in contrast to the adiabatic case in which isotropic pressure and homogeneous energy density lead to the stability of the vanishing complexity condition   [28]. It has been noted that the condition of zero complexity is affected by heat dissipation [29] and initial investigations into non-vanishing complexity for self-gravitating systems have recently been made [12].

Conclusion
We have generated a model for the collapse of a selfgravitating, shear-free and isotropic fluid with heat flux, based only on the constraint of vanishing complexity. No initial static configuration was used apart from the specification of the initial mass and radius. Vanishing complexity restricts the gravitational potentials to that of a Maiti type metric. The potentials were assumed to be separable in radial and time coordinates and a temporal behaviour similar to that of Tewari was followed. In other studies, temporal dependence is only included on the spatial part of the metric which can lead to models without horizon formation. In our model, temporal dependence on both metric functions ensured the formation of a horizon, albeit at a point where most of the mass has been radiated. We also concluded a relationship between the metric functions for the anisotropic case with vanishing complexity. Specification of one of the metric functions is then required, allowing for the inclusion of an initial state configuration. The modelling process then depends on the successful determination of the other metric function. It is expected that anisotropy would assist in improved stability. Lastly, future work on the complexity factor might result in non-vanishing trends in the factor during certain stages of the collapse process. The resulting differential equation in the metric functions would no doubt be much more complex and likely require numerical techniques.
Acknowledgements RB and MG acknowledge support from the office of the Deputy Vice-Chancellor for Research and Innovation at the Durban University of Technology.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data was obtained using the formulae explicitly given in the 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 copy-