Pith. sign in

REVIEW 3 major objections 4 minor 26 references

Neutrino spin oscillations near a black hole

T0 review · 3 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A high-precision simulation finds that a purely toroidal magnetic field around a black hole flips neutrino spins, contradicting earlier null results.

desk verdict Claims to overturn the earlier no-spin-flip result for toroidal fields, but the stated Hamiltonian coefficients are 76–86 orders of magnitude smaller than the physical phases, so the central result is unverified and likely mis-normalized. read the letter →

arxiv 2502.05238 v2 pith:5SBNWYCH submitted 2025-02-07 hep-ph astro-ph.HEgr-qc

Mridupawan Deka , Maxim Dvornikov This is my paper · ORCID
classification hep-phastro-ph.HEgr-qc PACS 95.30.Sf97.60.Lf14.60.Lm
keywords neutrinospinoscillationsblackholeKerrmetricaccretiondisktoroidalmagneticfieldmomentgravitationalscatteringspin-flip
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that neutrinos gravitationally scattered by a rotating black hole can undergo sizable spin-flip—conversion from left-handed active to right-handed sterile neutrinos—even when the accretion disk's magnetic field is purely toroidal. Earlier studies concluded that a poloidal field component was required for the effect; that conclusion, the authors say, came from a numerical precision failure in the old code, which rounded the tiny dimensionless coefficients $V_m \sim 10^{-87}$ and $V_B \sim 10^{-76}$ down to zero. Using a new high-precision code and millions of test neutrinos per black-hole spin, they find regions of the observer's sky where $P_{LL}<1$ for both a slow and a fast rotating black hole. If correct, the result changes how observed neutrino fluxes from black-hole accretion disks should be interpreted, since spin-flip removes neutrinos from detection and could be used to probe magnetic fields near black holes.

What carries the argument

The machinery is the spin-evolution equation $i\,d\psi/dx = \hat H_x\psi$, with $\hat H_x = -U_2(\boldsymbol\sigma\cdot\boldsymbol\Omega_x)U_2^\dagger$, where $\boldsymbol\Omega = \boldsymbol\Omega_g + \boldsymbol\Omega_{\rm em} + \boldsymbol\Omega_{\rm matt}$ collects the gravitational, electromagnetic (toroidal magnetic field), and electroweak matter contributions to the precession; the polarization vector obeys $d\boldsymbol\zeta/dt = 2(\boldsymbol\zeta\times\boldsymbol\Omega)$. Trajectories are exact Kerr geodesics parameterized by elliptic integrals, and each neutrino's spin is evolved along its own geodesic. The decisive numerical feature is that the new C++ code respects the smallness of the dimensionless coefficients $V_m$ and $V_B$ that the earlier MATLAB code rounded to zero.

What would settle it

Run an independent high-precision recomputation (e.g., arbitrary-precision arithmetic or a different integration scheme) of the same scattering setup—$M=10^8M_\odot$, $a=0.02M$ and $0.98M$, $B_{\max}=320$ G, $n_e=10^{18}$ cm$^{-3}$, $\mu=10^{-13}\mu_B$—for the toroidal-only field; if $P_{LL}=1$ everywhere (no spin-flip), the central claim is refuted. A simpler decisive check is rerunning the old MATLAB code with variable-precision arithmetic to see whether the earlier null result persists.

Watch

Extended reading notes

Core claim

The central claim is that a sizable neutrino spin-flip takes place near a Kerr black hole even in the presence of a toroidal magnetic field alone, contrary to Refs. [13–17], where such spin-flip was attributed to a poloidal component. The authors simulate 2.6 million (for $a=0.02M$) and 3.3 million (for $a=0.98M$) incoming left-handed neutrinos on exact Kerr geodesics, evolving each spin through the effective Schrödinger equation with gravitational, electromagnetic, and electroweak contributions. They find contour regions in the $(\theta_{\rm obs},\phi_{\rm obs})$ plane where the probability $P_{LL}$ of remaining left-handed is noticeably less than unity, for both black-hole spins. They attribute the previous null result to a precision issue in the older MATLAB-based code: the dimensionless coefficients $V_m\sim10^{-87}$ and $V_B\sim10^{-76}$ were below MATLAB's default precision, so the plasma density and magnetic pressure terms were evaluated as zeros.

