PLB140981

140981

S0370-2693(26)00831-210.1016/j.physletb.2026.140981

Fig. 1

Radiation entropy Srad(t, N)/N as a function of rescaled time t/tPage for different values of the Bekenstein–Hawking entropy N. The dashed black curve shows the large-N result (sharp kink at t=tPage). Coloured curves show the finite-N prediction from Eq. (29): the kink is progressively smoothed as N decreases. All curves were computed analytically using the log-sum-exp formula with α(N)=logN/(2N).

Fig. 2

Entanglement entropy (in nats) from the random unitary circuit simulation with nq=6 qubits (40 steps, 1 realisation shown). The simulation begins in a product state |0⟩⊗6 and evolves under Haar-random unitaries. After a brief equilibration the entropy saturates near the Page value SPage=nq2log2≈2.08 nats (dashed red), confirming that the finite system never exceeds this bound. The uncorrected Hawking entropy (dotted orange) grows unboundedly, illustrating the information paradox. The saturation at finite N is the key finite-dimensional effect captured by our analytic formula.

Algorithm 1

Finite-N page curve simulation via random unitary circuits.

Table 1

Summary of parameters and results for the random unitary simulation.

Parameternq=4nq=6nq=8nq=10Hilbert space dim d16642561024Neff=nq2log21.3862.0792.7733.466SPage (nats)1.3862.0792.7733.466α(Neff)0.2500.1820.1500.128Δtsmooth (steps)4.05.56.77.8ReviewNon-perturbative transitions in the island formula: A finite-N resolution of the page curve kink

Non-perturbative transitions in the island formula: A finite-N resolution of the page curve kink

Snigdha

Saha

Snigdhasaha2938@gmail.com

Independent resarcher, BangladeshIndependent resarcher Bangladesh

Independent resarcher, Bangladesh

Editor: G. Mandal

Abstract

The island formula, derived from replica wormhole saddles in the gravitational path integral, successfully reproduces the unitary Page curve for black hole radiation entropy in the large-N (large central charge) limit. However, the resulting entropy profile contains an unphysical sharp kink at the Page time tPage—a mathematical artefact of the min operation that dominates in that limit. Any finite quantum system must exhibit a smooth, analytic transition. In this paper we work within Jackiw–Teitelboim (JT) gravity coupled to a holographic CFT bath and derive systematically the leading finite-N correction to the radiation entropy,Srad(t,N)=SHawk(t)+δS(N,t),δS∼O(e−SBH),originating from non-perturbative interference between the Hawking and replica wormhole saddles. Using a soft-mode resummation of the gravitational path integral we show that these corrections replace the sharp kink with a log-sum-exp interpolation controlled by 1/N. We further study an interference term arising when both saddles contribute at comparable weight, providing the first analytic description of a smooth cross-over in the island formula. Our predictions are confirmed by a random-unitary-circuit simulation of an n-qubit Page transition and by direct numerical evaluation of the partition function in a matrix-model truncation of JT gravity. The results imply that information begins to leak gradually from the black hole interior long before tPage, consistent with unitarity but without invoking a firewall.

Keywords

Black hole information paradoxIsland formulaReplica wormholesFinite-N correctionsPage curveJT gravity

Data availability

Data will be made available on request.

Introduction

The question of whether black hole evaporation is a unitary process has occupied theoretical physicists for five decades [1]. Hawking’s semiclassical calculation showed that the radiation emitted by a black hole is thermal: the reduced density matrix of the radiation is a Gibbs state, implying that the von Neumann entropy S=−Tr(ρlogρ) grows monotonically without bound [1,2]. This behaviour directly contradicts the Page curve [3,4], which stipulates that the radiation entropy of a unitary evaporation process must mirror the black hole entropy—rising until the Page time and falling thereafter.1.1

The emergence of the island formula

A major breakthrough arrived when Almheiri, Mahajan, Maldacena, and Zhao [5] and, independently, Penington [6] showed that gravitational effects—specifically replica wormhole saddle-point contributions to the gravitational path integral—could correct the Hawking entropy calculation and reproduce a Page curve consistent with unitarity. This led to the island formula [7]:(1)Srad=minI[Area(∂I)4GN+SQFT(rad∪I)],where the minimisation is over all candidate island regions I in the black hole interior, ∂I denotes the boundary of I (the quantum extremal surface), and SQFT is the entanglement entropy of quantum fields on the region rad∪I. At early times the no-island saddle dominates, giving Hawking’s growing entropy. At late times the island saddle dominates, giving a decreasing entropy that mirrors the Bekenstein–Hawking entropy of the shrinking black hole.1.2

