{'fraction_non_alphanumeric': 0.0428904011461318, 'fraction_numerical': 0.021639207258834767, 'mean_word_length': 4.286367939413064, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Self-similar accretion in thin disks around near-extremal black holes

Near-maximally spinning black holes display conformal symmetry in their near-horizon region, which is therefore the locus of critical phenomena. In this paper, we revisit the Novikov-Thorne accretion thin disk model and find a new self-similar radiation-dominated solution in the extremely high spin regime. Motivated by the self-consistency of the model, we require that matter flows at the sound speed at the innermost stable circular orbit (ISCO). We observe that, when the disk pressure is dominated by radiation at the ISCO, which occurs for the best-fitting Novikov-Thorne model of GRS 1915+105, the Shakura-Sunyaev viscosity parameter can be expressed in terms of the spin, mass accretion rate and radiative efficiency. We quantitatively describe how the exact thin disk solution approaches the self-similar solution in the vicinity of the ISCO and for increasing spins.

One important boundary condition imposed in the original Novikov-Thorne model is that the torque vanishes at the ISCO, which leads to an inconsistency of the model (the fluid velocity diverges at the ISCO even though observables remain finite). This hypothesis has also been challenged by external considerations (Krolik (1999); Gammie (1999); Li (2000); Noble et al. (2010)). Modelling physical boundary conditions at the ISCO is crucial in particular for highly spinning black holes in order to accurately calibrate their spin estimate (Li et al. (2005)). In Penna et al. (2012), a self-consistent boundary condition hereafter called "sonic-ISCO" was proposed, which consists in equating the radial comoving fluid velocity at the ISCO with the sound speed. Two arguments were presented in its favour. First, in slim disk models around the Schwarzschild black hole analysed in Abramowicz et al. (2010), the sonic point asymptotes to the ISCO in the thin disk regime (where the accretion rate is sub-Eddington), independently of α. Secondly, the boundary condition implies that advection terms are negligible with respect to the stresses at the ISCO in the energy balance, in the limit where the disk height is negligible with respect to α. We will check that the latter assumption is self-consistent for typical stellar-mass black holes and luminous AGN parameters. A crucial consistency condition of thin disk is therefore satisfied. As we will discuss, the boundary condition also implies that the specific internal energy is negligible, which automatically enforces another hypothesis of the Novikov-Thorne model. In addition, we will show that the sonic-ISCO boundary condition allows for a well-defined near-horizon near-extremal scaling behavior, in contrast to the no-torque boundary condition. Our first objective is to derive the complete thin disk model with the sonic-ISCO boundary condition from first principles. We took particular care to correct various algebraic errors and typos found in the literature.
A new feature emerges from our analysis. For highly spinning stellar-mass black holes, gas pressure at the ISCO becomes negligible with respect to radiation pressure. The sonic-ISCO boundary condition then implies that the accretion rate is not a free parameter of the model. Instead, the disk height at the ISCO is the fourth independent parameter. It allows us to infer a best-fitting value for the Shakura-Sunyaev parameter α in terms of the total luminosity, radiative efficiency and spin of the source. For example in the case of the X-ray binary source GRS 1915+105, recent estimates of the mass and spin, using the Very Long Baseline Array, give M ∼ 12 ± 2 M and J/M 2 = 0.98 ± 0.01 (Reid et al. (2014)). Also, the best-fitting model of the NuSTAR observation of GRS 1915+105 in the plateau state gives a disk luminosity at 23% ± 4% of the Eddington rate (Miller et al. (2013)). Taking into account the torque contribution to the radiative efficiency of thin accretion disks (Li (2002)), we will infer that the best-fitting Novikov-Thorne model is a disk dominated by radiation around the ISCO with α = 0.43 and M = 1.6 × 10 18 g sec −1 .
It is well established that radiation-dominated phases in thin disks are both thermally and viscously unstable (Lightman & Eardley (1974); Shakura & Sunyaev (1976)). Without a stabilization mechanism, such disk models are not expected to occur as steady flows in nature or be visible in numerical simulations (see the thermal instability in Mishra et al. (2016)). Instead, the radiation pressure instability has been argued to drive the variability of radiation-dominated objects such as GRS 1915+105 (Janiuk et al. (2000)). The thermal and viscous instabilities are the main open issues to be understood in the Novikov-Thorne model. In this paper, we will further explore the radiation-dominated phase with the sonic-ISCO boundary condition in the extremely high spin regime, where new features arise.
According to the cosmic censorship conjecture, the spin of the Kerr black hole is bounded by its value at extremality. Extreme spinning black holes admit a near-horizon region with global conformal symmetry SO(2, 1) (Bardeen & Horowitz (1999)) and probe fields can therefore be described using the language of conformal representation theory (Porfyriadis & Strominger (2014) (2016)) or critical phenomena (Gralla et al. (2016b)). Approaching the extremality bound by realistic accreting processes is, however, limited by the absorption cross-section of retrograde photons as shown by Thorne (1974), leading to the bound J/M 2 < 99.8%. This bound is commonly accepted (see e.g. recent work by Kesden et al. (2010)), even though it can be lowered by magnetic fields (Gammie et al. (2004)) or instead challenged (Sadowski et al. (2011)). More fundamentally, for accretion disk models where the inner edge approaches the marginally bound orbit instead of the innermost stable orbit, the capture of retrograde photons asymptotes to zero, which allows spinning black hole to further approach extremality (Abramowicz & Lasota (1980)). Such scenarios have been shown to occur in slim disk models (Sadowski et al. (2011)). It is therefore worthwhile to explore which specific signatures may arise in the extremely high spin regime, where conformal symmetry is only slightly broken. This programme has already been carried out for gravitational wave signatures (Hadar et al. (2014(Hadar et al. ( , 2015; Gralla et al. (2016a)) and we initiate this programme for accretion disks in the simplest case of thin disk models. (for earlier related work see also Grumiller & Piso (2009)).
The paper is organized as follows. In Sec. 2, we first show how the emergent conformal symmetry constraints finite observables in the near-extremal near-horizon limit of Kerr spacetime, which completes the analysis of Gralla et al. (2016b). Sec. 3 is devoted to revisit the thin disk equations and the main assumptions of the Novikov-Thorne model. We impose the sonic-ISCO boundary condition (Penna et al. (2012)), and we explain its physical implications. In Sec. 4, we analyse the features of the Novikov-Thorne model and we construct piece-wisely the global solution. We focus on parameters for two notable cases of interest: stellar-mass and supermassive black holes. Moreover, we study the phase diagram of thin accretion disks in either cases in terms of the spin parameter. Remarkably, thin accretion disks reach a critical spin beyond which a phase transition occurs, close to the ISCO, from gas pressure dominated to radiation pressure dominated. This happens in the high spin regime. As an example, we consider the Novikov-Thorne model parameters to be those of the best-fitting values of GRS1915+105. Finally, in Sec. 5, we present a first explicit example of critical (self-similar) disk in the near-horizon region of near-extremal Kerr. In Appendix A, the reader can find all the necessary properties of the Kerr spacetime and the main definitions we used throughout the paper. Appendix B contains the local solutions to the thin disk equations.
Since it is obtained as a scaling limit of Kerr, it obeys the vacuum Einstein's equations and it admits enhanced scaling symmetry along the Killing vector H 0 ≡ T ∂ T − R∂ R . It turns out that the exact symmetry is further enhanced to the conformal group times the abelian group of azimuthal rotations, SO(2, 1) × U(1).
Conformal symmetry also appears for infinitesimal deviations from extremality (Amsel et al. (2009);Bredberg et al. (2010)), which implies that near-extremal black holes admit critical phenomena in their near-horizon region. In order to make that statement precise in the presence of matter accretion, it is fundamental to control the location of the ISCO. 1 It has been already noticed in the original paper of Bardeen et al. (1972) (see their Fig. 2) that the BL coordinate radius of the ISCO (r 0 ), as well as photon and marginally bound orbits, asymptotically merges with the BL coordinate radius of the horizon r + towards extremality, Therefore, in order to obtain a finite coordinate location for the ISCO in the near-horizon limit, it is necessary to scale the near-extremality parameter as σ =σλ 3/2 , whereσ is fixed, as described in Gralla et al. (2015). The λ → 0 limit still gives the NHEK spacetime in Eq. (2.3) and the ISCO is now located at the finite radial position R = R 0 ≡ 2 1/3σ2/3 . The near-horizon near-extremal limit above is invariant under the scaling λ → cλ for arbitrary c. Therefore, all physical quantities, which are finite in this limit, will be invariant under a simultaneous rescaling of R → cR, T → T/c and R 0 → cR 0 . In other words, finite observables in the above limit are zero eigenfunctions of the total dilatation operator T ∂ T − R∂ R − R 0 ∂ R 0 . In particular, stationary observables O ought to be of the form (R 0 /R) h , where h is the critical exponent, which equals the conformal weight of the observable in the near-horizon region (2.5) Thanks to an additional length-scale, namely R 0 , near extremal configurations allow a richer variety of critical exponents as compared with the extremal case discussed in Gralla et al. (2016b).
1 Local conformal symmetry also arises (Guica et al. (2009)) and plays an important role in understanding the quantum physics of black holes (see the reviews Bredberg et al. (2011);Compère (2017)). Conformal symmetry also appears for specific probes away from extremality (Castro et al. (2010)). We will, however, not use these additional structures for the present considerations.
In Sec. 5, we will provide an explicit example of such a critical behavior within the Novikov-Thorne accretion disk model with critical exponents Here, u R is the near-horizon radial velocity, h the disk opening angle, F the vertical energy flux, T the temperature, p the pressure and ρ the rest-mass density.

THIN DISK EQUATIONS
In this section, we revisit the accretion thin disk model around Kerr spacetime, originally presented in Novikov & Thorne (1973) and Page & Thorne (1974), but instead of imposing the "no-torque" boundary condition at the ISCO, we will impose the "sonic-ISCO" boundary condition introduced in Penna et al. (2012). The familiar reader might skip sub-sections 3.1-3.3 and directly reach sub-section 3.4, where new features arise.

Fundamental equations
We assume a single component relativistic viscous fluid, flowing along the four-velocity u µ in the fixed Kerr background. The stress-tensor can then be split uniquely with respect to u µ as where ρ is the rest-mass density, Π the specific internal energy, p the total isotropic pressure, h µν = u µ u ν + g µν the projector into the local rest frame (LRF), S µν the symmetric, transverse and traceless anisotropic stress-tensor that takes into account viscous stresses, and q µ the transverse energy flux. All these quantities are relative to the LRF and all indices are raised with g µν . The disk is assumed not to selfirradiate. Magnetic fields are ignored except for their contribution to the viscous stresses. Neutrinos and dark matter are ignored. The fundamental equations governing the dynamics of the accretion disk are respectively, the rest-mass conservation law, the relativistic Navier-Stokes equations and the energy conservation equation. These fundamental equations must be supplemented by the equation of state, the radiative energy transport law and prescriptions about the nature of the viscous effects and opacity. Conservation of the rest-mass along the fluid fourvelocity is valid for energies much below 2m b , where m b is the rest-mass of the baryon species, or equivalently, for temperatures much below 10 10 K for electrons and 10 13 K for protons/neutrons. We will see that the last hypothesis is self-consistent after presenting the solution to the model. Assuming thermal equilibrium and the existence of an equation of state of the form Π = Π(p, v) where v = 1/ρ is the specific volume of the fluid, one can define the temperature T and specific entropy s by the first law T ds = dΠ + pdv. The energy conservation equation might then be written in terms of local entropy production (Ellis (2009) ;α is the shear tensor and a µ = u ν u µ ;ν is the acceleration. The four-entropy flux is defined as S µ = sρu µ + q µ /T and the local form of the second law of thermodynamics S µ ;µ ≥ 0 implies the two constitutive equations with η, λ ≥ 0 being the viscosity and heat conduction coefficients, respectively.

Thin disk approximation
The fundamental equations (3.2) require additional simplifying assumptions in order to construct analytical accretion disk models. The original thin disk model provides a set of working assumptions, carefully listed in Page & Thorne (1974), to transform the full system of partial differential equations to an algebraic non-linear system of equations and obtain local solutions in analytical form.
In what follows, we will use cylindrical BL coordinates (t, r, z, φ) where z = r cos(θ) (see Appendix A for details and notation). We assume stationarity and axisymmetry. The half-thickness H of the disk is defined as the vertical distance between its upper surface and the equatorial plane at z = 0. A geometrical thin disk has an opening angle h(r) ≡ H/r 1. We suppose that all radiation is vertical, q µ ∼ δ µ z . We are only interested in vertically integrated quantities between z = −H and z = +H. For example the surface density of the disk is defined as (3.5) In order to directly obtain the equations for vertically integrated quantities from Eqs. (3.2) as the leading order of a Taylor expansion at the equator, one can simply ignore all z dependence of all physical quantities, except for the pressure p and radiation flux q z , which can be assumed to be (see discussion in Abramowicz et al. (1997)): The function F is then the radiation flux emitted from the upper or lower side of the disk. The orbital motion of the fluid is taken to follow nearly circular equatorial geodesics with a small radial (non-geodesic) component u r produced by viscous stresses and responsible for accretion on to the black hole.
Since the shear tensor σμν only admits the nonvanishing component σrφ = σφr < 0 in the local rest frame (see Eq. (A16)), and given the first constitutive equation (3.4), one can define the vertically integrated shear tensor as (3.7) We assume the α-viscosity prescription of Shakura & Sunyaev (1973) where α is a free parameter and p is total pressure. 2 Two further assumptions were used in the original analysis of Novikov & Thorne (1973), namely negligible specific internal energy density Π = 0 and negligible advection. We will see that these hypotheses are obeyed as a result of the sonic-ISCO boundary conditions and h α, which will be obeyed in turn by explicit check on the solutions.

Thin disk equations
The dynamics in the radial direction is governed by the restmass, energy and angular momentum conservation laws. The radial Navier-Stokes equation is trivially satisfied with the additional assumptions that velocity and pressure gradients are negligible. Standard manipulations of Eqs. (3.2) lead to is the conservation of rest mass, Eq. (3.9b) conservation of energy and Eq. (3.9c) is a combination of energy and angular momentum conservation laws. The constant of integration M is the accretion rate, while E, L, Ω and σrφ are kinematic quantities of circular equatorial geodesics (see Appendix A1 for their expressions). The right-hand side of Eq. (3.9c) is nothing else than M P defined in (A10). The dimensionless integration constant P 0 is fixed by the boundary conditions, which are discussed in sub-section 3.4. The vertical Navier-Stokes equation describes the pressure balance along the cylindrical vertical coordinate z and reads as as derived in Abramowicz et al. (1997) 3 . The total pressure is the sum of the radiation pressure and the gas pressure p = p (gas) + p (r ad) , (3.11a) where b = 4σ SB /c is the radiation constant density, k B Boltzmann's constant, σ SB Stefan-Boltzmann's constant and m p the rest-mass of the proton. We impose the energy transport law bT 4 =κΣF (3.12) 2 Some specific MHD turbulent disks might be modelled by this prescription (Balbus & Papaloizou (1999)). The modification of αp to αp g a s was argued in Lightman & Eardley (1974) and further developed in Bisnovatyi- Kogan & Blinnikov (1977) and Sakimoto & Coroniti (1981). Such models, later called βp models, do not suffer from thermal instabilities. However, MHD simulations do not conclude on their validity ( whereκ is the optical opacity of the disk originating from free-free (ff) absorption and electron scattering (es).

Sonic-ISCO boundary condition
A physical system is determined by its equations of motion and its boundary conditions. In the original analysis of Novikov & Thorne (1973), it was assumed that there is no torque at the ISCO (located at r = r 0 ), which is equivalent to assume that there is no radiation at that point, which is also equivalent to fixing the integration constant P 0 = 0 in Eq. (3.9c). Following Penna et al. (2012), we impose instead that the (purely radial) fluid velocity in the frame corotating with the geodesic flow equals (minus) the sound speed at the ISCO, c s (r 0 ) = −ur (r 0 ).
( 3.15) This boundary condition has three important advantages that we will describe in what follows. So far, we can rewrite the energy conservation equation as where Q diss = −S µν σ µν is the dissipation function, Q cool = q µ ;µ is the cooling function and Q adv = ρTu µ s ,µ is the advection function that takes into account the rate of change of the specific entropy along the four-velocity. Now, as shown in Penna et al. (2012), the boundary condition (3.15) implies the following scaling relations: Q adv ∼ h 2 , Q cool ∼ Q diss ∼ αh; so if h α, advection can indeed be neglected. Thin disk models usually assume that the specific internal energy density is negligible Π = 0. This hypothesis is justified if the sound speed is nonrelativistic, c s 1. Due to the gravitational potential, the sound speed is highest in the near-horizon region of the disk, where we will impose the boundary condition (3.15). The hypothesis c s 1 will therefore be obeyed as long as |ur | 1 at the ISCO, which is already part of the hypotheses since we assumed that the fluid follows nearly circular geodesics. We will check that the solutions indeed obey |ur | 1 and c s 1. Finally, with the no-torque boundary condition, the disk model has no regular limit at the ISCO as we will review below, while the boundary condition (3.15) regularizes the model.
Let us now fix the remaining integration constant P 0 for the sonic-ISCO boundary condition. It will be fixed as a function of the disk height at the ISCO, h 0 , as follows. The integral in the right-hand side of (3.9c) is zero when evaluated at the ISCO. Thus, Eq. (3.9c) reads as ( 3.17) On the other hand, dividing the first two equations (3.9a) and (3.9b), we get In the second step, we have expressed the radial velocity component in the LRF, u r = (∆ 1/2 /r)ur , we have substituted the definitions of Σ and W, we have used the Shakura-Sunyaev prescription and assumed the sonic-ISCO boundary condition (3.15), ur 0 = − p 0 /ρ 0 . Equating (3.17) and (3.18), we find after some algebra involving the functions defined in Appendix A. Our final formula for M P 0 disagrees with the constant C derived in Penna et al. (2012), which is easily seen to be incorrect since it has the wrong dimension of length. The physical meaning of the integration constant P 0 is to introduce a torque at the ISCO. More precisely, the torque might be derived by comparison of Eq. (3.9c) and Eq. (12) of Li (2002). The torque g 0 is then The total energy radiated per unit time as measured by an observer at infinity is therefore given by both accretion and torque contribution (Li (2002)) where η 0 = 1−E 0 is the specific conserved energy of a particle orbiting along the ISCO (A9) and η is the radiative efficiency of the disk Here, we are ignoring the capture of radiation from the hole, which decreases the efficiency for high spin (Thorne (1974)).

FEATURES OF THE GENERAL SOLUTION
We introduce the dimensionless mass and mass accretion rate where M is the mass of the Sun. The global solution can be approximated by a piece-wise construction of three local solutions described in Appendix B, which are patched according to their range of validity. The qualitative features of the global solution depend upon the region that dominates at the ISCO. There are three possibilities depending which of the three relevant local solution is valid around the ISCO.

Gas-pressure-dominated ISCO
In the "standard" first two cases, the ISCO lies in the region dominated by gas pressure, either Region In Fig. 1, the first region in which the ISCO is located is called the edge region. There, the gas pressure overwhelms the radiation pressure and the opacity due to electron scattering is dominant over the free-free absorption. The disk height is given by Eq. (4.2). The transition to the inner region occurs when radiation pressure starts becoming predominant over the gas pressure. A new transition occurs when the gas pressure becomes dominant again over radi-ation pressure and the resulting region is called the middle region. For very low spins, the inner region is absent and the edge and middle regions merge. The outer region is still gas-pressure-dominated, but the main mechanism responsible for the opacity is the free-free absorption. We checked that the height satisfies h α = 0.2 for r 10 13 M where the model breaks down for other reasons (because self-gravitation is not negligible). The temperature of the accretion disk is hotter in the region near the ISCO, and in general the temperature is higher for faster spin. By evaluating the temperature at the ISCO for a = 0.999M, we found T ∼ 10 7 K. This implies that the assumption of conservation of mass is valid. We also checked that the sound speed is negligible with respect to light speed, c s = p/ρ c, in the entire disk.
The original Novikov-Thorne model did not contain the edge region, but it should contain it by consistency: for standard spins, the radiation pressure in the inner region goes to zero at the ISCO and therefore the gas pressure needs to dominate at low enough radius. Yet, if one assumes the no-torque boundary condition, the Novikov-Thorne model is singular at the ISCO in the edge region. The reason is readily seen because in the [Gas-es] solution (B2) the function P defined in (A10) vanishes at the ISCO if P 0 = 0, and the radial velocity is then divergent at the ISCO. The non-zero torque introduced by the sonic-ISCO boundary condition allows us to regulate the model as claimed earlier.
In Fig. 2, the disk is always dominated by gas pressure, but the dominant contribution to opacity varies with radius. In the edge region, where the ISCO lies, the free-free absorptions are dominant. The disk height is therefore given by Eq. (4.3). The middle region is dominated by electron scattering. The outer region is again dominated by free-free absorption.
The assumption h α = 0.2 is obeyed for r 10 21 M where the model breaks down. The accretion disk for supermassive black holes is colder with respect to the stellar-mass black holes; it never exceeds T ∼ 10 4 K and the conservation of mass is obeyed. We also checked that the sound speed is negligible with respect to light speed, c s c, in the entire disk.

Radiation-pressure-dominated ISCO
Let us now discuss the configurations where the ISCO lies in the region dominated by radiation pressure and electron scattering (Region [Rad-es] of Appendix B). This scenario happens for very high spins as we will discuss below. A new feature arises as a result of the sonic-ISCO boundary condition: a constraint relates the accretion rate M, the mass, spin and α parameter, while the opening angle at the ISCO is unconstrained.
Indeed, the disk height (B4d) in geometric units is given by At the ISCO, P is given by P 0 (see Eq. (A10)) that can be evaluated using the expression for the sonic-ISCO boundary condition in Eq. (3.19). Therefore, h 0 appears linearly in both sides and we are left with a constraint among the parameters of the model given by The function f a M monotonically decreases and vanishes at extremality.
It is instructive to compare the accretion rate in Eq. (4.5) with the Eddington accretion rate M E dd = 4πM/(κ es η), where η is the radiative efficiency defined in Eq. (3.22). We find the reduced accretion rate (4.6) Since the α parameter is usually hard to estimate, it is useful to solve the relation (4.6) for α in terms of m using Eq. (3.22). The accretion rate is then determined from Eq. (4.5) and we obtain The free parameters of the model where the ISCO lies in a radiation-dominated region can be finally taken to be (M, a, m, h 0 ).
In the phase diagrams of both typical stellar-mass and supermassive black holes displayed in Figs 1 and 2, the gas pressure dominates at the ISCO for all standard spins. However, for sufficiently high spins, radiation pressure dominates as we will now show. Let us first discuss configurations where the ISCO lies in the Region [Gas-es] for standard spins such as the case studied in Fig. 1. If one (wrongly) assumes that the ISCO lies in Region [Gas-es] for very high spins, one deduces from Eqs. (B3a) and (4.2) that the ratio of pressures at the ISCO is where in the last step we took the near-extremal scaling a/M = √ 1 − σ 2 . For σ 1, one obtains that radiation pressure will instead dominate. In the example of M = M = 1 and α = 0.2, the transition occurs (in the sense that p r ad = p gas ) at a/M = 0.99996 which is much above the Thorne bound of 0.998 (Thorne (1974)) and therefore much probably unrealistic.
If instead the ISCO lies in the Region [Gas-ff] for standard spins such as the case studied in Fig. 2 which is 1 for σ 1. Again, radiation pressure dominates for sufficiently high spins and a new near-ISCO region opens up. However, for typical parameters M = 10 7 , M = 10 5 and α = 0.2, the transition to the near region occurs at a/M = 1 − 10 −18 which is unreasonably high to be realistic.
However, there are more interesting parameters to consider. Let us take an accretion rate at 23% Eddington ( m = 0.23) as a model for the plateau state of GRS 1915+105 (Miller et al. (2013)) with the spin estimate J/M 2 = 0.98 (Reid et al. (2014)). Assuming a radiation-dominated ISCO, we can derive the α parameter using Eq. (4.7) after fix-  ing an estimate for h 0 . We checked that for any value 0 < h 0 < 0.01, the resulting values of α = 0.43 and M = 16.5 differ by 1% or less. The continuous transition between the gas-pressure-and radiation-dominated phases occurs at a/M = 0.980 ± 0.001. For definiteness, we choose h 0 = 0.002, so that the transition exactly occurs at a/M = 0.98.
We conclude that GRS 1915+105 could be modelled by a high spin a/M > 0.98 and fastly accreting thin disk, which is radiation-dominated at the ISCO. As an example, we depict in Fig. 3 the phase diagram for the parameters M = 4, m = 0.23, h 0 = 0.002 and spins 0.98 < a/M < 1. The inner region is radiation-dominated. At finite radius away from the ISCO, there is a transition to a middle region, which is gas pressure-dominated, but whose opacity is still dominated by electron scattering. There is also an outer region further away, where the main mechanism for opacity is freefree absorption.

NEAR-HORIZON NEAR-EXTREMAL SOLUTION
In the near-extremal regime, a new spacetime region, the NHEK region, opens up as explained in Section 2, which is characterized by conformal symmetry. Since the radiationdominated solution is the only one relevant in the limit of extremely high spins, we will perform the near-extremal nearhorizon limit of the solution (B4) in order to exhibit its conformal properties. We perform the scaling (2.2) with a/M = √ 1 −σ 2 λ 3 and let λ → 0. The ISCO is located at R 0 = 2 1/3σ2/3 . The accretion rate M (in terms of asymptotic time) is constrained as in Eq. (4.5) which implies the scaling M ∼ λ. In terms of near-horizon time T, the accretion rate M is finite, Therefore, contrary to the asymptotically flat observer, the NHEK observer measures a finite non-zero accretion rate M . The ratio M /R 0 is fixed by the parameters of the model. The near-horizon behaviour of the solution (B4) can be obtained by performing the near-extremal near-horizon scaling and trading M for M . We get the following expressions at leading order (O(λ 0 )): All these quantities are defined for R ≥ R 0 . Note that all quantities, except p and T, depend on the ratio which is independent of the choice of the constantσ. The only exception is the radial component of the four-velocity, which has an additional power of R 0 , and therefore depends on the position of the ISCO. The radial scaling exponents reproduce the table of conformal weights in Eq. (2.6), announced in Section 2. Note one unusual property of the selfsimilar solution: starting from h 0 at the ISCO the disk height h decreases with the radius; it increases again outside of the range of validity of the self-similar solution.
The self-similar solution is an approximate solution of the disk around the ISCO. It is important to discuss its range of validity. The critical temperature does not depend upon the accretion rate M or the height of the disk h 0 . Since α or M are overall factors of the temperature profile (B4h), the relative error between the actual temperature profile and the constant critical temperature only depends upon the radius and the spin. We find that for near-extremal spins 0.96 ≤ a/M ≤ 1, the actual temperature profile deviates from the critical temperature by less than 25% only in the range r 0 ≤ r ≤ 2.1−2.2M, where the upper bound is nearly independent of the spin. This very limited near-ISCO region is the region where the disk is approximately described by the self-similar solution. The region is biggest when the ISCO approaches M, which occurs closest to extremality. If a higher precision    is required, the region of validity shrinks accordingly. We plot in Fig. 4 the range of validity of the temperature of the self-similar solution with 25%, 15% and 10% relative precision.
Other physical quantities can be analysed similarly. Unfortunately, the relative precision of the pressure requires a spin higher than the Thorne bound a/M = 0.998 and a narrower region around the ISCO, as plotted in Fig. 5. Another important physical quantity is the radiation flux F. Fig. 6 shows that the range of validity of the self-similar solution is limited to a narrow region around the ISCO. The physical relevance of the self-similar solution (5.2) is therefore uncertain.
In conclusion, we found a first example of critical accretion around a near-extremal black hole, which admits a scale invariance in a region close to the ISCO. More elaborated models are required to find more realistic solutions.   x ≥ x + = 1 + (1 − a 2 s ) 1/2 1/2 describes the Kerr spacetime outside the event horizon. The ISCO is located at Following the tradition, we define the dimensionless functions 4 : All functions (A6) go to unity far from the black hole x 1. The functions A, B and R are monotonically decreasing, while the functions C, D, F and G are monotonically increasing in their domains of definition. The radial coordinates where C = 0 are given by x 2 = 2 cos 1 3 arccos(a s ) + π 3 , x 3 = −2 cos 1 3 arccos(a s ) , with x 1 ≥ x 2 ≥ x 3 . The last photon orbit is given by x = x 1 and is placed between the ISCO and the event horizon, x 0 ≥ x 1 ≥ x + . The function D is positive outside the event horizon.
Here, we list the kinematic quantities characterizing the geodesic motion: respectively, the time component u t , the azimuthal component u φ , the conserved specific energy E, the conserved specific angular momentum L, the corotating Keplerian angular velocity Ω with respect to a stationary observer and a conserved quantity along the geodesic motion L 2 , which appears in the thin disk vertical equation.
It is useful to introduce the function P defined by Here, P 0 is a real constant to be fixed by imposing a boundary condition, as explained in Sec. 3.4. In the above sum, we identify x i = x i+3 . The function Q of Page & Thorne (1974) can be expressed as Q = x −1 BC −1/2 P but we find more economical and natural to use P.

A2 Observer frames
We define three privileged observers.