Microscopic State of BHs and an Exact One Body Method for Binary Dynamics in General Relativity

In gravitational collapses, the horizon and singularity's realisation in the finite future of the proper time used co-moving observer happens in the future of infinitely far away future of the normal time used outside probe. To the latter the horizon and singularity of a black hole formed through gravitational collapse are physical realities only in the sense of uncertainty principle and ensemble interpretation. We provide two exact time dependent solution families to the Einstein equation and show that they form a pair of complementarity description for the microscopic state of black holes by showing that the Bekenstein-Hawking entropy formula follows properly from their canonical wave function's degeneracy. We also develop an eXact One Body method for general relativity two-body dynamics whose conservative part requires no perturbative input from post newtonian approximation and applies to the full three stages of black hole binary merger events. By this method, we analytically calculate the gravitational wave forms following from such merger processes. In the case black holes carry exact and apriori horizon and singularity our wave forms agree with those following from conventional effective one body method but exhibit more consistent late time behaviour. In the case the black holes carry only asymptotic horizon and extended inner structure thus experiencing banana shape deformation as the merger progresses, our wave forms exhibit all features especially the late time quasi-normal mode type oscillation seen in real observations.

IX.An exact one body method 22 X.The full gravitational wave forms 25

XI. Conclusion 28 I. INTRODUCTION
The existence of Bekenstein-Hawking entropy [1][2][3][4][5][6] is a strong hint that black holes are inner-structured [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24].The current observational data [25][26][27] on gravitational waves arising from Black Hole Binary (BHB) merger events may brought us information about such innerstructures already.The purpose of this work is to provide a thorough understanding of such inner-structures in the framework of standard general relativity and canonical quantum mechanics.At the same time an exact one body method will be developed to calculate the gravitational wave form of BHB merger process so that the effects of their inner-structures on observational signals can be exhibited directly.Both these two points seem to be impossible purposes.We will show that the first impossible arises from our wrong taking of the horizon and singularity as physical realities independent of probe and physical interactions.The second impossible arises from our wrong choosing of the effective background geometry which is aimed to account for the relative motion between the two bodies as static.This choice not only makes exact solutions to the problem impossible, but also the known effective one body method inconsistent at the principle level.
By physical reality, we refer to measurement outputs of any single type of physical quantity.In contrasts, we will use physical law to refer to relations bridging at least two types of physical realities.Physical reality is invariant under general coordinate transformation but observer dependent.Physical law is not only general coordinate invariant but also observer independent.Taking Einstein Equation R µν − 1 2 g µν R = 8πGT µν as an example, both sides of it are physical realities.We will give concrete examples to show how the physical law's observer inde-arXiv:2311.11764v2[gr-qc] 11 Dec 2023 pendence is confused with the physical reality's observer dependence thus leading to the wrong belief of horizon and singularity's probe independence in section VII.But to understand the analysis there firmly we need five sections to build the right description for the microscopic state of black holes.In a sense, all sections of this work are original and innovative.However, for those who only want to know why our two purposes are possible, they can jump to section VII and IX directly.We also have a short version of this long paper which is available on email requsting.The following is a highlight of the key point of each sections of this work.
However, just as any existence question in philosophies do, the form of matter's existence in physics is also related with the definition of time.In general relativity, this is especially so.Because of time delay effects caused by gravitation, the outside fixed position observers cannot see the horizon's formation in a to be black hole collapsing star in any finite future quantifiable with their physical time definition.These observers are not required to be human being or any bio-wisdom.They can be any outside probe or detectors with fixed position or approximately fixed position.For example, either participant involved in a merger event of BHB system works as probes of its inspiral partner.These probes will not see the horizon's formation in their partners' evolution in any finite future quantifiable with their physical time definition.So to them, the form of matter's existence inside their merger partners is just that of the frozen star [60][61][62][63][64].But the reason for being frozen here is the infinite gravitational time delay due to the horizon to form, instead of any exclusive force such as that arising from the neutron degeneration pressure or any other unknown physics [65][66][67].We will provide exact metric and Penrose-Carter diagram description for this frozen star in section II of this work.
By the freely falling observer's time definition, materials consisting of a to be black hole collapsing star would indeed contract into the mathematical surface defined by r = 2GM in the finite future.However, both the materials being observed and the observer themselves experience not any exotics as that surface is being walked through, all of them are freely falling.As local observer and observable, the spacetime region around them are completely flatten.The non-zero curvature is physically relevant only to the global viewpoint holden investigators.This is also the case even when the r = 0 point is concerned.The observer and observables coming from the east world of the central point just normally walk over that point and get to the west world afterwards.Tidal forces affect only non-local or finite sized observer and observables.However, general relativity is a time-reverse invariant theory.No matter how strongly is an extended object tidally stretched as it walking to the r = 0 point, its shape and inner structure will be retrieved when it walking away from that point just as the world is time reversed [68][69][70][71][72]. Penrose and Hawking's singularity theorem [28][29][30][31] only requires that all matter contents get into the horizon onto the singularity (an equal time hyper surface) in finite proper time.It does not prohibit them from walking across that epoch.Exact metric description and PC-diagram representation for this picture will be provided in section III.
With section V supplemented as a working example, section IV discusses the quantisation of inner-structures described by the exact metric families of sections II and III.We will decompose the matter content of spherically symmetric black holes into many concentric shells characterised by each shell's mass and radial quantum number and write the wave functional of the whole system as direct product of all compositing shells.By considering all possible shell decomposition schemes allowed by the total mass summation rule and the radial quantum number assigning schemes allowed by the asymptotic horizon condition, we analytically show that the degeneracy of the system's wave functional is exactly that required by Bekenstein Hawking entropy formula, except for some logarithmic type correction.Numerically determining this correction is possible but challenging as the mass of the black hole becomes large.In contrast with the euclidean path integration method based on saddle point approximation [73][74][75] and string theory calculation [76][77][78][79][80][81][82][83][84][85][86][87][88][89], our microscopic state definition and number counting for black holes are purely general relativity based and resort not any hyper physical concepts such as supersymmetry or extra spatial dimension.But our black holes will exhibit similar observational features as string theory fuzz balls [90,91] if their shadow images are taken [92,93] or gravitational wave echos [94][95][96][97][98][99] are measured.
The black hole complementarity principle [100][101][102][103][104] believes that both the outside fixed position observer and the freely falling observer's description of a collapsing star's evolution are equally valid.We will elaborate in section VI that, for the outside fixed position observers, the uncertainty principle or quantum fluctuation's ubiquity implies that parallel universe or ensemble description of the initial status of materials consisting of the black hole is necessary.While for the freely falling co-moving observers, this necessity can be interpreted as the ergodicity of the system's evolution on the proper time axis.So when the initial status' uncertainty or the evolution dynamics' ergodicity is considered, not only are the two types of observer's description both valid, but also are they complete.The fire wall paradoxes [105][106][107] will be analysed in this section to show that the truly rational resolution to the information missing puzzle [108][109][110][111][112][113][114][115][116][117][118] is to give up the popular vacuum fluctuation and partial escaping mechanism for Hawking radiation instead of the idea of black hole complementarity.This section mainly consists of logic reasoning instead of mathematical derivation.So we put it in the middle part of the paper to avoid preaching.Mathematical derivations underlying can be found in our earlier works [16,17].
Through three concrete examples, section VII will analyse how the belief that gravitational collapse will cause the formation of observer independent horizon and singularity makes the wrong of (A) taking the horizon and singularity of global viewpoint holders as the horizon and singularity detectible to local viewpoint holders and (B) taking the physics law's general coordinate invariance as the physical reality's observer or probe independence.The physical reality's observer or probe independence is not the fact at all.Our three examples are, (i) the frequency of light signals emitted from the surface of a collapsing star and measured by an outside fixed position observer and a freely falling observer; (ii) the pressure of a dust type test mass shell freely falling towards the horizon of a pre-existing Schwarzschild black hole measured by a far away fixed position observer and a freely falling observer co-moving with the shell itself and (iii) the horizon of two cosmological sized super giants separated double hubble radius away in our real universe.This section is the key part of this work but we burry it here because it is conceptually so bold that some people may think it non-sense at all.Section VIII will show that when the horizon is taken as a probe dependent physical reality, the black holes involved in the binary merger event will experience banana shape deformation under the inhomogeneous back reaction force of gravitational wave radiation.We will illustrate such a deformation's occurrence in the newtonian gravitation theory and argue that the physical picture would be similar in the fully general relativistic treatments of the system's evolution.We then provide a static mechanical analysis for the banana shape deformation of black holes and show that the radiation activity of such BHBs will be suppressed by a factor of sin 4GM z/a(t) 4GM z/a(t) when they inspiral and merge.Here z is the system's stretching mulriple relative to the diameter of the undeformed Schwarzschild horizon.z is time dependent whose concrete form reflects the black hole's deformation sensitivity or susceptibility under the back reaction force of gravitational wave radiation.
Section IX will propose an exact one body (XOB) method to calculate the relative motion orbit of general relativity two-body systems.This method will avoid the self-contradicting ingredients of conventional effective one body (EOB) method [119][120][121].The conservative part of XOB hamiltonian involves no post-newtonian approximation [122][123][124][125][126] as inputs.Taking the most simple quadrupole radiation of gravitational waves as the source of dissipation, this method will allow us to trace the relative motion of binary systems even when their evolution develops to the ring-down stage.By this method we will calculate the relative motion orbit and gravitational wave forms for binary systems consisting of standard Schwarzschild black holes to make comparisons with the results of conventional EOB method.Our results will coincide with that of EOB at early times but behaves more rationally as the very late time stage is arrived.By looking the black hole as an arc of fixed length, this method will yield gravitational wave forms highly similar with those obtained in numeric relativity method [148][149][150][151][152][153][154][155] and black hole perturbation theories [156][157][158][159][160][161].
Section X will compare the relative motion orbits and gravitational wave forms following from our XOB method for three different time variation modes z(t) of banana shape deformation.These three modes describe black holes with three typical shape deformation sensitivity or susceptibilities under the inhomogeneous back reaction force of gravitational wave radiation.By exhibiting the resultant gravitational wave forms' difference, we will show that the inner structure and the shape deformation sensitivity of black holes are observable experimentally instead of purely theoretical arguments.By assuming that gravitational wave carries no energy away, we will derive a universal upper bound for the late time quasinormal mode's frequency as functions of the symmetric mass ratio of the binary system irrespective of their inner structure details.While by setting the conservative part of the XOB hamiltonian to the minimal possible value, we will get a lower bound for such quasi normal frequency.We will use existing data to test this two bounds quantitatively.
Section XI is the conclusion of the whole paper.

