{"id":"f1abec16-36f4-42f9-b3b8-d9cb85205f1b","arxiv_id":"1908.04855","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new first-order perturbation theory for solar-scale neutrino oscillations with non-unitary mixing shows that δ and the non-unitarity α parameters form physical, convention-independent phase correlations.","lead":"This paper works out how neutrinos oscillate if their mixing is slightly non-unitary, focusing on the energy range where the solar mass splitting dominates. It finds that the CP phase and the non-unitarity parameters form correlations that cannot be removed by changing phase conventions, which future low-energy neutrino experiments will need to handle.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SOL-convention δ-α correlation is proven only at first order, while the paper's exact solar-region numerics (Sec. 7.4) show patterns the analytic framework cannot explain; the physical-reality claim consequently lacks direct support in the displayed regime.","rationale":"The reader's weakest assumption — that the correlation is derived only in first-order perturbation theory and is not reproduced by the exact numerics in the displayed solar-region regime — is precisely the most load-bearing point. The paper is honest about this limitation in Sec. 7.4, conceding that the coexistence of vertical and circular phase-correlation patterns is not understood by the analytic framework. That concession directly undermines the strongest positive claim: the specific form of the correlation (e^{±iδ} with α-blobs) has not been verified against exact calculations in the solar region. The paper's argument that the correlation 'prevails to higher order' applies to the S-matrix level and does not automatically survive the |S|^2 and (1-α) projections needed for the physical probability. Nevertheless, the paper contains a detailed, internally consistent first-order derivation, symmetry checks, and an explicit translation rule between U_MNS conventions, all of which support the conditional value of the work. No fraud or circular reasoning is present. A single numerical test at α = 10^-3 would settle whether the concern is real or whether the first-order formula is simply a victim of the large-α display choice. The CONDITIONAL verdict is therefore appropriate and unchanged.","tokens_in":37424,"tokens_out":4332,"duration_ms":44838,"concrete_test":"Compute the exact numerical ΔP_μe in the φ_μe-δ and φ_τe-δ planes at E = 200 MeV, L = 3000 km and L = 12000 km, using α_μe = α_τe = 10^-3 (instead of 0.1), with all other α set to zero and constant density ρ = 3.2 g/cm^3. Overlay the first-order analytic prediction from Eqs. (5.7), (5.8), (6.3) using the same parameters. If the contour shapes coincide, the first-order correlation is the physical one in the perturbative regime; if vertical/circular structures persist at α = 10^-3, the analytic identification fails and the existence claim needs a nonperturbative proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central existence claim in Sec. 6.2 — that the δ-α correlation appears in the SOL convention at solar-scale enhancement — is derived from the first-order perturbative probabilities of Sec. 5. The statement in Sec. 6.2 that the correlation 'prevails to higher order' because it resides in the Φ matrix elements covers only the unitary-evolution S-matrix elements; the physical flavor-basis probability (5.1) includes the non-unitary projectors (1-α) and interference of orders, so higher-order survival is not demonstrated. More importantly, the exact numerical phase-correlation plots in Sec. 7.4 are made with α = 0.1, which the paper itself notes is outside the perturbative regime, and the observed vertical and circular φ-δ patterns are explicitly stated (Sec. 7.4) to be 'not understood' by the analytic framework. Thus the precise form of the correlation — e^{±iδ} times α-blobs — has not been confirmed by exact calculations in the solar region, and the physical-reality conclusion (Sec. 6.3) rests on an unverified truncation. If the exact correlation at realistic α differs in its δ-dependence, the claim that no UMNS convention removes the correlation in both oscillation regions would need qualification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the authors' previous 'helio-UV' perturbation theory to the solar-resonance region for three-neutrino evolution with a non-unitary mixing matrix in the SOL convention of UMNS, in which e^{±iδ} is attached to s12. It derives first-order expressions for the νμ→νe oscillation probability, decomposing the non-unitary correction into a unitary evolution part and a genuine non-unitary part. The central physics claim is that a δ−α parameter correlation does exist in the SOL convention at solar-scale enhanced oscillations, taking the form of e^{±iδ} correlated with 'blobs' of α parameters (K12 and K23), rather than the chiral combinations α̃_βγ e^{±iδ} found in the atmospheric region in the PDG convention. The authors conclude that no UMNS convention removes the correlation in both the atmospheric and solar regions, and hence that the δ−α correlation is physical rather than a convention artifact. Sections 7.2–7.4 add exact numerical studies of the UV contribution ΔPμe, showing cancellations between the EV and UV parts and among different α parameters, and presenting φβγ−δ correlation plots in the solar and atmospheric regions.","tokens_in":37609,"tokens_out":3063,"duration_ms":36188,"significance":"If the central claim is correct, the paper resolves a genuine open question from the companion work [38]: whether the δ−α phase correlation is a phase-convention artifact. It also provides a new analytic tool, the UV-extended solar-resonance perturbation theory, which is likely to be useful for future low-energy atmospheric-neutrino unitarity tests. The manuscript is technically careful: the first-order derivation is presented with full appendices, the ϕ→ϕ+π/2 dynamical symmetry is used as a nontrivial consistency check, and exact numerical integration is used to cross-check qualitative features. The authors also deserve credit for explicitly disclosing where their analytic framework fails, notably the vertical and circular φ−δ patterns in the solar region that are stated in Sec. 7.4 to be 'not understood, regrettably, by our analytic framework.' These strengths make the paper a serious contribution even though the physical-reality conclusion needs additional support at the quantitative level.","major_comments":[{"comment":"The physical-reality conclusion in Sec. 6.3 rests on the first-order perturbative formula, but the exact numerical phase-correlation plots that directly probe the solar region are made with αβγ = 0.1, which footnote 20 and the surrounding discussion admit is outside the perturbative regime. The observed vertical and circular φ−δ patterns are explicitly stated in Sec. 7.4 to be 'not understood' by the analytic framework, and the figures are computed in the PDG convention, not the SOL convention in which the central claim of Sec. 6.2 is formulated. The paper should either (i) show that the first-order e±iδ−K12/K23 correlation reproduces the exact numerical ΔPμe for α values within the existing bounds (e.g., the values used in Fig. 1), or (ii) explicitly qualify the claim that the correlation is physical as a first-order statement and discuss how higher-order corrections could modify the δ-dependence.","section":"Sec. 7.4"},{"comment":"The statement that the δ−(blob of α) correlation 'prevails to higher order in perturbation theory in the unitary evolution part' is demonstrated only for the Φ-matrix elements that enter the tilde-basis S-matrix elements. The physical flavor-basis probability (5.1) is built from S_flavor = (1−α)S_prop(1−α)†, whose modulus squared contains the non-unitary projectors and interference between S(0), S(1)_EV, and the αS(0)S(0)α† terms. Showing that Φ contains the correlation does not, by itself, show that the probability to higher order preserves exactly the same e±iδ-blob form; the authors should either prove this for the probability or soften the higher-order claim.","section":"Sec. 6.2"},{"comment":"The conclusion that 'there is no UMNS convention in which the phase correlation is absent both at around the atmospheric- and the solar-scale enhanced oscillations' is proven only within the first-order perturbative framework. Since the exact solar-region numerics display features not reproduced by this framework, the universal conclusion over all conventions would be on firmer ground if the authors demonstrated, for at least one representative realistic α value, that the exact probability in the SOL convention contains a δ-dependent correlation of the predicted sign and approximate magnitude.","section":"Sec. 6.3 / Sec. 7.4"}],"minor_comments":[{"comment":"There is a typographical error in Eq. (5.6): the terms 'P_EV|OD1' and 'P_EV|OD2' are printed without the intervening plus sign, making the decomposition hard to read.","section":"Eq. (5.6)"},{"comment":"In the sentence beginning 'One may wonder why the features of the correlation between α and the α parameters are so different...', the first 'α' should presumably be 'δ'; please correct the wording.","section":"Sec. 6.3"},{"comment":"The figure captions for Figs. 3–5 do not state that the computations use the PDG convention, even though Sec. 7 states this at the start of the section; adding the convention to the captions would prevent confusion for readers who jump directly to the figures.","section":"Section 7"},{"comment":"The target-sensitivity discussion would benefit from a one-sentence reminder that the first-order UV expression is being used for the '10^{-2}' accuracy estimate, since the 10^{-4} target is also mentioned and the second-order α² terms are discussed in the same paragraph.","section":"Sec. 4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JHEP and is a legitimate continuation of the authors' prior work. The main issue is not the derivation but the support for the 'physical reality' conclusion: the exact numerics at α = 0.1 are admittedly outside the perturbative regime and are not reproduced by the analytics, while the realistic-α comparison is not shown. This is fixable by adding numerical checks at α values within current bounds or by making the scope of the conclusion explicit. I would not reject the paper; the central analytic result is substantial even if the stronger interpretive claim needs qualification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a serious referee. The new piece is the extension of the solar-resonance perturbation theory to non-unitary mixing, and the demonstration that the δ–α phase correlation appears in the SOL convention of UMNS in the solar region. That directly answers the authors' own earlier puzzle from the atmospheric-region paper, where the SOL convention looked like the one place the correlation vanished. If it survives at solar-scale enhancement, the 'convention artifact' interpretation is dead, and the correlation is physical in the only sense that matters: no rephasing can remove it from both regions. That is a real step forward for the unitarity-violation program.\n\nThe calculation itself is solid. The first-order probability is derived explicitly, with F and K and Φ matrix bookkeeping, appendices, and a useful symmetry check under φ → φ + π/2. The decomposition into the unitary-evolution and genuine non-unitary parts is clean, and the observation that these two parts tend to cancel — both here and in the earlier atmospheric paper — is a genuine phenomenological insight with implications for how future solar-scale data (DUNE sub-GeV, Hyper-K, JUNO) should be fit. The paper is also honest about its limits: it says plainly in Sec. 7.4 that the vertical and circular phase-correlation patterns in the exact numerics are not understood by the analytic framework, and it flags the missing physical explanation for the EV–UV cancellation.\n\nWhere are the soft spots? First, the central existence claim is first-order. Section 6.2 says the correlation 'prevails to higher order' because it lives in the Φ matrix elements — but that argument covers the unitary S-matrix, not the full flavor-basis probability, which includes the (1−α) projectors and interference with the zeroth-order piece. So the higher-order survival is asserted, not proven. Second, the exact numerical checks in Figs. 3–5 are at α = 0.1, which the paper itself admits is outside the perturbative regime. Those plots show patterns that the first-order formula does not reproduce. That does not refute the first-order result — the numerics are consistent with a small-α expansion — but it does mean the displayed exact evidence does not directly confirm the e^{±iδ}-blob form in the regime where the analytic claim is controlled. Third, uniform matter density is assumed throughout; the paper notes adiabatic extension is possible, but the cancellations could shift with a realistic density profile. None of this is a load-bearing flaw; the first-order result stands on its own, and the limitations are disclosed. But the 'physical reality' conclusion is slightly stronger than the evidence displayed.\n\nWho is this for? Neutrino phenomenologists working on non-unitarity, parameter correlations, and future unitarity tests with solar-scale oscillations. They will want this as a reference and as a spur to higher-order and varying-density checks. I would send it to peer review with a request that the authors either prove the higher-order survival claim properly or soften it, and ideally show at least one exact-numerics panel at a genuinely perturbative α, like 0.01, to confirm the blob structure inside the regime of validity.","headline":"A genuinely new analytic result — the δ–α correlation survives in the SOL convention at solar-scale enhancement — but the 'physical reality' claim is proven only at first order, and the paper's own exact numerics at α = 0.1 show unexplained patterns.","tokens_in":38223,"tokens_out":1194,"would_cite":true,"duration_ms":15286,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq"],"model":"deepseek-v4-flash","headline":"This paper claims that the correlation between the Standard Model CP phase δ and the complex α parameters encoding unitarity violation in neutrino oscillations is physical, not a convention artifact, because it reappears in the SOL…","keywords":["neutrino oscillations","unitarity violation","non-unitary mixing matrix","CP phase correlation","solar resonance perturbation theory","alpha parameters","matter effects"],"falsifier":"Evaluate the exact all-orders appearance probability P(νμ→νe) in the SOL convention with αμe = 0.1 at E = 200 MeV and baselines 3000 km and 12000 km, as in the paper's figures, and check whether the δ dependence of ΔPμe survives when the second-order UV Hamiltonian and realistic varying matter density are included; if the δ dependence vanishes in that exact evaluation, the first-order correlation would be shown to be an artifact.","tokens_in":37147,"feed_emoji":"⚛️","tokens_out":6396,"duration_ms":62957,"temperature":0.7,"pith_summary":"The paper sets out to settle whether the phase correlation between the Standard Model CP phase δ and the complex α parameters that encode unitarity violation in neutrino oscillations is physically real or merely an artifact of the mixing-matrix convention. Working in the SOL convention of the MNS matrix, where $e^{{±iδ}}$ is attached to the solar angle s12, the authors build a first-order perturbative framework around solar-scale enhanced oscillations, extending their earlier solar-resonance perturbation theory to non-unitary mixing. They find that in the solar region δ and the α parameters are correlated after all, but through $e^{{±iδ}}$ multiplying composite blobs of α parameters rather than through simple chiral combinations. Because the correlation survives in the one convention where it had previously appeared absent, the paper concludes that no mixing-matrix convention wipes it out, so the correlation is physical. Numerically, the paper also finds that non-unitary effects tend to cancel between the unitary-evolution part and the genuine non-unitary part of the oscillation probability.","feed_headline":"Solar oscillations confirm the neutrino phase correlation","feed_subtitle":"A solar-resonance calculation ties the CP phase δ to unitarity-violation parameters in every mixing convention.","key_machinery":"The machinery is the solar-resonance perturbation theory extended to non-unitarity. It starts from the flavor-basis Hamiltonian with non-unitary mixing matrix N = (1 − \\tilde{α})U_SOL, transforms through tilde and hat bases to diagonalize the zeroth-order νSM Hamiltonian in matter, and treats both the matter-dressed effective parameter A_exp ≈ c13 s13 a/Δm²31 ~ $10^{{-3}}$ and the α parameters as small expansion parameters. The load-bearing objects are the F and K matrices, which repackage the α parameters, and the Φ matrix elements built from K and the diagonalized evolution phases; the correlated combinations K_{12}$e^{{-iδ}}$ and K_{23}$e^{{iδ}}$ emerge from Φ and carry the δ-(blob of α) correlation. The framework also produces a dynamical symmetry under φ → φ + π/2 that serves as a consistency check.","core_discovery":"In the solar-scale enhanced oscillation region, the νSM CP phase δ correlates with the unitarity-violating α parameters even in the SOL convention of UMNS, where $e^{{±iδ}}$ multiplies s12. The correlation is not of the 'chiral' form seen in the atmospheric region under the PDG convention ([$e^{{-iδ}}$\\bar{α}_{μe}, $e^{{-iδ}}$\\bar{α}_{τe}, \\bar{α}_{τμ}]); instead the first-order amplitudes contain K_{12}$e^{{-iδ}}$ and K_{23}$e^{{iδ}}$, where K_{12} and K_{23} are blobs built from the SOL-convention α parameters. Since no UMNS phase convention makes the correlation vanish in both the atmospheric and solar regions at once, the paper concludes that the δ-α correlation is physical rather than a convention artifact. The paper also reports that the exact numerical phase-correlation patterns in the solar region include vertical and circular contours that the first-order analytic framework does not capture, a gap it explicitly acknowledges.","pith_inferences":["If the δ-blob correlation survives exact all-orders computation, low-energy atmospheric neutrino detectors near the solar resonance could serve as independent probes of the CP phases of non-unitarity, complementing long-baseline experiments that probe the atmospheric region.","The vertical and circular φ-δ contours in the exact solar-region plots may be a second-order-in-α effect; computing the second-order UV correction would be a direct test of whether the first-order blob structure is the full story.","Since NSI parameters cluster into collective variables in a different pattern, a joint analysis of solar- and atmospheric-region appearance data could in principle separate non-unitarity from non-standard interactions by looking for these different correlation fingerprints."],"forward_implications":["A unitarity-violation fit that uses solar-scale enhanced oscillations must treat the Standard Model phase δ and the α parameters as correlated, because the correlation persists in the SOL convention.","The non-unitary contribution to P(νμ→νe) splits into a unitary-evolution part and a genuine non-unitary part that tend to cancel; a measurement of the appearance probability alone can therefore hide non-unitarity, so the departure-from-unitarity sum rule P(νμ→νe)+P(νμ→νμ)+P(νμ→ντ)≠1 becomes the more direct diagnostic.","Different α parameters also cancel against each other in the full ΔPμe, so bounds obtained by turning on one α at a time may be artificially strong compared with a global marginalization over all α parameters.","The form of the correlation changes from the atmospheric region to the solar region: only in the solar region does the e^{±iδ}-blob combination appear, so constraints and degeneracies derived for long-baseline experiments cannot be assumed to hold for low-energy atmospheric data."],"supporting_citations":[{"why":"Companion paper that found the chiral δ-α correlation in the atmospheric region and raised the SOL-convention skepticism that this paper answers.","marker":"[38]"},{"why":"The solar-resonance perturbation theory for the νSM sector that this paper extends by adding the α expansion.","marker":"[41]"},{"why":"Source of the α-parametrization of the non-unitary mixing matrix used throughout the paper.","marker":"[23]"},{"why":"Provides the α bounds used in the numerical figures and the NSI-UV correspondence that motivates the parameter-correlation discussion.","marker":"[26]"},{"why":"Gives the exact non-unitary evolution in constant matter density, used as the benchmark for exact numerical results and to justify the vacuum-mass-eigenstate evolution equation.","marker":"[25]"},{"why":"NSI perturbation theory whose collective variables Θ13 and Θ12 establish the precedent for clustering of νSM and new-physics parameters.","marker":"[55]"}],"fun_headline_variants":["New framework shows solar-scale CP-α correlation is physical","Solar oscillations tie CP phase to non-unitarity in all conventions","Unitarity-violation parameters correlate with CP phase in solar region","Neutrino CP-α correlation survives every phase convention","Solar-scale enhancement reveals physical δ-α correlation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analytic derivation of the correlation is first order in the small parameters α and A_exp and assumes constant matter density; the paper's exact numerical checks use α = 0.1, a value the paper acknowledges is outside the perturbative regime, and show φ-δ patterns the first-order formula cannot reproduce.","fun_headline_variants_meta":{"raw":{"variants":["New framework shows solar-scale CP-α correlation is physical","Solar oscillations tie CP phase to non-unitarity in all conventions","Unitarity-violation parameters correlate with CP phase in solar region","Neutrino CP-α correlation survives every phase convention","Solar-scale enhancement reveals physical δ-α correlation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1684,"prompt_tokens":1051,"completion_tokens":633,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":550}},"tokens_in":667,"tokens_out":633,"duration_ms":5910,"temperature":1.0,"reasoning_tokens":550,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:30:44.364241+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the exact all-orders appearance probability P(νμ→νe) in the SOL convention with αμe = 0.1 at E = 200 MeV and baselines 3000 km and 12000 km, as in the paper's figures, and check whether the δ dependence of ΔPμe survives when the second-order UV Hamiltonian and realistic varying matter density are included; if the δ dependence vanishes in that exact evaluation, the first-order correlation would be shown to be an artifact.","supporting_citations":[],"review_version":1}