{"id":"304fc4bb-df6a-4ad8-87b2-f66f907dfd8d","arxiv_id":"2412.13896","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a lattice simulation of the custodial two-Higgs-doublet model, an extra BSM scalar can be as light as about 0.2 times the W boson mass while Standard Model physics is held fixed.","lead":"This paper uses large-scale computer simulations to study a model with two Higgs doublets, one of which plays the role of the Standard Model Higgs. The main finding is that an extra scalar can be as light as about 0.2 times the W boson mass, roughly 16 GeV, while Standard Model physics stays fixed up to energies of 600 GeV.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.2 m_W 'plateau' may be a finite-volume floor: for Table 1 volumes m_H L ≈ 0.9–1.7, and no spectrum volume test is shown, so cutoff independence is not established.","rationale":"The reader's weakest assumption concerns whether the LPCP remains valid when kappa1 is scanned, because only R and S are monitored. That concern is legitimate, but the paper's own Fig. 4(left) provides at least a coarse check, and small quartic couplings make large retuning effects unlikely at the level of the claim. The more directly load-bearing issue is that the numerical value of the BSM mass is read off from two-point functions on lattices where the H correlation length is comparable to the box size. m_H L ranges from about 0.94 to 1.7, and at the beta = 8.56 point a m_H is within 10% of 1/L. Near the H12 transition the light scalar's mass in a finite box is governed by finite-size scaling, not by the infinite-volume spectrum, so a plateau is exactly what one would expect from the box floor. Since the quoted 0.2 m_W is close to 1/(m_W L) for two of the five points, the 'cutoff independence' could be a consequence of all boxes having m_W L in a similar range. This is checkable by one additional simulation per beta. I therefore keep the reader's CONDITIONAL verdict unchanged, but for a different reason: the manuscript must demonstrate volume independence of the BSM masses before the central bound can be accepted. The absence of released code or data also prevents independent checks, but the volume test is the specific missing analysis that would settle the most load-bearing concern.","tokens_in":8972,"tokens_out":11692,"duration_ms":107562,"concrete_test":"Repeat the zero-temperature spectrum determination at beta = 8.56 on a 48^4 or 64^4 lattice and at beta = 8.64 on a 64^4 lattice, keeping all bare couplings of Table 1 fixed and using the same S^4_12 and W^4_12 correlators of Eq. (3). If the resulting m_H/m_W or m_A/m_W shifts by more than ~10% from the L = 32 and L = 48 values, or if the plateau in Fig. 4 (right) moves downward, the claimed 0.2 m_W lower bound and its cutoff independence are finite-volume artifacts rather than continuum physics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract and §5) that the BSM scalar H has a lower mass m_H ~ 0.2 m_W independent of the cutoff rests on the plateau in Fig. 4 (right) as kappa1 approaches kappa_c1. No finite-volume study of the BSM spectrum is reported. From Table 1, m_H L ≈ (m_H/m_W)(m_W L) ≈ 1.7, 1.5, 1.3, 0.94 and 1.2 for the five LPCP points; at beta = 8.56 (L = 32, m_W L = 4.694) the H correlation length is a m_H ≈ 0.029, close to 1/L = 0.031. In a finite periodic box the extracted ground-state mass cannot decrease much below ~1/L, so as kappa1 approaches the H12 transition the effective mass will appear to saturate even if the infinite-volume mass continues to fall. The observation that m_A/m_W 'seems to keep decreasing' while m_H/m_W plateaus is consistent with the lighter H hitting the finite-volume floor first. The apparent cutoff independence could therefore be an artifact of boxes with physical extents of only 4.7 to 8.5/m_W. A volume test of the two-point functions in Eq. (3) is a prerequisite for the 0.2 m_W bound.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports exploratory lattice simulations of the custodial/inert two-Higgs-doublet model with SU(2) gauge fields. It defines a line of partially constant physics (LPCP) by fixing R = m_h/m_W ≈ 1.5 and g_GF^2(m_W) = 0.5, with cutoffs Λ ranging from about 300 to 630 GeV. Along this line the authors compute the SM spectrum, the gradient-flow running gauge coupling, the BSM mass ratios m_H/m_W and m_A/m_W as functions of the unbroken hopping parameter κ1, and the finite-temperature susceptibility of the Higgs angular variable. The central claim, stated in the abstract and Section 5, is that for small BSM quartic couplings the lighter BSM scalar has m_H ~ 0.2 m_W and that this value is independent of the cutoff.","tokens_in":9228,"tokens_out":6579,"duration_ms":65185,"significance":"If established, the main result would be a non-perturbative counterpoint to typical tree-level expectations: a weakly coupled 2HDM could contain a BSM scalar considerably lighter than the W boson, which is relevant for model building and for future searches. The paper also provides useful first lattice data on the running coupling and on the screening mass in this model. Strengths include the use of two independent SM quantities to define the LPCP, the gradient-flow definition of the coupling, and the clear description of the parameter space sectors. The principal limitations are the absence of a quoted uncertainty for the central m_H ~ 0.2 m_W value and the lack of any finite-volume test for the BSM spectrum; as a result the claim is plausible but not yet demonstrated to the standard needed for a quantitative mass bound.","major_comments":[{"comment":"The central claim m_H ~ 0.2 m_W is quoted without any uncertainty in the abstract and conclusion, and Fig. 4 (right) is not accompanied by a table of the plateau values or fit errors. Since the LPCP tuning quantities R and S in Table 1 carry relative errors of several percent, and the mass ratios are derived from correlator fits on finite statistics, the statement that the value is 'independent of the cutoff' needs numerical values with errors, e.g. m_H/m_W = 0.21(3) for each β. Without an explicit error budget and a χ²/dof for the plateau, the abstract overstates the precision of the result.","section":"Abstract; Section 4, Fig. 4 (right); Section 5"},{"comment":"No finite-volume study of the BSM two-point functions is reported. For the claimed m_H/m_W ~ 0.2, the five LPCP points have m_H L ~ 1.7, 1.5, 1.3, 0.94 and 1.2 (using m_W L from Table 1). These values are in a regime where finite-volume corrections to a ground-state mass extracted from a zero-momentum correlator can be substantial, and the volumes change together with β (L = 28, 28, 32, 32, 48). A volume dependence test at one or two β values, or at least an estimate of the expected finite-volume shift, is a prerequisite for the claim that the plateau in Fig. 4 (right) is a cutoff-independent physical mass rather than a finite-volume artifact.","section":"Section 4, Table 1 and Eq. (3)"},{"comment":"The LPCP fixes only R and g_GF, while the BSM bare couplings are held fixed, and the paper assumes that the BSM sector is 'mostly insensitive' to SM physics. The evidence shown is a coarse monitoring of R and S under the κ1 scan. Table 1 shows that R varies from 1.462(53) to 1.527(80), so the line is not tuned exactly to R = 1.5. If m_H/m_W depends on the residual deviations from the LPCP or on R itself, the plateau could be an artifact of an incompletely tuned constant-physics line. The authors should quantify the insensitivity, for example by varying η3 or R at one β and demonstrating that m_H/m_W shifts by less than the quoted statistical error.","section":"Section 3 and Section 4, Fig. 4 (left)"}],"minor_comments":[{"comment":"Several typographical errors should be corrected: 'bellow' in the abstract should be 'below', 'cutodial' in the Fig. 1 caption should be 'custodial', 'scaning' in the conclusion should be 'scanning', and 'LCPC' in the discussion of Fig. 4 should be 'LPCP'.","section":"Abstract; Fig. 1; Conclusion"},{"comment":"Reference [12] has a formatting error: '1604 [1410.2740]' should be '1604 (2016) [arXiv:1410.2740]' or similar.","section":"References"},{"comment":"The left panel caption should explicitly define which quantities are plotted (R and S) and explain the marker types and the meaning of the shaded bands, so that the figure is self-contained.","section":"Fig. 4 caption"},{"comment":"The sentence reporting m_screen ≈ 50 GeV should carry a caveat in the main text (not only in footnote 3) that the screening mass depends on the assumed Yukawa form and on the scheme, since the text later uses this value without further qualification.","section":"Section 4, screening mass discussion"},{"comment":"The notation 2μ^2 Tr(Φ1†Φ2) for the bare mass-mixing term is easily confused with the renormalization scale μ used later; renaming this bare parameter, e.g. as μ12^2, and stating its relation to the parameters in Eq. (1) would improve clarity.","section":"Section 2, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style contribution whose central quantitative claim is bolder than the supporting evidence. The main missing pieces are error bars on the BSM mass ratios and a finite-volume check; both are local additions rather than conceptual changes, so I would not reject the manuscript. I recommend that the editor require the authors to either provide these items or temper the abstract and conclusion to report the result as an exploratory observation with the current systematic uncertainties stated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look. This proceedings does something no lattice 2HDM paper has done: builds a line of constant SM physics (R = m_h/m_W ~ 1.5 and g^2_GF(m_W) = 0.5), measures gradient-flow running, and extracts BSM mass ratios across cutoffs from 300 to 630 GeV. The central observation is a plateau in m_H/m_W ~ 0.2 as kappa1 approaches the H12 transition, while m_A/m_W keeps falling. If true, a custodial/inert 2HDM with weak quartics can realize an extra neutral scalar near 16 GeV. That is a concrete, falsifiable map of parameter space, and the paper is honest that it is a proceedings with limited scope.\n\nCredit where due: the LPCP construction is a genuine technical step; the gradient-flow running collapses onto one curve with a perturbative match at high energy and a Yukawa/screening fit in the IR; and the finite-T crossover argument, based on absence of volume dependence in the susceptibility peak, is reasonable even if limited to the two finest lattices. The authors also flag their own limitations: small BSM couplings are a choice, not a necessity, and only two beta values give a clean susceptibility peak.\n\nThe soft spots are real. First, the headline number has no quoted uncertainty anywhere; Fig. 4 clearly has errors, but the abstract and conclusion state m_H ~ 0.2 m_W as a result. A lattice number without an error bar is not yet a result. Second, the LPCP tuning is approximate: R ranges 1.46(5) to 1.53(8), a few sigma off nominal, and the BSM mass ratios could shift if the SM line were retuned more tightly. Third and most important, the stress-test concern lands: from Table 1, m_H L ranges from about 0.94 to 1.7. The lightest BSM state has a correlation length comparable to the box size at the finest point (beta = 8.56). Without a volume test, the apparent plateau in m_H/m_W could be the finite-volume floor: the extracted ground-state mass cannot fall much below 1/L. The observation that m_A/m_W keeps decreasing while m_H/m_W flattens is exactly what one would see if the lighter H hits the floor first. So the cutoff-independence claim is not established; it is an invitation to run a volume study. Finally, no code or data are released, so the numerics cannot be checked from the manuscript.\n\nNone of this is disqualifying for a proceedings. All the issues are addressable in a full paper: quote errors, document the correlator fits, add a volume scan of the BSM spectrum, and publish the data. The paper deserves a serious referee, but the referee should insist on the finite-volume check before the 0.2 m_W claim is accepted. I would not cite it as a result yet, but I would follow the full paper when it appears.","headline":"First lattice line-of-constant-physics study of the inert/custodial 2HDM, with a striking low-mass claim that needs a finite-volume check before I'd trust it.","tokens_in":9920,"tokens_out":2021,"would_cite":false,"duration_ms":20479,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T25","81T13"],"pacs":["11.15.Ha","12.60.Fr"],"model":"deepseek-v4-flash","headline":"On a line of constant Standard-Model physics, the custodial two-Higgs-doublet model with weak quartic couplings allows an extra neutral scalar H with mass about 0.2 times the W boson mass, and this ratio does not change as the lattice…","keywords":["two Higgs doublet model","lattice gauge theory","custodial symmetry","inert doublet model","line of constant physics","gradient flow coupling","BSM scalar spectrum","finite temperature phase transition"],"falsifier":"On the same five lattice spacings, re-tune $\\kappa_2$ and $\\eta_2$ for each $\\kappa_1$ value so that $R$ and $S$ match the target values to much higher precision, then check whether $m_H/m_W$ still approaches $0.2$ and remains cutoff independent. If the plateau shifts or disappears, the reported lower bound is an artifact of the LPCP tuning rather than a property of the model with constant SM physics.","tokens_in":8695,"feed_emoji":"⚛️","tokens_out":5311,"duration_ms":46265,"temperature":0.7,"pith_summary":"The paper asks whether the custodial (inert) two-Higgs-doublet model can hide additional scalars far below the W boson mass while still reproducing Standard-Model physics. Using lattice simulations on a line of constant SM physics with cutoff from 300 to 630 GeV, the authors report that, for the small quartic couplings they choose, the neutral BSM state H reaches a mass of about $0.2\\,m_W$ (near 16 GeV), and this ratio stays flat as the cutoff changes. They also compute the running of the weak gauge coupling and the finite-temperature transition, finding a smooth crossover at weak couplings. A sympathetic reader would take the result as a non-perturbative, cutoff-independent lower-bound estimate for a realizable H mass in this regime, and as a benchmark for later scans at larger couplings.","feed_headline":"Extra Higgs scalar can be as light as 16 GeV on the lattice","feed_subtitle":"In a custodial two-Higgs-doublet model, the H mass ratio plateaus at 0.2 while A and H± keep falling toward the W.","key_machinery":"The argument is carried by a line of partially constant physics (LPCP) in the bare parameter space of the lattice action. The SM sector is tuned at each bare gauge coupling $\\beta$ by adjusting $\\kappa_2,\\eta_2$ so that $R\\equiv m_h/m_W\\approx1.5$ and the gradient-flow renormalized gauge coupling satisfies $g^2_{GF}(\\mu=m_W)=0.5$ (equivalently $\\sqrt{8t_0}\\,m_W=1.0$); all BSM couplings $\\kappa_1,\\eta_1,\\eta_3,\\eta_4,\\eta_5$ are held fixed at small values. The spectrum is read from two-point functions of composite operators that separate the SM Higgs and W boson from the BSM states $H,A,H^\\pm$. The BSM scan then varies $\\kappa_1$ within the $\\mathbf{H}_2$ phase while monitoring $R$ and $S$ to check that SM physics is unchanged.","core_discovery":"The central discovery is that the custodial inert two-Higgs-doublet model on the lattice admits a BSM neutral scalar whose mass ratio $m_H/m_W$ saturates at roughly $0.2$ before the transition into the phase where both doublets condense, and this saturation is independent of the lattice cutoff over $300$--$630$ GeV. Along the same line of constant SM physics the heavier states $A$ and $H^\\pm$ behave differently: $m_A/m_W = m_{H^\\pm}/m_W$ keeps decreasing as the hopping parameter $\\kappa_1$ is raised, so the mass gap between $H$ and $A/H^\\pm$ widens. Because the SM conditions $R=m_h/m_W\\approx1.5$ and $g^2_{GF}(m_W)=0.5$ are fixed at each cutoff, the light $H$ is a property of a theory that still looks like the SM in the Higgs-gauge sector.","pith_inferences":[],"forward_implications":["A neutral inert-doublet scalar $H$ with $m_H\\simeq 0.2\\,m_W$ can coexist with SM-like Higgs and W masses in a weakly coupled custodial 2HDM.","The cutoff independence of the ratio justifies calling $\\sim0.2\\,m_W$ a non-perturbative lower bound of the H mass for the chosen couplings.","$m_A/m_W=m_{H^\\pm}/m_W$ keeps falling toward the W as $\\kappa_1$ approaches the $\\mathbf{H}_{12}$ transition, so the H--A mass splitting grows along the LPCP.","The gradient-flow running gauge coupling collapses onto one curve over a wide energy range, matching one-loop massless running at high energy and a screened Yukawa form at low energy.","At small quartic couplings the electroweak transition is a crossover, with no volume dependence in the susceptibility peak; the transition sits near $m_W/T_c\\sim0.5$.","If the $\\sim0.2\\,m_W$ plateau survives a fully retuned LPCP, the model predicts a $\\sim16$ GeV neutral scalar that lies within reach of low-energy collider and Higgs-decay searches, a target the paper does not discuss.","A natural next step is to repeat the $\\kappa_1$ scan with larger $\\eta_3$ or $\\eta_4$, since the paper notes these couplings must be $O(1)$ for a strong first-order transition; the light-H plateau may move or disappear in that regime, which would change the phenomenology.","The partially constant physics line fixes only SM quantities; an extension that also matches a BSM renormalized quantity would test whether the H-mass bound is an accident of the tuning strategy."],"supporting_citations":[{"why":"Defines the inert and custodial 2HDM limits and the four phases H0, H1, H2, H12 used to place the lattice parameters.","marker":"[10]"},{"why":"Supplies the Wilson flow used to define the flowed action density and hence the renormalized gauge coupling.","marker":"[19]"},{"why":"Gives the perturbative relation between the flowed action density and the running gauge coupling used in Eq. (5).","marker":"[20]"},{"why":"Provides the physical W boson mass used to set the lattice scale and convert lattice masses to GeV.","marker":"[21]"},{"why":"Provides the SM crossover result used as reference for interpreting the finite-temperature transition in the 2HDM.","marker":"[2]"}],"fun_headline_variants":["Lattice study: extra Higgs can be as light as 16 GeV","Lightest extra Higgs in 2HDM lattice: 0.2 m_W, cutoff-independent","Two-Higgs-doublet lattice: new scalar dips to 0.2 m_W","16 GeV extra Higgs state emerges in lattice 2HDM","Cutoff-independent light H in custodial 2HDM on lattice"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction rests on the assumption that the BSM sector is 'mostly insensitive' to the SM sector, so that scanning $\\kappa_1$ while keeping all other BSM couplings fixed leaves the line of constant SM physics intact; the paper verifies only that $R$ and $S$ stay roughly constant, not that small mistunings would leave the BSM mass ratios unchanged.","fun_headline_variants_meta":{"raw":{"variants":["Lattice study: extra Higgs can be as light as 16 GeV","Lightest extra Higgs in 2HDM lattice: 0.2 m_W, cutoff-independent","Two-Higgs-doublet lattice: new scalar dips to 0.2 m_W","16 GeV extra Higgs state emerges in lattice 2HDM","Cutoff-independent light H in custodial 2HDM on lattice"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001197,"raw_usage":{"total_tokens":4918,"prompt_tokens":907,"completion_tokens":4011,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":3906}},"tokens_in":523,"tokens_out":4011,"duration_ms":28417,"temperature":1.0,"reasoning_tokens":3906,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:41:11.114466+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On the same five lattice spacings, re-tune $\\kappa_2$ and $\\eta_2$ for each $\\kappa_1$ value so that $R$ and $S$ match the target values to much higher precision, then check whether $m_H/m_W$ still approaches $0.2$ and remains cutoff independent. If the plateau shifts or disappears, the reported lower bound is an artifact of the LPCP tuning rather than a property of the model with constant SM physics.","supporting_citations":[{"cited_title":"Branco, P","cited_arxiv_id":null,"evidence_quote":"Defines the inert and custodial 2HDM limits and the four phases H0, H1, H2, H12 used to place the lattice parameters."},{"cited_title":"Lüscher,Properties and uses of the wilson flow in lattice qcd,Journal of High Energy Physics 2010 (2010)","cited_arxiv_id":null,"evidence_quote":"Supplies the Wilson flow used to define the flowed action density and hence the renormalized gauge coupling."},{"cited_title":"Lüscher and P","cited_arxiv_id":null,"evidence_quote":"Gives the perturbative relation between the flowed action density and the running gauge coupling used in Eq. (5)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the physical W boson mass used to set the lattice scale and convert lattice masses to GeV."},{"cited_title":"D’Onofrio and K","cited_arxiv_id":null,"evidence_quote":"Provides the SM crossover result used as reference for interpreting the finite-temperature transition in the 2HDM."}],"review_version":1}