The kink problem and research gap

The island formula as written in Eq. (1) involves an explicit min function. In the large-N (or equivalently large central charge c → ∞) limit the two saddles are sharply separated in weight, and the partition function is dominated by precisely one saddle at any given time. The resulting Page curve therefore has a sharp kink at t=tPage—a C0 but not C1 transition. This is mathematically the origin of the kink: it is not physical, but rather an artefact of the saddle-point approximation.From the perspective of a finite quantum system with Hilbert space dimension eN (where N is the Bekenstein–Hawking entropy in appropriate units), the true entropy curve must be:

(i)Smooth: the von Neumann entropy of any finite quantum system is an analytic function of time whenever the Hamiltonian is time-independent or adiabatically varying.

(ii)Bounded: Srad≤logdimHrad, which saturates only for a maximally mixed state.

(iii)Monotone-bounded: after the Page time the entropy is not strictly decreasing; fluctuations of order e−N/2 persist.

None of these properties are captured by the large-N island formula. The present paper addresses this gap.1.3

Main results and outline

Our principal results are as follows:

•We derive the finite-N corrected radiation entropy in JT gravity, showing that the Hawking and island saddles are joined by a smooth log-sum-exp interpolation with width Δt∼1/(NS˙BH).

•We identify a novel saddle interference term δSint that is exponentially small in N individually but becomes the leading correction near the Page time.

•We demonstrate via a random-unitary-circuit simulation that the finite-N curve matches our analytic prediction with no free parameters (after rescaling).

•We show that the smooth transition implies gradual information leakage beginning at t ≪ tPage, consistent with unitarity and with a smooth horizon (no firewall).

The paper is organised as follows. Secion 2 sets up the JT gravity model. Secion 3 derives the replica partition function and its saddles. Secion 4 derives the finite-N correction. Secion 5 presents numerical results. Secion 6 discusses physical implications, and Secion 8 concludes.2

Theoretical framework: JT gravity and the evaporating black hole

2.1

JT gravity action

We work with Jackiw–Teitelboim gravity in two spacetime dimensions coupled to a holographic CFT with central charge c [7–9]. The Lorentzian JT action is(2)IJT=ϕ016πGN∫MR+116πGN∫Mϕ(R+2)+ϕ08πGN∫∂MK+18πGN∫∂Mϕ(K−1),where ϕ is the dilaton field, ϕ0 is its extremal value, R is the Ricci scalar, and K is the extrinsic curvature of the boundary ∂M. The Euclidean path integral is(3)Z=∫DgDϕe−IE[g,ϕ],where IE is the Euclidean continuation of Eq. (2).On-shell, the bulk equations of motion set R=−2, fixing the metric to be locally AdS2. The only remaining dynamical degree of freedom is the boundary dilaton value ϕb=ϕ|∂, which in the Schwarzian description satisfies(4)SSchw=−ϕr16πGN∫0βdτ[{f,τ}+2π2β2(f′)2],where {f, τ} denotes the Schwarzian derivative and ϕr is the renormalised dilaton coupling.2.2

Black hole thermodynamics in JT gravity

The Bekenstein–Hawking entropy and temperature of the JT black hole are(5)SBH=ϕh4GN=ϕ04GN+2π2ϕrT,(6)TH=1β=ϕrϕhπ,where ϕh is the dilaton value at the horizon. We define the dimensionless parameter(7)N≡SBH(t=0)=ϕh(0)4GN,which plays the role of the large parameter in our expansion. Throughout the paper N ≫ 1 (semiclassical regime) but finite.2.3

Setup: evaporating black hole coupled to a bath

We couple the JT black hole to a non-gravitating CFT bath at temperature Tbath < TH, following [7]. The radiation region R consists of the exterior half-line in the bath. The island I=[b−,b+] is a connected interval inside the horizon, with endpoints (quantum extremal surfaces) b ±  determined by extremising the generalised entropy:(8)∂Sgen∂b±=0,Sgen=ϕ(b+)+ϕ(b−)4GN+SQFT(R∪I).For the linear-dilaton background at late times (t ≫ β) we have(9)SHawk(t)=c6logtϵ≈c6·ttPage·N,(10)Sisl(t)=2N−c6·ttPage·N=2N−SHawk(t),where ϵ is a UV cutoff and we normalise so that tPage=1 in units where c/6=N. The island saddle coincides with SBH(t) up to subleading terms.3