II. THE OUTSIDE FIXED POSITION OBSERVER'S DESCRIPTION
The physical reality's probe dependence is a very widely accepted concept in high energy physics.For example, the hadron structures measured in the deep inelastic scattering experiments depends on the probe particle's momentum transfer.One purpose of this work is to show that, the horizon and singularity's observer independent is just a counter part of this concept in general relativity and black hole physics.
By the word of outside fixed position observer, we refer to all out side probe or detectors that are affected by the gravitation emanate from the to be black hole collapsing stars.These observers can only probe or detect physics related with the collapsing stars during their life time, i.e. existence duration of the probe themselves.Typical examples for these observers include, either member of a BHB system which are experiencing merger process, and/or all imaging light signals that are emitted from the background accretion disk and enter the detector of the event horizon telescope on earth.These observer or probe's positions are not truly fixed, but during their whole existence duration and in the spherical coordinate system of the being probed black holes, their radial positions are only slowly varying functions of the time coordinate which is understood as physic times by the infinitely far away investigators.
To these observers, the general form of spacetime metric of the collapsing star can be written as The function form of B is determined by the normalisation of matter sources' four velocity We choose the scaling of t so that hdt 2 = dτ 2 , where τ is the proper time measured by observers fixed on the collapsing matter, up to a regular r-dependent scale factor.In the case of gravity dominating over all other interactions, neglecting the pressure is a rational doing.In this case, the Einstein equation where we used short notations m = GM [t, r].Equation (II.4) follows from the condition that G 0 0 G 1 1 − G 0 1 G 1 0 = 0. Equation (II.5) follows from the condition of G θ θ = G ϕ ϕ = 0 supplemented by (II.4).
Eq(II.4) can be written into the form of a first order differential equation and integrated formally y ′ (t, r) − p(t, r)y(t, r) = q(t, r), y ≡ ṁ2 (II.6) where c(t) is an arbitrarily tunable function of tcoordinate, while r 1 and r 2 are two arbitrary reference point.Taking m/r m ′ −1/2 as a slowly varying function of r, the integration p(t, z)dz in (II.8) can be done approximately so that This means that given initial mass function m(0, r), we can always tune the form of c(t) so that on the boundary of the matter occupation region ṁ m ′ | m.o.r bndry = dr dt | m.o.r bndry matches with the geodesic motion of test particles freely released in the Schwarzschild metric caused by the total mass of the collapsing star.But according to equation (II.4), the initial value of ṁ2 (0, r) is completely determined by the form of m(0, r) thus allows no artificial tuning to implement a given ṁ m ′ (0, r m.o.r bndry ).This is because we neglect all pressure effects as well as the external force/potentials exerted on the matter contents so that their initial speed equals zero.In our numerical examples, we will simply set c(t) = 0 and take r 2 as a small parameter to assure integral convergence and accuracy.In this case equation (II.10) implies that as m → r 2 , ṁ2 m ′2 ∼ (r−2m) → 0. In this limit, equation (II.4) leads to This explains the origin of asymptotic expression (II.3), which is in fact a result independent with respect to the matter source's equation of state.The Ricci scalar of spacetime metric (II.1)-(II.2) when simplified using (II.4) and (II.5) has the following form This is regular as t → ∞, because the asymptotic behaviour of (II.11).Divergent happens only when we try to go beyond the Future of infinite Far Future (FiFF), that is, when we use alternative time coordinate which covers range outside that t does.Eqs(II.4)-(II.5)allows us to express all the higher order time derivatives of m in terms of its spatial derivatives, e.g.
These expressions tell us that, as t → ∞ the n-th order time derivative of m are proportional to (r−2m) −n ∼ ∞.This makes numeric integration of equations (II.4) or (II.8) more and more difficult as we approach the configuration m(t, r) ≈ r 2 .However, the basic feature of → 0 is very robust.This causes the outside probes not being able to probe m(t, r)' evolution into profiles of r − 2m < 0 in any finite t time.
Physically this is just the fact that, as outside observers or probes, they cannot see the formation of horizons in any finite t durations.Of course, if the proper time description is adopted, the co-moving observers can see matter's collapse into the horizon in finite durations, just as was illustrated in references [12].However, since those observers or detectors are freely falling, they will experience not any exotics as they pass through the horizon surface.Consider two simple initial mass function (i) m(0, r < 1) = 1 4 r 3 , (ii) m(0, r < 1) = 1 4 r 1 and m(0, 1 < r) = 1.The former has central sparser initial mass density profile ρ ∼ (1 − r 2 4 ) −1 , the latter has a central denser and singular profile ρ ∼ r −1 but no horizon apriori.For these two initial mass functions the t = 0 integration (II.8) can be done exactly as well as analytically The physical black holes form through gravitational collapses are defined by their asymptotic horizon.The asymptotic horizon is only mathematical surface which is infinitely approachable but never reached in any finite physic time.It provides us a boundary to look the physical black holes as statistic systems and introduce parameters such as internal energy and entropy to describe their macroscopic features.Just as all normal statistic systems do, the microscopic state of the physical black holes can be simply attributed to the mechanical motion of their matter contents.At here, this motion is the freely falling of concentric matter shells consisting of the collapsing star under randomly specified initial radius and radial speed.The randomness of initial conditions arise from the uncertainty principle.We will discuss the quantisation of this mechanical motion and show that its wave functionals are indeed degenerated in the number of ways required by the area law formula of Bekenstein Hawking entropy.But before this doing, let us consider this motion's description from the inside co-moving observer's aspect.

III. THE INSIDE CO-MOVING OBSERVER'S DESCRIPTION
By the term inside co-moving observer, we refer to all local detectors or probes which are freely falling together with the matters consisting of the collapsing star and use proper time as time coordinate to describe physic evolutions accessible to them.As local probes, these observers are not necessarily point particles.As long as their size is much smaller than the space-time region's scale they are detecting, things will be okay.Their detection are just interactions between them and the to be detected spacetime region itself.The interaction between different parts of the system is also a detection which has the potential of putting the microscopic state of the whole system onto an ergodic evolution orbit.To these observers, the most general metric ansatz for an isotropic but inhomogeneous collapsing star can be written as where as was done in the previous section, we have neglected the pressures so that the inner structures of the system are characterised by a single function M [ϱ] exclusively.This is reasonable when the gravitation dominates all non-gravitational exclusive interactions; τ and ϱ are so called Lemaitre coordinate and r s ≡ 2GM tot .By the conventional Boyer-Linderquist coordinate, ds 2 out can be equivalently written as The advantage of Lemaitre coordinate is that, we can more directly see the smoothness of the inside-outside metric's connection.The advantage of Boyer-Linderquist coordinate is that, we can more easily understand what we are describing is a simple neutral and spherical symmetry black hole in the standard general relativity.Since we neglect all exclusive pressures, the concrete form of a[τ, ϱ] can be obtained exactly from the dust sourced Einstein equation R µν − 1 2 g µν = 8πG{ρ, 0, 0, 0} Since a(τ, ϱ) is negative allowable, the sign of the coefficient of the dϱ 2 term in (III.1) must be defined in such a way that the whole term is positive definitely.
For each given ϱ, the behavior of a[τ, ϱ] involved in (III.1) or (III.4) is just the oscillation driven by a square inverse force or linearly inverse potential, The full period of this oscillation is During each first quater period, the oscillation equation can be integrated implicitly and the full oscillation function can be written as   4πϱ 2 dϱ = ϱ 1 .The matters' migration on the co-moving grid, marked by cyan arrows, is from the superficial to central region so that the co-moving ϱ-coordinate covers smaller and smaller regions.The right has a central-denser initial mass density ρ[0, ϱ] ∼ ϱ −1 .The matters migration on the co-moving grid in this case is from the central to the outer region, marked by cyan-arrows.
uniform initial density case, the evolution of the system is completely determined by the contraction and expansion of the co-moving grid as a whole.But in the other two cases, the co-moving grid's contraction and the matter's migration on the grid are asynchronous.In the central sparser case of the middle column, the matters' migration on the co-moving grid is from the superficial layers to the central ones, so that the matters are distributed on smaller and smaller co-moving grids.While in the central-denser case of right most column, the matters' migration on the co-moving grid is from the central region to the superficial layers, so that they cover larger and larger co-moving grids.In both latter cases, the matters' migration on the co-moving grid is towards a more and more homogeneous configure until they reach an exactly uniform distribution and over-cross each other.At the over-cross epoch, the physical mass density dM/(a 3 ϱ 2 dϱ) is divergent, but the density dM/ϱ 2 dϱ on the co-moving grid is finite.After that epoch, the mass-energy distribution will retrieve its initial configuration but with the anti-podal point interchanged.So the matters' gathering onto the central point and causing divergence there is not the terminal of their motion inside the horizon but a normal epoch of a periodical over-cross oscillation.Of course this oscillation happens outside the domain of time definition of the outside observers.To them this oscillation especially the matter shells' expansion away from the gravitational central happens only in a quantum parallel universe or statistic ensemble spacetime which is causally in-connected with the one they own live in and is understandable only as the results of the quantum fluctuation implied by the uncertainty principle.
Relative to other mono-directional collapsing star solutions to the sourceful Einstein equation such as Oppenhiemer-Snyder [133,134], Yodzis-Seifert-Muller [135,136] and Vaidya [137,138] in the history, the most important feature of our oscillation solutions family above is its not taking the dusts' collapse to a[τ, ϱ] = 0 and hitting on the singularity as the terminal of their motion inside the horizon.Instead it takes such event only as an intermediate epoch of the dusts' motion inside the horizon.The next step and following on motion is their passing through each other and periodic oscillation across the central point.Their motion over cross the central point is an analogue of two classic wave's superposition and independent propagation after encountering each other.This over cross oscillation mechanism is firstly noted and quantised with a functional Schrodinger equation to resolve the Schwarzschild singularity in reference [12] and then quantised with the method of shelldecomposition in references [13,14] to interpret the origin of Bekenstein-Hawking entropy and resolve the black hole information missing puzzle caused by hawking radiation.While in reference [15], exact solutions for the simple AdS 2+1 -Schwarzshild are provided and the corresponding Bekenstein-Hawking entropy formula is analytically proved by the integer number partition formula.
The existence of our solutions family (III.1)-(III.6)does not contradict with the singularity theorem of Penrose and Hawking [28][29][30], because all matters described in this solutions family indeed fall onto the central point and cause a spatial-like singularity in finite proper time as expected.To the co-moving observer, such a singularity is nothing but an equal-τ hyper-surface.Just like the horizon can be trivially traversed, the spatial-singularity which occupies only a future point on the co-moving observer's world line is also trivially traversable in almost no τ -time, after which all kinematic parameters of the system get their τ -time reversed value routinely.Tidal forces affects only the finite sized objects.However, just as we pointed out in the introduction section, such effects are time inverse symmetric.So even when the finite sized objects are involved in a physic collapsing process, the status information encoded in such objects will be retrieved as the singularity is traversed.This means that we can define the microscopic state of the black hole with the initial status of its consisting matters.For an exactly spherical symmetric Schwarzschild black holes, the only possible initial status information is its matter contents' layering structure and each layer's radial collapsing speed.
To the outside fixed position observer or physical probes, the asymptotic horizon provides an ideal excuse for them to say that the system has to be considered as an ensemble of many parallel universe, each contains a collapsing star with equal mass but different initial radial mass profile.To the freely falling co-moving observers, this excuse cannot be used.However, the fact that the domain of proper time τ 's definition covers that of the coordinate time t many times brings us another interpretation, ergodicity.That is, during the infinitely long τ -time evolution, the point representing the microscopic state of the black hole will walk through each point of the phase space that represents all possible microscopic states of the collapsing star.See FIG. 6 for references.To display the oscillations inside the horizon more conveniently, we extended the rule of conventional Penrose Carter diagram so that the east world and west world of the spacetime are simultaneously displayed in the {τ, ρ} plane.We introduced this diagram method firstly in reference [16].We will show in the following that the microscopic state definition based on the ergodicity of co-moving observers and that based on the uncertainty principle of the outside fixed position observers are equivalent to each other.
In both FIG. 3 and FIG.6, the horizon represented by red straight lines and the singularity represented by red zigzags are physical realities only to the investigators which can see the Future of infinitely Far away Future (FiFF), we will call the viewpoint of such observers as global of God-type.While the viewpoint adopted by the true physic probes will be called local.Both the outside fixed position observer and the freely falling co-moving observers are local.To them the horizon and singularity are only mathematical realities instead of physical ones, because they have not any physical ways to detect such realities' existence.