Load-bearing premise

As the paper itself states in Sec. III, the contrast with earlier work rests on a precision issue in the old MATLAB code; the load-bearing premise is therefore that the new C++ code's handling of $V_m\sim10^{-87}$ and $V_B\sim10^{-76}$ is correct and that the computed spin-flip is not itself a numerical artifact, a premise the paper does not independently validate.

Editorial extensions

If this is right

  • Observed fluxes of neutrinos scattered off a black-hole accretion disk are depleted by the factor $P_{LL}<1$; future telescopes must account for this spin-flip loss when inferring source luminosities.
  • A poloidal field is not required for astrophysically relevant spin-flip; toroidal fields, which are naturally produced in accretion disks, suffice.
  • The effect persists for both slow ($a=0.02M$) and near-extremal ($a=0.98M$) black-hole spins, so it is robust to spin magnitude within the simulated cases.
  • The spatial pattern of $P_{LL}$ over the observer's sky encodes the geometry of the disk's magnetic field, supporting neutrino tomography of black-hole magnetic fields.
  • Precision handling of exponentially small coefficients is mandatory in these simulations; codes that silently underflow these terms will report null spin-flip.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • My inference: the same precision mechanism predicts that the earlier null results should be reproduced as nonzero once recomputed with arbitrary precision, so a direct side-by-side rerun would settle the discrepancy without new astrophysics.
  • My inference: if the spin-flip is real, the magnitude of the predicted flux depletion depends sensitively on the assumed neutrino magnetic moment $\mu=10^{-13}\mu_B$; varying $\mu$ and comparing with observed event rates could turn this calculation into a constraint on $\mu$.
  • My inference: the trajectory-based spin evolution could be extended to include flavor oscillations or multiple neutrino generations, which would produce energy-dependent modifications of $P_{LL}$ rather than a simple overall suppression.
  • My inference: because the result rests on numerical integration of tiny coefficients, an independent semi-analytic estimate of the spin-flip probability (e.g., in a simplified field model) would both validate the code and identify which parts of the sky are most promising for observation.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies spin oscillations of ultra-relativistic neutrinos gravitationally scattered by a Kerr black hole surrounded by a thick magnetized accretion disk. The authors consider only a toroidal magnetic field, in contrast to earlier work in Refs. [13–17] that included a poloidal component, and claim that a sizable spin-flip (P_LL < 1) occurs. The calculation uses exact Kerr geodesics for millions of test neutrinos and integrates a spin-evolution Hamiltonian along each trajectory with a newly written C++ code. The authors attribute the absence of spin-flip in the earlier MATLAB-based studies to the smallness of two dimensionless coefficients, V_m ~ 1e-87 and V_B ~ 1e-76, which they say produced zero plasma density and magnetic pressure in the old code. The results are presented as color contour maps of P_LL over the observer angle plane for two black hole spins.

Significance. If the claimed result is correct, it overturns the previous conclusion that a poloidal magnetic-field component is required for sizable neutrino spin-flip near a black hole, and it would have implications for interpreting observations of astrophysical neutrinos and for neutrino tomography of black-hole magnetospheres. The paper's strengths include the use of exact geodesic motion in Kerr spacetime, a high-statistics parallel simulation with millions of test neutrinos, and a clear geometric argument (Fig. 2) that a toroidal field can be transverse to the neutrino momentum. However, the numerical code is not released, no quantitative P_LL values or error estimates are given, and no validation against an analytic or previously benchmarked case is provided. The central claim therefore currently rests on an unverifiable simulation, and the manuscript's own explanation of the discrepancy in terms of the tiny coefficients V_m and V_B is, as written, internally difficult to assess.

