{'fraction_non_alphanumeric': 0.037999127888868126, 'fraction_numerical': 0.008783404971033451, 'mean_word_length': 4.299933971607792, 'pattern_counts': {'":': 0, '<': 18, '<?xml version=': 0, '>': 12, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Charged Vaidya Solution Satisfies Weak Energy Condition

The external matter stress-tensor supporting charged Vaidya solution appears to violate weak energy condition in certain region of the spacetime. Motivated by this, a new interpretation of charged Vaidya solution was proposed by Ori [1] in which the energy condition continues to be satisfied. In this construction, one glues an outgoing Vaidya solution to the original ingoing Vaidya solution provided the surface where the external stress-tensor vanishes is spacelike. We revisit this study and extend it to higher-dimensions, to AdS settings, and to higher-derivative f(R) theories. In asymptotically flat space context, we explore in detail the case when the mass function m(v) is proportional to the charge function q(v). When the proportionality constant \nu = q(v)/m(v) lies in between zero and one, we show that the surface where the external stress-tensor vanishes is spacelike and lies in between the inner and outer apparent horizons.


Introduction
Evolution of charged null shells in spherically symmetric settings has been a source of confusion in the past [2,3,4,5]. Despite the simple form of the charged Vaidya solution and its seemingly simple interpretation, it suffers from a rather serious conceptual difficulty: it seems to violate the weak (as well as null) energy condition in certain region of the spacetime. In four-dimensional asymptotically flat context the region where weak energy condition is violated is, where m(v) and q(v) are respectively the mass and charge functions of the ingoing charged Vaidya solution, and v is the advanced null coordinate. In modern literature this feature has been argued to be related to the possible violation of strong sub-additivity of holographic entanglement entropy [6]. For the ingoing charged Vaidya solution, the spherical shells of "charged photons" all start at infinity and flow inwards. The key physical question is therefore whether these shells enter the region (1.1) where the energy conditions can be violated. Ori [1] in 1991 did a systematic and thorough analysis of this question. In this paper we revisit that study and extend it in a number of ways.
The key confusion in the earlier literature was the assumption that the Lorentz force does not matter and the charged null fluid flows along the radial null geodesics. In the Reissner-Nordström geometry the radial null geodesics go all the way from r = ∞ to r = 0. Assuming that the Lorentz force is irrelevant, earlier authors concluded that nothing can prevent charged null shells from entering the region where energy conditions are violated. Ori, on the other hand, showed that the assumption that the Lorentz force can be ignored is incorrect (in fact inconsistent, as we review below). By properly taking it into account he showed that the weak energy condition continues to be satisfied.
The most important outcome of this rather intricate analysis was as follows: locally the Lorentz orbits are tangent to one of the two (ingoing or outgoing) radial null geodesics, however, at r = r c the tangent direction flips. The null vector k a along the Lorentz orbit vanishes at r = r c . The physical continuation of the initially ingoing orbit at r = r c is to switch to the outgoing one.
Intuitive explanation of this result is rather straightforward. Charged particle moving in charged background looses its kinetic energy due to electromagnetic repulsion at it falls inwards. At r = r c its kinetic energy is reduced to zero. At this point the particle starts moving outwards.
Precisely the same intuitive picture applies to the charged Vaidya spacetime. At r = r c the charged null shells switch from ingoing to outgoing orbits, and there is no violation of the energy conditions. In this paper, after providing a concise summary of Ori's work, we extend the analysis in the following ways: 1. We adapt Ori's arguments to the AdS and higher-dimensional settings. This extension brings in some minor technical changes in the details; the physical picture remains unchanged.
2. We adapt the same set of arguments to f (R) theories coupled to a Maxwell field. Again this extension work straightforwardly. It suggests that the picture suggested by Ori is valid in other settings as well, e.g., Einstein-Gauss-Bonnet gravity coupled to a Maxwell field.
3. In the asymptotically flat context, we analyse in some detail the case when the mass function m(v) is proportional to the charge function q(v). For the physically interesting case where the proportionality constant ν = q(v)/m(v) is in between zero and one, we show that the surface where the external stress-tensor vanishes is spacelike and lies in between the inner and the outer apparent horizons. Since this surface is spacelike, to its future is an outgoing Vaidya solution. We discuss this gluing construction in some detail.
The rest of the paper is organised as follows. In section 2 we start with a review of the charged Vaidya solution and its standard interpretation. We review the fact that the interpretation of charged null fluid as made up of discrete charged massless particles is inconsistent with the geodesic motion assumption. We highlight how the null Lorentz force equation appears from such an analysis.
In section 3 we analyse radial orbits of massless charged particles in higher-dimensional AdS Reissner-Nordström background. In this analysis we focus our attention on the behaviour of orbits near the vanishing point k a = 0. This probe analysis already shows how initially ingoing orbits switch to outgoing ones. We also analyse the behaviour of orbits near a generic vanishing point.
In section 4 we return to the interpretation of the charged Vaidya solution and apply the results of the former sections. The evolution beyond the hypersurface of the vanishing points is discussed. f (R)-Maxwell set-up is discussed in section 5. We end with a brief summary of our results and some speculative discussion in section 6. Power law profiles for m(v) and q(v) functions in asymptotically flat context are further discussed in appendix A.

Preliminaries: Charged Vaidya solution and the Lorentz force equation
In section 2.1 we start with a brief review of the charged Vaidya solution and its standard interpretation. Then, following Ori [1], we show in section 2.2 that this interpretation is inconsistent with the physical picture that the shell is made up of discrete charged null particles. Consistency of the equations of motion requires us to replace the geodesic equation used in the standard interpretation with the null Lorentz force equation.

Charged Vaidya solution
Let us consider the charged Vaidya solution for the following Lagrangian in d spacetime dimensions, where L matter is the Lagrangian for the external matter source. The equations of motion are and T

(m)
ab is the stress-tensor of the external matter source. The cosmological constant parameter Λ is related to a characteristic length scale l as Λ = − (d−1)(d−2) 2l 2 . The length scale l is a more convenient quantity to work with. We take the cosmological constant term to be negative; we are only interested in asymptotically anti-de-Sitter and flat (obtained by sending l → ∞) settings. The ingoing charged Vaidya metric is described by the line element where the metric function f (r, v) is and dΩ 2 d−2 is the line element on the unit round (d − 2) sphere. The Maxwell potential A a supporting the solution is, The external matter stress tensor needed to support this solution is In particular, it is of the form where the vector k a is (ingoing) null vector tangent to the v = constant lines. Normalisations of ρ and the null vector k a are not fixed by the form (2.8). Further physical inputs are required to fix these normalisations. In the standard interpretation, which has been a source of confusion in the past, the choice is made. With this choice k a is the null vector generating the ingoing null geodesics. Therefore, the charged null fluid is thought of as flowing along these null geodesics. Under this assumption, charged null particles starting from r = ∞ run all the way to r = 0 along v = constant lines.
Null energy condition: Generators of radial null orbits l a = (l v , l r , 0, · · · , 0) that intersect the null fluid satisfies f (r, v)l v = 2l r . The local energy density K = T ab l a l b "experienced" by these trajectories intersecting the charged null fluid simplifies to, Clearly, for K(v, r) becomes negative; the null energy condition is violated.
Weak energy condition: In a similar manner the energy density E = T ab u a u b measured by a timelike observer with normalised velocity u a = (u v , u r , 0, · · · , 0) simplifies to, (2.12) The first term corresponds to the "Maxwellian" energy density due to the presence of electric charge. The second term is the energy density of the particles constituting the collapsing charged null shell. Once again, in the region (2.11) the second term becomes negative. Moreover, for sufficiently large u r this term dominates, and the total E(v, r) becomes negative; weak energy condition is violated.

Lorentz force equation
Now we demonstrate that the geodesic motion assumption is inconsistent. Above we distinguished between the matter and electromagnetic stress tensor. Let us look at the equation arising from Bianchi identity, i.e., the conservation of the total stress-tensor. We get This expression upon using the Maxwell equation We adopt the hydrodynamic point of view and think of the charged null fluid as a stream of discrete charged massless particles. We associate a conserved null current associated with the flow a := ρ e k a is thought of as conserved electric current; ρ e being the electric charge density. The ratio e := J (e) a /J a = ρ e /ρ is the charge per particle. The conservation of e follows, and it reflects the fact that each particle carries its own electric charge. Using these inputs, is the key equation for the consideration of this paper. It is not the geodesic equation. It is the massless analog of the Lorentz force equation for a massive particle; for simplicity we continue to call it the Lorentz force equation. We call the parameter λ along the orbits satisfying (2.18) the Lorentz parameter. In general, the Lorentz parameter is different from the affine parameter, however, when F b c k c = 0, it becomes identical to the affine parameter. In that case the Lorentz orbit becomes a null geodesic.

Probe computation
To understand the full implications of the Lorentz force equation (2.18), it is instructive to start with a probe computation. Due to the electromagnetic repulsion the magnitude of k a for infalling "charged photons" gradually decreases and at some spacetime point k a vanishes. Our focus in this analysis is on the behaviour of orbits near the vanishing point k a = 0. In section 3.1 we regard the higher-dimensional AdS Reissner-Nordström metric and the associated electromagnetic field as a given background, and study the Lorentz orbits in this background. This probe analysis already shows how initially ingoing orbits switch to outgoing ones. We analyse the behaviour of orbits near a generic vanishing point in higher-dimensions in section 3.2. Other studies of the motion of charged particles in the field of (rotating) charged black holes include [7,8].

Vanishing point in the spherically symmetric context
The AdS Reissner-Nordström background in Schwarzschild coordinates in d dimensions is where the metric function f is The only non-vanishing components of field strength tensor are To solve the Lorentz force equation (2.18) for radial trajectories in this background we take k a of the form k a = (k t (r), k r (r), 0, . . . , 0). With this ansatz equations can be readily solved. For ingoing orbits we obtain The parameter κ is interpreted as the energy of the infalling particle measured by a static observer at infinity. The vanishing point is located at We can parametrise such an orbit with a real parameter λ such that k a = dx a dλ . We find by direct integration, where c is an integration constant. It is fixed by the initial conditions; its numerical value is not important for the arguments below. For simplicity we set c = 0. The function F is the hypergeometric function. For d = 4 as r approaches r c the Lorentz parameter λ diverges to +∞. For d > 4 at r = r c , the parameter λ becomes which also diverges to +∞. This can be easily seen by the identity At the vanishing point the momentum vector k a vanishes. How to continue the orbit beyond the vanishing point? To answer this question further physical inputs are required. Mathematically there are two possibilities: 1. Non-bouncing or ingoing continuation: This corresponds to k r = k t . To begin with k t < 0, see equation (3.5). Across the vanishing point k t (and hence k r ) changes sign. Thus, for this case, preserves its sign as the orbit crosses the vanishing point. As a result the orbit continues along the same (ingoing) null direction. That is why we call it the non-bouncing or the ingoing continuation.
It is instructive to analyse the non-bouncing continuation in two separate cases depending on whether the metric function f is positive or negative at the vanishing point r = r c .
(a) f > 0: In this case dr dt < 0 to begin with at r = r c . It continues to remain negative beyond the vanishing point. However, the null vector k a changes sign -from future directed it becomes past directed -so the Killing energy carried by the infalling particles changes sign.
(b) f < 0: In this case dr dt > 0 to begin with at r = r c . That is to say that for an ingoing particle in between the inner and the outer horizon, both r and t decrease. Beyond the vanishing point dr dt does not change sign. However, once again, the null vector k a changes sign so the Killing energy carried by the infalling particles changes sign.
Since the null vector k a changes from future directed to past directed, this is our first hint that this continuation is not the physical one.
2. Bouncing or outgoing continuation: This corresponds to taking k r = −k t beyond the vanishing point. For this case dr dt changes sign. The Lorentz orbit switches from ingoing to outgoing. For this reason we call this the bouncing or the outgoing continuation.
Once again it is instructive to analyse the bouncing continuation in two separate cases depending on whether the metric function f is positive or negative at the vanishing point r = r c .
(a) f > 0: To begin with k r = k t < 0, and therefore dr dt = −f < 0, i.e., r decreases with time. Beyond the vanishing point, Therefore, r starts to increase with time. The orbit switches from the ingoing to outgoing. From equation (3.5) we note that k t continues to be negative beyond the vanishing point. Hence, k a continues to be future directed null vector; the Killing energy carried by the particle continues to be positive.
(b) f < 0: In this situation the particle is in between the inner and the outer horizon.
For future directed causal curves in this region the coordinate r always decreases. To began with dr dt = −f is positive. Beyond the vanishing point dr dt is negative dr dt = f < 0. To begin with k t < 0. From equation (3.5) we note that beyond the vanishing point, in the present case, k t changes sign, as the coordinate r being timelike continues to decrease beyond r c . However, k r does not change sign. As a result k a remains future directed.
Since the null vector k a continues to be future directed, the bouncing continuation seems to a more physical one.
To see that the bouncing continuation is the physical continuation we introduce a small conserved angular momentum k φ = L > 0. The introduction of the parameter L prevents the vanishing of the vector k a and therefore removes the above mentioned mathematical ambiguity. In the limit L → 0 these orbits converge to bouncing radial orbits, thus establishing that the bouncing continuation is indeed the physical one [1].
Let us take k a = (k t , k r , 0, . . . , 0, L), and consider Lorentz null orbits in an appropriate equatorial plane where g φφ = r 2 . The null condition gives The Lorentz force equation (2.18) can also be written as where Γ a bc is the Christoffel symbol. Once again the k t equation can be readily integrated with solution exactly the same expression as (3.5). The k r equation can be written as We are interested in the location of the turning point in the radial coordinate, i.e., the point where k r vanishes. In that case equation (3.15) simplifies to (3.16) Substituting k r = 0 in (3.13) gives Once again it is useful to analyse f > 0 and f < 0 cases for the location of the turning point k r = 0 separately.
1. f > 0: Substituting (3.5) in (3.17) we note that the turning point is located at 1 the solution of the equation For small values of L, we immediately see that the solution of this equation is at r > r c . Moreover, for small L, equation (3.16) becomes Recall that initially kt < 0.
Therefore, the particle bounces with positive k r . Since the turning point is at r > r c , k t does not vanish at the turning point. It continues to maintain its sign beyond the turning point; k a remains future directed.
Taking L → 0 limit in this analysis we conclude that the bouncing continuation of the radial null orbits is the physical one.
2. f < 0: In this case the particle is in between the inner and outer horizon. For future directed causal curves in this region the coordinate r always decreases. Therefore r = r c is reached with negative k r . The null condition (3.13) becomes From equation (3.5) we note that k t vanishes at r = r c and changes sign as r continues to decrease. From (3.20) we see that k r never vanishes. It continues to be negative and hence future directed.
Once again, taking L → 0 limit in this analysis we conclude that the bouncing continuation of the radial null orbits is the physical one.
To summarise: the above probe analysis allows us to conclude that the initially ingoing radial null Lorentz orbits switch to the outgoing ones.

Properties of the generic vanishing point
In this section we analyse the local behaviour of null Lorentz orbits near a generic vanishing point. We perform this study in an infinitesimal small neighbourhood of the vanishing point; spacetime curvature and variations of the electromagnetic field are not important [1]. We can obtain relevant properties of vanishing points and behaviour of orbits near them by working in flat spacetime with constant electromagnetic field. Equation (3.14) simplifies to the form In Cartesian coordinates (t, x 1 , x 2 , . . . , x d−1 ) a canonical form of the electromagnetic fieldstrength tensor F ab in d dimensions [9] For ease of writing we use the notation x 1 = x, x 2 = y and x d−1 = z. The equations of motion (3.21) are These equations are identical to the ones analysed by Ori [1]. In the following we make only brief comments about the solutions to these equations, referring the reader to that reference for a more detailed discussion. We are mostly interested in the motion of a massless particle in the (t, z) plane. The solution for equations (3.23) is: From these expressions we see that the vanishing point is located at finite spacetime point with coordinates (t = t c , z = z c ). Now we are in position to discuss the motion of the orbit beyond the vanishing point. For concreteness we assume that the Lorentz parameter increases as the vanishing point is approached, i.e, case (i) applies. Under this assumption the orbit approaching the vanishing point (t = t c , z = z c ) has parameterisation In particular, it has dz dt = −1, i.e., the z coordinate decreases as the orbit approaches z = z c . Once the vanishing point is reached either λ decreases from +∞ or it increases from −∞ depending on whether case (i) or case (ii) is realised. If beyond the vanishing point case (i) is realised then dz dt = −1, i.e., it maintains its sign. Therefore case (i) is the non-bouncing continuation. On the other hand, in case (ii) dz dt = +1, i.e., it changes sign across the vanishing point. Therefore case (ii) is the bouncing continuation.
To see which continuation is physical, one can do a perturbative analysis, analogous to the one performed in the previous section, by introducing a small mass parameter and studying timelike Lorentz orbits in flat spacetime with constant electromagnetic field strength. Such an analysis is performed in detail in reference [1]. Indeed, one finds that the bouncing continuation is the correct continuation beyond the vanishing point. That analysis applies as it is to the higher-dimensional case.

Charged Vaidya solution: bouncing continuation
In this section we apply the results of the previous sections to get a new interpretation of the charged Vaidya solution.

Interpretation
Recall that the charged Vaidya solution is supported by the external matter stress tensor (2.7). The external matter stress tensor is indeed of the form (cf. (2.8)), however, the normalisation of neither ρ nor that of the null vector k a is fixed by this expression alone. The extra physics input required to fix these normalisations is provided by the continuity equation (2.17). Assuming k a ∝ −∂ a v, i.e., k a = S(v, r)δ a r implies J a = ρk a = ρ(v, r)S(v, r)δ a r . Conservation of J a fixes its form to be where φ(v) is the flux of particles at retarded time v. A convenient numerical factor is taken out from definition of φ(v) in equation (4.2) . For the electromagnetic field (2.6), Maxwell equation (2.14) gives the form of electric current to be Comparing (4.2) and (4.3) we find charge per particle, .
In principle there is nothing wrong in taking the charge per particle e(v) also a function of the retarded time v. However, if one imagines the incoming particles to be of the same type at all times, then e(v) is simply a constant. The expressions for the external energy momentum tensor (2.7) and matter current (4.2) fix both S(v, r) and ρ(v, r), .
From the expression of S(v, r) we immediately see that where κ ∞ is interpreted as the energy of the infalling particle measured by a static observer at infinity, and Comparing (1.1), (3.6), and (4.9) we see that the critical surface is precisely at the place where the vanishing point of the null vector field k a is located. Equation (4.7) tells us that the kinetic energy of the thin shell at advanced time v vanishes at r = r c (v). The fact that ρ diverges at r = r c (v) is not of any consequence; it has no direct physical interpretation 3 . Moreover, one can easily check that the vector k a is consistent with the Lorentz force equation (2.18). We conclude that the ingoing charged Vaidya solution is only valid till the surface r = r c (v). In particular, we are gluing the spacetime M − with metric g − given by ingoing charged Vaidya solution

Gluing construction
and M + with metric g + given by outgoing charged Vaidya solution 4 3 Note that the number density of these particles measured by an observer is n ob = Jau a ob = ρkau a ob . This quantity is finite throughout. Since kv vanishes at r = rc the problematic factor in ρ is cancelled. 4 In this line element, coordinate u is thought of as minus the retarded null coordinate. on the surface Σ, where we have identified radial coordinate r as well as all angular coordinates across Σ, cf. figure 1. The function f − and f + are given by, and respectively; m − and q − are monotonically increasing positive functions of v. The bouncing surface Σ is represented in g − and g + by, 14) respectively, where, For r c + the dot denotes derivative with respect to the u coordinate and for r c − it denotes derivative with respect to the v coordinate.
The preliminary junction condition is then the agreement of the induced metric on Σ, where,f − (v) ≡ f − (v, r = r c − (v)) andf + (u) ≡ f + (u, r = r c + (u)). Then gluing amounts to finding a function, u = ψ(v), (4.18) defined on the surface Σ, such that the preliminary junction condition (4.17) holds and the extrinsic curvature remains continuous 5,6 across Σ. Identifying on Σ leads to r c + (u) = r c − (v). Inserting this into (4.17) we get, The solutions are,ψ The first solution,ψ = 1, corresponds to matching to the non-bouncing continuation and is thus discarded. Choosing ψ(0) = 0 we obtain a unique matching given by,

q(v) = ν m(v)
Let us now consider in asymptotically flat settings, a case in which the mass and the charge functions are proportional to each other, i.e., charge to mass ratio is constant as a function of v, where ν ∈ [0, 1) is a free parameter. A slightly more general situation is commented on in appendix A. With this choice of charge distribution the outer and inner apparent horizons are at while the radial-coordinate r on the bouncing surface Σ depends on v through, Clearly, the bouncing surface lies in between the inner and outer apparent horizon. The norm of the normal to the surfaces We noted above that the surface r = r c − (v) lies in between the inner and outer apparent horizons. Therefore, the expression in brackets in equation (4.29) is negative (assumingṁ − (v) > 0), i.e., the normal to the surface r = r c − (v) is timelike, We conclude that the bouncing surface is spacelike. However, since the inner apparent horizon lies in the future of this surface, its location must be given in terms of the u coordinate. For a simple situation we provide such an expression now. The matching condition reads, where µ is a parameter and v 0 is the width of the shell. Then the ingoing and outgoing coordinates are related on Σ by, The mass function of outgoing null fluid reads, where,μ and the charge function q + (u) = ν m + (u) remains proportional to mass function by the same factor ν. The inner apparent horizon is at,

f (R)-Maxwell set-up
One of the most popular higher-curvature extension of general relativity is the so-called f (R) gravity [13,14,15]. A general framework regarding spherically symmetric backgrounds in f (R) theories is well developed [16,17,18]. A large class of vacuum solutions are obtained requiring constant Ricci scalar curvature. For non-vacuum solutions, a requirement to obtain constant Ricci scalar curvature solutions is that the trace of the external stress-tensor should be zero.
Since we are interested in charged Vaidya solutions, the requirement of zero trace of the external stress-tensor forces us to work only in four-dimensions. Fortunately, non-static generalisations, in particular, the four-dimensional charged Vaidya solution with constant Ricci scalar has been worked out [19]. The Lagrangian for this theory is The equation of motion for the metric is [16,17,19] where T The only non-vanishing components of F ab are The external matter stress tensor required to support the solution is This stress tensor vanishes at the hypersurface Following the analysis of section 4.1 it is straightforward to show that in this case also, the critical surface is the vanishing surface of the null vector field satisfying the Lorentz force equation. To this end let k a = S(v, r)δ v a . Once again the continuity equation for the matter and electric current fixes the form of the S(v, r) and ρ(v, r). We find .
Using these expression we have (5.10) It can be readily checked that k a satisfies the Lorentz force equation, Qualitatively these expressions are very similar to the expressions in section 4.1. Therefore, similar considerations apply. The f (R) set-up we have worked with is essentially equivalent to general relativity with a cosmological constant [19]. That is why the analysis of this section is similar to that of the section 4.1. Nonetheless, we hope that our analysis is a small step that can be viewed as motivation for a more detailed analysis into other higher-curvature theories.

Discussion
We want to emphasise that in the usual interpretation of the charged Vaidya solution null geodesics do enter the r < r c region; this was precisely the source of confusion in the past literature [2,3,4,5] (see also [6]) . The pathology that there is a region in the spacetime where null and weak energy conditions are violated is tied to the interpretation of the charged Vaidya solution itself. In the correct interpretation, as first discussed by Ori [1], and further explored in this paper, the pathological region is not there. In this paper we highlighted the fact that the assumption that "charged photons" move on null geodesics is inconsistent with the equations of motion. We showed that at r = r c the Lorentz orbits bounce. The continuation of the Vaidya solution beyond r = r c surface is unfortunately not a simple matter. Only in the case when r = r c is a spacelike surface the question can be adequately addressed by gluing an outgoing charged Vaidya spacetime in the future of the r = r c surface. When r = r c surface is timelike the correct interpretation will involve cross-flowing Vaidya solutions, which to the best of our knowledge has no simple metric representation.
In AdS/CFT correspondence the importance of null energy condition has been realised in several situations. In an attempt to sharpen this connection further, reference [6] explored the strong sub-additivity property of the holographic entanglement entropy in the context of planar charged Vaidya spacetime; a set-up closely related to what we discussed above. They worked on the assumption that the vanishing surface r = r c separates the entire spacetime in two regions: (i) r < r c where null energy condition is violated, and (ii) r > r c where it is satisfied. They argued that spacelike minimal surfaces can enter and explore the r < r c region at least for certain choices of the mass and charge functions. Although intentionally focusing on an unphysical region of a spacetime is not justified from a General Relativity perspective, it can still be useful for other reasons. Knowing that the charged Vaidya metric for (0 < r < ∞) leads to an unphysical region where the null energy condition is violated, reference [6] argued that spacelike minimal surfaces entering the r < r c region possibly translates into a violation of strong sub-additivity of entanglement entropy in the dual field theory.
Kaminaga [5] in 1990 suggested a model for charged black hole evaporation where a negative energy flux of charged null fluid flows into the black hole. The negative energy causes the event horizon to shrink to zero size. Kaminaga studied the model from the classical and semiclassical perspectives. In this model it was found that the inner structure of geometry is very different from that of the Reissner-Nordström geometry. In this context, the effect of bouncing continuation of charged null orbits has been studied by Ori and Levin [20]. They found that the inner structure of black hole suggested by Kaminaga is further drastically altered when this effect is taken into account.
In the mass inflation scenario [21], and in particular in the Ori model of mass inflation [22], the charged Vaidya solution plays a very important role. In that context we recall that only the mass is taken to be a function of the advanced time v, not the charge, i.e., q(v) = constant. Therefore, the above details do not alter any of the discussions there.
In a different context, Frolov and Vilkovisky [23] explored effects of higher curvature corrections, as quantum gravity effects, and argued that these effects cause an incoming null shell to bounce back. Hayward [24] has explored similar ideas as a model for formation and evaporation of a non-singular black hole. Implications of these ideas for information paradox are discussed recently by Frolov [25]. There are intriguing similarities between our purely classical considerations and ideas put forward in these papers. It will be interesting to explore this link further. where, The metric functionf on the critical surface reads, Sincef ∼ v −2s as v → 0 for s > 0, the metric at the critical surface is singular. On the other hand for s < 0 there exist a v = v c given by, v |s| c = ν below which r ± are not real.