IV. QUANTISATION
As the definition function for black hole microscopic state, both the asymptotic M [∞, r] in (II.3) and the comoving mass profile M [ϱ] in (III.1)-(III.6)at the classic level are arbitrary thus un-countable.However, at the quantum level, things are rather different.We have two methods to quantise the system canonically.The first is to use the hamiltonian formulation of general relativity to build a Wheeler-Dewitt like functional Schrodinger equation which takes the instantaneous mass function m(t, r) as independent variable and the radial mass-energy density profile ρ[t, r] as eigenvalues.
where h ≡ 1 − 2m(r) r and m(r) is the value of m(t, r) at any given time t.The 3 + 1 decomposition or foliation of the spacetime itself defines the relation between m(t, r) and the instantaneous energy density ρ(t, r).
In the classic hamiltonian constraint m should be understood as the function of t and r simultaneously, while in its quantum version m is understood as the function of r only.This quantisation method is a generalisation of the minimal super space quantum cosmology.We firstly proposed this generalisation in reference [12] but with an error happens in the functional Schrodinger equation corrected in [16].In their counter parts of references [12,16], the zero term positions of (IV.1) and (IV.2) are occupied by terms of ρ 2 and 1 respectively.This difference is due to the time coordinate adopted in this work to define the hamiltonian constraint is the outside fixed position observers time t instead of the proper time τ of the inside co-moving observers.
To directly discretise the independent variable m(r) and resolve the functional Schrodinger equation (IV.1) is very difficult, see reference [12] for our earlier attempts.Beginning from reference [13] and through out references [14][15][16][17], we proposed a new discretisation stratagem.In this new stratagem, we decompose the whole collapsing star into many concentric shells {m i } and quantise each of them canonically.The shell decomposition here is nonapriori.All schemes allowed by the following two conditions would be considered, The former follows from the total mass sum rule, while the latter follows from the square integrability of m i 's wave function; tum number of shell m i 's wave function and classically they are determined by the shell's initial radius.Quantum mechanically the mass function of the whole collapsing star can be written as since ψ i (r) here is the eigenfunction of the hamiltonian whose conjugating time is that of the outside fixed position observers, this quantum mass function corresponds to the asymptotic function M [∞, r] of eq(II. 3).This can also be taken as the initial mass function on the comoving grid since the initial time epoch on the co-moving grid is arbitrarily choosable.Obviously, the mass function's quantum definition (IV.5) tells us that it is the combination of the shell decomposition and initial radius assigning scheme that defines the microscopic state of the collapsing star, neither one of them alone is enough to do that work.So our number counting of the microscopic state of the system should be made on the partition here is order-concerned, e.g.{m i } = {1, 0.7, 0.3} and {m ′ i } = {1, 0.3, 0.7} are two different partitions of M tot = 2. Two attentions must be payed with the number counting.The first is, the global maximal value position of each shell's wave function modular far away from it.r i h is the classic horizon defined by the masses the shell m i wraps."not too far" means that the global maximal position of the outmost shell's wave function should take the lowest value which lies outside the classic horizon r h = 2GM tot .The second is, any radial excitation scheme should not make the global maximal value of shell m i 's wave function |ψ 2 i | happens outside that of the shell m i+1 , for all i.The position of these maximal values corresponds to the classic asymptotic radius of the shells.This implies that the the radial excitation number {n i } cannot be set arbitrarily.It can be numerically shown that given shell mass decomposition {m i }, the number of possible {n i }-values' assigning schemes which satisfy these two requirements grows at most polynomially with the total mass M tot increases.This polynomial non-uniqueness is the origin of logarithmic corrections to the area law formula in our microscopic state definition for black hole.
From the viewpoint of the outside fixed position observers, the classic motion of each shell m i is freely falling in an effective Schwarzschild geometry determined by the total mass of shells inside m i (including m i itself) ) Using x and x µ (λ) ≡ {t(λ), x(λ), 0, 0} to denote the radius of shell m i and the world line of a representative point's moving on it, the standard geodesic equation and four velocity normalisation of the representative point tell us where Γ µν and γ i are the usual Christoffel symbol corresponding to the effective metric (IV.7) and an integration constant of the geodesic equations respectively.γ i = 0 corresponds that the shell is released from r ini = 2GM i and ṙini = 0; γ i = 1 corresponds that the shell is released from r ini = ∞ and ṙini = 0.The former implies that the surface r = r ini happens to be a horizon.The latter implies the masses inside r = r ini is far from forming horizon.From the two equations of (IV.9) and (IV.10), we can easily derive out Multiplying both sides of this equation with the shell mass m i and manipulate the resultant equation routinely, we will get This is nothing but the hamiltonian constraint of the shell's motion.Now replace the classic momentum m ẋ into its operator form iℏ∂ x and introduce a wave function ψ i (x) to denote the probability amplitude the shell m i be measured at radius being x status, we will get the quantum version of the Hamiltonian constraint (IV.12) This quantises the motion of shell m i canonically.All other shells can be quantised similarly.Directly multiplying their wave-functions together, we will get the wave functional of the whole collapsing star as follows Equation (IV.13) is almost the standard eigenstate Schrodinger equation with a coulomb like potential except the normalisation condition which involves no angular space factor 4πx 2 because ψ i (x) is the probability amplitude of the spherical shell as a whole to lie on the state of radius taking x value instead of the probability amplitude a point particle to be measured at position (x, θ, ϕ).The eigenfunction and eigenvalue of (IV.13) can be written down as (IV.17) where is the associated Lagurre polynomial; N i is the normalisation factor following from (IV.15).Equation (IV.17) explains the origin of our microscopic state number counting constraint (IV.4).
Although the radial quantisation condition (IV.17) or (IV.4) for each component shell provides countability for the microscopic state of the system, the number of state counted basing on equations (IV.3)-(IV.4) is still infinite.This is because both the two conditions do not impose discreteness on the masses {m i } of the component shells, but only on the quotients { GMimi (2−2γ 2 i ) 1/2 ℏ } and γ i are not required to be zero exactly.As a regularisation for this countable infinite, we introduce the following concept of partitions' distinguishability two partitions of equal number of shells {m i }& {m ′ i } are indistinguishable on precision ε if (IV.18) two partitions of unequal shell numbers are always distinguishable By this standard, we can quantitatively state that to precision ε = 1, the finest way of decomposing the collapsing star into concentric shells is that defined by the condition that for all i, GMimi (2−2γ 2 i ) 1/2 ℏ = 1 and γ i = 0. We will call this decomposition scheme as fundamental.By this fundamental scheme, the number of shells the collapsing star can be partitioned into is k For the fundamental shell decomposition above, to assure that the global maximal value of each shell m i 's wave function modular square happens outside the classic horizon r i h ≡ 2GM i , the lowest level radial quantum number assigning scheme is to set n i = 2 for all i.This is because in this case, the global maximal value of |ψ 2 ni (x)| happens on x If we set any of the n i s to 3, then the corresponding x ψn i max = 6.525GM i will go outside its nearest outer neighbour.If we set all of the n i s to 3, then the corresponding wave functional of the whole system would have a too fat non-zero tail outside the classic horizon r h = 2GM tot .The object described by the corresponding wave functional would have matter contents distributed in a spherical spatial region of radius r = 6.525GM tot instead of being localised inside the classic horizon.So we will not take such highly excited quantum state as the microscopic state of black holes defined by the classic horizon r h = 2GM tot .We will call the microscopic state defined by the fundamental shell decomposition and pure 2 radial quantum number as the fundamental state of the classic black hole.For latter usages, we introduce the following diagram to represent this fundamental state where k − ≡k−1 and GM i m i = √ 2 for all i.By this representation, the most general radial excitation scheme can be written as However, by the explicit form of each shell's wave function (IV.16), it can be easily proved that one cannot find any radial quantum numbers' combination {n i } which is different from that of the fundamental state (IV.21) and for all i, the position of |ψ 2 i |'s maximal value does not go outside that of the shell m i+1 's |ψ 2 i+1 |.This means that the fundamental state of the system is unique.
By choosing two or more inner shells and merging them into some new single shell, we have ways to excite the inner structure of the fundamental state without change the total ADM-mass/energy of the system.The most simple choice is selecting two neighbouring shells and merging them into a single one.For example, merging the outmost two shells and set ≈4 means that as the two shells are merged, n k − and n i can be set to 4 if we tune the value of γ k − or γ i to some value which is larger than their value in the fundamental state.This is rational because when we merge two shells, the resultant single shell must be put onto a state whose radial quantum number is roughly the summation of the merger-before members' quantum number to assure the global maximal value of the resultant shell's wave function happens outside the horizon r i ′ h = 2GM i ′ .Before the two shells' recombination, all n i equals to 2. The radial quantum number of shells inside the recombined shell can now be set to values greater than 2 because as the latter takes now a larger radial quantum number, the formers have the possibility to be excited to higher radial quantum number level and still keep their global maximal value inside the that of the latter.This means that a given recombined shell decomposition scheme could have more than one radial excitation ways.
The most general shell decomposition and radial excitation scheme can be generated from the fundamental state by shell recombination and radial quantum numbers' reassigning.The resultant microscopic state of the system will have equal total mass as but different inner profiles from the fundamental one.They can be represented as where n min ℓ takes the minimal possible value so that |ψ 2 ℓ |'s global maximal happens outside the horizon r h = 2GM tot .The number of shell recombination ways from the fundamental state is the number of ordered partition of integer k = GM 2 tot 2 √ 2 , which is 2 k−1 exactly, including the fundamental state itself.So the total number of microscopic state of the system can be written as where X(k) is a correction factor which arises from the non-uniqueness of the radial quantum number assigning scheme in the specified shell decomposition scheme.In the case of neglecting the requirement that |ψ 2 i |'s global maximal value happens inside that of |ψ 2 i+1 |, references [16] proves that X(k), which is named F (k) there, grows exponentially as k increases.However, when this requirement is considered, it can be verified that F (k)'s growing trend is at most polynomial way.See the next section for a concrete working example.
So by the distinguishability definition of (IV.18) and the result (IV.26), we can now claim that to precision ε, the number of all distinguishable ways of decomposing a collapsing star with asymptotic horizon into concentric shells {m i } and setting each shell's radial quantum number {n i } is where M tot is only required to be mildly larger than G − 1 2 ≡M pl .Since the decomposition precision parameter ε linearly affects the number of shells of the fundamental decomposition, c(ε) is a linear function of ε.This means that we can easily get the logarithmically corrected Bekenstein-Hawking formula with c(ε) = 1 by simply setting ε ≈ 1 8 √ 2π .To get the explicit formula of the logarithmic correction X(A), we need quantitatively work out the relation between the number of radial quantum number assigning scheme {n i } and the collapsing star's total mass M tot for specified shell decomposition {m i } under the constraint that for all i, the global maximal value of |ψ 2 i | happens not outside that of |ψ 2 i+1 |, and the global maximal value of all |ψ 2 i | happens outside the classic horizon of r i h = 2GM i , with n i−outmost takes the minimal allowed value.This may be a challenging work but operable in principle.
Equation (IV.17) or (IV.4) implies that the outmost shell has a minimal allowed mass GM tot (IV.29) Since the shell with γ outmost =0 has the maximal probability be found on the horizon and that with γ outmost = 1 has the maximal probability be found on the infinity r region.Exact horizon boundary condition implies that γ outmost =0.This means that This is also the typical mass of most of the outer layered shells.This explains why we have chances to get the area law entropy by counting the bulk motion degrees of freedom of the black hole matter contents.Because the number of such shells reads Note that for all macroscopic black holes, the value of m min outmost is far less than the standard model particles which consists of the collapsing star.So what the Bekenstein-Hawking entropy reflects is not the motion degrees of freedom of these particles, but the collection motion modes of them which we call as composite shells.This fact is pointed out and emphasised in our series of works [12][13][14][15][16][17].
Equation (IV.30) also tells us that m min outmost has the same order as the hawking temperature of a mass M tot black hole.This means that this shell can be easily radiated away as hawking radiations.FIG.7 displays the wave function modular square of the fundamental microscopic state of a 3M pl black hole caused by a dust ball's gravitational collapse.From the figure we easily see the out most shell has very large probability be found outside the horizon.Of course this outmost shell's appearance outside the horizon itself is not hawking radiation.The true hawking radiation arises from the microscopic state's change of the system and can be described as the spontaneous radiation of usual atoms with explicitly hermitian hamiltonian, which is named gravitation induced spontaneous radiation in [17], see references [13,14,16] for more earlier and detailed investigations.
The physical picture and calculations up to here provide a quantum mechanics based general relativity interpretation for the Bekenstein-Hawking entropy and area law formula.This interpretation contrasts with the common belief that Bekenstein-Hawking entropy has pure geometric origin, i.e. the fundamental degrees of freedom of black holes is carried by their horizon, each planck unit of area element carries one bit of information.However, our calculation shows that such fundamental degrees of freedom may completely be carried by the bulk motion modes of the matter contents.As long as we find the right carrier the area law formula flows logically and naturally.Since in a collapsing star consists of dust only, no other forces except gravitation is involved, the matter's motion and fluctuation is an exact and faithful reflection of the geometric feature of the spacetime underground.As results, the quantisation of these matters' motion is also the quantisation of spacetime itself.In a sense what we provide in this work is a quantum gravitation theory indeed.Since it focuses on the origin of Bekenstein Hawking entropy exclusively, we consider it a mini-version of some underlying full quantum gravitation theories.

