{"id":"99db0cec-bdc8-4b0b-a099-752de5bfdcfb","arxiv_id":"1908.04335","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A parametrized quasiparticle model fitted to finite-magnetic-field lattice QCD thermodynamics predicts that magnetic fields and QCD interactions suppress the shear viscosity and electrical conductivity of the quark-gluon plasma.","lead":"Physicists built a simple parametrized model that mimics lattice QCD results for quark-gluon plasma in a magnetic field, then used it to estimate viscosity and electrical conductivity. The model predicts that both the magnetic field and quark-gluon interactions reduce these transport coefficients.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quasi-particle shear-viscosity expression Eq. (29) does not reduce to the quoted non-interacting limit Eq. (23): with v_av=f_k=g=1 it gives coefficient 3/(5π^2), not 4/(5π^2). The Fig.","rationale":"The reader's CONDITIONAL verdict is reasonable, and I keep it. The weakest assumption identified by the reader, the linear condensate-to-mass proxy M=(M_N/3)⟨qq⟩, is a genuine heuristic limitation, but it affects the size of g and hence the overall magnitude of suppression. The normalization mismatch between Eq. (29) and Eq. (23) is a more direct, internal defect in exactly the quantity the paper headlines: the interaction/non-interaction ratios for shear viscosity. It can be demonstrated algebraically without any external physics assumptions. The paper is transparent about the heuristic, rough-estimation character of the model and does not claim machine-checked proofs, so outright rejection is not warranted. However, the quantitative reduction percentages should be corrected or re-quoted with a consistent baseline before the results are used in phenomenological applications. The proposed check settles whether the effect is material or only cosmetic.","tokens_in":22953,"tokens_out":18852,"duration_ms":189096,"concrete_test":"Evaluate Eq. (29) at v_av=1, f_k=1, g=1 and require equality with Eq. (23). Since equality fails, insert the compensating factor 4/3 (or replace k_av^2 by ⟨k^2⟩=12T^2 f_k^2 in Eq. (25)) and recompute the η/(τ_c T^4) curves and the lower-panel Int./Non-Int. ratios in Fig. 3. If the high-T eta reduction changes from the quoted 50% by more than about 10 percentage points, the paper must report corrected percentages; if the change is smaller, the inconsistency is cosmetic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the interaction-driven reduction of η and σ. For η, the chain from Eq. (25)-(29) is internally inconsistent. Eq. (25) approximates η = v_av^2 k_av^2/(15) τ χ and then uses χ = g(gg+gQ)T^2/π^2 with k_av=3T f_k. Setting f_k=v_av=g=1 gives η = (gg+gQ) 3 T^4/(5π^2) τ. But Eq. (23), the massless spinless result used as the non-interacting baseline, is η = (gg+gQ) 4 T^4/(5π^2) τ. These differ by a factor 4/3, so the text's statement that Eq. (29) reduces to Eq. (23) is not correct. The origin is that k_av^2=(3T)^2 is used in place of the required second moment ⟨k^2⟩=12T^2 in the ∫ k^4/E^2 weight. Because Eq. (31) multiplies all η components by the same prefactor, η‖, η⊥ and η/s inherit the mismatch. The lower panels of Fig. 3 plot interacting values from Eq. (29) against a non-interacting baseline from Eq. (23), so the ratio is not a pure interaction-suppression factor. For σ, Eq. (29) does reduce to Eq. (24), so the flaw is specific to the shear-viscosity part of the headline. This does not remove the physical suppression from g<1, v_av<1 and the magnetic denominators, but it means the quoted percentages are not trustworthy until the normalization is fixed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a quasi-particle model of the quark-gluon plasma in a finite magnetic field, in which the interaction is encoded through temperature- and magnetic-field-dependent factors: an effective degeneracy factor g(T,B), an average-energy rescaling f_E(T,B), an average-momentum rescaling f_k(T,B), and an average velocity v_av(T,B). These factors are tuned so that the model reproduces LQCD results for the number density, pressure, energy density, and entropy density at eB = 0, 0.2, and 0.4 GeV^2. The tuned factors are then inserted into relaxation-time-approximation expressions for the shear viscosity and electrical conductivity, producing parallel and perpendicular components in a magnetic field. The paper's central quantitative claim is that both the magnetic field and the QCD interaction are dominant sources of reduction of these transport coefficients, with reported reductions of 50-98.5% for shear viscosity and 30-90% for electrical conductivity in the temperature range studied.","tokens_in":23319,"tokens_out":7452,"duration_ms":77363,"significance":"If the central claim survived scrutiny, the paper would offer a simple, analytic way to translate LQCD-matched thermodynamics into anisotropic transport coefficients for magnetized QGP, which could be useful for magnetohydrodynamic simulations and for quick phenomenological estimates. The authors are transparent about the heuristic nature of the model, and they compare their results with HRG and NJL estimates, which helps place the work in context. However, the quantitative claim is currently weakened by a normalization mismatch in the shear-viscosity expression and by the fact that the transport suppression is largely inherited from the fitted thermodynamic factors rather than emerging from an independent transport calculation. The paper therefore has value as a rough parametric tool, but its headline percentages and the statement that interaction is a 'dominating source' need to be reformulated or supported by additional analysis.","major_comments":[{"comment":"The claim that Eq. (29) reduces to the non-interacting baseline Eq. (23) is not correct for shear viscosity. Setting v_av = f_k = g = 1 in Eq. (29) gives η = (g_g + g_Q) 3 T^4 / (5π^2) τ_c, whereas Eq. (23) gives η = (g_g + g_Q) 4 T^4 / (5π^2) τ_c. The discrepancy is a factor 3/4 and originates from replacing the second moment ⟨k^2⟩ = 12T^2 with k_av^2 = (3T)^2 in the k^4/E^2 weight. Because Eq. (31) multiplies η_parallel, η_perp, and hence η/s by the same prefactor, the 'Int./Non-Int.' ratios in the lower panels of Fig. 3 do not isolate the interaction effect; they contain a spurious normalization factor. The quoted 50-98.5% reductions are therefore not trustworthy until this normalization is corrected.","section":"§3, Eqs. (23), (25), (29), and Fig. 3"},{"comment":"The factors g(T,B), v_av(T,B), and f_k(T,B) are obtained by requiring that Eqs. (14) and (20) reproduce n_LQCD, P_LQCD, and ε_LQCD point by point. Inserting these same fitted functions into Eqs. (29) and (31) means that the reported suppression of η and σ relative to the Stefan-Boltzmann limits is largely a restatement of the fitted suppression of the thermodynamic quantities. The abstract's statement that 'magnetic field and interaction both are two dominating sources' reducing transport coefficients is therefore, at present, partly a consequence of the construction rather than a falsifiable prediction of the transport calculation. The authors should either demonstrate that some feature of the result is not determined by the fitted factors (for example, the specific interplay between the τ_B denominators and the f_k, v_av factors) or explicitly reframe the claim as an estimate inherited from LQCD-matched thermodynamics.","section":"§2, Eqs. (14), (18), (20) and §3, Eqs. (29), (31)"},{"comment":"The conversion from the LQCD quark condensate to a constituent quark mass, M(T,B) = (M_N/3)⟨q̄q⟩_T, is a heuristic assumption motivated only by analogy to an NJL gap equation with zero current quark mass. No normalization of the condensate or derivation of the factor M_N/3 is provided. This relation determines g(T,B) through Eq. (14), and g(T,B) enters every transport coefficient in Eq. (31), so the quantitative reductions reported in the paper hinge on this assumption. The authors should test the sensitivity of their results to this choice, for example by varying the proportionality constant over a range typical of NJL/PNJL models, and should state clearly that the condensate-to-mass mapping is an unvalidated assumption rather than a QCD-derived relation.","section":"§2, Eqs. (13)-(14)"}],"minor_comments":[{"comment":"The sentence 'Eav(T, B = 0) < 1' is dimensionally inconsistent; the text should say E_av(T, B = 0) < 3T, as is stated correctly for k_av.","section":"§3.2"},{"comment":"The reference 'as shown in the right panel of Fig. (2)' appears to point to the wrong figure; the comparison with HRG and NJL results is displayed in Fig. (4), not Fig. (2).","section":"§3.1"},{"comment":"The notation η1, η2, σ0 is used without explicit identification with η_parallel, η_perp, and σ_parallel/σ_perp from Eq. (31); the caption should define these symbols to avoid confusion.","section":"Fig. 3 caption and legend"},{"comment":"There are several typographical errors, including 'coppied' in §2 and 'Completely vanished' in the Summary; the manuscript would benefit from a careful proofreading pass.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"This is a heuristic quasi-particle paper whose authors are admirably explicit about the limitations of their parametric approach. The main quantitative claim, however, is currently undermined by the normalization mismatch in Eq. (29) and by the circularity of using fitted thermodynamic factors as the sole source of transport suppression. Neither issue is fatal in principle: the normalization can be corrected, and the claims can be reframed as 'LQCD-matched estimates' rather than independent predictions. I therefore recommend major revision rather than rejection, provided the authors address the three major comments above."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read it. First, it does something new: it takes the same group's B=0 quasi-particle mapping of LQCD thermodynamics, extends it to finite B, and then feeds the fitted g(T,B), v_av(T,B), and f_k(T,B) into RTA transport expressions for shear viscosity and electrical conductivity. Second, the shear-viscosity formula has a normalization bug. With all interaction factors set to 1, Eq. (29) gives 3/(5π^2)(g_g+g_Q) τ_c T^4, while the non-interacting baseline Eq. (23) is 4/(5π^2)(g_g+g_Q) τ_c T^4. The paper says Eq. (29) reduces to Eq. (23); it does not. The source is the replacement of ⟨k^2⟩=12T^2 by (k_av)^2=9T^2. This affects η_parallel, η_perp, and η/s, so the ratios in Fig. 3 are not a clean interaction-suppression factor. The electrical-conductivity expression does reduce correctly, so the flaw is specific to shear viscosity.\n\nWhat the paper does well: the parametric mapping of the finite-B LQCD thermodynamics is transparent and honest. The authors state the limitations themselves: the condensate-to-mass conversion M=(M_N/3)⟨q̄q⟩ is a heuristic guided by an NJL gap equation, Landau quantization is neglected, and the model is meant only for order-of-magnitude estimates. The appendix gives a self-contained derivation of the anisotropic transport coefficients in RTA. The comparison with HRG and NJL results is useful context.\n\nThe soft spots beyond the normalization issue: the weakest link is the chain from LQCD condensate to a constituent mass to a number density, which then fixes g(T,B) that later controls transport suppression. This is acknowledged as heuristic, but it means the quantitative reductions are essentially propagated from the fit rather than independently predicted. The relaxation time τ_c is a free parameter, so absolute values of η/s are not predictions; only the ratios of interacting to non-interacting values carry the model's message, and those are distorted by the 4/3 error. There are also no error bars on the LQCD fits, which matters for claims about 50-98% reductions.\n\nThe qualitative conclusion—both magnetic field and interaction suppress transport coefficients—is likely robust, since g<1, v_av<1 and the magnetic denominators all act in the same direction. But the specific percentages should not be quoted until the normalization is fixed.\n\nWho is this for? Heavy-ion phenomenologists who need quick analytic estimates of magnetized QGP transport inputs for MHD simulations. It deserves a serious referee; the idea is useful and the problems are addressable. I would send it to review with a request to fix the normalization, add error estimates, and clarify the strong-field domain of validity.","headline":"Useful heuristic extension of the authors' B=0 quasi-particle mapping to finite B, but a 4/3 normalization error in the shear-viscosity expression undermines the reported suppression percentages.","tokens_in":23895,"tokens_out":4400,"would_cite":false,"duration_ms":40924,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper attempts to map lattice-QCD thermodynamics and quark condensate data at finite temperature and magnetic field onto a quasi-particle parametrization, then claims that magnetic field and QCD interaction are the two dominant…","keywords":["quark-gluon plasma","magnetic field","shear viscosity","electrical conductivity","quasi-particle model","lattice QCD thermodynamics","inverse magnetic catalysis","relaxation time approximation"],"falsifier":"A direct lattice calculation of the total number density $n/T^3$ at finite $B$ would settle the mapping without the condensate proxy: if the $g(T,B)$ extracted from the condensate differs from the one obtained directly from $n$, the transport chain collapses. A second check is to test the magnetic-field functional form itself: for fixed $T$, $\\eta_\\perp$ should fall as $1/[1+(\\tau_c/\\tau_B)^2]$ and $\\eta_\\parallel$ as $1/[1+4(\\tau_c/\\tau_B)^2]$ when $B$ is varied, which a kinetic simulation can verify.","tokens_in":22707,"feed_emoji":"🧲","tokens_out":7799,"duration_ms":75073,"temperature":0.7,"pith_summary":"This paper tries to establish that the quark-gluon plasma's interaction and its exposure to a magnetic field can both be captured by a small set of temperature- and field-dependent quasi-particle parameters, and that those same parameters control the plasma's fluid response. Concretely, it claims that interaction and magnetic field are the two dominant sources that reduce the shear viscosity and electrical conductivity of QGP, with interaction alone cutting $\\eta$ by 50% to 98.5% and $\\sigma$ by 30% to 90% over the temperature range considered. The motivation is practical: once the parameters are fitted to lattice QCD thermodynamics, the model gives quick analytical estimates of any phenomenological quantity in a magnetized, interacting QGP without a full microscopic transport calculation.","feed_headline":"Magnetic field and QCD interaction cut QGP viscosity by ~98 percent","feed_subtitle":"A quasiparticle fit to lattice thermodynamics shows both effects suppress transport coefficients by tens of percent.","key_machinery":"The load-bearing mechanism is the fitted set of quasi-particle parameters $g(T,B)$, $f_E(T,B)$, $f_k(T,B)$, and $v_{\\rm av}(T,B)$, defined so that average energy and momentum are $E_{\\rm av}=3T f_E$ and $k_{\\rm av}=3T f_k$ with $k_{\\rm av}=v_{\\rm av}E_{\\rm av}$. These parameters map the non-interacting massless Stefan-Boltzmann picture onto the interacting lattice picture by matching pressure, energy density, entropy density, and a number density built from the lattice quark condensate via $M=(M_N/3)\\langle\\bar q q\\rangle$. Eq. (31) then carries these parameters into the transport coefficients, multiplying the massless $\\eta$ and $\\sigma$ by the interaction factors and by the magnetic denominators, so the equilibrium fit directly determines how much transport is suppressed in the magnetized plasma.","core_discovery":"The paper's central claim is that the same factors used to map lattice QCD thermodynamics—a degeneracy factor $g(T,B)$, average energy $E_{\\rm av}=3T f_E(T,B)$, average momentum $k_{\\rm av}=3T f_k(T,B)$ and average velocity $v_{\\rm av}=f_k/f_E$—can be inserted into kinetic-theory expressions for the shear viscosity and electrical conductivity. In the resulting formulas, Eqs. (31), interaction enters as multiplicative factors $v_{\\rm av}^2 f_k^2 g$ for $\\eta$ and $v_{\\rm av}^2 g$ for $\\sigma$, while the magnetic field enters through the denominators $1/[1+(\\tau_c/\\tau_B)^2]$ and $1/[1+4(\\tau_c/\\tau_B)^2]$, where $\\tau_B=E_{\\rm av}/(qB)$ is the inverse synchrotron frequency. The two suppression mechanisms act independently: the magnetic field shortens the effective relaxation time, and the interaction shrinks the thermodynamic phase space. With both included, the anisotropic components satisfy $\\eta_\\parallel>\\eta_\\perp$ and $\\sigma_\\parallel>\\sigma_\\perp$, all suppressed relative to the free massless values.","pith_inferences":["If the condensate-to-mass proxy is trustworthy, then $g(T,B)$ can be read as an effective number of available degrees of freedom; a natural cross-check would be to extract $g(T,B)$ from the lattice entropy density alone and compare it with the value obtained from the condensate route.","Since the paper leaves $\\tau_c$ free, one could turn its $\\eta/s$ curves into a prediction for a temperature-dependent $\\tau_c(T,B)$ by fixing $\\eta/s$ near the KSS bound, a concrete input that magneto-hydrodynamic simulations could test.","The authors note that inverse magnetic catalysis is faint in $\\eta$ and $\\sigma$ because the thermal distribution suppresses the mass effect; the natural place to look for the IMC signal is the interaction measure $\\epsilon-3P$ or the bulk viscosity, where a peak shift with $eB$ should appear."],"forward_implications":["Heavy-ion fluid simulations that take the model at face value should use strongly reduced, anisotropic viscosities and conductivities, especially near the transition temperature where suppression is largest.","The inequality $\\eta_\\parallel>\\eta_\\perp$ and $\\sigma_\\parallel>\\sigma_\\perp$ means charge and momentum transport are faster along the magnetic field than across it, so magneto-hydrodynamic descriptions require separate longitudinal and transverse coefficients rather than one scalar value.","Because the suppression grows with magnetic field strength, the model predicts that the plasma's fluid behavior becomes more dissipative in the direction transverse to the field at large $eB$.","The same $T,B$-dependent parameters can be reused to estimate other QGP observables—strangeness enhancement, thermal dilepton and photon rates, heavy-quark diffusion, jet quenching—giving quick non-perturbative approximations."],"supporting_citations":[{"why":"Supplies the finite-$B$ lattice quark condensate data that the model converts into $M(T,B)$ and $n_{\\rm LQCD}$.","marker":"[5]"},{"why":"Supplies the finite-$B$ lattice pressure, energy density and entropy density that the fitted parameters $g,f_E,v_{\\rm av}$ must reproduce.","marker":"[6]"},{"why":"Provides the tensor decomposition of shear viscosity in a magnetic field used to split $\\eta$ into parallel, perpendicular and Hall components.","marker":"[16]"},{"why":"Gives the treatment of viscous flow in a strong magnetic field that supports the $\\eta_1,\\eta_2$ component structure adopted in the appendix.","marker":"[17]"},{"why":"Supplies the massless finite-$B$ relaxation-time formulas for shear viscosity and conductivity that the interacting model modifies.","marker":"[20]"},{"why":"Provides an NJL-model comparison for anisotropic shear viscosity in a magnetic field.","marker":"[21]"},{"why":"Provides the hadron-resonance-gas comparison curves for anisotropic transport at finite $B$.","marker":"[23]"},{"why":"Introduces the zero-field degeneracy-factor parametrization that this paper extends to finite magnetic field.","marker":"[66]"}],"fun_headline_variants":["Magnetic field and QCD coupling reduce QGP viscosity and conductivity","B-field and quark interactions shrink QGP transport coefficients","Both B-field and QCD interaction suppress QGP shear viscosity","Magnetic field plus interaction cut QGP transport coefficients","Magnetic field and QCD interaction lower QGP transport coefficients"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire suppression chain rests on treating the lattice quark condensate as a constituent quark mass through $M=(M_N/3)\\langle\\bar q q\\rangle$ and then requiring the resulting number density to equal $g(T,B)\\times 5.23\\,T^3$; if that linear proxy misrepresents the finite-$B$ lattice data, every reported reduction of $\\eta$ and $\\sigma$ changes.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic field and QCD coupling reduce QGP viscosity and conductivity","B-field and quark interactions shrink QGP transport coefficients","Both B-field and QCD interaction suppress QGP shear viscosity","Magnetic field plus interaction cut QGP transport coefficients","Magnetic field and QCD interaction lower QGP transport coefficients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000632,"raw_usage":{"total_tokens":2923,"prompt_tokens":952,"completion_tokens":1971,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":1888}},"tokens_in":568,"tokens_out":1971,"duration_ms":15299,"temperature":1.0,"reasoning_tokens":1888,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:46:47.586444+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct lattice calculation of the total number density $n/T^3$ at finite $B$ would settle the mapping without the condensate proxy: if the $g(T,B)$ extracted from the condensate differs from the one obtained directly from $n$, the transport chain collapses. A second check is to test the magnetic-field functional form itself: for fixed $T$, $\\eta_\\perp$ should fall as $1/[1+(\\tau_c/\\tau_B)^2]$ and $\\eta_\\parallel$ as $1/[1+4(\\tau_c/\\tau_B)^2]$ when $B$ is varied, which a kinetic simulation can verify.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-$B$ lattice quark condensate data that the model converts into $M(T,B)$ and $n_{\\rm LQCD}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-$B$ lattice pressure, energy density and entropy density that the fitted parameters $g,f_E,v_{\\rm av}$ must reproduce."},{"cited_title":"Lifshitz and L.P","cited_arxiv_id":null,"evidence_quote":"Provides the tensor decomposition of shear viscosity in a magnetic field used to split $\\eta$ into parallel, perpendicular and Hall components."},{"cited_title":"Tuchin, On viscous flow and azimuthal anisotropy of the quark–gluon p lasma in a strong magnetic field J","cited_arxiv_id":null,"evidence_quote":"Gives the treatment of viscous flow in a strong magnetic field that supports the $\\eta_1,\\eta_2$ component structure adopted in the appendix."},{"cited_title":"Shear viscosity and electrical conductivity of relativistic fluid in presence of magnetic field: a massless case","cited_arxiv_id":"1907.11164","evidence_quote":"Supplies the massless finite-$B$ relaxation-time formulas for shear viscosity and conductivity that the interacting model modifies."},{"cited_title":"Ghosh, B","cited_arxiv_id":null,"evidence_quote":"Provides an NJL-model comparison for anisotropic shear viscosity in a magnetic field."},{"cited_title":"Satapathy, S","cited_arxiv_id":null,"evidence_quote":"Introduces the zero-field degeneracy-factor parametrization that this paper extends to finite magnetic field."}],"review_version":1}