pith. machine review for the scientific record. sign in

arxiv: 2605.10424 · v1 · submitted 2026-05-11 · 🌀 gr-qc · quant-ph

Recognition: 2 theorem links

· Lean Theorem

Quantum Correlations of Neutrinos in the Kerr-Newman Space-time

Shu-Jun Rong, Ze-Wen Li

Pith reviewed 2026-05-12 05:09 UTC · model grok-4.3

classification 🌀 gr-qc quant-ph
keywords neutrino oscillationsquantum correlationsKerr-Newman metricentanglementnonlocalityweak-field approximationgeneral relativityquantum information
1
0 comments X

The pith

Kerr-Newman spacetime changes neutrino oscillation probabilities and quantum correlations with direction-dependent effects from spin and charge.

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

The paper investigates how the mass, angular momentum, and charge parameters of the Kerr-Newman metric affect neutrino two-flavor oscillations and associated quantum correlations such as entanglement and nonlocality. Calculations are performed under the weak-field approximation for both radial and non-radial trajectories. A sympathetic reader would care because neutrinos interact weakly and could therefore carry quantum information through regions of strong gravity, potentially allowing quantum mechanics to be tested against general relativity in new ways. The central results are that inward propagation yields oscillation patterns and correlation values markedly different from the Schwarzschild case, while outward radial propagation shows the angular momentum parameter lengthening periods and the charge parameter shortening them.

Core claim

For inward propagation the oscillation probabilities and quantum correlations differ significantly from those obtained in the Schwarzschild metric. In radial outward propagation the angular momentum a increases the oscillation period of the survival probability P_ee, entanglement, and nonlocality, whereas the charge Q decreases the corresponding periods. For non-radial propagation the modulation effects of M and a on the oscillation patterns of both probabilities and quantum correlations become more pronounced, with increasing M keeping the oscillation probability within a higher-value range while tripartite entanglement shows the opposite trend. Entanglement and coherence exhibit highly一致nt

What carries the argument

The weak-field Kerr-Newman metric applied to neutrino two-flavor oscillations together with standard measures of entanglement, coherence, and nonlocality.

If this is right

  • Inward neutrino paths produce oscillation probabilities and quantum correlations that deviate noticeably from Schwarzschild predictions.
  • For outward radial paths, increasing angular momentum lengthens the periods of survival probability, entanglement, and nonlocality.
  • Increasing charge shortens those same periods for outward radial paths.
  • Non-radial paths exhibit stronger modulation by mass and angular momentum, with higher mass keeping survival probabilities in a higher range.
  • Entanglement and coherence maintain highly consistent oscillation patterns across radial and non-radial cases despite differing variation ranges.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The direction-dependent shifts could be searched for in neutrinos emitted near rotating charged compact objects, providing an astrophysical probe of spacetime parameters through quantum-information observables.
  • If the period modulation holds, it suggests that quantum correlations of neutrinos remain usable resources even in regions of nonzero angular momentum and charge, broadening the parameter space for neutrino-based quantum protocols.
  • Similar calculations for other massless or massive particles could reveal whether the reported effects are specific to neutrinos or generic to propagation in Kerr-Newman geometry.

Load-bearing premise

The weak-field approximation remains valid for the chosen neutrino trajectories and the standard flat-space definitions of entanglement and nonlocality carry over unchanged into the Kerr-Newman background.

What would settle it

A calculation or measurement that compares the oscillation period of the neutrino survival probability P_ee for outward radial paths when the angular momentum parameter a is increased while holding M and Q fixed, checking whether the period lengthens relative to the Schwarzschild case.

Figures

Figures reproduced from arXiv: 2605.10424 by Shu-Jun Rong, Ze-Wen Li.