V. A WORKING EXAMPLE FOR THE QUANTISATION
Let us consider the inner structure's quantisation of an M tot = 3M pl black hole as a concrete example for the microscopic state definition ideas above.For this black hole, to precision ε = 1, its matter contents can be divided into 4 distinguishable concentric shells at most, under the condition that n i =2, γ i =0 for all i.The fundamental state of the whole system is defined by the following shell mass assigning and radial quantum level number's setting scheme, {m i } = {1.25118,0.718127, 0.559288, 0.471405}(V.1){M i } = {1.25118,1.96931, 2.5286, 3} (V.2) 2 } = {0.5, 0.5, 0.5, 0.5} (V.4) x.By the representation method of (IV.21), this ground state can be represented as The modular square of the ground state wave function of an M = 3M pl black hole.
FIG. 7 displays the modular square of the wave function of each sub-shell.In the figure, three things are worthy of noticing.The first is, the positions of the global maximal value of each sub-shell's wave function are ordered from inside to the outside.The second is, the global maximal value of all shell m i 's wave function happens outside r i h .This is nothing but the requirement that classically each shell lies outside the horizon r i h = 2G i i ′ =0 m i ′ determined by the matter mass they wrap.The third is, the probability of each shells' being found outside the global horizon defined by r h = 2GM tot is rather remarkable.Keeping the condition x ψi max < x ψi+1 max , ∀i and exciting all shells m i to their higher level quantum state uniformly by setting n i = 3 or more larger is possible, but that would make the probability of these shell's being found outside the global horizon becomes too larger than the ground state characterised by n i = 2, ∀i.Now if we let some but not all of the inner shells such as m 0 to be excited, we will find that even only let n 0 = 3, the global maximal value of ψ 0 (x) would happen outside that of ψ 1 (x).So under the requirement that all other shell's quantisation conditions are kept invariant, i.e.GMimi √ 2 = 2(i̸ =0&1), the only possible way of exciting m 0 is letting it be bounded with m 1 as a single shell.Similarly, if we want excite other one or more shells independently, the best way is also to let them be bounded with their closest outer neighbour into a single shell.For the M = 3M pl black hole whose fundamental state is defined by the four n i = 2(i=0, 1, 2, 3) shells, we have only 7 = 2 4−1 − 1 possible ways to recombine those shells.Let us consider them one by one.
Recombination A, x.This is the same as (V.6) and similarly in the other recombination schemes.
To satisfy the constraint 2GM ′ 0 =3.94⩽x 62, we find that n ′ 0 = 5 are allowed, other values are not.So this shell recombination scheme contributes only one microscopic state to the system.By the representation method of (IV.21), this state can be represented as To satisfy the condition 2GM ′ 2 =6⩽x ψ ′ 2 max , we find that n ′ 2 = 4 is the lowest requirement.Other greater than 4 values for n ′ 2 could also satisfy this constraint.But the corresponding wave function will have a very fat non-zero tail outside the global horizon thus contradict with definition of black hole in general relativity.We will exclude such state as the microscopic state of the black holes characterised by the classic horizon r h = 2GM tot .This means that the shell recombination scheme C also contributes only one microscopic state to the system.The diagram representation of this state can be plotted as follows To satisfy the condition 2GM ′ 0 =5.06⩽x 85, we find that n ′ 0 = 7, 8 are both allowed.Other values are not.
So this recombination contributes 2 microscopic states to the system.The diagram representation of these two states can be plotted as follows To satisfy the constraints 2GM ′ 0 =3.94⩽x , we find that n ′ 0 = 5, n ′ 1 = 4 are allowed.Other values for n ′ 0 and n ′ 1 which satisfy this three constraints are also possible, but all of them have the problem of giving the corresponding wave function a too fat nonzero tail outside the horizon thus contradict with the definition of black holes characterised by their classic horizon.So such highly excited state will not be considered as the microscopic state of the class black holes.The diagram representation of this state can be plotted as follows To satisfy the condition 2GM ′ 1 =6⩽x ψ ′ 1 max , we find that n ′ 1 = 6 is the lowest level required.Other greater than 6 levels for n ′ 1 which satisfy this condition are also possible, but has the problem of giving the corresponding wave function a too fat tail outside the classic horizon.So they will still not be taken as the microscopic state of class black holes.The diagram representation of this state can be plotted as follows To satisfy the condition 2GM ′ 0 =6⩽x ψ ′ 1 max , we find that n ′ 0 = 10 is the lowest level required.Other greater than 10 levels for n ′ 0 satisfying this condition are also possible, but has the problem of giving the corresponding wave function a too fat nonzero tail outside the classic horizon.So still can not be taken as microscopic states of the classic black holes.The diagram representation of this state can be plotted as follows From the calculations above, we can see that to precision ε = 1, an M = 3M pl black hole can be divided into 4 distinguishable concentric shells at most under the condition that n i =2 and γ i =0 for all i.Letting all this four shells lie on the quantum state of n i = 2 yields the fundamental state of the system.In this fundamental state, the global maximal value of each shells' wave function modular square happens outside the classic horizon determined by the masses the shell wraps To excite each shell independently and preserve the ordering of the excited shells is equivalent to recombine shells from the fundamental state and tuning the resultant shell's radial quantum number.Including the fundamental state itself, we have 2 4−1 = 8 ways to to such a recombination.Except the recombination scheme D, the number of ways to excite all shells in all the recombinations is unique.Generalising to the large mass case, this is our conclusion proved analytically in the previous section, the matter contents of a spherical symmetric black hole of mass M (mildly larger than M pl ) with asymptotic horizon r = 2GM can be divided into k = GM 2 2 √ 2 distinguishable concentric shells to precision ε = 1.These fundamental shells can be recombined and excited at 2 k−1 X(k) number of ways.X(k) is at most a polynomial function of k.

