{"id":"64cc0d54-c2b6-40a4-811e-1c5971bc8b93","arxiv_id":"2504.20236","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"Gravity-induced neutrino-antineutrino oscillations near a rotating primordial black hole are computed in Kerr-Schild polar coordinates, with quantum speed limit and entanglement entropy estimated.","lead":"This paper calculates neutrino-antineutrino flavor oscillations near a spinning primordial black hole by adding a gravity-induced potential to the neutrino mass matrix, then uses those oscillations to estimate the quantum speed limit and entanglement entropy. The application to Kerr-Schild polar coordinates is new, but the central QSLT equation appears to be dimensionally wrong, calling the paper's main quantitative result into question.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (28)/(32) use δ in eV while t from Eq. (31) is dimensionless, dropping the rg conversion; the resulting oscillation phases, QSLT ratio Eq. (57), and mass dependence in Fig. 10 are dimensionally inconsistent and likely overstate gravitational effects.","rationale":"The paper's central advertised result is quantitative: a gravitational potential in eV modifies neutrino-antineutrino transition probabilities and the QSLT ratio. A unit inconsistency in the phase is therefore load-bearing for the entire quantitative analysis. The reader's formal reason for rejection includes the dimensional inconsistency of Eq. (57); I agree with that symptom, but the same missing rg factor already appears in Eq. (28)/(32), so it also compromises the survival probabilities and entanglement results, not only the QSLT ratio. The WKBJ issue raised by the reader is a genuine robustness concern, but it is secondary: if WKBJ holds, the phase is still mis-scaled. A single recomputation with t=rg∫f dr would settle whether the claimed oscillations and mass dependence are physical or a unit artifact. Because this concern only strengthens the reader's rejection, the verdict should remain unchanged.","tokens_in":29754,"tokens_out":22765,"duration_ms":252100,"concrete_test":"Evaluate Eq. (32) with the coordinate-time conversion restored: t(r)=rg ∫_{r+}^{r} f(r',θ,a) dr', where rg=(G MPBH)/(ℏc) in eV^-1, for MPBH=10^16, 10^17, and 10^18 kg at a=0.998, θ=π/4, and compare with Fig. 10. If the three curves collapse onto one mass-independent curve and the sin^2 arguments remain much less than 1 over r≤100rg, the plotted oscillations and the empirical relation in Eq. (58) are artifacts of the missing rg factor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV fixes distances in units of rg=GM/c^2=1, so r in Eq. (31) is dimensionless and t=∫f dr is dimensionless. The same section converts Bd to eV through Eq. (16) with a factor ℏc/rg, so δ≈B0−|B| in Eq. (29) has units of eV. In Eq. (28) the phase δt is then not dimensionless; Eq. (32) uses δ∫f dr, omitting the required factor rg expressed in eV^-1. For MPBH=10^18 kg, rg≈3.8×10^-3 eV^-1, so the plotted phases are too large by roughly 260. This affects the survival probabilities and their mass dependence (Fig. 10), the QSLT ratio in Eq. (57) — which as printed has numerator S0√Q with energy units divided by a dimensionless integral — and the entanglement entropy curves that track PS. The central claims that gravity significantly changes transition probabilities within r∼100rg and that T_QSLT/T is suppressed for larger a therefore do not follow from the equations as written. This is more foundational than the WKBJ concern: even granting WKBJ, the phase is mis-scaled.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two-flavor neutrino-antineutrino oscillations near a spinning primordial black hole modeled by the Kerr metric in Kerr-Schild polar coordinates. It derives (or claims to derive) the four-vector gravitational potential Bd from the spin connection, embeds its temporal and spatial components into a Hermitian effective mass matrix for Majorana neutrinos, and then computes survival probabilities, the quantum speed limit time ratio, and entanglement entropy as functions of radial distance. The central claims are that the gravitational potential significantly modifies the flavor transition probabilities within roughly 100 gravitational radii, that a larger black hole spin suppresses the QSLT bound ratio, and that a smaller PBH mass produces faster dynamical evolution and more oscillatory entanglement entropy. The paper is an application of an established spin-curvature coupling formalism, with the new elements being the Kerr-Schild polar-coordinate choice and the QSLT/entanglement analysis.","tokens_in":30052,"tokens_out":11662,"duration_ms":131360,"significance":"If the results were correct, the paper would provide a concrete quantum-information signature of gravitational coupling to neutrinos near black holes, extending earlier work on the gravitational Zeeman effect. The use of PDG data for neutrino parameters and the transparent mass-matrix setup are positive features, and the QSLT and entanglement diagnostics are interesting tools for this problem. However, the quantitative content is not reliable: the oscillation phase in Eq. (32) is dimensionally inconsistent, the printed QSLT formula in Eq. (57) has the wrong algebraic structure and units, and the central numerical claims rest on these equations. Because the main physical conclusions are drawn from figures and empirical relations generated from these inconsistent formulas, the paper in its present form does not establish its claims. The topic may be worth revisiting after a careful revision that fixes the units, provides the closed-form Bd, validates the WKBJ assumption, and recomputes all figures.","major_comments":[{"comment":"The oscillation phase in Eq. (32) is dimensionally inconsistent. Section IV fixes r in units of rg = 1, so the integral in Eq. (31) is dimensionless, while delta in Eq. (29) has units of eV. For the phase delta * integral to be dimensionless in units with hbar = c = 1, the integral must be multiplied by rg expressed in eV^-1 (or delta must be converted to inverse length). For MPBH = 10^18 kg, rg is about 3.8 x 10^-3 eV^-1, so the plotted phases are too large by a factor of roughly 260. This directly affects the survival probabilities in Figs. 6 and 9-10. In particular, the mass dependence in Fig. 10 is suspect: since Bd scales as 1/rg through Eq. (16) and the coordinate time t from Eq. (31) carries a factor rg, the corrected phase becomes largely mass-independent, so the claim that lower PBH mass leads to highly oscillatory survival probability does not follow from the equations as written.","section":"Section VI, Eqs. (28)-(32), Fig. 10"},{"comment":"The QSLT ratio formula in Eq. (57) is algebraically and dimensionally wrong. Equation (2) gives T_QSLT = S0 / Delta H, and Eq. (55) defines Delta H = sqrt(Q). The correct ratio is therefore T_QSLT/T = S0 / (sqrt(Q) * T), with T proportional to the integral in Eq. (31). Equation (57) instead places sqrt(Q) in the numerator, which is the reciprocal of the correct expression and has units of eV divided by a dimensionless integral. This error propagates into the bottom panel of Fig. 11, the panels of Fig. 12, and the empirical relations (58)-(59). The claims that larger a suppresses T_QSLT/T and that smaller PBH mass speeds up evolution are therefore not supported by the printed formulas.","section":"Section VIII, Eqs. (55) and (57), Figs. 11-12"},{"comment":"The WKBJ approximation is not quantitatively justified. The condition as typeset in Eq. (30) is not a well-formed dimensionless inequality, and the statement that it is satisfied 'from Figs. 1 and 2' is not a quantitative check. The entire r-dependent treatment, including the t-to-r mapping in Eq. (31) and the survival probabilities derived from locally constant Bd, requires that the fractional change of Bd over a neutrino wavelength be small, i.e., lambda |dB/dr| << |B|. The paper does not verify this along the specific trajectories used in the figures. Since Bd varies rapidly near the PBH, this is a load-bearing gap for the central quantitative claims.","section":"Section VI, Eq. (30) and surrounding text"},{"comment":"The abstract and Section IV claim an analytical expression for the four-vector gravitational potential in Kerr-Schild polar coordinates, but no closed-form expression is ever displayed. The reader is given only the definition in Eq. (16), the tetrad choice in Eq. (15), and plots of B0 and |B|. Without the explicit algebraic form of Bd(r, theta, a), the claimed derivation cannot be independently checked, and the coordinate dependence that drives all subsequent results is not actually demonstrated.","section":"Section IV, after Eq. (15)"}],"minor_comments":[{"comment":"The notation for the survival probability alternates between Ps, PS, and P_S; please standardize it to a single symbol.","section":"Throughout"},{"comment":"The WKBJ condition contains typographical errors involving p and p c; it should be rewritten as a dimensionless inequality, preferably in the form lambda |dB/dr| << |B| or with an explicit characteristic length scale.","section":"Eq. (30)"},{"comment":"The Bures angle expression has a missing ket bracket: |psi_e(r)> should appear inside the inner product, not |psi_e(r>.","section":"Eq. (54)"},{"comment":"The panel labels for Delta H and T_QSLT/T are garbled in the typeset version; the axis labels should state the correct variable names and units.","section":"Fig. 11"},{"comment":"The data availability statement cites Ref. [59], which is the PDG review; it does not describe the data generated for the figures in this paper. Please state explicitly whether any numerical data or code is being released.","section":"Data availability"}],"recommendation":"reject","confidential_remarks":"The paper is heavily based on the authors' previous work on neutrino oscillations in Kerr spacetime, and the incremental elements are the Kerr-Schild polar coordinates and the QSLT/entanglement analysis. The dimensional inconsistency in Eqs. (32) and (57) is decisive: it invalidates the main quantitative claims and all figures that depend on them. A revision would require redoing the numerical analysis with correct units and providing the missing closed-form Bd, so I do not see a path to publication for the manuscript in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this paper has a load-bearing unit error, so the central quantitative claims do not follow from the equations as written. The one genuinely new thing—explicit use of Kerr-Schild polar coordinates for a spinning PBH—is undermined because the advertised analytical expression for B_d is never written down.\n\nWhat the paper does well: the structure is clear, the PDG neutrino parameters are handled properly, and the WKBJ approximation is at least discussed, not just assumed silently. The parametric scans over θ, a, and M_PBH are straightforward and complete.\n\nThe softness is in the middle. Section IV measures distances in units of r_g = 1, so r and t from Eq. (31) are dimensionless. But B^d is converted to eV via Eq. (16), making δ in Eq. (29) have units of eV. The phase δt in Eq. (28) is then not dimensionless; the correct phase is δ r_g ∫ f dr, with r_g in eV^-1. For M_PBH = 10^18 kg, r_g ≈ 3.8×10^-3 eV^-1, so the plotted phases are too large by roughly a factor of 260. This propagates into the survival probabilities, the entanglement entropy, and the mass-dependence in Fig. 10. Eq. (57) shows the same problem in different clothes: the numerator has units of eV, the denominator is dimensionless, so T_QSLT/T as printed is not a ratio at all. The correct expression is S0/(√Q r_g ∫ f dr). The empirical relations in Eqs. (58)-(59) therefore do not follow from the written equations.\n\nSeparately, the claimed analytical expression for the four-vector gravitational potential in KSP coordinates never appears. The paper shows plots, but the formulas for B0, B1, B2, B3 are not displayed. That makes the central new result unreproducible. I'd want those expressions in an appendix before taking the new-coordinate claim seriously.\n\nThe WKBJ issue the reader flagged is real but secondary. Eq. (30) is asserted, not demonstrated, and the paper gives no numbers to back the claim that the coupling is constant over local patches. Still, even granting WKBJ, the phase is mis-scaled.\n\nWho is this for? People working on gravity-induced neutrino-antineutrino oscillations and quantum speed limits in curved spacetime. The framework comes from the authors' earlier papers, but the KSP application is a legitimate extension. As written, I would not cite it. I would, however, send it to peer review: the error is fixable, the framework is established, and a referee can push the authors to write down the potentials and fix the units. Reject as written, but not because the idea is dead.","headline":"A load-bearing unit error in the oscillation phase and QSLT ratio guts the central quantitative claims; the KSP coordinate application is new but the advertised analytical potential is never shown.","tokens_in":30577,"tokens_out":4370,"would_cite":false,"duration_ms":44845,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.70.-s","14.60.Pq","03.65.-w","03.67.-a"],"model":"deepseek-v4-flash","headline":"A spinning black hole's gravity modifies neutrino-antineutrino flavor oscillations and the quantum speed limit of the flavor change.","keywords":["neutrino oscillations","neutrino-antineutrino oscillations","Kerr-Schild polar coordinates","four-vector gravitational potential","quantum speed limit time","entanglement entropy","primordial black holes","gravitational Zeeman effect"],"falsifier":"Compute the neutrino survival probability near a $10^{18}$ kg spinning primordial black hole by numerically integrating the Dirac equation in the Kerr background without assuming a locally constant $B^d$, for $a = 0.998$ and $\\theta = \\pi/3$. If the exact result differs materially from the paper's Eq. (32) for radii between about $3$ and $6$ gravitational radii, the WKBJ stitching assumption fails and the claimed suppression of oscillations near the hole is not established.","tokens_in":29570,"feed_emoji":"🌀","tokens_out":7194,"duration_ms":71032,"temperature":0.7,"pith_summary":"This paper tries to establish that the gravitational field of a spinning primordial black hole, described by the Kerr metric in Kerr-Schild polar coordinates, generates a four-vector gravitational potential that enters the Dirac equation as an axial term and changes how neutrinos oscillate into antineutrinos and into other flavors. That potential depends on the polar angle of the neutrino's position relative to the spin axis, the distance from the hole, and the dimensionless spin parameter a, and it vanishes for a non-spinning hole. If correct, the result means that near a spinning black hole the flavor survival probability, the quantum speed limit time ratio, and the entanglement entropy of the neutrino state are all altered by gravity within roughly 100 gravitational radii, after which standard vacuum oscillations return. The authors connect the same mechanism to a gravitational Zeeman splitting of neutrino and antineutrino energies and find that the resulting neutrino-number asymmetry at the Sun's surface is many orders of magnitude below current observations.","feed_headline":"Black hole spin speeds up neutrino flavor transitions near the horizon","feed_subtitle":"Gravity from a spinning black hole changes neutrino-antineutrino mixing, speed limits, and entanglement near the hole.","key_machinery":"The load-bearing object is the four-vector gravitational potential $B^d$ built from the spin connection, evaluated in Kerr-Schild polar coordinates. It produces an axial-vector term $B_d\\gamma^d\\gamma^5$ in the Dirac Lagrangian, which for Majorana neutrinos becomes an effective mass matrix that couples neutrino and antineutrino two-spinors. The computations then rest on three tools: the WKBJ approximation that the neutrino wavelength is small enough for $B^d$ to be locally constant, the ultrarelativistic approximation for the energy eigenvalues, and the null radial geodesic mapping from time to radial distance. The quantum speed limit ratio is built from the Bures angle $S_0 = \\cos^{-1}\\sqrt{P_S}$ and the energy variance $\\Delta H$.","core_discovery":"The central claim is that in the Kerr-Schild polar coordinate form of the Kerr metric the spin connection yields a nonzero four-vector gravitational potential $B^d = \\epsilon^{abcd}\\omega_{bac}$, and that this potential changes the effective mass matrix of a Majorana neutrino-antineutrino system. Diagonalizing that matrix gives a mixing angle and an energy splitting $\\delta \\simeq B_0 - |\\mathbf{B}|$ that drive gravity-induced oscillations between neutrino and antineutrino states. Using the WKBJ approximation and a null radial geodesic to convert propagation time to radial distance, the survival probability of an initial flavor state is computed; the same unitary evolution is used to obtain the quantum speed limit bound ratio $T_{\\mathrm{QSLT}}/T$. The paper finds that this ratio is suppressed for larger black hole spin $a$, and similarly suppressed for lower primordial black hole mass, meaning stronger spin-curvature coupling allows the neutrino flavor state to evolve faster. The entanglement entropy of the four-qubit flavor state is suppressed in the strong-coupling region and reaches Bell-state-like maxima where the survival probability is near $1/2$.","pith_inferences":["One could test the WKBJ assumption directly by numerically solving the Dirac equation in the Kerr background without the local-constant approximation; the predictions near $r \\sim 3\\,r_g$, where $B^d$ varies fastest, are the most vulnerable.","Because the gravity-induced effects are confined to roughly $100$ gravitational radii, any realistic astrophysical probe would need neutrinos produced or detected very close to the horizon; distant detectors would see ordinary vacuum oscillations.","The same Kerr-Schild polar-coordinate four-vector potential could also be used to derive gravitational geometric phases or spin precession for other fermions around spinning black holes, not just neutrino flavor oscillations.","A future measurement of solar neutrino-number asymmetry at the predicted $10^{-39}$ level would confirm the gravitational Zeeman mechanism, while any detection near the observed $10^{-10}$ level from the Sun would rule it out as the dominant cause."],"forward_implications":["Near a fast-spinning primordial black hole the neutrino-antineutrino oscillation length is short and the survival probability of the initial flavor state is high close to the horizon, so oscillations are suppressed in the strong-coupling region and recover farther out.","Larger black hole spin $a$ lowers the quantum speed limit ratio $T_{\\mathrm{QSLT}}/T$, so the flavor state can evolve faster close to the hole; lower primordial black hole mass has the same effect.","Entanglement entropy of the neutrino-antineutrino flavor mode is strongly suppressed near the hole and returns to Bell-state-like maxima where the survival probability equals $1/2$, then becomes independent of angle and spin beyond about $100$ gravitational radii.","For extremely long distances the gravitational effect dies out and the two-flavor oscillation approaches ordinary vacuum mixing, independent of black hole spin and angle.","The same gravitational Zeeman effect generates a neutrino-number asymmetry; at the solar surface the paper estimates $\\Delta n/n \\sim 10^{-39}$, far below the observed cosmic value near $10^{-10}$, so a dedicated solar probe would be needed to see it."],"supporting_citations":[{"why":"Introduces gravity-induced neutrino-antineutrino oscillation and the gravitational Zeeman effect that this paper builds on.","marker":"[21]"},{"why":"Supplies the modified mass matrix formalism for neutrinos coupled to gravity, which Section V adapts to the two-flavor case.","marker":"[22]"},{"why":"Provides the Kerr metric in Kerr-Schild polar coordinates and the tetrad choice used to derive the four-vector potential.","marker":"[33]"},{"why":"Defines the generalized quantum speed limit ratio used throughout the QSLT analysis.","marker":"[35]"},{"why":"Earlier study of gravity's influence on the quantum speed limit in neutrino oscillations, which this paper extends to Kerr-Schild polar coordinates.","marker":"[31]"},{"why":"Establishes the occupation-number four-qubit representation and entanglement entropy technique used in Section IX.","marker":"[54]"},{"why":"Supplies the WKBJ condition that justifies treating the gravitational potential as locally constant.","marker":"[58]"},{"why":"Provides the neutrino mixing parameters and mass splittings used to fix the Majorana masses and mixing angles.","marker":"[59]"}],"fun_headline_variants":["Black hole spin speeds neutrino flavor flips near horizon","Spinning black holes alter neutrino oscillation speed limits","Curved spacetime boosts neutrino quantum speed limit","Kerr black hole gravity tunes neutrino entanglement","Spin-curvature coupling revs neutrino flavor evolution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the WKBJ approximation: the neutrino wavelength is so small compared with the scale on which the gravitational potential varies that the potential can be treated as locally constant and the global solution built by stitching local plane waves; if this fails near the black hole, the derived survival probabilities and speed limits break down.","fun_headline_variants_meta":{"raw":{"variants":["Black hole spin speeds neutrino flavor flips near horizon","Spinning black holes alter neutrino oscillation speed limits","Curved spacetime boosts neutrino quantum speed limit","Kerr black hole gravity tunes neutrino entanglement","Spin-curvature coupling revs neutrino flavor evolution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000158,"raw_usage":{"total_tokens":1252,"prompt_tokens":1000,"completion_tokens":252,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":181}},"tokens_in":616,"tokens_out":252,"duration_ms":2965,"temperature":1.0,"reasoning_tokens":181,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:34:10.741167+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the neutrino survival probability near a $10^{18}$ kg spinning primordial black hole by numerically integrating the Dirac equation in the Kerr background without assuming a locally constant $B^d$, for $a = 0.998$ and $\\theta = \\pi/3$. If the exact result differs materially from the paper's Eq. (32) for radii between about $3$ and $6$ gravitational radii, the WKBJ stitching assumption fails and the claimed suppression of oscillations near the hole is not established.","supporting_citations":[{"cited_title":"Influence of gravity on the quantum speed limit in neutrino oscillations","cited_arxiv_id":"2411.18558","evidence_quote":"Provides the Kerr metric in Kerr-Schild polar coordinates and the tetrad choice used to derive the four-vector potential."},{"cited_title":"Space-time curvature coupling of spinors in early universe: Neutrino asymmetry and a possible source of baryogenesis","cited_arxiv_id":"hep-ph/0510351","evidence_quote":"Earlier study of gravity's influence on the quantum speed limit in neutrino oscillations, which this paper extends to Kerr-Schild polar coordinates."},{"cited_title":"Complete complementarity relations for quantum correlations in neutrino oscillations","cited_arxiv_id":"2205.01601","evidence_quote":"Establishes the occupation-number four-qubit representation and entanglement entropy technique used in Section IX."},{"cited_title":"Schwinger, Particles, Sources, And Fields, V olume 1 (Taylor and Francis, 2019)","cited_arxiv_id":null,"evidence_quote":"Provides the neutrino mixing parameters and mass splittings used to fix the Majorana masses and mixing angles."}],"review_version":1}