{"id":"88bc5f21-cca5-4ffa-ad43-35da1a82596b","arxiv_id":"2501.11262","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A continuum-estimated lattice QCD calculation shows that pi0 and K0 screening masses respond non-monotonically to magnetic fields while eta_s(s-bar) decreases monotonically, alongside magnetic catalysis and inverse magnetic catalysis in the quark condensates.","lead":"This lattice QCD study computes how quark condensates and the screening masses of neutral pseudoscalar mesons change with temperature and magnetic field at physical quark masses. It provides continuum-estimated data that can serve as benchmarks for effective models of QCD in strong magnetic fields.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Disconnected contributions to the π0 screening correlator are neglected without quantification; since the magnetic field breaks u-d degeneracy, the B=0 cancellation no longer holds and the central non-monotonic claim may shift.","rationale":"The reader's weakest assumption is the neglect of disconnected quark-line contributions to neutral pseudoscalar screening correlators, especially for π0. My analysis agrees and strengthens the concern: the single most load-bearing point is the quantitative effect on the π0 screening mass. The paper's own Figure 1 verifies the WT identity only for the integrated correlator (susceptibility), not for the asymptotic exponential decay that defines the screening mass. The disconnected contribution is exactly zero at B=0 by u-d symmetry but becomes nonzero once magnetic fields split the u and d propagators; Figure 4 shows the u-d condensate difference grows with eB, so the neglected term is not obviously negligible in the parameter range studied. The cited justification [34,55] comes from zero-temperature or heavier-mass contexts, not from the physical-point near-crossover regime explored here. The continuum-extrapolation concern (three spacings, two ansätze) is real but affects the quoted error bars and benchmark status, not the qualitative non-monotonic shape as directly. The abstract's WT-identity phrasing is overstated but not load-bearing for the measured screening masses. The paper is otherwise solid: it uses HISQ at physical quark masses, three lattice spacings, public data, and multiple-state fitting with model selection. The conditional verdict is appropriate because a targeted check of the disconnected term could settle the issue without invalidating the work.","tokens_in":15629,"tokens_out":9188,"duration_ms":90607,"concrete_test":"On existing Nτ=8 configurations at T≈151 MeV and eB≈0.8 GeV², compute the disconnected π0 screening correlator with stochastic sources and extract the effective screening mass from the full connected+disconnected correlator; compare with the connected-only result. If the mass shifts by more than the quoted error, or the non-monotonic dip in eB is reduced, the central claim needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III states: 'The contribution of the disconnected quark lines to the screening mass was neglected in this analysis, as their impact is expected to be small [34,55]' and Section IV repeats this for the neutral pion. For the π0 interpolator (uū − dd̄)/√2, the disconnected part of the correlator is proportional to ([Tr γ5 S_u] − [Tr γ5 S_d])(x)([Tr γ5 S_u] − [Tr γ5 S_d])(y). At B=0, u and d propagators are identical and the term vanishes; in a magnetic field, the u/d charge asymmetry (quantified by ΔΣ_{u−d} in Figure 4) makes it nonzero and growing with eB. Figure 1 checks the WT identity for the integrated susceptibility, not for the asymptotic screening mass; a small integrated disconnected term does not imply a negligible long-distance component. If the disconnected hairpin correlator has a smaller screening mass than the connected π0 at high eB, where inverse catalysis raises the connected pion mass, the neglected term could dominate at large z and flatten or remove the rise in the π0 screening mass, directly weakening the central claim. The cited references [34,55] do not quantify this near Tpc at the physical point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents lattice QCD results for chiral condensates and neutral pseudoscalar screening masses in a background magnetic field, using (2+1)-flavor HISQ fermions at physical quark masses. Simulations are performed on N_tau = 8, 12, and 16 lattices with aspect ratio 4, at five temperatures between about 145 and 166 MeV and eight magnetic field values up to eB ~ 0.8 GeV^2. The paper reports continuum estimates obtained by interpolating in the T–eB plane and extrapolating in 1/N_tau^2 using two ansätze. The main physical results are that the light and strange-light chiral condensates show magnetic catalysis at small eB followed by inverse magnetic catalysis at larger eB, that the pi0 and K0 screening masses are non-monotonic in eB and track the corresponding condensates, and that the eta_{s sbar} screening mass decreases monotonically with eB. The analysis includes a Ward-Takahashi identity check relating the condensates to the integrated pseudoscalar correlators.","tokens_in":15878,"tokens_out":7393,"duration_ms":78311,"significance":"If the results are correct, this is the first continuum-estimate study of these observables at physical quark masses in a thermomagnetic medium near the chiral crossover. The paper uses a standard and defensible lattice setup, with three lattice spacings, multiple correlator fit states selected by AICc, a B-spline interpolation procedure, two continuum-extrapolation ansätze, and a public dataset deposit. These are genuine strengths that make the paper a useful benchmark for effective models. The central caveat, acknowledged in the text, is the neglect of disconnected quark-line contributions to the neutral pion and kaon screening correlators; because the main quantitative claims concern precisely those channels, this issue is load-bearing and needs to be addressed before the results can be fully accepted.","major_comments":[{"comment":"The neglect of disconnected quark-line contributions to the neutral pseudoscalar screening correlators is not quantified. At zero magnetic field the disconnected contribution to the pi0 correlator cancels by u–d degeneracy, but at eB != 0 that cancellation is broken by the charge asymmetry, which the paper itself shows grows with eB in Figure 4 (Delta Sigma_{u-d}). The references cited for the expectation that the contribution is small, [34] and [55], are zero-temperature or different-setup studies and do not provide a bound at T near Tpc with physical quark masses. The Ward-Takahashi check in Figure 1 is performed on the integrated susceptibility, which is dominated by short-distance correlator contributions; it does not constrain the asymptotic, long-distance screening mass that defines M_{pi0}. If the disconnected hairpin correlator has a smaller screening mass than the connected pi0 channel at large eB, it could dominate the large-z behavior and change the extracted non-monotonic dependence. The authors should compute the disconnected contribution on at least a subset of ensembles, or provide an explicit model-based upper bound and demonstrate that the screening-mass results are stable under including it.","section":"Section III (lattice observables)"},{"comment":"The pion correlator is defined as G_{pi0} = (G_{uu} + G_{dd})/2 with equal weights for the up and down quark contributions. This is presented as an assumption in a footnote, but it is not a symmetry statement in the presence of a magnetic field, where the u/d charge splitting breaks SU(2)_V. The physical neutral pion is a field-dependent combination of the light flavors, and the choice of equal weights defines the interpolating operator rather than the mass eigenstate. The paper does not estimate the systematic uncertainty that this operator choice introduces. The assumption should be moved into the main text and justified, or the extraction should be interpreted as the screening mass of this particular operator rather than the physical pi0 mass.","section":"Figure 5 caption"}],"minor_comments":[{"comment":"The abstract states the temperature range as 145 MeV to 166 MeV, while Section III says 145 MeV to 165 MeV; the tables list temperatures such as 165.98 and 166.03 MeV. Please harmonize the quoted range.","section":"Abstract / Section III"},{"comment":"The caption says the ratios are 'normalized with respect to their corresponding quark masses and susceptibility' but does not define the plotted quantity explicitly. Please give the explicit ratio, for example (mu+md) chi_{pi0} / ( <psi-bar psi>_u + <psi-bar psi>_d ).","section":"Figure 1"},{"comment":"The sums over i and j in the multi-state fit ansatz are not defined. Please state explicitly that they run over excited states and that A and M are fit parameters for the non-oscillating and oscillating contributions.","section":"Section II"},{"comment":"For the N_tau = 12 ensemble at beta = 6.712, the screening-mass configuration counts for Nb = 1 through 6 are all listed as 3157. This repetition looks like a possible placeholder or transcription error; please verify the entries.","section":"Appendix C"},{"comment":"The B-spline interpolation procedure requires a smoothing factor, but its value is not reported. A sentence giving the chosen smoothing factor and its sensitivity would improve reproducibility.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid and useful contribution, and the main results are likely correct, but the unquantified disconnected contribution is a genuine load-bearing issue for the central pi0/K0 claims. I would encourage the authors to address it with a direct calculation or a conservative estimate on a subset of ensembles. The paper is otherwise well within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the physical-mass, three-lattice-spacing follow-up to Ding et al. [26], and it mostly delivers what it promises: continuum estimates of DeltaSigma_ud/ds/s and neutral pseudoscalar screening masses up to eB ~ 0.8 GeV^2 near Tpc. The setup is standard and well documented, the data are public via the Zenodo DOI in Ref. [58], and the new content is real. Prior continuum work on condensates in magnetic fields did not use physical masses, and there was no continuum estimate of pi0/K0/eta_s screening masses at the physical point. That alone is worth having.\n\nThe paper does several things well. Three lattice spacings with fixed aspect ratio, multi-state correlator fits with AICc selection, B-spline interpolation in the T-eB plane, and a two-ansatz continuum average are all reasonable choices. The WT identity check on lattice data is a useful consistency test, and the separate u/d condensates plus DeltaSigma_{u-d} are a nice addition. The qualitative picture of competing magnetic catalysis and inverse magnetic catalysis, with nonmonotonic pi0 and K0 screening masses and a monotonically decreasing eta_s, is plausible and consistent with earlier heavier-pion work and effective models.\n\nThe real soft spot is the neglected disconnected quark-line contribution to the pi0 screening correlator. The authors say it is expected to be small, citing Refs. [34,55], but at nonzero eB the u and d propagators are not degenerate, and Figure 4 itself shows the u-d condensate difference growing with eB. The WT identity check in Figure 1 is for the integrated susceptibility, not for the asymptotic screening mass; a small integrated disconnected term does not guarantee the long-distance part is small. If the disconnected hairpin channel has a lower screening mass than the connected pi0, it could flatten or pull down the extracted pi0 mass at large eB and soften the central non-monotonic claim. This should be quantified, at least at one lattice spacing and the largest eB, before the pi0 result is treated as a benchmark.\n\nTwo minor issues: the abstract's wording that condensates are \"intrinsically connected\" to screening masses via WT identities overstates the case, since the identities hold for integrated susceptibilities and the paper itself acknowledges screening masses probe long distances. And the continuum estimate rests on only three spacings and two ansatze; acceptable, but the systematic error could be larger than the bands suggest.\n\nThis is a solid lattice calculation with a legitimate, addressable caveat. I would send it to a serious referee and expect it to be publishable after the disconnected contribution is quantified or its potential impact discussed. I would cite it for the continuum benchmark values.","headline":"A solid, public-data lattice QCD benchmark at physical quark masses for chiral condensates and neutral pseudoscalar screening masses in magnetic fields, with one real caveat about neglected disconnected diagrams in the pi0 channel.","tokens_in":16433,"tokens_out":2108,"would_cite":true,"duration_ms":22878,"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":"Lattice QCD with physical quark masses yields continuum estimates of how magnetic fields reshape chiral condensates and neutral meson screening masses near the QCD crossover.","keywords":["lattice QCD","chiral condensate","screening mass","neutral pseudoscalar mesons","magnetic catalysis","inverse magnetic catalysis","continuum limit","HISQ fermions"],"falsifier":"Compute the disconnected contribution to the π⁰ screening correlator on the finest lattice (Nτ = 16) at T ≈ 157 MeV and eB ≈ 0.2–0.4 GeV², where the non-monotonic dip is most visible; if including it moves the extracted screening mass by more than the quoted uncertainty, the claimed non-monotonic eB dependence of m_{π⁰} is not robust.","tokens_in":15405,"feed_emoji":"🧲","tokens_out":6103,"duration_ms":53054,"temperature":0.7,"pith_summary":"This paper uses (2+1)-flavor lattice QCD with physical quark masses to determine, for the first time in the continuum limit, how strong magnetic fields change the chiral condensates and the screening masses of the neutral pseudoscalar mesons π⁰, K⁰, and η⁰_{s\\bar{s}} near the chiral crossover. The central result is that the screening masses of π⁰ and K⁰ depend non-monotonically on the magnetic field strength eB, first decreasing and then rising, closely tracking the non-monotonic behavior of the light and strange-light chiral condensates; the screening mass of the fictitious strange eta, η⁰_{s\\bar{s}}, instead decreases monotonically. These continuum estimates, covering temperatures 145–166 MeV and fields up to eB ≈ 0.8 GeV², quantify the competition between magnetic catalysis and inverse magnetic catalysis in spatial correlation lengths, and provide benchmarks that low-energy QCD models and effective theories can be tested against.","feed_headline":"π⁰ and K⁰ screening masses dip, then rise, with magnetic field","feed_subtitle":"Continuum estimates at physical quark masses show inverse magnetic catalysis in spatial correlation lengths — benchmarks for QCD models.","key_machinery":"The load-bearing object is the Ward–Takahashi identity (m_u+m_d)χ_{π⁰}=⟨ψ̄ψ⟩_u+⟨ψ̄ψ⟩_d (and its K⁰ and η⁰_{s\\bar{s}} analogues), which ties each chiral condensate combination to the space-time integral of the corresponding pseudoscalar correlation function; the screening mass is then the inverse correlation length extracted from the same spatial correlator. On the lattice, the correlators are computed with highly improved staggered (HISQ) fermions in a quantized magnetic flux background, and the screening masses are obtained from multi-state $\\cosh$ fits with an oscillating parity-partner term; continuum estimates come from averaging linear and quadratic 1/Nτ² extrapolations of the Nτ = 8, 12, 16 data.","core_discovery":"The paper establishes that in QCD with physical up, down, and strange quark masses, the magnetic-field dependence of the neutral pseudoscalar screening masses is governed by the same competition between magnetic catalysis and inverse magnetic catalysis that shapes the chiral condensates. Using the Ward–Takahashi identities that tie each condensate combination to a pseudoscalar susceptibility, the authors extract screening masses from the exponential decay of spatial correlation functions on Nτ = 8, 12, and 16 lattices and extrapolate to the continuum with linear and quadratic ansätze in 1/Nτ². The resulting continuum estimates show that ΔΣ_ud and ΔΣ_ds rise to a peak and then fall as eB grows at low temperature, with the peak shifting to smaller eB as T increases, while ΔΣ_s keeps rising in the covered window; at higher temperatures all three develop more intricate rise–fall–rise patterns. Correspondingly, the screening masses of π⁰ and K⁰ first decrease and then increase with eB, mirroring their condensates, whereas m_{η⁰_{s\\bar{s}}} decreases monotonically. All screening masses increase with T, with steeper slopes at larger eB, and the crossing of constant-eB curves is interpreted as the magnetic-field-induced reduction of the pseudocritical temperature.","pith_inferences":["If the disconnected-diagram contribution to the π⁰ correlator is not negligible, the reported non-monotonic dip in m_{π⁰}(eB) could be partially an artifact; a dedicated computation of the disconnected part on the Nτ = 16 ensemble would settle this without a full new simulation campaign.","The close tracking between screening masses and condensates suggests that at these temperatures the pseudoscalar screening masses might be expressible through a generalized Gell-Mann–Oakes–Renner relation with eB-dependent decay constants, which could be tested by measuring the amplitudes A_H of the correlators.","The crossing of constant-eB curves in m_{π⁰}(T) could be used as an alternative, correlation-based definition of T_pc(eB); comparing it with the inflection-point definition would show whether the field-induced T_pc reduction is observable-independent."],"forward_implications":["The non-monotonic screening masses mean that inverse magnetic catalysis is not only a short-distance/integrated effect: it extends to the long-distance spatial correlation lengths that govern how mesonic excitations screen color fields in the medium.","The peak of ΔΣ moving to smaller eB with temperature implies a T-dependent boundary between magnetic catalysis and inverse magnetic catalysis, consistent with a falling T_pc(eB).","Continuum estimates of ΔΣ_ud, ΔΣ_ds, and ΔΣ_s provide direct targets for NJL-type models, the linear sigma model, and holographic AdS/QCD constructions that currently disagree on whether inverse magnetic catalysis appears.","The monotonic decrease of the η⁰_{s\\bar{s}} screening mass with eB shows that strange-quark pseudoscalars remain in the magnetic-catalysis regime across the whole temperature window, so model comparisons should treat light and strange channels separately.","The growing up–down condensate asymmetry ΔΣ_{u−d} with eB, which shrinks as T rises, offers an observable signature of charge-dependent chiral symmetry breaking that could be probed by future simulations on larger volumes."],"supporting_citations":[{"why":"The previous heavier-pion study that established the connection between condensates and screening masses and reported similar non-monotonic behavior; this work extends it to physical quark masses and the continuum.","marker":"[26]"},{"why":"Provides the lattice determination of the QCD quark condensate in external magnetic fields whose continuum estimates are compared with here.","marker":"[21]"},{"why":"Supplies the benchmark result for the eB-dependence of the pseudocritical temperature and phase diagram that this paper's crossing behavior is compared with.","marker":"[22]"},{"why":"The zero-temperature (2+1)-flavor chiral-property study that supplies the Ward-Takahashi identities, magnetic-field implementation, and the strange-quark mass tuning procedure used here.","marker":"[34]"},{"why":"The (2+1)-flavor meson screening-mass study that provides the correlator fit form (oscillating cosh ansatz) and scale-setting parametrization.","marker":"[41]"},{"why":"The effective-model calculation that explained how sea-quark inverse magnetic catalysis makes the neutral pion screening mass rise with eB, which the lattice results now confirm.","marker":"[32]"}],"fun_headline_variants":["π⁰ and K⁰ screening masses fall then rise with magnetic field","Magnetic field makes π⁰ and K⁰ screening masses dip and rebound","Neutral pion and kaon screening lengths dip then climb with B","Lattice QCD: π⁰ and K⁰ screening masses non-monotonic in eB","Non-monotonic screening: π⁰ and K⁰ in a magnetic field"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The extraction of the π⁰ and K⁰ screening masses assumes that disconnected quark-line contributions to the neutral pseudoscalar spatial correlators are negligibly small, an assumption the paper states but does not quantitatively verify at physical quark masses.","fun_headline_variants_meta":{"raw":{"variants":["π⁰ and K⁰ screening masses fall then rise with magnetic field","Magnetic field makes π⁰ and K⁰ screening masses dip and rebound","Neutral pion and kaon screening lengths dip then climb with B","Lattice QCD: π⁰ and K⁰ screening masses non-monotonic in eB","Non-monotonic screening: π⁰ and K⁰ in a magnetic field"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001361,"raw_usage":{"total_tokens":5632,"prompt_tokens":1167,"completion_tokens":4465,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":783,"completion_tokens_details":{"reasoning_tokens":4357}},"tokens_in":783,"tokens_out":4465,"duration_ms":33079,"temperature":1.0,"reasoning_tokens":4357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:28:06.306212+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the disconnected contribution to the π⁰ screening correlator on the finest lattice (Nτ = 16) at T ≈ 157 MeV and eB ≈ 0.2–0.4 GeV², where the non-monotonic dip is most visible; if including it moves the extracted screening mass by more than the quoted uncertainty, the claimed non-monotonic eB dependence of m_{π⁰} is not robust.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The previous heavier-pion study that established the connection between condensates and screening masses and reported similar non-monotonic behavior; this work extends it to physical quark masses and the continuum."},{"cited_title":"Detar and J.B","cited_arxiv_id":null,"evidence_quote":"The (2+1)-flavor meson screening-mass study that provides the correlator fit form (oscillating cosh ansatz) and scale-setting parametrization."}],"review_version":1}