major comments (3)
  1. [Sec. III, Eq. (11)] The role of the coefficients V_m and V_B is not defined, and the explanation of the old MATLAB failure is not self-consistent as stated. The text reports V_m = G_F/(sqrt(2) m_p r_g^3) ~ 1e-87 and V_B = (mu/r_g)^2 ~ 1e-76 and says these were too small for MATLAB's default precision, causing the old code to yield zero plasma density and magnetic pressure. However, with the paper's own inputs (n_e(max)=1e18 cm^-3, B_max=320 G, mu=1e-13 mu_B, M=1e8 M_sun), the physical spin-precession phases are sqrt(2) G_F n_e r_g ~ 0.19 and mu B r_g ~ 0.28, i.e., of order unity and 76-87 orders of magnitude larger than the quoted values. If V_m and V_B are the prefactors multiplying dimensionless density and magnetic-field profile functions in H_x, then the matter and magnetic terms are utterly negligible and the claimed sizable P_LL<1 cannot follow from these interactions. If they are not the prefactors, then the stated reason for the old null result is incomplete or incorrect. The manuscript never displays the explicit matrix elements of H_x or the definitions of the dimensionless density and field variables, so the reader cannot determine which case holds. The authors must present H_x explicitly, define the dimensionless variables, and show how the physical O(1) phases emerge from the stated inputs.
  2. [Sec. IV, Fig. 1] The central quantitative claim is not supported by any numbers. Figures 1(a) and 1(b) are color contour projections of P_LL, but the paper gives no colorbar, no numeric P_LL values, and no statistical uncertainties. The statement that a 'sizable' spin-flip occurs is therefore unquantified. Because the paper's thesis is that a previous numerical code produced a qualitatively wrong null result, the new code must be validated: at least one analytic limit (e.g., flat spacetime with constant magnetic field and matter density, where the Rabi-like spin-flip probability is known), a convergence test for the Adams-Bashforth-Moulton integrator, and a side-by-side comparison with the old MATLAB code for identical inputs are needed. Without such validation, the new spin-flip signal could itself be a numerical artifact.
  3. [Sec. III, Sec. II B] Several inputs required for reproducibility are missing. The neutrino energy E never appears in the paper, although the matter and magnetic terms in the spin Hamiltonian generally depend on E (e.g., through the 1/E factor in the quasi-classical Hamiltonian). The density profile n_e(r,theta) and the magnetic-field profile B(r,theta) of the Polish doughnut are not given; only maximal values are quoted. The sampling of initial conditions (impact parameters or the distribution of L and Q) is not described. Without these details, the simulation cannot be reproduced, and the dependence of the claimed spin-flip on the neutrino energy and on the accretion-disk model cannot be assessed.
minor comments (4)
  1. [Eq. (9)] In the sentence preceding Eq. (9), 'outgoing neutrions' should be 'outgoing neutrinos'.
  2. [Sec. III] The predictor-corrector methods are called 'Adam–Bashforth' and 'Adam–Moulton'; the correct spellings are 'Adams–Bashforth' and 'Adams–Moulton'.
  3. [Table I] The entry 'T otal' contains a spacing artifact and should read 'Total'.
  4. [Sec. V] The phrase 'there are non-zero probabilities that neutrino spin-flip happens' should be rephrased, e.g., 'there is a nonzero probability of neutrino spin-flip'.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the spin-flip result is a simulation output, not a fitted parameter or a self-referential construction.

full rationale