Multi-saddle replica analysis

3.1

Replica trick and the Rényi entropy

The von Neumann entropy is obtained via the replica trick:(11)SvN(ρR)=−limn→1∂nTrρRn=−limn→1∂nlogZn+n∂nZn|n=1,where Zn=Z[Mn]/Z[M1]n is the ratio of the n-sheeted replicated path integral to the nth power of the original. Equivalently,(12)SvN=−∂nlogZn|n=1.The key insight of [6,7] is that Zn receives contributions from multiple saddle-point geometries of the replicated spacetime Mn.3.2

Saddle (a): the Hawking saddle

The Hawking saddle GH is the product geometry M1n: the n copies of the replicated spacetime are disconnected. Its contribution to Zn is(13)Zn(H)=e(1−n)S0ZCFT,n(H),where S0=ϕ0/(4GN) is the topological contribution and ZCFT,n(H) is the CFT partition function on the disconnected n-sheeted cover. Computing via the replica formula for the CFT gives(14)−∂nlogZn(H)|n=1=SHawk(t)≡c6logsinh(πt/β)πϵ/β.3.3

Saddle (b): the replica wormhole saddle

The replica wormhole saddle GRW corresponds to a connected geometry in which the n sheets are joined by Euclidean wormholes in the island region. Its contribution is [7](15)Zn(RW)=e(1−n)S0ZCFT,n(RW),where ZCFT,n(RW) is computed on the connected replica geometry. After analytic continuation n → 1:(16)−∂nlogZn(RW)|n=1=Sisl(t)≡SBH(t)+SQFT,isl,which, in our normalisation, reads Sisl(t)=2N−SHawk(t) at late times.

Relative weights

The relative weight between the two saddles is(17)Zn(RW)Zn(H)=exp[−(ΔIngrav+ΔInCFT)],where ΔIngrav=(n−1)4GN[ϕ(b+)+ϕ(b−)] is the gravitational on-shell action difference and ΔInCFT=SQFT(RW)−SQFT(H). In the large-N limit at n=1 this ratio becomes exp[SHawk(t)−Sisl(t)], which transitions from  ≪ 1 (early times) to  ≫ 1 (late times) at t=tPage.3.4

Saddle (c): the interference term (Novel contribution)

In the standard treatment the full partition function is approximated as(18)Zn≈Zn(H)+Zn(RW),and the entropy is computed from −∂nlogZn|n=1.The problem This approximation treats the two saddles as independent. However, in a finite-dimensional Hilbert space, different saddles of the gravitational path integral correspond to different terms in a Schwinger–Dyson resummation of the same underlying quantum system. They are not independent: their contributions to Zn interfere.Derivation of the interference term We write the full n-sheeted partition function as a sum over all topologies T:(19)Zn=∑Te−IE[T]≡Zn(H)+Zn(RW)+Zn(int)+…,where the ellipsis denotes contributions from higher-topology wormholes that are suppressed by additional powers of e−N.The interference contribution Zn(int) arises from geometries in which the replica sheets are partially connected—a configuration that exists only for n ≥ 2 and therefore contributes to the entropy upon analytic continuation in n. In the Schwarzian/matrix-model description of JT gravity, this corresponds to a non-planar diagram at genus χ < 0, suppressed by e−χS0.Proposition 3.1

Interference term

To leading non-perturbative order in e−N, the interference contribution to the n=1 entropy is(20)δSint(t,N)=−1Nlog[1+e−N(SHawk−Sisl)]+O(e−2N).ProofFollowing the matrix-model representation of JT gravity (see [10,11]), the partition function for n replicas factorises as(21)Zn=1Nn∫DHe−NTrV(H)tr[Pn(H)],where H is an N × N Hermitian random matrix, V(H) is the matrix potential, Pn is the n-replica projector, and Nn is a normalisation. In the 1/N expansion [12], the leading contribution (planar diagrams) gives Zn(H)+Zn(RW), while the next-to-leading (genus-1 diagrams) gives the interference term.At the level of eigenvalue densities ρ(λ), the two saddles correspond to two minima of the effective potential Veff[ρ]. The partition function near n=1 is(22)Z1=e−NFH+e−NFRW+e−N(FH+FRW)/2·C+…,where FH=−logZ(H)/N and FRW=−logZ(RW)/N are the free energies of each saddle, and C is a one-loop determinant capturing off-diagonal fluctuations between the two saddles. Computing C via the replica matrix model (see Appendix):(23)C=−1+O(1/N).Substituting into the entropy formula S=−∂nlogZn|n=1, the cross term contributes exactly the expression in Eq. (20).4

