{"id":"7160988e-ea82-4d92-b9b8-4fae50ae1484","arxiv_id":"1908.08409","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A modular S4 radiative seesaw model fits normal-hierarchy neutrino data and predicts a total neutrino mass of roughly 58-62 meV and a double-beta decay mass of 1-4 meV.","lead":"This paper builds a neutrino mass and dark matter model using a modular S4 symmetry, with neutrino mass generated in a one-loop radiative seesaw. It shows the model can match normal-hierarchy neutrino oscillation data and predicts a narrow range for the total neutrino mass and neutrinoless double beta decay.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"DM mass window 534±8.5 GeV is taken from a kinetic-term-only relic calculation; the model's extra Yukawa and quartic couplings could shift it, and the numerical scan depends on it.","rationale":"The reader's weakest_assumption identifies exactly the gap I consider most load-bearing: the input DM mass comes from an external calculation that does not include this model's interactions. I agree with that identification. The model's central numerical output is a parameter scan whose scalar mass spectrum is anchored to the quoted DM window; if the relic-density constraint moves m_I, the allowed region and derived predictions shift. This is not a disagreement with the modular S4 framework or an ad hominem; it is a missing computation in the argument as written. A single independent relic-density calculation would settle whether the concern lands: if the full model still requires m_I near 534 GeV, the numerical results stand; if not, the scan must be redone. The CONDITIONAL verdict is appropriate because the construction is coherent and the neutrino sector can be checked, but the paper should not be accepted as fully predictive until the DM input is justified within the model or the scan is shown to be insensitive to it.","tokens_in":10690,"tokens_out":22228,"duration_ms":209837,"concrete_test":"Use micrOMEGAs or MadDM to compute the thermal relic density of the full scalar sector (η_I, η_R, η±) for the two Table III sample points and for a grid of m_I in [500,570] GeV, including the Yukawa matrix yη from Eq. (II.9) and the scalar potential from Eq. (II.10) with λHη and λ'Hη scanned over [-1,1]. Compare the m_I value that reproduces Ωh² = 0.12 with the assumed 534±8.5 GeV input. If the required mass shifts by more than ~10 GeV, the scan in Section II.A must be redone and the reported prediction ranges are not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section II.A the scan restricts m_R to 525.5-542.5 GeV on the basis of the quoted relic-density window 534±8.5 GeV from Ref. [48]. That reference computes the thermal abundance for an inert scalar doublet whose interactions are only the gauge/kinetic terms. The present model has additional Yukawa couplings yη in Eq. (II.9) and the scalar quartic couplings in Eq. (II.10); the latter control (co)annihilation of η_I with η_R and η±, and the former add t-channel lepton/neutrino final states. The quartic couplings are not fixed by the paper and could be O(1). The paper neither computes Ωh² in this model nor shows that the [48] window remains valid when λHη, λ'Hη, and yη are included. Since m_R is an input to the neutrino loop function in Eq. (II.18) and to the scan, an O(1) shift in the required DM mass would force a rerun of the χ² scan. The quoted 58-62 meV and ⟨mee⟩ ranges are therefore contingent on an unverified relic-density assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a lepton model with a modular S4 symmetry in which neutrino masses arise from a radiative seesaw at one-loop order. An inert scalar doublet eta and heavy Majorana fermions are introduced; well-assigned modular weights forbid a tree-level neutrino mass and leave a remnant Z2 symmetry that stabilizes the lightest inert scalar component as dark matter. The authors specify the charged-lepton, Dirac Yukawa, and heavy-neutrino mass matrices, derive the one-loop neutrino mass formula, and perform a chi-square scan of the parameter space. They report that the normal hierarchy of neutrino masses is reproduced, that the inverted hierarchy is disfavored, and that the model predicts the sum of neutrino masses (58-62 meV), the neutrinoless double beta decay mass (1.2-4 meV), and a Dirac CP phase localized near pi/2, while satisfying charged-lepton flavor violation bounds.","tokens_in":11089,"tokens_out":4630,"duration_ms":52730,"significance":"If the numerical results are robust, the model would provide an interesting connection between modular flavor symmetry, radiative neutrino mass, and dark matter stability, with concrete testable predictions for the neutrino mass sum, 0nu beta beta decay, and the Dirac CP phase. The manuscript is transparent in specifying the field content, the Lagrangian, the mass matrices, and the chi-square procedure, and it explicitly checks lepton flavor violation constraints. Its main limitations are that the dark matter mass range is imported from a simplified relic-density calculation, and the solar mass-squared difference is used as an input to fix the quartic coupling rather than being a prediction, which reduces the claimed predictive content.","major_comments":[{"comment":"The scan restricts m_R to 525.5-542.5 GeV using the dark matter mass window 534 +/- 8.5 GeV from Ref. [48]. That reference computes the thermal relic density for an inert scalar doublet whose interactions consist only of the gauge and kinetic terms. The present model contains additional couplings, notably the Yukawa matrix y_eta in Eq. (II.9) and the quartic scalar couplings lambda_Heta, lambda'_Heta, lambda''_Heta in Eq. (II.10), which provide extra annihilation and coannihilation channels for the DM candidate eta_I with eta_R, eta+- and lepton final states. The paper neither computes Omega h^2 in this full model nor demonstrates that the quoted window remains valid once these couplings are switched on. Since m_R enters the neutrino mass loop function in Eq. (II.18) and the chi-square scan, a shift in the required DM mass would force a re-evaluation of the numerical predictions. This issue is load-bearing for the central claim, and should be addressed either by computing the relic density in the model or by treating m_R as an independent parameter and showing the quoted ranges for sum m_i and <m_ee> are stable.","section":"II.A"},{"comment":"The quartic coupling lambda''_Heta is fixed by Eq. (II.21) to reproduce the solar mass-squared difference Delta m^2_sol. Consequently, the solar splitting is an input to the fit, not a prediction of the model. The chi-square analysis in Sec. II.A includes Delta m^2_sol as one of the five observables, so the effective number of independently predicted observables is smaller than the nominal five. The abstract and conclusion should be qualified to state clearly which quantities are predictions (mixing angles, phases, mass ratios, and the correlated ranges for sum m_i and <m_ee>) and which are used as inputs (the charged-lepton masses and the solar splitting).","section":"II, Eq. (II.21)"},{"comment":"There is an inconsistency in the sign of the mass-splitting relation. From Eqs. (II.12) and (II.13), m_R^2 - m_I^2 = lambda''_Heta v_H^2. With a positive lambda''_Heta, the real scalar is heavier than the imaginary one, which is the correct hierarchy for eta_I to be the dark matter candidate. However, the text in Sec. II.A states m_I = sqrt(m_R^2 + 2 lambda''_Heta v_H^2), which implies m_I^2 - m_R^2 = 2 lambda''_Heta v_H^2 and would make m_I heavier for positive lambda''. This sign error affects the expression used in the numerical scan and, through the sign of the loop function in Eq. (II.18), the neutrino mass matrix. Please correct the formula and verify that the numerical results in Table III and Figs. 1-3 correspond to the correct sign convention.","section":"II, Eqs. (II.12)-(II.13) and Sec. II.A"}],"minor_comments":[{"comment":"In the list of input parameters, 'alpha_D' appears twice; presumably one is 'alpha_eta' or another Yukawa coupling. Please correct this typo.","section":"Sec. II.A"},{"comment":"The phrase 'Mojorana phases' should be 'Majorana phases'.","section":"Fig. 3 caption and Sec. III"},{"comment":"The caption reads 'at the sample points of NH and IH', but the text and the values indicate that both sample points are for the normal hierarchy. Since the paper concludes that IH is disfavored, the caption should be changed to 'NH'.","section":"Table III caption"},{"comment":"The probability density function is written as f(x,nu) = x^{nu/2 - 1} exp(x/2) / (2^{nu/2} Gamma(x/2)); the Gamma function argument should be nu/2, not x/2.","section":"Appendix"},{"comment":"The scan region for tau is Re(tau) in [-1.5, 1.5], which extends beyond the standard fundamental domain of the modular group. Since modular transformations identify points outside the fundamental domain, the authors should either restrict tau to the fundamental domain or explain why the larger scan region is justified.","section":"Sec. II.A and Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a typical modular flavor model paper with a radiative seesaw, and the central mechanism is clearly presented. The main concern is not the modular S4 construction itself but the unverified import of the DM mass window from a kinetic-term-only relic calculation. If the authors were to compute the relic density in the full model or otherwise justify the mass window, the paper could become acceptable. The sign inconsistency in the m_I formula is also concerning and should be corrected before any further consideration."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a modular S4 version of the Ma one-loop radiative seesaw. The genuinely new piece is the specific embedding — three Majorana families as singlet plus doublet, charged leptons as one singlet plus one doublet, left-handed doublets as triplet — and the way modular weights leave a remnant Z2 that stabilizes the inert doublet. That combination is not in the earlier modular S4 papers [18-20], so it is a real if modest addition to the modular flavor program.\n\nWhat is well done: the mass matrices are written explicitly, the χ² scan over the modulus and the couplings is transparent, two sample points are given, and the LFV bounds are checked. The NH fit is coherent and yields concrete, narrow ranges: sum m_i in 58–62 meV, <m_ee> in 1.2–4 meV, and Dirac CP near π/2. The IH exclusion is stated as a no-solution; reasonable, but it would be stronger if the scan details were published.\n\nThe main soft spot is the DM mass input. Section II.A takes m_I ≈ 534 ± 8.5 GeV from Ref. [48], which computed the relic density of an inert doublet with essentially kinetic-term-only interactions. The present model adds the Yukawa y_η and the quartic couplings of Eq. (II.10); these affect DM annihilation and coannihilation, and none of them is fixed by the paper. The authors are explicit about the assumption, but they do not compute Ωh² in this model or show the window survives. Because m_R is scanned in the narrow 525.5–542.5 GeV range and feeds the neutrino loop, a change in the relic-density-required mass would require a rerun of the whole scan. So the quoted neutrino predictions are contingent, not robust.\n\nA related but lesser point: λ''_Hη is fixed by Eq. (II.21) to reproduce Δm²_sol, and the charged-lepton couplings are fitted to the masses. Those are inputs, not outputs. The genuinely predicted quantities are the mixing angles, the CP phases, and the correlations in Figs. 2–3, which is still worth having but is less than the abstract's 'predictive' suggests.\n\nNo catastrophic misunderstanding at the core: the model is consistent on its own terms, and the authors flag the fragile assumptions. It is one of many modular flavor models, so the audience is specialists who will use it as a benchmark example. I would send it to a serious referee, asking for an in-model relic density calculation or a convincing argument that the [48] window holds, and for the scan details. With that revision it becomes a solid model-building paper; without it, the DM and mass predictions are on borrowed footing.","headline":"A clean but routine modular S4 Ma-model extension that fits NH data; the DM mass is borrowed, not computed, so the predictions are less robust than the abstract suggests.","tokens_in":11593,"tokens_out":3743,"would_cite":false,"duration_ms":58812,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq","95.35.+d"],"model":"deepseek-v4-flash","headline":"A modular $S_4$ symmetry whose weights forbid a tree-level neutrino mass generates the masses radiatively and predicts a normal hierarchy, $\\sum m_i$ near 58–62 meV, $\\langle m_{ee}\\rangle$ 1.2–4 meV, and a $Z_2$-stabilized dark matter…","keywords":["modular S4 symmetry","radiative seesaw","neutrino masses","normal hierarchy","dark matter stability","remnant Z2 symmetry","neutrinoless double beta decay","Dirac CP phase"],"falsifier":"Recompute the thermal relic density of the full model with the new Yukawa couplings $y_\\eta$ and the quartic scalar couplings included; if the dark matter mass required to match the observed abundance moves outside the assumed $534\\pm8.5$ GeV window, the input scalar masses used for the neutrino fit are invalid and the predicted $\\sum m_i$ and $\\langle m_{ee}\\rangle$ windows shift. Separately, a measurement of $\\sum m_i$ outside 58–62 meV, or of $\\langle m_{ee}\\rangle$ outside 1.2–4 meV, at $5\\sigma$ would rule out the normal-hierarchy solution.","tokens_in":10487,"feed_emoji":"⚛️","tokens_out":20373,"duration_ms":329543,"temperature":0.7,"pith_summary":"This paper proposes that a single modular $S_4$ flavor symmetry, through carefully chosen modular weights, can forbid a tree-level neutrino mass and at the same time leave a remnant $Z_2$ parity that stabilizes a dark matter candidate. With the lepton fields placed in specific $S_4$ representations, the neutrino mass matrix arises from a one-loop radiative seesaw and fits the observed normal hierarchy. The payoff is a compact set of predictions: $\\sum m_i$ in a narrow window around 58–62 meV, $\\langle m_{ee}\\rangle$ between 1.2 and 4 meV, and a Dirac CP phase that clusters near $\\pi/2$. If these survive experimental scrutiny, the model would connect the origin of neutrino flavor, the dark matter mass, and an unbroken $Z_2$ stability in one parameter $\\tau$.","feed_headline":"58–62 meV: modular symmetry fixes the neutrino mass sum","feed_subtitle":"The model also predicts a 534 GeV dark matter particle and a testable double-beta decay signal.","key_machinery":"The load-bearing mechanism is the modular $S_4$ symmetry with fixed modular weights, integers that control how fields transform under the modular group and which couplings are allowed. The lowest-weight modular forms $Y_2^{(2)}$, expressed through the eta function and its derivative, generate the higher-weight triplets $Y_3^{(4)}$ and $Y_{3'}^{(4)}$ that enter the Yukawa matrices. These matrices feed the one-loop radiative seesaw formula for the neutrino mass, where the inert scalar doublet $\\eta$ runs in the loop; the same assignments that forbid the tree-level mass produce a remnant $Z_2$ parity that stabilizes the lightest $\\eta$ component as dark matter.","core_discovery":"The central discovery is that the modular weights themselves do the structural work: assigning weight $-2$ to the lepton fields and weight $-1$ to the Majorana fields makes the tree-level neutrino mass matrix vanish, so the leading contribution is the one-loop diagram with the inert doublet $\\eta$. After modular symmetry breaking a $Z_2$ subgroup survives, odd parity for odd-weight fields and even for even-weight fields, and the imaginary component of $\\eta^0$ is the stable dark matter. The Yukawa sector is then structured by modular forms built from the eta function and its derivative, and a $\\chi^2$ fit to current neutrino oscillation data selects a narrow region of the modular parameter $\\tau$ near $-0.5+1.35i$. In that region the model reproduces the normal hierarchy and predicts $58\\,\\text{meV}\\lesssim\\sum m_i\\lesssim62\\,\\text{meV}$, $1.2\\,\\text{meV}\\lesssim\\langle m_{ee}\\rangle\\lesssim4\\,\\text{meV}$, and a Dirac phase that clusters near $\\pi/2$; it finds no solution for the inverted hierarchy. Two sample points are exhibited, one minimizing $\\Delta\\chi^2$ and one reproducing the best-fit $\\delta_{CP}=195^\\circ$.","pith_inferences":["Inference: recomputing the relic density with the full coupling set would likely move the allowed dark matter mass, and because $m_R$ and $m_I$ enter the one-loop neutrino mass formula, the predicted $\\sum m_i$ and $\\langle m_{ee}\\rangle$ windows could shift accordingly.","Inference: the same modular parameter $\\tau$ that sets the flavor structure also determines the scalar mass spectrum through the potential, so a search for the charged inert scalar $\\eta^\\pm$ could test the neutrino and dark matter regions in a way the paper does not spell out.","Inference: the allowed region sits near $\\tau\\approx -0.5+1.35i$ rather than at the fixed point $\\tau=i$, suggesting the predictions are not a pure symmetry-forced accident; variants that impose exact fixed-point values would give different correlations between $\\delta_{CP}$ and the Majorana phases."],"forward_implications":["The model selects the normal hierarchy: no parameter choices were found that fit the inverted hierarchy.","The sum of neutrino masses is confined to 58–62 meV, within reach of upcoming cosmological and neutrino-mass experiments.","The neutrinoless double beta decay effective mass sits between 1.2 and 4 meV, below current limits but potentially observable by next-generation experiments.","The Dirac CP phase clusters near $\\pi/2$ in the $5\\sigma$ allowed region, while the Majorana phases show a linear correlation rather than being free.","Lepton flavor violating branching ratios stay far below current bounds, so the model introduces no observable flavor problem."],"supporting_citations":[{"why":"supplies the one-loop radiative seesaw structure with an inert scalar doublet and heavy Majorana fermions that the model builds on.","marker":"[1]"},{"why":"establishes the modular flavor symmetry formalism that the paper applies to the group $S_4$.","marker":"[2, 3]"},{"why":"provides the explicit $S_4$ modular forms, written through the eta function and its derivative, used to build the Yukawa couplings.","marker":"[39]"},{"why":"gives the lepton-flavor-violating branching ratio formulas used to constrain the model.","marker":"[40, 41]"},{"why":"sets the experimental upper bounds on $\\mu\\to e\\gamma$, $\\tau\\to e\\gamma$, and $\\tau\\to\\mu\\gamma$ imposed in the scan.","marker":"[42–44]"},{"why":"supplies the cosmological upper bound on the sum of neutrino masses that the predicted $\\sum m_i$ must satisfy.","marker":"[45, 46]"},{"why":"defines the future neutrinoless double beta decay sensitivity against which the predicted $\\langle m_{ee}\\rangle$ is compared.","marker":"[47]"},{"why":"is the relic-density calculation whose dark matter mass window $534\\pm8.5$ GeV is adopted as input.","marker":"[48]"},{"why":"provides the global fit to neutrino oscillation data used for the $\\chi^2$ analysis and confidence levels.","marker":"[49]"}],"fun_headline_variants":["Modular S4 symmetry predicts 58–62 meV neutrino mass sum","One-loop radiative seesaw from modular S4 yields dark matter","Modular weights kill tree-level mass, one-loop fixes neutrinos","Modular S4 model: normal hierarchy, 534 GeV dark matter","Radiative seesaw with modular S4 predicts 58–62 meV sum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dark matter mass is imported from a relic-density calculation that only includes the scalar kinetic term, even though the model's new Yukawa and quartic couplings would change the annihilation rate; if the full calculation moves the required mass, all the numerical predictions built on it shift.","fun_headline_variants_meta":{"raw":{"variants":["Modular S4 symmetry predicts 58–62 meV neutrino mass sum","One-loop radiative seesaw from modular S4 yields dark matter","Modular weights kill tree-level mass, one-loop fixes neutrinos","Modular S4 model: normal hierarchy, 534 GeV dark matter","Radiative seesaw with modular S4 predicts 58–62 meV sum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1553,"prompt_tokens":971,"completion_tokens":582,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":485}},"tokens_in":587,"tokens_out":582,"duration_ms":5739,"temperature":1.0,"reasoning_tokens":485,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:58:24.252534+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the thermal relic density of the full model with the new Yukawa couplings $y_\\eta$ and the quartic scalar couplings included; if the dark matter mass required to match the observed abundance moves outside the assumed $534\\pm8.5$ GeV window, the input scalar masses used for the neutrino fit are invalid and the predicted $\\sum m_i$ and $\\langle m_{ee}\\rangle$ windows shift. Separately, a measurement of $\\sum m_i$ outside 58–62 meV, or of $\\langle m_{ee}\\rangle$ outside 1.2–4 meV, at $5\\sigma$ would rule out the normal-hierarchy solution.","supporting_citations":[{"cited_title":"The ranges of DN are respectively 40 TeV 9 ≤DN1≤ 80 TeV, 130 TeV≤DN2≤ 420 TeV, and 160 TeV≤DN3≤ 480 TeV","cited_arxiv_id":null,"evidence_quote":"supplies the one-loop radiative seesaw structure with an inert scalar doublet and heavy Majorana fermions that the model builds on."},{"cited_title":"Ishimori, T","cited_arxiv_id":null,"evidence_quote":"provides the explicit $S_4$ modular forms, written through the eta function and its derivative, used to build the Yukawa couplings."},{"cited_title":"Simple A4 models for dark matter stability with texture zeros","cited_arxiv_id":"1811.10619","evidence_quote":"is the relic-density calculation whose dark matter mass window $534\\pm8.5$ GeV is adopted as input."}],"review_version":1}