{"id":"a4a7415a-c802-4d4b-8a55-272b81d1bb8a","arxiv_id":"2608.11144","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In chiral perturbation theory, the QCD vacuum responds nonlocally to an inhomogeneous magnetic field, yielding parameter-free NLO predictions for the free energy, chiral condensate, and an induced vacuum current.","lead":"This paper calculates how the vacuum of the strong nuclear force reacts to a magnetic field that varies across space, using the low-energy effective theory of QCD. The reaction is nonlocal and includes an induced current that vanishes when the field is uniform.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gradient-operator concern is answered by standard ChPT power counting; no load-bearing error identified, though the abstract's 'without undetermined parameters' overstates the body.","rationale":"The reader correctly identifies the absence of gradient operators at NLO as the most load-bearing assumption for the parameter-free claim. I examined this assumption against the standard chiral power counting and found it to be sound: f_+ is O(p^2), a derivative is O(p), and no independent parity-even O(p^4) invariant with a gradient of f_+ survives reduction to the known NLO basis. The paper's assertion that gradient effects enter only at NNLO is therefore consistent with ChPT, and the single counter-term h is sufficient for renormalization. I also checked the internal consistency of the main results: the integrated free energy, condensate, and current are properly renormalized; the local quantities are consistent with the integrated ones up to the admittedly unproven identity in Appendix B, which is verified numerically and analytically to O(b^6). The reader's three conditions for the CONDITIONAL verdict are reasonable but none is load-bearing: the Appendix B identity would affect only the local free-energy representation, the abstract overclaim is corrected by the body, and the absence of code/data is a reproducibility issue rather than a scientific one. My pass therefore does not change the verdict: the calculation is careful, internally consistent, and delivers parameter-free NLO predictions for the width-dependent deviations from uniform-field behavior, as long as the reader interprets 'without undetermined parameters' as referring to those differences and not to the overall B^2 terms.","tokens_in":27637,"tokens_out":18425,"duration_ms":239337,"concrete_test":"Independently enumerate the independent O(p^4) chiral invariants containing one f_+ and one covariant derivative (e.g., <f_+^{mu nu} nabla_mu u_nu>), and check whether integration by parts and the leading-order equations of motion reduce all such terms to the standard L9/L10/H1/H2 basis used in Eq. (2.7). If an independent invariant remains, the paper's counter-term set would be incomplete; if none remain, the gradient-operator assumption is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central parameter-free NLO program relies on the claim that no gradient operators of the magnetic field enter the chiral Lagrangian at NLO. In standard two-flavor ChPT, the external field strength f_+ counts as O(p^2), so f_+^2 terms are O(p^4) and appear in the NLO Lagrangian. A derivative of f_+ is O(p^3); a parity-even O(p^4) invariant would require pairing it with an O(p) building block such as u_mu. Such combinations are not independent: integration by parts and the leading-order equations of motion reduce them to the standard invariants of the Gasser-Leutwyler basis (L9, L10, H1, H2) appearing in Eq. (2.7). Gradient effects therefore first enter at NNLO, O(p^6), exactly as the abstract asserts. The h counter-term found from the uniform-field limit is sufficient to cancel the logarithmic divergence in the inhomogeneous case (Eq. (3.19) vs. Eq. (2.34)). No internally inconsistent step was found in the derivation of the integrated or local observables. The only real weakness is the abstract's phrasing: the total free energy in Eq. (3.20) and the full vacuum current in Eq. (4.31) retain the unknown LEC h in the B^2/dBdx contact terms; the parameter-free predictions are the width-dependent differences (renormalized susceptibility, condensate ratios, and the pion-loop contribution to the current). This is a wording issue, not a mathematical flaw.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the next-to-leading-order (NLO) chiral perturbation theory response of the QCD vacuum to a static, localized magnetic field with profile B(x)=B sech^2(x/lambda). Using the exact resolvent of the associated one-dimensional Schrödinger problem, it derives ultraviolet-finite expressions for the integrated and local free energy, the chiral condensate, and an induced vacuum current. The uniform-field limit is checked against known results, and the local free-energy density is checked against the integrated result in Appendix B. The paper demonstrates that the locally constant approximation fails when the field width is of order the pion Compton wavelength and that a gradient-induced vacuum current appears. The advertised parameter-free statements apply to width-dependent, inhomogeneity-induced parts of the observables after subtraction of LEC-dependent contact terms.","tokens_in":27895,"tokens_out":10998,"duration_ms":119367,"significance":"If the results hold, this is a useful and nontrivial extension of ChPT to spatially varying external fields. The solvable profile enables closed-form NLO expressions, and the renormalization is treated carefully, with cross-checks including the uniform-field limit, numerical integration, and series expansions through O(b^6). The paper also identifies a factor-of-two discrepancy with Ref. [13] and provides supporting consistency checks. The predicted failure of the locally constant approximation and the inhomogeneity-induced vacuum current are concrete, potentially falsifiable signatures that can be compared with lattice QCD in non-uniform magnetic backgrounds. The main caveat is that the abstract's claim of 'without undetermined parameters' is stronger than what the body actually establishes; the parameter-free statements hold for subtracted, width-dependent quantities rather than for all NLO observables as literally written.","major_comments":[{"comment":"The abstract's claim that equilibrium vacuum observables can be determined at NLO 'without undetermined parameters' overstates the body of the paper. The full integrated free energy in Eq. (3.20) and the full vacuum current in Eq. (4.31) retain the LEC h in the Maxwell-type contact terms, through chi_B and the dB/dx coefficient respectively. The parameter-free statements apply to the renormalized integrated susceptibility chi_B^r(lambda) in Eq. (3.22), to the condensate ratios in Eqs. (3.27)-(3.29) and Eq. (4.23), and to the pion-loop current after the linear response is subtracted in Eqs. (4.38)-(4.39). Please revise the abstract and conclusions to state this qualification explicitly.","section":"Abstract; §3.2 Eq. (3.20); §4.3 Eq. (4.31)"},{"comment":"The trace of the Green's function is the input to all integrated observables. Appendix A states that the result is a factor of two smaller than the corresponding result in Ref. [13] and attributes the discrepancy to a typo in that reference, supported only by a statement that the endpoint behavior matches and that a d=3 check agrees. Because this factor is load-bearing for Eqs. (3.13), (3.17), and every integrated quantity that follows, the derivation of Eq. (A.7) should be given in enough detail to verify the factor of two, or the d=3 comparison should be shown explicitly.","section":"Appendix A, Eq. (A.7), and Eq. (3.13)"},{"comment":"The paper explicitly states that it has not been able to prove analytically that the local free-energy density integrates to the closed-form integrated result; Eq. (B.4) is verified only numerically and via series expansions through O(b^6). Since this identity is asserted to be required for consistency, its status should be clarified: either supply an analytic proof or clearly label the local-integrated agreement as a numerically supported check rather than a derived relation. If the latter, the statements in §4.1 presenting Eq. (4.15) as the local free-energy density should not rely on this unproven identity.","section":"Appendix B, Eqs. (B.2)-(B.4)"}],"minor_comments":[{"comment":"The text says 'Using the renormalized magnetization in a uniform magnetic field Eq. (3.30)', but Eq. (3.30) is the integrated magnetization in the inhomogeneous field; the uniform-field magnetization is defined via Eqs. (2.45) and (2.32). Please correct the cross-reference.","section":"§4.3, after Eq. (4.40)"},{"comment":"The caption states that the susceptibility is plotted as a function of the width lambda, while the horizontal axis appears to be mu = lambda m_pi. Please make the axis definition explicit.","section":"Figure 1 caption"},{"comment":"The local free-energy density is representation-dependent up to a total derivative, as acknowledged in footnote 1, yet the discussion of Fig. 4 treats the spatial profile of the quadratic response as a meaningful local probe. Please add a sentence emphasizing that this profile is a property of the chosen resolvent representation, not a uniquely defined local observable.","section":"§4.1, footnote 1 and Eq. (4.19)"},{"comment":"The conclusion states that the quantities are determined 'non-perturbatively in both the width lambda ~ m_pi^{-1} and strength of the magnetic field eB ~ m_pi^2'. The strength statement should be qualified, since the NLO chiral expansion still requires eB/(4 pi F_pi)^2 << 1, as the paper itself notes in §3.3.","section":"§5, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically solid and the advertised physics is interesting. The recommended revision is driven by the need to make the factor-of-two correction in Appendix A and the Appendix B identity either fully derived or explicitly qualified, and to fix the overstatement in the abstract. Once these points are addressed, the manuscript would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First ChPT treatment of a spatially varying magnetic field, and it is a careful job. The genuinely new content: closed-form NLO expressions for the integrated and local free energy, the chiral condensate, and an induced vacuum current for the solvable B(x)=B sech^2(x/lambda) profile. The current is a clean nonlocality observable: it vanishes in a uniform field, and the locally constant approximation fails once the width is around the pion Compton wavelength. For the lattice community this gives an analytic benchmark for non-uniform backgrounds; for ChPT it demonstrates that gradient effects first enter at NNLO.\n\nThe execution is solid. Dimensional regularization is handled transparently: the logarithmic divergence is cancelled by the known h counterterm, the power-law divergences are isolated by add-subtract plus analytic continuation, and the uniform-field limit in Sec. 2.2 reproduces known results. The hardest internal check, equality of the one- and two-dimensional forms of the integrated free energy (Appendix B), is verified numerically and by series through O(b^6); the authors openly state they could not find an analytic proof. That is an honest limitation, minor rather than fatal—it is a consistency check, not a load-bearing claim. They also flag the ambiguity in defining a local free-energy density (footnote 1) and note a factor-of-two discrepancy with Ref. [13] in Appendix A with a plausible explanation. A referee should spot-check that last point, but nothing I read suggests a real problem.\n\nSoft spots, in proportion. The abstract's 'without undetermined parameters' is stronger than the body: the total free energy and the full vacuum current retain the LEC h in the B^2 and dB/dx contact terms. What is actually parameter-free is the set of width-dependent differences—renormalized susceptibility, condensate ratios, pion-loop current. This is a wording problem in the headline claim, not a mathematical flaw. Second, if a skeptic worries that gradient operators of the magnetic field enter at NLO, standard power counting answers the worry: f_+ is O(p^2), any derivative-of-f_+ combination needs an O(p) partner, and IBP plus the LO equations of motion reduce such terms to the standard Gasser–Leutwyler basis. The gradient effects really are NNLO.\n\nI found no load-bearing error, no circular fitting, and the use of the known LEC h is legitimate. The paper deserves a serious referee; with the abstract toned down and Appendix B's numeric check maybe strengthened, it should appear in JHEP largely as is. Send it to review.","headline":"A careful, internally consistent first ChPT treatment of inhomogeneous magnetic fields, with a new induced-current observable — just tone down the abstract's parameter-free claim.","tokens_in":28448,"tokens_out":5241,"would_cite":true,"duration_ms":50222,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that the QCD vacuum's response to a spatially varying magnetic field can be computed at next-to-leading order in chiral perturbation theory with no undetermined parameters, and that this response is genuinely nonlocal.","keywords":["chiral perturbation theory","inhomogeneous magnetic field","QCD vacuum","chiral condensate","vacuum current","magnetic susceptibility","dimensional regularization","nonlocal response"],"falsifier":"One concrete test: compute the induced vacuum current for the same sech-squared field profile in lattice QCD at widths around one inverse pion mass. If matching the chiral prediction requires a new low-energy constant beyond the one fixed by uniform fields, the claim that gradients enter only one order higher is wrong. A cheaper calculation is to evaluate the one-loop charged-pion determinant with an independent regulator; a surviving divergence proportional to the square of the field gradient would also signal a missing term.","tokens_in":27383,"feed_emoji":"🧲","tokens_out":10955,"duration_ms":98989,"temperature":0.7,"pith_summary":"This paper asks how the QCD vacuum responds when the external magnetic field is not uniform but varies in space, and it argues that this response can be computed at next-to-leading order in chiral perturbation theory with no free parameters. The concrete setting is a magnetic field localized along one spatial direction with the profile $B(x)=B\\,\\mathrm{sech}^2(x/\\lambda)$, a shape chosen because the underlying one-dimensional quantum problem is exactly solvable. The key structural claim is that gradient operators of the field enter the chiral Lagrangian only at next-to-next-to-leading order, so the uniform-field counterterm is enough to renormalize all inhomogeneous observables at the order considered. Using this, the paper obtains ultraviolet-finite expressions for the integrated and local free energy, the chiral condensate, and an induced vacuum current, and shows that pion fluctuations respond to the field profile over the pion Compton wavelength rather than to the local field value. A reader should care because this gives a parameter-free, first-principles effective-field-theory account of a nonlocal vacuum response that uniform-field studies cannot see.","feed_headline":"Magnetic-field shape exposes the QCD vacuum's nonlocal response","feed_subtitle":"Chiral perturbation theory predicts the QCD vacuum's response to a non-uniform magnetic field without free parameters.","key_machinery":"The carrying object is the resolvent Green's function $G(x,x|p_2,E)$ of the effective one-dimensional Schr\\\"odinger operator $-d^2/dx^2+V(p_2,x)+E$, where $V(p_2,x)=(p_2-eA(x))^2+m^2$ and $A(x)=\\lambda B\\tanh(x/\\lambda)$. Because this potential has a sech-squared (P\\\"oschl\\u2013Teller-type) shape, its eigenfunctions are hypergeometric functions and the coincident Green's function can be written in closed form, so the dimensionally regulated trace-log of the charged-pion operator can be evaluated explicitly. Three steps carry the argument: the uniform-field limit fixes the counterterm $h$; the same $h$ cancels the ultraviolet divergences of the inhomogeneous free energy and current, because gradient operators are postponed to next-to-next-to-leading order; and local observables are evaluated by changing variables to $(z_0,\\Delta)$, which isolates the divergent regions and leaves finite integrals for numerical evaluation.","core_discovery":"The paper establishes that, for the solvable profile $B(x)=B\\,\\mathrm{sech}^2(x/\\lambda)$, equilibrium vacuum observables at next-to-leading order are determined without undetermined parameters: the magnetic-field gradient does not generate new chiral Lagrangian operators until next-to-next-to-leading order, so the single low-energy constant $h$ fixed in the uniform-field problem renormalizes the inhomogeneous problem as well. Concretely, the charged-pion contribution to the effective action is evaluated through the resolvent Green's function of an effective one-dimensional Schr\\\"odinger equation whose potential is built from the gauge potential $A(x)=\\lambda B\\tanh(x/\\lambda)$; the resulting renormalized free-energy density, chiral condensate, and induced vacuum current are ultraviolet finite and are compared with locally constant approximations. The central physical finding is that for widths $\\lambda$ of order the inverse pion mass the locally constant approximation fails: the condensate is suppressed where the field is strongest and develops tails, and a nonzero vacuum current along the $y$-direction appears, driven by the magnetic-field gradient, which vanishes identically in a uniform field.","pith_inferences":["Inference: a measurement or lattice computation of the induced vacuum current at $\\lambda\\sim m_\\pi^{-1}$ would be a direct test of the gradient-counting argument, because a new next-to-leading-order counterterm would appear as a current component not captured by the uniform-field $h$.","Inference: the same exactly solvable resolvent method could be pushed to time-dependent backgrounds, turning the induced current into a real-time response and connecting the framework to dynamical electromagnetic fields in heavy-ion collisions.","Inference: the nonlocal magnetic susceptibility kernel $\\chi_B(x,x')$ implicit in the local quadratic free-energy response could be extracted from these expressions, yielding a first-principles estimate of the spatial correlation length of the magnetized QCD vacuum.","Inference: because the ratio of local to locally constant condensate is independent of the leading unknown constant, that ratio is a particularly clean observable for future lattice-QCD tests of the chiral prediction."],"forward_implications":["For widths $\\lambda$ of order the pion Compton wavelength, the locally constant approximation fails, so chiral perturbation theory predicts a nonlocal contribution to the condensate and free energy in non-uniform fields.","A vacuum current along the direction perpendicular to both the field and its gradient is induced whenever the magnetic field is inhomogeneous; it vanishes for uniform fields, giving a clean signature of inhomogeneity.","All next-to-leading-order results share the single counterterm fixed by the uniform-field problem, so the predictions are parameter-free and can be falsified by lattice or model calculations.","The renormalized integrated magnetic susceptibility is scale and scheme independent and returns to the uniform-field value as the width grows, quantifying how the response becomes local.","The linear part of the induced current is governed by the QCD vacuum polarization tensor in an inhomogeneous background, which the paper identifies as the next target for a full nonlocal characterization."],"supporting_citations":[{"why":"supplies the exactly solvable semilocalized magnetic-field Green's function and the closed-form integration technique used throughout.","marker":"[13]"},{"why":"provides the next-to-leading-order chiral Lagrangian and the counterterm structure from which $h$ is built.","marker":"[11]"},{"why":"establishes the proper-time effective-action method for vacuum polarization in background fields.","marker":"[1]"},{"why":"gives the uniform-field chiral condensate result that must be reproduced and serves as the locally constant benchmark.","marker":"[23]"},{"why":"supplies the isoscalar contact operator needed for the magnetic-field Lagrangian in Eq. (2.7).","marker":"[21]"},{"why":"provides the vacuum polarization tensor in inhomogeneous magnetic fields used to interpret the linear current response.","marker":"[29]"},{"why":"gives the lattice-QCD study of thermal QCD in a non-uniform magnetic background that this effective-theory calculation complements.","marker":"[19]"}],"fun_headline_variants":["Bumpy B-field: QCD vacuum response without new parameters","Magnetic-field gradient drives vacuum current in QCD","Inhomogeneous fields expose QCD vacuum's nonlocal side","Parameter-free QCD vacuum response to shaped magnetic fields","Non-uniform B-field predicts QCD vacuum's hidden current"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation rests on the assumption that a magnetic field's spatial variation does not introduce any new terms in the low-energy theory until one order higher than the order computed here; if a term involving the field's gradient were needed at the computed order, the one number fixed by uniform fields would no longer be enough.","fun_headline_variants_meta":{"raw":{"variants":["Bumpy B-field: QCD vacuum response without new parameters","Magnetic-field gradient drives vacuum current in QCD","Inhomogeneous fields expose QCD vacuum's nonlocal side","Parameter-free QCD vacuum response to shaped magnetic fields","Non-uniform B-field predicts QCD vacuum's hidden current"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1339,"prompt_tokens":939,"completion_tokens":400,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":333}},"tokens_in":555,"tokens_out":400,"duration_ms":5111,"temperature":1.0,"reasoning_tokens":333,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:17:57.547542+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete test: compute the induced vacuum current for the same sech-squared field profile in lattice QCD at widths around one inverse pion mass. If matching the chiral prediction requires a new low-energy constant beyond the one fixed by uniform fields, the claim that gradients enter only one order higher is wrong. A cheaper calculation is to evaluate the one-loop charged-pion determinant with an independent regulator; a surviving divergence proportional to the square of the field gradient would also signal a missing term.","supporting_citations":[{"cited_title":"Gasser and H","cited_arxiv_id":null,"evidence_quote":"provides the next-to-leading-order chiral Lagrangian and the counterterm structure from which $h$ is built."},{"cited_title":"Anomalies and WZW-term of two-flavour QCD","cited_arxiv_id":"hep-ph/0011377","evidence_quote":"supplies the isoscalar contact operator needed for the magnetic-field Lagrangian in Eq. (2.7)."},{"cited_title":"Vacuum polarization tensor in inhomogeneous magnetic fields","cited_arxiv_id":"1107.0286","evidence_quote":"provides the vacuum polarization tensor in inhomogeneous magnetic fields used to interpret the linear current response."}],"review_version":1}