VI. BLACK HOLE COMPLEMENTARITY PRINCIPLE
By our proof and calculations above, one can easily understand that in the co-moving observer's description, the gravitational collapse passes through the horizon and Schwarzschild singularity fluently in finite and seeable proper time.The resulting black hole owns an quasiperiodically oscillatory matter core and multiply possible mass profiles which is ergodically accessed by the evolution of system.While in the outside fixed position observers' time definition domain, a collapsing star can infinitely approach but will never contract into the horizon defined by its total mass; the system forms only an asymptotically implementable horizon and asymptotically linear mass function, the microscopic state of the system is characterised by its mass function's deviation from the exact linear profile M [r] = r 2G .The gap between this two observers' description is filled up as long as we note that, the horizon and singularity encountered by the co-moving observer happen only in the Future of infinitely Far away Future (FiFF) by the outside fixed position observers' physical time definition.They are realities only in the parallel universe or ensemble spacetime of these outside fixed position observer or detectors.
This fact naturally brings us to the idea of black hole complementarity principle [100][101][102][103][104], which conjectures that both the co-moving observer and the outside fixed position observer's description of a collapsing star's evolution are equally true.What we want to complement here is that, when the ergodicity feature of the co-moving observer and its observables' interaction dynamics and the uncertainty principle implied parallel universe or ensemble spacetime interpretation for the outside fixed position observers' measurements are considered, the two types of observer's definition for the black hole microscopic state will not only be equally true, but also be equally complete.In another word, the two observers' description are equivalent.By the extended Penrose-Carter diagram, we display this equivalence in FIG. 8.
In the quantisation stratagem of references [13][14][15][16] we assumed γ 2 i < 0 for each single shell when do the microscopic state number counting for black holes caused by gravitational collapse.Technically, this stratagem is nothing but the requirement that classically all shells consisting of the collapsing star oscillate inside the horizon defined by r < r i h ≡ 2GM i so that for these shells both ṫ and ẋ take imaginary values but the four velocity normalisation condition h ṫ2 − h −1 ṙ2 = 1 holds just as the shell is moving outside the horizon r i h .While in our proofs of this work, section IV, this stratatem is given up and the same area law is obtained successfully.Physically this means that we do not need any shell m i to move inside its horizon r i h = 2GM i , we only need to consider their motion of falling towards there.The fact that the two proofs lead to the same area law entropy formula implies that, there is an equivalence between the two pictures of (i) the motion of matters consisting of a black hole is their oscillation over-crossing the central point and (ii) the motion of matters consisting of a black hole is its composite shells' falling towards their own's horizon.In the latter picture, the quantum fluctuation implied by the uncertainty principle has the power of causing various shells to escape away from their asymptotical horizon thus making their inside horizon space looks like something just be dug out from the whole spacetime and antipodal point identified [162][163][164][165][166][167][168][169][170][171][172].
Although our proof or illustration of the equivalence relation between the two observers' definition of the black hole microscopic states is based on the physical black holes with non-singular inner mass distribution.We emphasise here that, even for the mathematical black holes defined as the exact solutions to the vacuum Einstein equation thus characterised by the apriori horizon and singularity.Similar complementarity principle must also be possible be understood from the two types observers' microscopic state definition.Because as long as such microscopic state are physical, they must be observable in principle.For such black holes, the idea of black hole complementarity principle is challenged by the firewall paradoxes in recent years [105,106] By definition, the black hole in the fire wall paradox is nothing but a horizon wrapped singularity and between the horizon and the singularity is a vacuum region with wrongly signatured space-time metric.By the fire wall paradox's reasoning, the hawking particle arises from the particle pair production and partial escaping from the vacuum fluctuation around the horizon.The escapers have nothing to do the microscopic state of the black hole.At late times any given escaper is entangled with both its partner falling into the horizon and the early escaping successful colleagues as a whole.This two side entanglement would break the monodromy of the quan-tum entanglement thus calls for protections from mechanisms such as firewalls.Obviously, for black holes caused by gravitational collapse, exact horizon and singularity never form to the outside observer or detectors.To them the black holes are inner-structured and characterised by special radial mass profiles.Referring FIG. 9, late time hawking particles can be recorded only when the in-going member of the early produced pairs get annihilated with matter particles consisting of the collapsing star because otherwise we will have vacuum decay inevitably.Effectively, this makes each hawking particle looks like arising from the microscopic state change of the black hole instead of completely random fluctuation of the vacuums around it.References [13,14,16] provide an explicitly hermitian hamiltonian description for this radiation mechanism which is named gravity induced spontaneous radiation in [17].
FIG. 9: The vacuum fluctuation and partial escaping mechanism for hawking radiation is a policy of deficit financing.The cost of the outside Gieger counter's recording hawking particles is the vacuum definition standard's lowering, an analog of the currency releasing standard's lowering in conventional financial systems.An event of particle pair's annihilation would uplift this vacuum standard, analog of the currency releasing standard's up-tuning.The potential curve V [ϕ]'s dash to solid style change denotes the vacuum definition standard's change.For a black hole caused by gravitational collapse, the horizon is only an asymptotically realisable but not a truly implemented physic surface so matter particles of the system can be found outside of it.The time separation ∆t between the Gieger counter's recording particle and the ingoing particle's annihilation with matter particles consisting of the collapsing star is very short so that ∆t•∆E ≈ ℏ 2 is holden very well.
For the black holes defined as exact solutions to the vacuum Einstein equation so carry absolute horizon and singularity apriori, hawking particles cannot be produced continuously as the results of particle pair production and partial escaping from the vacuum fluctuation around the horizon.This is because the hawking radiation recorded FIG.10: For black holes defined as exact solutions to the vacuum Einstein equation, the horizon is absolute and apriori.The in-going particles arising from vacuum pair production needs infinite t-time to go across the horizon and annihilate with matters inside to uplift the standard of vacuum definition.During this period, if the outside fixed position Gieger counter records hawking particles continuously, the vacuum definition standard would be observed to be lowered remarkably.This would cause vacuum decay.Particles would be observable everywhere, just like currency inflations in the conventional financial system.So the deficit financing policy interpretion for hawking radiation's production in these black holes does not work at all.That is, the outside fixed position observer could not see the black hole's evaporation as the result of vacuum fluctuation caused particle production and partial escaping.
by the outside fixed position observers comes from the vacuum definition standard's lowering.The particle pairs annihilation take the responsibility of uplifting this back thus retrieving such standard.To those observers, each member of the particle pairs falling towards the horizon needs infinite time to get across the horizon and reach the central point which accumulates all matter contents of the black hole to annihilate with particles there so that the vacuum definition standard can be uplifted back, see FIG. 10 for intuitions.If an outside fixed position observer records hawking particles continuously in his/her detector, then his/her vacuum definition standard would be lowered continuously without compensation.As results, he/she would see infinite number of particles everywhere as the results of vacuum decaying.So, for black holes defined as exact solutions to the vacuum Einstein equation, the situation imagined in the firewall paradoxes where the black hole radiated half of its mass will not appear at all.In fact the fraction of this black hole's mass change due to hawking radiation is strictly constrained by the uncertainty principle ∆M •∆t ≈ ℏ 2 .So the complementarity principle is right enough to be taken as a basis for understanding the black hole inner structures and their microscopic state definition.What the firewall paradox challenges is not the complementarity principle, but the underlying mechanism of hawking particle's generation mechanism.The right mechanisms such as gravity induced spontaneous radiation must contain some ingredients which couples the hawking particle's observational feature such as interaction amplitude with the microscopic state change of the radiating black holes, so that the information encoded in the initial state of the black hole can be transferred into observational features of the hawking particles.See references [13,14,16,17] for more concrete and detailed implementation of this idea.

VII. THE PHYSICAL REALITY'S OBSERVER DEPENDENCE
Why is the singularity theorem chosen by so many people as the basis for understanding the inner structure and microscopic state of black holes in general relativity?Why is it not the complementarity principle interpretation?Pushing aside the establishing time reason, the other one important reason for this status may be people's mistaking of the global viewpoint as local viewpoints and negligence of the physical reality's observer dependence, under the name of the physic law's general coordinate invariance.
It is well known that in a black hole formed through gravitational collapse, the outside fixed position observers will not see the formation of horizon and singularity in their any finite future.These observers can be any outside test particle or probes whose motion is affected by the gravitational field of the to be black hole collapsing star, such as the other partner in a binary merger event and the light rays skimming over the surface of a single black hole to be imaged.During the whole lifetime of these test particle's existence, their motion is affected by the gravitational field provided by the matter contents of the collapsing star being measured.So they see no black holes defined by the horizon and singularity at all, they see just a to be black hole collapsing star.
Due to the co-moving observers existence and the singularity theorem's being proven, many people believe that the fate of all heavy enough collapsing stars in the nature is the horizon wrapped singularity definitely.The logic behind this belief is, the physics laws are invariant under the general coordinate transformation.In that the singularity and horizon form in the finite future of the collapsing matter themselves' time definition, their effects on the motion of the outside fixed position observers must also happen in the form of horizon and singularity naturally.The outside fixed position observers cannot see such objects appearance only because they adopt a bad coordinate system or time definition.
However, this popular belief makes two mistakes rather hiddely, (i) taking the global viewpoint of the horizon's existence as the viewpoint of local observers and (ii) tak- TABLE I: Two types of observer's measurement of a light signal emitted from the surface of a freely collapsing star whose outside geometry is exactly Schwarzschild type.The upper is in the standard Boyer-Linderquist coordinate, the lower is in the Eddington-Finckelstein coordinate.As a physical reality, the light signal's frequency is invariant under general coordinate transformation but observer dependent.h0 is the value of h at the light signal's starting point.h fx is the value of h at the point the fixed position observer sits on.k µ is the wave number of the light signal, whose components are determined by k0 + Γ 0 µν k µ k ν = 0 and k•k = 0. u µ fx is the four velocity of the fixed position observer far away, u µ fr is the four velocity of the freely falling observer outside the star.The components of u µ are determined by u0 + Γ 0 µν u µ u ν = 0 and u•u = −1 ing the physics law's general coordinate invariance as the physical reality's observer independence.The physics law here refers to relations like G µν = 8πG N T µν which bridge observables following from at least two types of physical measurement.The physical reality refers to any physic quantity or physical features which is measured through only one type of physical measurement.These are two totally different concepts, the invariance of the physics laws under the general coordinate transformation does not imply the physical reality's observer independence in any sense.In fact, the physical reality's observer independence is not a fact at all.Let us elaborate this point through three examples.
The first is a concrete measurement example.Consider two observer's measuring of a same light signal's frequency emitted from an atom fixed on a freely collapsing star whose outside geometry is written in two different coordinate system.The first observer is a fixed position one, the second is a freely falling one.TABLE I lists these two observer's definition and measurement mathematics.From the Table, we easily see that, the frequencies measured by the two observers ω fx and ω fr are obviously unequal, although they are invariant before and after the coordinate is transformed from the Boyer-Lindrquist to the Edington-Finkelstein type, ω fx = ωfx , ω fr = ωfr .In this example, the frequency of the light signal, and even the redshift's dependence on the radial coordinate r is what we called physical realities.These realities are obviously observer dependent, because their value follows from the inner product of the observer's four velocity and the tensor representing the observable

itself.
The second is still a concrete measurement example.But this time let us consider two observers measurement of a same test shell's pressure when it is freely falling towards the horizon of a pre-existing Schwarzschild black hole characterised by the exact horizon and central singularity.The first observer is a fixed position observer which is sitting far away from the horizon of the black hole.The second is a freely falling observer co-moving with the test shell.We display these two observers' four velocity and measurements of the test mass shell's pressure in TABLE II.In this example, the pressures of the test mass shell is our physical reality.From the table, we can easily see that this physical reality is also observer dependent p fx ̸ = p co , although coordinate invariant, p fx = pfx and p co = pco .p fx has the property of approaching to infinite as the shell falls infinitely close to the horizon.So to the far away fixed position observers, the asymptotic mathematical horizon's existence makes the black hole look like an incompressible but finitely sized objects.p co has the property of being zero constantly, which means that to the co-moving observer, the horizon exhibit not any observational effect at al.
The third is a conceptual example.Let us consider the thought experiment cartooned in FIG.11   8πG N , H0 denotes just the Hubble constant of the current universe.Although divorced for almost 1.37 × 10 10 years already, the gravitation between them and their own ex-friend still makes Alice and Bob recede from each other everyday.From Alice and Bob, this is nothing but the usual freely falling.Charle, their common friend or a global observer, also thinks this is nothing but the conventional big bang.The question is, why is it possible Alice and Bob get known to each other and begin to date?Because each of them lives inside the horizon caused by her or his ex-friend but outside that of the other's.
a paradox we called F2FCD, see the figure and captions carefully.This paradox magnifies the error of the common belief that the horizon and singularity are observer independent physical realities to all black holes.Charle's confusion arises from two factors.(I) believing that in Alice's freely falling reference frame, Robert has a horizon and similarly in Bob and Walker's case.This is an error of taking the global viewpoint of Charle himself as the viewpoint of local observers such as Alice and Bob.(II) believing that the horizon-carrying feature of Robert to Alice is also the case to Bob and vice versa.This is an error of taking the physic laws' general coordinate invariance as the physical reality's observer independence.In fact, as an outside observer, in any finite time of his own, Bob can not see or detect Robert's horizon at all.So he can get known with Alice without any difficulty.This F2FCD paradox is firstly introduced in reference [14] to necessitate the horizon's fuzziness to resolve the information missing puzzle.But its value may be more right here to help understanding the physical reality's observer dependence.
As long as we know that the horizon-carrying is not a physical reality of black holes independent of observers, we can easily understand that the microscopic state picture of black holes is observer dependent.By our exposition in section II,III and IV, to the outside fixed position observer, the matter contents of physical black holes live inside a spatial region which is infinitely close to that wrapped by the mathematic surface r h = 2GM tot .But they live there just like some in-compressible and static fluid.Because inside that region, matters inside each sphere of smaller radius has the same feature.This is nothing but the concept of the old frozen star [60][61][62][63][64], see references [65][66][67] for recent investigations.Except that the in-compressibility now is understood as a result of the infinite gravitational time delaying instead of any physical exclusive forces.To the co-moving observer, these matters just freely falling towards the central point, over cross each other there and then oscillate periodically.But these inside horizon motion and oscillation happen in the future of infinitely far away future by the outside fixed position observers' physic time definition.They are physical realities only in the sense of parallel universe or statistic ensembles.Just as we showed above and in our earlier works [12][13][14][15][16][17], quantisation of these matters' fluctuation on the frozen star background or oscillations across the central point both lead to the right area law featured Bekenstein Hawking entropy, thus form equivalent and complementarity definitions for the microscopic state of black holes.

