{"id":"2610bfb3-9568-4337-92f1-cfeca7fbe6e7","arxiv_id":"2502.06216","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Continuum-estimated lattice QCD screening masses of neutral pseudoscalar mesons show a magnetic-field-induced minimum for pi and K, but not for the strange eta, near the QCD transition temperature.","lead":"Lattice QCD simulations with physical quark masses were used to measure how the screening masses of neutral pseudoscalar mesons change with temperature and magnetic field strength. The results show a non-monotonic behavior for pions and kaons, with a minimum at intermediate magnetic fields, while the strange eta meson decreases steadily.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The continuum estimate is under-validated: Eqs. (9)-(10) have zero residual degrees of freedom, so the claimed minima and crossing rest on an assumed O(a^2) scaling that the data cannot test.","rationale":"I read the paper as a proceedings summary whose main contribution is a qualitative, continuum-estimated picture of neutral pseudoscalar screening masses in a thermomagnetic medium. The workflow is credible: HISQ/tree action at physical quark masses, three lattice spacings, AICc-based model selection, plateau determination, and bootstrap errors. The weakest step is the continuum extrapolation. The reader's weakest_assumption identifies exactly this: the linear and quadratic fits in Eqs. (9)-(10) have no residual degrees of freedom, so the O(a^2) scaling is an assumption rather than a demonstrated property. I agree with that identification. My concern is not that the authors did anything improper—this is a standard way to present a continuum estimate in a proceedings—but that the strength of the conclusions (presence and shift of minima, crossing of temperature curves) exceeds what a zero-degree-of-freedom extrapolation can establish on its own. A concrete check using only the two finer lattice spacings, or an additional Ntau=20 ensemble, would settle whether the qualitative features survive. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":7265,"tokens_out":7038,"duration_ms":68903,"concrete_test":"Recompute the continuum estimate for pi0 and K0 at T = 145, 157, and 166 MeV using only the Ntau = 12 and 16 data with the linear ansatz (Eq. 9), and compare the eB-dependence with the published quadratic result that includes Ntau = 8. If the location of the minimum in Fig. 4 shifts by more than about 0.1 GeV^2, or if the minimum disappears, then the O(a^2) assumption is not supported by the two finer lattices and the non-monotonic claim should be regarded as not yet continuum-stable. This test uses only data already shown in the paper (or in companion Ref. [11]) and isolates the scaling assumption from the interpolation and fitting details.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central qualitative claims—the non-monotonic pi0 and K0 screening masses with eB and the crossing of constant-eB curves in temperature—are read off the continuum bands in Figs. 4 and 5. These bands are produced by Eqs. (9)-(10): a linear ansatz using only Ntau=12 and 16, and a quadratic ansatz using Ntau=8, 12 and 16. Neither fit has residual degrees of freedom, so the O(a^2) scaling assumption is not validated by the data. The situation is made worse by the fact that the continuum extrapolation is applied to a 2D B-spline interpolation in (T,eB), not to independent raw points; smoothing at Ntau=8 can propagate into the continuum band. If the discretization error contains an O(a^2 alpha_s) or a taste-violating term with a different Ntau dependence, the location of the pi0/K0 minima in Fig. 4 and the crossing temperatures in Fig. 5 could shift. Because the conclusions are qualitative statements about turning points and crossings, a shift in those features, not just in the central values, is enough to weaken the paper's message. The paper explicitly depends on companion Ref. [11] and provides no raw screening masses or fit parameters, so the extrapolation cannot be independently audited from this text.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper reports a lattice QCD study of the screening masses of neutral pseudoscalar mesons (π0, K0, and the fictitious η_sbar) in (2+1)-flavor HISQ/tree QCD at physical quark masses, for temperatures between 145 and 166 MeV and magnetic fields up to eB = 1 GeV^2. The screening masses are extracted from spatial correlators using multi-state fits with AICc model selection, interpolated in the (T, eB) plane via B-splines, and extrapolated to the continuum using linear (Nτ = 12, 16) and quadratic (Nτ = 8, 12, 16) fits in 1/Nτ^2. The central claims are that the continuum-estimated π0 and K0 screening masses are non-monotonic in eB with minima that shift to smaller eB as T increases, that the η_sbar screening mass decreases monotonically with eB, and that constant-eB curves in T cross, which is interpreted as evidence for a reduction of T_pc with magnetic field.","tokens_in":7534,"tokens_out":6918,"duration_ms":62434,"significance":"If substantiated, these results would provide continuum-estimated, physical-quark-mass evidence for inverse magnetic catalysis in the screening spectrum and would support the picture of T_pc reduction in a magnetic field. The paper uses a sensible lattice setup: a tuned HISQ/tree action, physical light and strange quark masses, AICc model selection with bootstrap errors, and three lattice spacings. The main caveats are that the continuum extrapolation relies on fits with zero residual degrees of freedom, the disconnected quark-line contributions are neglected without finite-temperature justification, and the reported AICc formula is nonstandard. Because the qualitative conclusions are statements about minima and crossings of continuum bands, these issues are directly relevant to the paper's central claims.","major_comments":[{"comment":"The continuum extrapolation has zero residual degrees of freedom: the linear ansatz uses only Nτ = 12 and 16 (two data points, two parameters) and the quadratic ansatz uses Nτ = 8, 12, and 16 (three data points, three parameters). Neither fit can test the assumed O(a^2) scaling. Since the minima in Fig. 4 and the crossings in Fig. 5 are features of the continuum bands constructed from these fits, an unquantified taste-violating or O(a^2 α_s) term could shift these features. In addition, the fits are applied to B-spline interpolated values rather than to independent raw lattice points, and the smoothing factor used at Nτ = 8 can feed into the continuum estimate in a way that is not tested. Please provide per-Nτ screening masses (or a table of fit parameters), a scaling plot, and/or an estimate of the systematic error from the extrapolation ansatz. If these checks are documented in companion Ref. [11], please state explicitly which results are taken from there rather than derived in this proceedings text.","section":"§4.3, Eqs. (9)-(10)"},{"comment":"The reported AICc formula is not the standard one: AICc = 2k - ln(hat L) + (2k^2 + 2k)/(n - k - 1) should contain -2 ln(hat L) rather than -ln(hat L). As written, the criterion changes the relative weight of goodness-of-fit and model complexity, which can alter the selected ansatz and hence the extracted screening masses. Please correct Eq. (8) and confirm that the quoted masses are unchanged under the standard definition of AICc.","section":"§4.1, Eq. (8)"},{"comment":"The neglect of disconnected quark-line contributions is justified by citing Ref. [7], a zero-temperature study. For the neutral π0 and η_sbar channels near T_pc, the disconnected part is expected to be more prominent, and the paper itself attributes the observed behavior to sea-quark effects (inverse magnetic catalysis). Without a finite-temperature estimate of the disconnected contribution, the quoted screening masses may not represent the full physical meson screening masses. Please quantify this systematic or cite a finite-temperature study that supports the claim that the disconnected part is small.","section":"§3.3"}],"minor_comments":[{"comment":"The text says 'The constant temperature curves diverge at higher temperature', but Fig. 5 shows screening mass versus T at fixed eB; this should read 'constant magnetic field curves'.","section":"§5.2"},{"comment":"There is a typo in 'to decrease reduce the fitting error': it should be 'to reduce'.","section":"§4.1"},{"comment":"The text contains 'InfigureFigure4' and 'FigureFigure5'; these should be 'In Fig. 4' and 'Fig. 5'.","section":"§5, first paragraphs"},{"comment":"The phrase 'The simulated temperatures ranges from 145 MeV to 166 MeV' should be 'The simulated temperatures range from ...' or 'The temperature ranges from ...'.","section":"Abstract"},{"comment":"The notation for the fictitious eta meson is inconsistent: it appears as η^0_{s\\bar s} in Eq. (3) but as η^0_{\\bar s s} or η^0_{\\bar{ss}} in figure captions and text; please unify the notation.","section":"§2.1 and figure captions"},{"comment":"Reference [11] is a preprint (arXiv:2501.11262); if it has been published or accepted, please update the citation.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings contribution that relies heavily on companion Ref. [11]. The editor may wish to verify that the companion paper contains the raw screening masses and a validation of the continuum extrapolation; if not, the present text should either include enough numerical detail or soften the qualitative claims. The topic is suitable for the venue and the core idea is worth publishing after the load-bearing concerns about the continuum extrapolation, the AICc formula, and the disconnected diagrams are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a LATTICE2024 proceedings that summarizes the same group's companion paper [11]. The new content is modest: physical quark masses and a continuum estimate for pi0, K0, and eta_s sbar screening masses in magnetic fields, extending their earlier heavier-quark study [5]. The paper says so itself, so I don't hold the lack of novelty against it as a proceedings, but it does mean the standalone value is mainly as a pointer to [11].\n\nWhat the paper does well: the lattice setup is standard HISQ/tree with physical pion and kaon masses, multiple Ntau, and a careful AICc-based model selection for the correlator fits with bootstrap uncertainties. The qualitative picture — pi0 and K0 screening masses dipping then rising with eB, eta_s decreasing monotonically, and crossing of constant-eB curves in T — is plausible and consistent with the group's earlier work at heavier quark masses. The neglect of disconnected diagrams is stated explicitly and attributed to earlier estimates, which is fine for a proceedings.\n\nThe soft spots are in the continuum extrapolation, and they are real but not disqualifying for a proceedings. The linear ansatz uses only Ntau=12 and 16, and the quadratic uses Ntau=8, 12, and 16; neither has residual degrees of freedom, so the O(a^2) scaling is assumed rather than tested. Taste violations or O(a^2 alpha_s) terms could shift the minima in Fig. 4 and the crossings in Fig. 5. The extrapolation is applied to a B-spline interpolation in (T,eB), which can smooth features, and no raw screening masses or fit parameters are given, so the continuum bands cannot be independently audited from this text. Those are real limitations, but the paper is honest about being directly based on [11], so the full systematic checks presumably live there.\n\nWho gets value from this: people working on heavy-ion phenomenology or effective models who want a quick summary of where the lattice numbers point, and lattice practitioners who want to see the group's latest analysis choices. It deserves a serious referee because the physics question is important and the simulation work is careful, even though the conclusions are preliminary. My recommendation: accept as a proceedings contribution, but make the authors state explicitly in the text that the continuum estimates are based on zero-degree-of-freedom fits and that the companion paper is the definitive reference. Strictly on its own merits, I would not cite this proceedings; I would cite [11] instead.","headline":"A useful proceedings summary of the group's physical-quark-mass screening mass project, but the continuum estimate is under-validated and the real results live in the companion paper.","tokens_in":8044,"tokens_out":1902,"would_cite":false,"duration_ms":19448,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc"],"model":"deepseek-v4-flash","headline":"In lattice QCD at physical quark masses, neutral pion and kaon screening masses dip to a minimum, then climb, as the magnetic field grows; the strange eta falls throughout — a pattern tied to inverse magnetic catalysis.","keywords":["lattice QCD","screening mass","pseudoscalar mesons","magnetic field","inverse magnetic catalysis","continuum extrapolation","HISQ action","chiral symmetry restoration"],"falsifier":"Run the same physical setup on a fourth, finer lattice — for instance $N_\\tau = 20$ at the same aspect ratio 4 — and redo the continuum estimate: if the $1/N_\\tau^2$ scaling is right, the new point must land within errors of both the $N_\\tau = 12/16$ linear extrapolation and the $N_\\tau = 8/12/16$ quadratic one. A fourth point that forces a different scaling law would move the continuum screening masses and with them the location and depth of the claimed $\\pi^0$ and $K^0$ minima. A complementary check on the same ensembles is to measure the pseudocritical temperature directly from the chiral condensate and compare it with the crossing temperatures of the fixed-field screening-mass curves.","tokens_in":7080,"feed_emoji":"🧲","tokens_out":23125,"duration_ms":160329,"temperature":0.7,"pith_summary":"This paper seeks to establish how the screening masses of neutral pseudoscalar mesons behave in a quark-gluon plasma that is both hot and threaded by a magnetic field, using (2+1)-flavor lattice QCD with physical quark masses and continuum estimates from three lattice spacings. Its central result is that the neutral pion and kaon screening masses are non-monotonic in the field: they fall as $eB$ grows from zero, pass through a minimum, and then rise, whereas the fictitious $\\eta_{s\\bar{s}}$ meson's screening mass decreases monotonically. In temperature, the constant-field screening-mass curves for the pion and kaon cross one another, a pattern the authors read as a magnetic-field-induced reduction of the pseudocritical temperature, consistent with inverse magnetic catalysis. These results matter because screening masses probe the long-distance structure of the medium near the QCD transition, giving a handle on chiral symmetry restoration that complements the short-distance chiral condensate.","feed_headline":"Magnetic field sends pion, kaon screening masses through a minimum","feed_subtitle":"Neutral meson masses dip then climb near the QCD transition — a sign magnetic fields ease chiral restoration.","key_machinery":"The object that carries the argument is the screening mass itself: the decay rate of the spatial meson correlation function at large separation, which sets the length scale over which a mesonic disturbance propagates in the hot, magnetized medium. On the lattice, the masses are extracted by fitting folded correlators to multi-state sums of non-oscillating and oscillating hyperbolic cosines, with the number of states chosen by the corrected Akaike information criterion; the continuum value is then obtained by assuming that the remaining discretization error scales as $1/N_\\tau^2$, combining a linear fit through $N_\\tau = 12$ and 16 with a quadratic fit through $N_\\tau = 8$, 12, and 16. Conceptually, the interpretation is carried by two competing mechanisms: inverse magnetic catalysis from sea quarks pulls the screening mass down as $eB$ grows, while magnetic catalysis from valence quarks pushes it up, and the paper reads the location of the minimum as the balance point of the two.","core_discovery":"Stated as the authors would state it to a fair reader: in the continuum limit of (2+1)-flavor QCD at physical quark masses and temperatures from 145 to 166 MeV, the screening masses of the neutral pion $\\pi^0$ and the neutral kaon $K^0$ are convex functions of the magnetic field. Starting from $eB = 0$ they decrease to a minimum and then increase, and the minimum moves to smaller $eB$ as the temperature rises; the fictitious strange eta $\\eta_{s\\bar{s}}$ instead falls monotonically with $eB$, with no minimum, a difference the authors attribute to strange valence quarks feeling mainly magnetic catalysis. In temperature, the fixed-field curves for $\\pi^0$ and $K^0$ cross at low temperatures, which the paper interprets as stronger fields lowering the pseudocritical temperature and promoting earlier chiral symmetry restoration, while the $\\eta_{s\\bar{s}}$ curves converge without crossing. The continuum estimates rest on lattices with temporal extents $N_\\tau = 8$, 12, and 16 at aspect ratio 4, using the HISQ/tree action with physical quark masses, magnetic fields up to $eB = 1$ GeV$^2$, and an extrapolation that combines linear and quadratic fits in $1/N_\\tau^2$.","pith_inferences":["A direct test of the balance-point reading would compute the $\\pi^0$ pseudoscalar susceptibility on the same ensembles: the Ward-Takahashi identity ties it to the chiral condensate, so the field strength at which the screening mass bottoms out should coincide with where the condensate's drop is steepest.","Because the magnetic field couples directly to charged states, computing the charged-pion screening mass alongside the neutral one would separate valence- and sea-quark effects more cleanly than neutral channels alone; this paper computes only neutral mesons.","Finer $eB$ spacing near zero and additional temperatures flanking the five used here would test whether the initial downward slope from the origin and the crossing temperatures both survive, or whether the B-spline interpolation is smoothing a sharper structure."],"forward_implications":["If the dip-then-rise is real, the light-meson channel's screening length in the plasma first lengthens and then shortens as the magnetic field grows, meaning the range of mesonic correlations is itself non-monotonic in $eB$ near the crossover.","The crossing of the fixed-field $\\pi^0$ and $K^0$ screening-mass curves implies that a stronger magnetic field brings the medium closer to chiral restoration at the same temperature, i.e. that the pseudocritical temperature decreases with $eB$.","The monotonic decline of the $\\eta_{s\\bar{s}}$ screening mass implies that the strange-quark sector responds to the magnetic field mainly through magnetic catalysis in the investigated window, in contrast to the light-quark sector.","Because the minimum shifts to smaller $eB$ as temperature increases, the non-monotonicity is tied to proximity to the crossover and should be most visible just below the pseudocritical temperature."],"supporting_citations":[{"why":"Supplies the quantized magnetic flux condition used in the setup and the field-induced lowering of the pseudocritical temperature that the crossing interpretation builds on.","marker":"[2]"},{"why":"Defines inverse magnetic catalysis by showing the light quark condensate drops with $eB$, the sea-quark effect the paper invokes for the dip.","marker":"[3]"},{"why":"Reviews magnetic catalysis, the valence-quark effect the paper takes to dominate the $\\eta_{s\\bar{s}}$ behavior.","marker":"[4]"},{"why":"Reports the earlier heavier-quark-mass study whose inverse chiral magnetic effect the paper uses to explain the minimum in the pion channel.","marker":"[5]"},{"why":"Provides the meson screening mass extraction method with multi-state cosh fits and the $f_K$ scale setting used throughout.","marker":"[6]"},{"why":"Gives the Ward-Takahashi identities linking pseudoscalar susceptibilities to chiral condensates and supports neglecting disconnected contributions.","marker":"[7]"},{"why":"Provides the lattice implementation of the external magnetic field through quantized flux.","marker":"[8]"},{"why":"Is the parent work this proceedings is directly based on, containing the full analysis from which the continuum estimates are taken.","marker":"[11]"}],"fun_headline_variants":["Magnetic fields send pion and kaon screening masses through a dip","Pion and kaon screening masses dip then climb as magnetic field grows","Strong magnetic fields make pion, kaon screening masses sag to a minimum","Continuum QCD: magnetic fields bend pion, kaon screening masses to a minimum","Magnetic fields cause pion and kaon screening masses to dip, then surge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the leftover lattice-spacing error falls as the square of the lattice spacing, so the screening masses at $N_\\tau = 8$, 12, and 16 lie on a single smooth curve in $1/N_\\tau^2$; the straight-line fit through the two coarser points and the curved fit through all three have no spare points to verify that scaling, and if taste-breaking or other discretization effects are not purely quadratic in the spacing, the continuum masses and the magnetic-field minima would shift.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic fields send pion and kaon screening masses through a dip","Pion and kaon screening masses dip then climb as magnetic field grows","Strong magnetic fields make pion, kaon screening masses sag to a minimum","Continuum QCD: magnetic fields bend pion, kaon screening masses to a minimum","Magnetic fields cause pion and kaon screening masses to dip, then surge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000633,"raw_usage":{"total_tokens":2988,"prompt_tokens":1079,"completion_tokens":1909,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":695,"completion_tokens_details":{"reasoning_tokens":1809}},"tokens_in":695,"tokens_out":1909,"duration_ms":12852,"temperature":1.0,"reasoning_tokens":1809,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T16:23:21.145180+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same physical setup on a fourth, finer lattice — for instance $N_\\tau = 20$ at the same aspect ratio 4 — and redo the continuum estimate: if the $1/N_\\tau^2$ scaling is right, the new point must land within errors of both the $N_\\tau = 12/16$ linear extrapolation and the $N_\\tau = 8/12/16$ quadratic one. A fourth point that forces a different scaling law would move the continuum screening masses and with them the location and depth of the claimed $\\pi^0$ and $K^0$ minima. A complementary check on the same ensembles is to measure the pseudocritical temperature directly from the chiral condensate and compare it with the crossing temperatures of the fixed-field screening-mass curves.","supporting_citations":[],"review_version":1}