Kerr Black Holes within a Modified Theory of Gravity

The Kerr black hole is studied within a modified theory of gravity, which adds the effects of vacuum fluctuations near a black hole. These vacuum fluctuations are treated as a dark energy. A parameter is introduced to account for these fluctuations. It is zero for the standard theory and acquires a maximal value, just before there would be no event horizon. The existence of an event horizon not only depends on the value of this parameter, but also on the spin of the black hole. In addition, we study the existence of a light-ring. We also elaborate on the relation of the appearance and vanishing of the event horizon and light-ring to phase transitions.


Introduction
The General Relativity (GR) is one of the best tested theories [1]. Recently, the detection of gravitational waves [2] and the first optical resolution of the black hole in M87 [3][4][5][6][7][8] are all consistent with GR, i.e., one observes a shadow of the black hole. The size of it can be understood within GR and its radius is about 3 m 0 . Unfortunately, the resolution is only of 24 µas, which will smooth detailed structures, which are predicted by us (see the discussion in the main body of the text).
There are some concerns that, for strong gravitational fields, compared to the solar system, modifications to GR probably have to be applied. This did lead to several attempts to modify GR, which will be all resumed briefly in the next section. Here, we will concentrate on the so-called pseudo-complex General Relativity (pc-GR).
A complete review on the pc-GR, its motivation and relations to other approaches can be found in the book by Hess, Schäfer and Greiner [9]. The mathematics of the pc-GR is explained in the first and last chapter of the book, using differential geometry. Information on this can also be found in [10]. The pc-GR was first proposed in [11]. A general formulation, besides in [9], can also be found in [12,13], where, in [12], circular stable orbits were investigated and the last stable orbits, compared to GR. One important result is that, for small a, the last stable orbit in pc-GR is a little bit closer to the black hole than in GR. From a limiting value of a on, stable orbits up to the surface of the star exist. In this review, we will elaborate on the consequences of it. The theory predicts a different behavior of an accretion disc near the black hole, showing a dark ring followed by a bright one further in. Studies on this were also published in [14,15] and more recently in [16,17]. The main feature is the ring structure, mentioned above, which is quite independent from the mass of the black hole and the type of the accretion disc (if it is thin, thick, a torus, . . . ), save that the positions of the structures scale with the mass. The intensities of the light emission are larger in pc-GR than in GR (see also the discussion in the main body of this review).
In addition, neutron stars were studied in [18][19][20], in addition to [9]. Including the effects of the pc-GR, stars up to six solar masses were obtained in [18,20] and up to 200 solar masses in [19].

