{"id":"6a62f6d9-cbfe-4c7b-8a44-315cd96d7ac6","arxiv_id":"1909.00752","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In the one-singlet seesaw plus two-Higgs-doublet model, the neutrino Yukawa couplings can be expressed analytically through the measured neutrino masses, the PMNS matrix, and only two free parameters, phi' and lambda_D.","lead":"This paper studies a minimal neutrino mass model that adds one heavy right-handed neutrino and a second Higgs doublet to the Standard Model, giving mass to two light neutrinos through a combination of the seesaw mechanism and a one-loop radiative correction. It shows how to rewrite the model's neutrino Yukawa couplings in terms of measured neutrino masses and mixing angles plus two free parameters, making the allowed parameter space much easier to scan.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central inversion rests on the unquantified tree-level-plus-delta-M_L approximation (Eqs. 2.21-2.22); if neglected one-loop blocks are non-negligible, the two-parameter description is not the full model prediction.","rationale":"The reader's weakest_assumption correctly identifies the Grimus-Lavoura one-loop formula, Eq. (2.22), as the foundation of the analytic inversion. My reading of the paper confirms that the derivation of a, b, c and all subsequent formulas depends entirely on that approximation, and the paper nowhere quantifies the neglected terms. This is a genuine correctness risk because the central claim of testability is conditional on the approximation being numerically accurate across the sampled parameter space. The paper is otherwise algebraically consistent: I checked the determinant computation (Eqs. 3.29-3.35), the orthogonality argument for the massless state (Eq. 3.8 and footnote 2), and the parameter counting; these hold. The consistency check in Sec. 4.1 uses the same approximate formulas, so it cannot reveal a breakdown of the approximation. Thus the reader's CONDITIONAL verdict is appropriate, and I see no reason to move it to ACCEPT or REJECT. The concrete test of computing complete one-loop corrections for a benchmark would resolve whether the concern actually lands.","tokens_in":36343,"tokens_out":15064,"duration_ms":160737,"concrete_test":"For benchmark point B1 (Table 3) and three values of m_4 (e.g., 10^5, 10^10, 10^12 GeV), compute the full one-loop correction to the complete 4x4 neutrino mass matrix, including delta-M_D and delta-M_R and the charged-Higgs and W-boson contributions, using the complete Grimus-Lavoura/Pilaftsis formulas, with the Yukawa couplings obtained from Eqs. (3.7), (3.9), (3.35), and (3.43). Diagonalize the resulting 4x4 mass matrix and compare the light neutrino masses and PMNS matrix with the input values. If the deviations exceed the experimental 1-sigma ranges for any point, the inversion omits non-negligible contributions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's key claim is that the six complex parameters in Delta_1 and Delta_2 can be replaced by measured neutrino data plus two real parameters (Sec. 3, Summary). The derivation of a, b, c (Eqs. 3.12-3.14), and hence d (Eq. 3.35), |d'| (Eq. 3.43), and the Yukawa couplings (Eqs. 3.7, 3.9), uses only the one-loop correction delta-M_L evaluated at zero external momentum, while Eq. (2.21) explicitly drops delta-M_D and delta-M_R and the text states that charge-changing contributions are subdominant (Sec. 2.4). These are approximations inherited from Grimus-Lavoura, not proven within this paper, and no estimate is given for their size in the sampled parameter space (m_4 from 10^2 to 10^12 GeV, Higgs masses up to 3 TeV, Yukawa couplings spanning several orders of magnitude). If delta-M_D, delta-M_R, charged-scalar/W loops, or momentum dependence are not negligible for some benchmark points, then the Yukawa couplings computed from the inversion are not the true model predictions, so the claim that only phi' and lambda_D remain undetermined would fail. The consistency check in Sec. 4.1 recomputes masses and mixing angles from the derived Yukawas using the same approximate formulas, so it validates the algebra but not the approximation itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the Grimus-Neufeld model: type-I seesaw with a single right-handed singlet and a second Higgs doublet, so that one light neutrino mass arises at tree level and another at one loop. Using the Grimus-Lavoura one-loop approximation for the effective light neutrino mass matrix, the authors derive an analytic inversion in which the six complex parameters of the neutrino Yukawa couplings Δ1 and Δ2 are re-expressed in terms of the measured neutrino mass squared differences, the PMNS matrix, the heavy neutrino mass m4, the Higgs-sector parameters, and two additional real parameters λD and φ'. The derivation writes the effective 3x3 mass matrix in the tree-level seesaw basis, solves for the coefficients d and |d'| from the determinant and trace conditions, and then constructs Δ1 and Δ2 from the PMNS columns. The paper presents numerical distributions of d, |d'|, and the Yukawa couplings for a CP-conserving 2HDM Higgs sector subject to stability, unitarity, and oblique-parameter constraints, and it compares the analytic method with a numerical χ2 fitting approach. The central claim is that the neutrino Yukawa parameter space can be mapped analytically rather than by expensive scans.","tokens_in":36644,"tokens_out":9168,"duration_ms":87508,"significance":"If the central derivation is correct, the analytic parameterization is a useful technical result: it reduces the neutrino Yukawa sector of the Grimus-Neufeld model to two continuous parameters and provides explicit formulas, Eqs. (3.35) and (3.43), that are much cheaper to evaluate than a global fit over the twelve real Yukawa parameters. The paper is honest that neutrino masses and mixings are inputs, not predictions, and the algebraic inversion is internally consistent; the numerical consistency check confirms that the constructed couplings reproduce the input masses and angles. The main value is in mapping allowed Yukawa couplings for future phenomenological studies. The strength of the 'only two non-physical parameters' claim is, however, conditional on the one-loop approximation in Eqs. (2.21)-(2.22) and on the assumption in Eq. (3.8), so the paper should either quantify those limitations or temper the final testability statement.","major_comments":[{"comment":"The entire inversion in Sec. 3 is built on the approximation that the one-loop corrected neutrino mass matrix is Mν ≈ ((δML, M_D^T), (M_D, M_R)), with δML evaluated at zero external momentum and with charge-changing (charged-scalar/W) contributions neglected. Eq. (2.21) explicitly drops δM_D and δM_R, and the text says the charge-changing contributions are 'subdominant' but provides no estimate of their size over the scan range (m4 from 10^2 to 10^12 GeV, Higgs masses up to 3 TeV, Yukawa couplings spanning several orders of magnitude). Since Eqs. (3.35) and (3.43) determine d and |d'| from the f1, f2, f3 built out of this δML, a sizable neglected contribution would invalidate the central claim that only λD and φ' remain undetermined. The authors should either quantify δM_D, δM_R, the momentum dependence, and the charge-changing terms for representative benchmark points, or explicitly state that the parameterization is a property of the approximate one-loop formula rather than of the full model.","section":"Sec. 2.4, Eqs. (2.21)-(2.22)"},{"comment":"The consistency check in Sec. 4.1 recomputes the neutrino masses and mixing angles from the derived Yukawa couplings using the same approximate formulas, Eqs. (2.21)-(2.22), that were used for the inversion. The agreement therefore validates the algebraic inversion but not the underlying approximation. The Summary's statement that 'a few measurements that restrict the neutrino Yukawa couplings can confirm or rule out our model' is stronger than what this check establishes; the analysis uses the measured masses and the PMNS matrix as input, so it does not provide an independent prediction of neutrino observables. I recommend softening the testability claim or adding a genuinely independent test, for example by computing a loop-induced observable that was not used as input.","section":"Sec. 4.1"},{"comment":"The reduction of Δ2 to two complex parameters via Δ2 = d V_r^† + d' V_s^† assumes that the massless state V_o is orthogonal to Δ2. The footnote argues that a vector orthogonal to both Δ1 and Δ2 always exists, but the identification of that vector with the physical massless neutrino and with a specific column of the PMNS matrix is not demonstrated. If the loop-corrected mass matrix in the general model can have a Δ2 component along V_o, or if the massless eigenstate is not exactly the orthogonal vector, then the parameterization covers only a subclass of the model. The paper should prove that Eq. (3.8) is not a loss of generality, or state explicitly that the analysis applies to the subclass of parameter space satisfying this condition.","section":"Sec. 3, Eq. (3.8)"}],"minor_comments":[{"comment":"The labels 'for NH' are repeated, and the second line uses Δm21^2 while the first uses |Δm31^2|; please clarify which scenario of Table 2 each formula refers to.","section":"Sec. 3.1, Eqs. (3.52)-(3.55)"},{"comment":"The polynomial coefficients a4, a3, ... are denoted with the same symbol a that was used for the matrix element in Eq. (3.11); renaming the coefficients would avoid confusion.","section":"Sec. 3, Eqs. (3.37)-(3.43)"},{"comment":"Reference [31] appears incomplete: it has no source, journal, or arXiv identifier, and should be updated or removed.","section":"References"},{"comment":"The caption and text contain typos such as 'a finer study of is shown'; please correct the wording.","section":"Sec. 4.2, Fig. 3"},{"comment":"The text says the sum over k runs over all neutral physical Higgses, but the formula has three terms plus a separate Z term; please clarify whether the Goldstone boson contribution is included in the Z term or omitted.","section":"Sec. 2.4, Eq. (2.22)"},{"comment":"The paper uses m4 and M_R almost interchangeably before defining their relation in Eq. (3.46); a short convention statement near Eq. (2.22) would improve readability.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a borderline major-revision case. The algebraic parameterization is sound and potentially useful, and the numerical consistency check confirms the inversion. The two load-bearing issues are the unquantified one-loop approximation of Eqs. (2.21)-(2.22) and the assumption in Eq. (3.8). Both are addressable in a revision, so I would not recommend rejection; however, the current text overstates the testability of the model because the neutrino observables are inputs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nPunchline: this is a solid, honest reparameterization of the Grimus-Neufeld model. The genuinely new thing is the analytic inversion: instead of fitting six complex Yukawa parameters, the authors express them directly in terms of the measured neutrino mass differences, the PMNS matrix, and just two real numbers (λ_D and φ'). The algebra is internally consistent—they get d² from the determinant of the 2×2 block and |d'| from a quartic trace equation—and they verify numerically that the constructed couplings reproduce the inputs. They also show their method is roughly four hundred times faster than a differential-evolution fit for the same model. That's a practical win.\n\nWhat's good: the paper is transparent about what is input and what is output—the scare quotes around \"predict\" are appropriate—and the statistical maps of Yukawa couplings across m4 and λ_D are clearly presented. The CP-conserving Higgs case is worked out fully, and the appendices on b-vectors and PMNS parametrization are useful. The citation of Grimus-Lavoura and Ibarra-Simonetto is fair; the claim that the inversion is new seems correct.\n\nSoft spots, in proportion. The entire construction rests on the Grimus-Lavoura one-loop formula (Eq. 2.22): zero external momentum, and dropped δM_D, δM_R, and charge-changing pieces. The paper never quantifies these neglected terms in the sampled region (m4 from 10^2 to 10^12 GeV, Higgs masses up to 3 TeV). If they are non-negligible for some benchmark points, the two-parameter description is not the full model prediction. This is a real caveat, but it's an inherited standard approximation and the authors are explicit about using it. I wouldn't reject over this, but a referee should ask for an estimate of the neglected blocks. Minor issues: no error propagation (best-fit values only), the λ_D range [1/2,2] is an educated guess, and no code or data are released.\n\nThe consistency check in Sec. 4.1 validates the algebra, not the approximation—the authors don't overstate it. Net: the central claim holds within the stated approximation.\n\nWho this is for: phenomenologists working on minimal seesaw plus 2HDM, especially anyone scanning this model's parameter space. It deserves a serious referee. The right referee will focus on the domain of validity of Eq. (2.22) and will ask for error estimates, but the paper is honest, well-organized, and the new analytic result is publishable.\n\nBest.","headline":"A solid, honest reparameterization of the Grimus-Neufeld model; the analytic inversion is genuinely new, but the one-loop approximation underpinning it is unquantified.","tokens_in":37188,"tokens_out":4990,"would_cite":false,"duration_ms":180805,"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":"One-loop corrections turn the seesaw model's six complex neutrino Yukawa couplings into functions of two free parameters.","keywords":["seesaw mechanism","two-Higgs-doublet model","radiative neutrino mass","neutrino Yukawa couplings","Grimus-Neufeld model","PMNS matrix","CP-conserving Higgs potential","one-loop self-energy"],"falsifier":"Compute the neutrino self-energies at nonzero external momentum $p^2$, including charged-Higgs and $W$-boson loops, at the benchmark point B1 with $m_4=10^{10}\\,\\mathrm{GeV}$ and $\\lambda_D=1$, then re-solve for $d$ and $|d'|$ from the resulting mass matrix; if the reconstructed neutrino mass-squared differences deviate from $\\Delta m^2_{21}$ and $|\\Delta m^2_{31}|$ by more than the experimental errors, the zero-momentum, neutral-only approximation that underpins Eqs. (3.12)-(3.14) is not subdominant and the parametrization fails.","tokens_in":36155,"feed_emoji":"⚛️","tokens_out":6630,"duration_ms":61476,"temperature":0.7,"pith_summary":"This paper studies the Grimus-Neufeld model, a seesaw extension of the Standard Model with one right-handed neutrino and a second Higgs doublet that generates a radiative neutrino mass at one loop. Its central claim is that, in the CP-conserving two-Higgs-doublet case, the six complex neutrino Yukawa couplings are not free parameters: they can be expressed analytically through the measured neutrino mass differences, the PMNS matrix, the heavy neutrino mass, and two remaining parameters, the phase $\\phi'$ and the ratio $\\lambda_D$. If correct, this converts the model's neutrino sector into a nearly closed prediction, reducing a twelve-parameter numerical fit to a two-parameter scan and making the model directly testable by any future measurement that constrains a neutrino Yukawa coupling. The paper demonstrates the approach numerically, including a consistency check that reproduces the input mass differences and mixing angles.","feed_headline":"Neutrino Yukawa couplings collapse to two free parameters","feed_subtitle":"Measured mass differences and PMNS angles fix all six couplings, leaving a phase and a mass ratio free.","key_machinery":"The load-bearing object is the effective $3\\times3$ light-neutrino mass matrix in the basis where the tree-level seesaw is diagonal: after the unitary matrix $V$ from the tree-level diagonalization, the one-loop corrected matrix has the rank-2 block form $\\mathrm{diag}(0,M_{2\\times2})$ with $M_{2\\times2}=\\begin{pmatrix}a&b\\\\b&c\\end{pmatrix}$. The entries $a,b,c$ are quadratic in the Yukawa coefficients $d,d'$ and the loop functions $f_1,f_2,f_3$ of Eqs. (3.12)-(3.19), where these functions encode the neutral-Higgs and $Z$ self-energies through $L(m^2)$ and the $b$-vectors of the Higgs mass eigenfields. The determinant condition fixes $d$, while the trace condition on $M_{2\\times2}^\\dagger M_{2\\times2}$ produces a fourth-order polynomial in $|d'|$; together with the Takagi factorization of the $2\\times2$ block, this turns the measured mass-squared differences and the PMNS matrix into inputs that determine the Yukawa couplings.","core_discovery":"The paper establishes an inversion of the one-loop seesaw formula. Using the Grimus-Lavoura approximation for the effective light neutrino mass matrix, the tree-level seesaw leaves one massless and one massive state; at one loop the $2\\times2$ block $\\begin{pmatrix}a&b\\\\b&c\\end{pmatrix}$ is built from Yukawa coefficients $d$, $d'$, and loop functions $f_1,f_2,f_3$. The determinant relation $ac-b^2=d^2(2m_D^2/v^2)(f_1f_3-f_2^2)$ gives $d^2$ directly, and the trace relation for $M_{2\\times2}^\\dagger M_{2\\times2}$ yields a fourth-order polynomial in $|d'|$ whose real positive roots are the allowed couplings. Interpreting the orthogonal vectors as columns of the PMNS matrix then expresses the neutrino Yukawa couplings $\\Delta_1$ and $\\Delta_2$ through Eqs. (3.7) and (3.9). With the lightest neutrino massless, the two measured mass-squared differences and the PMNS matrix fix the Yukawa couplings up to the phase $\\phi'$ and the scaling parameter $\\lambda_D$.","pith_inferences":["A single future measurement that pins down one neutrino Yukawa coupling, for example a charged-lepton-flavour-violating rate or a collider signature involving the second Higgs doublet, would overconstrain the two-parameter plane and could by itself rule the model in or out.","The determinant-plus-trace trick is not specific to $n_R=1$: for two right-handed singlets one would instead solve a small polynomial system for the Yukawa coefficients, and it may be possible to reduce similarly large parameter spaces to a handful of physical inputs.","The observed speed-up suggests that the same inversion strategy could turn global fits of other minimal radiative-seesaw models into scans over one or two physically meaningful parameters."],"forward_implications":["If the parametrization is correct, the neutrino Yukawa couplings are no longer scan parameters: choosing $\\phi'$ and $\\lambda_D$ together with the Higgs-sector masses and mixing angle produces a specific coupling matrix, so the model can be confronted directly with any observable sensitive to these couplings.","A scan over the two remaining parameters replaces a twelve-parameter numerical minimization; on the benchmark point the analytical method was about 430 times faster than the differential-evolution fit while giving the same distribution of $|\\Delta_{23}|$.","The remaining freedom is mild: $\\lambda_D$ is restricted to $[1/2,2]$ by loop-expansion sanity, and $\\phi'$ is often confined to narrow intervals where the quartic has real positive roots.","The lightest neutrino stays massless at one loop, so the two measured mass-squared differences fix the masses directly for both hierarchies, keeping the sum of neutrino masses in agreement with the Planck bound.","The procedure extends in principle to a CP-violating Higgs potential by allowing $f_1$ to become complex, with the same qualitative structure."],"supporting_citations":[{"why":"Introduces the model with one right-handed singlet and a second Higgs doublet, whose radiative neutrino masses this paper parametrizes.","marker":"[6]"},{"why":"Supplies the one-loop effective neutrino mass formula (Eq. 2.22) that the analytic inversion starts from.","marker":"[7]"},{"why":"Provides the companion formulation of soft lepton flavour violation in the multi-Higgs seesaw framework used for the Yukawa structure.","marker":"[8]"},{"why":"Provides the experimental neutrino mass-squared differences and oscillation angles used as input to fix the PMNS matrix.","marker":"[17]"},{"why":"Supplies the basis-independent CP-conserving 2HDM parameterization and the constraints on mixing angles used in the numerical analysis.","marker":"[19]"},{"why":"Provides the Takagi factorization routine used to diagonalize the complex symmetric light neutrino mass matrix.","marker":"[25]"},{"why":"Establishes that the lightest neutrino remains massless at one loop in this model, justifying the use of two measured mass-squared differences as input.","marker":"[40]"}],"fun_headline_variants":["One-loop seesaw leaves neutrinos with just two free parameters","Seesaw with one singlet and second Higgs: Yukawas pinned by data","Two parameters control all neutrino Yukawa couplings","Radiative seesaw shrinks neutrino parameter space to two angles","Extra Higgs doublet makes neutrino masses mostly determined"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inversion stands on the Grimus-Lavoura one-loop formula Eq. (2.22), which evaluates the neutral-Higgs and $Z$ self-energies at zero external momentum and drops charged-current contributions; if those neglected terms are not small, the derived couplings are not the model's true predictions.","fun_headline_variants_meta":{"raw":{"variants":["One-loop seesaw leaves neutrinos with just two free parameters","Seesaw with one singlet and second Higgs: Yukawas pinned by data","Two parameters control all neutrino Yukawa couplings","Radiative seesaw shrinks neutrino parameter space to two angles","Extra Higgs doublet makes neutrino masses mostly determined"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000185,"raw_usage":{"total_tokens":1304,"prompt_tokens":913,"completion_tokens":391,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":529,"completion_tokens_details":{"reasoning_tokens":306}},"tokens_in":529,"tokens_out":391,"duration_ms":4654,"temperature":1.0,"reasoning_tokens":306,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:38:20.473842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the neutrino self-energies at nonzero external momentum $p^2$, including charged-Higgs and $W$-boson loops, at the benchmark point B1 with $m_4=10^{10}\\,\\mathrm{GeV}$ and $\\lambda_D=1$, then re-solve for $d$ and $|d'|$ from the resulting mass matrix; if the reconstructed neutrino mass-squared differences deviate from $\\Delta m^2_{21}$ and $|\\Delta m^2_{31}|$ by more than the experimental errors, the zero-momentum, neutral-only approximation that underpins Eqs. (3.12)-(3.14) is not subdominant and the parametrization fails.","supporting_citations":[{"cited_title":"Grimus and H","cited_arxiv_id":null,"evidence_quote":"Introduces the model with one right-handed singlet and a second Higgs doublet, whose radiative neutrino masses this paper parametrizes."},{"cited_title":"Gauge dependence of tadpole and mass renormalization for a seesaw extended 2HDM","cited_arxiv_id":"1806.04675","evidence_quote":"Establishes that the lightest neutrino remains massless at one loop in this model, justifying the use of two measured mass-squared differences as input."}],"review_version":1}