{"id":"aa4a2fcd-498c-454a-8ae3-df8f12b52292","arxiv_id":"2411.11398","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"New analytic and numerical solutions for 1+1D viscous magnetohydrodynamics show that viscosity and a decaying magnetic field heat the fluid and produce an early temperature peak.","lead":"This paper derives new mathematical formulas for how the temperature of the quark-gluon plasma created in heavy-ion collisions changes over time when both viscosity and a magnetic field are present. The formulas could serve as simple benchmark solutions for complex computer simulations of these collisions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (50) contains a spurious homogeneous term: it violates the initial condition T(τ0)=1 by ε1 and is not a systematic O(ε1) perturbative solution of Eq. (41).","rationale":"The paper's central contribution is the pair of analytical solutions Eqs. (40) and (50) for first-order viscous non-resistive MHD. The reader's conditional verdict already focuses on Eq. (50)'s crude approximation and limited regime. My concern is more specific: Eq. (50) does not satisfy the initial condition at O(ε1), because Eq. (49) contains an unsourced homogeneous mode. This makes Eq. (50) not a systematic perturbative solution of the initial-value problem, independent of how small ε1 is. The numerical comparisons in Fig. 5 are therefore contaminated by a known initialization error, which strengthens the reader's call for revision but does not invalidate the paper entirely: Eq. (40) appears to track the numerics, and the paper itself concedes Eq. (50)'s limitations. The appropriate verdict remains CONDITIONAL, as the reader already recommended, so no change to the verdict is needed. I partially agree with the reader because they identified the approximation but not the explicit initial-condition violation and the order-mixing problem in Eq. (47).","tokens_in":18934,"tokens_out":26378,"duration_ms":284843,"concrete_test":"Set ε2=0 and a=2 so that T0=(τ0/τ)^{1/3} exactly. Solve T1' + T1/(3τ) = (τ0/τ)^3 / τ with T1(τ0)=0; the unique solution is T1 = (3/8)(τ0/τ)^{1/3} − (3/8)(τ0/τ)^3. Compare with Eq. (49), which gives T1 = (11/8)(τ0/τ)^{1/3} − (3/8)(τ0/τ)^3; the difference is (τ0/τ)^{1/3}. Then plot Eq. (50) with this homogeneous term removed against the numerical solution of Eq. (53) for ε1=0.05, ε2=1, a=2. If the corrected curve removes the initial offset and reduces the late-time deviation, the flaw is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The second perturbative solution, Eq. (50), is not a legitimate solution of the stated initial-value problem. In Sec. III D, the first-order equation Eq. (48) has general solution T1 = 3/(4−6a) u^{2a−1} + C u^{1/3}, with u = τ0/τ. Requiring T1(τ0)=0 fixes C = −3/(4−6a). Equation (49) instead has C = (1−6a)/(4−6a), which gives T1(τ0)=1 and hence T_NS-2(τ0)=1+ε1. Equivalently, Eq. (50) is the correct zero-initial-value solution plus the unsourced homogeneous term ε1 u^{1/3}. This is not a small-regime subtlety: at the initial time the error is exactly ε1, independent of τ, so the 5% deviation reported at ε1=0.05 and the larger deviations in Fig. 5 are partly an initialization artifact of Eq. (50), not purely a truncation error. The derivation also mixes orders: Eq. (47) keeps ε1T1 in the denominator and then approximates the denominator by τ0/τ, instead of dropping the ε1T1 term before linearization. Consequently, the paper's central claim that Eqs. (40) and (50) are two valid perturbative NS-MHD solutions is overstated; Eq. (50) should be corrected to preserve the initial condition before it can be used as a benchmark.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies 1+1 dimensional relativistic non-resistive magnetohydrodynamics with longitudinal boost invariance, a power-law decaying transverse magnetic field, and shear viscosity. It derives analytical solutions for the ideal MHD (Victor-Bjorken) case, the first-order Navier-Stokes viscous case with and without a magnetic field, and presents two perturbative analytical solutions for the combined case (Eq. (40) and Eq. (50)). These are compared with direct numerical integration of the energy-conservation equation. The paper also solves the second-order Israel-Stewart equations numerically and studies the effects of initial shear stress on the temperature evolution. The central physical claims are that both the magnetic field and shear viscosity heat the fluid, producing an early temperature peak, and that in the Israel-Stewart theory the initial shear stress significantly affects the cooling rate.","tokens_in":19213,"tokens_out":20623,"duration_ms":148237,"significance":"If the results are correct, the paper provides useful closed-form benchmarks for validating numerical codes in a simplified setting relevant to heavy-ion collisions. The derivations are mostly explicit and the numerical comparisons are straightforward, which makes the paper valuable as a pedagogical or reference contribution. However, the significance is moderated by two issues: the second perturbative solution Eq. (50) does not satisfy the stated initial condition, and the Israel-Stewart energy-conservation equation (59) contains an algebraic error in its shear-stress coupling. The first perturbative solution Eq. (40) is the strongest result: it appears correct, reduces properly to known limits, and agrees well with numerics for small magnetic-field parameters. The paper also honestly documents the limited accuracy of Eq. (50), though the root cause (an initialization artifact rather than purely truncation error) is not diagnosed.","major_comments":[{"comment":"The second perturbative solution does not satisfy the initial condition T̃(τ0)=1. The general solution of Eq. (48) is T̃1 = 3/(4−6a) u^{2a−1} + C u^{1/3} with u = τ0/τ. Requiring T̃1(τ0)=0 fixes C = −3/(4−6a), whereas Eq. (49) uses C = (1−6a)/(4−6a), giving T̃1(τ0)=1 and hence T̃NS-2(τ0)=1+ϵ1. The difference is exactly the unsourced homogeneous term ϵ1(τ0/τ)^{1/3}. This is not a truncation error: the deviation at τ=τ0 is exactly ϵ1, independent of τ, so the 5% discrepancy at ϵ1=0.05 in Fig. 5(b) is partly an initialization artifact. The derivation of Eq. (47) also mixes orders by retaining ϵ1T̃1 in the denominator and then approximating it away; a consistent first-order truncation would give the corrected T̃1 = 3/(4−6a)[(τ0/τ)^{2a−1} − (τ0/τ)^{1/3}]. Please correct Eq. (50) and the limits in Eqs. (51)–(52) accordingly, or explicitly state that this solution does not obey the initial condition.","section":"§III D, Eq. (50)"},{"comment":"The energy-conservation equation for the Israel-Stewart case is inconsistent with the preceding derivation. Starting from Eq. (54), using ε=3p, p=a1T^4, B²=σT0^4(τ0/τ)^{2a}, and π=a1b3T0^4π̃, one obtains after dividing by 12a1T0^4T̃^3: ∂T̃/∂τ = −T̃/(3τ) + ϵ1/T̃³(τ0/τ)^{2a}/τ + b3π̃/(12τT̃³). Equation (59) prints the last term as + b3π̃T̃/(12τ), which is a different dynamical term. Since Eqs. (59)–(60) are the basis of the numerical IS results in Figs. 6 and 7, either the printed equation is misprinted or the numerics solved a different equation; please state which and correct the equation.","section":"§IV B, Eq. (59)"}],"minor_comments":[{"comment":"The text after Eq. (40) states \"this solution is stable where ϵ1 is small. We will prove this in the next section IV A,\" but Sec. IV A supplies only a numerical comparison, not a stability proof. The stability claim in the abstract and conclusions is therefore supported by numerics, not by a formal analysis; please rephrase to avoid the word \"prove\" or add a genuine perturbative stability argument.","section":"§III C, §IV A"},{"comment":"The derivation of Eq. (47) is not a systematic first-order expansion: the denominator (T̃0³+3ϵ1T̃0²T̃1) is kept before being approximated away, which obscures the order counting. The final result is equivalent to the standard first-order truncation, but the presentation should be made more transparent.","section":"§III D, Eq. (47)"},{"comment":"For a < 2/3 the exponent 2a−4/3 is negative, so the factor (τ0/τ)^{2a−4/3} grows with time; the solution may develop a negative argument in the fourth root for sufficiently large τ. The paper should state the domain of validity of Eq. (19) for these parameter values, which are used in Fig. 1.","section":"§III A, Eq. (19)"},{"comment":"The expression for x1 in Eq. (39) is very long and involves hypergeometric functions, making it difficult to verify. Since the figures only use a=2, reporting the simplified a=2 limit would improve readability and facilitate checks.","section":"§III C, Eq. (39)"},{"comment":"There are several language and typographical errors, e.g., \"quluon\" for \"gluon\" in Sec. I, \"fouce\" for \"focus\" in Sec. II, \"follow\" for \"follows\", \"regress\" for \"revert\" in Case-C of §III D, and \"consistented\" for \"consistent\" in Sec. V. A careful proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The primary contribution of this paper is the first perturbative solution (Eq. (40)) and the explicit ideal/viscous limits, which appear correct and could serve as useful benchmarks. The errors in Eq. (50) and Eq. (59) are local and fixable, so rejection is not warranted. However, as written, the abstract's claim of \"two perturbative analytical solutions\" is overstated, and the IS numerical results may be based on an incorrect equation. The authors should be asked to correct these points and to report updated comparisons with numerics after the fix."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new part of this paper is the first-order-in-ε1 solution, Eq. (40): a closed-form temperature profile for 1+1D boost-invariant Navier-Stokes MHD with a power-law magnetic field. The construction (nonconserved-charge method, with the magnetic term treated as a perturbation on Muronga's viscous Bjorken flow) is transparent, and the numerical comparison in Fig. 5 shows it tracks the ODE solution well for ε1 up to about 0.1. The ideal MHD limit, the a → 2/3 and a → ∞ limits, and the no-magnetic-field limit are recovered correctly. If I were testing a viscous MHD code, I would use this curve as a benchmark.\n\nThe soft spot is the second perturbative solution, Eq. (50). It does not satisfy the initial condition T(τ0)=1. For the first-order equation (48), requiring T1(τ0)=0 fixes the homogeneous coefficient to -3/(4-6a); Eq. (49) instead uses (1-6a)/(4-6a), which gives T1(τ0)=1, hence T_NS-2(τ0)=1+ε1. The 5% deviation seen at ε1=0.05 in Fig. 5 is therefore an initialization artifact, not a truncation error. The derivation also mixes orders by replacing T0^3(1+3ε1 T1/T0) with T0^3 before linearizing. The paper itself flags the approximation as \"very significant,\" but it still presents Eq. (50) as a valid solution, which understates the problem.\n\nThe stability claim in Sec. IV A is supported only by numerical evidence, not a proof; the wording \"we will prove this\" overstates it. The Israel-Stewart part is purely numerical and adds less.\n\nNet: this is a useful contribution but not ready as is. The first solution is solid and novel, the second should either be corrected (choose the right homogeneous term) or withdrawn. I would accept this for peer review, with the expectation of a major revision that fixes Eq. (50), softens the stability language, and ideally provides the code or data for the numerics. If those changes land, it becomes a good benchmark paper. As it stands, I would not quote Eq. (50).","headline":"First perturbative solution (NS-1) is a solid new benchmark; the second (NS-2) violates its initial condition and needs fixing before use.","tokens_in":19763,"tokens_out":4304,"would_cite":true,"duration_ms":38415,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives analytic temperature solutions for a 1+1 dimensional, longitudinally boost-invariant viscous magnetohydrodynamic flow and shows that a decaying magnetic field and shear viscosity both heat the quark-gluon plasma…","keywords":["relativistic magnetohydrodynamics","Bjorken flow","shear viscosity","boost invariance","quark-gluon plasma","heavy-ion collisions","analytical solutions"],"falsifier":"Numerically integrate Eq. (53) with $a=2$, $\\tau_0=0.6$ fm/c, $T_0=0.65$ GeV, $\\epsilon_2=1$, and $\\epsilon_1 = 0.02$, $0.05$, $0.1$, and $0.5$, then compare the output with Eq. (50) and Eq. (40). The paper itself reports a five-percent gap for Eq. (50) at $\\epsilon_1=0.05$; observing that the gap grows faster than linearly with $\\epsilon_1$, or that Eq. (40) departs from numerics at small $\\epsilon_1$, would show how much of the claimed early heating is an artifact of the perturbative truncation.","tokens_in":2077,"feed_emoji":"🧲","tokens_out":4203,"duration_ms":90729,"temperature":0.7,"pith_summary":"This paper asks how a quark-gluon fluid cools when it is stretched along one spatial direction, carries a decaying transverse magnetic field, and resists shearing. It derives analytic formulas for the fluid temperature as a function of proper time in 1+1 dimensional, longitudinally boost-invariant relativistic magnetohydrodynamics with a conformal equation of state. The two new perturbative solutions for the first-order Navier-Stokes limit show that both the magnetic field and shear viscosity deposit energy into the fluid, producing an early temperature peak. In the second-order Israel-Stewart theory the paper's numerical solutions show the same heating effect and a slower cooling when the initial shear stress is larger. These formulas give a simple benchmark for how magnetic-field and viscous effects modify the standard Bjorken cooling law in heavy-ion collisions.","feed_headline":"Magnetic fields and viscosity slow quark-gluon cooling","feed_subtitle":"Analytic temperature formulas for boost-invariant relativistic MHD match numerics when the magnetic field is small.","key_machinery":"The load-bearing object is the ordinary differential equation for the normalized temperature $\\tilde T$, Eq. (28), which combines ideal Bjorken cooling $-\\tilde T/(3\\tau)$, magnetic heating $\\epsilon_1 \\tilde T^3 (\\tau_0/\\tau)^{2a}/\\tau$, and viscous heating $\\epsilon_2/\\tau^2$. The first perturbative solution follows the nonconserved-charge method of rewriting the equation as $d\\tilde T/d\\tau + \\tilde T/(3\\tau) = \\tilde T\\, d\\lambda/d\\tau$, then solving the auxiliary function $x(\\tau)=\\exp[\\lambda(\\tau)-\\lambda(\\tau_0)]$ order by order in $\\epsilon_1$, which introduces hypergeometric functions. The second solution instead writes $\\tilde T = \\tilde T_0 + \\epsilon_1 \\tilde T_1$ and replaces $\\tilde T_0^3(1+3\\epsilon_1 \\tilde T_1/\\tilde T_0)$ with $\\tilde T_0^3 \\approx \\tau_0/\\tau$; that replacement is the step that limits its range of validity. The Israel-Stewart treatment adds a relaxation equation for the normalized shear stress $\\tilde\\pi$, making the system two coupled differential equations that are solved numerically.","core_discovery":"The central claim is that the energy-conservation equation for a boost-invariant, non-resistive viscous fluid in a transverse power-law magnetic field reduces to the single first-order ordinary differential equation $$\\partial_\\tau \\tilde T + \\frac{\\tilde T}{3\\tau} - \\epsilon_1 \\frac{\\tilde $T^{3}$}{\\tau}\\left(\\frac{\\tau_0}{\\tau}\\right)^{2a} - \\frac{\\epsilon_2}{\\$tau^{2}$} = 0,$$ where $\\tilde T = T/T_0$, $\\epsilon_1 = (a-1)\\sigma/(12a_1)$ packages the initial field strength and the power-law decay exponent $a$, and $\\epsilon_2$ packages shear viscosity. The paper constructs two perturbative analytic solutions of this equation. The first, Eq. (40), is obtained through the nonconserved-charge method and tracks the numerical solution closely for small $\\epsilon_1$ even when the shear viscosity is not tiny; the second, Eq. (50), is a linear-in-$\\epsilon_1$ expansion that matches numerics only for very small magnetic fields and runs about five percent high at $\\epsilon_1 = 0.05$. Both solutions interpolate among the Bjorken, Victor-Bjorken, and Azwinndini-Bjorken limits, and the Israel-Stewart numerical solutions show that a stronger initial shear stress slows the temperature drop.","pith_inferences":["The early temperature peak occurs for proper times below the initial time $\\tau_0 = 0.6$ fm/c, outside the plotted window, so the formulas make a concrete prediction about the pre-thermal regime that a full 3+1 dimensional MHD code could check.","A direct testable extension would replace the power-law magnetic decay with a decay driven by temperature-dependent electrical conductivity; the same reduction to a single ODE would show whether the heating peak survives a more realistic field evolution.","Because the second solution loses accuracy for $\\epsilon_1 > 0.05$, a resummed or matched-asymptotic version of Eq. (50) could extend the analytic description toward the paper's stated stability bound $B_0^2 < 6 a_1 T_0^4$.","The same perturbative machinery could be applied to anomalous magnetohydrodynamics with a chiral magnetic current, where the field and the chiral transport coefficient couple in an additional term."],"forward_implications":["For small magnetic fields with decay exponent $a > 1$, the fluid temperature declines more slowly than ideal Bjorken cooling and develops an early peak whose height grows with both $\\epsilon_1$ and $\\epsilon_2$.","The two perturbative solutions reduce to the Bjorken, Victor-Bjorken, and Azwinndini-Bjorken limits, so they serve as interpolation formulas connecting ideal MHD, purely viscous flow, and combined viscous magnetohydrodynamics.","The first Navier-Stokes solution, Eq. (40), remains close to the numerical solution for $\\epsilon_1$ up to at least $0.1$ when $\\epsilon_2 = 1$, giving a practical analytic benchmark for code comparisons in the boost-invariant setting.","In the Israel-Stewart theory, increasing the initial shear stress from $0.1 p_0$ to $2 p_0$ visibly lowers the cooling rate, meaning the initial viscous stress is a controlling input for the temperature history.","Adding the magnetic field to a viscous flow produces slower cooling than viscosity alone, so neglecting the field in dissipative flow models would misestimate the temperature evolution."],"supporting_citations":[{"why":"Provides the ideal MHD Victor-Bjorken flow solution that the paper extends to include shear viscosity.","marker":"[49]"},{"why":"Introduces the Bjorken-flow MHD setup with magnetization and the temperature and magnetic-field parameterization used here.","marker":"[50]"},{"why":"Supplies the nonconserved-charge method that produces the first perturbative solution, Eq. (40).","marker":"[51]"},{"why":"Gives the non-resistive viscous MHD energy-momentum tensor and Israel-Stewart equations in a magnetic field that the paper adopts.","marker":"[57]"},{"why":"Defines the second-order dissipative fluid dynamics and the shear-stress relaxation structure used in the Israel-Stewart section.","marker":"[62]"},{"why":"Provides the Azwinndini-Bjorken viscous flow solution, thermodynamic relations, and transport coefficients that the new solutions reduce to in the viscous-only limit.","marker":"[63]"}],"fun_headline_variants":["Viscous MHD: B-field and shear stress slow QGP cooling","Magnetic fields and viscosity delay quark-gluon plasma cooling","New MHD solutions: B-field and viscosity control QGP temperature","Relativistic MHD: viscosity and B-field slow QGP cooling","Boost-invariant MHD: B-field and viscosity slow QGP temperature drop"],"cache_read_input_tokens":21760,"weakest_assumption_plain":"The calculations go through only if the magnetic-field parameter $\\epsilon_1$ is small enough that the perturbative series can be truncated at first order; for the second analytic solution the additional replacement $\\tilde T_0^3(1+3\\epsilon_1 \\tilde T_1/\\tilde T_0) \\approx \\tilde T_0^3 \\approx \\tau_0/\\tau$ is an even cruder assumption, and the paper itself labels its limitations as very significant.","fun_headline_variants_meta":{"raw":{"variants":["Viscous MHD: B-field and shear stress slow QGP cooling","Magnetic fields and viscosity delay quark-gluon plasma cooling","New MHD solutions: B-field and viscosity control QGP temperature","Relativistic MHD: viscosity and B-field slow QGP cooling","Boost-invariant MHD: B-field and viscosity slow QGP temperature drop"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000631,"raw_usage":{"total_tokens":2940,"prompt_tokens":999,"completion_tokens":1941,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":1846}},"tokens_in":615,"tokens_out":1941,"duration_ms":16708,"temperature":1.0,"reasoning_tokens":1846,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:33:54.471327+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate Eq. (53) with $a=2$, $\\tau_0=0.6$ fm/c, $T_0=0.65$ GeV, $\\epsilon_2=1$, and $\\epsilon_1 = 0.02$, $0.05$, $0.1$, and $0.5$, then compare the output with Eq. (50) and Eq. (40). The paper itself reports a five-percent gap for Eq. (50) at $\\epsilon_1=0.05$; observing that the gap grows faster than linearly with $\\epsilon_1$, or that Eq. (40) departs from numerics at small $\\epsilon_1$, would show how much of the claimed early heating is an artifact of the perturbative truncation.","supporting_citations":[{"cited_title":"External-magnetic-field-induced paramagnetic squeez- ing effect in heavy-ion collisions at energies available at the CERN Large Hadron Collider","cited_arxiv_id":null,"evidence_quote":"Provides the ideal MHD Victor-Bjorken flow solution that the paper extends to include shear viscosity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Bjorken-flow MHD setup with magnetization and the temperature and magnetic-field parameterization used here."},{"cited_title":"An- alytic Bjorken flow in one-dimensional relativistic magnetohy- drodynamics","cited_arxiv_id":null,"evidence_quote":"Supplies the nonconserved-charge method that produces the first perturbative solution, Eq. (40)."},{"cited_title":"Haddadi Moghaddam, W","cited_arxiv_id":null,"evidence_quote":"Gives the non-resistive viscous MHD energy-momentum tensor and Israel-Stewart equations in a magnetic field that the paper adopts."},{"cited_title":"Second order dissipative fluid dynam- ics for ultrarelativistic nuclear collisions","cited_arxiv_id":null,"evidence_quote":"Defines the second-order dissipative fluid dynamics and the shear-stress relaxation structure used in the Israel-Stewart section."},{"cited_title":"Most and Jorge Noronha","cited_arxiv_id":null,"evidence_quote":"Provides the Azwinndini-Bjorken viscous flow solution, thermodynamic relations, and transport coefficients that the new solutions reduce to in the viscous-only limit."}],"review_version":1}