The Modified Theory: Pseudo-Complex General Relativity
At first, a historical overview on attempts to extend the GR is given, in addition to the reasons for it and why we decided for the pc-GR: Extensions of the GR have been discussed several times in the past: Einstein extended the metric to a complex one [30,31], in an attempt to unify GR and the Electro-Magnetism. He defined a complex metric where the real part is the standard metric of GR and the imaginary part is the electromagnetic tensor. The real part is symmetric while interchanging the indices and the imaginary part is anti-symmetric. This can be seen as follows: Comparing both sides leads to g µν = g µν and F µν = − F νµ .
Why a complex extension did not work will become obvious in a moment. The motivations of Born [32,33] were quite different. His concern is that, in Quantum Mechanics, the coordinates and momenta are treated on an equal footing. Canonical transformations of all kinds are allowed, which leave the commutation relations invariant, and one can even interchange coordinates and momenta. On the contrary, in GR, the coordinates play a singular role-in the length element square, only coordinates appear. He tried to recover the symmetry between the coordinates and momenta, leading to a modification of the length element, including momentum dependent terms, however, with the price of a mass dependence on a particle. The attempt by Born was retaken by Caianiello [34], who introduced in the length element squared an infinitesimal quadratic 4-velocity term, without the reference to a particle, implying an infinitesimal length scale, equivalent to a maximal acceleration: The l is a length parameter and is not subject to a Lorentz (or Poincaré) transformation. Thus, Lorentz invariance is guaranteed with an infinitesimal length scale in the model! Extracting an eigentime dτ 2 , taking into account that du µ dτ is an acceleration a µ and using that −η µν a µ a ν = a 0 a 0 − a i a i = −a 2 (η µν is the Minkowski metric), we obtain dω 2 = 1 − l 2 a 2 g µν , i.e., it corresponds to a new metric G µν = 1 − l 2 a 2 g µν . The acceleration is then limited by a ≤ 1 l 2 . We will see that these modifications are included in pc-GR in a particular limit. A more detailed resume on former intentions to extend the GR is given in [9].
In [35], all kinds of algebraic extensions of the coordinates were considered and the field equations for weak gravitational fields were obtained. It was shown that nearly all algebraic extensions contain solutions for tachyons or ghosts, which shows that nonphysical solutions appear. Only real and pseudo-complex (called in [35] hyper-complex) coordinates did not have this problem. This is the reason why we discuss only this particular extension: In pc-GR, the coordinates of GR are extended to these so-called pseudo-complex (pc) variables: (µ = 0, 1, 2, 3), with I 2 = 1, which justifies its name (though the name is not universal in the literature, where these variables are denoted as para-complex, etc.). The complex conjugate is defined as X * µ = x µ − Iy µ . In [34], the variable y µ is proportional to a minimal length scale, l, multiplied by the 4-velocity component u µ . It is important to stress that the parameter l is not a physical length and, thus, is not effected by a Lorentz transformation. Unfortunately, the relation of y µ to the 4-velocity is only correct in flat space, where the constraint used further below in (9) reduces to the standard dispersion relation g µν y µ dy ν = 0, with the solution y µ ∼ u µ . The minimal length parameter l is required as a factor by dimensional considerations. When the space is not flat, no easy solution has been found yet, but, in principle, can be found by solving the constraint (9).
An important nature of these variables is revealed when we change to the basis It implies that there are variables of the type vσ ± that have a zero norm. Thus, the variables form a ring and not a field. The components of σ ± are called zero divisor components.
The division into the zero divisor components can be done for any pc-function F(X) = F(X + )σ + + F(X − )σ − . Mathematical manipulations are independent from each other (see also [9]).
The fact that σ + σ − = 0 allows for formulating in each zero-divisor component a theory of General Relativity! In order to get a consistent theory, both zero-divisor components have to be connected. In [9,11], this is done introducing a modified variational principle, where the variation of the action is proportional to a "general zero", i.e., a function with a zero norm. However, this is not necessary: in [36] and the last chapter of [11], it is shown that a constraint can be implemented and a usual variational principle leads to modified Einstein equations. In what follows, we will resume the main steps.
The constraint requires that the pc length element is real. We used the definitions In addition, the representation of d 2 ω in the zero-divisor is indicated and in the original basis ((1, I)-basis). (7) coincides with the one of [34], when y µ is substituted by lu µ . In each component of the zero divisor basis, only Riemannian manifolds are considered, thus other types of manifolds are not included yet.
Setting the pseudo-imaginary part to zero leads to the constraint The action, without the constraint, is given by where R is the Riemann scalar. The last term in the action integral allows for introducing the cosmological constant in cosmological models, where α has to be constant in order not to violate the Lorentz symmetry. This, however, changes when a system with a uniquely defined center is considered, which has spherical (Schwarzschild) or axial (Kerr) symmetry. In these cases, the α is allowed to be a function in r, for the Schwarzschild solution, and a function in r and ϑ, for the Kerr solution.
The variation is independent in each zero-divisor basis. The result is the set of the modified Einstein equations (see [36]) with The energy-momentum tensor corresponds to the one of an asymmetric ideal fluid. The y µ is related to the appearance of a minimal scale, as explained in the case of the proposal given in [34]. Because these effects are practically impossible to measure due to their smallness, it suffices to restrict to the real part of the equations, which leads to the only one equation: The energy-momentum tensor acquires under this approximation the form The relations to the parameter in the action are (see [36]) where k µ are the components of a space-like vector, orthogonal to the 4-velocity.
Because we do not know the exact solution for y µ , we cannot derive the exact form of the energymomentum tensor and, thus, we are left to propose a phenomenological model. The energy-momentum tensor describes an ideal anisotropic fluid. That it has to be anisotropic is explained in [9]: First of all, as explained further below, the energy-momentum tensor on the right-hand side of the Einstein equations is assumed to be related to the contribution of vacuum fluctuations, which represents a dark energy. This dark energy can be calculated within the one-loop approximation of GR as described in [37] and done in [38]. The equation of state of a dark energy in the radial pressure is p Λ r = −ρ Λ , where p Λ r is its pressure and ρ Λ is the dark energy density. Within the Tolmann-Oppenheimer-Volkoff (TOV) equations, the radial derivative of the pressure is proportional to the sum of the radial pressure (which is equal to the tangential one in case of an isotropic fluid) and the density, i.e., it is zero. An anisotropic fluid is characterized by an addition term, proportional to the difference of the radial and tangential pressure (p Λ θ ). Without this term, the radial derivative of the pressure is zero and thus constant in r. Because of the proportionality of the energy density to the pressure, this also implies a constant energy density. Requiring that the density has to vanish at infinity results in a zero density, contrary to the requirement that the density is building up toward smaller distances. Only an anisotropic fluid can resolve this contradiction (see [9]).
In the absence of a quantized theory of gravity, it is impossible to deduce, for example, the radial dependence of the dark-energy density. Here, we rely on information from one-loop calculations in GR [37], which where performed in [38] for a static Schwarzschild back-ground metric. In [38], the energy density rises proportional to 1/ r 6 1 − 2m 0 r 2 toward the center, i.e., it is singular at the Schwarszschild radius. Clearly, the assumption of a static metric fails when the gravitational field increases too much. In this case, one has to include back-reaction effects, which is very difficult to do. This tells us that the fluctuations increase toward the center and they can become large. To avoid the singularity at the Schwarzschild radius, we propose a phenomenological model, where the dark-energy density is treated classically and behaves as where, in a first attempt, n = 3 was taken. It is strong enough in order not to contribute to the known observations within the solar system and other systems with not too strong gravitational fields [1]. There may be other dependencies with n > 3 and, in fact, in [27,28], it was shown that the fall-off of the dark-energy density has to be stronger. We will discuss other cases of the fall-off, ordered by a number n as in (16); however, the main conclusions and structural predictions of the theory remain similar.
In [12], the pc-Kerr solution was derived and here we list it for any n The a is the spin parameter in units of m 0 . The solution is identical to the standard Kerr solution, except for the term − B n (n−1)(n−2)r n−2 . The m 0 is the mass value observed at infinite distance and occasionally it will be abbreviated in the figures simply by m, clear from the context. The B n is given by For the equal sign, an event horizon is located at In this contribution, some examples are calculated with n = 3 or n = 4 and we will vary the parameter b n in B n = b n m n 0 from zero to a maximal value 64 27 for n = 3 and 81 8 for n = 4. For a larger value, there is no event horizon present anymore. We will also vary the spin-parameter a from zero to m 0 and show that, already for b n = 0, there is a region of large a values where this horizon disappears. For example, the question is, if there is or is not an event horizon already becomes relevant for tiny admixtures of vacuum fluctuations! This is similar, but not equal, to the so-called naked singularities: First, there is no singularity within the pc-GR and, second, what is exposed is the surface of a star. However, seeing the surface will be very difficult because the red-shift tends to infinity.
Another topic will be the calculation of the position of the light-ring. A light-ring is defined as the geodesic of a photon in a circular orbit (r = const). In GR, there is, for a = 0, only one light-ring at r = 3m 0 and the value is slightly lowered for larger a. In contrast, in pc-GR and, for values b n > 0, the structure becomes quite involved. We will investigate this property and relate it to phase transitions, using the Catastrophe Theory [29].
The application of pc-GR is not only limited to the region outside of a mass distribution, but also was applied to the interior of the star [18,19,[39][40][41][42][43]. The main problem is to propose a coupling of the dark energy to the mass distribution, which, due to not knowing it from first principles, is always charged with phenomenological assumptions. In [18], a linear relation was assumed, which, however, leads to an upper limit of six solar masses for a star. The reason is that, near the surface, the repulsion is strong enough that the star sheds its upper parts of the surface. In [19], this was resolved partially, applying a calculation within the one-loop approximation of quantized gravity [37] and using the monopole approximation. As a result, the distribution of the dark energy turns out to fall off stronger near the surface. With that, stars with up to 200 solar masses were found to be stable. This already indicated that any mass of a star can probably be stabilized, i.e., even the large masses in the center of galaxies of billions of solar masses are stable and rather dark stars than black holes. There are different, alternative, approaches, e.g., in [39][40][41][42][43] compact and dense objects were investigated within the pc-GR and maximal masses were also deduced.

Circular Orbital Motion
The infinitesimal length square element dω 2 in terms of the metric components is given by with the metric components listed in (17). Circular orbital motion is defined byṙ = dr dτ = 0, where dτ is either the infinitesimal difference of the eigentime or the length element dω. The particle moves on a geodesic path and the deduced frequency is given by [17] In [46], the orbital motion frequency of a particle in a prograde circular orbit (assumed here) was deduced, using the method of an effective potential. In addition, the motion of a photon in a circular orbit is deduced, which satisfies, in addition, dω 2 = 0.
The approach presented here is different from [47], where the path of a particle following a gravitational collapse of a star is investigated, including a semi-quantum mechanical description of the particle. Although it is very interesting, it would lead here to a too extensive discussion, out of the scope of this article, also requiring adding a complete new description to the problem. It should, however, be useful to discuss a collapsing star within pc-GR. One first step is to change the metric used in [47] to the modified metric in pc-GR.
The orbital frequency of a particle in its prograde orbit is shown in Figure 1, for a = 0.9m 0 and n = 3, n = 4, where the b n values are 64 27 for n = 3 and 81 8 for n = 4. The upper curve is the result of GR while the two lower ones are for n = 3 and n = 4 in pc-GR. The maximum is at independent from a. In case b = b n max , there is a zero of the orbital frequency, coinciding with a horizon (19), and it is also independent on a. r 0 is at the estimated position of the star's surface. For 0 ≤ b n < b n max , the lower curves slowly merge into the upper one for b → 0. The stable orbits are defined in terms of a saddle point of the effective potential, see [9,14], and they are plotted in Figure 2. The left panel shows the region of stable orbits in GR and pc-GR for n = 3, while the right shows the region of stable orbits for GR and pc-GR for n = 4 (see also the figure caption for explanation). When n increases, for small a, both curves, GR and pc-GR, approach each other. However, the curve in pc-GR starts at a = 0 at smaller distances to the black hole than GR. The light emission from a disc is calculated using the 2014 version of the GYOTO routines [48], which uses the raytracing method and also allows for changing the metric input. Because the last stable orbit reaches further in, more gravitational energy is released and light emission increases. This can be verified in Figure 3, where simulations of thin accretion discs are presented for a = 0.9 m 0 , using the model of [49]. The left panel is pc-GR for a = 0, being similar in structure to GR but with a greater intensity, and the right panel for a = 0.9 m 0 , which now shows a dark ring followed further in by a bright ring. In Figure 4, the simulations are shown for the same type of disc, but for different n. The left figure refers to n = 3 and the right one to n = 4. When n is increased, the position of the dark and bright rings are shifted to larger r values. In fact, the parameter n can be determined through (21) and measuring the position of the dark ring, which gives the maximum of the frequency.   Figure 3. The left panel shows the simulation for pc-GR and n = 3 (the same as the right panel in Figure 3), while the right panel is a simulation for n = 4. As noted, the position of the dark and bright rings are shifted slightly to larger radial distances. The figures were obtained using the open accessible 2014 version of the GYOTO routines [48]. For n = 4, the modified C++ rountines can be retrieved in [45]. This behavior is understood by inspecting Figure 1: at the maximum of the orbital frequency, neighboring orbitals are very similar and friction is low, such that the disc gets less excited and a dark ring forms. Above and below the maximum, the change in orbital frequency for neighboring orbitals increases and friction is high, which results in a larger excitation and the light emission increases. In this sample calculation, the orientation of the accretion disc relative to the observer is 80 • . However, these studies require a resolution of at least 0.5 µas, not yet reached by the EHT. In the upper row of Figure 5, simulations of an accretion disc are shown, with a large resolution (left) and a low resolution (right). The inclination of the disc with respect to the observer is now 70 • and a = 0.6. Below these figures, the EHT result is copied from [3] and rotated by 90 • .
In Figure 6, we compare the simulation of an accretion disc in GR (left panel) with the one in pc-GR (right panel), assuming a resolution of 20 µas, as reported by the EHT-collaboration. As can be seen, the resolution is too low for detecting significant structural differences. This probably also affects other extensions of GR, i.e., the EHT results make it difficult to discriminate alternative theories. The flux of light emitted by a thin disc, discussed in the above simulations, is determined as explained in [14,49] (please consult these references for more detailed information). The formula used is where ω |r is the derivative in r of the circular orbital frequency, L z is the orbital angular momentum, with L z|r its derivative in r and E is the orbital energy. The lower limit rms = r ω max is the position of the maximum of the orbital frequency. As shown in [14], the integral yields only positive values. The origin of the flux is in the friction between neighboring orbitals and it is distributed from orbitals of larger frequencies to orbitals with lower frequencies. The intensity in pc-GR is always larger. The peak shifts to the left (lower distances) as a increases. For a = 0.9 m 0 , the peak is around r = 0.4 m 0 , while the peak of the pc-GR curve is always at lower r. The intensities are obtained using (23). The relation between L z and E can be retrieved from [14].

Event Horizons
In this section, we discuss the dependence on the existence of an event horizon as a function on the parameters b n and a. In addition, we take into account the range of b n from 0 to b n max , for n = 4. Furthermore, we relate the disappearance of the event horizon with a phase transition, whose physical significance is not yet clear to us and will be investigated in future. The relation to phase transitions is completely new.
The surface of infinite redshift is defined by the zero of g 00 in Equation (17). It depends on the azimuth angle ϑ. As in GR, this defines an outer surface that engulfs the event horizon.
To obtain an equation, which determines the surface of the event horizon, one has to look for time independent axially symmetric surfaces with a null-norm [9,51]. A surface is described by u(r, ϑ) = 0. The normal vector is given by satisfying n µ n µ = 1. The steps to find the final equation is identical to what is presented in [9,51]. We define new variables y = r m 0 and a = a m 0 and the equation for the event horizon is now which is equivalent to where we redefined the r to y and a to a, as given above. At the pole (ϑ = 0), the infinite redshift surface joins the even-horizon, while the infinite redshift surface is further out at the equator (ϑ = π 2 ). The space in between the infinite red-shift surface and the event horizon is called the ergosphere [51].
In order to discuss phase transitions, we will follow a similar path as applied in [52]: let us integrate (26), which leads to the auxiliary "potential" This potential has, to our knowledge, no physical meaning, but the discussion on phase transitions is best explained with a potential [29]. The extremal points of the potential are determined via Equation (26). We also allow values of y < 3 2 , where, for n = 4, the surface of the star is supposed to be. Though these values do not have a physical meaning, it is better for the illustration of the phase transitions.
As noted, at b n = 0, we have the known solutions of y versus a: For a = 0, the only physical solution is y = 2, which approaches y = 1 for maximal rotation a = 1 m 0 . The second positive solution is always lower than 1 and, thus, in GR always below the event horizon at a = 1 m 0 . For non-zero b n , a novel feature appears, namely that, from a certain a on, there is no event horizon! For maximal b n , there is only one event horizon at a = 0 and, for all a > 0, the event horizon disappears! Thus, even for b n smaller than the maximal value, an event horizon does not always exist, even for tiny contributions of vacuum fluctuations! In order to construct the surface of the critical points, i.e., where the derivative of the potential is zero, we follow for the case n = 4 closely the method of the Catastrophe Theory [29] as presented in [53]. First, we construct in the three-dimensional space with the coordinates (y, b n , a 2 ) the surface of critical points where (26) is satisfied. The critical b n (b cr ) is obtained resolving (26) and renaming y cr by λ 1 , Further below, we will also need the derivative of b cr with respect to λ 1 , which is ∂b cr ∂λ 1 = −12a 2 λ 1 + 36λ 2 1 − 24λ 3 1 . Now, we consider the projection map of the two-dimensional critical surface onto the two-dimensional sub-space (b n , a 2 ): This defines a singular mapping, if the Jacobian of the transformation satisfies, redefining a 2 as λ 2 , This gives the condition for a 2 = a 2 cr : Substituting this into (28) for a 2 and taking the square root of a 2 (which is always positive), we construct the parametric curve which defines the separatrix. The separatrix is depicted in Figure 8 in the plane (b n , a), together with the surface of critical points. It is clearly seen that the separatrix is the projection of the points where the slope of the critical surface is infinite. The left panel shows the case n = 3 and the right one for n = 4. This figure illustrates that for a given b n the range of a is finite where an event horizon exists. i.e., as soon as there is a little accumulation of dark energy near a black hole, it loses its event horizon from a given a on! In Figure 9, we plot the position of extrema of the potential, one for a = 0 and the other for a = 0.5 m 0 . The upper part in each pair of curves corresponds to a maximum, while the lower one corresponds to the deformed minimum. Both join at a certain value of b n , after which no extrema exists. This point is just at the position of the separatrix (denoted in the figure by circles). The larger a is, the point of the separatrix (where the two curves meet) moves further to the left.
In Figure 10, we show a sequence of potential, all for a = 0.5 m 0 , as a function in b n . While for b n = 2 the global minimum of the potential is at negative values, i.e., it is lower than U at y = r m 0 = 0. For b n = 3.3, the deformed minimum is approximately at the same height and, at b n = 4, it is at positive values. At b n = 7.5 it is close to the point when maximum joins the minimum, i.e., forming a saddle point. This type of behavior, the relative position of minima and the change of it, is typical for a phase transition.
Finally, we will discuss the order of phase transition when reaching the separatrix: For that, one has to calculate the derivatives of the potential as a function of the parameters (b cr , a cr , y cr ), where y cr is the critical value for the distance y = r m 0 . The best way to do so is to expand the function in a Taylor series around the critical point, also called the germ of the critical point. Expanding solely in y, we obtain up to the 5th order the expansion 1 5 (y − y cr ) 5 + (y − y cr ) 4 (− 1 2 + y cr ) + (y − y cr ) 3 ( a 2 cr 3 − 2y cr + 2y 2 cr ), +(y − y cr ) 2 (a 2 y cr − 3y 2 cr + 2y 3 cr ) + (y − y cr )( b 6 + a 2 y 2 cr − 2y 3 cr + y 4 cr ).
When the critical values for a cr and b cr are substituted, as given above and renaming y cr by λ 1 , the terms of lower order in (y − y cr ) p (p ≤ 2) vanish. The factor of the first non-vanishing term, which is proportional to (y − λ 1 ) 3 ,with y cr = λ 1 , is given by Because the first non-vanishing term in the expansion is proportional to (y − y cr ) 3 , we can conclude that it corresponds to a third order phase transition. The meaning of this is not clear yet but will be investigated in future.
That we find for the critical surface, for a given a and b n , two solutions at positive y can be understood readily: for that, let us restrict to a = 0 (Schwarzschild). Then, the zero of g 00 = 1 − 2m 0 r + B n (n−1)(n−2)r n coincides with the position of the event horizon. Due to the positive third term in g 00 , there are two-solutions. For B n = 0, there is only one at r = 2 m 0 . When B n starts to be different from zero, for very small B n , the lower solution is at small positive r, while the upper one starts to decrease. Increasing B n , little by little, the two solutions approach each other until they meet at the point of the separatrix.

Light-Rings
As mentioned above, the condition for a light-ring in a circular orbital motion is that it has to be geodesic and simultaneously dω 2 = 0. The geodesic equation leads to the prograde frequency [12] ω = − g 03 g 33 + g 03 2 − g 00 g 33 which, in turn, leads to Equation (21). The prime denotes the derivative with respect to r. For the light-cone propagation (dω 2 = 0), the frequency for a circular orbit is [12] ω = − g 03 g 33 + (g 03 ) 2 − g 00 g 33 It has the same form as in (35), except that no derivatives of the metric components appear.
The two frequencies deduced in (35) and (36) have to be set equal, which leads to the following equation, deduced in [27,28] (we changed the variable r into y = r m 0 and defined a = a m ) √ ∆ y 3 − a 2 F + a 2y 2 m(y) + y 2 + a 2 F − y yFg 22 with ∆ = y 2 + a 2 − 2ym(y), F(y) = m(y) − m (y)y and m(y) = 1 − b 4y 2 . Resolving (37) for r gives the position of the (or more) light-ring(s). This equation and, thus, the position of the light-ring depends on the azimuth angle ϑ via g 22 . The expression is the same as in [27,28].
Here, we will discuss only the motion in the orbital plane with ϑ = π 2 and n = 4 (the other cases have the same structures). The equation for the light-ring, using the above described path, is given by We proceed in the same manner as explained in the last sub-section on the event horizon, constructing the surface of critical points, given by the solution of (38) that corresponds to the light-ring. A light-ring represents a close contact of neighboring light trajectories, and, thus, an envelope where light focusing takes place (a caustic). Caustics are extremely important in wave theory [54]. The light intensity grows to infinity in a caustic, limited by diffraction, and the light-wave field structure is dominated by these caustics singularities. In this way, we proceed again by means of Catastrophe Theory in order to study the regions of stability and its structure of the light-ring. As it shall be demonstrated, within pc-GR theory, there exists a further singularity of light rays, of a greater order: a double-light-ring, predicted within a limited region of the parameter space as a two-fold critical surface, experiences a further coalescence establishing an even brighter caustic of caustics, and signaling a second stability separatrix in the parameters' space that limits the possibility to any further existence of light-ring singularities at the outer region of parameter values.
The numerical solution of the critical surface is depicted in Figure 11, together with the critical surface for the event horizon, identified as the inner surface. The y-variable starts at 3 2 , where the surface of the star is assumed to be (n = 4). As can be noted, increasing b n , there is an increasing range of large a-values, where there is no light-ring anymore. Reaching b 4 = 81 8 , the light-ring ceases to exist from about a = 0.35 m 0 on. Even for a tiny contribution of the vacuum fluctuations, there is no light-ring near a = 1 m 0 ! This result puts some doubt on the assumption that a light ring still exists [55], which is based on the use of pure GR and not including the fact that vacuum fluctuations are building up. An interesting feature can be observed for large b n values: the light ring exhibits two real solutions. Figure 11. The critical surfaces for the horizon (inner surface) and the light-ring (outer surface). The critical surface of the light-ring follows the one of the horizon, only further out. Note the region where the light-ring does not exist. The b value is equal to b n . For the construction of the light ring surface, (38) was used. The corresponding MATHEMATICA routine can can be retrieved from [13].

Conclusions
We presented a report on the current status of the pseudo-complec General Relativity (pc-GR). This theory requires that, around a mass, dark energy has to accumulate, for which we constructed a phenomenological ansatz, i.e., that it falls off proportional to 1 r n+2 . This is multiplied by a parameter B n = b n m n that describes the coupling of the dark energy to the central mass. In most applications, the parameter b n is chosen such that there is no event horizon for any a-value. For this case, there is a lower value for b n , which is 64 27 for n = 3 and 81 8 for n = 4. In this contribution, we also investigated the range of b n -values from zero on, i.e., we considered small vacuum fluctuation too, which by force is there [37,38]. We found that, even a tiny amount of vacuum fluctuations erase the event horizon near a = 1 m 0 and also the light-ring ceases to exit for small values of b n and a near 1 m 0 , which is a new interesting feature. This has nothing to do with the so-called naked singularities because there is no singularity in pc-GR and what is exposed most is the surface of a star. However, the red-shift tends to infinity and the surface can not be seen.
One of the most important predictions is that, when an accretion disc is present and for a rotation large enough, described by the parameter a, it shows a dark ring followed further in by a bright one. This is an effect of the dependence of a particle in a circular orbit, whose frequency shows a maximum and falls off for smaller r. At the maximum, neighboring orbitals have approximately the same orbital frequency, thus friction is low and no light emission is produced. Below and above the position of the maximum in the orbital frequency, it shows a change for neighboring orbital and friction is large, resulting in a stronger light emission. We also compared the GR simulation with the pc-GR one, taking into account the low resolution of 20 µas, similar to the EHT observation. From this we can conclude that there is, for now, little possibility to see any structural differences. Only an increase in the resolution by a factor of 5 can probably show differences between GR and pc-GR. Until such a resolution can be obtained, a couple of decades have to pass because it will imply to put radio-telescopes in the orbit of the moon and beyond. Nevertheless, the pc-GR can be tested with these telescopes in future. We are currently in contact with members of the EHT on how to obtain a Fourier transform from our pictures, which is what the EHT observes. We plan to look for details in the intensity distribution and its Fourier transform. Unfortunately, we cannot present any results on this yet.
This phenomenon of rings depends on the rotational parameter a: For low a, the last stable orbit follows the one of GR, only further in. The result is a similar pattern of light emission in pc-GR, but, at a higher value, because the last stable orbit reaches further in and, thus, more energy is released which is distributed within the disc. Note that, for a = 1 m 0 also in GR, no light ring can be observed because the inner edge of the accretion disk is at r = m 0 , while the light ring is at larger r. Thus, the light ring is hidden by the accretion disk.
As a new contribution, we showed that the disappearance of the event horizon can be related to a phase transition, which is of the order of 3 when we reach the points where the surface of critical points has a singularity in its derivative. What consequences it may have is not clear to us and we are investigating it further.
Finally, we discussed the light-ring in the orbital plane (θ = π 2 ). We also found that, at the moment, some dark energy is accumulated around a black hole, and there is a region where the light-ring ceases to exist from a certain a-value on! The range increases with b n . We found that the critical surface of the light-ring engulfs the one of the event horizon.
Thus, adding even a small amount of dark energy around a large mass, leads for a certain a-value on to the disappearance of the event horizon and also eliminates the existence of a light-ring! One has to not assume a sufficient large value of b n in order for no event horizon or light-ring to exist, though the limiting value of b n = b n max serves that, for any a, no horizon exists.
Though a fluid identified with a dark energy is introduced, this fluid is treated classically, as all of the approaches presented in this review. We hope that the results may shed some light on how to include quantum mechanical effects directly. Often, classical treatments can shed some light on the quantum mechanical extension of the theory. The connection to a quantum mechanical description may be the reproduction of the dark energy behavior near a large mass, requiring first the observational confirmation of the structures predicted by the pc-GR.