Computational Algorithms for the Study of Distributions with Electric Charge and Radiation Flux in General Relativity

: The following research contributes to the understanding of how electric charge inﬂuences the unwrapping of spherically symmetric distributions using a relativistic approach. With the recent detection of gravitational waves and the location of black holes, where a strong emission of electromagnetic radiation prevails, interest has arisen to consider the existence of compact charged objects. Thus, the appearance of charge in small quantities during gravitational collapse, the process by which black holes originate, is not ruled out; this also includes the emission of electromagnetic waves from them. This article intends to establish algorithms and write ﬁeld equations for a charged ﬂuid as those corresponding to an anisotropic ﬂuid with radiation ﬂux. Using an appropriate deﬁnition of the mass function, considering self-similar symmetry and Bondian observers, dynamical results are obtained for the Einstein–Maxwell electromagnetic system with added gravity. Imposing an additional homothetic symmetry, the ﬁeld equations are solved, and the most relevant conclusions are drawn about the inﬂuence of the electric charge during the collapse and subsequent changes in the physical variables.


Review of the Field Equations Posed by Einstein
Field equations developed by Einstein, which form the basis of the theory of general relativity (GR), relate the geometry of space-time (Einstein Tensor) to the distribution of matter (energy-momentum tensor) (Equation (1)) GR replaced the theory of universal gravitation by explaining effects that cannot be accounted for by Newton's laws, such as the anomaly in the orbit of Mercury and other planets. It also makes numerous predictions, some already formed about the effects of gravity, such as the bending of light and time contraction. In addition, it predicts space-time distortions, a phenomenon known as gravitational waves.
The study of a system composed of a radiating fluid is possible through these equations. The solution of these equations allows, on the one hand, to obtain a set of physical variables, density, pressure, fluid velocity, and energy flow, and, on the other hand, information about the geometrical structure of space-time. Following a literature-based approach, the problem can be analyzed by considering that the surface of the distribution separates space-time into two different sections: the section inside the distribution, referred to as the inner zone, where there are high concentrations of both matter and energy, and the section outside the surface, referred to as the outer zone, with radiation freely escaping to infinity. It is, therefore, of great importance to conceptualize the inner and outer regions properly from a geometrical point of view, as well as the characteristics of their energetic material content.
Some physical characteristics of collapsing charged matter distributions may justify the interest in considering whether the electric charge affects gravitational collapse. For example, at certain times of collapse, electric charge in reduced amounts can change the final state of an astrophysical object.
The study of spherically symmetric charged distributions begins with Eddington [1], who showed that, if an isolated point charge has an oscillatory potential, the Lorentz gauge [2] allows all components of the metric to be radiative, even if they are static. Although some refer to static situations [3], significant efforts have been made to describe dynamic situations. Recently, great interest has arisen in studying the effect of large amounts of charge in spherically symmetric systems, although by an approximation in pre-relativistic physics it can be proved that a physical object possesses a quantity of charge greater than one hundred coulomb over M Θ (solar mass) [4].
It is important to highlight the possibility that strong electric fields may exist in strange stars [5]. Some authors [6] describe models of collapsing newly charged spheres including dissipative processes and emphasize the importance of electric charge in the interactions present in these systems. Specifically, they show that the electric charge does not cause regeneration effects (it does not add more activity to the present field).
To solve problems in GR, it is necessary to introduce certain initial conditions, such as thermodynamic principles, a correlation between measurable variables, additional symmetries to space-time [7][8][9][10], or to use the post-quasistatic approach. For example, Ray et al. [11] study charged fluid spheres in an Einstein-Maxwell space-time under the imposition of conformal symmetry. In a study on relativistic stars [12], it was observed that the event horizon occurs when the radius reaches its limit; it is at that time that the gravitational attraction can be balanced in the presence of a large amount of electric charge [2,13,14]. Bonnor and Wikramasuriya [15] investigated the effects of electric charge on the matter, whose electrostatic repulsive force balances the gravitational attractive force [15,16].
Other research on electric charge includes that of Barreto and Da Silva [13,14], who study self-similar charged spherical distributions in the diffusion approximation. Ray et al. [11] use classical arguments to check the electric neutrality of stars and show that, if the pressure, density, and gravitational field are large, it means that high values of charge and intense electric fields could exist in such distributions. Ray et al. [11] analyze the influence of electric charge on compact objects (in that case, stars) and observe how the equilibrium of forces allows a charge in them of about 10 20 coulomb and electric fields around 10 20 coulomb and electric fields around 10 21 V/cm [11,17,18].
Barreto et al. [17,18] noted that solving the differential equations defining these problems can be used to find the physical characteristics of collapsing stars. These studies show that the exact solutions of electromagnetic systems with coupled gravity are important in the representation of compact objects, such as quark stars or preon stars. Some researchers [19] found solutions to these systems where charged spheres are present with a specific type of field, both electric and gravitational. These researchers [20] showed the solutions they obtained for these types of arrangements where the electric and gravitational fields have a specific shape.
In studies of compact relativistic objects, strange stars have been proposed as a different option to the well-known study of neutron stars, and anisotropy (intrinsic, electric charge-induced, or by any other mechanism) plays an important role when analyzing these objects. If the anisotropy is induced by electric charge, naturally, the Tolman-Oppenheimer-Volkoff Equation (1) [21] should be modified to consider the charge within the distribution. Among other possibilities, one of the reasons for considering anisotropy within stars could be the presence of strong electric fields as recently suggested by Usov [5].
Among other works related to electric charge, we can mention that of Rosales et al. [17], where they study static fluid spheres in isotropic coordinates and find a central equation of state for the distribution. Ivanov [15] demonstrates the relationship between anisotropy and charged fluids.
A comprehensive study of the influence of electric charge on anisotropic strange star aspirants takes into account an interior measure of space-time with spherical symmetry. To solve the resulting field equations, it is observed that the cumulative electric charge distribution is maximal over the surface possessing an estimated electric charge in the range of 10 20 coulomb and an electric field in the range of 10 21 V/cm a 10 22 V/cm [22].
In [17,18], new ways were obtained to solve the field equations where matter possesses both anisotropy and electric charge. When a spherically symmetric and static space-time is assumed, an equation of state of quadratic type (Equation (2)) is used to know the different gravitational potentials and the electric field strength, through which the abovementioned new exact solutions for isotropic graphical charts are obtained. In [13,14,16,18], an algorithmic method was disclosed that could model the behavior of all fluid distributions possessing charge with spherical symmetry, considering other different function, and a set of physically acceptable solutions was found, and representative compact stellar models were reported.
The study of a charged compact object with anisotropic pressures and core governed by a quadratic equation of state was analyzed in [23], and it was concluded that the masses, radio, and other compactness factors of stellar objects are affected by the presence of an electric field when these values agree with those of previous studies [24].
This work attempts to delve into the study of charged, radiative, and spherically symmetric fluids. An important aspect to consider is that the question of whether charged relativistic spheres form black holes remains unanswered. In this work, we make an interesting contribution in this sense, solving the equations of the governing electromagnetic systems that also include gravity (Einstein-Maxwell systems) and considering both the equation describing the conservation of charge and the Tolman-Oppenheimer-Volkov equation [21].