VIII. BANANA SHAPE DEFORMATION
It is very important to emphasise that the observers referred to in this work are not limited to human beings but are proxies of any probe or detectors, especially those living outside the matter occupation region of the being probed objects.They can be the imaging photons skimming over a single black hole or the other partner involved in a BHB merger event.These test objects' motion is determined by the gravitational field produced by the matter contents of the being tested black holes which never fully collapse and cause singularity in any finite future by the time definition of the outside fixed position observers.These being probed objects carry approximate horizon which is formed so long that during the short duration of the imaging process or binary merger process the imaging photon's perturbation or the inspiral partner's tidal force has no chances to cause disintegration.Because longer formation time means smaller matter occupation region radius thus more remarkable time delaying from the matter occupation region's boundary to the test particle's position, so that more time consuming for their matter contents to be affected by the imaging photon or merger partner's gravitation, and vice versa.In the BHB meger events, this provides us chances to treat dynamics of the binary system with the adiabatic approximation.By this approximation, we look the motion of the binary system as a combination of each participant's mass central's motion and their different parts' relative motion.The latter arises from the inhomogeneous radiation back reaction forces and will cause banana shape deformation of FIG.12.
The key reason for the banana shape deformation is the inhomogeneity of the gravitational wave radiation's back reaction.Referring to FIG. 12, let us consider the newtonian motion of a representative point such as ⃗ x i on the participant B and its counter partner on participant A as an illustration for this phenomena.By the standard quadrupole formula for the back reaction force arising from the system's inspiral motion and gravitational wave radiation, we can write down the dynamic equation controlling the inspiral orbit's decay of these two representative point and get solutions 3) The synchronicity of different composite points' inspiral motion an approximate constant for different point i.So as the inspiral motion progresses on those points with larger initial x i will reduce their orbit radius more remarkably ∆x In the case of all points inspiral at the same angular frequency, this will cause the compression of the system along the ⃗ x i direction, the net effect is a banana shape deformation.The true case in general relativity would be much more complicated than the newtonian illustrations here and can only be simulated numerically.But the qualitative physical picture will be the same.Obviously, when the two participating black holes experience banana shape deformation as in FIG. 12, the radiation activity of the binary system will decrease inevitably, because in this case the rotational symmetry of the system enhances and the radiation relevant quadrupole of the system decreases correspondingly.The extremal case is when the two ends of the banana shape deformed black holes contact each other.In such cases the gravitational wave radiation will stop completely as the result of rotating symmetry's retrieving.For physic black holes originating from gravitational collapse, nonsingular inner mass profile will make this banana shape deformation inevitable and will become very important as the two participants inspiral to the merger and ring down stage.From equation (VIII.5),we can see that this effect is of O[G 0 ], so it is non-perturbative and not accounted for in the higher multipole moments appearing as the response of the participating black holes to the tidal force caused by their partners and written into the tidal love numbers perturbatively.The tidal love numbers are well known to be zero for non-spinning black holes defined as exact solutions to the vacuum Einstein equation [143][144][145][146]. See FIG. 13 for intuitive pictures.
Since both of the merger participants are outside observers of their inspiral partner, neither of them would see the other's horizon.Both of them would see the other as some frozen star occupying invariant spatial volume.So their compression along ⃗ x i direction implies stretching along the ⃗ v insp direction inevitably.In a sense, the spatial volume contained in a spherical symmetric Schwarzshild horizon is the minimal volume we can compress matter contents into.This can be seen from the following fact that a spherically symmetrical collapsings star with linear radial mass profile M (r) ≈ r 2G implement the lowest possible self gravitating potential energy without causing the physical horizon's formation Note that M (r) ≈ r 2G implies ρ sph (r) ≈ 1 8πGr 2 necessarily.Compressing this star and tuning its mass value so that the surface of the resulting ellipsoid has equal newtonian potential [147] as the spherical symmetric star − GM a = − 1 2 , the self gravitation energy of the resultant object will become where F (k, φ) is the elliptic integral of the first kind and dv = abcλ 2 dλ sin θdθdϕ.In this calculation, the star is considered as many ellipsoid shape shells of equal axis length ratio {aλ, bλ, cλ} and radial mass function By our convention, the star is compressed along y and z direction but kept invariant along the x-direction.Letting the surface potential of the resulting ellipsoid take the same value as that of the Schwarzschild black hole GM F (k,φ) with the same x-directional size, we will find that M ell < M sph and E ell /E sph < 1. See FIG. 14 for concrete numerics.Conversely, it is very natural to understand that if we keep the total mass invariant when compressing the star along one direction, then the star must be stretched along the other so that the horizon formation condition will not be realised as any outside test body cannot detect it in the finite future by their physical time definition.
Considering the banana shape deformation, we will get similar conclusion.In the case the ellipsoid star are parameterised as, x=aλ sin θ cos ϕ, y=bλ sin θ sin ϕ, z=cλ cos θ(VIII.10)λ ∈ (0, 1), θ ∈ (0, π), ϕ ∈ (0, 2π) (VIII.11) the banana shape deformation can be written as whose own volume has the form 4πa bnn b bnn c bnn λ 2 dλ and wrap mass M • λ inside, then the self-gravitational potential energy of the system can be calculated approximately as ) where ψ is the bending angle of the banana, referring to FIG 14 for physical pictures; subscripts bnn means banana shape deformation.In the case of ψ ≪ π, equation (VIII.16)provides rather good approximation to the underlying exact result.For large ψ, equation (VIII.16) is not a good approximate, but still catches the qualitatively feature of banana shape deformation costing extra energy relative to the ellipsoid deformation.
Referring to FIG. 15, when we consider a banana shape deformed BHB system's rotation and gravitational wave radiation, the time dependent quadrupole will become FIG. 15: The radiation active quadrupole of a banana shape deformed BHB system is suppressed by a factor of ζ = sin α α .4GMi and 4GMizi are the deformation before and after diameter/length of the black holes respectively.≡ α 2 = α with respect to their inspiral central, the elongation factor z i is required to satisfy z1M1 M2 = z2M2 M1 ≡ z.As results, the radiation active quadrupole of the system can be written as Just as we pointed out previously, during the inspiral and merger process the banana shape deformed black hole would be elongated under the inhomogeneous radiation back reaction force, so z is a growing function of the time coordinate t.In this work, we will not try to determine the function form of z(t) dynamically but choose to parameterise it as various power functions of a(t) As time passes by, a(t) would decrease due to the relative motion orbit decay caused by gravitational wave radiation.So a(t) −n (0<n)s are increasing functions of time.z = 0 corresponds to the case the black hole has point like singular structure.z = 1 corresponds to the case the black hole has extended inner structure but such structures experience only banana shape bending but no elongation as the system inspirals and emits gravitational wave radiation.Black holes with z∝a(t) −1 experience more banana shape elongation during the early inspiral stage than those with z∝a(t) −2 so are more soft than the latter.It is necessary to note that although the black holes with non-singular inner-structure experience banana shape deformation, the protection provided by their asymptotic horizon would prevent them from tidal disintegration during the most duration of their binary merger process, from the early inspiral to the very late ring-down stage.However, to show this inner structure truly exists and happens there, we need to develop new theoretical tools to calculate their observational signal analytically.