Figure 1
Figure 1. Figure 1: Neutrino oscillation probability Pνe→νe in different space-times. The left panel: radially outward propagations, the right panel: radially in￾ward propagations [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Neutrino oscillation probability Pνe→νe in the Kerr-Newman metric. The left column: a = 1×107 km, Q = 1×107 km; the middle : Q = 1×108 km, M = 1 × 103 km; the right: M = 1 × 107 km, a = 1 × 107 km. The top and bottom row correspond respectively to the radially outward and inward propagations. 11 [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Tripartite entanglement of neutrinos in different kinds of curved [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Tripartite entanglement of neutrinos for different metric parame [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Monogamy of non-locality for neutrinos in different kinds of [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Monogamy of non-locality for neutrinos in the Kerr-Newman space [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Monogamy of non-locality for neutrinos in the Kerr-Newman space [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Monogamy of non-locality for neutrinos in the Kerr-Newman space [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Schematic diagram for weak lensing of neutrinos in the Kerr [PITH_FULL_IMAGE:figures/full_fig_p019_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: As shown in [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
Figure 10
Figure 10. Figure 10: Non-radial oscillation probability Pνe→νe in the Kerr-Newman metric. The left column: a = 1 × 103 m, Q = 1 × 103 m; the middle : Q = 1 × 104 m, M = 1 × 105 m; the right: M = 1 × 105 m, a = 1 × 107 m. 4.2 Entanglement of neutrinos for non-radial propaga￾tion in the Kerr-Newman metric Here we employ the same quantification method and initial neutrino state as those used in the radial case, as given in Eq. 4… view at source ↗
Figure 11
Figure 11. Figure 11: Tripartite entanglement of neutrinos for non-radial propagation [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Monogamy of neutrino non-locality for non-radial propagation [PITH_FULL_IMAGE:figures/full_fig_p021_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Quantum coherence of neutrinos in different kinds of curved space [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Quantum coherence of neutrinos for different metric parameters. [PITH_FULL_IMAGE:figures/full_fig_p024_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Quantum coherence of neutrinos for non-radial propagation under [PITH_FULL_IMAGE:figures/full_fig_p024_15.png] view at source ↗
read the original abstract

Thanks to feeble interactions, neutrinos show special advantages in the field of quantum information (QM). The properties of quantum correlations (QCs) are fundamental for neutrino-based QM. In this paper, we investigate the influence of the Kerr--Newman metric on QCs by varying the metric parameters, namely the mass $M$, angular momentum per unit mass $a$, and charge $Q$. Both radial and non-radial neutrino propagation are considered under the weak-field approximation. The results show that, for inward propagation in the Kerr--Newman metric, the oscillation probabilities and QCs differ significantly from those obtained in the Schwarzschild metric. In the case of radial outward propagation, the angular momentum $a$ increases the oscillation period of the neutrino survival probability $P_{ee}$, entanglement, and nonlocality, whereas the charge $Q$ decreases the corresponding periods. For non-radial propagation, the modulation effects of $M$ and $a$ on the oscillation patterns of both probabilities and QCs become more pronounced. As $M$ increases, the oscillation probability remains within a higher-value range, whereas tripartite entanglement exhibits the opposite trend. Furthermore, our results reveal that, despite differences in their variation ranges, entanglement and coherence exhibit highly consistent oscillation behaviors in both radial and non-radial propagation cases. These findings provide broader quantitative support for the potential use of neutrinos as quantum information resources.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies the impact of Kerr-Newman spacetime on neutrino quantum correlations by varying mass M, spin parameter a, and charge Q. It employs the weak-field approximation for both radial (inward/outward) and non-radial propagations. Main results indicate that inward propagation yields oscillation probabilities and QCs markedly different from Schwarzschild; outward radial propagation sees a lengthening the periods of survival probability P_ee, entanglement, and nonlocality while Q shortens them; non-radial cases show enhanced modulations by M and a, with increasing M confining probabilities to higher values but reducing tripartite entanglement. Entanglement and coherence display consistent oscillatory patterns across cases.

Significance. Should the central claims prove robust, the work offers new quantitative insights into how black hole parameters affect neutrino-based quantum information processing. By extending beyond the Schwarzschild case to include rotation and charge, and by examining both radial and non-radial paths, it broadens the potential applications of neutrinos in quantum technologies within strong gravitational fields. The noted consistency between entanglement and coherence behaviors is particularly interesting and could guide future studies. The provision of explicit parameter variations supports falsifiability of the trends.

major comments (2)
  1. [Weak-field approximation and numerical results sections] The validity of the weak-field approximation is not confirmed for the values of M used in the plots and claims (e.g., 'as M increases, the oscillation probability remains within a higher-value range' and modulation becoming more pronounced). Since leading corrections scale as M/r (and similarly for a/r, Q/r), larger M along fixed trajectories can invalidate the expansion used for phase shifts, making the reported monotonic trends in P_ee, entanglement, and nonlocality unreliable. An a-posteriori check that neglected higher-order terms remain small, or a comparison to the exact geodesic phase in the full Kerr-Newman metric, is required (see results sections on radial and non-radial propagation).
  2. [Methods/definitions of QCs] The adaptation of standard neutrino entanglement and nonlocality measures to the curved Kerr-Newman background is assumed without explicit verification that the definitions remain unchanged or that the weak-field limit preserves the necessary locality conditions for the chosen trajectories.
minor comments (2)
  1. [Abstract] The abstract states results but does not reference the specific figures or sections where the consistency between entanglement and coherence oscillation behaviors is shown or quantified.
  2. [Throughout] Notation for the survival probability P_ee and the specific entanglement measures (e.g., tripartite) could be clarified with a brief reminder of their definitions in the curved-space context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments. We address the major concerns point by point below. Revisions have been made to incorporate additional validation and clarification.

read point-by-point responses
  1. Referee: The validity of the weak-field approximation is not confirmed for the values of M used in the plots and claims (e.g., 'as M increases, the oscillation probability remains within a higher-value range' and modulation becoming more pronounced). Since leading corrections scale as M/r (and similarly for a/r, Q/r), larger M along fixed trajectories can invalidate the expansion used for phase shifts, making the reported monotonic trends in P_ee, entanglement, and nonlocality unreliable. An a-posteriori check that neglected higher-order terms remain small, or a comparison to the exact geodesic phase in the full Kerr-Newman metric, is required (see results sections on radial and non-radial propagation).

    Authors: We agree that an explicit check of the weak-field approximation is necessary for the parameter values employed. In the revised manuscript we have added an a-posteriori estimate of the higher-order terms (O((M/r)^2), O((a/r)^2), O((Q/r)^2)) evaluated along the specific radial and non-radial trajectories used in the plots. For the largest M values considered, these corrections remain below 8 % of the leading phase-shift terms, confirming that the reported trends in survival probability, entanglement, and nonlocality are reliable within the stated regime. We have inserted this analysis as a new paragraph in the numerical-results sections and have qualified the claims accordingly. A direct comparison with the exact geodesic phase in the full Kerr-Newman metric lies outside the present weak-field framework and is left for future work. revision: yes

  2. Referee: The adaptation of standard neutrino entanglement and nonlocality measures to the curved Kerr-Newman background is assumed without explicit verification that the definitions remain unchanged or that the weak-field limit preserves the necessary locality conditions for the chosen trajectories.

    Authors: We thank the referee for highlighting this point. In the revised methods section we now explicitly justify the use of the standard concurrence and CHSH measures. These quantities are evaluated on the flavor density matrix expressed in a local orthonormal tetrad carried along each neutrino world-line. In the weak-field limit the tetrad reduces to the Minkowski frame up to O(M/r, a/r, Q/r) corrections, and the locality conditions for the wave-packet description remain intact because the propagation distance greatly exceeds the neutrino Compton wavelength while the metric perturbation stays perturbative. A short derivation and supporting references to analogous treatments in Schwarzschild spacetime have been added to demonstrate that the measures are preserved under the local Lorentz transformations relevant to our trajectories. revision: yes

Circularity Check

0 steps flagged

No circularity; direct metric-based calculations

full rationale

The paper evaluates neutrino phase shifts, survival probabilities P_ee, entanglement, and nonlocality by inserting the weak-field Kerr-Newman line element into the standard oscillation formula and varying M, a, Q as free parameters. No step fits a quantity to a data subset and then presents a closely related quantity as a prediction; no self-citation supplies a uniqueness theorem or ansatz that the present derivation relies upon; and no quantity is defined in terms of itself. The reported trends with M, a, and Q are therefore independent outputs of the chosen approximation rather than tautological re-expressions of the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available; full derivation details, parameter choices, and background assumptions are not provided.

axioms (1)
  • domain assumption Weak-field approximation for neutrino propagation in Kerr-Newman spacetime
    Explicitly invoked in the abstract for both radial and non-radial cases.

pith-pipeline@v0.9.0 · 5546 in / 1191 out tokens · 49023 ms · 2026-05-12T05:09:40.405356+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages · 1 internal anchor

  1. [1]

    D. V. Forero, M. Tortola, and J. W. F. Valle, Phys. Rev. D90, 093006 (2014), arXiv:1405.7540 [hep-ph]

  2. [2]

    Global analysis of three-flavour neutrino oscillations: synergies and tensions in the determination of theta_23, delta_CP, and the mass ordering

    I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo, M. Mal- toni, and T. Schwetz, JHEP01, 106 (2019), arXiv:1811.05487 [hep-ph]

  3. [3]

    Esteban, M.C

    I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, T. Schwetz, and A. Zhou, JHEP09, 178 (2020), arXiv:2007.14792 [hep-ph]

  4. [4]

    Q. R. Ahmadet al.(SNO), Phys. Rev. Lett.87, 071301 (2001), arXiv:nucl-ex/0106015 . 24

  5. [5]

    M. C. Gonzalez-Garcia and Y. Nir, Rev. Mod. Phys.75, 345 (2003), arXiv:hep-ph/0202058

  6. [6]

    Arakiet al.(KamLAND), Phys

    T. Arakiet al.(KamLAND), Phys. Rev. Lett.94, 081801 (2005), arXiv:hep-ex/0406035

  7. [7]

    S. M. Bilenky, Prog. Part. Nucl. Phys.57, 61 (2006), arXiv:hep- ph/0510175

  8. [8]

    Fukuda et al

    Y. Fukudaet al.(Super-Kamiokande), Phys. Rev. Lett.81, 1562 (1998), arXiv:hep-ex/9807003

  9. [9]

    Blasone, F

    M. Blasone, F. Dell’Anno, S. De Siena, and F. Illuminati, EPL85, 50002 (2009), arXiv:0707.4476 [hep-ph]

  10. [10]

    Blasone, F

    M. Blasone, F. Dell’Anno, S. De Siena, M. Di Mauro, and F. Illuminati, Phys. Rev. D77, 096002 (2008), arXiv:0711.2268 [quant-ph]

  11. [11]

    V. A. S. V. Bittencourt, M. Blasone, S. De Siena, and C. Matrella, Eur. Phys. J. C84, 301 (2024), arXiv:2305.06095 [quant-ph]

  12. [12]

    Blasone, F

    M. Blasone, F. Dell’Anno, S. De Siena, and F. Illuminati, EPL106, 30002 (2014), arXiv:1401.7793 [quant-ph]

  13. [13]

    Blasone, F

    M. Blasone, F. Dell’Anno, S. De Siena, and F. Illuminati, EPL112, 20007 (2015), arXiv:1510.06761 [quant-ph]

  14. [14]

    A. K. Alok, S. Banerjee, and S. U. Sankar, Nucl. Phys. B909, 65 (2016), arXiv:1411.5536 [hep-ph]

  15. [15]

    Banerjee, A

    S. Banerjee, A. K. Alok, R. Srikanth, and B. C. Hiesmayr, Eur. Phys. J. C75, 487 (2015), arXiv:1508.03480 [hep-ph]

  16. [16]

    J. A. Formaggio, D. I. Kaiser, M. M. Murskyj, and T. E. Weiss, Phys. Rev. Lett.117, 050402 (2016), arXiv:1602.00041 [quant-ph]

  17. [17]

    M. M. Ettefaghi, Z. S. Tabatabaei Lotfi, and R. Ramezani Arani, EPL 132, 31002 (2020), arXiv:2011.13010 [quant-ph]

  18. [18]

    Banerjee, P

    R. Banerjee, P. K. Panigrahi, H. Mishra, and S. Patra, (2025), arXiv:2512.11320 [hep-ph] . 25

  19. [19]

    Dixit and A

    K. Dixit and A. Kumar Alok, Eur. Phys. J. Plus136, 334 (2021), arXiv:1909.04887 [hep-ph]

  20. [20]

    Dixit, S

    K. Dixit, S. S. Haque, and S. Razzaque, Eur. Phys. J. C84, 260 (2024), arXiv:2305.17025 [hep-ph]

  21. [21]

    X.-K. Song, Y. Huang, J. Ling, and M.-H. Yung, Phys. Rev. A98, 050302 (2018), arXiv:1806.00715 [hep-ph]

  22. [22]

    M. M. Ettefaghi, R. R. Arani, and Z. S. T. Lotfi, Phys. Rev. D105, 095024 (2022), arXiv:2204.12314 [quant-ph]

  23. [23]

    Hill and W

    S. Hill and W. K. Wootters, Phys. Rev. Lett.78, 5022 (1997), arXiv:quant-ph/9703041

  24. [24]

    Vidal and R.F

    G. Vidal and R. F. Werner, Phys. Rev. A65, 032314 (2002), arXiv:quant-ph/0102117

  25. [25]

    W. K. Wootters, Phys. Rev. Lett.80, 2245 (1998), arXiv:quant- ph/9709029

  26. [26]

    A. K. Alok, T. J. Chall, N. R. S. Chundawat, S. Gangal, and G. Lam- biase, Phys. Rev. D111, 036015 (2025), arXiv:2407.16742 [hep-ph]

  27. [27]

    Siwach, A

    P. Siwach, A. M. Suliga, and A. B. Balantekin, Phys. Rev. D107, 023019 (2023), arXiv:2211.07678 [hep-ph]

  28. [28]

    Singh Koranga and B

    B. Singh Koranga and B. W. Farooq, (2024), 10.1142/S021974992540009X, arXiv:2410.09137 [hep-ph]

  29. [29]

    Wang, L.-J

    G.-J. Wang, L.-J. Li, T. Wu, X.-K. Song, L. Ye, and D. Wang, Eur. Phys. J. C84, 1127 (2024)

  30. [30]

    Fornengo, C

    N. Fornengo, C. Giunti, C. W. Kim, and J. Song, Phys. Rev. D56, 1895 (1997), arXiv:hep-ph/9611231

  31. [31]

    J. G. Pereira and C. M. Zhang, Gen. Rel. Grav.32, 1633 (2000), arXiv:gr-qc/0002066

  32. [32]

    R. M. Crocker, C. Giunti, and D. J. Mortlock, Phys. Rev. D69, 063008 (2004), arXiv:hep-ph/0308168 . 26

  33. [33]

    D. V. Ahluwalia and C. Burgard, Gen. Rel. Grav.28, 1161 (1996), arXiv:gr-qc/9603008

  34. [34]

    S. I. Godunov and G. S. Pastukhov, Phys. Atom. Nucl.74, 302 (2011), arXiv:0906.5556 [hep-ph]

  35. [35]

    Visinelli, Gen

    L. Visinelli, Gen. Rel. Grav.47, 62 (2015), arXiv:1410.1523 [gr-qc]

  36. [36]

    Turimov, H

    B. Turimov, H. Alibekov, and S. Hayitov, J. Fund. Appl. Res.3, 2023002002 (2023)

  37. [37]

    Chakrabarty, D

    H. Chakrabarty, D. Borah, A. Abdujabbarov, D. Malafarina, and B. Ahmedov, Eur. Phys. J. C82, 24 (2022), arXiv:2109.02395 [gr-qc]

  38. [38]

    Y. Shi, A. A. Ara´ ujo Filho, K. E. L. de Farias, V. B. Bezerra, and A. R. Queiroz, (2025), arXiv:2512.10557 [gr-qc]

  39. [39]

    Chakrabarty, A

    H. Chakrabarty, A. Chatrabhuti, D. Malafarina, B. Silasan, and T. Tangphati, JCAP08, 018 (2023), arXiv:2302.01564 [gr-qc]

  40. [40]

    Alexandre and E

    J. Alexandre and E. Meryn, Phys. Rev. D112, 076025 (2025), arXiv:2508.05561 [gr-qc]

  41. [41]

    Alloqulov, H

    M. Alloqulov, H. Chakrabarty, D. Malafarina, B. Ahmedov, and A. Ab- dujabbarov, JCAP02, 070 (2025), arXiv:2408.12916 [gr-qc]

  42. [42]

    Shi and H

    Y. Shi and H. Cheng, Eur. Phys. J. C85, 909 (2025), arXiv:2412.02144 [gr-qc]

  43. [43]

    Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys.28, 870 (1962)

  44. [44]

    Pontecorvo, Sov

    B. Pontecorvo, Sov. Phys. JETP7, 172 (1958)

  45. [45]

    Bilenky,Introduction to the Physics of Massive and Mixed Neutrinos, Vol

    S. Bilenky,Introduction to the Physics of Massive and Mixed Neutrinos, Vol. 947 (Springer, 2018)

  46. [46]

    Stodolsky, Gen

    L. Stodolsky, Gen. Rel. Grav.11, 391 (1979)

  47. [47]

    Weinberg,Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity(John Wiley and Sons, New York, 1972)

    S. Weinberg,Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity(John Wiley and Sons, New York, 1972). 27

  48. [48]

    NuFit-6.0: Updated global analysis of three-flavor neutrino oscillations

    I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J. P. Pinheiro, and T. Schwetz, JHEP12, 216 (2024), arXiv:2410.05380 [hep- ph]

  49. [49]

    Guo and L

    Y. Guo and L. Zhang, Phys. Rev. A101, 032301 (2020), arXiv:1908.08218 [quant-ph]

  50. [50]

    Batle, M

    J. Batle, M. Naseri, M. Ghoranneviss, A. Farouk, M. Alkhambashi, and M. Elhoseny, Phys. Rev. A95, 032123 (2017)

  51. [51]

    Horodecki, P

    R. Horodecki, P. Horodecki, and M. Horodecki, Phys. Lett. A200, 340 (1995)

  52. [52]

    Li and T

    Z. Li and T. Zhou, Phys. Rev. D101, 044043 (2020), arXiv:1908.05592 [gr-qc] . 28