The Metric, Energy-Momentum Tensor, and Field Equations
The metric for defining the shape of the distribution on which the profiles of geometrical and physical variables in general relativity will be evaluated; the choice of a suitable coordinate system greatly simplifies the analysis of the system under study. In particular, in the study of radioactive systems, the description provided by the so-called radiation coordinates or Bondi metrics [25] is adequate. Therefore, for the description of the geometry of the interior of the distribution, the Bondi metrics [25] are used. Considering a spherically symmetric matter distribution, which is given by a charged fluid, with values of pressure, energy, electric energy density, and radiation energy flux given by p, ρ, ǫ, and , respectively, and measured from the perspective of an observer in a Minkowski representation moving at par (comoving) with the fluid (with velocity −ω), the metric is expressed in a radiation coordinate system [26][27][28] as Equation (3): of which β is a function u, and V is a function of r, the Bondi time being u = x 0 [25] at null infinity, where u is constant. Where r = x 1 is the radial coordinate tangent to the null cones. Moreover, both θ = x 2 as well as ϕ = x 3 correspond to angular coordinates. Geometric units will be used throughout this work. G = c = 1. From the metric function V, the function m (u, r) is defined, which represents the mass term in the interior of the distribution. This quantity was introduced by Bondi [29] as a generalization of the mass in the static case to the radiation case and is related to β and for the case of charged distributions by the following expression: where C represents the charge inside the distribution. Knowing that, in spherical symmetry, the only components with nonzero values of the electromagnetic field tensor are F ur = −F ru (Equation (4)). Then, a function C(u, r) can be defined by the correspondence observed in [30][31][32].
The nonhomogeneous Maxwell equations are Equations (6) and (7) C ,r = 4πr 2 J u e 2β (6) and The partial derivative as a function of the indicated coordinate is denoted by the comma symbol (,); then, the time component is represented by J u and the radial component by J r , both belonging to the electric current density quadrivector; then, the idea is that C (u, r) shows a function of charge at radius r at a given time u. In this way, the conservation of charge in the inhomogeneous electromagnetic equations is expressed as well as observed in Equation (8) U where the quadravelocity of the fluid is given by (9).