IX. AN EXACT ONE BODY METHOD
The full process of a BHB system's evolution can be divided into three stages, i.e. inspiral, merger and ringdown.The post newtonian approximation and effective one body (EOB) method can be used to describe evolutions of the system during the first and early part of the second stage.Numeric relativity can be used to model the system's evolution during the whole process but is usually applied only for the second stage due to computational resource contstraint.The practical signal wave form used in observation are compositions of the three method's output, with the joining condition tuned carefully by hand.
Even one day, our computational resource gets rich enough to calculate the wave form of the merger process with numeric method exclusively, the popular belief that Schwarzschild black hole has singular inner-mass distribution is still not a conclusion extractable from such calculations.The reason is, Kerr black holes have extended inner structure embodied in their singular ring.If such a structure follows from the merger events of two Schwarzschild black holes, then the inner mass distribution of the parent Schwarzschild black hole must not be point-like singular.Instead it must be extended and will experience stretching and bending during the binary merger process.Signals of such structures would be hidden inside the boundary data [152][153][154][155] of numeric relativity, what we see in the inner-motion irrelevance is only an illusion.To understand the existence of such an inner-structure picture observationally, what we should do is not to escalate the precision of the existing numeric relativity's simulation by finding more precise initial data basing on higher order post newtonian expansions or something else.What we should do is to develop some more self consistent and full process applicable analytical description or analytical wave form generating techniques which directly translate the inner-structure of black holes to the gravitational wave forms measurable from the real BHBs' merger events.
The basic idea of conventional EOB method [119][120][121] is to map the general relativity two bodies' motion onto those of a total mass particle's shifting and a reduced mass particle's inspiral.Due to the non-linearity of general relativity, the gravitational field controlling the reduced mass particle's inspiral is not a Schwarzschild metric sourced by the total mass particle, but an effective and corrected but still static and Schwarzschild-like one.This effective metric must be constructed by requiring that the reduced mass particle's motion in it given by the world line action coincide with the two particle's relative motion in the original general relativity two-body problem under the post newtonian approximation.For the inspiral and early merger stage, this method can yield gravitational wave forms coincide with those following from numeric relativities very precisely.However, as the late merger and ring-down stage is coming, this method's inconsistence manifests.Most manifestly, the relative inspiral motion's frequency following from ω≡ φ and ω∼p ϕ /µa 2 exhibits opposite trends of variation as the time passes by, see FIG. 16 for illustrations.This means that this method cannot be used to tell us wheter the two Schwarzschild black holes involved in the binary merger process own inner-structure and experience shape deformation or not.
FIG. 16: The relative motion orbit of general relativity two particle system and the corresponding gravitational wave form following from the conventional EOB method.This method yields manifestly opposite variation trends for the system's inspiral frequency defined through ω≡ φ and ω∼p ϕ /µa 2 as the system evolves to the very late stage.FIG.17: The spacetime is dynamically partitioned into three regions, red, green and their inter-connection region.Each region has its own static geometry whose application region rotates synchronously with the merger participants.For our purposes, we only need to know the metric functions' form in the red and green region and in the equatorial plane.
The conventional EOB method's inconsistency arises from its decomposition of the essentially time dependent radiation geometry into "a static background gravitation + dissipative test particle".This doing is not self-consistent because a test particle's motion in a fixed background spacetime is well known to be conservative so no radiation caused dissipation is allowed at all.To avoid this inconsistency, we come back to the original general relativity two-body action and write the time dependent radiation geometry as the combination of three independent static patches 1 with time-dependent application region, see FIG. 17 for intuitive pictures.On the motion plane, the two patches' accommodating the two merger participants and rotating synchronously with them can be written as, , ϕ 2 (τ )} (IX.4)By this application-region wise metric, the action of the general relativistic two-body system becomes Since the orbit decay caused by gravitational wave radiation is always slow relative to the rotation of the system, we can always neglect the variation of the two participants' radial coordinate in a single circular orbit motion period thus setting ȧ1 = ȧ2 = 0.So for the circular orbit motion, the variation principle for a 1 and a 2 will yield that For the non-circular orbit motion, as long as we understand a 1 and a 2 as the length of the semi-major axis of the two participants' approximately periodic elliptic orbit, equation (IX.6) will also be valid.At the same time, to assure the unshifting of the two merger participants' central of mass position, we impose As results, the synchronicity of the two participants' inspiral motion implies that 1 To determine the binary system's inspiral motion, what we need is only the space-time metric on the motion plane instead of the whole space.So exactly speaking, the metric expressions (IX.2)-(IX.3)are valid only on the motion plane.
It can be easily verified that, by this choice of E 1 and E 2 , equations (IX.6) will lead to a general relativistic version of Keppler's third law for the reduced mass particle's motion in the total mass particle's gravitational field.Equations (IX.7) and (IX.8) transform the two-body problem (IX.5) into a single body one M 2 is always less than 1.The logic of the one bodylisation method (IX.9) is to capture the non-linearity of general relativity gravitation non-perturbatively through the non-perturbative central fixing condition (IX.7) and the relative motions' synchronicity condition (IX.8).In contrasts, the logic of conventional EOB method of [119][120][121] is to capture the non-linearity of general relativity gravitation through the perturbative match of the relative motion hamiltonian with that of a test body's motion in the background of a static and corrected Schwarzschild-like field.Since our one bodylisation conditions (IX.7)-(IX.8)are exact and call no post-newtonian approximation as input, we name it as an eXact One Body (XOB) method.
For those who can't help from evaluating how far is our one body-lisation idea to the conventional EOB method, the post newton expansion of the hamiltonian following from the conservative part of (IX.9) is as follows (IX.10) Except H 0 , all the higher order terms do not coincide with the input of the conventional EOB method [119][120][121].However, just as we pointed out above and illustrated in FIG.16, including the dissipation term in the hamiltonian of the conventional EOB method is not selfconsistent.The inspiral frequency ω ≡ φ and ω∼p ϕ /µa 2 following from such a method exhibit manifestly opposite variation trends at the late time of the merger process.In contrast, adding dissipation term in our hamiltonian (IX.10) is self consistent, because the frequency ω it involves plays the role of the merger participants' inspiral speed and the resultant spacetime's rotational speed simultaneously.So the whole system allows radiation caused dissipation consistently.
For the above reason, we will not take the higher order post newton expansion of (IX.12)-(IX.14)differs from the input of reference [119][120][121] as a negative diagnosis for our idea's rationality.Instead we will take the direct gravitational wave form output of our method as the justification for its validity.Taking the most simple quadrupole radiation power as source of dissipation, we will get evolution equations for a(t) From this equation we can solve the function a(t) and calculate the gravitational wave through We compare in FIG.18 the relative motion orbit and gravitational wave forms following from the EOB method and our XOB method.From the figure, we can easily see that the two orbits and wave forms highly agree with each other till the inner most circular orbit (ISCO) epoch of EOB.However, XOB exhibits much rational or selfconsistent behaviour after that epoch.For example by XOB, the late time inspiral speed defined through ω = φ and ω ∼ p ϕ /µa 2 (p ϕ ≡ δL δω ) are both growing functions of time.While the variation trends of those predicted by EOB are contrary with each other.
FIG. 18: The relative motion orbit and gravitational wave form of an equal mass BHB system following from the standard EOB method and our XOB method.Both members of the binary system are not inner-structured and experience not any shape deformation during the merger process.The dotted magenta line marks the ISCO radius defined by the EOB method.The dashed red line is the horizon radius of the total mass black hole r h = 2GMtot.
Two very impressive features of the gravitational wave forms following from XOB displayed in FIG.18 are (i) the decay of the gravitational wave forms is invisible from the time interval displayed and (ii) the asymptotic frequency ω∼p ϕ /µa 2 of the inspiral motion is logarithmically divergent.Both features imply that the binary system has infinite energy to be carried away by the gravitational wave radiation.This is provided by the infinite decreasable space of the gravitational potential energy between the two point particle type participating black holes.In the language of newtonian mechanics, the is just the reflection of ε= − GM1M2 a a→0 − −− → −∞.While in the language of general relativity, ε∝ M } is the asymptotical separation the two black holes can get close to, see equation (X.3) of next section for references.In the case of M 1 = M 2 , this asymptotic value a f = 3GM 4 is less than the mathematic horizon radius r h = 2GM .
From FIG. 18, careful readers may see that XOB method has the power of traceing relative motions of the binary system and apply the quadrupole formula (IX.16) for the gravitational wave forms even when a f <2GM tot , i.e. when the two participating black holes enter into the horizon determined by their total mass.This is rational because to each of them, their partner's horizon is only 2GM i which is less than a f .So all of them are outside observers of the other during the whole inspiral, merger and ring-down process.The horizon defined by r h = 2GM tot is never implemented at all.So, tracing relative motions of the binary system inside its mathematical horizon does not violate any physical law such as causality at all.In a sense, this is a conceptual revolution.If one cannot understand this properly, he/she will have difficulties to understand why the outside fixed position observers can see the gravitational wave forms reflecting the relative motion of the two black holes inside the horizon defined by their total mass r h ≡ 2GM tot .This will be even more clear when we consider the nonsingular inner structure and banana shape deformation of the participating black holes.In such cases, XOB's ability to trace the system's evolution inside the mathematic horizon will bring us the full gravitational wave forms with well behaved late time quasi-normal mode feature and decodable information about the inner structure of the participating black holes.
By our XOB method, during the BHB merger process the event of exiting from innermost stable circular orbit (ISCO) does not happen on the epoch defined by the EOB method a = 5.718GM tot .Because, by XOB method even after that epoch, the relative motion orbit of the system is stable and the system's evolution after the horizon entrance r < 2GM tot is still traceable analytically until the system enters the ring down phase.However, for comparing and referring convenience, we will preserve this nomenclature and take 5.718GM tot as a ISCO , and mark it out in most of the relative motion orbit and gravitational wave form picture below.

X. THE FULL GRAVITATIONAL WAVE FORMS
Consider the non-singular inner-structure of black holes and their radiation activity suppression caused by banana shape deformation, we have to write the relative motion orbit decay equation as where ζ is the radiation activity factor of the binary system and z grows with time slowly, see equation (VIII.20)and comments there for their definition and concrete expressions.Focusing on the circular orbit or understanding a only as the semi-major axis length of the quasiperiodic elliptic orbit, In the special case of In the general case of ν⩽1/4, the function form of dH da will be more complicated.But its basic feature of containing a zero point somewhere inside the region 32 ⩽ GM a ⩽ 6 is very robust, referring to FIG. 19 for illustration.Denoting this value of a as By equation (X.1), the zero value of dH/da will manifest as a divergent point on da dt .However such a divergence will not bring us disconnection on the variation of a(t).So it can be cured by hand in the numeric integration of eq(X.1).For latter uses, we define another two extra special values of a through where a ∞ H is the point on which H itself diverge and a 0 ζ is the maximal zero point of ζ.In the special case of ν = 1 4 , according to equation (X.3), we can easily see that a ∞ H = 3GM 4 .On the evolution line of a(t), whether the divergence implied by dH da = 0 will be implemented depends on the relative relation between a 0 H ′ and a 0 ζ , see FIG. 19 for references.
FIG. 20: The relative motion orbit and gravitational wave form of binary systems consisting of equal mass members following from XOB method.The participating black holes in the upper row have point like singular inner-structure thus are described by the standard Schwarzschild metric exactly.Their shape deformation parameter z≡0.The participants in the lower row have extended inner-structure but fixed length of 4GM .So the shape deformation parameter z≡1.This inner-structure will experience bending but no elongation as the merger progresses.The dotted magenta line marks the ISCO radius defined by the EOB method.The dashed purple line is the mathematic horizon radius r h = 2GMtot.
Since ζ= sin 4GM z/a 4GM z/a , the concrete value of a 0 ζ depends on the form of z(t) inevitably.The standard Schwarzschild black hole is defined by the exact horizon and singularity.By the language of newonian mechanics, its matter distribution is singular and experiences no stretching or bending.By our definitions of banana shape deformation, such black holes have z ≡ 0 and the corresponding radiation activity factor ζ ≡ 1.In the case the two members have equal mass, equation (X.1) allows us to trace their merger evolution down to a(t) The former equals to half of the real part of the late time quasi-mormal mode's frequency ω re 022 .The latter has the explanation of the final black hole's total mass or energy.Its divergence and negativity is a reflection of the irrationality of the point particle type inner-structure picture of black holes.The upper part of FIG.20 displays the full three stage gravitational wave forms correspondingly.To make the wave forms readable, we reduced the late time frequency and amplitude of the gravitational wave by a factor of 0.3.The decay feature of late time gravitational wave forms is non-visible in the time interval displayed.The logarithmic divergence feature of ω ′ ≡ p ϕ /µa 3 is due to the point particle type inner structure of the black holes involved.So this figure tells us that, if black holes have point like singular structure, then we cannot get quasi-normal mode type late time feature on the gravitational wave form of their binary merger event.Now let us consider the extended inner structure and non-singular matter distribution of black holes.As first step, we let the black hole have fixed length which equals to the diameter of the corresponding Schwarzschild horizon so that 4GM z=4GM , or z≡1 constantly.Although the length is fixed, we let the extended black hole be bendable as the inspiral and merger progresses, so that the radiation activity factor varies with time in the most simple but still non-trivial way ζ = sin 4GM/a 4GM/a , see equation (VIII.20)for the derivation this expression.When this radiation activity factor comes in, the asymptotic value of a(t) is no longer the singular point a ∞ H of the hamiltonian, but a 0 ζ .In the case the two participating black holes have equal mass, a The former still corresponds to half of the real part of the late time quasi-normal mode's frequency ω re 022 .The latter equals to the final black hole's total mass or energy.Combining with the full process gravitational wave form displayed in the lower part of FIG. 20, we easily see that, relative to those following from the exact horizon wrapped point particle type inner-structure picture, the gravitational wave form following from inner-structure picture characterised by extended inner-mass distribution with fixed length but banana shape bendable and wrapped only by approximate or asymptotic horizons are much more reasonable, due to its late time quasi-normal mode type oscillation.
In the practical BHB merger process, the degree of the participating black holes' banana shape deformation varies with time.Just as we pointed out in section VIII, when such banana shape deformation happens, the participating black holes will be elongated along the inspiral direction.So during the whole merger process, we expect their length varies from zero to some final length value 2πa f .Substituting the value of this final length and the corresponding ω f = GM a 3 f 1/2 into the hamiltonian (IX.10), we should get an energy less than the total mass of the initial black holes This is because the gravitational wave radiation dissipates energy of the system inevitably.From this inequality, we can solve out the lower limit of a low f,lim and a corresponding frequency ω f , which equals to half of the real part of the merged black hole's quasi-normal mode's frequency f,lim On the other hand, from the hamiltonian expression (X.3) we note that given ν and M , H is a lower bounded function of a, i.e.H ν,M min <H ν,M (a).So we cannot set the mass of the final merged black hole to arbitrarily small value.Setting H[a, ν, ( GM a 3 ) 1 2 ] to this minimal value H ν,M min , we will get an upper limit for the value value a upp f,lim and a corresponding lower limit on the value of ω re However, since we have not any physical principle to prohibit a black hole be banana shape deformed and elongated longer than 2πa upp f,lim , so that the corresponding lower limit ω re022 lolim can only be considered a qualitatively estimation for the value of ω re 022 .The upper part of FIG.21 displays H[a, ν] and ω re022 uplim 's dependence on a.The lower part of the figure displays ω re022 lolim 's dependence on the symmetric mass ratio ν.Also displayed on the downer part of the figure is the currently available observational data of [25][26][27].From the figure, we can easily see that, our definite prediction for ω re022 uplim is perfectly satisfied by the observational data except one point.While ω re022 lolim gives rather precise estimation for the value of ω re 022 .By this estimation, when ν = 1 4 , the value of a upp f,lim = 3 2 GM .This means that the two participating black holes are elongated to πa upp f,lim = 3π 2 GM as the merger finishes.When ν = 0, the value of a upp f,lim = 6GM .That is, the final black hole in this case should be considered a rotating ring of perimeter 6πGM .In principle, by tracing up the matter distributions' evolution under the radiation back reaction force's driving, it is possible to determine the function form of z(t) theoretically.However, such a work is very difficult to accomplish in practice.So we will illustrate the observational feature of the black hole deformability's difference by comparing the gravitational wave forms following from two shape deformation modes of The corresponding radiation activity factors read The black holes characterised by z 1 (t) experience the banana shape deformation more intensively at the early stages of the merger process, while those characterised by z 2 (t) experience the shape deformation more intensively at the late stages of the merger process.
FIG. 22: The relative motion orbit and gravitational wave form of BHB systems consisting of equal mass members with inner structure and experience banana shape deformation.
The banana shape deformation parameter in the upper black holes varies with time at the way z∝a f a(t) −1 , the lower varies at the way z∝a 2 f a(t) −2 .The the gravitational wave amplitude in the lower figure not only grows more suddenly, but also decays more suddenly around the a ≈ 2GM region.The dotted magenta line marks the ISCO radius defined by the conventional EOB method.The dashed red line is the horizon radius of the total mass black hole r h = 2GMtot.FIG.22 displays the relative motion orbit and the corresponding gravitational wave forms of these two types of black holes.From the figure, we can see that the z∼a(t) −1 black holes experience banana shape deformation more mildly and more intensively during the early stages of the inspiral motion, so their relative motion orbit radius evolves to the asymptotic value a 0 f in longer time durations.Numerically this is t 1 isco = 222, t 1 h = 302 while t 2 isco = 220, t 2 h = 293.We will call this type of black holes softer than those characterised by z∼a(t) −2 .The most general form of the banana shape deformation function can be written as Except a f , all other parameters k i in this function are dimensionless.For these dimensionless parameters, we conjecture that they are independent of the black hole sizes.By carefully measuring of the late time gravitational wave form following from various BHB merger process, we have chances to determine these parameters observationally.So wether the black holes in the nature have inner structures and how hard such inner structures are are a question answerable experimentally.

