{"id":"78ae7ecd-72ae-406e-afff-96c2e4eedc8d","arxiv_id":"hep-lat/9412091","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A 3d effective-theory framework, with 2-loop gauge-invariant condensates and the lattice-to-continuum conversion constants, for determining electroweak phase transition parameters by Monte Carlo.","lead":"The authors lay out a two-step recipe for pinning down the electroweak phase transition: reduce the hot 4d Standard Model to an effective 3d theory analytically, then simulate that 3d theory on the lattice. They supply the conversion constants that translate lattice numbers into physical temperature and Higgs mass, plus a set of gauge-invariant order parameters measurable in Monte Carlo.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The rho calibration is the load-bearing soft spot: a five-parameter fit with an unconstrained 4-loop coefficient and a 2-sigma tension between the direct (2.18(6)) and A0-integrated (2.03(5)) routes means the 0.3% temperature-accuracy claim is not demonstrated.","rationale":"The reader correctly identifies the rho extraction as the weakest link. The analytic content (nu_rho, omega_rho) is solid, and the framework is well-motivated. The concern is not that the framework is wrong, but that the numerical calibration of rho is not demonstrated to the claimed accuracy. The 2-sigma tension and the unconstrained psi-tilde are specific, acknowledged red flags. The paper is honest about them, but honesty about a limitation does not remove it. The verdict CONDITIONAL is appropriate. I would not adjust it: the paper's own estimate that rho is known to 0.12 at 2 sigma gives 0.3% accuracy in temperature, but the 2-sigma tension between 2.18(6) and 2.03(5) is 0.15, already larger than the quoted statistical error, so the systematic reach is less well controlled than the statistical reach. The concrete test above would settle whether the flat directions in the fit are biasing rho or whether the tension is a real discrepancy.","tokens_in":59014,"tokens_out":9539,"duration_ms":173623,"concrete_test":"Re-fit the published MC data with (i) an added cubic finite-a term C x^3 and (ii) the coarsest lattice (beta_G = 12) excluded, holding psi-tilde fixed at -700, 0, and +700. If rho shifts by more than the quoted statistical error 0.06 under any of these variations, the flat directions are not controlled and the 0.3% claim is unsupported. Alternatively, compute the 2-loop Debye mass constants c1, c2 (eq. 30) analytically and check whether including them resolves the 2.18 vs 2.03 tension.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the constant physics curves (eqs. 91, 92) and condensate relations (eqs. 138, 140) are exact in the continuum limit, with rho, nu_rho, omega_rho determined well enough to give T_c to 0.3%. nu_rho and omega_rho are computed analytically and appear solid. rho, however, is extracted from a five-parameter fit (rho, A, B, phi-tilde, psi-tilde) to the scalar condensate at beta_G = 12, 20, 24, 32, 40. The fit assumes a specific form for the finite-a corrections (A x + B x^2, eq. 146) and a specific form for the 3- and 4-loop continuum corrections (eq. 147). The paper reports psi-tilde unconstrained over -700 < psi-tilde < 700, meaning the 4-loop correction is essentially a free function that could absorb residual finite-a or higher-order effects. The independent cross-check via the A0-integrated-out theory gives rho = 2.03(5), about 2 sigma below the direct value 2.18(6). The authors acknowledge this tension and attribute it possibly to the uncomputed 2-loop Debye mass constants (c1, c2 in eq. 30), which are bounded by an assumed inequality |zeta_g + ...| < 14 (eq. 156) rather than evaluated. The 'exact in the continuum limit' claim is therefore not established for rho: the continuum extrapolation is a fit with a free 4-loop parameter, and the two routes disagree at the level of the quoted uncertainty. This is load-bearing because rho is the single number that calibrates the lattice-to-continuum map for the Higgs mass relation, and the 0.3% temperature accuracy claim rests on it.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"This is a technical framework paper for lattice Monte Carlo study of the electroweak phase transition via a 3d dimensionally reduced effective theory. The authors (i) motivate dimensional reduction and estimate its accuracy and its computational advantage over 4d lattice simulations (§2), (ii) summarize the 2-loop 3d effective potential and the 4d→3d matching (§3), (iii) introduce gauge-invariant condensates defined as derivatives of the renormalized vacuum energy with respect to the 3d parameters, and give their 2-loop MSbar values (§4–5), (iv) write down the lattice action and derive the constant-physics curves relating lattice couplings to 3d continuum parameters, which are exact in the continuum limit up to three pure numbers rho, nu_rho, omega_rho (§6), and (v) compute nu_rho and omega_rho analytically in lattice perturbation theory and extract rho from a five-parameter fit to Monte Carlo data for the scalar condensate (§7), with a cross-check from the A0-integrated-out theory (§8). The quoted results are rho = 2.18(6), nu_rho = 1.01, omega_rho = 0.44, with a 2-sigma tension against the A0-integrated route rho = 2.03(5).","tokens_in":59354,"tokens_out":6321,"duration_ms":126953,"significance":"The paper's structural results are genuinely valuable and are stated with unusual care. The 2-loop effective potential (Eq. 34) and the sunset function (Eq. 36) are explicit; the lattice tadpole, one-loop vacuum energy, and sunset expansions (Eqs. 169–174) are given with numerical coefficients; and the two non-trivial 2-loop integrals delta = 1.94 (Eq. 175) and zeta = -0.314 (Eq. 176) are written as definite integrals with numerical values, which makes them independently checkable. The derivation of omega_rho in pure scalar theory (Eqs. 110–117) is a clean step-by-step computation. The claim that the constant-physics curves have no higher-loop contributions follows from 3d super-renormalizability and is correct as a statement about the structure of the matching. The paper also ships an explicit erratum to its own published Letter (the correction of rho and nu_rho in ref. [4], §6.3), which is the kind of self-correction that should be credited. The framework is falsifiable: the fitted rho must be Higgs-mass independent, and the paper tests this at m_H = 160 GeV (Eq. 149).","major_comments":[{"comment":"The extraction of rho = 2.18(6) rests on a five-parameter fit (rho, A, B, phi-tilde, psi-tilde) to the scalar condensate at beta_G = 12, 20, 24, 32, 40, with the finite-a correction assumed to have the form A x + B x^2 (Eq. 146) and the 3- and 4-loop continuum corrections assumed to have the form of Eq. (147). The paper itself reports that psi-tilde is unconstrained over -700 < psi-tilde < 700, i.e. the 4-loop coefficient is effectively a free function that can absorb residual finite-a and higher-order effects. The quoted error 2.18(6) is statistical only; the systematic from truncating f(x) at two terms is not quantified, and the paper offers no third-term stability check. This is load-bearing because rho is the single number that calibrates the lattice-to-continuum map for the Higgs mass relation (Eq. 99), and the 0.3% temperature-accuracy claim in §8 depends on it. The authors should either add a third term to f(x) and quote the shift in rho, or state explicitly that the quoted error is statistical and that the systematic from the fit form is not bounded.","section":"§7, Eq. (148) and surrounding text","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The independent cross-check via the A0-integrated-out theory gives rho = 2.03(5), which the authors state is about 2 standard deviations below the direct value 2.18(6). They attribute the discrepancy to the uncomputed 2-loop Debye mass constants c1, c2 in Eq. (30), bounded by the assumed inequality |zeta_g + (lambda_3/g_3^2) zeta_{g lambda}| < 14 (Eq. 156). That bound is asserted, not derived, and the paper does not show where the number 14 comes from. If the true value sits at the edge of the bound, the discrepancy disappears but so does the 0.3% accuracy claim; if it sits near zero, the two routes genuinely disagree at the quoted precision. The paper should either derive the bound or present the 0.3% figure as conditional on it, and should state what value of the combination would be required to reconcile 2.18(6) with 2.03(5).","section":"§8, Eq. (165) versus Eq. (148)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract (abstract:s3) states that the relations between lattice couplings and continuum parameters are 'exact in the continuum limit', and §6.2 says the constant-physics curves 'do not have any higher loop contributions'. Both statements are correct for the structure of the matching, but the numerical implementation involves truncations that are numerically, not structurally, motivated. Footnote 14 states that f_2D is formally multiplied by log^2(aT), which is singular in the continuum limit, and is dropped because 'in real lattice simulations with beta_G < 40, the value of this log never exceeds 2'. Similarly P_2D = 0 and the O(g^4) terms in Eq. (30) are dropped by power counting, and the paper itself notes that 'the formal suppression of the higher-order terms does not necessarily imply numerical suppression'. The 'exact' language should be qualified to 'exact up to the stated truncations', or the truncations should be shown to vanish in the continuum limit. As written, a reader could take 'exact' to cover the numerical values of rho, nu_rho, omega_rho, which it does not.","section":"§6.3, footnote 14, and the abstract claim of exactness","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The text of §6.2 quotes sigma = 3.17159, while Appendix A, Eq. (170) gives sigma = 3.1759114. These should be reconciled; the value enters Eq. (91) and hence the constant-physics curve.","section":"§6.2 versus Appendix A, Eq. (170)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The symbol rho is used both for the constant-physics parameter rho = 2.18(6) (Eq. 102) and for the lattice tadpole coefficient rho = 0.15281 (Eq. 169). A different symbol for one of the two would help the reader, especially since both appear in the same section.","section":"Notation, Eq. (102) versus Eq. (169)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The erratum corrects ref. [4] to rho = 2.0 and nu_rho = 0.66, whereas this paper quotes rho = 2.18(6) and nu_rho = 1.01. The reader would benefit from a sentence explaining whether the difference is due to improved fits, a different data set, or a different truncation.","section":"§6.3, erratum to ref. [4]","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The vertical axis is labelled 'Delta R^2 beta_H beta_G / x' and the horizontal axis 'x = phi/(g T beta_G)'. Since the collapse of points at different beta_G onto a single curve is the visual evidence for the fit, the caption should state explicitly what normalization is applied and what the units of the vertical axis are.","section":"Figs. 5 and 6","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The estimates of '10^3 larger 3d volumes' and '3 to 4 orders of magnitude gain in computer time' rest on assumed values m/T ~ 0.2 and p = 10%. These assumptions should be stated in the text at the point where the numbers are quoted, not only in the surrounding prose.","section":"§2.2.2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The numerical constants rho = 0.15281, zeta = -0.6, delta = 1.94, and zeta = -0.314 are quoted without an indication of the integration method or the numerical precision. Since these enter the central results, a sentence on how they were obtained and to what accuracy would be appropriate.","section":"Appendix A","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract's 'will be able to reliably determine' and the 0.3% temperature-accuracy figure in §8 are projections; the actual determination of T_c is deferred to ref. [3], listed as 'in preparation'. This is appropriate for a framework paper, but the abstract should make clear that the 0.3% figure is an estimate of the accuracy the framework will deliver, not a demonstrated result of this paper.","section":"Abstract and §9","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a careful framework paper and the structural claims are sound. The soft spot is the rho extraction: the five-parameter fit has a flat direction (psi-tilde unconstrained over a factor of 1400) and the two independent routes to rho disagree at 2 sigma. The authors are transparent about both, which is to their credit, but the 'exact in the continuum limit' phrasing in the abstract and the 0.3% accuracy claim in §8 should be qualified. The fixes are local (qualification of claims, derivation or explicit statement of the bound in Eq. 156, a third-term stability check on f(x)), so minor revision seems right rather than major. The paper is otherwise well suited to the journal."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"This is the infrastructure paper for the 3d lattice electroweak program, and it does the unglamorous work well. The genuinely new pieces are the gauge-invariant condensates defined as derivatives of the vacuum energy with respect to the five theory parameters (section 4), and their 2-loop evaluation in the broken phase (eqs. 71, 75-78). That construction is clean and it is the right one: it sidesteps the gauge-dependence of the effective potential without hand-waving. The lattice perturbation theory is also real work, not decoration. omega_rho = 0.44 comes out of pure scalar theory (eqs. 110-117), nu_rho = 1.01 out of the SSV diagrams, and both reduce to definite integrals (eqs. 175-176) with quoted numerical values. Anyone can check those. The paper also corrects its own earlier published numbers from ref [4], which is the kind of thing that buys credibility rather than costing it. The distinction between the linear and logarithmic counterterms, and the demonstration that the 1/m^2 poles cancel between one- and two-loop pieces, are handled properly. Where I part company with the harshest reading: the 'exact in the continuum limit' claim is about the form of eqs. 91, 92, 138, 140, and that form is exact by super-renormalizability. rho is a scheme-conversion constant that has to be fixed nonperturbatively; the paper says so. The soft spot is narrower and more specific: the numerical value of rho. It comes from a five-parameter fit to the condensate, with A and B absorbing finite-a corrections and psi-tilde unconstrained over -700 to 700. Nothing in the data separates rho from those except curve collapse across beta_G. The A0-integrated cross-check giving 2.03(5) against 2.18(6) is the honest internal signal, and the authors report it rather than bury it. The systematic from the uncomputed Debye-mass constants is bounded by an assumed inequality, not evaluated. So the 0.3% temperature figure is a statistical reach, not a demonstrated accuracy. Read it as a framework with a calibrated constant, not a finished prediction. Worth a serious referee, and worth citing if you use the condensate relations.","headline":"A solid, honest methods paper: the condensate construction and the analytic lattice perturbation theory hold up, but rho is a calibration constant, not a computation, and the 0.3% temperature claim leans on it harder than the data supports.","tokens_in":60081,"tokens_out":970,"would_cite":true,"duration_ms":18160,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["11.15.Ha","11.10.Wx","98.80.Cq"],"model":"selfhosted/deepseek-v4.1-flash","headline":"The paper argues that the electroweak phase transition can be computed reliably by simulating a three-dimensional effective theory on the lattice, because the map between lattice couplings and continuum couplings — the constant physics曲线 —,","keywords":["dimensional reduction","electroweak phase transition","3d effective theory","lattice Monte Carlo","constant physics curves","gauge-invariant condensates","lattice perturbation theory","Debye screening"],"falsifier":"An independent determination of $rho$ would settle it — either a direct lattice perturbation theory computation of the three- and four-gluon contributions to the two-loop Higgs mass operator, or substantially better statistics that separate $rho$ from $A$ and $B$; separately, pinning down the two unknown constants in the two-loop Debye mass would show whether the $2.18(6)$ versus $2.03(5)$ discrepancy is statistical or a real systematic.","tokens_in":60293,"feed_emoji":"⚛️","tokens_out":3890,"duration_ms":80173,"temperature":0.7,"pith_summary":"The paper asks how to compute, from first principles, the numbers that characterise the electroweak phase transition: the temperature at which it happens, the size of the jump in the order parameter, and the latent heat. Its answer is a two-stage strategy. First, replace the full four-dimensional finite-temperature Standard Model with a three-dimensional effective theory in which the heavy thermal modes have been integrated out analytically; second, simulate that 3d theory on the lattice, where the strongly coupled infrared physics is handled numerically. The bridge between the two stages is a set of relations — the constant physics curves — connecting lattice couplings to continuum couplings, and the paper's claim is that these relations, together with the relations between lattice and $overline{mathrm{MS}}$ condensates, are exact in the continuum limit. If that holds, a lattice simulation stops being a model of a model and becomes a numerically controlled statement about the hot Standard Model, with the residual error budget set by ordinary perturbation theory rather than by the lattice.","feed_headline":"Three pure numbers link lattice to hot electroweak theory","feed_subtitle":"Dimensional reduction plus two-loop condensates pin the transition temperature to roughly 0.3 percent.","key_machinery":"The load-bearing mechanism is super-renormalizability of the 3d theory: because only one- and two-loop graphs diverge, the difference between lattice and continuum regularization can only be a linear counterterm, a logarithmic counterterm, and a constant. That reduces the entire matching problem to three numbers multiplying $g_3^4$, $lambda_3 g_3^2$ and $lambda_3^2$ in the mass polynomial, which is what makes an exact-in-the-continuum-limit relation possible at all. The second piece of machinery is the gauge-invariant condensate: differentiating the finite vacuum energy density with respect to $m_3^2$, $m_D^2$, $lambda_3$, $lambda_A$ and $g_3^2$ produces order parameters that are computable at two loops and measurable on the lattice with small statistical error, sidestepping the gauge dependence of the effective potential itself.","core_discovery":"The central claim is that the 3d gauge-Higgs theory is super-renormalizable, so its ultraviolet divergences occur only at one and two loops and all scheme dependence collapses into a small set of constants. Consequently the lattice-to-continuum map can be written exactly in the continuum limit and is parametrized by just three pure numbers, denoted $rho$, $nu_rho$ and $omega_rho$. The authors compute $nu_rho=1.01$ and $omega_rho=0.44$ analytically in lattice perturbation theory, and extract $rho=2.18(6)$ by matching Monte Carlo measurements of the scalar condensate $langle R_L^2 rangle$ against the two-loop continuum prediction, fitting simultaneously for the finite-lattice-spacing coefficients $A=5.09(25)$, $B=-3.40(23)$ and a higher-order coefficient $tilde{phi}=-162(22)$. They also define gauge-invariant condensates as derivatives of the renormalized vacuum energy density with respect to the couplings, so that the latent heat of a first-order transition is just the jump of the scalar condensate between phases. As a cross-check they integrate out the $A_0$ field, obtain $rho'=1.62(7)$ from a separate simulation, and convert it analytically to $rho=2.03(5)$ — a value the authors themselves note is about two standard deviations below the direct determination.","pith_inferences":["The roughly two-sigma gap between $rho=2.18(6)$ and the $A_0$-eliminated value $2.03(5)$ is the most informative number in the paper for anyone judging the programme: if it survives better statistics, it points to the neglected $g^4$ contributions to the Debye screening mass being numerically larger than power counting suggests, which would be a cheap diagnostic for the whole dimensional-reduction","Because finite-spacing effects in 3d start at order $a$ rather than $a^2$, the practical reach of this method is set by how tightly the linear coefficient $A$ can be pinned — a generic feature of simulating any super-renormalizable 3d theory, not something specific to the electroweak case.","The trick of defining order parameters as derivatives of the vacuum energy with respect to couplings is portable: any super-renormalizable 3d theory one wishes to simulate, hot QCD somewhat above its transition being the obvious neighbour, admits the same construction.","The framework buys precision only on the broken side; the paper concedes that the symmetric phase is strongly coupled and inaccessible to perturbation theory, so the most cosmologically interesting regime is precisely where this machinery hands off rather than delivers."],"forward_implications":["The transition temperature can be located to roughly 0.3 percent, which at $m_H=80$ GeV means an uncertainty of about 0.5 GeV — fine enough to decide whether electroweak baryogenesis survives for a given Higgs mass.","The latent heat and the order-parameter jump follow directly from the difference of the scalar condensate between the broken and symmetric phases, with no need for a four-dimensional lattice simulation.","A 3d lattice with $beta_G$ between 12 and 40 reaches the same finite-spacing accuracy as a 4d simulation at a cost roughly three to four orders of magnitude smaller in computer time.","Comparing measured condensates with the two-loop prediction yields an empirical estimate of the neglected three-loop contribution, quoted as $phi=-49(22)$, which is a direct check on how well perturbation theory works in the broken phase.","The same constant physics curves can be written for the theory with the $A_0$ field integrated out, giving an independent route to the same physics and a consistency test of the matching procedure."],"supporting_citations":[{"why":"Supplies the two-loop dimensional reduction, the 3d effective Lagrangian, the two-loop effective potential and the mass counterterms on which every subsequent computation in this paper rests.","marker":"[2]"},{"why":"The earlier paper that initiated the programme and produced the first qualitative lattice results for the light-Higgs case, establishing the setting this framework formalizes.","marker":"[1]"},{"why":"Introduces the gauge-invariant effective potential for $langle phi^dagger phi rangle$ that Appendix B clarifies and extends to two loops here.","marker":"[13]"},{"why":"Provides the prototype for relating a lattice coupling to a continuum coupling scheme, the model the constant-physics-curve derivation follows.","marker":"[11]"},{"why":"The continuum side of that same lattice-to-continuum coupling relation, used as the template for how the matching constants are defined.","marker":"[12]"},{"why":"Establishes that the $1/m^2$ poles cancel between one- and two-loop contributions for any gauge-fixing parameter, which underwrites the claim that the extracted condensate is gauge invariant.","marker":"[27]"},{"why":"Numerical study of the pure-gauge triplet sector, supporting the assertion that the $A_0$ field develops no vacuum expectation value and can be treated as dynamically inessential.","marker":"[29]"},{"why":"Four-dimensional lattice work with a temperature-dependent subtraction scheme, the alternative route against which the efficiency argument for 3d simulations is measured.","marker":"[9]"}],"fun_headline_variants":["Three pure numbers tie lattice couplings to electroweak theory","Lattice-to-continuum map needs just three pure constants","Two-loop condensates pin the electroweak lattice map","Three constants map the lattice; cross-check sits two sigma off","Super-renormalizable 3d theory reduces scheme dependence to three numbers"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The number $rho$ is not measured directly: it is inferred from a five-parameter fit in which the finite-lattice-spacing coefficients $A$ and $B$ and an unconstrained higher-order coefficient float alongside it, so $rho$ is identified only by requiring data at different lattice spacings to collapse onto a single curve.","fun_headline_variants_meta":{"raw":{"variants":["Three pure numbers tie lattice couplings to electroweak theory","Lattice-to-continuum map needs just three pure constants","Two-loop condensates pin the electroweak lattice map","Three constants map the lattice; cross-check sits two sigma off","Super-renormalizable 3d theory reduces scheme dependence to three numbers"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1407,"prompt_tokens":935,"completion_tokens":472,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":400}},"tokens_in":935,"tokens_out":472,"duration_ms":10961,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":400,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-11T13:17:36.083028+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"An independent determination of $rho$ would settle it — either a direct lattice perturbation theory computation of the three- and four-gluon contributions to the two-loop Higgs mass operator, or substantially better statistics that separate $rho$ from $A$ and $B$; separately, pinning down the two unknown constants in the two-loop Debye mass would show whether the $2.18(6)$ versus $2.03(5)$ discrepancy is statistical or a real systematic.","supporting_citations":[],"review_version":1}