Page Curves of Reissner-Nordstr\"om Black Hole in HD Gravity

We obtain the Page curves of an eternal Reissner-Nordstr\"om black hole in the presence of higher derivative terms in four dimensions. We consider two cases: gravitational action with general ${\cal O}(R^2)$ terms plus Maxwell term and Einstein-Gauss-Bonnet gravity plus Maxwell term. In both the cases entanglement entropy of the Hawking radiation in the absence of island surface is increasing linearly with time. After including contribution from the island surface, we find that after the Page time, entanglement entropy of the Hawking radiation in both the cases reaches a constant value which is the twice of the Bekenstein-Hawking entropy of the black hole and we obtain the Page curves. We find that Page curves appear at later or earlier time when the Gauss-Bonnet coupling increases or decreases. Further, scrambling time of Reissner-Nordstr\"om is increasing or decreasing depending upon whether the correction term (coming from ${\cal O}(R^2)$ terms in the gravitational action) is increasing or decreasing in the first case whereas scrambling time remains unaffected in the second case (Einstein-Gauss-Bonnet gravity case). As a consistency check, in the limit of vanishing GB coupling we obtain the Page curve of the Reissner-Nordstr\"om black hole obtained in arXiv:2101.06867.


Introduction and Motivation
The black hole information paradox [2] has to do with the unitary evolution of the black holes. Page proposed that entropy of Hawking radiation of an evaporating black hole should fall after the Page time [3]. For eternal black holes entropy of the Hawking radiation should reach a constant value which is Bekenstein-Hawking entropy [4] of the black holes. To calculate the Page curve of a black hole one is required to couple the black hole to bath where we can collect the Hawking radiation. For eternal black holes we take two copies of the aforementioned setup to calculate the Page curve.
With the discovery of AdS/CFT correspondence [5] as a duality between conformal field theory and anti-de Sitter background in one higher dimensions, it become easy to resolve the black hole information paradox [2]. We can obtain the entanglement entropy of conformal field theory from gravity dual background using Ryu-Takayanagi formula [6] for static background. For time dependent background one is required to use the prescription given in [7]. Quantum corrections to all order in to the Ryu-Takayanagi formula was incorporated in [8] where one is required to extremize the generalised entropy. Surfaces which extremize the generalised entropy are known as quantum extremal surfaces(QES). If there are more than one quantum extremal surfaces then we need to consider the one with minimal area. In [9], authors generalised the QES prescription to island surfaces where we are required to extremize the generalised entropy like functional which includes contribution from the island surfaces. In this case extremal surfaces are known as quantum extremal islands.
We can calculate the Page curves of eternal black holes using Island proposal [9]. Island formula is given below: S gen (r) = Min I Ext I Area(∂I) 4G N + S matter (R ∪ I) , where R is the radiation region, G N is the Newton constant and I is the island surface. First term in the above formula is the area of the island surface and second term is the matter contribution to the entanglement entropy of the Hawking radiation which includes contribution from the island surface. Since initially there is no island surface therefore (1) reduces to the entanglement entropy of only radiation part which increases linearly with time. At late times island surface comes into the picture and saturates the entanglement entropy growth of the Hawking radiation at the Page time and produces the Page curve. In general higher derivative gravity we can write the above formula in most general form as, where S gravity and S matter are the gravitational and matter contributions to the total entanglement entropy. In higher derivative gravity theories first term in equation (2) can be calculated using Dong formula [10] and S matter can be calculated using Cardy formula [11,12]. Then we are required to extremize the total entanglement entropy with respect to location of the island surfaces. If there are more than one surface then we need to choose the surface with the minimal area among those surfaces. By following this approach Page curves of the Reissner-Nordström black hole, charged dilaton black hole, Schwarzschild black hole and hyperscaling violating black branes were calculated in [1,13,15,16]. Page curve in charged linear dilaton model for nonextremal black hole and the extremal black hole have been studied in [17] 1 . Recently islands in Kerr-de Sitter spacetime and generalized dilaton theories have been studied in [18,19]. Page curve and the information paradox for the flat space black holes was studied by authors in [20] 2 . The role played by mutual information of subsystems on the Page curve was studied in [21] 3 . Page curve calculations in higher derivative gravity theories can be found in [14,22]. In [14], authors calculated Page curves of the Schwarzschild black holes in the presence of higher derivative terms which are O(R 2 ) terms. By following [14], we are calculating the Page curves of an eternal Reissner-Nordström black hole in the presence of O(R 2 ) terms as considered in [14] and in Einstein-Gauss-Bonnet gravity [23] in four dimensions in this paper.
Charged black hole in Einstein-Gauss-Bonnet gravity in four dimensions were studied in [24][25][26][27][28][29][30]. In [24], authors studied superradiance, stability and quasinormal modes of regularised charged Einstein-Gauss-Bonnet black hole. In [25], authors have discussed thermodynamics of the charged Einstein-Gauss-Bonnet black hole via van der Walls equation. In [26], thermodynamics and phase transitions of charged AdS black holes was studied, additionally quasinormal modes for the massless scalar perturbation was also studied in the same paper. In [27], correspondence between shadow and test field was studied for the charged Einstein-Gauss-Bonnet black hole. In [28], charge of the Einstein-Gauss-Bonnet black hole was studied by considering negative and positive charge of the particle-antiparticle pair on cauchy horizon and physical horizon. In [29], instability of charged Einstein-Gauss-Bonnet de-Sitter black holes was studied by authors under charged scalar perturbations. In [30], authors have discussed motion of the charged and spinning particles and photons in the vicinity of charged Einstein-Gauss-Bonnet black hole in four dimensions. Further radius of the innermost circular orbit, gravitational deflection angle and some collision processes were also studied by the authors. For review on Einstein-Gauss-Bonnet gravity in four dimensions, see [31].
As summarised above literatures on the charged Einstein-Gauss-Bonnet black hole in four dimensions. We find that study of effect of the Gauss-Bonnet coupling on the Page curve of 1 We thank H. S. Jeong to bring their work to our attention. 2 We thank C. Krishnan to bring their work to our attention. 3 We thank A. Saha to bring their work to our attention.
charged Einstein-Gauss-Bonnet black hole was missing. Therefore it will be interesting to study that how the Page curves of Reissner-Nordström black hole will be modified with the presence of Gauss-Bonnet term. With this motivation we have calculated the Page curves of charged black hole in higher derivative gravity with O(R 2 ) terms and for the charged Einstein-Gauss-Bonnet black hole in this paper.
Structure of the paper is as follow. Section 2 has been divided into two subsections 2.1 and 2.2. In subsection 2.1 we review the charged black hole solution in Einstein-Gauss-Bonnet gravity in four dimensions [23] and in subsection 2.2 we review the Page curve calculation of the Reissner-Nordström black hole [1]. In section 3 we have calculated the Page curves of an eternal Reissner-Nordström black hole in the presence of O(R 2 ) terms in the gravitational action. In section 4 we have calculated the Page curves of the same black hole in Einstein-Gauss-Bonnet gravity. In section 5 we have some discussion about the island rule which also includes the summary of the results obtained in this paper.