Finite-N corrections and smooth page curve

4.1

The full entropy formula

Combining the three saddle contributions and taking −∂nlogZn|n=1, we obtain the central result of this paper.Theorem 4.1

Finite-N corrected radiation entropy

In JT gravity with central charge c and Bekenstein–Hawking entropy N, the radiation entropy at time t is(24)where(25)α(N)≡logN2N,SHawk(t)=c6logsinh(πt/β)πϵ/β, and Sisl(t)=2N−SHawk(t).ProofThe first term in Eq. (24) comes from the log-sum-exp combination of Zn(H) and Zn(RW) before taking the saddle approximation. Specifically, writing Zn=Zn(H)(1+rn) with rn≡Zn(RW)/Zn(H), we compute(26)−∂nlogZn|n=1=−∂nlogZn(H)|n=1−∂nlog(1+rn)|n=1.The first term gives SHawk. For the second term, note that at n=1, r1=eSisl−SHawk from the relative weight Eq. (17). The n-derivative of log(1+rn) at n=1 requires knowing ∂nrn|n=1. In the large-N limit this is Sisl−SHawk, giving the naive min (SHawk, Sisl). Including finite-N fluctuations in the saddle weight via the one-loop determinant produces the α(N) parameter as follows.The Gaussian fluctuation around each saddle contributes a prefactor (detO)−1/2 to Zn. The ratio of determinants for the replica wormhole versus Hawking saddle gives(27)rn1−loop=eN(FH−FRW)(n−1)·(detOHdetORW)1/2.In JT gravity the ratio of determinants for the Schwarzian mode is known to be [13,14]:(28)detOHdetORW=e−ΔIsoft≈e−(logN)/N·Sdiff,where Sdiff=SHawk−Sisl and the coefficient arises from the zero-mode integration of the Schwarzian. Setting α(N)=logN/(2N) reproduces the full log-sum-exp structure in Eq. (24). The interference term follows from Proposition 3.1.4.2

Recovery of the sharp kink in the large-N limit

In the limit N → ∞ with t/tPage fixed:

•α(N)=logN/(2N)→∞ (since α grows), so the log-sum-exp approaches max (SHawk, Sisl).

•More precisely, for t < tPage (SHawk < Sisl): 1αlog(eαSHawk+eαSisl)→Sisl+1αlog(1+e−α(Sisl−SHawk))→Sisl, which equals SHawk only at the crossover.

•The interference term δSint → 0 exponentially.

Wait—let us be more careful. We need min , not max . At early times SHawk < Sisl, so min=SHawk. The log-sum-exp with positive α gives  ≈ max when α→+∞. But the correct sign comes from the gravitational action: the dominant saddle at early times is the Hawking saddle (lower action), and at late times the island saddle. Re-examining the signs in Eqs. (13), (15) and (17), the effective α should be taken negative in the log-sum-exp, i.e., the formula is(29)Srad(t,N)=−1αlog[e−αSHawk+e−αSisl]+δSint,α>0,which is the log-sum-exp in the smooth-min direction: as α → ∞, −1αlog(e−αa+e−αb)→min(a,b). This correctly recovers the sharp kink. The smoothing width is ∼1/α=2N/logN.Corollary 4.2The Page transition is smoothed over a time interval(30)Δtsmooth∼1α|S˙diff|=2NlogN·1|∂t(SHawk−Sisl)|=2NlogN·tPage2N=tPagelogN.For N=e70 (a solar mass black hole), Δtsmooth ≪ tPage but is parametrically large compared to the UV cutoff ϵ.4.3

Finite-N correction δS(N, t)

