{"id":"693c8f56-cb2a-4aa1-b38f-b8d66df1e3de","arxiv_id":"2502.04611","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using 3+1D resistive magnetohydrodynamics, the authors show that the slope of proton-antiproton charge-dependent directed flow in Au+Au collisions at 200 GeV varies with the QGP's electrical conductivity and can change sign when an initial positive charge density is included.","lead":"A numerical model of quark-gluon plasma (QGP) with finite electrical resistivity shows that the difference in flow between protons and antiprotons in gold-gold collisions depends on the plasma's electrical conductivity and on collision centrality. The result suggests that measurements of charge-dependent directed flow could be used to learn about electrical transport in the QGP, provided initial field and charge conditions are better constrained.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The σ-scan in Fig. 1 compares runs with identical initial EM fields, but Eq. A1 makes those fields σ-dependent; the claimed conductivity sensitivity may be a normalization artifact.","rationale":"The reader's weakest_assumption focuses on the known overestimation of the initial EM fields (App. C), which affects absolute magnitudes. I agree that this is a limitation, but the more load-bearing problem for the central claim is the internal consistency of the σ scan. The central claim is a comparative statement: Δv1 depends on σ across centralities. For that comparison to be meaningful, σ must be the only changed parameter. Eq. (A1) gives initial fields that depend strongly on σ; the text of Sec. V asserts the initial conditions are identical. These two statements cannot both be true unless a fixed reference σ was used to set initial conditions for all runs, in which case large-σ runs are initialized with unrealistically strong fields. This directly contaminates the slopes in Figs. 5 and 6 and the conclusion that larger σ gives larger |slope|. The paper is otherwise careful: it repeatedly labels the work as qualitative, provides a useful appendix on field-strength sensitivity (App. C), and reports the centrality trend correctly. A request for the authors to specify and justify the initialization convention, and to re-run the σ scan with σ-consistent initial fields, is the appropriate level of conditionality. Hence the reader's CONDITIONAL verdict stands, with the focus shifted to this confound.","tokens_in":15105,"tokens_out":8642,"duration_ms":82893,"concrete_test":"Rerun or instrument the setup for Fig. 1: extract the numerical values of the initial By(x=0, y=0, η=0, τ0=0.4 fm/c) for each σ run. (i) If these values are identical across σ, identify the common σ used to evaluate Eq. (A1), then recompute Fig. 6 with initial fields evaluated at each run's own σ; if the slopes change sign or ordering, the reported σ-sensitivity is an artifact. (ii) If the values differ across σ, verify whether the 'earliest point is equal' statement holds; if it does not, correct Fig. 1 and re-quantify the centrality dependence. In either case, report the initial-field values in the paper or supplement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Δv1 is sensitive to σ rests on a controlled comparison of runs with different σ. That comparison is undermined by an unspecified initialization of the EM fields. Sec. V and Fig. 1 state 'we have assumed the initial conditions are the same, the point at the earliest time in Fig. 1 is equal.' But the analytic initial conditions of Eq. (A1) depend on σ: the field scales with σ and decays as exp(−b²σ/(4τ0)). At b=10 fm and τ0=0.4 fm/c, increasing σ from 0.0294 to 0.294 fm⁻¹ reduces the initial By by roughly seven orders of magnitude. If the code instead uses a common reference σ for the initial fields of all runs, then the large-σ runs start with a field far above the physical in-medium value for that conductivity, inflating both the magnetic-field lifetime (Fig. 1) and the Δv1 slopes (Fig. 6). If the code uses σ-dependent initial fields, then the statement that the earliest points are equal is false. Sec. IV and Appendix A do not specify which choice was made. Because the abstract claims sensitivity to resistivity across centralities, and Figs. 5 and 6 are the evidence, this ambiguity means the central quantitative result is not cleanly isolated from the initialization convention.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a 3+1D relativistic resistive magnetohydrodynamic (RRMHD) model for Au+Au collisions at sqrt(s_NN)=200 GeV. It initializes the QGP with an optical Glauber profile and initial electromagnetic fields from Tuchin's analytic in-medium solution, evolves the system with a constant scalar conductivity, and computes directed flow via Cooper-Frye freezeout with a finite electric chemical potential on the hypersurface. The central results are: (i) the time evolution of the magnetic field is sensitive to sigma, with larger sigma giving longer field lifetimes; (ii) the slope d(delta v1)/dy for protons and antiprotons is negative and grows in magnitude with sigma; (iii) the slope decreases toward central collisions; and (iv) adding an initial positive electric charge density (Zq=0.4) can flip the slope positive, as suggested in Ref. [15]. The authors explicitly list simplifying assumptions: constant scalar sigma, ideal equation of state, overestimated initial fields, no spectator currents during the evolution, and no afterburner.","tokens_in":15366,"tokens_out":9567,"duration_ms":97125,"significance":"The paper addresses an active experimental question and provides a concrete mechanism connecting the electric conductivity of QGP to a charge-dependent flow observable. Its strengths are the honest treatment of limitations, the consistency of the negative-slope sign with Ref. [15], and the explicit scaling test in Appendix C. If the conductivity-sensitivity claim were quantitatively robust, the model would yield falsifiable centrality trends that can be compared with STAR and ALICE data. However, because the slopes are roughly proportional to the initial field strength and the sigma-scan uses an ambiguous initialization convention, the present results establish a proof of mechanism rather than a quantitative extraction. The sign-flip demonstration with an initial charge density is illustrative, not a quantitative prediction.","major_comments":[{"comment":"The controlled sigma-scan is undermined by an unspecified initialization convention for the electromagnetic fields. Equation (A1) makes the initial fields explicitly sigma-dependent through both a prefactor sigma and the exponential exp[-b^2 sigma/(4(t +/- vz))]. At a fixed initial time tau0=0.4 fm/c, the initial By values for different sigma are not equal; for a field point at transverse distance b=10 fm from a single charge, the exponential factor alone changes by about seven orders of magnitude when sigma is increased from 0.0294 to 0.294 fm^-1. The statement in Sec. V that 'we have assumed the initial conditions are the same, the point at the earliest time in Fig. 1 is equal' is therefore either inconsistent with Eq. (A1) or implies that a common reference conductivity was used to generate the initial fields for all runs. The manuscript never states which convention was adopted. If a common reference is used, the high-sigma runs start with a field that is unphysically large for their conductivity, inflating the field lifetime in Fig. 1 and the slopes in Fig. 6; if sigma-dependent initial fields are used, the earliest points in Fig. 1 cannot be equal. This ambiguity directly affects the abstract's claim of sensitivity across centralities, so it must be resolved and the sigma-scan rerun with a stated, physically motivated convention.","section":"Sec. V, Fig. 1, Appendix A"},{"comment":"The quantitative slopes are not yet isolated from the known overestimation of the initial fields. Appendix C demonstrates that reducing the initial field strength to 10% reduces the delta v1 slope by an order of magnitude. Therefore, the statement in Sec. VI that 'an order increase in the electric conductivity sigma brought a similar increase in the slope' does not by itself establish a sigma-specific effect, because the same order-of-magnitude change can be generated by renormalizing the initial fields. The authors acknowledge this degeneracy in Appendix C, but the main-text presentation in Fig. 6 and the abstract should be conditioned on it. A two-dimensional scan in (sigma, initial field strength), or at least a normalization of the initial fields to a common physical reference, is needed before the centrality dependence in Fig. 6 can be attributed to sigma.","section":"Sec. VI and Appendix C"}],"minor_comments":[{"comment":"Both equations are written as □M(tau B2), but the surrounding text and Eq. (18) indicate that the third component B3 is intended.","section":"Eqs. (11) and (15)"},{"comment":"The caption spells the model as 'RRHMD'; it should be 'RRMHD'.","section":"Fig. 6 caption"},{"comment":"The functions f±(eta), H(eta), and the normalization used in Eq. (B1) are not defined in the text; please provide explicit forms or precise references to Ref. [31].","section":"Appendix B"},{"comment":"The sentence 'an order increase in the electric conductivity sigma brought a similar increase in the slope' should specify the sigma range and centrality to which it applies, since Fig. 6 does not show a single power-law relation across all centralities.","section":"Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and the model is a useful step, but the initialization ambiguity is serious enough that the sigma-scan should be redone. I would not reject because the mechanism and numerical setup are plausible and the ambiguity can be fixed by clarifying the convention and rerunning the comparison. The paper's claim of being 'the first work to establish such behavior' is somewhat overstated given earlier resistive-model studies, but this can be tempered with appropriate citations and wording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's new pieces are genuine: proton/antiproton Δv1 in symmetric Au+Au, a centrality scan, and a first pass at an initial positive charge density. The model setup and freezeout are described carefully, and the authors are open about their simplifications. But the central claim—that Δv1 is sensitive to σ—rests on a comparison that is not cleanly defined. The initial EM fields come from the in-medium analytic solution of Eq. A1, which depends on σ, while the text says the earliest points in Fig. 1 are equal because the initial conditions are the same. Those two statements are in tension. At b=10 fm and τ0=0.4 fm/c, the analytic By changes by roughly seven orders of magnitude when σ goes from 0.0294 to 0.294 fm⁻¹. If the code uses a common initial field for all runs, the high-σ runs start with a field far above the physical value for that conductivity, which inflates the field lifetime and the resulting slope. If it uses σ-dependent fields, the earliest points would not be equal. The paper does not say which. Since Appendix C shows the slope scales roughly with the initial field strength, this ambiguity means Fig. 6 does not isolate the dynamical role of σ. That is a load-bearing problem, not a footnote.\n\nThe rest of the paper is more solid. The µQ-on-freezeout mechanism is laid out clearly, and the negative-slope result is consistent with Ref. [15]. The five-significant-figure slopes in Eqs. 25 and 26 are overpressed given the model's known sensitivity, and the 'first work' phrasing in Sec. VII is a bit strong next to Refs. [15,21]. Both are minor next to the initialization issue.\n\nThis is worth refereeing. The model extension to protons and the centrality pattern are potentially useful for the EM-field-in-heavy-ion community. But the referee should ask the authors to state explicitly what initial fields they used for each σ and, if they used a common field, to justify it physically or rerun with σ-consistent initial fields. Without that, the σ-sensitivity claim is not established. I'd bring it to a reading group as a case study in how a simulation comparison can be quietly contaminated by initial conditions.","headline":"A useful RRMHD extension to proton/antiproton Δv1 with a centrality scan, but the σ-scan anchoring the main claim is contaminated by an inconsistent EM-field initialization.","tokens_in":15923,"tokens_out":4600,"would_cite":false,"duration_ms":44766,"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 claims that the proton-antiproton directed-flow difference in symmetric heavy-ion collisions is a conductivity-sensitive observable, and that an initial positive charge can reverse its slope.","keywords":["relativistic resistive magnetohydrodynamics","charge-dependent directed flow","electric conductivity","quark-gluon plasma","heavy-ion collisions","electromagnetic fields","Au+Au collisions","freezeout chemical potential"],"falsifier":"A calculation of pre-equilibrium EM fields with a first-principles kinetic model (e.g., CGC/glasma) that gives fields an order of magnitude smaller than the analytic solution would also reduce the predicted $\\Delta v_1$ slopes by roughly an order of magnitude. If that reduced slope falls below the sensitivity of STAR or ALICE data while conductivity values stay at lattice estimates, the paper's quantitative predictions fail; the qualitative conductivity dependence would still be testable in the centrality trend.","tokens_in":14873,"feed_emoji":"⚛️","tokens_out":3865,"duration_ms":34062,"temperature":0.7,"pith_summary":"This paper argues that the difference in directed flow between protons and antiprotons in symmetric gold-gold collisions at RHIC energies carries information about the electric conductivity of the quark-gluon plasma. Using a 3+1D relativistic resistive magnetohydrodynamic model, it shows that the slope of this charge-dependent flow difference changes systematically with conductivity across collision centralities. The model connects the conductivity to the lifetime of the magnetic field, which drives Faraday and Lorentz currents that leave an electric chemical potential imprinted on the freezeout surface. A positive initial electric charge in the plasma can flip the slope positive, offering an alternative explanation for the sign seen in data.","feed_headline":"Conductivity shows up in proton flow split","feed_subtitle":"Proton vs antiproton flow slope at RHIC energy tracks the quark-gluon plasma's electric conductivity.","key_machinery":"The central object is the relativistic resistive magnetohydrodynamic (RRMHD) evolution coupled to Ohm's law $J^\\mu = q u^\\mu + \\sigma F^{\\mu\\nu}u_\\nu$. The conductivity $\\sigma$ controls how the magnetic field decays like a non-linearly damped oscillator, and the resulting Lorentz and Faraday currents accumulate an electric chemical potential $\\mu_Q$ on the freezeout hypersurface. The Cooper-Frye formula then converts $\\mu_Q$ into a charge-dependent split of proton and antiproton directed flow, whose slope is the observable $\\Delta v_1$.","core_discovery":"The paper's central claim is that the slope of $\\Delta v_1 = v_1^+ - v_1^-$ for protons and antiprotons in symmetric Au+Au collisions at $\\sqrt{s_{NN}}=200$ GeV is sensitive to the electric conductivity $\\sigma$ of the quark-gluon plasma. In the RRMHD model, a finite conductivity controls the decay of the magnetic field; larger $\\sigma$ prolongs the field, strengthens the Lorentz and Faraday currents, and produces a larger negative slope of $\\Delta v_1$ across centralities. The paper also shows that if the plasma initially carries a net positive charge, the slope can become positive, so the sign of the measured slope need not come solely from transported baryon charge.","pith_inferences":["Because the sign of the slope can be flipped by the initial charge density, measurements of $\\Delta v_1$ slope alone cannot uniquely fix $\\sigma$ without independent control of the initial charge distribution; joint fits to centrality trends might separate the two.","The same mechanism should generate charge-dependent flow for other identified hadrons (kaons, pions) with magnitudes set by their electric charge, which could be tested by comparing STAR and ALICE data across energies.","If the conductivity is shown to be tensor-valued in strong magnetic fields, the predicted $\\Delta v_1$ slope could develop a dependence on the orientation of the field relative to the reaction plane, an effect this scalar-conductivity model cannot produce."],"forward_implications":["The slope of proton-antiproton $\\Delta v_1$ is a probe of QGP electric conductivity even in symmetric collisions, not only asymmetric ones.","Larger conductivity implies longer-lived magnetic fields and larger negative slopes of $\\Delta v_1$, with Faraday current dominating Lorentz current.","An initial net positive electric charge in the plasma can reverse the slope to positive, providing an alternative to transported-baryon explanations of STAR data.","Conductivity dependence varies non-trivially with centrality, with lattice-like $\\sigma$ showing the least change across centralities.","The result is qualitative; matching data quantitatively requires better initial conditions, a $\\mu_Q$-dependent equation of state, and an afterburner."],"supporting_citations":[{"why":"Supplies the analytic in-medium Maxwell solutions used to initialize the electromagnetic fields, overestimation of which is flagged in the paper.","marker":"[25]"},{"why":"Establishes the Faraday-versus-Lorentz current competition and the expected negative slope of $\\Delta v_1$ that this model reproduces.","marker":"[15]"},{"why":"Previous RRMHD study of pions that this work extends to protons, antiprotons, and multiple centralities.","marker":"[21]"},{"why":"Describes the RRMHD code that solves the coupled hydrodynamic and Maxwell equations used here.","marker":"[22]"},{"why":"Earlier application of the RRMHD model to directed flow in symmetric and asymmetric collisions, providing the baseline setup.","marker":"[23]"},{"why":"STAR data on charge-dependent directed flow that motivates the proton-antiproton observable and provides the experimental comparison.","marker":"[16]"},{"why":"Lattice QCD estimates of the electric conductivity that justify the scalar conductivity values and the restriction $\\sigma<5$ fm$^{-1}$.","marker":"[7]"},{"why":"Phenomenological ansatz for initial baryon charge density that the paper adapts to estimate the initial electric charge density.","marker":"[40]"}],"fun_headline_variants":["Charge flow split reveals quark-gluon conductivity","Proton-antiproton flow bent by plasma conductivity","How conductivity steers proton vs antiproton flow","Plasma conductivity leaves fingerprint in flow split","Directed flow slope tracks QGP electric conductivity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The initial electromagnetic fields are taken from an analytic in-medium solution that the paper itself notes overestimates the field strength, and reducing that input to 10% shrinks the predicted flow slope by an order of magnitude, so the quantitative size (and the sign-flip demonstration) rests on a known-overestimated input.","fun_headline_variants_meta":{"raw":{"variants":["Charge flow split reveals quark-gluon conductivity","Proton-antiproton flow bent by plasma conductivity","How conductivity steers proton vs antiproton flow","Plasma conductivity leaves fingerprint in flow split","Directed flow slope tracks QGP electric conductivity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1413,"prompt_tokens":831,"completion_tokens":582,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":511}},"tokens_in":447,"tokens_out":582,"duration_ms":6316,"temperature":1.0,"reasoning_tokens":511,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T22:07:57.463266+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation of pre-equilibrium EM fields with a first-principles kinetic model (e.g., CGC/glasma) that gives fields an order of magnitude smaller than the analytic solution would also reduce the predicted $\\Delta v_1$ slopes by roughly an order of magnitude. If that reduced slope falls below the sensitivity of STAR or ALICE data while conductivity values stay at lattice estimates, the paper's quantitative predictions fail; the qualitative conductivity dependence would still be testable in the centrality trend.","supporting_citations":[{"cited_title":"Charge-dependent anisotropic flow in high-energy heavy-ion collisions from relativistic resistive magneto-hydrodynamic expansion","cited_arxiv_id":"2212.02124","evidence_quote":"Earlier application of the RRMHD model to directed flow in symmetric and asymmetric collisions, providing the baseline setup."},{"cited_title":"Conductivities of magnetic quark-gluon plasma at strong coupling","cited_arxiv_id":"1809.02178","evidence_quote":"Lattice QCD estimates of the electric conductivity that justify the scalar conductivity values and the restriction $\\sigma<5$ fm$^{-1}$."},{"cited_title":"Resistive relativistic MHD simulations of astrophysical jets","cited_arxiv_id":"2308.09477","evidence_quote":"Phenomenological ansatz for initial baryon charge density that the paper adapts to estimate the initial electric charge density."}],"review_version":1}