The paper's derivation chain is: Kerr geodesic trajectories for test neutrinos, a spin-precession Hamiltonian taken from the authors' earlier work, and a numerical integration of the resulting Schrödinger equation for millions of incoming neutrinos. The central output, P_LL, is computed directly from these inputs; no parameter is fitted to the target result. The Hamiltonian in Eq. (11) is attributed to Refs. [15,16], which are the authors' own prior papers, but this is a normal physical input rather than a circular construction: the cited works derive the precession terms from standard electroweak and electromagnetic interactions, and the present paper does not define the spin-flip probability in terms of the Hamiltonian coefficients by fiat. The parameters M = 10^8 M_sun, a = 2e-2 M or 0.98 M, n_e(max) = 10^18 cm^-3, B_max = 320 G, and mu = 10^-13 mu_B are fixed external values from the literature, not adjusted to produce the claimed spin-flip. The contrast with Refs. [13-17] is attributed to a claimed precision issue in the old MATLAB code; this is a numerical-reproducibility claim, not a reduction of the prediction to its inputs. The quoted coefficients V_m ~ 1e-87 and V_B ~ 1e-76 raise a serious correctness concern about normalization, since the physical spin-precession phases estimated from the paper's own inputs are O(0.1-0.3), but that is a numerical or physical consistency issue rather than circularity: even if the Hamiltonian is mis-scaled, the output is still a genuine consequence of the stated equations. A score of 2 reflects the presence of self-citations in the Hamiltonian source while recognizing that the central claim has independent computational content and is not a self-fulfilling fit.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The simulation depends on a standard Kerr geodesic framework, a spin evolution Hamiltonian borrowed from the authors' prior work, physical disk parameters from the literature, and several modeling approximations. None of these are fitted to the target result, but the chosen field strength and neutrino magnetic moment set the overall scale of the effect.

free parameters (5)
  • Black hole mass M = 1e8 solar masses
    Fixed as an input parameter; determines the Schwarzschild radius and the overall scale of the gravitational field.
  • Black hole spin a = 0.02 M and 0.98 M
    Two chosen values to represent slowly and rapidly rotating black holes; the spin affects the geodesics and hence the scattering angles.
  • Maximal electron density n_e^(max) = 1e18 cm^-3
    Taken from a cited astrophysical observation; sets the electroweak interaction strength with the disk matter.
  • Maximal toroidal magnetic field B_max = 320 G (about 1% of Eddington limit)
    Chosen from an astrophysical estimate; the spin-flip probability scales with the magnetic field strength, so this value is critical for the 'sizable' claim.
  • Neutrino magnetic moment mu = 1e-13 Bohr magnetons
    Assumed below the astrophysical upper bound; sets the electromagnetic interaction strength and directly controls the spin precession rate.
assumptions (6)
  • standard math Kerr geodesic equations for ultra-relativistic test particles (Eqs. 1-9) from Ref. [20].
    The paper uses known Kerr geodesic results to compute neutrino trajectories; these are not derived in the paper.
  • domain assumption Spin evolution equation (Eq. 10) and effective Schrodinger equation (Eq. 11) with Omega from Refs. [15,16].
    The Hamiltonian for spin precession, including gravitational, electromagnetic, and electroweak contributions, is borrowed from the authors' prior work without derivation.
  • domain assumption Forward scattering approximation for the neutrino-matter electroweak interaction [25].
    Assumes the matter interaction can be treated in the forward scattering limit, as is standard in neutrino spin oscillation studies.
  • domain assumption One neutrino generation with only a diagonal magnetic moment.
    The paper restricts to spin oscillations within a single flavor and ignores off-diagonal magnetic moments and flavor mixing.
  • domain assumption The accretion disk is a 'Polish doughnut' with a purely toroidal magnetic field [18,19].
    The disk geometry and field configuration are adopted from the literature; the paper explicitly excludes poloidal fields.
  • domain assumption Neutrinos are ultra-relativistic test particles with negligible back-reaction on the spacetime.
    The paper treats neutrinos as test particles moving on Kerr geodesics, ignoring any gravitational influence of the neutrinos themselves.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Neutrino spin oscillations near a black hole." pith.science (2026). https://pith.science/paper/5SBNWYCH

