{"id":"b39fa79b-14df-4abe-8a4e-ae7cb0d7f2c3","arxiv_id":"2607.21346","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For generalized domain-wall fermions, the O(e^2) electromagnetic expansion requires new local seagull and anti-quark contact vertices, derived here for the first time.","lead":"Lattice QCD calculations of electromagnetic corrections to quark processes need special 'seagull' terms to remain gauge-consistent; this paper derives those missing terms for generalized domain-wall fermions. The result gives lattice groups a complete O(e^2) vertex set for computing radiative corrections to hadronic decays and precision Standard Model tests.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sign error in Eq. (90): the anti-quark contact term \\bar Q_2 enters with the wrong sign; Eq. (85) implies +1/2⟨q\\bar Q_2⟩, not -1/2.","rationale":"The paper’s core construction—the derivative expansion of the generalized DWF propagator in a background QED field—is coherent and the algebraic steps leading to Eq. (85) check out. The finite-difference validation in Sec. VII provides independent support for Eq. (85). However, the central presentation of the result, Eq. (90), contains a sign error in the \\bar Q_2 contact term. The reader’s identified weakness (the convention for the physical anti-quark field) is not the most load-bearing issue: the paper derives qbar from the overlap/DWF relation rather than assuming it arbitrarily, and different field normalizations would not invalidate the derivation. The sign error is a concrete, localized flaw in the stated main result. Since Eq. (85) is correct and the error is likely typographical rather than conceptual, rejection is too strong; the paper should be accepted only after Eq. (90) is corrected (and the missing repository URL in Ref. [34] is supplied). The verdict therefore remains CONDITIONAL/UNCHANGED.","tokens_in":15159,"tokens_out":53115,"duration_ms":510162,"concrete_test":"Use the released Grid implementation on the 8^4 single-configuration test of Fig. 1. Compute the finite-difference second derivative δ2φ of the physical 4D propagator (Eq. 131) and separately evaluate Eq. (90) with the \\bar Q_2 term as printed (−) and with the sign corrected to +. The version that matches δ2φ is correct; Eq. (85) predicts the + sign. Alternatively, independently re-derive Eq. (90) from Eq. (85) by substituting S = S_dw D_- and collecting the coefficient of S_dw Γ_2 C; it is -1/2, which equals +1/2⟨q\\bar Q_2⟩ because of the minus in Eq. (89).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s main operator-insertion formula, Eq. (90), is internally inconsistent with Eq. (85). Starting from Eq. (85) and the paper’s own definitions \\bar Q_1^\\mu = -\\bar\\psi \\Gamma_1^\\mu C R_5 P (Eq. 46) and \\bar Q_2^\\mu = -\\bar\\psi \\Gamma_2^\\mu C R_5 P (Eq. 89), the O(e^2) physical propagator is\n1/2 ∂²_e S_q = -⟨VV⟩ + ⟨V\\bar Q_1⟩ - 1/2⟨T⟩ + 1/2⟨\\bar Q_2⟩.\nThe key term is -1/2 S_dw Γ_2 A^2 C in Eq. (85). Since \\bar Q_2 contains an explicit minus sign, S_dw Γ_2 C = -⟨q\\bar Q_2⟩, so -1/2 S_dw Γ_2 C = +1/2⟨q\\bar Q_2⟩. Equation (90), however, writes the last bracket as -1/2(⟨T⟩+⟨\\bar Q_2⟩), giving -1/2⟨\\bar Q_2⟩. The sign of the anti-quark contact term is therefore opposite to what Eq. (85) dictates. This is not a convention ambiguity: it follows from the definitions stated in the paper. The numerical validation in Sec. VII tests Eq. (85), not Eq. (90), so the discrepancy is not caught by the reported finite-difference checks. Any implementation based on Eq. (90) would insert the O(α_em) contact interaction with the wrong sign.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the perturbative expansion of the physical quark propagator in generalized domain-wall QCD+QED with respect to the electromagnetic charge e up to O(e^2). Starting from the U(3)-link Wilson operator, the authors obtain the first-order conserved vector current and, at second order, a new local seagull current T^mu plus anti-quark contact operators Qbar1^mu and Qbar2^mu, the latter arising because the physical anti-quark field qbar depends on the QCD+QED links through the operator D_- in Eq. (26). The derivation is performed with a background-field method, and the resulting Ward-Takahashi identity for the DWF currents is used to justify the seagull vertex. The authors also provide a numerical finite-difference validation on an 8^4 lattice for Shamir, scaled-Shamir, and zMobius domain-wall actions, with code released in a public repository.","tokens_in":1181,"tokens_out":1598,"duration_ms":293008,"significance":"If the claimed result is correct, it fills a genuine gap: until now the O(e^2) seagull and contact terms for generalized domain-wall fermions were missing, forcing RM123-style calculations with Mobius/zMobius actions to use a non-gauge-invariant local current. The derivation is explicit, self-contained, and parameter-free, and the numerical validation against finite-difference ground truth is a strong check. The public implementation in Grid is also a practical strength. However, the displayed operator-insertion formula Eq. (90) contains a sign error in the anti-quark contact term, which must be corrected before the paper can be accepted as stated.","major_comments":[{"comment":"Equation (90) is internally inconsistent with Eqs. (85) and (89). The last operator in Eq. (85) is -S_dw (Gamma2 dot A^2) C. Using the paper's own definition Qbar2^mu = -psi_bar Gamma2^mu C R5 P (Eq. (89)), the sandwich [P^-1(-S_dw Gamma2 C R5 P)]_11 equals +<q Qbar2>, so 1/2 partial^2_e S_q contains +1/2<q Qbar2>, not -1/2 as printed in Eq. (90). The same sign follows from the path-integral derivative: partial^2_e<q qbar> contains <q partial^2_e qbar> = <q Qbar2>, and the overall factor 1/2 gives +1/2. The Feynman-diagram expression in Eq. (95) already displays the +1/2 sign, so Eq. (90) is the outlier. This is not a convention ambiguity; it follows from the definitions stated in the paper. Since Eq. (90) is the main operator-insertion formula, an implementation based on it would insert the anti-quark contact term with the wrong sign. The numerical validation in Sec. VII tests Eq. (85)","section":"Sec. IV, Eq. (90)"}],"minor_comments":[{"comment":"The validation is performed on a single gauge configuration and a single source. This is acceptable for a consistency check, but should be stated clearly in the caption. In addition, reporting the fitted slopes or numerical values of epsilon1 and epsilon2 at several e^2 would make the claim that the errors scale linearly with e^2 more quantitative than a visual inspection of points.","section":"Sec. VII, Fig. 1"},{"comment":"The function F is used without an explicit definition. Please define F(D) = D^-1 partial D / partial mu (or state that it denotes the integrand in Eq. (111)), so that the trace identity tr F(D_ov) = tr F(D_dw(m)) - tr F(D_dw(1)) is unambiguous.","section":"Sec. VI, Eq. (114)"},{"comment":"The analytic continuation is written as eA_nu -> i mu delta_nu4. It would help to state explicitly that derivatives are taken with respect to the variable eA_nu, not A_nu, to avoid apparent missing factors of 1/e.","section":"Sec. VI, Eq. (115)"}],"recommendation":"major_revision","confidential_remarks":"The central derivation, the numerical validation, and the Ward identity all appear sound. The problem is the sign in Eq. (90): the anti-quark contact term Qbar2 is printed with the wrong sign, contradicting Eqs. (85), (89), and (95). This is a local, fixable error, but it is load-bearing because Eq. (90) is one of the main results that practitioners would implement. I expect that after the authors correct the sign and add a short consistency check, the paper could be acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWorth a referee's time, but Eq. (90) has a sign error that should be fixed before the paper is used as a reference.\n\nThe new content is real: the O(e^2) seagull current T^mu and the anti-quark contact terms Qbar1/Qbar2 for c_s != 0. The derivation from the U(3)-link Wilson operator is explicit and the signs and factors in Eq. (85) check out. The Ward identity (97) tying the seagull to the vector current is clean, and the numerical validation on an 8^4 lattice shows the expected e^2 scaling for Shamir, scaled Shamir, and zMöbius DWF. That validation tests Eq. (85), not Eq. (90).\n\nThe soft spot is the Qbar2 sign. From Eq. (85) and the definitions Qbar1 = -psi bar Gamma1 C R5P and Qbar2 = -psi bar Gamma2 C R5P, the O(e^2) correction to the physical propagator should be -<VV> + <V Qbar1> - 1/2<T> + 1/2<Qbar2>. Eq. (90) writes the last term as -1/2<Qbar2>. The sign flip is not a convention choice; it follows from the paper's own definitions. Any implementation based on Eq. (90) would insert the Qbar2 contact interaction with the wrong sign. I checked this by differentiating Eq. (85) and contracting the projector; the minus in Qbar2 flips the sign.\n\nTwo smaller issues. Ref. [34] has no URL, so the claimed public code cannot be located. And the paper leans on the standard qbar = [(psi bar D_-) R5P]_1 identification, which is fine — that is the established overlap/DWF field definition, not a weakness.\n\nFor the lattice QCD audience computing radiative corrections with Möbius/zMöbius DWF, this is a necessary ingredient. The derivation in Eq. (85) is the actual technical result and it looks right. Send it to peer review; a referee should ask for the sign correction in Eq. (90) and a locatable repository. Once those are fixed, it is a solid contribution.","headline":"Solid, useful derivation of O(e^2) DWF seagull vertices — but Eq. (90) carries a sign error in the Qbar2 contact term that must be fixed.","tokens_in":16105,"tokens_out":17486,"would_cite":true,"duration_ms":140508,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.20.-m"],"model":"deepseek-v4-flash","headline":"Generalized domain-wall fermions require new O(e^2) seagull and anti-quark contact terms for gauge-covariant QED corrections.","keywords":["domain-wall fermions","lattice QCD","QED corrections","seagull vertex","Ward-Takahashi identity","conserved vector current","background-field expansion","radiative corrections"],"falsifier":"Compute the left and right sides of the lattice Ward-Takahashi identity for the two-current correlator on a small volume with a fixed U(1) background and check whether the difference vanishes when the seagull term is dropped. Alternatively, replace the physical anti-quark definition by D_- evaluated at e=0 and compare the O(e^2) propagator derivative with a direct finite-difference derivative; any mismatch would show the contact terms are artifacts of the chosen convention.","tokens_in":15070,"feed_emoji":"⚛️","tokens_out":5714,"duration_ms":59924,"temperature":0.7,"pith_summary":"This paper establishes that when physical quark fields in generalized domain-wall fermions are defined through the gauge-link-dependent operator D_-, the expansion of the quark propagator in the electromagnetic charge e acquires new local insertions at second order: a seagull current and two anti-quark contact terms. These terms are not artifacts; they are required for the lattice U(1) Ward-Takahashi identity to hold and for O(alpha_em) radiative corrections obtained by the standard perturbative expansion around e=0 to be gauge covariant. The paper rederives the known first-order conserved current from a background-field method, derives the previously missing second-order seagull vertices, validates them numerically against direct derivatives of the propagator, and extends the pattern to all orders.","feed_headline":"Two contact terms complete QED on domain-wall lattices","feed_subtitle":"Without them, O(alpha_em) lattice calculations violate the U(1) Ward identity.","key_machinery":"The load-bearing object is the definition of the physical anti-quark field, qbar(x) = [(psibar D_-)(x) R5 P]_1, where D_- = diag(1 - c_s D_w) contains the QCD+QED Wilson gauge links. Because D_- carries the U(3) links, derivatives of the field generate contact terms not present in formulations with c_s = 0. The second-order seagull current T^mu(z) = sum_{x,y} psibar(x) Gamma_2^mu(x,z,y) Omega(m) psi(y), built from a point-split kernel Gamma_2^mu, and the anti-quark contact currents Qbar_1^mu and Qbar_2^mu are the new local insertions that carry the argument; together they enforce the Ward-Takahashi identity and make the O(e^2) propagator gauge covariant.","core_discovery":"The central claim is that the physical anti-quark field in generalized domain-wall fermions, qbar = (psibar D_-) R5 P with D_- = 1 - c_s D_w, depends on the QCD+QED gauge links. Expanding this definition in e produces operator insertions that are absent in the standard c_s=0 case. At O(e) the expansion yields the conserved vector current V^mu plus a contact term Qbar_1^mu. At O(e^2) it yields a double-current insertion, a single second-order insertion involving the seagull current T^mu(z) = sum_{x,y} psibar(x) Gamma_2^mu(x,z,y) Omega(m) psi(y), and anti-quark contact terms Qbar_1^mu and Qbar_2^mu. With these terms, the lattice Ward-Takahashi identity for the two-current correlator gives the","pith_inferences":["Because the extra terms come entirely from the gauge-link dependence of D_-, any future reformulation of domain-wall fermions that changes this field definition will need its own seagull derivation; the present calculation provides the template but not a uniqueness proof.","Existing generalized-domain-wall results that used local currents plus renormalization factors are not saved by renormalization at the level of gauge invariance; their systematics should be re-evaluated observable-by-observable before being used for sub-percent precision.","A direct test is to measure the Ward identity residual on a lattice with a non-trivial U(1) background; the residual should vanish only when both the seagull and anti-quark terms are included, giving a check independent of numerical derivative validation.","The all-order pattern suggests that seagull-like terms appear at every even order, so a resummed all-order expression for the propagator could be constructed and might simplify future automated perturbation theory."],"forward_implications":["Any O(alpha_em) calculation with generalized domain-wall fermions that drops the seagull or anti-quark contact terms will violate the lattice U(1) Ward-Takahashi identity and lose gauge covariance.","The conserved vector current and the first-order Ward identity are recovered from the same background-field expansion, unifying the first-order description.","The vacuum polarization tensor remains transverse, qhat_mu Pi_mu_nu = 0, only when the seagull term is included in the two-current correlator.","The pattern alternates: at every odd order the point-split kernel is Gamma_1^mu and at every even order Gamma_2^mu, so higher-order radiative corrections can be generated systematically.","The same contact structures apply to imaginary chemical potential through analytic continuation, so quark number susceptibilities in finite-temperature QCD need the same corrections."],"fun_headline_variants":["Seagull vertices restore QED Ward identity","Gauge covariant QED needs new seagull terms","Second-order QED corrections demand extra insertions","Domain-wall QED: missing terms violate Ward identity","New seagull operators complete lattice QED"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire result rests on the chosen definition qbar(x) = [(psibar D_-)(x) R5 P]_1, with D_- containing the QCD+QED links; if a different field identification were used, every second-order contact term would change.","fun_headline_variants_meta":{"raw":{"variants":["Seagull vertices restore QED Ward identity","Gauge covariant QED needs new seagull terms","Second-order QED corrections demand extra insertions","Domain-wall QED: missing terms violate Ward identity","New seagull operators complete lattice QED"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1350,"prompt_tokens":693,"completion_tokens":657,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":437,"completion_tokens_details":{"reasoning_tokens":582}},"tokens_in":437,"tokens_out":657,"duration_ms":7573,"temperature":1.0,"reasoning_tokens":582,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:43:38.336688+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the left and right sides of the lattice Ward-Takahashi identity for the two-current correlator on a small volume with a fixed U(1) background and check whether the difference vanishes when the seagull term is dropped. Alternatively, replace the physical anti-quark definition by D_- evaluated at e=0 and compare the O(e^2) propagator derivative with a direct finite-difference derivative; any mismatch would show the contact terms are artifacts of the chosen convention.","supporting_citations":[],"review_version":1}