{"id":"bb4182cc-2b0c-4390-ad0a-7120e745f47f","arxiv_id":"2505.12564","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The bilinear formalism of the two-Higgs-doublet model is extended to one-loop order, yielding gauge-invariant expressions for the quantum-corrected scalar masses and vacuum shift.","lead":"This paper develops a gauge-invariant way to compute one-loop quantum corrections to the scalar masses in the two-Higgs-doublet model, using bilinear field combinations and a perturbative expansion in powers of the loop factor. It provides compact formulas that avoid gauge-dependent artifacts and apply to any two-Higgs-doublet model.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Goldstone-IR divergence cancellation (Eq. 5.82) is argued via a general reference, not demonstrated for this formalism; the one-loop mass formulas (7.5)-(7.6) need an explicit check of that cancellation and an independent numerical cross-check.","rationale":"The reader's weakest assumption (perturbative hbar-expansion about the tree-level vacuum) is structurally important, but the more pressing concern after reading the full text is the explicit cancellation of Goldstone infrared divergences, which the paper only argues for by reference to a generic result. The numerical example is not cross-checked against any existing THDM one-loop computation. These are not fatal flaws, but they justify a CONDITIONAL verdict: the framework is plausible and detailed, but the cancellation of the spurious Goldstone IR divergences and the numerical implementation must be verified before the final formulas (7.5), (7.6), and (7.19) are used as black boxes. I partially agree with the reader: the perturbative expansion around the tree-level vacuum is a genuine structural assumption, but the specific and testable issue of Goldstone IR cancellation is more central to the main result as presented, so my concern overlaps only partially with the reader's. The recommendation is to keep the CONDITIONAL verdict, not to strengthen it, because the known cancellations in this context are plausible and the burden is a concrete computational check rather than a demonstrated error.","tokens_in":40419,"tokens_out":2110,"duration_ms":21196,"concrete_test":"Evaluate the Goldstone-loop contributions in Eq. (5.82) for the type-I CP-conserving example of Sec. 8 with a small Goldstone mass xi as a regulator: compute the rotation matrices, the couplings in Eq. (7.20), and the sum of all B(0,0)-type terms that arise in the xi -> 0 limit. Check whether the coefficient of log(xi^2/mu^2) vanishes identically when the complete one-loop mass expression (7.6) is assembled, including the kappa(1) shifts and the combination entering the secular equation (6.22). If the divergent coefficient does not vanish, the spurious Goldstone IR divergence is not cancelled, contradicting the paper's claim. In addition, recompute the one-loop masses of Fig. 1 for tan(beta)=1 with a standard THDM one-loop code to verify the numerical formulas.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central computational claim is that one-loop pole masses can be obtained from the hbar-expanded effective potential at the tree-level minimum. The most fragile step is the treatment of the Goldstone-boson infrared-divergent terms in the scalar contributions. In Eq. (5.82), the terms proportional to B(0,0) (arising from Goldstone loops) are logarithmically divergent. In Sec. 6, the paper invokes the generic argument of [34] about cancellation in the secular equation, but it does not explicitly demonstrate, for the derived couplings and the hbar-expanded mass matrix, that the divergent B(0,0) terms cancel when the one-loop scalar mass matrix is combined with the self-energy at p^2=0. If this cancellation fails for the specific hbar-expansion formula (7.6), the one-loop scalar masses would be scheme-dependent, and the central claim of providing a gauge-invariant one-loop formalism would not hold. A second, related issue is that the numerical example in Sec. 8 is not compared to any independent one-loop THDM computation, so a systematic error in the lengthy algebraic expressions would not be detected. The strongest claim is not settled by internal consistency alone; it requires an explicit check of the Goldstone cancellation and an independent numerical verification.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a gauge-invariant, bilinear-formalism treatment of one-loop corrections in the general Two-Higgs-Doublet Model (THDM). The authors review the bilinear formulation of the THDM potential, gauge sector, and Yukawa sector, then combine it with the ℏ-expansion (Sec. 4) to solve the one-loop stationary-point equations iteratively around the tree-level vacuum. Sections 5 and 6 derive the gauge, fermionic, and scalar contributions to the one-loop effective potential and its first and second bilinear-field derivatives, and express the one-loop shifts of the charged and neutral scalar masses in terms of these quantities. The central results are Eqs. (7.5) and (7.6), which give m_H±^2 and m_a^2 at next-to-leading order, together with the collection of all required derivatives and couplings in Sec. 7. The method is applied in Sec. 8 to a CP-conserving type-I THDM, with a numerical spectrum shown in Fig. 1 as a function of tanβ.","tokens_in":1932,"tokens_out":2277,"duration_ms":197896,"significance":"If the central formulas (7.5)–(7.6) are correct, the paper provides a compact and manifestly gauge-invariant route to one-loop THDM scalar masses that avoids gauge-fixing of the Higgs fields and Goldstone-sector mixing, and it is directly applicable to any THDM. The manuscript has genuine strengths: the derivations are detailed and systematic; the summary section (Sec. 7) is self-contained; the approach is parameter-free in the sense that no parameter is fitted to the one-loop output (tree-level masses, tanβ, and cos(α−β) are used as inputs); and the appendices provide the nontrivial second-derivative formulae. The main unresolved points concern the infrared behaviour of the final mass formula and the absence of an external numerical validation, both of which are load-bearing for the paper's central claim.","major_comments":[{"comment":"The cancellation of the spurious Goldstone-boson divergences is asserted but not demonstrated for the specific formula that is the paper's main result. In Eq. (5.82) (reproduced in Eq. (7.19)) the terms proportional to B(0,0) have coefficients (1/2)(λ̄_{aG0G0}λ̄_{bG0G0} + λ̄_{aG±G±}λ̄_{bG±G±}); using the couplings of Eq. (7.20), λ̄_{aG0G0} and λ̄_{aG±G±} are each proportional to m_a^2 k̄_a / √(2K0), so this coefficient is generically nonzero at a charge-conserving vacuum. The remaining contributions to (m_a^2)^(1) in Eq. (7.6) — the vacuum-shift terms (K_0^(1)/K_0^(0))(m_a^2)^(0) and f̄^±_a δ̄^(1)_a — depend on D̄_a V^(1) through Eq. (7.1) and involve only the A-functions, so they cannot cancel the B(0,0) divergence. The discussion in Sec. 6 (Eqs. (6.22)–(6.27)) is the standard argument, following Ref. [34], that on-shell pole masses are IR finite once the momentum-dependent self-energy combination Π(p^2)−Π(0) is included; it does not show that the zero-momentum object (7.6) is that pole mass. In the computational flow the divergence is instead disposed of by Eq. (7.28), which sets Bs(0,0)=0 — a prescription that conflicts with the definition B(x,x)=log(x/µ^2) in (5.18) and that changes the numerical result by the finite part of the Goldstone loop unless an independent cancellation is proved. As written, (7.6) is an MS-bar zero-momentum curvature mass, and its IR finiteness and gauge invariance do not follow from the arguments given. The authors should either provide the explicit pole-mass formula, including the Π(p^2)−Π(0) contributions that replace B(0,0) by the finite B(p^2,0,0), or prove that the B(0,0) coefficients in the final combination (7.6) vanish identically.","section":"Secs. 5.3, 6, 7 (Eqs. (5.82), (6.22)-(6.27), (7.6), (7.28))"},{"comment":"The numerical application is not validated against any independent one-loop THDM calculation. The one-loop masses in Fig. 1 follow from many pages of algebraic expressions (Secs. 5–7 and Appendices A–B), and internal consistency cannot detect a global sign error or a missing factor in, for example, the scalar couplings of Eq. (7.20) or the fermionic coefficients of Eq. (5.70). I recommend adding a comparison with an existing THDM one-loop spectrum for a benchmark point — for instance by recomputing the same type-I scenario with a public code or with the one-loop formulae of an independent calculation (see Ref. [5] and references therein) — and reporting the numerical one-loop shifts, not only the plot, so that the central claim of Eqs. (7.5)–(7.6) can be checked.","section":"Sec. 8.2, Fig. 1"}],"minor_comments":[{"comment":"The third displayed equation of (7.19) has unbalanced braces and a stray period, and the quantities λ0aa and λ0H±H± are written without the bar used in Eq. (5.85); please align the notation with (5.85)–(5.86).","section":"Sec. 7, Eq. (7.19)"},{"comment":"The value Bs(0,0)=0 should be introduced as an explicit regularization prescription with a justification, rather than as a value that follows from Eqs. (5.18)–(5.19), since the limit x→0 of log(x/µ^2) is not defined there.","section":"Sec. 7, Eq. (7.28)"},{"comment":"The constraint |α|=|β|=π/2 combined with cos(α−β)=1 is confusing, since the latter relation already suggests α=β; moreover the barred four-vector f̄K in Eq. (8.22) is used before it is defined. Please clarify the alignment conditions and the rotation that defines f̄K.","section":"Sec. 8.1, Eqs. (8.20) and (8.22)"},{"comment":"The manuscript explains that results are first obtained in the MS-bar scheme and that one should switch to the on-shell scheme for pole masses, but Fig. 1 does not state which scheme the plotted curves correspond to; please specify the scheme, the renormalization scale, and the scale dependence of the one-loop shifts.","section":"Secs. 6–8, Fig. 1"},{"comment":"The running title contains 'T wo-Higgs' with an extra space, and there are a few similar formatting artifacts; these should be cleaned up before submission.","section":"Title and general text"}],"recommendation":"major_revision","confidential_remarks":"The paper's central claim rests on the IR-divergence resolution discussed in Sec. 6, which is currently argued at the level of principle (following Ref. [34]) rather than demonstrated for the formulas actually implemented in Eqs. (7.5)–(7.6) and used in Fig. 1. The heavy overlap with Refs. [7,19] by the same group is not circular, since those results are independent and previous, but an independent numerical cross-check would substantially strengthen the case. If the authors can supply the explicit pole-mass formula and a benchmark comparison against an independent calculation, I expect the paper to be publishable in JHEP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it says: it extends the bilinear formalism to one loop using the ℏ-expansion and packages the results as compact, ready-to-use formulas. If you work in the THDM, having the scalar mass corrections in gauge-invariant bilinear form is a real convenience, and the authors deserve credit for carrying the algebra through and collecting it cleanly in Sec. 7. The tree-level masses are used as inputs, not fitted to the output, so there is no circularity problem.\n\nThe structure is coherent and the internal derivations are careful. The weak points are exactly where the stress-test note lands. First, the cancellation of the B(0,0) Goldstone divergences is argued by citing the general secular-equation logic of Denner, but the paper never shows, for its own ℏ-expanded mass matrix and couplings, that the divergent terms in Eq. (5.82) actually drop out of the physical masses. That is not a fatal flaw -- the argument is plausible and the mechanism is standard -- but it is a gap, and an explicit check in this formalism would close it.\n\nSecond, the numerical example in Sec. 8 has no comparison to any existing one-loop THDM calculation. Since the algebraic expressions are long, a systematic sign or factor error could easily hide. A single benchmark compared against, say, a conventional diagrammatic one-loop calculation would substantially raise confidence. Without that, the central claim rests on internal consistency, which is good but not conclusive.\n\nOne more thing: the paper overstates the novelty a little. The components -- bilinears, ℏ-expansion, gauge-invariant effective potential -- all exist. What is new is the combination and the explicit formulas, which is enough to justify publication but not to call it a breakthrough.\n\nOverall: the framework is plausible, the formulas are likely correct, and the paper deserves refereeing. A serious referee should ask for two additions: an explicit demonstration of the Goldstone cancellation and an independent numerical comparison. If those are provided, I would be comfortable with the result. Who is this for? THDM phenomenologists doing precision studies -- they will want the formulas in Sec. 7. I would not cite it in my own work right now, but I would bring it to the attention of people who might.","headline":"A genuinely useful gauge-invariant one-loop THDM toolkit, but the Goldstone-IR cancellation is asserted rather than demonstrated and the numerical example has no independent cross-check.","tokens_in":41274,"tokens_out":1048,"would_cite":false,"duration_ms":13903,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.60.Fr","11.30.Qc","14.80.Cp"],"model":"deepseek-v4-flash","headline":"One-loop corrections to the two-Higgs-doublet scalar masses can be computed in a manifestly gauge-invariant way, using bilinears and the $\\hbar$-expansion, without fixing the Higgs gauge or dealing with Goldstone mixing.","keywords":["two-Higgs-doublet model","bilinear formalism","gauge invariance","one-loop corrections","hbar-expansion","effective potential","Goldstone infrared divergences","scalar mass spectrum"],"falsifier":"For a fixed benchmark of a CP-conserving THDM, compute the one-loop pole masses through the secular equation (6.22) with explicit self-energies in a general $R_\\xi$ gauge and compare with equations (7.5)--(7.6). Residual dependence of the physical masses on the gauge-fixing parameter, or non-vanishing Goldstone masses, would falsify the paper's claim.","tokens_in":40163,"feed_emoji":"⚛️","tokens_out":8216,"duration_ms":78019,"temperature":0.7,"pith_summary":"The paper establishes that one-loop corrections to the scalar mass spectrum of the general two-Higgs-doublet model (THDM) can be computed while keeping electroweak gauge invariance manifest at every step. Its method combines the bilinear formalism, four real gauge-invariant combinations of the two Higgs doublets, with the $\\hbar$-expansion, an iterative solution of the stationarity conditions around the tree-level vacuum. The payoff is that no gauge choice for the Higgs fields is needed, and the problematic Goldstone-sector mixings and infrared divergences are avoided or cancel explicitly. If correct, this gives a compact set of analytic formulas, equations (7.5) and (7.6), applicable to any THDM and useful for precision tests such as the electroweak $\\rho$ parameter.","feed_headline":"One-loop THDM masses, gauge-invariant at every step","feed_subtitle":"Four bilinears plus the hbar-expansion yield closed-form corrections and keep Goldstone bosons massless.","key_machinery":"The central object is the $2\\times 2$ positive-semidefinite bilinear matrix $K = \\psi\\psi^\\dagger$, whose four real components $K_0, K_1, K_2, K_3$ are gauge-invariant bilinears. At a charge-conserving minimum $K$ has rank 1, which encodes correct electroweak symmetry breaking and makes the neutral mass matrix expressible as $\\gamma_3(M - 2u\\,\\tilde{g})\\gamma_3^T$. The second ingredient is the $\\hbar$-expansion with $\\kappa = 1/(16\\pi^2)$: the potential and vacuum are expanded around their tree-level values, and the stationarity conditions are solved order by order, with derivatives $D_a = \\gamma_{a\\nu}\\partial_\\nu$ projecting onto the neutral mass eigenstates. These two ingredients turn the one-loop calculation into a small set of algebraic formulas involving tree-level masses, Yukawa bilinears, and gauge-boson bilinears.","core_discovery":"The central claim is that the one-loop corrected charged and neutral Higgs masses are given by $m^2_{H^\\pm} = (m^2_{H^\\pm})^{(0)} + \\kappa (m^2_{H^\\pm})^{(1)} + O(\\kappa^2)$ and $m_a^2 = (m_a^2)^{(0)} + \\kappa (m_a^2)^{(1)} + O(\\kappa^2)$, where the first-order shifts follow from gauge-invariant derivatives of the one-loop effective potential evaluated at the tree-level minimum. Writing the full THDM Lagrangian in terms of bilinears makes the all-order scalar mass matrix block-diagonal in a canonical basis, with the charged-pair mass equal to $4uK_0$ at every order; the $\\hbar$-expansion then converts the one-loop stationarity condition into simple algebraic equations for the vacuum shifts and mass corrections. The paper provides explicit closed-form expressions for the gauge, fermion, and scalar contributions to these shifts in Sec. 7.","pith_inferences":["Beyond the paper: the same vacuum-shift formulas should feed directly into gauge-invariant one-loop predictions for the electroweak oblique parameters $S$, $T$, $U$, which are not computed here.","Beyond the paper: the rank-1 vacuum condition gives a built-in consistency check for numerical implementations, since physical masses must stay gauge-independent and Goldstone modes must remain massless.","Beyond the paper: because the $\\hbar$-expansion is iterative, the structure carries to two loops whenever a two-loop effective potential in bilinears becomes available, with no new gauge-fixing step required.","Beyond the paper: one could test the method's reach by benchmarking it against conventional diagrammatic one-loop calculations for non-CP-conserving or type-II Yukawa scenarios, where the $|\\xi_{ud}|^2$ terms are active."],"forward_implications":["The one-loop shifted masses in (7.5) and (7.6) are direct analytic predictions for any THDM once the tree-level vacuum and masses are known.","The charged-Higgs mass relation $m^2_{H^\\pm} = 4uK_0$ holds at every perturbative order, so its one-loop correction is fixed entirely by the shifts in $u$ and $K_0$.","Goldstone bosons stay massless order by order and never mix with physical scalars, and the usual IR-divergent Goldstone contributions cancel when physical pole masses are extracted through the secular equation.","The gauge, fermion, and scalar contributions to the shifts separate into modular closed-form expressions involving the integrals $A$ and $B$, making phenomenological scans straightforward.","The method applies to any THDM, including models with explicit or softly broken $Z_2$ symmetries and different Yukawa types, once the parameters are expressed in bilinear form."],"supporting_citations":[{"why":"Introduces the bilinear formalism, the charge-conserving rank condition, and the tree-level stability and electroweak-symmetry-breaking analysis used throughout.","marker":"[7]"},{"why":"Derives the all-order gauge-invariant scalar mass matrices in the canonical basis, including the charged-mass relation m^2_Hpm = 4uK0, which this paper extends to one loop.","marker":"[19]"},{"why":"Introduces the hbar-expansion method for iteratively solving stationarity conditions around the tree-level solution.","marker":"[30]"},{"why":"Provides the one-loop effective potential formula that all gauge, fermion, and scalar contributions start from.","marker":"[33]"},{"why":"Supplies the secular equation and self-energy conventions used to show that infrared Goldstone divergences cancel in physical pole masses.","marker":"[34]"},{"why":"Provides the bilinear expression for the fermion mass matrices squared used in the fermionic one-loop contributions.","marker":"[15]"}],"fun_headline_variants":["Gauge-invariant one-loop THDM masses","Bilinear formalism yields one-loop THDM corrections","No gauge fixing for one-loop THDM masses","Closed-form one-loop THDM corrections, gauge-invariant"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that quantum corrections shift the vacuum only by small perturbative steps away from the tree-level, electrically neutral minimum; a true one-loop minimum of a different character, such as charge-breaking or far away, would be missed.","fun_headline_variants_meta":{"raw":{"variants":["Gauge-invariant one-loop THDM masses","Bilinear formalism yields one-loop THDM corrections","No gauge fixing for one-loop THDM masses","Closed-form one-loop THDM corrections, gauge-invariant"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000435,"raw_usage":{"total_tokens":2211,"prompt_tokens":937,"completion_tokens":1274,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":1213}},"tokens_in":553,"tokens_out":1274,"duration_ms":9091,"temperature":1.0,"reasoning_tokens":1213,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:31:48.136433+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a fixed benchmark of a CP-conserving THDM, compute the one-loop pole masses through the secular equation (6.22) with explicit self-energies in a general $R_\\xi$ gauge and compare with equations (7.5)--(7.6). Residual dependence of the physical masses on the gauge-fixing parameter, or non-vanishing Goldstone masses, would falsify the paper's claim.","supporting_citations":[{"cited_title":"The hbar Expansion in Quantum Field Theory","cited_arxiv_id":"1009.2313","evidence_quote":"Introduces the hbar-expansion method for iteratively solving stationarity conditions around the tree-level solution."},{"cited_title":"Coleman and E.J","cited_arxiv_id":null,"evidence_quote":"Provides the one-loop effective potential formula that all gauge, fermion, and scalar contributions start from."},{"cited_title":"Multiple point principle in the general Two-Higgs-Doublet model","cited_arxiv_id":"2001.10541","evidence_quote":"Provides the bilinear expression for the fermion mass matrices squared used in the fermionic one-loop contributions."}],"review_version":1}