{"id":"1fa92bfa-3c7c-4264-9fad-86f08f2d0503","arxiv_id":"1908.00800","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Analytic expressions for quark pressure and magnetic susceptibility of dense quark-gluon plasma in strong magnetic fields are derived within the Field Correlator Method, predicting strong paramagnetism.","lead":"This paper derives formulas for the pressure and magnetic susceptibility of hot, dense quark-gluon plasma in a strong magnetic field. It predicts that this plasma is strongly paramagnetic, more so at higher temperature, and that the results line up with lattice simulations of quantum chromodynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central prediction of a rapidly rising magnetic susceptibility rests on Polyakov-loop and Debye-mass inputs taken from zero-field lattice fits; a 7% B-dependence of L(T), if propagated into Eq.","rationale":"I read the paper in good faith and verified the internal consistency of the main formulas. Equation (21) correctly reduces to the zero-field pressure (9) in the B to 0 limit: the apparent linear term in eB cancels through the Macdonald-function identity z K2'(z) = -z K1(z) - 2 K2(z), and the (eB)^2 expansion of Eq. (21) yields exactly the claimed Eq. (34) for chi_q, including the correct T/n factor after the Landau-level sum is handled. The numerical check in Fig. 1, though not reproduced here, is the right validation for the Euler-Maclaurin approximation behind Eq. (21). The reader's weakest assumption is the one I also find most load-bearing: the inputs L(T) and m_D(T) are taken from zero-field lattice fits and assumed independent of B and mu. The paper is transparent about this in Section VI and cites lattice bounds, but the 15% error estimate is not derived by propagating those uncertainties through Eq. (34). Because the susceptibility claim is quantitative, this gap is enough to keep the verdict CONDITIONAL rather than ACCEPT. I do not see a more severe problem: the derivation is internally consistent, the model limitations are stated, and the lattice agreement is honestly described as qualitative. The verdict should remain CONDITIONAL with the concrete check above as the natural path to confirmation or refutation.","tokens_in":9680,"tokens_out":30758,"duration_ms":292147,"concrete_test":"Recompute chi_q(T) from Eq. (34) at mu_B=0 and mu_B=0.4 GeV for T in [0.16,0.5] GeV, using the central L(T), and then with L(T) shifted by +7% and -7% and m_D(T) shifted by +5% and -5%, as permitted by Refs. [44,45]. If the band of chi_q(T) curves has width less than 15% and preserves a monotone rapid increase, the concern is resolved. If the band is wider or the temperature slope changes sign over part of the range, the quantitative claims in the abstract should be scaled back to qualitative. A complementary check is to compare the mu_B=0, eB=0.2 and 0.4 GeV^2 pressure shift from Eq. (21) point-by-point with lattice data of Refs. [29,31,33]; deviations larger than 20% at T below 0.2 GeV would indicate that the B-independent input assumption fails just where susceptibility grows fastest.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—the pressure shift Delta P in Eq. (21) and the susceptibility chi_q in Eq. (34)—are obtained from the zero-field FCM inputs L(T) and m_D(T), with the explicit assumption (Section VI) that L(T,mu_B,B) approximately L(T), m_D(T,B) approximately m_D(T), and m_D(T,mu) approximately m_D(T). The paper cites lattice bounds of 7%, 5%, and 2% for these approximations and adds them linearly to estimate a total 15% error. This is the weakest step for two reasons. First, chi_q in Eq. (34) has L^n inside the sum, so a relative uncertainty in L is amplified by the thermal average over n; near T_c, where L is small but rapidly rising, a 7% shift in L can change chi_q by more than 7%. Second, the central claim that chi_q increases rapidly with temperature is in practice driven by the T-dependence of L(T) and K0(nMbar/T); if L(T,B) differs from L(T) by the quoted 7%, the slope of chi_q(T) can change enough to weaken a quantitative lattice comparison, which is already only qualitative. The paper's own error estimate is an order-of-magnitude bound, not a propagation through Eq. (34). The result is a valid model prediction, not a confirmed QCD prediction, and the conditional verdict is appropriate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytic Field Correlator Method (FCM) treatment of quark pressure and magnetic susceptibility in deconfined QCD at nonzero baryon density and in a uniform external magnetic field. Starting from the zero-field pressure formula, the authors introduce Landau quantization via Eqs. (15)-(16), sum over Landau levels to obtain Eq. (21), and expand it to obtain the magnetic susceptibility in Eq. (34). The numerical section tests the Landau-level summation against the integral form, studies the magnetic shift of the pressure, and compares the combination Delta of Eq. (38) qualitatively with lattice data. The central claims are that the quark pressure has a simple analytic form in a magnetic field and that the QGP is strongly paramagnetic, with the susceptibility increasing rapidly with temperature and slowly with density.","tokens_in":9926,"tokens_out":4887,"duration_ms":55850,"significance":"If correct, these are the first analytic nonperturbative expressions for dense QGP thermodynamics in a magnetic field, and they give a useful model benchmark for lattice and heavy-ion phenomenology. The manuscript has genuine strengths: the Landau-level summation in Eq. (21) is checked numerically against the partial sum of Eq. (20); the B-to-0 limit of the integral form reproduces the zero-field pressure; and the authors are explicit about their input assumptions and error budget. At the same time, the predictive content comes entirely from zero-field Polyakov-loop and Debye-mass inputs, so the magnetic-field dependence and the susceptibility are model predictions that require a sensitivity analysis before the quantitative claims can be accepted.","major_comments":[{"comment":"The derivation of Eq. (34) from Eq. (21) is not shown. The susceptibility is the central new result, and a reader cannot verify the coefficient Nc/(3 pi^2), the sign, or the cancellation of the 1/(eB) and finite terms that must occur in the quadratic expansion. Please display the expansion explicitly, or provide it in an appendix, including the intermediate terms that cancel.","section":"IV, Eq. (34)"},{"comment":"The assumption that L and m_D are independent of B and mu is load-bearing for the quantitative claim that the susceptibility increases rapidly with temperature. The paper cites 7%, 5%, and 2% lattice bounds and adds them linearly, but Eq. (34) contains L^n inside a thermal sum, so a relative uncertainty in L can be amplified near T_c. A sensitivity scan or an error propagation through Eq. (34) is needed; without it the statement that chi_q rises rapidly with T is not yet robust.","section":"VI and Eq. (34)"},{"comment":"The B-to-0 limit of Eq. (27) is a useful internal consistency check, but it is not an independent verification of the magnetic-field dependence, because the parameters V1(infinity,T) and m_D(T) are tuned to reproduce the zero-field lattice pressure. The paper should state explicitly that the predictive content resides in the B-dependent terms, and should not present the B-to-0 agreement as validation of the magnetic-field effects.","section":"II and III, Eqs. (28) and (21)"}],"minor_comments":[{"comment":"The text states that the difference between Eq. (20) and Eq. (21) is negligible, but Fig. 1 would be much more informative with a relative-difference panel or axis; as printed it is difficult to quantify the accuracy of the Landau-level summation.","section":"V, Fig. 1"},{"comment":"Equation (8) quotes c_D ~ 2 from Ref. [23], while the numerical analysis sets c_D = 1.6; please state the reason for this difference or adjust the notation.","section":"V and Eq. (8)"},{"comment":"The phrase 'Polyakov loops interaction' should read 'Polyakov loop interactions', and the text around Eq. (38) should clarify that Delta is the trace anomaly only if the pressure is treated as isotropic, as the authors themselves note.","section":"VI"},{"comment":"Reference [36] is incomplete, with no article number or journal volume; please provide the full citation.","section":"References"},{"comment":"The definition of the magnetic susceptibility in Eq. (30) contains the conventional factor 1/2; it would be helpful to state explicitly that the same definition is used in the comparison with lattice values, to avoid a factor-of-two ambiguity.","section":"IV, Eq. (30)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a model calculation with clear assumptions and useful analytic results. The main risk is that the central quantitative claims are presented a bit more strongly than the sensitivity analysis supports; a careful expansion of Eq. (34) and an uncertainty propagation would address this. The manuscript is within the scope of the journal if model-based predictions are acceptable, but the authors should avoid implying that the magnetic-field dependence is a direct QCD prediction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nMy take on 1908.00800: the genuinely new piece is combining nonzero baryon density with a magnetic field inside the Field Correlator Method. Earlier FCM papers did B at zero mu or mu at zero B; here you get closed-form analytic expressions for the quark pressure, Eq. (21), and the magnetic susceptibility, Eq. (34). The derivation is mostly transparent, and it passes sensible internal checks: the B to 0 limit of Eq. (21) reproduces their earlier zero-field pressure, and Fig. 1 shows the Landau-level summation approximation matches the integral form. That is real work, and the paper deserves credit for checking those limits.\n\nThe main soft spot is not hidden. Section VI states plainly that L(T, mu_B, B), m_D(T, B), and m_D(T, mu) are replaced by their zero-field values, with lattice-based error bounds of 7%, 5%, and 2%, added linearly to a 15% total error. That is honest. But the estimate is an order-of-magnitude bound, not a propagation through Eq. (34). Since chi_q in Eq. (34) sums over n with L^n inside, a 7% uncertainty in L can be amplified, especially near T_c where L is small and rising. So the central prediction—chi_q increasing rapidly with temperature—is model-dependent at the level of its slope, not just its normalization. The paper's own language ('agreement with lattice data' for the zero-density limit) is appropriately cautious; they call the comparison qualitative, which is right.\n\nA smaller issue: the expansion of Eq. (21) that yields Eq. (34) is asserted rather than shown. I would want the algebra checked in review, but nothing suggests it is wrong.\n\nThe citation pattern is heavily self-referential, but that is typical for a method paper, and they do cite lattice data for the parameter ranges and error bounds. No red flag there.\n\nBottom line: this is a solid, well-scoped extension, not a breakthrough. It gives the community manageable analytic expressions for the magnetic equation of state and susceptibility of dense QGP that can be tested on the lattice. The conditional verdict from the report is appropriate. I would send it to peer review; a good referee should ask for the omitted expansion details and a sensitivity analysis of chi_q to L(B) variation near T_c, but the paper deserves referee time.","headline":"A competent FCM extension to dense QGP in a magnetic field, giving analytic pressure and susceptibility, with honest but under-propagated error estimates; worth refereeing.","tokens_in":10514,"tokens_out":2060,"would_cite":false,"duration_ms":21850,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Mh","25.75.-q"],"model":"deepseek-v4-flash","headline":"This paper derives analytic expressions for the quark pressure and magnetic susceptibility of a dense quark-gluon plasma in a uniform magnetic field, and finds the plasma is strongly paramagnetic, with susceptibility rising rapidly with…","keywords":["quark-gluon plasma","magnetic field","magnetic susceptibility","Polyakov loop","Debye mass","Landau levels","baryon chemical potential","field correlator method"],"falsifier":"A lattice computation of the magnetic susceptibility at fixed nonzero baryon density, using the same zero-field values of the Polyakov loop and Debye mass as inputs, would settle the claim: if the analytically predicted rapid temperature growth and slow density growth disagree with the lattice by more than the estimated roughly 15% at $eB\\sim0.5\\,\\mathrm{GeV}^2$, the independence assumption fails.","tokens_in":9444,"feed_emoji":"🧲","tokens_out":6102,"duration_ms":55413,"temperature":0.7,"pith_summary":"Above the deconfinement transition, at temperatures between about one and three times $T_c$ ($T_c\\sim 0.16\\,\\mathrm{GeV}$), with baryon chemical potential up to $0.5\\,\\mathrm{GeV}$ and magnetic field strength $eB$ up to $0.5\\,\\mathrm{GeV}^2$, this paper claims that the quark pressure of the quark-gluon plasma can be written analytically. Starting from the Field Correlator Method, with the nonperturbative physics carried by the Polyakov loop and a Debye mass that grows linearly with temperature, it replaces quark transverse motion by Landau levels and obtains closed forms for the pressure along the field, Eq. (21), and for the magnetic susceptibility, Eq. (34). The susceptibility is claimed to grow rapidly with temperature and only slowly with baryon density, and the zero-density limit is claimed to agree with lattice data. The result matters because magnetized quark matter appears in non-central heavy-ion collisions and in neutron-star interiors.","feed_headline":"A new formula captures dense quark-gluon plasma magnetism","feed_subtitle":"Closed-form expressions cover strong magnetic fields and nonzero baryon density, and match lattice results at zero density.","key_machinery":"The central objects are the Polyakov loop $L(T)=\\exp(-V_1(\\infty,T)/(2T))$, which encodes color-electric deconfinement dynamics, and the Debye mass $m_D(T)$ from color-magnetic confinement, defined by $m_D=c_D\\sqrt{\\sigma_s(T)}$ with $\\sigma_s\\sim g^4 T^2$. The magnetic field enters through the replacement of transverse quark energy by Landau levels, Eqs. (15)-(16), and the summation over Landau levels converts the series over $n$ into the closed forms (21) and (34). The Macdonald functions $K_0,K_1,K_2$ carry the thermal sums, while the same nonperturbative inputs, $L$ and $m_D$, are taken from the zero-field, zero-density fit to lattice QGP pressure.","core_discovery":"The paper's central claim is that, in the temperature range $1<T/T_c<3$, the quark contribution to the pressure of quark-gluon plasma at baryon chemical potential $\\mu_B<0.5\\,\\mathrm{GeV}$ in a uniform external magnetic field $eB<0.5\\,\\mathrm{GeV}^2$ is given by Eq. (21): a sum over Matsubara modes of Bessel functions $K_0,K_1,K_2$, with the field entering only through Landau-quantized quark energies. Expanding this expression in powers of $eB$ yields a magnetic susceptibility, Eq. (34), proportional to a sum of $K_0(n\\bar M/T)$ weighted by Polyakov-loop powers $L^n$ and $\\cosh(\\mu n/T)$. The resulting QGP is strongly paramagnetic: the pressure grows with the field, the susceptibility increases rapidly with temperature and slowly with density, and at zero density the results match lattice data.","pith_inferences":["If the same analytic susceptibility continues to match lattice data at finite density, it could be used to constrain the equation of state of magnetized neutron-star matter at finite isospin or strangeness chemical potentials.","The paper's main limitation is explicit: the Polyakov loop and Debye mass are assumed independent of $B$ and $\\mu_B$; a lattice measurement of these quantities at $eB\\sim0.5\\,\\mathrm{GeV}^2$ as functions of both variables would show how far the formulas can be pushed.","Because the magnetic field enters only through Landau-level energy shifts, the same summation technique could be applied to other magnetized-plasma observables, such as quark-number susceptibilities or conductivities.","The anisotropic-pressure relation noted as a future topic has measurable consequences for flow patterns in non-central heavy-ion collisions if the magnetization is large enough."],"forward_implications":["At zero baryon density, Eq. (21) reproduces the previously known FCM pressure and agrees with lattice data, anchoring the new finite-density extension.","The magnetic susceptibility (34) is an analytic function of $T$, $\\mu_B$, and the nonperturbative inputs, so it can be evaluated without approximate Landau-level partial sums.","The combination $\\Delta=(\\epsilon-3P_z+\\mu n)/T^4$ reproduces the qualitative lattice features of Ref. [33], including the dependence of the peak position and width on $eB$.","The pressure computed is the longitudinal pressure $P_z$; the paper notes the anisotropic relation $P_x=P_y=P_z-\\boldsymbol{M}\\cdot\\boldsymbol{B}$ and leaves its detailed study to future work.","The results support strong paramagnetism of QGP, in line with the lattice study cited as Ref. [29]."],"supporting_citations":[{"why":"Supplies the FCM quark-gluon thermodynamics with Polyakov loop and color-magnetic confinement that the paper extends to finite density and magnetic field.","marker":"[17]"},{"why":"Derives the dense QGP equation of state and integral representations that the magnetic-field generalization builds on.","marker":"[21]"},{"why":"Gives the Debye mass $m_D=c_D\\sqrt{\\sigma_s(T)}$ and the high-temperature string tension behavior used as the nonperturbative input.","marker":"[23]"},{"why":"Provides the nonperturbative inputs $V_1(\\infty,T)$ and $m_D(T)$ that were adjusted to match lattice pressure at zero field and density.","marker":"[26]"},{"why":"Lattice study reporting strong paramagnetism of QGP that the present susceptibility result is compared with.","marker":"[29]"},{"why":"Lattice computation of the $\\Delta$ combination whose qualitative temperature and field dependence is reproduced.","marker":"[33]"},{"why":"Earlier FCM treatment of QGP thermodynamics in a magnetic field at zero density, giving integral forms and equation-of-state expressions used here.","marker":"[39]"}],"fun_headline_variants":["Quantum chromodynamics plasma becomes strongly paramagnetic in B-fields","New formulas for dense QGP in strong magnetic fields match lattice","Magnetic susceptibility of quark-gluon plasma rises fast with temperature","Field-correlator theory predicts strong magnetism in hot QGP","Dense QGP turns paramagnetic in strong magnetic fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The nonperturbative inputs—the Polyakov loop and the Debye mass—are taken from the zero-field, zero-density lattice fit and assumed not to change with magnetic field or baryon density; if they do change, the predicted pressure shift and susceptibility are altered.","fun_headline_variants_meta":{"raw":{"variants":["Quantum chromodynamics plasma becomes strongly paramagnetic in B-fields","New formulas for dense QGP in strong magnetic fields match lattice","Magnetic susceptibility of quark-gluon plasma rises fast with temperature","Field-correlator theory predicts strong magnetism in hot QGP","Dense QGP turns paramagnetic in strong magnetic fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00103,"raw_usage":{"total_tokens":4311,"prompt_tokens":889,"completion_tokens":3422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":3336}},"tokens_in":505,"tokens_out":3422,"duration_ms":22904,"temperature":1.0,"reasoning_tokens":3336,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:32:29.822513+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice computation of the magnetic susceptibility at fixed nonzero baryon density, using the same zero-field values of the Polyakov loop and Debye mass as inputs, would settle the claim: if the analytically predicted rapid temperature growth and slow density growth disagree with the lattice by more than the estimated roughly 15% at $eB\\sim0.5\\,\\mathrm{GeV}^2$, the independence assumption fails.","supporting_citations":[{"cited_title":"Field Correlator Method for the confinement in QCD","cited_arxiv_id":"1804.08946","evidence_quote":"Supplies the FCM quark-gluon thermodynamics with Polyakov loop and color-magnetic confinement that the paper extends to finite density and magnetic field."},{"cited_title":"Dynamical role of Polyakov loops in the QCD thermodynamics","cited_arxiv_id":"1610.01472","evidence_quote":"Gives the Debye mass $m_D=c_D\\sqrt{\\sigma_s(T)}$ and the high-temperature string tension behavior used as the nonperturbative input."},{"cited_title":"New nonperturbative approach to the Debye mass in hot QCD","cited_arxiv_id":"hep-ph/0604004","evidence_quote":"Provides the nonperturbative inputs $V_1(\\infty,T)$ and $m_D(T)$ that were adjusted to match lattice pressure at zero field and density."},{"cited_title":"Thermodynamics of quark-gluon plasma at finite baryon density","cited_arxiv_id":"1906.08677","evidence_quote":"Lattice study reporting strong paramagnetism of QGP that the present susceptibility result is compared with."},{"cited_title":"The magnetic susceptibility in QCD","cited_arxiv_id":"1312.5070","evidence_quote":"Lattice computation of the $\\Delta$ combination whose qualitative temperature and field dependence is reproduced."},{"cited_title":"Borsanyi, G","cited_arxiv_id":null,"evidence_quote":"Earlier FCM treatment of QGP thermodynamics in a magnetic field at zero density, giving integral forms and equation-of-state expressions used here."}],"review_version":1}