Review
We have divided this section into two subsections. In subsections 2.1 and 2.2, we are briefly reviewing the charged black hole in Einstein-Guass-Bonnet gravity in the vanishing cosmological constant limit [23] and Page curve calculation of Reissner-Nordström black hole [1] in four dimensions.

Brief Review of Charged Black Hole in Einstein-Gauss-Bonnet Gravity
In this section we are reviewing spherically symmetric charged black hole solution in Einstein-Gauss-Bonnet gravity in four dimensions [23]. Authors in [32] have shown that we can make the Gauss-Bonnet term in four dimensions dynamical if we rescale the Gauss-Bonnet coupling as α → α D−4 . Consistent theory of Gauss-Bonnet gravity was constructed by authors in [33] 4 in (d + 1)-dimensions where all the relevant degrees of freedom is present. Action of the Aoki, Gorji and Mukhohyama (AGM) theory is: where the Gauss-Bonnet term. If we take the d → 3 limit of equation (3) then we obtain the consistent theory of Gauss-Bonnet gravity in 4 We thank S. Siwach to bring [33] to our attention.
four dimensions for the neutral black hole. Author in [23] constructed charged Einstein-Gauss-Bonnet black hole solution in four dimensions by rescaling of the Gauss-Bonnet coupling. We are discussing this solution in the vanishing cosmological constant limit. Action of the Einstein-Gauss-Bonnet gravity including Maxwell term is: where, R[g] is the Ricci scalar, α is the Gauss-Bonnet coupling, L GB is the Gauss-Bonnet term and F µν is the field strength tensor corresponding to A µ gauge field which are defined below: Black hole solution of the action (4) is: where, If we perform small expansion in α with negative sign chosen in (7) then we obtain: Black hole horizons are given by solving F (r) = 0 for the negative sign in equation (7), which are given below: where r + is the physical horizon of the charged black hole and r − is the Cauchy horizon. Now let us see how the black hole horizon r + is varying with the Gauss-Bonnet coupling α. We have plotted it in figure 1 for M = 1 and various values of the black hole charges.
From the figure 1 and equation (9), we can see that when Gauss-Bonnet coupling (α) is zero then we have usual Reissner-Nordström black hole. When Gauss-Bonnet coupling is decreasing towards negative values then black horizon is increasing and when Gauss-Bonnet coupling is increasing towards positive values then black horizon is decreasing.
Authors in [24,27] have discussed that for M = 1 following are the allowed values of the Gauss-Bonnet coupling α.
As discussed in [24,27], A region contain both the horizons (r ± ) whereas B region contain only the physical horizon (r + ). Since we are interested in the non-extremal black holes therefore relevant region for us is the region A. Numerically in region A, α ∈ (−6.66, 0.36) for Q = 0.8. Authors in [34,35] studied the effect of Gauss-Bonnet term on the shear viscosity to entropy density ratio η s , causality violation and instability of the black brane in five dimensions and arbitrary D dimensions from the perspective of gauge-gravity duality. They found that there is an upper bound on the Gauss-Bonnet coupling, λ ≤ 1 : α appear as Gauss-Bonnet coupling in [34,35]) when D → ∞ and, λ ≤ 0.09 when D → 5. There is instability of RN-AdS black brane in five dimensions when 0 < λ ≤ 0.09.
Charged black hole in Einstein-Gauss-Bonnet gravity were studied in detail in various papers mentioned earlier in the introduction section 1. In this paper we will calculate the Page curves of the charged black hole in the presence of higher derivative terms which are O(R 2 ) terms and then we focus explicitly on the Einstein-Gauss-Bonnet gravity [23]. Since we are working with higher derivative gravity therefore we will use [10] to calculate the entanglement entropy wherever required for higher derivative terms in this paper.
In this paper we are focusing only on the non-extremal black holes. Page curve of the extremal black holes were studied in [17,36,37]. Authors in [17,36] studied the Page curve of eternal charged linear dilaton black holes and Reissner-Nordstrom black holes for the non-extremal and extremal both cases and found that island formulation is ill-defined for the extremal black holes, which implies that we may not obtain the Page curve for the extremal black holes. Authors in [37] considered rotating BTZ(Bañados-Teitelboim-Zanelli) black holes in three dimensions and found that for the non-extremal rotating BTZ black holes one obtains the Page curve using island formulation where as for the extremal rotating BTZ black holes, Page time and scrambling time turn out to be divergent and this problem can be removed by considering superradiance phenomenon in the system. With the inclusion of superradiance, Page time and scrambling time were decreasing as the angular momentum of the black holes were increasing and interestingly one obtains the Page curve after finishing of superradiance.

