REVIEW 4 major objections 5 minor 30 references
Charged Taub-NUT type Black Holes in Einstein-Weyl-Maxwell Theory
T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Charged Taub-NUT black holes in Einstein-Weyl-Maxwell theory split into three disconnected branches, with two meeting at a minimal horizon radius.
desk verdict A plausible new numerical solution family whose central branch-structure claim is unreproducible as written because the Weyl coupling α is never stated. 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 carrying mechanism is the numerical shooting method applied to the three coupled second-order ordinary differential equations (2.9) for the metric functions $h(r)$, $f(r)$ and the electric potential $a(r)$ in the NUT ansatz (2.8). Near the horizon $r_0$, the functions are expanded in Taylor series, and the expansion coefficients are fixed by four parameters; the paper normalizes $a_0=0$, $h_1=1/r_0$, $f_1=1/r_0+\delta$, leaving $\delta$ as the shooting parameter tuned so that the functions approach the required asymptotic behavior at infinity. The sign of $\delta$ organizes the branches: positive $\delta$ gives one branch, negative $\delta$ gives two branches connected at a minimal $r_0$. Temperature follows from surface gravity, $T=\sqrt{h_1 f_1}/(4\pi)$, and entropy from the Wald formula evaluated on the horizon, which simplifies to $S=\pi(r_0^2+n^2)-4\pi\alpha r_0\delta$.
What would settle it
Recompute the solutions with an independent numerical method (for example, a spectral or relaxation solver, or a different shooting implementation with far smaller error tolerance) at $a_1=1.0$, $n=0.5$ and at $a_1=0.02$, $n=0.5$, and count the $\delta$ values satisfying the boundary conditions at infinity. The central claim fails if the count is not three for the charged cases, if the two negative-$\delta$ branches do not meet at a single minimal horizon radius, or if the charged branches intersect for some horizon radius.
Extended reading notes
Core claim
The central discovery is a three-branch structure for charged Taub-NUT black holes in Einstein-Weyl-Maxwell theory. For fixed $n$ and $a_1$, the shooting method finds three values of the slope parameter $\delta$ (defined by $f_1 = 1/r_0+\delta$) that satisfy the asymptotic boundary conditions at infinity: one positive-$\delta$ solution, whose metric function shows a 'bump' and whose mass parameter is negative, and two negative-$\delta$ solutions with positive mass. The two negative-$\delta$ branches have a minimal horizon radius and connect smoothly there, producing 'hook'-shaped temperature and entropy curves. The positive-$\delta$ branch has consistently higher temperature and lower entropy than the other two. The paper emphasizes that the charged solutions are disconnected while the neutral ones ($a_1=0$) intersect: introducing charge separates the two groups, with the separation widening as $a_1$ increases. In the limit $n\to 0$ the solutions reduce to the static spherical charged solutions with two branches, and in the limit $a_1\to 0$ they reduce to the neutral NUT solutions.
Load-bearing premise
The load-bearing premise is that the numerical integrations from the horizon to infinity reliably find every solution branch for the chosen parameter values and the fixed value of the Weyl coupling $\alpha$, since the paper reports no convergence tests, residual error estimates, or an explicit numeric value for $\alpha$.
Editorial extensions
If this is right
- If the central claim is correct, Einstein-Weyl-Maxwell theory provides a concrete example where electric charge changes the solution topology: the two intersecting groups of neutral NUT solutions become two disconnected groups of charged solutions, and the static spherical two-branch structure grows to three branches once the NUT parameter is switched on.
- The existence of a common minimal horizon radius for the two negative-$\delta$ branches implies a limiting smallest black hole for given $(n,a_1)$, with temperature and entropy curves forming a cusp-like 'hook' at that radius.
- Because the positive-$\delta$ branch has negative mass and higher temperature, any complete thermodynamic account of these solutions must include negative-mass configurations that are absent in the static charged case.
- The disconnectedness of the charged branches means there is no continuous family of charged NUT black holes interpolating between the two groups; transitions between them, if they occur, would be discontinuous in the space of solutions.
- These solutions give a concrete setting in which to test definitions of mass, NUT charge, and the first law in higher-derivative gravity, since those definitions are ambiguous for NUT spacetimes and are further complicated by higher-curvature terms.
Reading between the lines
- If the numerical claim holds, the small-charge limit should show the neutral intersection point splitting into a gap whose size scales as a power of $a_1$; computing that scaling perturbatively would give an analytic check of the disconnection claim.
- The three-branch structure may be generic to higher-derivative theories with a massive spin-2 mode rather than specific to Weyl gravity, so similar three-branch NUT solutions could exist in other quadratic-curvature theories; the paper does not consider this.
- A natural testable extension is to compute the free energy and heat capacities of the three branches to see whether phase transitions (for example, Hawking-Page-like transitions) occur among them; the paper leaves the phase diagram to future work.
- The negative mass of the positive-$\delta$ branch combined with the minimal radius of the negative-$\delta$ branches suggests a charge- and NUT-dependent extremality bound that could be mapped numerically.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper numerically constructs charged Taub-NUT-type black hole solutions in four-dimensional Einstein-Weyl-Maxwell theory. Adopting the metric ansatz (2.8) and gauge potential, it derives three second-order ordinary differential equations (2.9), performs near-horizon Taylor expansions (2.10)-(2.11), and uses a shooting method in which the parameter delta in (3.1) is adjusted to satisfy unspecified boundary conditions at infinity. The central claim, stated in the abstract and Section 4, is that for fixed NUT parameter n and electric parameter a1 there are exactly three branches of regular asymptotically flat solutions: one with delta>0 and two with delta<0 that meet at a minimal horizon radius. The paper further claims that the charged branches are disconnected, whereas the neutral branches intersect, and reports temperature and Wald entropy for the branches in Figures 4-16.
Significance. If the reported solutions exist as claimed, they are genuinely new and provide a concrete example of how the NUT parameter changes the branch structure of higher-derivative gravity black holes relative to both the static spherical case and the neutral NUT case. The paper is commendably explicit about the field equations, the horizon expansions, and the Wald entropy calculation, and the plotted metric functions are smooth and appear asymptotically flat. However, the central existence and branch-count claims are entirely numerical and currently rest on an underspecified computation: the Weyl coupling alpha is never given, and the shooting method's boundary conditions, outer radius, and tolerances are not stated. The three-branch claim is therefore not yet established as a property of a definite theory, and the apparent disconnection of the charged branches is not verifiable from the manuscript as it stands.
major comments (4)
- [Sec. 3 with Eq. (2.1), (2.9), (3.1)] The value of the Weyl coupling alpha is never stated, although alpha appears in the field equations (2.9) and in the denominators of the near-horizon expansion coefficients (2.11). The initial data used for the shooting are therefore not reproducible, and the massive spin-2 mode, whose asymptotic falloff is exp(-r/sqrt(2 alpha)), depends on alpha. Please state the value of alpha used in all integrations and provide at least one table of the shooting parameter delta for representative parameters (e.g., n=0.5, a1=1.0, r0=0.95) so that the three-branch claim is attached to a definite theory.
- [Sec. 3, shooting boundary conditions] The asymptotic selection of delta is described only as 'adjust delta to satisfy the boundary conditions at infinity' (paragraph before Figure 1). The paper does not state the explicit conditions imposed on h, f, and a at the outer boundary, the outer radius rmax, or the integration tolerance and precision. The plotted integrations extend only to roughly r=30-35; if alpha is not small, the Yukawa tail exp(-r/sqrt(2 alpha)) may not have converged at that radius. Please specify the asymptotic boundary conditions, rmax, and tolerances, and provide residual or convergence tests demonstrating that the number of delta roots and their locations are stable under changes of rmax and integration accuracy.
- [Secs. 3.1, 3.2 and Figures 4-16] The claim of exactly three branches and the claimed smooth connection of the two delta<0 branches at a minimal horizon radius are inferred from plotted discrete data, but no root-finding details are given: no initial guesses, no bracketing intervals, no delta resolution, and no number of data points per branch. The connection at the minimum is a turning point in the shooting parameter and is precisely the type of feature that can be a numerical artifact if roots are missed or if the shooting procedure does not resolve a fold. Please add a delta-versus-r0 diagram or a table for representative parameters and demonstrate convergence of the branch curves under mesh refinement, changes in rmax, and changes in the asymptotic matching point.
- [Eq. (2.11)] The displayed formulas for h2, f2, and a2 are typeset ambiguously, with unclear denominator grouping: for example, the h2 expression combines a fraction, a subtracted term, and a separate term without showing whether the latter two share the same denominator as the first. As written, these equations cannot be used to reproduce the near-horizon initial data even in principle. Please rewrite them with unambiguous parentheses and common denominators, or provide the recurrence relations in a supplementary file or code.
minor comments (5)
- [Eq. (2.7)] In the definition F_mu nu = grad_mu A_nu - grad_nu A_nu, the second term should read grad_nu A_mu; this is presumably a typographical error.
- [Sec. 2, below Eq. (2.10)] The phrase 'can be impressed in terms of h1, f1, a0, a1' should read 'can be expressed in terms of'; similarly, 'Talor expansion' before Eq. (2.12) should be 'Taylor expansion'.
- [Figure 1 and captions] The horizontal axis in Figure 1 is labeled r0 in the captions and in some panels, but the plotted variable is the radial coordinate r; please correct the labels and captions.
- [Sec. 3, parameter terminology] The parameter a1 is variously called the 'electric charge parameter' and the 'coupling parameter'; since a1 is the near-horizon coefficient of the gauge potential, please clarify its relation to the physical electric charge of the solution.
- [Sec. 3, mass extraction] The asymptotic extraction of the mass through 'g_tt = -1 + 2m/r + ...' is not justified for the NUT ansatz (2.8), where g_tt = -h(r) and h(r) may approach a constant with subleading corrections; please define the mass parameter m precisely in terms of h(r) or f(r).
Circularity Check
No circularity: the shooting parameter delta is a solution label fixed by asymptotic boundary conditions, and temperature/entropy are computed from the resulting solutions rather than used as fitting targets.
full rationale
The paper's derivation chain is self-contained. The field equations (2.9) are integrated from horizon expansions (2.10)-(2.11) with the shooting parameter delta introduced in (3.1) as f1 = 1/r0 + delta. The parameter delta is adjusted to satisfy prescribed boundary conditions at infinity, not to reproduce any particular temperature or entropy value. Temperature (2.12) and Wald entropy (2.17) are evaluated only after a solution is obtained; the entropy expression S = pi(r0^2+n^2) - 4*pi*alpha*r0*delta (3.2) is a derived consequence, not an input to the shooting. The central claim of three branches for fixed n and a1 arises from the number of delta roots found at each r0, which is a direct numerical outcome of the integration, not an ansatz or a fitted target. The connection of the two negative-delta branches at a minimal horizon radius is read off from the computed temperature/entropy versus r0 curves. The neutral comparison uses the same group's prior work [19], but only as a benchmark for contrast; it does not enter the charged construction or the branch counting. Similarly, the static two-branch result [26] is cited for comparison only. While the numerical setup lacks explicit convergence tests and a stated value of alpha, those are reproducibility concerns, not circularity: no reported prediction is equivalent by construction to a fitted parameter or to a self-cited result.
Assumptions & free parameters
free parameters (5)
- Weyl coupling α =
not specified in text
- shooting parameter δ =
one positive and two negative roots per (r0, n, a1)
- NUT parameter n =
0.5, 0.7, 1.0, 1.5
- electric charge parameter a1 =
0.02, 0.1, 0.5, 1.0
- horizon radius r0 =
0.85 to 3.5 (scanned)
assumptions (4)
- domain assumption The cohomogeneity-one ansatz (2.8) with round S^2 and r-dependent h, f, a covers all relevant charged Taub-NUT black holes.
- domain assumption R=0 holds on the NUT ansatz, reducing the field equations to three second-order ODEs.
- ad hoc to paper The shooting procedure with near-horizon Taylor data and δ adjustments yields all and only the regular asymptotically flat solutions.
- standard math Wald entropy formula (2.13) with binormal (2.14) gives the correct entropy (2.17).
Cite this review
Pith. "Pith review of Charged Taub-NUT type Black Holes in Einstein-Weyl-Maxwell Theory." pith.science (2026). https://pith.science/paper/ETSNEUP4
@misc{pith2026250801790,
author = {Pith},
title = {Pith review of: Charged Taub-NUT type Black Holes in Einstein-Weyl-Maxwell Theory},
year = {2026},
howpublished = {\url{https://pith.science/paper/ETSNEUP4}},
note = {Machine review of arXiv:2508.01790}
}
read the original abstract
We numerically construct new charged Taub-NUT-type black hole solutions in four-dimensional Einstein-Weyl-Maxwell theory. We explore the effects of the NUT parameter and the electric charge parameter on the black hole solutions in great detail. Compared with static spherical black holes in the same theory, there is one major difference: there are three branches of NUT solutions, whereas there are only two branches of static spherical black hole solutions. Compared with neutral NUT solutions in Weyl gravity, where the solutions intersect with each other, the charged NUT solutions we find are disconnected.
Reference graph
Works this paper leans on
-
[1]
LIGO Scientific, Virgo collaboration, B. P. Abbott et al.,GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral, Phys. Rev. Lett.119 (2017) 161101, [1710.05832]
arXiv 2017
-
[2]
LIGO Scientific, Virgo, Fermi-GBM, INTEGRAL collaboration, B. P. Abbott et al., Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A, Astrophys. J. Lett.848 (2017) L13, [1710.05834]
arXiv 2017
-
[3]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole, Astrophys. J. Lett.875 (2019) L1, [1906.11238]
arXiv 2019
-
[4]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. II. Array and Instrumentation, Astrophys. J. Lett.875 (2019) L2, [1906.11239]
arXiv 2019
-
[5]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. III. Data Processing and Calibration, Astrophys. J. Lett.875 (2019) L3, [1906.11240]
arXiv 2019
-
[6]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole, Astrophys. J. Lett.875 (2019) L4, [1906.11241]
arXiv 2019
-
[7]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring, Astrophys. J. Lett.875 (2019) L5, [1906.11242]
arXiv 2019
-
[8]
Akiyama et al.,First M87 Event Horizon Telescope Results
Event Horizon Telescope collaboration, K. Akiyama et al.,First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole, Astrophys. J. Lett. 875 (2019) L6, [1906.11243]
arXiv 2019
Show all 30 references
-
[9]
K. S. Stelle,Renormalization of Higher Derivative Quantum Gravity, Phys. Rev. D16 (1977) 953–969
1977
-
[10]
W. Li, W. Song and A. Strominger,Chiral Gravity in Three Dimensions, JHEP 04 (2008) 082, [0801.4566]
2008 arXiv
-
[11]
Deser, R
S. Deser, R. Jackiw and S. Templeton,Topologically Massive Gauge Theories, Annals Phys. 140 (1982) 372–411
1982
-
[12]
E. A. Bergshoeff, O. Hohm and P. K. Townsend,More on Massive 3D Gravity, Phys. Rev. D 79 (2009) 124042, [0905.1259]
2009 arXiv
-
[13]
Lu and C
H. Lu and C. N. Pope,Critical Gravity in Four Dimensions, Phys. Rev. Lett.106 (2011) 181302, [1101.1971]. – 16 –
2011 arXiv
-
[14]
Deser, H
S. Deser, H. Liu, H. Lu, C. N. Pope, T. C. Sisman and B. Tekin,Critical Points of D-Dimensional Extended Gravities, Phys. Rev. D83 (2011) 061502, [1101.4009]
2011 arXiv
-
[15]
H. Lu, A. Perkins, C. N. Pope and K. S. Stelle,Black Holes in Higher-Derivative Gravity, Phys. Rev. Lett.114 (2015) 171601, [1502.01028]
2015 arXiv
-
[16]
Huang, D.-J
Y. Huang, D.-J. Liu and H. Zhang,Novel black holes in higher derivative gravity, JHEP 02 (2023) 057, [2212.13357]
2023 arXiv
-
[17]
K. Lin, A. B. Pavan, G. Flores-Hidalgo and E. Abdalla,New Electrically Charged Black Hole in Higher Derivative Gravity, Braz. J. Phys.47 (2017) 419–425, [1605.04562]
2017 arXiv
-
[18]
Wu, D.-C
C. Wu, D.-C. Zou and M. Zhang,Charged black holes in the Einstein-Maxwell-Weyl gravity, Nucl. Phys. B 952 (2020) 114942, [1904.10193]
2020 arXiv
-
[19]
Chen and H.-S
Y.-Q. Chen and H.-S. Liu,New Taub-NUT Black Holes with Massive Spin-2 Hair, 2406.03692
-
[20]
A. V. Smilga,Supersymmetric field theory with benign ghosts, Journal of Physics A: Mathematical and Theoretical47 (Jan., 2014) 052001
2014
-
[21]
H. Lü, A. Perkins, C. Pope and K. Stelle,Lichnerowicz modes and black hole families in ricci quadratic gravity, Physical Review D96 (Aug., 2017)
2017
-
[22]
J. F. Plebanski and M. Demianski,Rotating, charged, and uniformly accelerating mass in general relativity, Annals Phys. 98 (1976) 98–127
1976
-
[23]
H.-S. Liu, H. Lu and L. Ma,Thermodynamics of Taub-NUT and Plebanski solutions, JHEP 10 (2022) 174, [2208.05494]
2022 arXiv
-
[24]
R. M. Wald,The thermodynamics of black holes, Communications in Mathematical Physics 87 (1983) 259–276
1983
-
[25]
Iyer and R
V. Iyer and R. M. Wald,Some properties of noether charge and a proposal for dynamical black hole entropy, Physical Review D50 (1994) 846–864
1994
-
[26]
Li and H.-S
Z. Li and H.-S. Liu,On Charged Black Holes in Einstein-Weyl-Maxwell Theory, 2505.02340
-
[27]
R. A. Hennigar, D. Kubizňák and R. B. Mann,Thermodynamics of Lorentzian Taub-NUT spacetimes, Phys. Rev. D100 (2019) 064055, [1903.08668]
2019 arXiv
-
[28]
Wu and D
S.-Q. Wu and D. Wu,Thermodynamical hairs of the four-dimensional Taub-Newman-Unti-Tamburino spacetimes, Phys. Rev. D100 (2019) 101501, [1909.07776]
2019 arXiv
-
[29]
Wu and S.-Q
D. Wu and S.-Q. Wu,Consistent mass formulas for the four-dimensional dyonic NUT-charged spacetimes, Phys. Rev. D105 (2022) 124013, [2202.09251]
2022 arXiv
-
[30]
Awad and S
A. Awad and S. Eissa,Lorentzian Taub-NUT spacetimes: Misner string charges and the first law, Phys. Rev. D105 (2022) 124034, [2206.09124]. – 17 –
2022 arXiv
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.