{"id":"fa23c306-05a6-445a-a823-721bb0819627","arxiv_id":"1908.10618","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Magnetized hadron resonance gas with van der Waals excluded volume predicts paramagnetic hadronic matter and interaction-induced suppression of baryon and electric charge susceptibilities.","lead":"Scientists combined a hadron resonance gas model with repulsive interactions and a constant magnetic field to calculate the thermodynamics of hot dense hadronic matter. They found that the matter is paramagnetic, and that repulsive interactions barely change the magnetization but strongly suppress conserved charge fluctuations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paramagnetic claim rests on the vacuum-pressure sign, and the printed MFIR derivation of that vacuum pressure is internally inconsistent.","rationale":"The reader's weakest assumption was the excluded-volume extension in Eq. (23), which is a legitimate concern for the quantitative susceptibility predictions. However, the central qualitative claim — paramagnetic hadronic matter with negligible interaction effect — is set by the vacuum part of the magnetization, since Fig. 3(a) is essentially T-independent and identical between HRG and EVHRG. Therefore the first, most load-bearing link is the regularized vacuum pressure. The printed derivation contains a sign mismatch between Eq. (10) and Eq. (11): the displayed square root has the wrong sign relative to Eq. (5), yet the zeta-function evaluation assumes the correct sign. This is exactly the type of internal inconsistency that should be resolved before accepting the headline conclusion. The paper has no machine-checked derivation, no external benchmark for Eq. (22), and no reproducible code, so an independent algebraic check is warranted. If the vacuum-pressure sign is confirmed correct, the qualitative paramagnetic claim likely stands; if it flips, the central claim fails. The excluded-volume anisotropy concern remains secondary and would affect quantitative susceptibilities rather than the sign of M.","tokens_in":17814,"tokens_out":26757,"duration_ms":285880,"concrete_test":"Independently redo the step from Eq. (9) to Eq. (11) using the plus-sign Landau dispersion of Eq. (5), then expand the resulting Eq. (22) for eB << m^2. Check whether the B^4 coefficient reproduces the standard Euler-Heisenberg vacuum limit and whether the sign of ΔP_vac at x = 2 matches Fig. 1(a). If Eq. (11) is not the correct regularization of Eq. (9), or if the sign of ΔP_vac changes, the magnetization sign in Fig. 3 is not established.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim — positive magnetization that is insensitive to excluded volume — is supported in Fig. 3 by a total magnetization almost completely dominated by the vacuum contribution. Hence the load-bearing calculation is the regularized vacuum pressure of Sec. 3, Eqs. (8)–(22). That derivation is internally inconsistent as printed: Eq. (10) writes sqrt(p_z^2 + m^2 − 2eBn), but the Landau dispersion for a spin-1/2 charged hadron in Eq. (5) is sqrt(p_z^2 + m^2 + 2eBn) for the lowest spin-up level, and Eq. (11) evaluates the integral with the plus-sign form, since its Hurwitz-zeta argument x = m^2/(2eB) would not arise from the printed minus sign. The final expression Eq. (22) inherits all sign choices made at this step. If any sign in ΔP_vac is wrong, M_vac = ∂ΔP_vac/∂B can change sign, which would remove the paramagnetic conclusion. An independent evaluation of Eq. (22) at moderate x (e.g., x = 2, corresponding to eB ≈ m^2/4) gives a negative ΔP_vac for spin-1/2, in tension with the statement that the vacuum pressure is positive for a wide range and with Fig. 1(a). This is not the excluded-volume issue: Eq. (23) changes only the small thermal part of M and cannot rescue a wrong vacuum sign.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the hadron resonance gas (HRG) model to a constant magnetic field, regularizes the charged-hadron vacuum pressure with the magnetic-field-independent regularization (MFIR) scheme, and adds van der Waals excluded-volume interactions with hard-core radii taken from an earlier fit to hadron yields. It computes the pressure, energy and entropy densities, magnetization, and baryon- and electric-charge susceptibilities of hadronic matter up to T ≈ 0.16 GeV at zero and finite baryon chemical potential, comparing the equation of state with lattice data of Ref. [77] up to T ≈ 140 MeV. The central claims are: (i) low-temperature hadronic matter is paramagnetic, with the sign of the magnetization set by the regularized vacuum contribution; (ii) repulsive excluded-volume interactions have negligible effect on the total magnetization; and (iii) both the magnetic field and the excluded-volume repulsion suppress the conserved-charge susceptibilities and modify the susceptibility ratios.","tokens_in":18081,"tokens_out":29792,"duration_ms":267523,"significance":"If the claims hold, the paramagnetism of the low-temperature hadronic phase and its insensitivity to repulsive interactions are clean, falsifiable predictions that can be confronted with lattice QCD results for the equation of state in magnetic fields; the same model reproduces the available lattice pressure data up to T ≈ 140 MeV both with and without B. The paper deserves credit for not fitting any parameter to its target observables: the hard-core radii are taken from Ref. [44], the gyromagnetic ratios are the tree-level values g_i = 2, and eB = 0.2 GeV^2 is a representative choice, so the results are genuine model predictions rather than fits. The vacuum-pressure formulas for spins 0, 1/2 and 1 are given explicitly in the appendices, and the susceptibility ratios χ4/χ2 and χ6/χ4 are concrete outputs that can be compared with future lattice data. The main caveats are that the printed vacuum-pressure derivation contains internal sign inconsistencies that must be corrected before the central claim can be trusted, and that the susceptibility conclusions rest on an assumed, rather than derived, extension of the excluded-volume scheme to magnetic fields.","major_comments":[{"comment":"The printed derivation of the vacuum pressure is internally inconsistent. Eq. (10) contains sqrt(p_z^2 + m^2 - 2eBn), which becomes imaginary for n > m^2/(2eB) and is inconsistent both with the Landau dispersion in Eq. (5), which has a plus sign for the lowest spin-up level, and with the Hurwitz-zeta evaluation in Eq. (11), whose argument x = m^2/(2eB) can only arise from the plus-sign form sqrt(p_z^2 + m^2 + 2eBn). In addition, Eq. (9) subtracts E_{p,0}/2 inside the sum while Eq. (10) subtracts E_0 once, so the lowest-Landau-level subtraction differs by a factor of two between the two equations. Since Eq. (22) determines the sign of the vacuum magnetization, which is the basis of the paper's headline paramagnetic claim, this step must be corrected and re-derived, or explicitly identified as typographical with the corrected intermediate steps supplied.","section":"Sec. 3, Eqs. (9)-(11)"},{"comment":"The sign of the spin-1/2 vacuum pressure given by the paper's own formula appears to contradict the text. Evaluating Eq. (22) at x = 2 gives Delta P_vac ≈ -0.468 (eB)^2/(2π^2) ≈ -0.024 (eB)^2, i.e., a negative value, and the value at the proton-relevant x ≈ 2.2 is likewise negative, whereas Sec. 5 states that the vacuum pressure is positive for a wide range of magnetic fields for all three spin channels. Because Fig. 1(a) uses a logarithmic vertical scale it cannot display negative values, so the figure should state what is actually plotted, and the sign of Delta P_vac, and of M_vac = ∂(Delta P_vac)/∂B, should be given explicitly over the plotted range of x; this is the quantity that supports the paramagnetic conclusion.","section":"Sec. 3, Eq. (22) and Fig. 1(a)"},{"comment":"The extension of the excluded-volume scheme to finite magnetic fields is assumed rather than derived. In a magnetic field the ideal-gas pressure is anisotropic, with P_parallel = -Ω and P_perpendicular = -Ω - MB, whereas Eq. (23) is a scalar transcendental equation using the isotropic P_id, and the excluded-volume parameter v is taken to be B-independent. This equation is the sole input generating the interaction effects on the susceptibilities in Figs. 5-8, so the paper should either justify the scalar form from a microscopic starting point or clearly state that Eq. (23) is an ansatz and estimate the resulting uncertainty; a concrete diagnostic would be to compare the solution of Eq. (23) obtained with P_parallel and with P_perpendicular in place of P_id.","section":"Sec. 4, Eq. (23)"},{"comment":"The Landau spectrum and the degeneracy factor are written with the signed charge e_i throughout, e.g., E_{i,c}^2 = p_z^2 + m_i^2 + 2 e_i B (n + 1/2 - S_z) and the prefactor eB/(2π)^2. For the negatively charged hadrons in Table 1 (π^-, K^-, ρ^-, K*^-, Σ^- and antiprotons) this would give imaginary energies at large n and a negative degeneracy factor; the physical spectrum and degeneracy involve |e_i|B. The text should state explicitly that |e_i|B is used throughout; as printed, the model equations are only defined for positively charged species.","section":"Sec. 2, Eqs. (2)-(6)"}],"minor_comments":[{"comment":"The text says n is 'any positive integer', but the lowest Landau level n = 0 is used throughout, for instance in Eq. (9); n should be described as a non-negative integer.","section":"Sec. 2, Eq. (5)"},{"comment":"The Boltzmann factors contain unmatched parentheses, with '(exp [Ej - µi)/T ]))' appearing twice; the intended expression is exp[(E_j - μ_i)/T], and the sign convention for fermions and bosons should be stated explicitly.","section":"Eq. (25)"},{"comment":"The caption contains a duplicated '(b)' label in the sub-caption for the energy-density panel.","section":"Fig. 4 caption"},{"comment":"The summary states that 'all the thermodynamical quantities are strongly suppressed due to non-zero background magnetic field', which is not consistent with Fig. 1(b), where the vacuum contribution makes the total pressure at low T larger for eB = 0.2 GeV^2 than for eB = 0; the statement should be restricted to the thermal parts of the thermodynamic quantities.","section":"Sec. 6"},{"comment":"Describing the positive sign of the magnetization as 'a fundamental characteristic of thermal QCD vacuum' overstates what an HRG model calculation can establish; a more cautious formulation is advised, since the result is a property of the MFIR-regularized hadron resonance gas model.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"The printed vacuum-pressure derivation in Sec. 3 appears to contain sign typos; the final expressions (22), (B.1) and (B.2) match the known MFIR results in Refs. [63,67,68], so I regard this as a fixable revision rather than grounds for rejection, but the authors must correct the intermediate equations and verify the sign statements in Fig. 1(a). The excluded-volume-in-magnetic-field assumption should be handled as a stated ansatz. The incremental contribution over Ref. [63], which already contains MFIR vacuum terms in the HRG context, lies in the excluded-volume extension and susceptibility analysis; the authors should make this novelty clearer. The citation concentration on the authors' own previous work is noticeable but not inappropriate given the continuity of the research line."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first combination I know of the excluded-volume HRG with a magnetic-field HRG, and the qualitative claims—positive magnetization that is insensitive to excluded volume, magnetic suppression of baryon and charge susceptibilities—are the sort of baseline the field can use. It is not a breakthrough, and the quantitative results should not be used until a sign problem at the center of the derivation is resolved.\n\nWhat's good: the model has no fitted parameters aimed at the target observables. The hard-core radii come from an earlier yield fit, g=2 is the tree-level value, and eB=0.2 GeV^2 is representative. The MFIR regularization follows Endrődi and earlier work, and the particle list is standard. The paper is also clear that no experimental comparison is possible yet because there is no parametrization of B with collision energy, and lattice results for the susceptibilities in a field are not available. That is honest.\n\nThe soft spots are real. The printed derivation of the vacuum pressure, Eqs. (10)–(11), contains a sign inconsistency: the dispersion in Eq. (5) for the lowest spin-up Landau level is p_z^2 + m^2 + 2eBn, but Eq. (10) writes sqrt(p_z^2 + m^2 − 2eBn). The Hurwitz-zeta form in Eq. (11) implicitly uses the plus sign, since x = m^2/(2eB) would not arise from the minus sign. So as printed the derivation cannot be followed. This matters because Fig. 3 shows the total magnetization is almost entirely the vacuum contribution; the thermal part is two orders of magnitude smaller. If the sign of ΔP_vac is wrong, the paramagnetic claim goes with it. My own quick evaluation of the printed Eq. (22) at x=2 gives a negative vacuum pressure, which conflicts with the text's claim that it is positive for a wide range and with the log plot in Fig. 1(a). This might be a typo or a sign convention issue, but it has to be fixed and checked against Endrődi's result before the physics claim can be trusted.\n\nThe second soft spot is the excluded-volume extension in Eq. (23), which uses the isotropic ideal pressure in a magnetic field where the pressure is anisotropic. The paper does not justify this or explain why the excluded-volume parameter v should be B-independent. That concern mostly affects the thermal part and hence the quantitative susceptibility suppression; it does not rescue or destroy the vacuum part.\n\nBottom line: this is a useful paper for people building hadronic equations of state in magnetic fields, but as written the key derivation is inconsistent. A serious referee could sort it out; I would send it out rather than desk reject. If the sign check fails, the paper is mainly a null/incomplete result; if it passes, the model is a solid baseline.","headline":"A useful hadronic EoS baseline in a magnetic field, but the vacuum-pressure sign that drives the paramagnetic claim is inconsistent as printed.","tokens_in":18635,"tokens_out":4943,"would_cite":false,"duration_ms":45273,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Mh","12.39.-x","11.30.Rd","11.30.Er"],"model":"deepseek-v4-flash","headline":"This paper claims that the low-temperature hadronic phase of QCD remains paramagnetic even when short-range repulsive interactions between hadrons are included, and that excluded-volume repulsion barely changes the magnetization.","keywords":["hadron resonance gas","excluded volume","magnetic field","magnetization","paramagnetic matter","conserved charge susceptibilities","Landau levels","QCD phase diagram"],"falsifier":"Compute the excluded-volume pressure using the parallel and perpendicular pressures in place of the isotropic ideal pressure in the self-consistent equation and check whether the magnetization remains positive; alternatively, measure the magnetization of QCD matter on the lattice at $eB\\approx 0.2$ GeV$^2$ and temperatures below 160 MeV and see whether it is positive and nearly independent of hadron hard-core radii.","tokens_in":1597,"feed_emoji":"🧲","tokens_out":4319,"duration_ms":87950,"temperature":0.7,"pith_summary":"This paper extends the hadron resonance gas model to a constant external magnetic field with short-range repulsive interactions between hadrons, implemented through an excluded-volume correction to the pressure. It asks whether the magnetic response of low-temperature hadronic matter survives once repulsion is included. The authors find that the magnetization of hadronic matter is positive, so the hadronic phase is paramagnetic, and that repulsive interactions have only a negligible effect on the magnetization. They also find that both the magnetic field and excluded-volume repulsion suppress the baryon and electric-charge susceptibilities, with the suppression stronger for higher-order fluctuations. If the paper is right, the hadronic phase below roughly 160 MeV is paramagnetic even when hadron sizes are taken into account.","feed_headline":"Hadron gas stays paramagnetic even with repulsive cores","feed_subtitle":"A magnetic-field resonance gas model finds positive magnetization survives excluded-volume repulsion below 160 MeV.","key_machinery":"The load-bearing object is the excluded-volume pressure equation $P_{\\mathrm{EV}}(T,\\mu,B) = P_{\\mathrm{id}}(T, \\mu - v P_{\\mathrm{EV}}, B)$, a self-consistent transcendental equation that shifts the chemical potential by the excluded-volume pressure. The ideal pressure for charged species is a sum over Landau levels and spin projections, while neutral species use the free-particle dispersion. Vacuum divergences are removed by magnetic-field-independent regularization, which separates the magnetic-field-dependent vacuum pressure from the zero-field vacuum pressure and yields renormalized expressions for spin-0, spin-1/2, and spin-1 hadrons. This regularized vacuum term supplies the dominant positive magnetization, while the thermal term carries the temperature dependence and the excluded-volume correction carries the repulsive interaction.","core_discovery":"The central claim is that an interacting hadron resonance gas in a constant magnetic field produces positive magnetization for low-temperature hadronic matter. The vacuum contribution, regularized by a magnetic-field-independent scheme, dominates and fixes the sign; the thermal contribution is negative at low temperature once pions populate but turns positive when spin-1 and spin-1/2 hadrons contribute. Adding excluded-volume repulsion suppresses the pressure, energy density, entropy density, and all studied baryon and electric-charge susceptibilities, but changes the magnetization by only a negligible amount, leaving the paramagnetic character intact. At finite baryon chemical potential, $\\mu_B = 300$ MeV, the same qualitative picture holds.","pith_inferences":["If the positive vacuum magnetization is dominant, hadron-gas magnetization is largely fixed by the vacuum sector, so a lattice QCD measurement of magnetization below the transition temperature could directly test whether the vacuum and thermal balance is as the model predicts.","A testable extension would replace the isotropic pressure in the excluded-volume equation with the parallel or perpendicular pressure of the magnetized system; this could change the size of the repulsive correction, though it would likely leave the sign of the magnetization unchanged.","The model's susceptibility ratios suggest an experimental route: comparing higher-order fluctuation ratios between central and peripheral heavy-ion collisions, where the magnetic field differs, could probe magnetic-field effects if freeze-out parameters are controlled.","The assumption that the excluded-volume parameter is independent of the magnetic field could be relaxed, and a field-dependent hard-core radius would introduce a new source of suppression that could be checked against future lattice results."],"forward_implications":["Below about 160 MeV, the hadronic phase should be paramagnetic, with positive magnetization, even when hadron hard-core repulsion is included.","Excluded-volume repulsion suppresses pressure, energy density, and entropy density, and the suppression grows with temperature and becomes stronger at finite baryon chemical potential.","Baryon and electric-charge susceptibilities of all studied orders are reduced by both the magnetic field and repulsive interactions, with higher-order susceptibilities suppressed more strongly.","Ratios such as $\\chi^6_B/\\chi^4_B$ drop sharply in the interacting model, making higher-order fluctuation ratios a sensitive probe of repulsion.","The magnetic field suppresses the thermal population of spin-0 hadrons by increasing their effective mass, while spin-1 hadrons become lighter, and this shapes the thermal contribution to magnetization."],"supporting_citations":[{"why":"Supplies the scattering-phase basis for treating the hadron spectrum as a noninteracting resonance gas.","marker":"[39]"},{"why":"Supplies the magnetic-field single-particle dispersion and the renormalization setup for charged hadrons in a magnetic field.","marker":"[63]"},{"why":"Introduce the magnetic-field-independent regularization scheme used to separate and renormalize the vacuum pressure.","marker":"[67, 68]"},{"why":"Supplies the thermodynamically consistent excluded-volume equation used to incorporate repulsive interactions.","marker":"[73]"},{"why":"Provides the hadron hard-core radii used for pions, kaons, mesons, and baryons.","marker":"[44]"},{"why":"Provides the lattice QCD pressure data against which the model's pressure is compared.","marker":"[77]"},{"why":"Establishes the excluded-volume hadron resonance gas baseline for baryon and charge susceptibilities that this paper extends to finite magnetic field.","marker":"[48]"}],"fun_headline_variants":["Paramagnetism survives hadron repulsion in magnetic fields","Repulsive hadron cores can't flip positive magnetization","Excluded volume barely dents hadron gas paramagnetism","Hadron gas magnetization stays positive under repulsion"],"cache_read_input_tokens":20736,"weakest_assumption_plain":"The argument hinges on assuming that the excluded-volume equation, built from the isotropic ideal pressure, remains valid inside a magnetic field with a single magnetic-field-independent excluded-volume parameter, even though the pressure in a magnetic field is anisotropic.","fun_headline_variants_meta":{"raw":{"variants":["Paramagnetism survives hadron repulsion in magnetic fields","Repulsive hadron cores can't flip positive magnetization","Excluded volume barely dents hadron gas paramagnetism","Hadron gas magnetization stays positive under repulsion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000141,"raw_usage":{"total_tokens":1160,"prompt_tokens":933,"completion_tokens":227,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":164}},"tokens_in":549,"tokens_out":227,"duration_ms":3018,"temperature":1.0,"reasoning_tokens":164,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:38:50.491907+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the excluded-volume pressure using the parallel and perpendicular pressures in place of the isotropic ideal pressure in the self-consistent equation and check whether the magnetization remains positive; alternatively, measure the magnetization of QCD matter on the lattice at $eB\\approx 0.2$ GeV$^2$ and temperatures below 160 MeV and see whether it is positive and nearly independent of hadron hard-core radii.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the magnetic-field single-particle dispersion and the renormalization setup for charged hadrons in a magnetic field."},{"cited_title":"Peskin and D","cited_arxiv_id":null,"evidence_quote":"Provides the lattice QCD pressure data against which the model's pressure is compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the excluded-volume hadron resonance gas baseline for baryon and charge susceptibilities that this paper extends to finite magnetic field."}],"review_version":1}