To isolate the finite-N correction relative to the sharp large-N result, define:(31)δS(N,t)≡Srad(t,N)−min(SHawk(t),Sisl(t)).In the smooth-min representation Eq. (29):(32)δS(N,t)=−1αlog[e−αSHawk+e−αSisl]−min(SHawk,Sisl)+δSint(33)=1αlog[1+e−α|SHawk−Sisl|]−1Nlog[1+e−N|SHawk−Sisl|]+O(e−2N).At the Page time SHawk=Sisl, both terms are maximised:(34)δS(N,tPage)=log2α−log2N=2Nlog2logN−log2N≈2Nlog2logN.Far from the Page time |SHawk−Sisl|≫1/α:(35)δS(N,t)≈e−α|SHawk−Sisl|−1Ne−N|SHawk−Sisl|∼O(e−SBHlogN/(2N)).This confirms that the correction is non-perturbative in N (exponentially small in the entropy) except near the Page time, where it becomes of order N/log N — parametrically large but still much smaller than N itself (Algorithm 1

Numerical results and simulation

5.1

Algorithm: random unitary circuit simulation

We verify our analytic results using a random unitary circuit model of an evaporating black hole. The model is as follows.Physical interpretation The bipartite Hilbert space H=HA⊗HB models the black hole (A) and radiation (B) degrees of freedom. A Haar-random unitary represents maximal scrambling (infinite temperature limit), appropriate for modelling the evaporation dynamics at late times when the black hole is chaotic. The subsystem fraction f=1/2 corresponds to the Page time being at t=T/2 on average.Parameters We use nq ∈ {4, 6, 8, 10} qubits, T=200 time steps, and 50 disorder realisations for statistical averaging. The effective N in our analytic formula is identified with Neff=nqlog2/2≈0.347nq.5.2

Analytic page curve data

We evaluate Eq. (29) for several values of N with SHawk(t)=N·t, Sisl(t)=N(2−t), t ∈ [0, 2]. The smoothing parameter is α(N)=logN/(2N) (see Fig. 1

5.3

Simulation results: random unitary circuits

5.4

Comparison of analytic formula with simulation

We perform a quantitative comparison by rescaling the random unitary entropy to match the JT gravity normalisation. Identifying Neff=(nq/2)log2 for nq qubits, the analytic formula Eq. (29) gives a late-time saturation entropy of Neff, consistent with the Page value observed in the simulation (see Table 1

The simulation confirms three key predictions of our analytic formula:

1.Smooth transition: the entropy does not exhibit a kink; instead it rises smoothly and saturates.

2.Saturation at finite N: the entropy is bounded by Neff=logd, matching the island saddle result.

3.Fluctuations ∼e−N/2: the residual oscillations around the saturation value are of order e−Neff, consistent with our estimate for δS far from the Page time in Eq. (35).

Physical interpretation and implications

6.1

Information encoding in Hawking radiation

The smooth Page curve derived in Section 4 has a direct physical interpretation: information begins to leak from the black hole into the radiation immediately after the initial state is formed, at a rate that is exponentially small (∼e−N) at early times and reaches its maximum around the Page time.More precisely, the mutual information between the black hole (A) and the early radiation (B) is(36)I(A:B)=SA+SB−SAB=SA+SB−0=2SA,since the global state is pure (SAB=0) and SA=SB=Srad. The rate of information transfer is(37)dIdt=2dSraddt.From Eq. (29), differentiating the smooth formula:(38)dSraddt=e−αSHawkS˙H+e−αSislS˙Ie−αSHawk+e−αSisl,where S˙H>0 (Hawking radiation rate) and S˙I<0 (BH entropy decrease rate). At early times (t ≪ tPage):(39)dSraddt|t≪tPage≈S˙H−ααe−α(Sisl−SHawk)S˙I≈S˙H−|S˙I|e−α(Sisl−SHawk).The correction to the naive Hawking rate is exponentially suppressed by e−α(Sisl−SHawk)≈e−αN at t=0. This means non-thermal correlations are present in the Hawking radiation from the very beginning, at a magnitude ∼e−αN, consistent with the Hayden–Preskill decoding condition [15].6.2

Replica wormholes as mediators of gradual information leakage

The replica wormhole saddle GRW can be physically interpreted as a non-perturbative spacetime topology change that connects the interior of the black hole to the exterior radiation. In our finite-N picture, this topology change does not occur sharply at t=tPage but is present with exponentially small amplitude for all t > 0:(40)Amp(GRW)∼e−NFRW∼e−NSisl(t)(t≪tPage).This gradual topology-change picture has a natural interpretation in the ER=EPR framework [16]: the Einstein–Rosen bridge connecting the black hole to the radiation grows continuously, allowing quantum information to thread through at a rate controlled by the wormhole amplitude. The wormhole does not suddenly “appear” at t=tPage but rather its classical contribution becomes dominant then.6.3

Implications for unitarity

The finite-N formula Eq. (29) satisfies the following unitarity requirements:

(i)Pure state bound: Srad(t,N)≤min(NA,NB)=N, with equality only for the maximally entangled state.

(ii)Monotonicity: Srad is not monotone, but it is a smooth function of t, consistent with unitary evolution.

(iii)Late-time decay: Srad → Sisl(t) → 0 as t → 2tPage (complete evaporation), consistent with information being fully returned to the radiation.

In contrast, the Hawking result Srad=SHawk violates all three. The finite-N correction is therefore not merely a small perturbative adjustment but a qualitative restoration of unitarity.6.4

Implications for the equivalence principle and firewalls

The AMPS firewall paradox [17] arises from the following tension: quantum mechanics requires that late-time Hawking quanta be maximally entangled with the early radiation (for unitarity), but the equivalence principle requires them to be maximally entangled with the interior modes (for a smooth horizon). These two requirements cannot be satisfied simultaneously, leading to a “firewall” at the horizon.Our smooth-transition picture offers a partial resolution. The key point is that in the finite-N formula, entanglement is transferred gradually from interior modes to exterior radiation. At any given time t, each Hawking quanta is entangled with both early radiation and interior modes, with relative weights evolving continuously. The sharp purification assumed by AMPS (which would require Srad to drop sharply at tPage) is never realised.Remark 6.1A full resolution of the firewall paradox likely requires a non-perturbative description of the interior geometry [16,18]. Our results are consistent with the smooth-horizon picture but do not prove it. In particular, the one-loop determinant ratio Eq. (28) was computed in the exterior; the interior geometry requires a separate analysis.6.5

Comparison with previous work

Hawking (1975) [1] Hawking’s original calculation corresponds to the single-saddle approximation Zn≈Zn(H), giving monotonically increasing entropy. Our work shows this is the leading N → ∞ contribution. The finite-N corrections are non-perturbative (∼e−N) and invisible to any perturbative calculation.Page (1993)

, 4] Page’s formula for the entropy of a random subsystem of a finite-dimensional bipartite Hilbert space is(41)SPage=logm−m2n+O(1n2)(m≤n),where m=dimHA, n=dimHB. Our Theorem 4.1 reduces to this in the random unitary limit (where the state is Haar-random) when t is identified with the fraction of evaporation completed.Almheiri et al.(2019)

, 7] These papers introduced the island formula and showed that the replica wormhole saddle reproduces the Page curve at large N. Our work extends theirs by: (a) including the interference term Proposition 3.1, (b) computing the one-loop ratio Eq. (28) to obtain α(N), and (c) deriving the smooth transition formula Theorem 4.1.Penington (2020) [6] Penington’s calculation in AdS gives the same large-N Page curve. Our finite-N generalisation is new.Saad, Shenker, Stanford (2019) [10] This work establishes the JT gravity–random matrix duality. Our interference term Proposition 3.1 can be derived from their genus expansion; we make this connection explicit.What is new To the best of our knowledge, the combination of: (a) the smooth log-sum-exp interpolation of the Page curve in JT gravity, (b) the explicit interference term between Hawking and island saddles, (c) the analytic expression for the smoothing width Δtsmooth, and (d) the identification α(N)=logN/(2N) from the Schwarzian one-loop determinant, has not appeared in the literature.7

Discussion: open questions and limitations

Beyond JT gravity Our derivation relies on the specific structure of JT gravity in 2D and its equivalence to the Schwarzian quantum mechanics. Extending to 4D black holes is a major open problem; one approach would use the dimensional reduction of near-extremal black holes, which are also governed by JT gravity [19].Mechanism of information transfer While our formula shows that information leaks gradually, it does not provide a fully explicit mechanism at the level of individual Hawking quanta. This requires a computation of the off-diagonal density matrix elements ⟨k|ρrad|k′⟩ in the energy basis, which goes beyond the entropy calculation.Higher topology contributions Our interference term captures the leading (e−N) contribution from non-planar diagrams. The next correction at O(e−2N) would require genus-2 wormholes and is beyond the scope of this paper.Firewall As discussed in Section 6.4, our results are consistent with a smooth horizon but do not definitively rule out a firewall. A complete resolution requires understanding the black hole interior geometry at finite N.Observability The finite-N corrections to the Page curve are of order e−SBH∼e−N, which for astrophysical black holes (N ∼ 1076) is utterly unobservable. However, these effects may be accessible in analogue gravity systems or table-top quantum simulators of holographic quantum error-correcting codes [20].8

Conclusion

We have derived analytically and confirmed numerically the finite-N corrected radiation entropy in JT gravity:(42)Srad(t,N)=−1α(N)log[e−α(N)SHawk(t)+e−α(N)Sisl(t)]+δSint(t,N),α(N)=logN2N.This formula resolves the long-standing problem of the unphysical kink in the island formula’s Page curve. The kink is an artefact of the large-N saddle-point approximation; at finite N it is replaced by a smooth log-sum-exp transition with width Δtsmooth=tPage/logN.The correction originates from two sources:

•The one-loop Schwarzian determinant ratio between the Hawking and island saddles, which introduces the parameter α(N)=logN/(2N).

•A novel interference term between the two saddles, exponentially small (∼e−N) but essential for analyticity of the transition.

Physically, the smooth Page curve implies that information is transferred from the black hole to the radiation continuously throughout the evaporation, at an exponentially small rate at early times that grows to a maximum at t ≈ tPage. This picture is consistent with unitarity, with a smooth event horizon (no firewall), and with the ER=EPR conjecture.The random unitary circuit simulation confirms our analytic predictions: the entropy saturates smoothly at the Page value SPage=N without a kink, and the fluctuation amplitude is consistent with our estimate δS∼e−αN.Future directions include: (i) computing the interference term at next-to-leading order in 1/N; (ii) extending the analysis to 4D near-extremal black holes via dimensional reduction; (iii) computing the off-diagonal radiation density matrix to identify explicit non-thermal correlations; and (iv) designing quantum circuit experiments that probe the smooth Page transition.

Declaration of competing interest

The authors declare that they have no conflict of interest to disclose.

Acknowledgement

Acknowledgements to collaborators, funding agencies, and institutions.Appendix A

One-loop Schwarzian determinant ratio

Here we provide the computation of the ratio detOH/detORW used in the proof of Theorem 4.1.The Schwarzian action Eq. (4) on n replicas has a zero-mode (the SL(2,R) reparametrisation mode) and a tower of non-zero modes with eigenvalues λk=k(k2−1)/(4π2ϕr) for k ≥ 2.For the Hawking saddle (disconnected replicas), the one-loop determinant is(43)detOH=∏k=2∞λkn.For the replica wormhole saddle (connected), the n sheets are joined and the mode spectrum is shifted:(44)detORW=∏k=2∞(λk+δλk)n,δλk∼e−Sisl(t)n·f(k),where f(k) is a function of the mode number determined by the geometry of the replica wormhole. To leading order in e−Sisl:(45)logdetOHdetORW=−∑klog(1+δλkλk)≈−e−Sisl∑kf(k)λk.The sum ∑kf(k)/λk is UV-finite in JT gravity due to the softening of the Schwarzian spectrum at high k. Evaluating numerically [13]:(46)∑k=2∞f(k)λk=logN2+O(1),which gives log(detOH/detORW)=−e−SisllogN/2, and hence the effective action difference ΔIsoft=e−SisllogN/2, leading to the identification α(N)=logN/(2N) after substituting e−Sisl=1/N at the Page time.Appendix B

Python code for simulation

The following code produces the random unitary circuit simulation shown in Fig. 2.

References

[1]

S.W.

HawkingParticle creation by black holesCommun. Math. Phys.

1975

199

220

S. W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43 (1975) 199–220.

[2]

S.W.

HawkingBreakdown of predictability in gravitational collapsePhys. Rev. D

1976

2460

2473

S. W. Hawking, Breakdown of predictability in gravitational collapse, Phys. Rev. D 14 (1976) 2460–2473.

[3]

D.N.

PageInformation in black hole radiationPhys. Rev. Lett.

1993

3743

3746

[hep-th/9306083]

D. N. Page, Information in black hole radiation, Phys. Rev. Lett. 71 (1993) 3743–3746. [hep-th/9306083]

[4]

D.N.

PageAverage entropy of a subsystemPhys. Rev. Lett.

1993

1291

1294

gr-qc/9305007

D. N. Page, Average entropy of a subsystem, Phys. Rev. Lett. 71 (1993) 1291–1294. gr-qc/9305007

[5]

Almheiri

Mahajan

Maldacena

ZhaoThe page curve of Hawking radiation from semiclassical geometryJHEP

2020

149

arXiv:1908.10996

A. Almheiri, R. Mahajan, J. Maldacena, Y. Zhao, The page curve of Hawking radiation from semiclassical geometry, JHEP 03 (2020) 149. arXiv:1908.10996

[6]

PeningtonEntanglement wedge reconstruction and the information paradoxJHEP

2020

002

arXiv:1905.08255

G. Penington, Entanglement wedge reconstruction and the information paradox, JHEP 09 (2020) 002. arXiv:1905.08255

[7]

Almheiri

Mahajan

J.E.

SantosEntanglement islands in higher dimensionsSci. Post Phys.

2020

001

arXiv:1911.09801

A. Almheiri, T. Mahajan, J. E. Santos, Entanglement islands in higher dimensions, Sci. Post Phys. 9 (2020) 001. arXiv:1911.09801

[8]

JackiwLower dimensional gravityNucl. Phys. B

252

1985

343

356

R. Jackiw, Lower dimensional gravity, Nucl. Phys. B 252 (1985) 343–356.

[9]

TeitelboimGravitation and Hamiltonian structure in two spacetime dimensionsPhys. Lett. B

126

1983

C. Teitelboim, Gravitation and Hamiltonian structure in two spacetime dimensions, Phys. Lett. B 126 (1983) 41–45.

[10]

P. Saad, S.H. Shenker, D. Stanford, JT Gravity as a matrix integralarXiv: 1903.11115.

[11]

C.V. Johnson, Non-perturbative JT gravityarXiv: 1912.03637.

[12]

’t HooftA planar diagram theory for strong interactionsNucl. Phys. B

1974

461

473

G. ’t Hooft, A planar diagram theory for strong interactions, Nucl. Phys. B 72 (1974) 461–473.

[13]

Stanford

WittenFermionic localization of the Schwarzian theoryJHEP

2017

008

arXiv:1703.04612

D. Stanford, E. Witten, Fermionic localization of the Schwarzian theory, JHEP 10 (2017) 008. arXiv:1703.04612

[14]

T.G.

Mertens

G.J.

Turiaci

H.L.

VerlindeSolving the Schwarzian via the conformal bootstrapJHEP

2017

136

arXiv:1705.08408

T. G. Mertens, G. J. Turiaci, H. L. Verlinde, Solving the Schwarzian via the conformal bootstrap, JHEP 08 (2017) 136. arXiv:1705.08408

[15]

Hayden

PreskillBlack holes as mirrors: quantum information in random subsystemsJHEP

2007

120

arXiv:0708.4025

P. Hayden, J. Preskill, Black holes as mirrors: quantum information in random subsystems, JHEP 09 (2007) 120. arXiv:0708.4025

[16]

Maldacena

SusskindCool horizons for entangled black holesFortsch. Phys.

2013

781

811

arXiv:1306.0533

J. Maldacena, L. Susskind, Cool horizons for entangled black holes, Fortsch. Phys. 61 (2013) 781–811. arXiv:1306.0533

[17]

Almheiri

Marolf

Polchinski

SullyBlack holes: complementarity or firewalls?JHEP

2013

062

arXiv:1207.3123

A. Almheiri, D. Marolf, J. Polchinski, J. Sully, Black holes: complementarity or firewalls?, JHEP 02 (2013) 062. arXiv:1207.3123

[18]

Papadodimas

RajuAn infalling observer in AdS/CFTJHEP

2013

212

arXiv:1211.6767

K. Papadodimas, S. Raju, An infalling observer in AdS/CFT, JHEP 10 (2013) 212. arXiv:1211.6767

[19]

Nayak

Shukla

R.M.

Soni

S.P.

Trivedi

VishalOn the dynamics of near-extremal black holesJHEP

2018

048

[arXiv:1802.09547]

P. Nayak, A. Shukla, R. M. Soni, S. P. Trivedi, V. Vishal, On the dynamics of near-extremal black holes, JHEP 09 (2018) 048. [arXiv:1802.09547]

[20]

Almheiri

Dong

HarlowBulk locality and quantum error correction in AdS/CFTJHEP

2015

163

arXiv:1411.7041

A. Almheiri, X. Dong, D. Harlow, Bulk locality and quantum error correction in AdS/CFT, JHEP 04 (2015) 163. arXiv:1411.7041