Review of Page Curve of Reissner-Nordström Black Hole
In this section we are reviewing the Page curve calculation of Reisnner-Nordström black hole done in [1]. Action for the Einstein-Maxwell system is: where R is the Ricci scalar, F µν is the field strength tensor corresponding to the A µ gauge field and I matter is the matter contribution to the action. Metric of the Reissner-Nordstörm black hole is: where dΩ 2 = dθ 2 + sin 2 θdφ 2 and, Black hole horizons will be given by solution to the F (r) = 0, and are given below: We can write metric (11) in Kruskal coordinates as given below [41]: where and r * is defined as: κ ± are the surface gravities corresponding to r + and r − and are defined as: Hawking temperature and Bekenstein-Hawking entropy of the charged black hole are defined below: Defining the conformal factor in the metric (14) as: Therefore we can write the metric (14) in the following form: Penrose diagrams of an eternal Reissner-Nordström black hole are given in figures 2 and 3 where R − and R + are the left and right wedges of the radiation regions, b − and b + are the boundaries of the R − and R + , and, a − and a + are the boundaries of the island surface in the left and right wedges. 5 .

Entanglement Entropy without Island
Initially there is no island surface therefore entanglement entropy of the Hawking radiation can be calculated using the following formula [11][12][13]: are the boundaries of the radiation regions in the right and left wedge of Reissner-Nordström black hole. Formula (21) is valid for two dimensional vacuum CFT. Since entanglement entropy formula for higher-dimensional spacetime is not known. But we can apply (21) for higher-dimensional spacetime in s-wave approximation [15] by assuming that location of the observer or cut-off surface (b ± ) is very far away from the black hole horizon (r + ), i.e., b ± r + . Therefore in s-wave approximation we can ignore the angular part of the metric and we will be left with asymptotically flat black hole for which (21) can be used . β is  the inverse of Hawking temperature and defined as, β = 2π/κ + . Geodesic distance between two points (l 1 , l 2 ) is given as: Now using equations (15), (16), (19), (21) and (22) entanglement entropy of the Hawking radiation in the absence of island surface simplifies to [1]: when t → ∞, i.e., at late times, writing cosh κ + t ∼ e κ + t , time dependent part of equation (23) simplifies to: Therefore from equation (24) we can see that entanglement entropy of the Hawking radiation in the absence of island surface is increasing linearly with time and becomes infinite at late times which leads to the information paradox for the Reissner-Nordström black hole. In the next subsection we will see that at late times island appears and entanglement entropy of the Hawking radiation in the presence of island surface will be constant and dominates after the Page time. On combining the both contributions, we obtain the Page curve.