@misc{pith2026250205238,
  author       = {Pith},
  title        = {Pith review of: Neutrino spin oscillations near a black hole},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5SBNWYCH}},
  note         = {Machine review of arXiv:2502.05238}
}
read the original abstract

In this work, we study neutrino spin oscillations in the case when they are gravitationally scattered off a rotating Kerr black hole surrounded by a thick magnetized accretion disk. We consider only toroidal magnetic field inside the disk. Neutrino spin precession is caused by the interaction of the neutrino magnetic moment with the magnetic field in the disk. Our treatment of the spin oscillations of the observed neutrino fluxes is based on numerical simulations of the propagation of a large number of incoming test neutrinos using High Performance Parallel Computing. We briefly discuss our results and their applications in the observations of astrophysical neutrinos.

Figures

Figures reproduced from arXiv: 2502.05238 by the authors.

Figure 1
Figure 1. FIG. 1. Contour projections of [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Schematic diagram showing neutrino trajectories which cross the accretion disk. The toroidal [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 14 canonical work pages

  1. [1]

    M. A. Aceroet al.[NOvA], Phys. Rev. D106, no.3, 032004 (2022) doi:10.1103/PhysRevD.106.032004. [arXiv:2108.08219 [hep-ex]]

    arXiv 2022

  2. [2]

    Giunti and C

    C. Giunti and C. W. Kim,Fundamentals of Neutrino Physics and Astrophysics(Oxford University Press, Norfolk, 2007), p. 245

    work page 2007

  3. [3]

    A. G. Beda, V. B. Brudanin, V. G. Egorov, D. V. Medvedev, V. S. Pogosov, M. V. Shirchenko and A. S. Starostin, Adv. High Energy Phys.2012, 350150 (2012) doi:10.1155/2012/350150

    doi:10.1155/2012/350150 2012

  4. [4]

    Dvornikov and A

    M. Dvornikov and A. Studenikin, Phys. Rev. D69, 073001 (2004) doi:10.1103/PhysRevD.69.073001 [arXiv:hep-ph/0305206 [hep-ph]]

    arXiv 2004

  5. [5]

    Broggini, C

    C. Broggini, C. Giunti and A. Studenikin, Adv. High Energy Phys.2012, 459526 (2012) doi:10.1155/2012/459526 [arXiv:1207.3980 [hep-ph]]

    arXiv 2012

  6. [6]

    Electromagnetic neutrinos in laboratory experiments and astrophysics

    C. Giunti, K. A. Kouzakov, Y.-F. Li, A. V. Lokhov, A. I. Studenikin, and S. Zhou, Ann. Phys. (Berlin) 528, 198 (2016); arXiv:1506.05387

    work page Pith review arXiv 2016

  7. [7]

    Fujikawa and R

    K. Fujikawa and R. Shrock, Phys. Rev. Lett.45, 963 (1980) doi:10.1103/PhysRevLett.45.963

    doi:10.1103/physrevlett.45.963 1980

  8. [8]

    O. L. Caballero, G. C. McLaughlin and R. Surman, Astrophys. J.745, 170 (2012) doi:10.1088/0004- 637X/745/2/170 [arXiv:1105.6371 [astro-ph.HE]]

    arXiv 2012

Show all 26 references
  1. [9]

    L. B. Okun, M. B. Voloshin and M. I. Vysotsky, Sov. Phys. JETP64, 446-452 (1986) ITEP-86-82

    1986

  2. [10]

    A. A. Pomeranski ˘ ı and I. B. Khriplovich, J. Exp. Theor. Phys.86, 839 (1998); gr-qc/9710098

    1998 arXiv

  3. [11]

    Dvornikov, Int

    M. Dvornikov, Int. J. Mod. Phys. D15, 1017-1034 (2006) doi:10.1142/S021827180600870X [arXiv:hep- ph/0601095 [hep-ph]]

    2006

  4. [12]

    Dvornikov, JCAP06, 015 (2013) doi:10.1088/1475-7516/2013/06/015 [arXiv:1306.2659 [hep-ph]]

    M. Dvornikov, JCAP06, 015 (2013) doi:10.1088/1475-7516/2013/06/015 [arXiv:1306.2659 [hep-ph]]

    2013 arXiv

  5. [13]

    Dvornikov, Class

    M. Dvornikov, Class. Quant. Grav.40, no.1, 015002 (2023) doi:10.1088/1361-6382/aca45a [arXiv:2206.00042 [hep-ph]]

    2023 arXiv

  6. [14]

    Dvornikov, Phys

    M. Dvornikov, Phys. Part. Nucl. Lett.20, no.3, 461-465 (2023) doi:10.1134/S154747712303024X [arXiv:2304.03622 [hep-ph]]

    2023 arXiv

  7. [15]

    Dvornikov, Int

    M. Dvornikov, Int. J. Mod. Phys. D33, no.14, 2340001 (2024) doi:10.1142/S0218271823400011 [arXiv:2212.03479 [hep-ph]]

    2024 arXiv

  8. [16]

    Dvornikov, JCAP09, 039 (2023) doi:10.1088/1475-7516/2023/09/039 [arXiv:2307.10126 [astro- ph.HE]]

    M. Dvornikov, JCAP09, 039 (2023) doi:10.1088/1475-7516/2023/09/039 [arXiv:2307.10126 [astro- ph.HE]]. 10

    2023 arXiv

  9. [17]

    Deka and M

    M. Deka and M. Dvornikov, Phys. Atom. Nucl.87, no.4, 483-488 (2024) doi:10.1134/S1063778824700327 [arXiv:2311.14475 [hep-ph]]

    2024 arXiv

  10. [18]

    Astronomy and Astrophysics63, 221 (1978)

    Abramowicz, M., Jaroszynski, M., Sikora, M. Astronomy and Astrophysics63, 221 (1978)

    1978

  11. [19]

    S. S. Komissarov, Mon. Not. R. Astron. Soc.368, 993 (2006); astro-ph/0601678

    2006 arXiv

  12. [20]

    S. E. Gralla, A. Lupsasca, and A. Strominger, Mon. Not. R. Astron. Soc.475, 3829 (2018); arXiv:1710.11112

    2018 arXiv

  13. [21]

    Abramowitz and I

    M. Abramowitz and I. Stegun,Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables(National Bureau of Standards, Washington D.C., 1964)

    1964

  14. [22]

    Jiang, A

    J. Jiang, A. C. Fabian, T. Dauser, L. Gallo, J. A. Garcia, E. Kara, M. L. Parker, J. A. Tomsick, D. J. Walton, and C. S. Reynolds, Mon. Not. R. Astron. Soc.489, 3436 (2019); arXiv:1908.07272

    2019 arXiv

  15. [23]

    V. S. Beskin,MHD Flows in Compact Astrophysical Objects: Accretion, Winds and Jets(Springer, Heidelberg 2010), p. 184

    2010

  16. [24]

    Viaux, M

    N. Viaux, M. Catelan, P. B. Stetson, G. G. Raffelt, J. Redondo, A. A. R. Valcarce, and A. Weiss, Astron. Astrophys.558, A12 (2013); arXiv:1308.4627

    2013 arXiv

  17. [25]

    Dvornikov and A

    M. Dvornikov and A. Studenikin, J. High Energy Phys.09, 016 (2002); hep-ph/0202113

    2002 arXiv

  18. [26]

    Abusleme et al

    A. Abusleme et al. (JUNO Collaboration), Prog. Part. Nucl. Phys.123, 103927 (2022); arXiv:2104.02565. 11

    2022

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.