{"id":"b7ef4917-e01b-47ad-98df-400bbb441e13","arxiv_id":"2505.03000","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An A4 modular inverse seesaw model with a U(1) B-L Z' fits neutrino data and produces the observed baryon asymmetry via resonant leptogenesis.","lead":"This paper builds a model of neutrino masses and mixing using A4 modular symmetry with an inverse seesaw mechanism and a Z' boson, and shows the same parameter points can generate the observed matter-antimatter asymmetry through resonant leptogenesis. The new ingredient is the Z' contribution to the leptogenesis equations, and the model yields a tight, testable correlation between two neutrino mixing angles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed sin^2(theta23) lower bound and sin^2(theta12)-sin^2(theta23) linear relation are scan-box artifacts until a dramatically wider scan shows they persist.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the scan ranges in Eq. (17) are hand-picked, and the claimed exclusions and correlations may be artifacts of that box. This is the correct central soft spot because the paper's distinctive claims are predictive statements about allowed mixing angles, not just existence of a fit. The model has several continuous parameters with no symmetry-enforced ranges; any random-scan result must be shown stable under enlarging the sampled region before it can be called a prediction of the A4 modular framework. The proposed test directly checks stability by broadening the scan over all free parameters while keeping the same oscillation acceptance criteria. A secondary concern noted by the reader is the typo in the Boltzmann equation (the missing derivative in Eq. (41) and the repeated Y_chi1 index in the j=2 term), but that is a presentation/numerical-implementation issue that can be fixed by a corrected integration, whereas the scan-prior problem affects the interpretation of the main phenomenological results. The verdict should remain CONDITIONAL: the paper's consistency analysis is plausible, but the predictive claims require the wider scan test before acceptance.","tokens_in":18765,"tokens_out":18969,"duration_ms":222288,"concrete_test":"Repeat the analysis with a substantially wider scan: alpha_p in [1e-6, 1e-1], beta_p in [1e-5, 1], mu_0 in [1e-3, 1e3] GeV, v_phi in [1e3, 1e9] GeV, and tau over the full fundamental domain (Im[tau] > 0, |tau| >= 1, Re[tau] in [-1/2, 1/2]), keeping the same 3-sigma acceptance criteria from Table III. Record the minimum sin^2(theta23) among accepted points and the slope and width of the sin^2(theta12)-sin^2(theta23) band. If the lower bound disappears or the band widens substantially, the paper's predictions are scan artifacts; if both survive, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.A and Eq. (17) restrict the scan to Re[tau] in [-0.5, 0.5], Im[tau] in [0.5, 1.5], alpha_p in [1e-4, 1e-3], beta_p in [1e-3, 1e-2], mu_0 in [0.1, 10] GeV, and v_phi in [1e5, 1e8] GeV. Nothing in the A4 modular structure or the U(1)_{B-L} charge assignment forces these parameters into those intervals: alpha_p and beta_p are unconstrained Yukawa coefficients, and mu_0 and v_phi are dimensionful parameters whose natural ranges are not fixed by the symmetry. The reported exclusion of sin^2(theta23) < 0.44 and the narrow linear band in Fig. 3 are therefore statements about the sampled box, not about the model. The paper does not report scan size, sampling density, or convergence, so one cannot distinguish a robust modular-texture prediction from a selection effect. This is load-bearing because the abstract's headline claims are precisely these predictions; the mere existence of some parameter point fitting oscillation data is trivial for a model with several free parameters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a minimal inverse seesaw model with A4 modular symmetry and a local U(1)_{B-L} symmetry. Neutrino masses and mixing are generated through a 7x7 neutral fermion mass matrix depending on the modulus tau and four continuous parameters (alpha_p, beta_p, mu_0, v_phi). The authors scan these parameters, retain points consistent with the 3 sigma neutrino oscillation data, and report a lower bound sin^2(theta23) >= 0.44 and a linear correlation between sin^2(theta12) and sin^2(theta23). They then check charged lepton flavor violating decays, neutrinoless double beta decay, and compute the baryon asymmetry via resonant leptogenesis including Z'-mediated lepton-number-conserving scatterings, obtaining Y_B ~ 8.6e-11.","tokens_in":19111,"tokens_out":6101,"duration_ms":66241,"significance":"If the claimed predictions are robust, the model offers a simultaneous, economical explanation of neutrino oscillation parameters, charged lepton flavor violation, and the baryon asymmetry within a modular-symmetry framework. The paper provides explicit modular forms and covers a broad set of observables, which is a strength. However, because the results are derived from a filtered random scan over a restricted parameter box, the advertised 'predictions' need stronger support to be credible; the current analysis does not yet establish that the lower bound on sin^2(theta23) and the theta12-theta23 correlation are features of the model rather than of the scan.","major_comments":[{"comment":"The claimed lower bound sin^2(theta23) > 0.44 and the linear correlation between sin^2(theta12) and sin^2(theta23) are derived from a random scan restricted to Re[tau] in [-0.5, 0.5], Im[tau] in [0.5, 1.5], alpha_p in [1e-4, 1e-3], beta_p in [1e-3, 1e-2], mu_0 in [0.1, 10] GeV, and v_phi in [1e5, 1e8] GeV. These intervals, especially for alpha_p, beta_p, mu_0, and v_phi, are not fixed by the A4 modular structure or the charge assignments; they are chosen for numerical convenience. The paper does not report the number of scan points, the sampling density, or any convergence test. Consequently, the reported exclusion of sin^2(theta23) < 0.44 and the narrow band in Fig. 3 could be artifacts of the chosen box. Since these are the headline claims of the abstract and conclusion, the authors must demonstrate robustness by substantially widening the scan, at least for the unconstrained continuous parameters, or provide an analytic argument that the features are inherent to the model.","section":"Sec. III.A, Eq. (17); Figs. 2 and 3"},{"comment":"The Boltzmann equation for Y_{B-L} contains a sum over j = 1, 2 with a factor (Y_chi1/Y_chi1^eq - 1) inside the sum; for the j = 2 term this factor should be (Y_chi2/Y_chi2^eq - 1). As written (or as likely implemented), the decay term of chi2 is weighted by the abundance of chi1, which would spoil the CP-asymmetry contribution from chi2 and make the computed Y_B unreliable. The authors should correct the equation and re-run the numerical integration, or clarify that the printed equation is a typographical error and that the code uses the correct abundance for each species.","section":"Sec. IV.B, Eq. (41)"},{"comment":"The paper states that Delta L = 1 and Delta L = 2 scattering processes 'can be safely neglected in our study as in our work K >> 1' with references [73,74]. In the standard leptogenesis literature, K >> 1 denotes the strong-washout regime, where inverse decays and Delta L = 1 scatterings are typically important and must be included. The reasoning as stated appears to be the opposite of the usual expectation. Please justify the neglect quantitatively, e.g., by comparing the relevant reaction densities to the decay density gamma_D for the scanned parameter points, or modify the Boltzmann equations to include these processes. This is directly relevant to the reliability of the reported baryon asymmetry.","section":"Sec. IV, text before Eq. (24)"},{"comment":"The conclusion states that the region Re[tau] in (-0.11, 0.98) is excluded by the model, but the scan in Sec. III.A only covers Re[tau] in [-0.5, 0.5]. Values of Re[tau] above 0.5 were never sampled, so the model cannot exclude them on the basis of the presented analysis. This claim should be corrected to the range actually scanned, or the scan must be extended beyond 0.5 to justify the exclusion.","section":"Sec. V, Conclusion"}],"minor_comments":[{"comment":"The left-hand side of the equation is printed as 'YB-L/dz' rather than 'dY_{B-L}/dz'; the missing differential operator should be corrected.","section":"Eq. (41)"},{"comment":"The text claims a lower limit sin^2(theta23) = 0.44, but the figure does not indicate this boundary explicitly; adding a horizontal guidance line would help the reader verify the claim.","section":"Fig. 2, bottom panel"},{"comment":"The table is labeled NuFIT 5.2 (2022), but the cited reference [55] is the 2020 JHEP paper by Esteban et al. Please update the reference to the NuFIT 5.2 publication or confirm that the numbers indeed correspond to the cited paper.","section":"Reference [55]"},{"comment":"The notation gamma_D and gamma_z' is introduced without explicit definitions in the text; the reader is directed to Appendix C, but the definitions of gamma_D and gamma_z' in terms of Y_eq and cross sections should be stated at first use for clarity.","section":"Sec. IV.B, after Eq. (43)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and addresses a timely topic. The main worries are (i) the scan-box dependence of the claimed neutrino-mixing predictions and (ii) an apparent error in the Boltzmann equation for leptogenesis. Both are fixable in revision. I recommend requiring a substantially wider scan and a corrected, re-run leptogenesis calculation before considering publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my honest read. This is a competent but incremental numerical study. The genuinely new piece is the inclusion of Z'-mediated, lepton-number-conserving scattering in the resonant leptogenesis Boltzmann equations; that is a real addition to the existing A4 modular inverse seesaw literature, including the same group's 2024 paper. The model is clearly laid out, the charged-lepton sector is diagonal by construction, and the scan consistently connects neutrino oscillation data to LFV bounds, 0νββ, and the final baryon asymmetry. That consistency is the paper's main strength. The claim that the Z' interaction delays but does not reduce the final asymmetry is physically reasonable and worth checking.\n\nThe soft spots are real. The scan ranges in Eq. (17) are not forced by the symmetry: αp and βp are unconstrained Yukawa coefficients, and μ0 and vφ are dimensionful parameters chosen by hand. The paper reports no scan size, sampling density, or convergence test. So the abstract's claims—that the model \"excludes\" sin²θ23 < 0.44 and produces a linear θ12–θ23 relation—are statements about the sampled box, not about the model. They may survive a wider scan, but the paper doesn't show it. This is load-bearing, because those correlations are the headline results.\n\nThere are also technical issues. Equation (41) has typos: the left side reads YB−L/dz, and the sum uses Yχ1/Yeqχ1 for both j=1,2 instead of Yχj/Yeqχj. The decay source term should depend on each heavy neutrino's own abundance. This is fixable but needs to be redone. The mZ' = 4 TeV value is cited to the PDG, which is not a real source for that number; they should cite a specific LHC constraint for the U(1)B−L model. The neglect of ΔL=1,2 scattering is justified with a brief K>>1 argument; that may be fine, but it deserves a bit more care given strong washout.\n\nThe citation pattern is honest: the overlap with the group's earlier work is visible and acknowledged. No code or data is shipped, but the equations are explicit enough to reproduce.\n\nWho is this for? BSM phenomenologists working on modular flavor symmetry or low-scale leptogenesis. It's a usable worked example, not a breakthrough. I'd send it to peer review, and I'd ask for a dramatically wider scan with reported sampling, a corrected Boltzmann equation, and a real collider bound for Z'. If those are fixed, the Z' effect on leptogenesis is a legitimate contribution.","headline":"Competent incremental A4 modular inverse seesaw paper; the Z' leptogenesis term is new, but the headline mixing-angle 'predictions' are scan-box artifacts until a wider scan proves otherwise.","tokens_in":19575,"tokens_out":4639,"would_cite":true,"duration_ms":55932,"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 minimal inverse seesaw model with A4 modular symmetry can fit neutrino mixing, evade charged-lepton flavor violation bounds, and generate the observed baryon asymmetry through resonant leptogenesis.","keywords":["inverse seesaw","A4 modular symmetry","neutrino mixing","lepton flavor violation","resonant leptogenesis","baryon asymmetry","Z' boson","neutrinoless double beta decay"],"falsifier":"A future measurement of sin² θ23 below 0.44 in the currently allowed 3σ range, or a dedicated scan covering the full fundamental domain of the modulus τ that finds oscillation-allowed points with sin² θ23 < 0.44, would falsify the model's lower-bound prediction. Similarly, a measured deviation from the predicted linear relation between sin² θ12 and sin² θ23 at the precision of upcoming experiments would rule out the correlation.","tokens_in":1592,"feed_emoji":"⚛️","tokens_out":2879,"duration_ms":75437,"temperature":0.7,"pith_summary":"This paper proposes a minimal inverse seesaw model with A4 modular symmetry and claims that one parameter regime can address three puzzles simultaneously: neutrino masses and mixing, the absence of charged lepton flavor violation, and the matter-antimatter asymmetry of the universe. Within the current 3-$\\sigma$ allowed region of oscillation data, the model excludes part of the parameter space for the atmospheric angle, predicting $sin^{2}$ θ23 > 0.44, and finds a clear linear decrease of $sin^{2}$ θ12 with $sin^{2}$ θ23. The same parameter points that fit oscillations satisfy the MEG and BaBar bounds on μ→eγ, τ→eγ, and τ→μγ, keep the neutrinoless double $\\beta$ decay rate below current sensitivity, and produce a baryon asymmetry near Y_B ~ 8.6×10⁻¹¹ through resonant leptogenesis of a nearly degenerate heavy neutrino pair, including Z′-mediated lepton number conserving scatterings. A sympathetic reader would care because the paper demonstrates that a single modular symmetry can tie low-energy flavor data to the cosmological baryon asymmetry.","feed_headline":"Neutrino mixing and matter asymmetry from one A4 modulus","feed_subtitle":"Model predicts sin² θ23 > 0.44 and a linear θ12–θ23 relation, all while staying below LFV limits.","key_machinery":"The central object is the 7×7 neutral fermion mass matrix in the flavor basis (ν_L, N_R^c, S), whose block structure gives the inverse seesaw light neutrino mass formula m_ν = M_D $M_R^{{-1}}$ μ (M_R^T)^{-1} M_D^T. The A4 modular symmetry, with a single complex modulus τ controlling all Yukawa couplings, provides the flavor structure without a large flavon sector; the heavy sector splits into two nearly degenerate pairs, and the lightest pair's decay in the resonance regime generates the CP asymmetry for leptogenesis. The same modulus τ sets the Dirac CP phase in the PMNS matrix, so the model connects low-energy CP violation to the CP asymmetry relevant for baryogenesis.","core_discovery":"The paper's central claim is that an inverse seesaw extension of the Standard Model based on A4 modular symmetry, with two right-handed neutrinos and two singlet fermions plus a local U(1)_{B−L} symmetry, is enough to reproduce all current neutrino oscillation observables while making three testable predictions: a lower bound of 0.44 on sin² θ23, a linear correlation between sin² θ12 and sin² θ23 in the 3σ region, and a strong link between the Dirac CP phase, the Jarlskog invariant, and the atmospheric angle. Using the same parameter space that fits oscillations, the model respects all three measured charged lepton flavor violating branching ratios, places the effective neutrinoless double beta decay mass |m_ee| below the KamLAND-Zen and projected nEXO sensitivities, and yields the observed baryon asymmetry through resonant leptogenesis from the decay of the lightest nearly degenerate heavy neutrino pair, with the Z′-mediated lepton number conserving scatterings included but not spoiling the final asymmetry.","pith_inferences":["If the linear θ12–θ23 relation survives a scan over the full fundamental domain of τ rather than the chosen parameter box, it would provide a rare cross-check between two independent oscillation observables that JUNO and DUNE could test at the sub-percent level.","The claimed insensitivity of the final baryon asymmetry to g_{B−L} suggests a general lesson: adding lepton number conserving Z′ interactions need not spoil low-scale resonant leptogenesis, a possibility worth testing in other inverse seesaw constructions.","The paper's use of a single complex modulus τ as the only flavor source points toward a geometric program: map which regions of the fundamental domain of τ produce each phenomenological signature, turning the current scan box into a full characterization of the model.","The combination of neutrino fit, LFV constraints, and leptogenesis within one parameter scan suggests that future data, especially a precise measurement of θ23 or the Dirac phase, could distinguish this A4 modular inverse seesaw from other TeV-scale seesaw frameworks."],"forward_implications":["The model predicts a lower bound sin² θ23 > 0.44, which current global fits already permit but which next-generation long-baseline experiments can directly test.","The linear relation between sin² θ12 and sin² θ23 in the 3σ region means a precise measurement of either angle sharpens the prediction for the other.","Parameter points that fit neutrino oscillations automatically satisfy the current LFV bounds from MEG and BaBar, so the model makes the absence of observed charged lepton flavor violation a consequence of the same structure that fixes neutrino mixing.","The final baryon asymmetry is essentially independent of the Z′ gauge coupling g_{B−L}, because stronger Z′-mediated scatterings delay the generation of asymmetry but do not suppress it, preserving Y_B ~ 8.6×10⁻¹¹.","The model predicts restricted ranges for the CP observables, including J_CP in [-0.07,0.07] and δ_CP in the intervals 0°–85° and 279°–359°, which future CP-violation measurements can check."],"supporting_citations":[{"why":"Supplies the (2,2) inverse seesaw realization and the resonant leptogenesis CP asymmetry formulas, including the decay widths and regulator, used throughout the leptogenesis analysis.","marker":"[34]"},{"why":"Supplies the minimal B−L extended Standard Model resonant leptogenesis framework, including the Boltzmann equations and the Z′-mediated scattering treatment adopted here.","marker":"[49]"},{"why":"Provides the Kadanoff-Baym-based regulator R_ij for the resonant CP asymmetry in the nearly degenerate heavy neutrino pair.","marker":"[50]"},{"why":"Supplies the NuFIT 5.2 global oscillation data, with 3σ ranges for mixing angles and mass-squared differences, used to filter the model parameter points.","marker":"[55]"},{"why":"Supplies the Planck 2018 observed baryon asymmetry and the bound on the sum of neutrino masses used as targets in the numerical analysis.","marker":"[37]"},{"why":"Provides the MEG-I experimental upper bound BR(μ→eγ) < 4.2×10⁻¹³ used to check LFV viability.","marker":"[65]"},{"why":"Provides the BaBar upper bound 3.3×10⁻⁸ on BR(τ→eγ) and BR(τ→μγ) used to check LFV viability.","marker":"[67]"},{"why":"Supplies the KamLAND-Zen current limit on the neutrinoless double beta decay effective mass, used to test the model's 0νββ prediction.","marker":"[57]"}],"fun_headline_variants":["A4 modular symmetry links neutrino mixing to matter asymmetry","One modulus yields neutrino angles and baryon asymmetry","Inverse seesaw with A4 predicts θ23 > 0.44 and leptogenesis","A4 model reproduces neutrino data and respects LFV limits","A4 modular symmetry: one modulus for neutrino mixing and baryogenesis"],"cache_read_input_tokens":21760,"weakest_assumption_plain":"The predictions depend on the randomly scanned parameter ranges for the modulus and couplings; if those ranges do not cover the model's full viable parameter space, the claimed exclusions and correlations could be artifacts of the scan rather than inherent predictions.","fun_headline_variants_meta":{"raw":{"variants":["A4 modular symmetry links neutrino mixing to matter asymmetry","One modulus yields neutrino angles and baryon asymmetry","Inverse seesaw with A4 predicts θ23 > 0.44 and leptogenesis","A4 model reproduces neutrino data and respects LFV limits","A4 modular symmetry: one modulus for neutrino mixing and baryogenesis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001031,"raw_usage":{"total_tokens":4336,"prompt_tokens":930,"completion_tokens":3406,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":3318}},"tokens_in":546,"tokens_out":3406,"duration_ms":24553,"temperature":1.0,"reasoning_tokens":3318,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:37:22.976403+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future measurement of sin² θ23 below 0.44 in the currently allowed 3σ range, or a dedicated scan covering the full fundamental domain of the modulus τ that finds oscillation-allowed points with sin² θ23 < 0.44, would falsify the model's lower-bound prediction. Similarly, a measured deviation from the predicted linear relation between sin² θ12 and sin² θ23 at the precision of upcoming experiments would rule out the correlation.","supporting_citations":[],"review_version":1}