Entanglement Entropy with Island
In this subsection we are reviewing calculation of the entanglement entropy of the Hawking radiation in the presence of island surface. We can calculate it using the following formula [11][12][13]: where a + (t a , a) and a − (−t a + ι β 2 , a) are the boundaries of the island surface in the right and left wedge of Reissner-Nordström geometry. Entanglement entropy formula (25) can be used for higher-dimensional spacetime in s-wave approximation where we can ignore the angular part of the metric for a distant observer (b ± r + ). At late times, t a , t b b > r + , from equations (15), (16), (19) and (22), equation (25) simplifies to [1]: Hence generalised entropy in the presence of island surface using (1) at late times is given by [1]: where first term in the above equation is coming from the area of the island surface. On varying equation (27) with respect to a, i.e., ∂S late times gen (a) ∂a = 0, we are required to solve the following equation: solution to the above equation is, Substituting the value of a from equation (29) in equation (27), we obtain the total entanglement entropy of the Hawking radiation at late times as: From the above equation we can see that total entanglement entropy is constant which dominates after the Page time and we obtain the Page curve of an eternal Reissner-Nordström black hole.

Page Curve
From equation (24) we saw that entanglement entropy in the absence of island surface is increasing linearly with time and from equation (30) we saw that entanglement entropy is constant at late times. Therefore initially entanglement entropy of the Hawking radiation will increase and after the Page time island surface emerges which saturates the linear growth of entanglement entropy and reaches a constant value, which is twice of the Bekenstein-Hawking entropy of the Reissner-Nordström black hole, and we obtain the Page curve of the Reissner-Nordström black hole.
We have plotted the Page curve of the Reissner-Nordström black hole for, M = 1, Q = 0.8, G N = 1 in figure 4 by considering only leading order term in G N in equation (30). In figure 4, blue line corresponds to the linear time growth of the entanglement entropy (equation (24)) and orange line corresponds to the entanglement entropy in the presence of island surface (equation (30)) which dominates after the Page time and leads to the Page curve.
Page time: Page time is defined as time at which entanglement entropy of the Hawking radiation start falling to zero for an evaporating black hole and reaches constant value for an eternal black hole. We can get the Page time by equating equations (24) and (30), which turn out to be: Scrambling time: Scrambling time is defined as time interval in which we recover the information thrown into the black hole in the form of Hawking radiation [38,39]. It was discussed in [38] that for an evaporating black hole information can be quickly retrieved in the form of Hawking radiation if the black hole has evaporated half away. If the black hole has not evaporated half away then one has to wait to evaporate half away, after that information can be quickly recovered. In the language of entanglement wedge reconstruction [40], scrambling time is defined as time at which information from the cut-off surface (r = b + ) reaches the boundary of the island surface (r = a + ). When information reaches to the boundary of island surface (r = a + ) then entanglement entropy of Hawking radiation includes contribution from the island surface because island degrees of freedom becomes part of the entanglement wedge of the Hawking radiation and we start recovering the information thrown into the black hole in the form of Hawking radiation.
If we want to send the information into the black hole from the cutoff surface (r = b) then time taken by information to reach the boundary of island surface (r = a) is given by: Substituting the value of a from equation (29) in the above equation, scrambling time of the Reissner-Nordström black hole turns out to be: where T RN is the Hawking temperature and S BH is the Bekenstein-Hawking entropy of the Reissner-Nordström black hole. From equation (33) we can see that scrambling time turns out to be logarithmic of black hole thermal entropy and this result is similar to [39] which implies that black holes are fastest scramblers.

