REVIEW 2 major objections 4 minor 2 cited by
Black holes in thermal bath live shorter: implications for primordial black holes
T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read A black hole immersed in a thermal bath emits more than Hawking's Planckian formula predicts, so black holes in the hot early universe evaporate faster and primordial black holes live shorter.
desk verdict The finite-temperature spectrum is a plausible detailed-balance result, but the paper's central claim that hot baths shorten PBH lifetimes is unsupported because the mass-loss equation drops the ingoing flux, which dominates in the regime they advertise. 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 engine of the derivation is the thermofield-double (TFD) formalism, which doubles the Fock space and defines a temperature-dependent pure state $|0,\beta\rangle$ whose reduced density matrix is the thermal ensemble at inverse temperature $\beta$. Thermal creation and annihilation operators act on the zero-temperature vacuum, so the field on $I^-$ is decomposed in thermal modes. A Bogoliubov transformation between the $I^-$ thermal modes and the $I^+$ outgoing modes, using the eikonal coefficient ratio $|\beta_{\omega\omega'}|^2 = e^{-2\pi\omega/\kappa}|\alpha_{\omega\omega'}|^2$, produces the extra term $(e^{2\pi\omega/\kappa}+1)/(e^{\beta\omega}-1)$ in the bosonic spectrum and the corresponding Fermi term. This correction term is what carries the bath temperature into the decay rate and ultimately into the black hole lifetime.
What would settle it
A concrete test is to compute the full mass-loss rate including the absorption cross-section of the ingoing thermal flux for $T_b \geq T_{\mathrm{BH}}$. If detailed balance forces the net mass change to become positive (accretion) in the hot-bath regime the paper emphasizes, or if a full emission-plus-absorption calculation shows the lifetime is not shortened, the central claim is falsified.
Extended reading notes
Core claim
Working with a massless scalar (and a spin-1/2 fermion) field prepared in a thermal state on past null infinity and evolved through a collapsing spherically symmetric geometry, the paper obtains the bosonic number density per mode $$n_\omega = \frac{\Gamma_\omega}{$e^{{2\pi\omega/\kappa}}$-1}\left[1+\frac{$e^{{2\pi\omega/\kappa}}$+1}{$e^{{\beta\omega}}$-1}\right],$$ with $\beta=1/T_b$, and a fermionic analogue with the signs in the exponentials flipped. In the limit $\beta\to\infty$ the bracket tends to 1 and the standard Hawking spectrum is recovered; at finite bath temperature the spectrum is a function of both temperatures and is not Planckian. Feeding this spectrum into the standard decay equation $dM/dt = -(M_p^4/M^2)\sum_i g_i \epsilon_i$, the enhanced emission makes bosonic contributions grow as $T_b$ when the bath is hot, so the lifetime scales as $\tau_{\mathrm{BH}}\propto 1/T_b$. For a 10 g primordial black hole at $T_b/T_{\mathrm{BH}}=10^2$ the lifetime drops from $\sim 10^{-24}$ s to $\sim 10^{-29}$ s, and in a two-stage reheating model with a time-dependent bath the lifetime of a $\sim 10^2$ g PBH is shortened by up to a factor of $10^2$ relative to the zero-temperature case.
Load-bearing premise
The load-bearing premise is that the thermal bath changes only the outgoing radiation spectrum; the paper explicitly neglects the ingoing flux of bath particles and its backreaction, so in a bath hotter than the black hole the neglected absorption could counter or outweigh the enhanced emission and the 'shorter lifetime' conclusion could fail.
Editorial extensions
If this is right
- Finite bath temperature must be treated as a first-order correction to black hole evaporation, alongside greybody factors and kinematic cut-offs.
- Very light primordial black holes formed right after inflation evaporate earlier than the Hawking-only estimate, shifting constraints and observable signatures tied to their evaporation such as gamma-ray bursts, relic particle production, and entropy injection.
- In the hot-bath limit the bosonic decay coefficient grows as $T_b$, so the lifetime scales as $\tau\propto 1/T_b$, making the enhancement strongest for the lightest PBHs with the highest Hawking temperature.
- In the two-stage reheating scenario there is a specific initial PBH mass that decays fastest, because the bath temperature falls while the Hawking temperature rises as the hole shrinks.
Reading between the lines
- The paper explicitly neglects ingoing bath flux and its backreaction; a full emission-plus-absorption calculation for $T_b \geq T_{\mathrm{BH}}$ could show accretion dominating evaporation in the hottest regime, which would reverse the 'shorter lifetime' conclusion there.
- If the corrected spectrum holds, applying the same thermofield-double construction to rotating or charged black holes would yield a bath-temperature dependence shaped by a different eikonal coefficient ratio.
- The boson-fermion asymmetry (bosonic $\epsilon_i \propto T_b$ while fermionic $\epsilon_i$ stays constant) is a distinctive fingerprint: in a hot bath, photon and graviton emission should be enhanced relative to neutrinos, which could show up in the evaporation products of light PBHs.
- A laboratory analogue — a radiating body or analogue horizon coupled to a tunable heat bath — could test the predicted non-Planckian correction factor directly, since the ratio to the zero-bath spectrum is known.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives a finite-temperature correction to the Hawking spectrum by quantizing a massless scalar or fermion in a thermofield-double state at inverse temperature β on I^- and evolving it through a collapsing black-hole spacetime. The result, Eqs. (15)-(16), is a non-Planckian outgoing spectrum containing a stimulated-emission factor that depends on the ratio of black-hole and bath temperatures. The authors then substitute this spectrum into the Page mass-loss formula, Eq. (17), and integrate the PBH evolution in a reheating cosmology to claim that PBHs in a thermal bath decay faster, reducing lifetimes by orders of magnitude (Figs. 1-2). The caveat that ingoing classical flux and backreaction are neglected is stated explicitly in the Summary.
Significance. The TFD calculation leading to Eq. (15) is a clean, parameter-free derivation and, if correct, would be a useful addition to the finite-temperature QFT-in-curved-spacetime literature. I also note that the cosmological parameters in Sec. 6 are scanned rather than fitted to the desired outcome, so there is no circularity in the numerical application. However, the paper's headline application to PBH lifetimes does not follow from the derived spectrum, because Eq. (17) omits the ingoing bath flux and its absorption by the black hole. In the regime T_b > T_BH emphasized by the authors, the omitted absorption term can dominate and even reverse the sign of dM/dt. The claimed shorter PBH lifetimes are therefore unsupported, and the paper's significance for PBH cosmology is not established.
major comments (2)
- [Eq. (17) and following paragraph; Summary] The mass-loss equation (17) is Page's emission formula and contains no term for energy absorbed from the thermal bath; the authors explicitly neglect 'any ingoing classical flux or its backreaction towards the black hole dynamics' in the Summary. In the regime highlighted in Fig. 1 (T_b/T_BH = 10^2), the ingoing thermal flux at the horizon scales as T_b^4, while the stimulated-emission part of Eq. (20) grows only as T_b for bosons and is constant for fermions. Energy conservation requires dM/dt = F_abs(T_b,M) - F_em(T_BH,T_b); for a sufficiently hot bath this quantity becomes positive, so the black hole accretes rather than evaporating faster. The advertised conclusion that black holes in a thermal bath live shorter therefore does not follow from Eqs. (17)-(20).
- [Eq. (15) and Eq. (17)] The spectrum (15) is the outgoing number density on I^+, but it includes bath quanta that are scattered by the black-hole potential barrier and never cross the horizon. Substituting this n_out into the Page equation (17) treats those reflected quanta as mass loss and double-counts the radiation that is not extracted from the black hole. The mass-loss rate should instead be obtained from the net energy flux across the horizon, with the absorption fraction Γω in Eq. (12) controlling the difference between ingoing and outgoing fluxes, rather than being applied only to the emission side.
minor comments (4)
- [Eq. (14)] Equation (14) defines Nω with an overall minus sign, while the final spectrum in Eq. (15) is positive; if this is only a typo, it should be removed for consistency.
- [Before Eq. (20)] The phrase 'geometrical optical limit (ω ≫ MADM)' is dimensionally inconsistent; the geometric optics limit should be stated as ωM ≫ 1 (or ω ≫ 1/M).
- [Sec. 6 (reheating model)] The two-stage reheating model with α1 = 0.15, α2 = 0.95 and T_re = 0.1 GeV is introduced without a sensitivity study, so the quantitative lifetime reductions in Fig. 2 should be presented with parameter scans or uncertainty bands before being used as predictions.
- [Abstract and Summary vs. Sec. 2] The Abstract and Summary describe the setup as emitted radiation being 'thermalised at the bath temperature,' but the calculation instead assumes a thermal state on I^- from the start; these two physical pictures should be reconciled in the text.
Circularity Check
No circularity: the modified spectrum follows from an assumed thermal initial state plus standard Hawking Bogoliubov coefficients, and the shorter-lifetime claim is a direct evaluation of the Page mass-loss equation with that spectrum.
full rationale
The paper's central result, Eqs. (15)-(16), is derived by evolving a TFD thermal state with inverse temperature β (an assumed initial condition) through the standard Bogoliubov transformation and evaluating the I^+ number operator; the Hawking-temperature dependence enters through the eikonal relation |β_ωω'|^2 = e^{-2πω/κ}|α_ωω'|^2 (Eq. 13), which is an external input from Hawking, not a fit to the target lifetime. The bath temperature T_b is an input parameter of the thermal state, and the nontrivial content is how it combines with κ; this is a derivation, not a renaming. The shorter-lifetime conclusion follows by substituting the derived spectrum into the independently cited Page mass-loss equation (Eq. 17) and integrating; no parameter is tuned to reproduce the claimed reduction from 10^-24 s to 10^-29 s or the τ_BH ∝ T_b^{-1} scaling. The self-citations present ([30,31] for reheating temperature behavior; [2] for the relation between outgoing flux and horizon decrease) are peripheral: the fixed-temperature result in Fig. 1 and the analytic scaling do not rest on them. The paper's explicit limitation, 'we neglect the effects of any ingoing classical flux or its backreaction towards the black hole dynamics,' is a physical-correctness caveat about the model, especially in the T_b > T_BH regime, but omitting absorption is an assumption rather than a circular reduction: the plotted lifetimes are consequences of the assumed evolution equation, and an added absorption term would modify the model, not reveal that the claimed prediction was secretly an input. No circular step can be exhibited from the paper's own equations.
Assumptions & free parameters
free parameters (5)
- M_in (initial PBH mass)
- T_in (bath temperature at PBH formation)
- Reheating temperature T_re =
0.1 GeV
- Reheating exponents alpha_1, alpha_2 =
0.15, 0.95
- Relativistic degrees of freedom g_* =
106.75
assumptions (5)
- domain assumption The thermal initial state on I^- is exactly the TFD thermal vacuum at inverse temperature beta (Eq. 1 and Eq. 3).
- domain assumption The eikonal Bogoliubov relation |beta_omega_omega'|^2 = e^(-2 pi omega / kappa) |alpha_omega_omega'|^2 (Eq. 13).
- ad hoc to paper Ingoing bath flux and its backreaction are neglected in the PBH mass evolution equation (Eq. 17).
- ad hoc to paper The bath temperature during reheating follows T_b = T_in (a/a_in)^-alpha with alpha_1 = 0.15, alpha_2 = 0.95 and T_re = 0.1 GeV.
- domain assumption All emitted Standard Model particles are treated as massless (z_i -> 0) with g_* approximately 106.75.
Cite this review
Pith. "Pith review of Black holes in thermal bath live shorter: implications for primordial black holes." pith.science (2026). https://pith.science/paper/57TAJTRB
@misc{pith2026250111925,
author = {Pith},
title = {Pith review of: Black holes in thermal bath live shorter: implications for primordial black holes},
year = {2026},
howpublished = {\url{https://pith.science/paper/57TAJTRB}},
note = {Machine review of arXiv:2501.11925}
}
read the original abstract
Hawking radiation from a non-extremal black hole is known to be approximately Planckian. The thermal spectrum receives multiple corrections including greybody factors and due to kinematical restrictions on the infrared and ultraviolet frequencies. We show that another significant correction to the spectrum arises if the black hole is assumed to live in a thermal bath and the emitted radiation gets thermalised at the bath temperature. This modification reshapes the thermal spectrum, and leads to appreciable deviation from standard results including modification in the decay rate of black holes. We argue that this altered decay rate has significance for cosmology and, in a realistic setting, show that it alters the life time of primordial black holes (PBHs) in the early universe. In particular, the very light PBHs formed right after the end of inflation decay faster which may have interesting phenomenological implications.
Figures
Forward citations
Cited by 2 Pith papers
-
Revisiting PBH accretion, evaporation and their cosmological consequences
Relativistic accretion onto Kerr primordial black holes gives roughly 4.5x mass growth and fast spin-down, strengthening BBN bounds, lowering the survival mass to ~2.7e14 g, and erasing the high-frequency stochastic g...
-
Thermal bath induced suppression of Hawking radiation
A Schwarzschild black hole coupled bilinearly to an environmental scalar field radiates less than Hawking's prediction, according to a semiclassical Keldysh calculation.
Reference graph
Works this paper leans on
-
[10]
J. D. Bekenstein and A. Meisels, Phys. Rev. D 15, 2775- 2781 (1977)
work page 1977
-
[1]
S. W. Hawking, Commun. Math. Phys. 43, 199-220 (1975), [erratum: Commun. Math. Phys. 46, 206 (1976)]
work page 1975
-
[2]
A. Chatterjee, B. Chatterjee and A. Ghosh, Phys. Rev. D 87, no.8, 084051 (2013)
work page 2013
-
[3]
Israel, Phys
W. Israel, Phys. Lett. A 57, 107-110 (1976)
1976
-
[4]
J. B. Hartle and S. W. Hawking, Phys. Rev. D 13, 2188- 2203 (1976)
work page 1976
-
[5]
G. W. Gibbons and M. J. Perry, Proc. Roy. Soc. Lond. A 358, 467-494 (1978)
work page 1978
-
[6]
R. M. Wald, University of Chicago Press, ISBN 978-0- 226-87027-4
-
[7]
D. N. Page, Phys. Rev. D 13, 198-206 (1976)
work page 1976
Show all 31 references
-
[8]
Cheek, L
A. Cheek, L. Heurtier, Y. F. Perez-Gonzalez and J. Turner, Phys. Rev. D 105, no.1, 015022 (2022) doi:10.1103/PhysRevD.105.015022
2022 doi
-
[9]
D. N. Page, Phys. Rev. D 14, 3260-3273 (1976)
1976
-
[11]
Boonserm and M
P. Boonserm and M. Visser, Phys. Rev. D 78, 101502 (2008). 6
2008
-
[12]
Barcelo, S
C. Barcelo, S. Liberati, S. Sonego and M. Visser, Phys. Rev. D 83, 041501 (2011)
2011
-
[13]
Visser, Int
M. Visser, Int. J. Mod. Phys. D 12, 649-661 (2003)
2003
-
[14]
Visser, JHEP 07, 009 (2015)
M. Visser, JHEP 07, 009 (2015)
2015
-
[15]
M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042- 5045 (2000)
2000
-
[16]
Takahashi and H
Y. Takahashi and H. Umezawa, Int. J. Mod. Phys. B 10, 1755-1805 (1996), Collect. Phenom. 2, 55-80 (1975)
1996
-
[17]
Das, Finite temperature field theory, World Scientific, (2023)
A. Das, Finite temperature field theory, World Scientific, (2023)
2023
-
[18]
G. F. Chapline, Nature 253, no.5489, 251-252 (1975) doi:10.1038/253251a0
1975 doi
-
[19]
B. Carr, F. Kuhnel and M. Sandstad, Phys. Rev. D 94, no.8, 083504 (2016) doi:10.1103/PhysRevD.94.083504
2016 doi
-
[20]
B. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, Rept. Prog. Phys. 84, no.11, 116902 (2021) doi:10.1088/1361- 6633/ac1e31
2021 doi
-
[21]
B. J. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, Phys. Rev. D 81 (2010), 104019 doi:10.1103/PhysRevD.81.104019
2010 doi
- [22]
- [23]
-
[24]
M. He, K. Kohri, K. Mukaida and M. Yamada, [arXiv:2407.15926 [hep-ph]]
-
[25]
M. He, K. Kohri, K. Mukaida and M. Yamada, JCAP 01, 027 (2023) doi:10.1088/1475-7516/2023/01/027
2023 doi
-
[26]
Altomonte, M
C. Altomonte, M. Fairbairn and L. Heurtier, [arXiv:2501.05531 [astro-ph.CO]]
-
[27]
B. J. Carr and S. W. Hawking, Mon. Not. Roy. Astron. Soc. 168, 399-415 (1974)
1974
-
[28]
Bousso and S
R. Bousso and S. W. Hawking, Phys. Rev. D 54, 6312- 6322 (1996)
1996
-
[29]
B. J. Carr, Astrophys. J. 201, 1-19 (1975)
1975
-
[30]
M. R. Haque, D. Maity and R. Mondal, JHEP 09, 012 (2023) doi:10.1007/JHEP09(2023)012
2023 doi
-
[31]
Chakraborty, M
A. Chakraborty, M. R. Haque, D. Maity and R. Mondal, Phys. Rev. D 108, no.2, 023515 (2023) doi:10.1103/PhysRevD.108.023515
2023 doi
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.