XI. CONCLUSION
In one sentence, our conclusion can be stated that black holes in the nature have microscopic states understandable in standard general relativity and innerstructures testable from the gravitational wave forms produced in the BHB merger events, which is calculable analytically with our exact one body method.Expansions of this sentence are as follows.
We provided two exact time dependent solution families to the Einstein equation sourced by dust fluids and showed how such solution families form basis of black hole microscopic state both from the viewpoint of outside fixed position observers and from the viewpoint of inside co-moving observers.By the canonic quantisation scheme, we proved that the degeneration degree of such solution families has right exponential-area law feature as required by the Bekenstein-Hawking entropy formula except some polynomial corrections.Singularity theorem judges impossible to understand black hole microscopic state in general relativity only because, in applying such a theorem, the horizon and singularity's observer dependence is neglected under the name of physical reality's general coordinate invariant.
A time dependent collapsing star solution indeed contains a future singularity and an asymptotically realisable horizon.But such a singularity and horizon would not become realities in any finite future of the outside fixed position observer by their physic time definition.They will be realised in the finite future only to the freely falling co-moving observers which use the proper time of the collapsing matter themselves as time definition.In another word, the horizon and singularity are physical realities only to the super investigators who adopt the global viewpoint or have the ability to see the future of infinitely far away future of the outside fixed position observers.To physical probe or observers such as the participants of BHB merger events or freely falling observers co-moving with the matter contents of a collapsing star, the horizon and singularity of their being observed objects are either non-implementable in any finite future or although reachable in finite future but non-detectable experimentally.The belief that black holes in the real world carry exact horizon and singularity makes two mistakes of (A) taking the horizon and singularity of the global viewpoint holders as the horizon and singularity of the local viewpoint holders and (B) taking the physical law's general coordinate invariance as the physical reality's observer independence, which is not the case at all.
Since either participant involved in a BHB's merger event can only see their merger partner's regular inner structures characterised by the approximately linear mass function and an asymptotic horizon, both of them will experience banana shape deformation under the inhomogeneous back reaction arising from gravitational wave radiations.We introduced the concept of radiation activity ζ to quantify this deformation caused suppression of the system ability to emit gravitational waves and derived out that ζ = sin 4GM z/a 4GM z/a , where z is the elongation times of the banana shape deformed black holes relative to their shape deformation before diameter.ζ will be multiplied on the quadrupole of the binary system when the quadrupole formula is used to calculate the gravitational wave radiations of the system.
We then proposed an exact one body-lisation method to transform the general relativity two-body problem into a single body one.This method gives up the static spherical symmetric background of EOB and the conservative part of the system's action calls no post newtonian approximation as input.It puts the two merger participants in a three partitied spacetime which rotates synchronously with the merger participants.The metric of the spacetime region that accommodate the merger participants can be determined non-perturbatively by requiring that the two participant's inspiral motion is synchronous and central fixed.This method avoids the intrinsic inconsistence of the conventional EOB and provides descriptions for the relative motion of BHB systems throughout the full inspiral, merger and ring-down stages.Taking the quadrupole radiation as the source of dissipation, and neglecting the inner structure of black holes by looking them as point particles wrapped by exact horizons, we calculated the relative motion orbits and the gravitational wave forms following from the BHB merger process and get results which are highly consistent with those from conventional EOB but behave more rationally at the very late time stages.
When considering the extended inner structure of black holes, even the most simple fixed length banana shape deformation model, the gravitational wave forms following from our XOB method are almost the same as those following from the combination of EOB+NR+BHPT methods and measured in the real observations.We derived out an exact upper bound on the real part of the final black holes' quasi-normal modes ω re022 upbnd irrespective of the black hole inner structure details as a concrete prediction of XOB.We also derived out a not so strong lower bound ω re022 lobnd for the late time quasi normal modes which can be used as an estimation for the typical value of of ω re 022 .Putting the currently available data on the theoretical prediction figure we see very good agreement between the two.We compared gravitational wave forms following from BHBs whose members are characterised by different shape deformation susceptibilities and point out that the differences between them may probably be measurable in the current and future observations.So the question on what inner structures the physic black hole has is answerable observationally.
Two possible directions for future works would be, to find observational evidences for the black hole inner structures we proposed here through other channel such as black hole images and gravitational wave echos, and to revise and apply our XOB method to BHBs or neutron star binaries whose members are spinning and carry nonzero angular momentum.

I. Introduction 1 II. The outside fixed position observer's description 3 III. 6 IV. Quantisation 8 V
The inside co-moving observer's description .A working example for the quantisation 12 VI.Black hole complementarity principle 15 VII.The physical reality's observer dependence 17 VIII.Banana shape deformation 19

FIG. 1 :
FIG. 1: The t-time evolution of the mass function of two collapsing stars.The upper has central sparser initial configuration, the downer has a central denser one.Different colour curves family denotes the mass function evolution inside three concentric spheres of different radius.

16 )FIG. 1 FIG. 2 :FIG. 3 :
FIG.1displays their latter time evolution following from equations (II.4)-(II.5).From the figure we can easily see that different radius concentric shells need different tduration to reach their horizon size.The duration dif- (III.8)-(III.9)will become the function (III.4)-(III.6)with fixed ϱ routinely.The mass function M [ϱ] involved in (III.1) and (III.4)-(III.6)are completely arbitrary at the classic level.Its physics meaning is the initial mass function of the collapsing star on the co-moving coordinate ϱ line.FIG.5 displays the τ -time evolution of three representative collapsing stars characterised by different M [ϱ].In the left most

FIG. 4 :
FIG.4:The upper part is the potential and phase diagram of oscillators driven by a linearly inverse potential or square inverse force, the lower part is its oscillation function.

FIG. 6 :
FIG. 6: To a co-moving local observer which uses proper time τ to describe physic evolutions, the Schwarzschild singularity caused by gravitational collapse is trivially traversable as the horizon is.So its world line in the τ, ρ plane A,B,C,D• • • ,A ′ ,B ′ ,C ′ ,D ′ is periodic.The Schwarzschild singularity represents not any terminal of physic evolution.But the definition domain of proper time τ covers that of the outside fixed position observers' time t infinite times.During the infinitely long τ -time evolution, the shell structure a freely falling co-moving observer see would not be exactly periodic due to the many-body essence of the interactions among the different shells.This would cause ergodicity of the shell structure's evolution.

FIG. 8 :
FIG. 8: Equivalence between two types of observers microscopic state definition for the black holes.Relative to FIG.3, the right hand part of this figure adopts the extended rule of conventional Penrose Carter diagram so that the east world and west world of the spacetime are simultaneously displayed in the {τ, ρ} plane.

ν +ū co µν 3 ≡0
TABLE II: Two types of observers' measurement of the pressure of a dust type mass m test spherical shell freely falling towards the horizon of a pre-existing Schwarszchild black hole of mass M .The upper is in the Boyer-Lindquist coordinate, the downer is in the Lemaitre coordinate.As a physical reality, the test mass shell's pressure is invariant under the general coordinate transformation, but observer dependent.T µν is the energy momentum tensor of the test mass shell; u µ fx and u µ fr are the four velocity of the fixed position observer and the freely falling observer respectively.T µν , ūµ fx and ūµ fr are their counter parts in the Lemaitre coordinate system.
FIG. 11: The Face to Face Communication Difficulty (F2FCD).In an spatial flat and homogeneous Friedman-Robertson universe, Alice and Bob are a couple who date frequently.Behind each of them is their ex-boy/girl friend Robert and Walker respectively, two spherical giants of radius H −1 0 and average mass density

FIG. 12 :
FIG.12: When two black holes caused by gravitational collapse merge, both participants are outside observers of the other.So none of them can see the other's horizon as physical realities, instead each of them will receive light signal and thus gravitation from the matter contents of the other.The inhomogeneous back reaction of gravitational wave radiation would cause non-uniform orbit decay thus leading to banana shape-deformation in their matter contents' spatial distribution.The resulting black holes' quadrupole can be written as that of the two arc shape strings.

FIG. 13 :
FIG.13: When two black holes defined as exact solutions to the vacuum Einstein equation merge, both participants have absolute and apriori horizon.None of them can see the other's inside horizon structure.The back reaction force of the gravitational wave radiation can only affect the orbit of each other as a whole.None of them will experience innermass distribution's spatial change.Only their horizons will experience droplet shape deformation and merge into a single one finally.Considering the inhomogeneity of the time delaying effect, the true horizon shape deformation[151] would be two comma like way of '' .

ℓ 2
FIG.14:The ratio of self gravitational potential energy of ellipsoid and banana shape deformed star over that of the Schwarzshild black hole with equal surface potential −GM/a = − 1 2 and x-directional size.The y,z-directional size are set to b = c = aη.The right is for banana shape deformation; ψ is the bending angle opened by the arc of the banana relative to the radius of inspiral motion.

FIG. 21 :
FIG. 21: The upper panel is the hamiltonian H[a, ν]'s dependence on a for four typical values of ν.During the merger process of two black holes, the relative obit radius a decreases continuously.Setting a = a low f,lim will leads to an upper bound on the real part of quasi normal frequency ω re 022 ; setting a = a upp f,lim will give a lower bound on ω re 022 Thismeans that the shell recombination scheme B also contributes only one microscopic state to the system.The diagram representation of this state can be plotted as follows 85′ 1 max < 7.85, we find that n ′ 1 = 4 are allowed, other values are not.