Page Curves of Charged Black Hole in HD Gravity up to O(R 2 )
In this section we are considering general higher derivative terms at O(R 2 ) in the gravitational action as considered in [14] and we are going to calculate the Page curves of an eternal Reissner-Nordström black hole in the presence of O(R 2 ) terms. In our case we have charged black hole whereas in [14] authors obtained the Page curve of Schwarzschild black hole. Action with inclusion of O(R 2 ) terms and Maxwell-term is: where R[g] is the Ricci scalar, R µν [g] is the Ricci tensor and L GB is the Gauss-Bonnet term, in four dimensions. Wald entropy of the four dimensional charged black hole in higher derivative gravity can be calculated using the formula given in [42,43] and used in [44] to calculate Wald entropy in the presence of O(R 4 ) terms in M-theory action [45]. The formula is: where h is the determinant of the induced metric on co-dimension two surface. Writing the metric (11) in the following form: From the action (34) we find that: where the terms written in the square bracket will also be used in calculation of Wald entropy of charged black hole in Einstein-Gauss-Bonnet gravity in section 4. For the metric (36), Wald entropy corresponding to the action (34) using equation (37) turns out to be: Since r + r − , therefore above equation simplifies to: In general higher derivative gravity theories holographic entanglement entropy can be calculated using the following formula [10]: where wher m (a) µ are unit orthogonal vectors along y a (tangential) directions, n (i) µ are the untit normal vectors along the normal directions, K λµν is the extrinsic curvature and q α in equation (40) is a number which is sum of the certain combination of Riemann tensor components and extrinsic curvatures which can be obtained by first calculating the second order derivative of Lagrangian and labelling each term by α and then performing expansion of Riemann tensor components along certain directions (for more details, see [10]). We can write the total entanglement entropy in the presence of higher derivative terms as: For the action (34), gravity contribution to the entanglement entropy turns out to be [14]: where i corresponds to the normal directions, K i being trace of extrinsic curvature, K i,µν = −h α µ h νβ , h µν is the induced metric on boundary of the island surface. Further, G N,ren , λ 1,ren , λ 2,ren and α ren are the renormalised Newton constant and coupling constants appearing in higher derivative gravity action (34). These are used to absorb the UV divergences of the von Neumann entropy of matter field discussed in [14]. Matter contribution to the entanglement entropy of the Hawking radiation in the presence of island surface can be obtained using equations (15), (16), (19), (22) and (25) and is given below [1] 6 : At late times, t a , t b b > r + , above equation simplifies to the following form: Therefore generalised entropy for the action (34)  +Smatter(R ∪ I) .

