{"id":"8dfa16d6-c0bb-4021-8d15-cec1ef2008c1","arxiv_id":"2607.17825","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A U(1) gauge model with six chiral singlet fermions produces naturally small Dirac neutrino masses plus smaller Majorana masses, giving a quasi-Dirac spectrum and a TeV-scale new boson.","lead":"This paper builds a particle physics model in which neutrinos get tiny masses from a new gauge force, making them a hybrid called quasi-Dirac neutrinos. It is a theory proposal with specific predictions for a heavy new boson, a dark matter candidate, and very slow neutrino decays, but no experimental data yet.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quasi-Dirac splitting R~10^-3 eV gives δm²~10^-4 eV² and appears excluded by solar neutrino data; the paper never confronts this, so 'realistic neutrino physics' is not established.","rationale":"I read the paper in good faith as a construction of quasi-Dirac neutrinos from a gauged U(1)Lμ−Lτ with anomaly-free chiral singlets. The anomaly cancellation, the Landau-pole estimate, and the neutrino mass-matrix algebra in Eqs. (8)–(22) are internally consistent, and the illustrative fits to oscillation data are plausible as existence arguments. The reader's weakest assumption, the hierarchical scalar VEV pattern, is real: the cross-quartic terms λk4|φk|²|φ4|² in Eq. (5) force a severe cancellation if vk∼GeV while v4∼tens of TeV, and the paper does not analyze the potential minimum. However, the most load-bearing soft spot is more direct: the model's natural quasi-Dirac splitting appears to violate solar neutrino data. With R3×3∼10^-3 eV and active masses ∼10^-2 eV, the mass-squared splitting within each quasi-Dirac pair is δm²∼10^-4 eV², far above the decoherence threshold for solar neutrinos, so half of the active flux would be lost to sterile components. The paper instead checks only N6 stability and collider/BBN constraints, not the active-neutrino survival probability. A smaller R could evade the bound, but then the 'natural' quasi-Dirac character is lost, so the central claim requires a quantitative benchmark that passes solar data. This reinforces the conditional verdict rather than overturning it.","tokens_in":12696,"tokens_out":28388,"duration_ms":278356,"concrete_test":"Using the quoted benchmark (D3×3 corresponding to Eqs. (14)–(16), R3×3 with entries yN v²/M∼10^-3 eV, yN∼O(1), v2∼1–10 GeV, M=10^14 GeV), construct the full 6×6 neutrino mass matrix of Eq. (18), compute the solar ν_e survival probability including the quasi-Dirac pairs with decoherent averaging, and compare with Borexino/SNO data. If the low-energy P_ee is roughly half the standard MSW value (≈0.27 instead of ≈0.55), the natural parameter point is excluded; then repeat with R reduced until P_ee is within 2σ of the data and check whether the required yN or v_i values remain natural.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing gap is the unaddressed phenomenology of the quasi-Dirac splitting itself. For the stated R3×3∼10^-3 eV (Eq. (20)) and Dirac masses m_i∼10^-2 eV (Eq. (19)), each active-sterile pair is split by δm²≈2m_i R_ii∼(2×10^-5–1×10^-4) eV². Since the active and sterile components enter with equal weight (1/√2), once this splitting is resolved the survival probability of solar ν_e is reduced by roughly half relative to the standard three-neutrino MSW expectation. The solar baseline and energies (L∼1.5×10^11 m, E∼0.1–10 MeV) average over δm²≫4E/L∼10^-12 eV², so the suppression is maximal. SNO and Borexino data determine the low-energy survival probability within the standard framework; halving it is excluded. The paper quotes R from natural yN∼O(1) and v2∼1–10 GeV, but never computes the resulting solar survival probability or applies quasi-Dirac oscillation bounds despite citing Ref. [9]. The scalar VEV hierarchy identified by the reader is a further gap, but this one directly attacks the abstract's claim that 'realistic neutrino physics can be produced' for the quoted natural parameters.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constructs an L_mu - L_tau gauge extension of the Standard Model with six chiral SM-singlet fermions and four dark Higgs fields. The new U(1) forbids renormalizable neutrino Yukawa and Majorana mass terms, so active neutrino Dirac masses arise from dimension-5 operators at a high scale M, giving m_nu ~ 10^-2 eV with O(1) couplings, while the right-handed neutrino Majorana masses from the same operators are R ~ 10^-3 eV. This produces a quasi-Dirac neutrino spectrum with approximate lepton number conservation. The paper illustrates that the Dirac mass matrix can accommodate the measured neutrino masses and mixings, and it discusses X-boson collider phenomenology, muon g-2, the BBN bound on m_X, and a long-lived sub-GeV dark matter candidate N6.","tokens_in":13115,"tokens_out":15185,"duration_ms":140501,"significance":"If the phenomenological issues were resolved, the model would be an elegant proof of principle: a gauged symmetry, rather than a global symmetry, explains the smallness of Dirac neutrino masses while also producing a calculable quasi-Dirac splitting and testable signals at a muon collider. The charge assignments and mass-matrix decompositions in Section II are internally consistent, and the model makes falsifiable predictions for m_X and for the decay channels of N6. However, the advertised claim that 'realistic neutrino physics can be produced' is not backed by a fit, the natural benchmark is in strong tension with solar neutrino data, and the scalar VEV hierarchy is assumed rather than demonstrated. These gaps are load-bearing, so the paper is not yet suitable for publication.","major_comments":[{"comment":"The advertised natural parameter choice is excluded by solar neutrino data, and the paper never confronts this. With D_3x3 ~ 10^-2 eV and R_3x3 ~ 10^-3 eV, each active-sterile pair has a mass splitting delta m^2 ~ 2 m_i R_ii between roughly 2 x 10^-5 and 1 x 10^-4 eV^2. Because each mass eigenstate contains active and sterile components with equal weight, once this splitting is resolved the solar nu_e survival probability is reduced from the standard three-neutrino value to approximately (1/4) Sigma_i |U_ei|^4 for vacuum-averaged propagation, and it is similarly suppressed in the MSW regime. This is in strong conflict with SNO and Borexino. The paper cites Ref. [9] on quasi-Dirac oscillations but does not apply the resulting bounds, so the abstract's claim that 'realistic neutrino physics can be produced' is not supported for the quoted parameters. A constraint on R_3x3, or an explicit demonstration that the splitting can be made small enough to evade solar bounds without destroying the naturalness of the model, is required.","section":"II.A, Eqs. (19)-(20)"},{"comment":"The hierarchical VEV pattern v1, v3 ~ 10 GeV, v2 ~ 1-10 GeV, and v4 of tens of TeV is assumed rather than derived. The scalar potential in Eq. (5) contains quartic and trilinear terms with arbitrary coefficients, and the paper does not show that a minimum with this hierarchy exists, nor that the hierarchy is stable against radiative corrections. This hierarchy sets the Dirac mass scale, the quasi-Dirac splitting, m_X, m_N6, and the BBN safety of the model, so it is load-bearing; without an existence or naturalness argument for this VEV pattern, the central mechanism is incomplete.","section":"II, Eq. (5) and Eq. (7)"},{"comment":"The claim that realistic neutrino physics can be produced is not backed by an explicit fit. The authors display target matrices for the combination U diag(m_i^2) U^T, but they never exhibit Yukawa matrices y^nu and VEV ratios that reproduce these matrices, nor do they quantify the required tuning in the row vectors y^nu_i. Since the right-hand side of Eq. (10) has the factorized structure v_i v_j y^nu_i y^nu_j^dagger, it is not self-evident that the measured PMNS matrix and mass splittings can be obtained without severe fine-tuning of the Yukawa vectors. An explicit numerical example, with all parameters specified, would substantiate the abstract's central claim.","section":"II.A, Eqs. (10), (14)-(16)"}],"minor_comments":[{"comment":"The (2,3) and (3,2) entries of the matrix in Eq. (17) are printed with inconsistent notation, 'y^N_56 v2' and 'y2_56 v2'; the latter should presumably be 'y^N_56 v2'.","section":"Eq. (17)"},{"comment":"The text first fixes m_X ~ (5-10) TeV in Section III and later requires m_X >= 43 TeV for BBN safety; the final adopted benchmark should be stated once and consistently.","section":"III.C"},{"comment":"There are several typos and incomplete reference entries, including 'namly' and 'betweem' in Section III.C, and missing publication years for Refs. [20], [21], and [29].","section":"Throughout"},{"comment":"The dark matter discussion would be strengthened by an estimate of the thermal relic abundance of N6; longevity alone does not establish that N6 has the observed dark matter density.","section":"III.C, dark matter"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and I found no problematic citation pattern. My main concern is that the authors may respond to the solar neutrino bound by lowering R without addressing the resulting fine-tuning in y^N, which would undermine the naturalness motivation of the model. I recommend that the revision be required to show a viable parameter region, including a solar-oscillation constraint, and to either justify the scalar VEV hierarchy or explicitly state it as a tuning assumption."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper gives a clean L_mu-L_tau gauge model with six chiral SM singlets, higher-dimensional operators producing small Dirac masses, and a quasi-Dirac spectrum. The construction is coherent and the anomaly math checks out. There is a novel combination here: the six-fermion charge assignment, the quasi-Dirac hierarchy, and the N6 sub-GeV dark matter candidate with suppressed decays.\n\nBut the central claim that 'realistic neutrino physics can be produced' is not established and, for the quoted natural parameters, looks wrong. With D ~ 1e-2 eV and R ~ 1e-3 eV, each active-sterile pair has mass-squared splitting about 4e-5 eV^2. That is huge compared to the solar oscillation baseline, so the electron neutrino survival probability is suppressed by roughly half relative to standard MSW once the splitting is resolved. The paper cites Anamiati et al. on quasi-Dirac bounds but never applies them to their own R. This is a load-bearing omission, not a side detail.\n\nThe paper also assumes the scalar potential has a minimum with v4 tens of TeV and v1,v3 ~ 10 GeV, v2 ~ 1-10 GeV, with no demonstration that such a hierarchy is natural or even exists. The entire mass pattern and the X boson mass depend on it. In addition, the neutrino fit is illustrative: they show matrices that give the right oscillation angles, but it is a tuned example, not a fit, so calling it 'realistic' is optimistic.\n\nOn the plus side, the paper is honest about what it cannot do: the muon g-2 prediction is far below the observed anomaly, and the BBN bound raises mX above 43 TeV, taking the collider reach beyond near-term muon colliders. The Landau pole analysis is clear, and the decay-length estimate for N6 is a sensible touch.\n\nSo: a well-constructed model with a serious phenomenological gap. The quasi-Dirac splitting needs to be reconciled with solar data; either R has to be much smaller (say <1e-12 eV) or the active-sterile mixing must be suppressed, which would require abandoning the naturalness argument that is the paper's main point. I would send it to peer review because the construction deserves scrutiny and a referee could ask for the fixing analysis. But as it stands, the abstract's promise is not met. I would not cite it in my own work right now. It's a decent reading-group candidate for a discussion of quasi-Dirac constraints.\n\nRecommended: peer review, with expectation of major revision or rejection if the solar bound cannot be answered.","headline":"A coherent quasi-Dirac neutrino model that fails to confront its own active-sterile splitting against solar data; the natural parameter choice is likely excluded.","tokens_in":13609,"tokens_out":5166,"would_cite":false,"duration_ms":45188,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq","12.60.-i","95.35.+d"],"model":"deepseek-v4-flash","headline":"A gauged $\\mathrm{U}(1)_{L_\\mu-L_\\tau}$ model with six chiral singlet fermions naturally produces quasi-Dirac neutrinos, with $10^{-2}$ eV Dirac masses from dimension-five operators and one-order-smaller Majorana masses.","keywords":["quasi-Dirac neutrinos","gauged L_mu-L_tau symmetry","anomaly-free chiral fermions","dimension-five operators","neutrino mass and mixing","sub-GeV dark matter","X boson phenomenology"],"falsifier":"Minimise the scalar potential of Eq. (5) numerically over its couplings: if no open region of parameter space yields $v_1\\sim v_3\\sim10$ GeV, $v_2\\sim1$–10 GeV, and $v_4\\sim30$ TeV, the hierarchy assumption is untenable. Independently, a positive observation of neutrinoless double-$\\beta$ decay at an effective Majorana mass above roughly $10^{-3}$ eV would contradict the paper's quasi-Dirac parameter region, since in that region lepton-number violation is bounded by $R_{3\\times3}\\sim10^{-3}$ eV.","tokens_in":12518,"feed_emoji":"⚛️","tokens_out":12659,"duration_ms":101524,"temperature":0.7,"pith_summary":"This paper argues that a specific gauge extension of the Standard Model can naturally produce quasi-Dirac neutrinos, without invoking global or discrete symmetries. The new ingredient is a $\\mathrm{U}(1)_{L_\\mu-L_\\tau}$ gauge symmetry under which six chiral Standard-Model-singlet fermions carry anomaly-free charges; three of them serve as right-handed neutrinos. The symmetry forbids renormalizable neutrino mass terms, so Dirac masses only appear through dimension-five operators suppressed by an assumed high scale $M\\sim10^{14}$ GeV, making the smallness of neutrino masses natural. After spontaneous breaking, the active-neutrino Dirac mass matrix is of order $10^{-2}$ eV and the right-handed Majorana mass matrix is about one order smaller, so the light spectrum forms three quasi-Dirac pairs with approximate lepton-number conservation. The paper shows that this structure can reproduce the measured neutrino mass splittings and lepton mixing, and that the accompanying $X$ boson, dark-matter candidate, and cosmology are consistent with current constraints when the new-physics scale is pushed to tens of TeV.","feed_headline":"Gauged symmetry makes neutrinos quasi-Dirac without tiny couplings","feed_subtitle":"Anomaly-free chiral singlets and dimension-five operators yield 0.01 eV Dirac masses with tiny lepton-number breaking.","key_machinery":"The load-bearing object is the $9\\times9$ neutral-fermion mass matrix, whose upper $6\\times6$ block is the quasi-Dirac matrix built from a Dirac block $D_{3\\times3}\\sim10^{-2}$ eV and a sterile-Majorana block $R_{3\\times3}\\sim10^{-3}$ eV. These scales are fixed by the assumed vacuum-expectation-value hierarchy $v_1\\sim v_3\\sim10$ GeV, $v_2\\sim1$–10 GeV, and $v_4$ at tens of TeV, together with the high suppression scale $M\\sim10^{14}$ GeV. The $\\mathrm{U}(1)$ gauge symmetry is the mechanism that makes the small entries natural: renormalizable mass terms are absent by charge conservation, and every physically relevant mass comes from a dimension-five operator, so the dimensionless couplings can be of order one while masses come out at the $10^{-2}$ eV scale.","core_discovery":"The paper's central claim is that a $\\mathrm{U}(1)_{L_\\mu-L_\\tau}$ gauge symmetry, broken by four dark Higgs fields with hierarchical vacuum expectation values, naturally yields quasi-Dirac neutrinos. Anomaly freedom fixes six chiral Standard-Model-singlet Weyl fermions with charges $N_{1,2,3}(2z)$, $N_{4,5}(-8z)$, and $N_6(10z)$; the same charges forbid tree-level Yukawa couplings of the Standard Model Higgs to right-handed neutrinos, so Dirac masses arise only from quadrilinear dimension-five operators such as $\\varphi_i\\tilde H^\\dagger L_i N_j/M$ with $M\\sim10^{14}$ GeV. With $v_1\\sim v_3\\sim10$ GeV and $v_2\\sim1$–10 GeV, the active Dirac mass matrix $D_{3\\times3}\\sim10^{-2}$ eV is large enough for realistic oscillations, while the right-handed Majorana matrix $R_{3\\times3}\\sim10^{-3}$ eV, generated by $\\varphi_2^2$ and $\\varphi_1\\varphi_3$ terms, is an order of magnitude smaller; the remaining $N_{4,5,6}$ sector is seesaw-suppressed and does not disturb the quasi-Dirac pairs. The paper demonstrates with explicit parameter choices that both neutrino mass orderings and the observed mixing angles can be fitted, identifies $N_6$ as a sub-GeV dark-matter candidate with a lifetime far exceeding the age of the Universe, and derives collider and Big-Bang-nucleosynthesis constraints on the $X$-boson mass.","pith_inferences":["Editorial extension: the same logic would work for any anomaly-free chiral $\\mathrm{U}(1)$ charge assignment over Standard-Model singlets, not just $L_\\mu-L_\\tau$, so the construction actually defines a family of quasi-Dirac models with different flavour patterns.","Editorial extension: if the high scale $M\\sim10^{14}$ GeV is identified with a unification or high-energy breaking scale, the $10^{-2}$ eV neutrino mass becomes a derived relation between that scale and the TeV-scale VEVs; the paper does not pursue this identification.","Editorial extension: the radiative decay $N_6\\to\\nu+\\gamma$ at energy around $0.1$ GeV is a concrete spectral signature whose non-observation could bound $\\varepsilon$ and the $N_6$ relic density, complementing the paper's discussion of its lifetime.","Editorial extension: a systematic numerical scan over the Yukawa matrices and VEVs could reveal how much of the lepton-mixing parameter space the model covers and whether the CP-violating phase is predicted; the paper gives two worked examples with $\\delta=1.19\\pi$ but no full scan."],"forward_implications":["Neutrinos behave as effectively Dirac particles: neutrinoless double-beta decay is unobservable at upcoming sensitivities because the lepton-number-violating Majorana block is only about $10^{-3}$ eV.","The model predicts a new $X$ gauge boson with mass at least around 43 TeV once Big-Bang nucleosynthesis constraints are enforced, together with new four-fermion interactions concentrated in the muon/tau/neutrino sector; these can be probed at a future muon collider.","The sub-GeV state $N_6$ is a viable dark-matter candidate with lifetime much longer than the age of the Universe, and its three-body decays can inject electrons and positrons in dense regions such as the Galactic center.","The heaviest neutrino's decay lifetime is about $6.6\\times10^{52}$ seconds, so neutrinos are effectively stable in all observable settings.","The model's $Z$--$X$ mixing and muon anomalous magnetic moment contributions are tiny, with $\\Delta a_\\mu\\simeq9.3\\times10^{-11}$, so precision electroweak and muon-anomaly measurements do not currently distinguish it from the Standard Model."],"supporting_citations":[{"why":"introduces the quasi-Dirac neutrino scenario as Dirac masses with small Majorana perturbations.","marker":"[8]"},{"why":"supplies the phenomenological framework of quasi-Dirac neutrino oscillations that the model aims to realise.","marker":"[9]"},{"why":"defines the Dirac seesaw naturalness problem that motivates suppressing Dirac masses by a high scale rather than tiny couplings.","marker":"[10]"},{"why":"provides the $\\mathrm{U}(1)_{L_\\mu-L_\\tau}$ gauge symmetry whose charge assignment is the model's starting point.","marker":"[15]"},{"why":"gives the general solution of the $\\mathrm{U}(1)$ anomaly equations, so the six chiral fermion charges are known to be consistent.","marker":"[19]"},{"why":"supplies the world-average neutrino mass splittings and mixing parameters used to test the Dirac mass matrix.","marker":"[24]"},{"why":"provides the Big-Bang-nucleosynthesis bound requiring the new gauge boson to be heavier than about 43 TeV, which sets $v_4$.","marker":"[29]"},{"why":"gives the lower bound on the new gauge boson mass from dilepton resonance searches, which the model's $X$ mass must satisfy.","marker":"[25]"}],"fun_headline_variants":["Gauged L_mu - L_tau symmetry yields quasi-Dirac neutrinos","Quasi-Dirac neutrinos from gauged L_mu - L_tau with dark matter","Anomaly-free gauge symmetry makes neutrinos quasi-Dirac","Natural quasi-Dirac neutrinos from L_mu - L_tau gauge symmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes, without demonstrating it, that the scalar potential in Eq. (5) has a minimum with the hierarchical vacuum expectation values $v_1\\sim v_3\\sim10$ GeV, $v_2\\sim1$–10 GeV, and $v_4$ of tens of TeV, with $v_4$ much larger than the others; this hierarchy sets every mass scale in the neutrino and dark-matter sectors, and if it cannot be realised the construction collapses.","fun_headline_variants_meta":{"raw":{"variants":["Gauged L_mu - L_tau symmetry yields quasi-Dirac neutrinos","Quasi-Dirac neutrinos from gauged L_mu - L_tau with dark matter","Anomaly-free gauge symmetry makes neutrinos quasi-Dirac","Natural quasi-Dirac neutrinos from L_mu - L_tau gauge symmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3493,"prompt_tokens":983,"completion_tokens":2510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":2419}},"tokens_in":599,"tokens_out":2510,"duration_ms":14406,"temperature":1.0,"reasoning_tokens":2419,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:34:39.315528+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Minimise the scalar potential of Eq. (5) numerically over its couplings: if no open region of parameter space yields $v_1\\sim v_3\\sim10$ GeV, $v_2\\sim1$–10 GeV, and $v_4\\sim30$ TeV, the hierarchy assumption is untenable. Independently, a positive observation of neutrinoless double-$\\beta$ decay at an effective Majorana mass above roughly $10^{-3}$ eV would contradict the paper's quasi-Dirac parameter region, since in that region lepton-number violation is bounded by $R_{3\\times3}\\sim10^{-3}$ eV.","supporting_citations":[{"cited_title":"Demonstration of Single Barium Ion Sensitivity for Neutrinoless Double Beta Decay using Single Molecule Fluorescence Imaging","cited_arxiv_id":"1711.04782","evidence_quote":"introduces the quasi-Dirac neutrino scenario as Dirac masses with small Majorana perturbations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the Dirac seesaw naturalness problem that motivates suppressing Dirac masses by a high scale rather than tiny couplings."},{"cited_title":"Dirac Neutrinos and Dark Matter within a Minimal Discrete Symmetry Model","cited_arxiv_id":"2408.14166","evidence_quote":"provides the $\\mathrm{U}(1)_{L_\\mu-L_\\tau}$ gauge symmetry whose charge assignment is the model's starting point."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the general solution of the $\\mathrm{U}(1)$ anomaly equations, so the six chiral fermion charges are known to be consistent."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the lower bound on the new gauge boson mass from dilepton resonance searches, which the model's $X$ mass must satisfy."}],"review_version":2}