REVIEW 4 major objections 5 minor 33 references
Unitary Black hole radiation: Schwarzschild-global monopole background
T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper argues that black hole radiation from a collapsing shell in a Schwarzschild–global monopole background is unitary, checked by both the density-matrix trace and probability conservation.
desk verdict The GM-shell calculation is real work, but the 'unitarity' claim is baked into the setup; the paper would benefit from a referee who asks what is actually being traced over. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a time-dependent simple harmonic oscillator equation for the scalar field mode $b$, obtained from the Wheeler–DeWitt minisuperspace quantization of the shell-plus-scalar system: $$\left[-\frac{1}{2\$\alpha$}\frac{\$partial^{2}$}{\partial $b^{2}$}+\frac{1}{2}\$\alpha$\$omega^{2}$(\tilde{\eta})$b^{2}$\right]\psi(b,\tilde{\eta})=i\frac{\partial\psi}{\partial\tilde{\eta}}.$$ The unitarity argument is carried by the auxiliary function $\zeta(\tilde{\eta})$, which solves $\zeta_{\tilde{\eta}\tilde{\eta}}+\omega^2(\tilde{\eta})\zeta=1/\zeta^3$; its Bessel-function solution fixes the mode coefficients and hence the density-matrix trace. The second line of proof uses the identity $\nabla_\mu J^\mu=\partial_t|\psi|^2=E\,\partial_{\tilde{\eta}}|\psi|^2$, which vanishes as $E\to0$ in the incipient limit.
What would settle it
Recompute $\operatorname{Tr}(\hat{\rho}_f)$ from eqn (70) by integrating the defining equation for $\zeta$ numerically over the full time dependence $E=e^{-\epsilon t/R_{GM}}$, without the Bessel-form ansatz; any deviation from 1 would falsify the unitarity claim. A quicker check is to verify the dimensions of $\tilde{\eta}$ in eqn (47) and to evaluate $u_0$ at late times where $\tilde{\eta}>1$: the square root becomes imaginary, which invalidates the reality assumption behind eqn (97).
Extended reading notes
Core claim
The paper demonstrates unitarity for a global-monopole black hole by two independent routes. In the incipient limit $R\to R_{GM}=2M/(1-\eta^2)$, the time-dependent frequency of the scalar wave functional is $\Omega(t)=e^{-\epsilon t/(2R_{GM})}\omega_0$ with $\epsilon=1-\eta^2$, and the exact Bessel-function solution for the auxiliary parameter $\zeta$ yields transition coefficients $c_n$ that vanish for odd $n$. Summing the even-$n$ probabilities gives $\operatorname{Tr}(\hat{\rho}_f)=1$; because the final density matrix is idempotent, $\hat{\rho}_f^2=\hat{\rho}_f$, this also gives $\operatorname{Tr}(\hat{\rho}_f^2)=1$. Independently, with the spatial probability current vanishing and $E=1-\eta^2-2M/R\to0$, the conservation law $\nabla_\mu J^\mu=0$ follows. Both results hold in the incipient limit, i.e., as the shell approaches the horizon without yet crossing it.
Load-bearing premise
The load-bearing premise is that $\tilde{\eta}$, defined by the time integral in eqn (47), can be used as a dimensionless variable in the Bessel argument $u_0=2\omega_0\sqrt{1-\tilde{\eta}}$, with $\zeta$ and $\zeta_{\tilde{\eta}}$ real; if these assumptions fail, the trace simplification $\operatorname{Tr}(\hat{\rho}_f^2)=1$ no longer follows.
Editorial extensions
If this is right
- If the central claim is correct, an asymptotic observer sees pure-state radiation throughout the collapse up to horizon formation, so no information is lost in this phase.
- Setting $\eta=0$ recovers the Schwarzschild unitarity result, making the analytic proof a strict generalization of earlier numerical work.
- Because the calculation uses only the general form of the metric, the same machinery can be applied to other spherically symmetric static line elements of the type $ds^2=-f(r)dt^2+f(r)^{-1}dr^2+r^2d\Omega^2$.
- The analytic rather than purely numerical confirmation strengthens the case that unitarity survives in non-asymptotically flat black hole spacetimes.
Reading between the lines
- Inference: the paper's unitarity result is established in the incipient limit before the horizon forms; extending it to the full evaporation phase is a separate step that the paper does not claim to make.
- Inference: since the proof depends only on the metric's functional form, a direct test would be to repeat the trace calculation with a different radial function $f(r)$ that still vanishes at a horizon, such as a general power-law deficit, and see whether $\operatorname{Tr}(\hat{\rho}_f^2)=1$ persists.
- Inference: the probability-conservation check is essentially a statement about the incipient limit $E\to0$; a sharper unitarity probe would be to compute the von Neumann entropy $S=-\operatorname{Tr}(\hat{\rho}\ln\hat{\rho})$ of the radiation state during the collapse.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that black hole radiation from an infinitesimally thin collapsing shell carrying a global monopole charge, in a Schwarzschild-global monopole background, is processed by unitary evolution. The authors use the Wheeler-DeWitt formalism to obtain a Schrödinger-like equation for a scalar field mode, solve it in the incipient limit R -> R_GM, and then test unitarity by two methods: showing that Tr(rho_f) = 1 and Tr(rho_f^2) = 1 for the final density matrix, and showing that the probability current is conserved, nabla_mu J^mu = 0, in the incipient limit. The paper concludes that black hole radiation is unitary in a non-asymptotically-flat spacetime and that this result, together with the authors' earlier work on Reissner-Nordstrom, settles unitarity for a class of spherically symmetric backgrounds.
Significance. If the claim were established, the paper would provide an analytic extension of previous unitarity results for black hole radiation to a non-asymptotically-flat background, complementing the numerical work of Saini and Stojkovic and recovering the Schwarzschild case as a limit. The manuscript is clearly organized, and the mode-expansion calculation is carried out in explicit detail. However, the advertised checks do not address the information-loss problem as posed in the literature. The density-matrix check establishes purity of a globally defined wavefunction, not unitarity of the reduced state of the outgoing radiation after tracing over unobserved degrees of freedom. The probability-conservation check is trivially satisfied in the limit considered. Because the central physical question is never directly confronted, the significance of the result, even if the algebra were correct, is considerably weaker than claimed.
major comments (4)
- [V.B, Eq. (77)] The demonstration that Tr(rho_f^2) = 1 is a consequence of the construction of rho_f, not of any property of black hole radiation. In Eqs. (67)-(68), rho_i and rho_f are projectors built from the expansion coefficients l_n and c_n of the global wavefunction psi(b,t). Because Eq. (63) is an expansion in a complete orthonormal set, the normalization sum_n |c_n|^2 = 1 used in Eq. (72) follows from the solution of the Schrodinger equation (50), which is unitary by assumption. Equation (72) is the identity |psi><psi|^2 = |psi><psi| for any normalized state. The physically relevant object for the information-loss problem is the reduced density matrix of the outgoing modes after tracing over interior modes or black-hole degrees of freedom; no such trace is computed anywhere in the manuscript. Therefore Eqs. (71)-(73) do not support the claim of unitary black hole radiation.
- [V.B, Eq. (77)] The probability-conservation check is not independent and is trivial in the limit analyzed. Since the spatial current is set to zero because b is independent of spatial coordinates, Eq. (76) reduces to partial_t |psi|^2. Using partial_t eta_tilde = E and taking E -> 0 as R -> R_GM, Eq. (77) merely states that the time derivative of the global wavefunction norm vanishes in the incipient limit. This is a consequence of the assumed unitary Schrodinger evolution (50), not a physical condition on the outgoing radiation. The apparent 'second independent line of approach' is therefore not independent and does not test unitarity of the radiation state.
- [Appendix, Eqs. (91)-(97)] The explicit computation of Tr(rho_f) rests on an invalid premise. The parameter eta_tilde is defined in Eq. (47) as an integral over dt and therefore has dimensions of time, yet it appears in u0 = 2 omega0 sqrt(1 - eta_tilde) in Eq. (93); the combination 1 - eta_tilde is dimensionally inconsistent. Moreover, for eta_tilde > 1, which occurs at late times in the incipient limit, u0 is imaginary, and the asserted realness of zeta and zeta_eta in Eqs. (91)-(96) is neither proved nor evident. Equation (97) is evaluated with Mathematica and reduces to 1 only under this realness assumption. Because the appendix's Bessel-function solution is the basis for the 'analytic' confirmation in Eq. (98), this is a load-bearing gap in the only direct evaluation of the density-matrix trace.
- [VI and Eq. (54)] The scope of the result is overstated. The calculation is confined to the incipient limit R -> R_GM with E -> 0, and the time dependence of E is imposed by hand in Eq. (54) with a cut-off time t_f. The statement in Section VI that 'if unitarity is preserved in this limit, it should be valid at every instant of time' is an unsupported extrapolation. No argument is given that the incipient limit is representative of the full collapse, nor that the chosen early-time and cut-off behavior does not affect the conclusion. The conclusion should be limited to the specific model and limit analyzed, and the extrapolation should be either removed or derived.
minor comments (5)
- [Eq. (55)] The exponent in Eq. (55) contains an undefined mass parameter m: the term e^{-m omega0 b^2/2} should likely read e^{-alpha omega0 b^2/2} to be consistent with Eq. (56). Please correct the typo.
- [Eq. (66) and Eq. (88)] The expressions for c_n in Eq. (66) and Eq. (88) disagree in the prefactor: one contains (Omega_f zeta^2)^{1/4}, while the other contains 1/(Omega_f zeta^2)^{1/4}. This discrepancy propagates into the trace formula (70) and should be reconciled before the computation can be considered consistent.
- [Eq. (47) and throughout] The notation eta_tilde for the new time parameter is easily confused with the global monopole charge eta, especially because both appear in the same equations. Please choose a distinct symbol, such as tau or s, for the integrated time variable.
- [Section IV.A] There is a literal placeholder in the text 'are[ ? ] just the simple harmonic oscillator ground states'; the citation is missing. Please supply the proper reference or remove the placeholder.
- [Section IV.A, Eq. (37)] The mode expansion (37) is introduced on the basis of separability of Eqs. (35) and (36), but the step is not shown. Since Eq. (35) contains T_tt/T_t^3 and Eq. (36) has r-dependent coefficients, the reader cannot verify that the same set of functions f_k(r) diagonalizes both kinetic and potential terms. Please provide the separation argument or state explicitly what assumptions on f_k are needed.
Circularity Check
The unitarity conclusion is tautological: ρ_f is constructed as the projector onto a normalized global wavefunction, so Tr(ρ_f²)=1 is guaranteed by construction; the probability-conservation 'check' vanishes only because E→0.
-
self definitional
[Section V.A, Eqs. (67)-(73), especially Eq. (72)]
"ρ̂f = Σ_{m,n} c_m c^*_n |φ_m⟩⟨φ_n| ... ρ̂f^2 = ... = ρ̂f (as, (Σ_n |c_n|^2)=1 by eqn(71)) ... Therefore, by eqn(72) we get, Tr(ρ̂f^2)=Tr(ρ̂f)=1."
The final density matrix is defined directly from the expansion coefficients c_n of the single global wavefunction ψ(b,t) (Eqs. 63-68). For any normalized pure state, ρ_f²=ρ_f and Tr(ρ_f²)=1 automatically; the only input used in Eq. (72) is Σ|c_n|²=1, i.e., wavefunction normalization. No trace is taken over interior modes, modes behind the horizon, or black-hole degrees of freedom, so the 'unitarity' test is exactly the assumption that a global Schrödinger wavefunction stays pure, not an independent demonstration that outgoing radiation is unitary.
-
other
[Section V.B, Eq. (77)]
"∇µJ µ = ∂|ψ|^2 /∂t = ∂|ψ|^2 /∂η~ ∂η~ /∂t =E∂|ψ|^2 /∂η~ ... For, R →RGM, ∇µJ µ = 0 ( as, E → 0) (77)."
This purported conservation check is not derived from the Schrödinger equation or from a continuity equation; it is simply the statement that ∂|ψ|²/∂t vanishes because the prefactor E=1−η²−2M/R tends to zero in the incipient limit. Any function of η~(t) would satisfy this, so the check tests the freezing of the time reparametrization, not unitarity. It is therefore a trivial consequence of the already-assumed limit, not an independent confirmation.
full rationale
The paper's central unitarity claim reduces to its own input by construction. The density-matrix calculation defines ρ_f from the coefficients of a single pure-state wavefunction; Eq. (72) then uses only Σ|c_n|²=1, which is wavefunction normalization, to prove ρ_f²=ρ_f and hence Tr(ρ_f²)=1. This would hold for any normalized ket and does not address the information-loss question, which requires a reduced density matrix after tracing over the interior or over black-hole modes; no such trace is computed. The second 'independent' check, ∇_μ J^μ=0, is also trivial in the incipient limit because the factor E multiplies the time derivative and E→0; it does not establish probability conservation in any dynamical sense. The authors' adoption of the formalism from their own Ref. [7] is present, but the main circularity is visible in the present paper's own equations without needing to rely on that citation. The appendix's Bessel-function calculation has a separate dimensional inconsistency (u0=2ω0√(1−η~) with η~ from Eq. (47) having dimensions of time), but that is a correctness defect rather than the source of the circularity; even if the appendix were corrected, Tr(ρ_f)=1 would remain just the normalization of the expansion coefficients. Overall, the result is forced by the definition of the pure-state density matrix, giving a circularity score of 8.
Assumptions & free parameters
free parameters (1)
- cut-off time t_f =
t_f → ∞
assumptions (5)
- domain assumption The Wheeler-DeWitt minisuperspace formalism yields a Schrödinger-like equation (50) for the scalar field modes.
- standard math The evolution equation (50) is unitary, i.e., the Hamiltonian is Hermitian with real frequency ω.
- domain assumption The spectral theorem applied to the infinite matrices A and B, and the reduction to a single eigenvector b, captures the full field behavior.
- ad hoc to paper The incipient limit (R→RGM) result extends to all times: 'if unitarity is preserved in this limit, it should be valid at every instant of time' (Section VI).
- ad hoc to paper The Bessel function combinations ζ and ζ_η are real for the entire range of the time parameter, including when u0 is imaginary.
Cite this review
Pith. "Pith review of Unitary Black hole radiation: Schwarzschild-global monopole background." pith.science (2026). https://pith.science/paper/2VWEZ7HY
@misc{pith2026190809616,
author = {Pith},
title = {Pith review of: Unitary Black hole radiation: Schwarzschild-global monopole background},
year = {2026},
howpublished = {\url{https://pith.science/paper/2VWEZ7HY}},
note = {Machine review of arXiv:1908.09616}
}
read the original abstract
Black hole radiation from an infinitesimally thin massive collapsing shell, possessing a global monopole charge, which in turn leads to a Schwarzschild black hole with a global monopole charge has been shown to be processed by a unitary evolution. The exterior metric of the collapsing shell is described by the global monopole (GM) metric. The analysis is performed using the Wheeler-deWitt formalism which gave rise to a Schr\"{o}dinger-like wave equation. Existence of unitarity is confirmed from two independent lines of approach. Firstly, by showing that the trace of the square of the density matrix, of the outgoing radiation, from a quantized massless scalar field, is unity. Secondly, by proving that the conservation of probability holds for the wave function of the system.
Reference graph
Works this paper leans on
-
[1]
are unrealistic and perhaps appear in some instances of cosmic phase transition [ 10]. 2 The surface gravity κ for the metric as given in eqn(1) is obtained by noting that the metric is of the form [ 15], ds2 = −fdt 2 +f −1dr2 +r2dΩ 2 2, (6) implying,κ = f ′ (r) 2 , (7) implying,κ GM = ( 1 −η2) 2 4M (8) ( as, f(r) = 1 −η2 − 2M r ) . whereκGM is the surfac...
-
[2]
was done in [ 11] and it was show that the outgoing Hawking radiation is thermal possessing a Planck spectrum, N = 1 e8πM ω/(1−η2)2 − 1 (9) where N is the number density of outgoing quanta of particles. The Hawking temperature is recovered to be, TGM = ( 1 −η2) 2 8πM , (10) which can also be obtained from eqn(8) using the Hawk- ing relation TH = κ 2π (whi...
-
[3]
Why Black Hole Information Loss is Paradoxical
D. Wallace, Why Black Hole Information Loss is Para- doxical, gr-qc/1710.03783v2 (2017)
work page Pith review arXiv 2017
-
[4]
S. D. Mathur, Class. Quant. Grav. 26, 224001 (2009)
work page 2009
-
[5]
Hawking, Commun
S.W. Hawking, Commun. Math. Phys. 43, 199 (1975)
1975
-
[6]
S. W. Hawking, Phys. Rev. D 14, 2460 (1976)
1976
- [7]
-
[8]
DeWitt, Phys
B.S. DeWitt, Phys. Rev. 160, 1113 (1967)
1967
Show all 33 references
-
[9]
Polchinski, The Black Hole Information Problem , hep- th/1609.04036v1 (2016)
J. Polchinski, The Black Hole Information Problem , hep- th/1609.04036v1 (2016)
2016 arXiv
-
[10]
D. N. Page, Phys. Rev. Lett. 71, 3743-3746 (1993)
1993
-
[11]
Dadhich, K
N. Dadhich, K. Narayan and U. A. Yajnik, Pramana - J Phys 50: 307 (1998)
1998
-
[12]
Saini and D
A. Saini and D. Stojkovic, Phys. Rev. D 97, 025020 (2018)
2018
-
[13]
Superspace and the nature of quantum geometrodynamics
J. A. Wheeler, “Superspace and the nature of quantum geometrodynamics” in Batelle Recontres, Benjamin, New York (1968). 10
1968
-
[14]
Barriola and A
M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63, 341 (1989)
1989
-
[15]
Poisson, A Relativist’s Toolkit: The Mathematics of Black-Hole Mechanics , Cambridge University Press, 2007
E. Poisson, A Relativist’s Toolkit: The Mathematics of Black-Hole Mechanics , Cambridge University Press, 2007
2007
-
[16]
Vachaspati, D
T. Vachaspati, D. Stojkovic and L.M. Krauss, Phys. Rev. D 76, 024005 (2007)
2007
-
[17]
Saini and D
A. Saini and D. Stojkovic, Phys. Rev. Lett. 114, 111301 (2015)
2015
-
[18]
K. G. Zloshchastiev, Phys. Rev. D. 57, 4812 (1998)
1998
-
[19]
C. A. L´ opez, Phys. Rev. D 38, 3662 (1988)
1988
-
[20]
Greenwood, JCAP 1001 002 (2010)
E. Greenwood, JCAP 1001 002 (2010)
2010
-
[21]
Israel, Nuovo Cimento 44B, 1 (1966)
W. Israel, Nuovo Cimento 44B, 1 (1966)
1966
-
[22]
Israel, Nuovo Cimento A 51: 744 (1967)
W. Israel, Nuovo Cimento A 51: 744 (1967)
1967
-
[23]
Then we have, Hshell = [(EΠ shell)2 +E(4πµR2)2]1/2 ≡ [q2 +m2]1/2, (27) where q2 := (EΠ shell)2 and m2 :=E(4πµR2)2
and eqn(24) we note that, as R →RGM , Π shell = 4πµR2Rt√ E √ E2 −R2 t , (25) Hshell = 4πE 3/2µR2 √ E2 −R2 t , (26) where,µ :=σ ( 1 − 2πσRGM − η2 8πσRGM ) . Then we have, Hshell = [(EΠ shell)2 +E(4πµR2)2]1/2 ≡ [q2 +m2]1/2, (27) where q2 := (EΠ shell)2 and m2 :=E(4πµR2)2. Hshell...
-
[24]
I. A. Pedrosa, J. Math. Phys. 28, 2662 (1987)
1987
-
[25]
C. M. A. Dantas, I. A. Pedrosa and B. Baseia, Phys. Rev. A 45, 1320 (1992)
1992
-
[26]
H. R. Lewis, J. Math. Phys. 9, 1976 (1968)
1968
-
[27]
H. R. Lewis. and W. B. Riesenfeld, J. Math. Phys. 10, 1458 (1969)
1969
-
[28]
Pal and N
S. Pal and N. Banerjee, Phys. Rev. D 91 044042 (2015)
2015
-
[29]
Kolopanis and T
M. Kolopanis and T. Vachaspati, Phys. Rev. D 87 085041 (2013)
2013
-
[30]
Sakurai and J
J.J. Sakurai and J. Napolitano, Modern Quantum Me- chanics, 2nd Edition , Addison-Wesley, 2011
2011
-
[31]
Pal and N
S. Pal and N. Banerjee, Phys. Rev. D 90, 104001 (2014)
2014
-
[33]
Pal and N
S. Pal and N. Banerjee, J. Math. Phys. 57, 122502 (2016)
2016
-
[70]
(71) eqn(71) shows that the necessary condition for the uni- tary evolution of states holds
to obtain (see appendix), Tr (ˆρf ) = 1. (71) eqn(71) shows that the necessary condition for the uni- tary evolution of states holds. For the sufficient condi- tion, we compute Tr (ˆρ2 f ). From eqn( 68), ˆρf = ∑ m,n cmc∗ n|φm⟩⟨φn| leading to, ˆρ2 f = (∑ m,n cmc∗ n|φm⟩⟨φn| ) ...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.