(46)
For the metric (11), we find that: , n r 2 = F (r),K 1 = 0 and i n ν i = 0, Hence gravitational contribution to the entanglement entropy of the Hawking radiation (43) simplifies to, Therefore total generalised entropy will be given by the sum of entanglement entropy contri- 6 Matter contribution to the entanglement entropy, Smatter(R ∪ I), up to some extent is same as [1]. Therefore we are writing the final result up to that point.
butions from gravity part (48) and matter (45) part both as: S gen (a) = S gravity (a) + S late times where r * (b), g(a) and g(b) can be substituted from equations (16) and (19) in equation (49). t b dependent term in equation (45) is small and hence has been ignored. Location of the island surface can be found by extremizing the above equation with respect to a, i.e.
∂S gen (a) ∂a ∼ 4π −2M r + λ 2,ren + 2Q 2 λ 2,ren + r 4 solution to the above equation is, From the above equation we can see that island lies outside the black hole horizon and this leads to causality paradox. It was shown in [46] that this result appears in all the two sided eternal black holes or black holes in Hartle-Hawking state and one can restore the causality by quantum focusing conjecture (QFC) [47]. Idea of [46] is that when we decouple the black hole from the bath then finite amount of energy flux will be produced and this energy flux pushes the black hole horizon outwards and therefore island always lies behind the horizon. It was discussed in [48] that finite amount of energy will also be produced even when we couple the black hole to bath and this energy flux pushes the horizon outwards which implies that island lies behind the horizon similar to the decoupling process and we can get rid of causality paradox.
Substituting value of a from equation (51) in equation (49), generalised entropy simplifies to the following form 7 : Therefore total entanglement entropy in the presence of island surface is twice of the Bekenstein-Hawking entropy of the black hole plus matter contribution. If we consider only leading order term in G N then we can see that total entanglement entropy of an eternal Reissner-Nordström black hole in the presence of O(R 2 ) terms reaches a constant value consistent with the literatures. It is interesting to notice that λ 2,ren dependence is appearing in higher order in G N similar to [14]. In the α → 0 then equation (54) reduces to, which is the almost same as in [1]. From equation (54), it is clear that total entanglement entropy at late times is constant. Combining equations (24) and (54)  ]. Now using equation (19). This term can be written as: Now substituting a ≈ r+ in the above equation we obtain: This will not affect the Page curve of the Reissner-Nordström black hole in the α → 0 limit because to obtain the Page curve we are considering only leading order term in equations (54) and (74) which in the α → 0 limit reduces to leading order term in equation (30). Same thing can also be verified by substituting the value of a from equation (51) in equation (49) and performing small GN expansion and retaining the terms up to O(G 0 N ). we obtain: where α in the above equation (and anywhere in this section) is α ren , for the sake of simplicity we have just written α.   (24) and (54) and is given below: From the above equation it is clear that in α → 0 limit, the Page time of the Reissner-Nordström black hole in the presence of O(R 2 ) terms reduces to Page time of the Reissner-Nordström black hole (31) without higher derivative terms.
Since in equations (54) and (57) , Gauss-Bonnet coupling (α) is coming with positive sign, therefore we obtain the Page curves at later or earlier time when α is increasing or decreasing.
Scrambling time: Substituting the value of a from equation (51) in equation (32), scrambling time of the charged black hole in the presence of O(R 2 ) terms turns out to be: In the above equation if we take λ 2,ren → 0 limit then we obtain the scrambling time of the Reissner-Nordström black hole (33) in the absence of higher derivative terms. Performing small λ 2,ren expansion of the equation (59), we obtain t HD scr ∼ From the above equation it is clear that scrambling time of the charged black hole in the presence of higher derivative terms will increase or decrease depending on whether second term in equation (60) is positive or negative.