Energy-Momentum Tensor
We now follow the approach of Bondi [29][30][31][32], reported as Bondian observers of the GR [9]: a local observer is located, according to the Minkowski perspective, in motion with the fluid charged with velocity -ω; for such an observer following the radial direction (the components of the energy-momentum tensor are given by the energy density, the electric energy density, the unpolarized energy density, and the pressure represented by the symbols ρ, µ, ǫ, and p, respectively), Equation (10) will represent the covariant energy-momentum tensor.
where density ρ, pressure p, metric tensor components g µν , radiation flux , null vector l µ , and electromagnetic tensor E µν are identified. For this paper, it is more useful to write Equation (10) as shown in Equation (11).
where they are identified as densityρ and effective pressurep, respectively, andμ = E 2 /8, with the electric field given by E = C/r 2 and p t = p +μ. As mentioned above, C(u, r), is the function representing the charge for any r at time instant u and satisfies the conservation of charge described in Equation (12) for numerical purposes.
As shown above, the comma (,) is the partial derivative concerning the indicated coordinate. Equation (13) represents the velocity of matter.
When the fluid is adiabatic, it possesses a diagonal energy tensor in covariant form (ρ, p r , p t , p t ) definingρ = ρ +μ, p r = p −μ, which is the same for an anisotropic fluid [16,17]. It is clear that, in the presence of electric charge, local anisotropy occurs. On the contrary, in the absence of electric charge, the Pascalian character of matter, i.e., isotropic fluid, is recovered. The variables P t and P r are the tangential and radial pressures, respectively.
In addition, by defining the mass function as shown in Equation (14), the so-called Missner-Sharp mass generalization [33,34] is obtained. When there is charge, it is expressed as
Making adjustments,ρ = ρ + , p r = p −μ ,μ = E 2 /8π, with E = C r 2 y p t = p +μ. The field equations can be written as Equations (20)-(23). where By knowing the metric β(u, r) and the metric m(u, r), together with their derivatives, it is possible to determine p, ρ, , and ω from the symmetry equations. Equations (20)- (23) are the same for anisotropic matter [18,[33][34][35], so that it can be interpreted as anisotropy of the fluid, to the electric charge. Such an interpretation is possible because of how the mass (Equation (14)) is defined, so that the electric charge can be assumed as a constituent of the metric. Assuming this position, the electric charge would then become another physical variable that contributes to both energy density and pressure. This model helps with the understanding of what has been studied due to the approximation of the Bondian observers to the GR (comoving frame of reference). Barreto et al. [32,36] showed that the observers used by Bondi for the treatment of matter in GR and those used by the standard ADM 3 + 1 of numerical relativity are the same, both have comoving and Minkowskian-type observers, as mentioned above. These observers are the same as those of the post-quasistatic approximation of GR for the study of collapsing spheres. The treatment is identical for an intrinsic anisotropic fluid, except for the mass function.
At this point, there are four field equations for seven unknowns: two geometric variables m and β and five physical variables p, rho, radiation flux, velocity, and charge.
The solution algorithm is as follows: - The geometric variables are obtained from the imposition of an additional homothetic (self-similar) symmetry, i.e., from the symmetry equations as a function of the surface variables. -Previously, to find the surface variables, the surface equations must be solved by the Runge-Kutta method of order 4 [36]. - The electric charge is determined from the numerical solution of the charge conservation equation. -Pressure, density, radiation flux, and velocity are determined from the field equations.

Results
The case study of a system composed of a radiant fluid with electric charge is possible through the field equations describing the Einstein-Maxwell system. The solution of these equations allows, on the one hand, for obtaining a set of physical variables, density, pressure, fluid velocity, energy flow, and electric charge, and, on the other hand, information about the geometrical structure of space-time. As already mentioned, in studying this problem, it is necessary, according to the theoretical approach, that space-time is divided in these configurations by the surface of the distribution into two regions: an inner one containing matter and radiation and an outer one where radiation escapes if difficulty of the system moves toward infinity. It is therefore necessary to show the geometry that defines both the inner and outer regions, as well as to characterize their energetic material content.

Reissner-Nordstrom-Vaidya Space-Time
The spherically symmetric charged distribution is considered to be bordered by the surface ∑. Whatever is not contained in it belongs to an RNV-type space-time described in Equation (26), according to [16,17].
In Equation (26), the total mass is given by M(u), and the total charge is C. Smoothly coupling Equations (2) and (26) on the hypersuperciliation r = ∑.
Once this moment is reached, which is the Bondi approximation for spherical symmetry, the equations obtained are of a general type. We identify five variables (ρ, p, ω, , C) that are physical, plus two (m, β) metric functions and four field equations in addition to the conservation equation. Other considerations need to be added in order to solve the characteristic initial value problem. Under these circumstances, space-time is assumed to be self-similar, which will allow us to illustrate how the transport mechanism influences the evolution of charged spheres.

Self-Similarity and Surface Equations
The inner zone, characterized by Equation (2) and the outer zone, described by the Reissner-Nordstrom-Vaidya metric [16,17], are separated by the hypersurface r = ∑. Therefore, the coupling of the inner and outer solutions of the field equations must be considered. Furthermore, to know the evolution of the physical variables in the interior of the distribution (p, ρ, C, , ω) in order to evaluate quantities over the surface, it is necessary to construct a system of ordinary differential equations relevant.
Having done the above, a system of three ordinary differential equations is obtained, with three functions of time µ associated with the radius of the distribution, the velocity at the edge of the distribution, and the function m evaluated at the surface.
Some authors [37] characterized covariantly the self-similarity when the matter is distributed spherically, noting that the self-similar solutions depend on whether they are invariant or not in the presence of scale transformations. Seen mathematically, the interesting thing about self-similarity is that it helps to reduce the field equations to a group of ordinary differential equations, and, on the other hand, because it helps to generalize the symmetries, homothetic, and Killing symmetry, according to [38][39][40].
To define self-similarity requires the existence of a homothetic Killing vector field. In a given coordinate system, self-similarity is evidenced by a simple scaling relation for the metric functions. Probably no problem in GR has been more exhaustively studied than that of spherically symmetric homothetic space-times [37]. Now, an additional homothetic (self-similar) symmetry is imposed, which mathematically is expressed in Equation (28) [16][17][18].
L ξ g µν = 2g µν (28) where g µν is the metric tensor, and the homothecal generator ξ is expressed as Putting Equations (28) and (29) together with (22) leads to and which can also be written as ζ(X(u, r)) = 0 (32) By solving Equations (30) and (31), it is possible to know the form of the metric variables β and m for all Y(u, r).

Proposed Solution to the Problem
When a problem in classical mechanics is very symmetric, it is required to simplify the partial differential equations to ordinary differential equations, where the dependent variables behave in essence as a single variable, which depends on the unstressed combination of the independent variables. Thus, when a problem has spherical symmetry, the independent variables are the radius at the center of symmetry r and the time t, and the dependent variables are a function of Equation (34) [16,17].
This method has been successfully employed in the nonrelativistic theory of selfgravitating fluids. This technique will be applied to Equations (30) and (31).
If new independent variables µ and r are defined, then X = X(ζ) and Y = Y(ζ) are solutions of the Equations (32) and (33) when the self-similar variable is expressed as Defined, Y = Λe 2β r and X = m r . The solution of Equations (32) and (33) is of the form and, making use of differential geometry for the continuity of I (first fundamental form), there are solutions for m and β, which are e 2β(u,r) = r r Σ l + 1 (37) and In the above equations, the values of k and l are constants, and the quantity is evaluated on the surface Σ, as indicated by the subscript, such that r Σ (µ) represents the radius of the distribution.

Equations on the Surface. Coupling Conditions
In general relativity, there are coupling conditions known as Darmois conditions [41]. The ultimate purpose of these conditions is to avoid the appearance of singular behavior of physical variables on the coupling surface [18,19]. According to Darmois [41], the necessary and sufficient condition for the coupling of two space-time regions on a hypersurface is that both the metrics induced on it (first fundamental form).
The inner region, characterized by the line element and the outer region, described by the Reissner-Nordstrom-Vaidya metric [16], are separated by the hypersurface r = Σ(µ). Therefore, the coupling of the inner and outer solutions of the field equations must be considered. Furthermore, to know the evolution of the physical variables inside the distribution, it is necessary to construct a system of ordinary differential equations for quantities evaluated on the surface. To obtain these equations, we proceed as follows:

•
The functions m and β are required to be continuous across the edge of the distribution.

•
The material velocity at the surface is evaluated. • The first field Equation (20) evaluated on the ∑ surface.
Once the above is done, two ordinary differential equations are obtained with two functions of time to be determined, corresponding to

•
The velocity of the edge of the distribution.
These are the surface radius r Σ , the distribution edge velocity, and the function evaluated on r = Σ. It is noticeable that the metric functions and, thus, the electric charge and physical variables are obtained from the time dependence of the variables at the surface m Σ y r Σ , for which the differential equations on the surface are required, which are obtained from Equations (40) and (41) .
where the point on the variable indicates time derivative. In addition, the radiation flux at the surface: The field Equation (22) evaluated on the surface gets to know the total mass of the distribution of matter where ξ = (l + 1) 2 and ψ = 2ξ + k(4 + 3l + k). This allows us to obtain Σ from Equations (22) and (23) and the surface equations.
Such that one has a single equation on the surface that is numerically integrated with the specification of the total mass m Σ and the radius of the distribution r Σ = R and the radius of the distribution rΣ = R. Once the dynamics on Σ have been obtained, the following procedure is followed to calculate the values of the physical variables:

•
Specifying the initial electric charge distribution allows for finding the physical variables at any radius (or any layer) of the sphere. • Integrate Equation (3) to advance the charge C in time.

•
We have the evolution of the physical variables for all time and all radius within the matter distribution.
Below is the routine in FORTRAN 77 for the solution of the surface equations with Runge-Kutta ( Figure 1).

Discussion of Results
A priori, to find solutions to electromagnetic systems with the presence of gravity (Einstein-Maxwell systems) as those proposed in the present work, there are no restrictions in the specification of which values the constants k, l, and the initial charge C can take. Only that the choice allows us to obtain physically acceptable behaviors, as illustrated in the results. Any choice of the above parameters must be consistent with −1 < Σ < 1 and Σ > 0 over the initial light cone (integration cone). At this point, we have a system of four nonlinear differential equations in partial derivatives, in which geometry (metric functions m and β with the material content of the system under study (pressure, density, electric charge, and radiation) and dynamics (velocity of matter in the comoving, ω) are mixed. All these variables correspond to the interior of the distribution. The time evolution of the surface of the distribution, how it collapses, and the speed at which it collapses must also be known.

Integration of the Conservation Equation
Once the integration of the equations on the studied surface has been achieved, it is necessary to proceed in the same way with the conservation equation to know the values of all the physical variables within the source. For this purpose, the spatial marker of comovility x = r/a. is used. Thus, the conservation equation can be written in the following way: which is defined as a wave equation that is integrated by using the numerical method of Lax-Friedrichs, with convergence conditions of Courant-Friedrichs-Levy (CFL). When treating the conservation equation in the described way and obtaining its solution, it must be remembered that this solution is conditioned by the evolution of the surface, so it was implemented as follows:

Discussion of Results
A priori, to find solutions to electromagnetic systems with the presence of gravity (Einstein-Maxwell systems) as those proposed in the present work, there are no restrictions in the specification of which values the constants k, l, and the initial charge C can take. Only that the choice allows us to obtain physically acceptable behaviors, as illustrated in the results. Any choice of the above parameters must be consistent with −1 < wΣ < 1 and mΣ > 0 over the initial light cone (integration cone).
At this point, we have a system of four nonlinear differential equations in partial derivatives, in which geometry (metric functions m and β with the material content of the system under study (pressure, density, electric charge, and radiation) and dynamics (velocity of matter in the comoving, ω) are mixed. All these variables correspond to the interior of the distribution. The time evolution of the surface of the distribution, how it collapses, and the speed at which it collapses must also be known.

Integration of the Conservation Equation
Once the integration of the equations on the studied surface has been achieved, it is necessary to proceed in the same way with the conservation equation to know the values of all the physical variables within the source. For this purpose, the spatial marker of comovility x = r/a. is used. Thus, the conservation equation can be written in the following way: which is defined as a wave equation that is integrated by using the numerical method of Lax-Friedrichs, with convergence conditions of Courant-Friedrichs-Levy (CFL). When treating the conservation equation in the described way and obtaining its solution, it must be remembered that this solution is conditioned by the evolution of the surface, so it was implemented as follows: In Equation (47), one has that the superscript n is associated with µ = nδµu, which is a hypersurface, and the subscript j is associated with the observer in comoving with x = nδx. Typically, δµ = 10 −2 with the CFL condition δµ = 2δx. The initial and boundary conditions need to be defined in order to proceed to integrate the conservation equation. For the case study of this paper, the boundary condition is given by Adding the initial condition results in the following: where the power p made it possible to investigate the level of effect of the initial conditions on the obtained results. The initial charge function is specified (its radial dependence) so that it satisfies the conservation equation. At this point, there are no restrictions to go inside the distribution and find the time evolution of the physical variables. For this, the specific charge function in relation to the radial coordinate r, at the time µ = 0, needs to be determined. Taking into account whether the time evolution does not spoil the physical requirements, the models can be obtained without difficulty. The sphere was divided into five layers, where the physical variables obtained from the field equations are evaluated. The collapse was imposed by imparting an initial velocity to the surface, which allows the evolution of the surface. The integration of the differential equations on the surface was done with the fourth-order Runge-Kutta method. In addition, the physical variables during the integration were kept bounded within physically acceptable values.
This is density greater than zero, pressure greater than density, and velocity bounded between −1 y 1, thus avoiding propagation velocities above the speed of light. This is why the CFL condition is used.
An important feature is the anisotropy factor p r − p t = 2μ, where it is evident that the anisotropy is induced by the charge. If the value of the charge is zero, of course, the anisotropy factor is zero. The radius of the distribution collapses inexorably. Unlike [2,13,18], the electric charge does not stop the collapse. Figure 2 shows the FORTRAN 77 routine to solve numerically the charge conservation equation. Figure 3 models the rate of change of the radius of the distribution over time, where it is observed that it collapses inexorably. In addition, the evolution of the physical variables is shown. Figure 4 shows the evolution when r/r Σ varies from 0.6 to 1.0, for k = 0.35; l = 0.5 initial charge and radius C = 0.292 and R (0) = 2.923, respectively. Figure 5 shows the evolution of radial pressure with time, for r/r Σ from 0.6 to 1.0, for the values k = 0.35; l = 0.5; C = 0.292; and R (0) = 2.923. The evolution is shown for different layers of the matter distribution. Figure 6 shows the radiation flux (multiplied by 102) for k = 0.35; l = 0.5 and initial conditions C = 0.292 and R(0) = 2.923. The evolution with time for different layers of the matter distribution is shown. Figure 7 shows the evolution of the charge function with time for constants k = 0.35; l = 0.5 and initial charge conditions and radius C = 0.292 and R(0) = 2.923.

Conclusions
An additional homothetic (self-similar) symmetry to space-time was assumed, and nonstatic solutions for the electromagnetic system with added gravity (Einstein-Maxwell systems) are obtained in a simple form. Coupling with the outer Reissner-Norstrom-Vaidya space-time and using it as a transport mechanism, the radiation flux freely escapes from the distribution.
It has been demonstrated, once again, as well as matter can be considered as an anisotropic fluid, when it is electrically charged. The approach is achieved by obtaining a dynamical model by means of the free escape approximation and self-similarity within the source. It is represented in the chosen example by manipulating only the initial conditions k, l, and CT and setting physical conditions: −1 < ω < 1, p < ρ y ρ > 0. A considerable set of initial conditions was selected for several runs of the programs, verifying the results

Conclusions
An additional homothetic (self-similar) symmetry to space-time was assumed, and nonstatic solutions for the electromagnetic system with added gravity (Einstein-Maxwell systems) are obtained in a simple form. Coupling with the outer Reissner-Norstrom-Vaidya space-time and using it as a transport mechanism, the radiation flux freely escapes from the distribution.
It has been demonstrated, once again, as well as matter can be considered as an anisotropic fluid, when it is electrically charged. The approach is achieved by obtaining a dynamical model by means of the free escape approximation and self-similarity within the source. It is represented in the chosen example by manipulating only the initial conditions k, l, and CT and setting physical conditions: −1 < ω < 1, p < ρ y ρ > 0. A considerable set of initial conditions was selected for several runs of the programs, verifying the results

Conclusions
An additional homothetic (self-similar) symmetry to space-time was assumed, and nonstatic solutions for the electromagnetic system with added gravity (Einstein-Maxwell systems) are obtained in a simple form. Coupling with the outer Reissner-Norstrom-Vaidya space-time and using it as a transport mechanism, the radiation flux freely escapes from the distribution.
It has been demonstrated, once again, as well as matter can be considered as an anisotropic fluid, when it is electrically charged. The approach is achieved by obtaining a dynamical model by means of the free escape approximation and self-similarity within the source. It is represented in the chosen example by manipulating only the initial conditions k, l, and CT and setting physical conditions: −1 < ω < 1, p < ρ y ρ > 0. A considerable set of initial conditions was selected for several runs of the programs, verifying the results shown. Several solutions were discarded, as they did not meet the above energy conditions.
More general solutions of the symmetry equations are needed to explore the generality of the results shown in this paper. Nevertheless, the anisotropic approach employed achieves results that are in agreement with those reported in other previously published work and which conclude that the final state is the steady-state [18].
An important feature is the anisotropy factor p r − p t = 2µ, where it is clear that the loading induces anisotropy. If the charge is zero, of course, the anisotropy factor is zero. This fact was used as a way to validate the code, because results for zero charges have been reported previously. Figure 3 shows the radius of the distribution, which collapses inexorably. Unlike [2,13,18,42,43], the electric charge does not stop the collapse.
In theory, while a black hole can possess a very large charge, the charge of the newly detected holes will be considered negligible. This assumption is based on the fact that a black hole is always surrounded by plasma. With the self-similar solutions proposed in this work, new models of black hole formation can be considered and also verify that the electric charge in the formation of a black hole can be significant. Models already reported establishing that, although the charge of a black hole is large, the charge/mass ratio is small and is not enough to significantly affect either the gravitational field of the matter distribution or the dynamics of the star or its collapse [44][45][46][47][48][49][50][51][52][53].