{"id":"f7fd021a-1450-4677-9851-c5338b86638d","arxiv_id":"2501.15748","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A radiative seesaw model with non-holomorphic modular S3 symmetry fits neutrino data and predicts ranges for the Dirac CP phase, Majorana phases, neutrinoless double beta decay, and dark matter mass.","lead":"A model of neutrino masses using the Ma radiative seesaw mechanism and a non-holomorphic modular S3 flavor symmetry is analyzed. It fits current neutrino oscillation data and makes predictions for neutrino CP phases, neutrinoless double beta decay, and dark matter mass.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inverted-hierarchy scan points lie above the DESI+CMB bound the paper itself quotes, so the IH predictions should not be presented as viable.","rationale":"The reader's weakest_assumption is that the explicit non-holomorphic modular forms from ref [22] are not displayed, which is a legitimate reproducibility concern but secondary to the IH inconsistency. The reader's rationale does mention the IH tension with the DESI+CMB bound, but they do not elevate it to the weakest or most load-bearing issue. In my reading, the IH contradiction is the single most load-bearing concern because it is a quantitative, internally verifiable conflict with a bound the paper itself adopts. It does not rely on external forms or fitting conventions; it follows from the model's two-RHN structure, which forces one neutrino mass to zero, and the quoted cosmological limit. The IH predictions at roughly 100 meV are robust and always above 72 meV, so the IH cases cannot be accepted as viable. The NH central claim remains plausible, and the reader's CONDITIONAL verdict is appropriate. I therefore keep the verdict unchanged, but I disagree with the reader's choice of weakest assumption: the IH cosmological bound is the more decisive issue. The modular-forms reproducibility problem, the missing chi-square details, and the two-RHN assumption are real but would collectively affect the strength of the NH claim, not the validity of the IH sections. The proposed check directly settles whether the IH sections violate the quoted bound: recompute PDnu for the IH scan points and compare to 72 meV, which requires no external code or modular-form conventions.","tokens_in":9733,"tokens_out":11992,"duration_ms":105582,"concrete_test":"For each IH scan point in Figs. 8 and 11, compute PDnu = Dnu1 + Dnu2 (with Dnu3 = 0) using kappa_IH = sqrt(delta m^2_atm / D~nu2^2) from Eq. (II.16) and NuFit 6.0 central values for the mass-squared differences. If the resulting PDnu always exceeds 72 meV, then every IH point violates the DESI+CMB bound quoted in the paper, and the IH predictions must be removed or reframed as excluded. Also verify that the published IH chi-square regions lie entirely to the right of the 72 meV vertical line, confirming the inconsistency.","verdict_should_be":"UNCHANGED","load_bearing_attack":"With two right-handed neutrinos, one active neutrino mass is exactly zero (the paper states the lightest neutrino mass eigenvalue is zero). In the inverted hierarchy, m3 = 0, so the sum of neutrino masses is PDnu = m1 + m2, which is about 2*sqrt(delta m^2_atm) = 100 meV. The paper explicitly quotes the DESI+CMB combination PDnu <= 72 meV (ref [32]) and draws this bound as a vertical line in Figs. 8 and 11. Yet the IH scan points in those figures are localized near PDnu = 100 meV, to the right of the bound, and the text does not acknowledge that these points are excluded. Because the abstract claims to demonstrate predictions for both normal and inverted hierarchy, this is an internal inconsistency in the central claim: the IH cases violate the very cosmological constraint the authors use in the same section. The concern does not depend on modular form conventions or fitting details; it follows directly from the two-RHN structure and the quoted bound. The NH predictions (57–60 meV) are consistent with the bound, so the NH central claim is not affected. The IH sections should either be withdrawn or explicitly labeled as excluded, otherwise the paper overstates its viable parameter space.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a radiative seesaw model based on the Ma model, augmented with a non-holomorphic modular S3 flavor symmetry. Active neutrino masses are generated at one loop through an inert scalar doublet and two right-handed neutrinos, which forces one active neutrino mass to be exactly zero. The authors carry out a chi-square scan over the model parameters to fit the three neutrino mixing angles and two mass-squared differences from NuFit 6.0, and they derive predictions for the Dirac and Majorana phases, the neutrinoless double beta decay effective mass, the sum of neutrino masses, the dark matter mass, and lepton flavor violation rates, in four cases: normal or inverted hierarchy with fermionic or bosonic dark matter. For the normal hierarchy, the paper reports rather sharp predictions, including a Dirac CP phase between about 140 and 240 degrees, m_ee between about 2 and 4.4 meV, a neutrino mass sum between 57 and 60 meV, and a fermionic dark matter mass between 10 and 3600 GeV, with BR(mu -> e gamma) reaching the current experimental upper bound. For the inverted hierarchy, the paper reports m_ee at about 50 meV and a neutrino mass sum at about 100 meV, without noting that this sum exceeds the DESI+CMB bound of 72 meV that the same section quotes and plots.","tokens_in":10011,"tokens_out":5848,"duration_ms":53499,"significance":"If the normal-hierarchy predictions survive scrutiny, the model is a complete and testable framework: it unifies a radiative seesaw with a modular flavor symmetry, treats the Dirac CP phase and mass sum as predictions rather than inputs, and makes falsifiable statements for upcoming experiments on neutrinoless double beta decay, dark matter, and charged lepton flavor violation. The paper also applies current LFV and DM constraints rather than ignoring them. However, the inverted-hierarchy scenario is internally inconsistent: because the model has only two right-handed neutrinos, the lightest neutrino mass is zero, and in the inverted hierarchy the predicted mass sum is about 100 meV, which is above the DESI+CMB bound of 72 meV that the authors themselves adopt. This discrepancy is not acknowledged, so the IH sections as written overstate the viable parameter space. The NH results are not affected by this issue, but their numerical values depend on explicit modular form expressions that are not displayed in the paper, which is a separate load-bearing gap.","major_comments":[{"comment":"The inverted-hierarchy scan points are excluded by the very DESI+CMB bound the paper quotes. The model has only two right-handed neutrinos, and the text states that the lightest neutrino mass eigenvalue is zero (Sec. II, after Eq. (II.11)); in the inverted hierarchy this gives m3 = 0, so the sum of neutrino masses is m1 + m2, which is about 2*sqrt(Delta m^2_atm) ≈ 100 meV. The paper quotes the DESI+CMB bound sum m_nu <= 72 meV in Sec. II and draws it as a vertical line in Figs. 8 and 11, yet the IH points in those figures are localized near 100 meV, to the right of the line, and the text reports these values without noting the exclusion. This is a load-bearing internal inconsistency in the claimed IH predictions. The IH sections should either be withdrawn or explicitly labeled as excluded by the quoted bound, and the abstract and conclusions should be adjusted accordingly.","section":"Sec. IV.B, Figs. 8 and 11"},{"comment":"The numerical results depend entirely on the explicit expressions of the non-holomorphic modular forms Y_2^(2), Y_2^(-2), and Y_2^(0) for S3, which are taken from ref. [22] but are not displayed or tested in this paper. Different normalizations or phase conventions of these modular forms could shift the fitted regions and all the predictions in Sec. IV. The authors should provide the explicit expressions (in the text or an appendix) and state the convention they use, so that the scan is reproducible and the basis matches that of ref. [22].","section":"Sec. II, Eqs. (II.3), (II.5), (II.7)"},{"comment":"The chi-square analysis is not fully specified. The text does not define the chi-square function, the covariance matrix used for the NuFit 6.0 data, the number of degrees of freedom, or the criterion that defines the colored sigma intervals in Figs. 1-12. Without this information, the quality of the fit cannot be assessed and the colored regions are not reproducible. Please provide the explicit chi-square definition and the mapping from chi-square values to the sigma labels.","section":"Sec. IV"}],"minor_comments":[{"comment":"In the paragraph following Eq. (II.1), the phrase 'is given by is given by' is duplicated; please remove the repetition.","section":"Sec. II"},{"comment":"The last term of the Dirac Yukawa Lagrangian is written with beta_l, but Eq. (II.7) and the text indicate that this coefficient should be beta_nu (or beta_tilde_nu) to match the definition in the Yukawa matrix.","section":"Eq. (II.6)"},{"comment":"The notation P Dnu for the sum of neutrino masses is nonstandard and easy to misread; please use sum m_nu and define it as sum_i D_nu_i at first occurrence.","section":"Sec. II after Eq. (II.11)"},{"comment":"The phrase 'v1/2 rel ≈ 0.3' appears to be a typo; the standard notation would be v_rel ≈ 0.3 for the relative velocity in the p-wave cross-section expansion.","section":"Sec. III.A"},{"comment":"The summation in the LFV branching ratio formula is written with an 'X' placeholder; please write out the sum over alpha = 1-3 explicitly.","section":"Eq. (II.8)"},{"comment":"The sentence 'We achieve chi-square analysis and demonstrate some predictions' is grammatically awkward; please rephrase to something like 'We perform a chi-square analysis and demonstrate predictions'.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the inverted-hierarchy inconsistency: the IH predictions violate the very cosmological bound the paper quotes, so the IH sections should be withdrawn or explicitly labeled as excluded. I would also ask the authors to provide the explicit modular form expressions, since all numerical results hinge on a convention taken from a reference that is not reproduced. Once these points are addressed, the NH part appears solid and the paper could be publishable. The paper would also benefit from a clearer statement of the chi-square methodology, as noted in the major comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the NH analysis is a usable, concrete application of non-holomorphic modular S3 to the Ma radiative seesaw; the IH analysis, as printed, is excluded by the very cosmological bound the paper itself quotes. That inconsistency is real and needs to be fixed before the IH sections can be taken seriously.\n\nWhat the paper does well: it takes the recent non-holomorphic modular S3 framework from Qu–Ding and applies it to a well-known radiative seesaw, performing a genuine chi-square scan against NuFit 6.0. The NH fermionic-DM case gives a compact set of predictions: δCP roughly 140–240°, m_ee between 2 and 4.4 meV, sum of neutrino masses 57–60 meV, and a DM mass in the 10–3600 GeV range, with μ→eγ able to touch the current upper bound. Those are concrete, falsifiable statements, and the paper is honest about treating δCP as a prediction rather than an input. The LFV and DM constraints are imposed, not hand-waved.\n\nThe soft spot that matters: with two right-handed neutrinos, the lightest active neutrino mass is exactly zero, as the paper states. In the inverted hierarchy this gives m3 = 0 and therefore Σmν = m1 + m2 ≈ 100 meV, which is well above the DESI+CMB bound Σmν ≤ 72 meV that the authors quote and draw in Figs. 8 and 11. Yet the IH scan points sit right at that 100 meV region, and the text presents them as viable predictions. This is not a matter of modular-form conventions or fitting details; it follows directly from the model's particle content and the paper's own constraint. The IH sections either need to be withdrawn or explicitly flagged as excluded. The NH predictions remain unaffected, so the central NH claim survives.\n\nSmaller issues: the explicit modular forms Y^(2)_2, Y^(-2)_2, and Y^(0)_2 are only referenced from Qu–Ding, not displayed, which makes the scan harder to reproduce. The chi-square definition, the scan size, and how the sigma-level color bands are assigned are also missing. For a model paper this is not fatal, but it is the kind of detail a referee should ask for. The DM treatment is approximate in places (e.g., assuming degenerate inert scalars for FDM, and simply rescaling a known result for BDM), but adequate for the paper's purpose.\n\nNovelty is modest: this is a direct extension of the same authors' earlier holomorphic S3 paper to the non-holomorphic version, with new numerical predictions. That is enough to make it citable within the modular-flavor community, and it deserves a serious referee rather than a desk rejection. The referee should insist on fixing the IH presentation and on providing enough detail to reproduce the NH scan.","headline":"A serviceable modular-flavor application whose NH predictions are worth a look, but the IH sections are internally inconsistent with the paper's own DESI+CMB bound.","tokens_in":10605,"tokens_out":2730,"would_cite":false,"duration_ms":27660,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A radiative seesaw driven by a non-holomorphic modular S3 flavor symmetry can fit all measured neutrino observables and, for normal mass ordering with fermionic dark matter, predicts narrow windows for the Dirac CP phase, the neutrinoless…","keywords":["non-holomorphic modular symmetry","S3 flavor symmetry","radiative seesaw","Ma model","neutrino mass","dark matter","lepton flavor violation","modular forms"],"falsifier":"A cosmological measurement of the neutrino mass sum outside the $57$–$60$ meV window, under the normal mass ordering, would falsify the model's normal-hierarchy predictions.","tokens_in":9478,"feed_emoji":"⚛️","tokens_out":13642,"duration_ms":108324,"temperature":0.7,"pith_summary":"The paper constructs a model in which neutrino masses arise at one loop via the Ma radiative seesaw, with a non-holomorphic modular $S_3$ flavor symmetry fixing the lepton Yukawa structure. The authors perform a chi-square fit to the measured neutrino mixing angles and mass-squared differences from NuFit 6.0 and show that the model can accommodate all five observables. In the normal mass ordering with a fermionic dark-matter candidate, the fit becomes predictive: the Dirac CP phase $\\delta_{\\mathrm{CP}}$ is confined to about $140^\\circ$–$240^\\circ$, the neutrinoless double-$\\beta$-decay mass $m_{ee}$ to $2$–$4.4$ meV, the sum of neutrino masses to $57$–$60$ meV, and the dark-matter mass to $10$–$3600$ GeV, with $\\mathrm{BR}(\\mu\\to e\\gamma)$ able to reach the current experimental upper bound. These are concrete, testable targets for upcoming neutrino, dark-matter, and flavor experiments.","feed_headline":"140-240 degree CP phase emerges from S3 flavor fit","feed_subtitle":"The model fits all oscillation data and predicts dark matter at 10-3600 GeV, within reach of upcoming experiments.","key_machinery":"The load-bearing object is the non-holomorphic modular $S_3$ symmetry and its three modular forms $Y_2^{(2)}$, $Y_2^{(-2)}$, and $Y_2^{(0)}$ of weight 2, which respectively enter the charged-lepton, right-handed-neutrino, and Dirac-Yukawa sectors. These forms are taken from the explicit construction in ref. [22] and are not reproduced in the paper; they replace the flavon fields and fix the flavor structure through the modulus $\\tau$. Neutrino masses are generated at one loop by the Ma radiative seesaw, with the lightest right-handed neutrino and the inert doublet $\\eta$ running in the loop; the overall scale is set by $\\kappa = M_2\\gamma_\\nu^2$. The $\\chi^2$ scan varies the complex parameter $\\widetilde{M}_1$, the Yukawa ratios $\\tilde{\\alpha}_\\nu$ and $\\tilde{\\beta}_\\nu$, and the coupling $\\gamma_\\nu$, then maps the allowed regions onto predictions for $\\delta_{\\mathrm{CP}}$, $m_{ee}$, $\\sum m_\\nu$, the dark-matter mass, and the flavor-violating branching ratios.","core_discovery":"The central claim is that a single non-holomorphic modular $S_3$ symmetry, applied to the Ma radiative-seesaw particle content, is enough to reproduce the observed lepton mixing and mass splittings while also providing viable dark matter and respecting lepton-flavor-violation bounds. The charged-lepton mass matrix, the right-handed-neutrino mass matrix, and the Dirac Yukawa coupling are all built from three weight-2 modular forms $Y_2^{(2)}$, $Y_2^{(-2)}$, and $Y_2^{(0)}$ together with a small set of free parameters. A scan over these parameters with a $\\chi^2$ statistic against the five oscillation observables yields four branches: normal or inverted ordering, each with a fermionic or bosonic dark-matter candidate. In the normal-ordering fermionic branch the fit localizes $\\delta_{\\mathrm{CP}}$ between about $140^\\circ$ and $240^\\circ$, $m_{ee}$ between $2$ and $4.4$ meV, the neutrino mass sum between $57$ and $60$ meV, and the dark-matter mass between $10$ and $3600$ GeV, while allowing $\\mathrm{BR}(\\mu\\to e\\gamma)$ to reach the current upper bound. In the inverted-ordering branches the mass sum is predicted near $100$ meV, which exceeds the DESI+CMB bound of $72$ meV that the paper itself cites, so those branches are in tension with that cosmological constraint.","pith_inferences":["Our reading of the quoted ranges: the inverted-ordering branches predict $\\sum m_\\nu \\sim 100$ meV, which would conflict with the DESI+CMB bound of $72$ meV; if that bound stands, the model effectively selects the normal ordering.","Because the predictions depend on the adopted modular-form expressions, an independent derivation or numerical evaluation of the $S_3$ non-holomorphic modular forms would provide a cross-check and might shift the allowed windows.","The model forces the lightest neutrino mass to zero by using only two right-handed neutrinos; a three-right-handed-neutrino extension could test whether the predictions persist when all three active masses are nonzero.","The sharply predicted $\\delta_{\\mathrm{CP}}$ window could serve as a handle for leptogenesis scenarios, since a non-zero CP phase is a necessary ingredient; the paper does not pursue this connection."],"forward_implications":["If the normal-ordering fermionic branch is realized, the Dirac CP phase lies in $140^\\circ$–$240^\\circ$, making it directly testable in long-baseline neutrino oscillation experiments.","The predicted $m_{ee}$ of $2$–$4.4$ meV is within the projected reach of next-generation neutrinoless double-beta-decay experiments such as KamLAND-Zen.","A neutrino mass sum of $57$–$60$ meV is a narrow, cosmologically testable band that upcoming DESI/CMB measurements can confirm or exclude.","The dark-matter mass window ($10$–$3600$ GeV for fermionic DM in the normal ordering, and about $534$ GeV for bosonic DM) gives a concrete target for direct and indirect dark-matter searches.","The $\\mu\\to e\\gamma$ branching ratio can approach the current experimental upper bound in the fermionic normal-ordering branch, so MEG II could observe the signal."],"supporting_citations":[{"why":"supplies the Ma radiative-seesaw field content and the one-loop neutrino mass mechanism","marker":"[10]"},{"why":"provides the explicit non-holomorphic modular forms for S3 that fix the flavor structure","marker":"[22]"},{"why":"the NuFit 6.0 data set for the oscillation observables used in the chi-square fit","marker":"[28]"},{"why":"gives the experimental upper bounds on mu->e gamma, tau->e gamma, tau->mu gamma that constrain the model","marker":"[31]"},{"why":"provides the DESI+CMB bound on the neutrino mass sum used to constrain the predictions","marker":"[32]"},{"why":"supplies the measured dark-matter relic density that the DM candidates must reproduce","marker":"[36]"},{"why":"computes the inert-scalar relic density that fixes the bosonic DM mass near 534 GeV","marker":"[38]"},{"why":"derives the lepton-flavor-violation formulas used for the branching-ratio predictions","marker":"[29]"}],"fun_headline_variants":["Radiative seesaw from S3 modular symmetry fits all data","Modular S3 predicts CP phase 140-240° and DM mass range","S3 flavor model gives neutrino mass and dark matter at loop level","Inverted ordering tensioned in S3 neutrino dark matter model","S3 symmetry yields viable neutrino sector and dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions stand or fall on the explicit expressions and normalization conventions of the three non-holomorphic modular forms taken from ref. [22]; if those expressions differ, the fitted regions and all quoted ranges would shift.","fun_headline_variants_meta":{"raw":{"variants":["Radiative seesaw from S3 modular symmetry fits all data","Modular S3 predicts CP phase 140-240° and DM mass range","S3 flavor model gives neutrino mass and dark matter at loop level","Inverted ordering tensioned in S3 neutrino dark matter model","S3 symmetry yields viable neutrino sector and dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000648,"raw_usage":{"total_tokens":2957,"prompt_tokens":912,"completion_tokens":2045,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":1955}},"tokens_in":528,"tokens_out":2045,"duration_ms":12627,"temperature":1.0,"reasoning_tokens":1955,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:58:15.101558+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A cosmological measurement of the neutrino mass sum outside the $57$–$60$ meV window, under the normal mass ordering, would falsify the model's normal-hierarchy predictions.","supporting_citations":[{"cited_title":"Neutrino phenomenology in the modular $S_3$ seesaw model","cited_arxiv_id":"2403.00593","evidence_quote":"provides the explicit non-holomorphic modular forms for S3 that fix the flavor structure"}],"review_version":1}