Page Curves of Charged Black Hole in Einstein-Gauss-Bonnet Gravity
In this section we are going to calculate the Page curves of charged black hole in Einstein-Gauss-Bonnet gravity using Dong formula (40). Working action with vanishing cosmological constant is [23]: where R[g] is the Ricci scalar, L GB is the Gauss-Bonnet term (5) and F µν is the field strength tensor corresponding to gauge field A µ in four dimensions. Wald entropy of the charged black hole in Einstein-Gauss-Bonnet gravity can be obtained from equations (35), (36), (37) and (61) and is given as: Holographic entanglement entropy for the action (61) was calculated in [10,49,50] and is given below: where h is the determinant of the following induced metric on constant t, r = a surface: Now we are going to simplify equation (63) using Gauss-Codazzi equation [14,50], From equations (63) and (65) entanglement entropy in Einstein-Gauss-Bonnet gravity simplifies to the following form (which is the gravitational contribution the generalised entropy): and matter contribution to the entanglement entropy is [1] 8 : Entanglement entropy at late times, t a , t b b > r + is given by the following expression: where r * (b) and g(a), g(b) can be obtained from equations (16) and (19). From equations (66) and (68), generalised entropy of the charged black hole in the Einstein-Gauss-Bonnet gravity is: where first term in the above equation is coming from the first term in equation (66). For the induce metric (64), we find that: R[∂I] = 2 a 2 . Therefore gravitational contribution to entanglement entropy (66) simplifies to: Now generalised entropy will be given by sum of equations (68) and (70), where we did not consider the small t b dependent term in equation (71). Location of the island surface can be found by extremizing the above equation with respect to a, i.e., equation (72) has the following solution, Now substituting value of a from equation (73) in equation (71), total entanglement entropy of the charged Einstein-Gauss-Bonnet black hole at late times (71) simplifies to the following form 9 : Since for the charged black hole in Einstein-Gauss-Bonnet gravity total entanglement entropy at late times (74) is same as total entanglement entropy for the charged black hole in higher derivative gravity with O(R 2 ) terms given in equation (54). Therefore Page curves of the charged black hole in Einstein-Gauss-Bonnet gravity will be same as the Page curves of charged black hole in higher derivative gravity with O(R 2 ) terms (figure 5).
Page Time: Since at t = t EGB Page , S WI EE = S EGB total , therefore on equating (24) and (71), we obtain the Page time of the charged Einstein-Gauss-Bonnet black hole as given below: From equations (71) and (75) we can see that in the α → 0 limit we obtain the total entanglement entropy at late times and Page time of charged black hole without higher derivative term as reviewed in section 2.2. Physical significance of the Gauss-Bonnet term is that when Gauss-Bonnet coupling is increasing then we are getting the Page curves at later time and when Gauss-Bonnet coupling is decreasing then we are getting the Page curves at earlier time with reference to Page curve of charged black hole without higher derivative terms, i.e., Page curve of the Reissner-Nordström black hole.
Scrambling time: In this case scrambling time is unaffected by the higher derivative terms and is same as scrambling time of the Reissner-Nordström black hole [1],

Conclusion and Discussion
In this paper we have obtained the Page curves of an eternal Reissner-Nordström black hole in the presence of higher derivative terms where entanglement entropy had been calculated using the formula given in [10].
Island rule was originated by coupling the evaporating JT(Jackiw Teitelboim) black hole plus conformal matter to two dimensional CFT bath [9]. Idea was to consider a black hole in JT gravity plus conformal matter (JT gravity + 2D CFT) to a CFT bath. There are three equivalent descriptions of this setup.
• 2D-Gravity: Black hole in JT gravity plus matter theory is coupled to 2D CFT bath where we are collecting the Hawking radiation.
• 3D-Gravity: The CFT bath has its own holographic dual which is AdS 3 .
• QM: Boundary of JT gravity plus matter theory is one dimensional QM system. Third description is 2D CFT bath with some boundary QM degrees of freedom.
Island rule for evaporating black holes in the above setup arises in a beautiful way that extra dimension(radial coordinate) in AdS 3 background connects the CFT bath to a region in the interior of the black hole which is known as, "Island", which appear at late times. Initially there is no island and one finds that entanglement entropy of Hawking radiation is increasing linearly with time and after the Page time island comes into the picture (which gives finite contribution to the entanglement entropy) and we obtain the Page curve. Above discussion depends strongly on holographic duality.
Island formula can be derived from gravitational path integral using replica trick for special JT black holes as done in [51,52]. It was discussed by authors that one obtains the Page curve of eternal black holes in the following way. There are two saddles: disconnected and connected. Below the Page time disconnected saddle dominates and after the Page time connected saddle dominates which is replica wormhole. Disconnected saddle is responsible for the linear time growth of entanglement entropy of the Hawking radiation and connected saddle gives the finite contribution. Combination of contributions from the disconnected and connected saddles to the entanglement entropy of Hawking radiation reproduces the Page curve. The argument of [51] is also applicable to n boundary wormholes which are known as replica wormholes.
In this paper we are focusing on non-holographic models therefore we should be careful when we will be discussing the higher dimensional asymptotically flat black holes. Since we are assuming that observer is very far away from the black hole therefore we can use s-wave approximation to calculate the entanglement entropy of Hawking radiation using formula for 2D CFT. Most of the papers existing in the literature discuss the application of island proposal to higher dimensional asymptotically flat black holes for the eternal black holes only. There is a paper [14] in which authors have discussed the Page curve in higher dimensional evaporating black hole in the presence of higher derivative terms in the gravitational action.
Effect of higher derivative terms on the Page curve had been explored in [14] for the neutral black hole. Authors considered the eternal black hole as well as evaporating black hole in their study. They focus only on the general O(R 2 ) terms in the gravitational action as higher derivative terms. We consider the non-extremal Reissner-Nordström black hole and we are focusing only on the eternal black hole. It is nice to study the effect of the Gauss-Bonnet coupling on the Page curves of Reissner-Nordström black hole. Presence of higher derivative terms in the gravitational action affect the Page time, scrambling time and Page curves of eternal black holes. Following are the main results obtained in this paper.
1. In the first case we considered the general O(R 2 ) terms in the gravitational action as discussed in [14] and calculated the Page curves. We find that, since initially there is no island surface, the entanglement entropy of the Hawking radiation increases linearly with time forever. Therefore we have information paradox for the charged black hole in the presence of O(R 2 ) terms. At late times island emerges and entanglement entropy of the Hawking radiation reaches a constant value which is twice of the Bekenstein-Hawking entropy of the black hole and we obtain the Page curves for fix values of the Gauss-Bonnet coupling (α). In this case we find that as Gauss-Bonnet coupling (α) increases Page curves shift towards later times and when Gauss-Bonnet coupling (α) decreases Page curves shift towards earlier times.
2. In the second case we considered only Gauss-Bonnet term as higher derivative term which is relevant to studying the charged black hole in Einstein-Gauss-Bonnet gravity [23]. Similar to the first case we calculated the Page curves in this setup, i.e., of charged Einstein-Gauss-Bonnet black hole [1]. Similar to the first case we have linear time growth of the entanglement entropy of the Hawking radiation at earlier times and island emerges at late times which saturates the linear time growth of the Hawking radiation and entanglement entropy of Hawking radiation becomes equal to twice of the Bekenstein-Hawking entropy of the black hole. Therefore we obtain the Page curves of the charged Einstein-Gauss-Bonnet black hole in four dimensions. Interestingly we find a similar behaviour of the Page curves with Gauss-Bonnet coupling as discussed earlier in the first case.
3. In the first case scrambling time of the charged black hole is affected by the higher derivative terms. It will increase if the correction term is positive and decrease if the correction term is negative. Further if we take the coupling appearing in scrambling time to zero then we recover the scrambling time of the Reissner-Nordström black hole. In the second case, i.e., when we consider only Gauss-Bonnet term as higher derivative term, scrambling time is unaffected by the higher derivative term and is same as scrambling time of the Reissner-Nordström black hole.
4. In the limit of vanishing Gauss-Bonnet coupling, we obtain the Page curve of the Reissner-Nordsrtöm black hole as obtained in [1].