{"id":"a3d587e2-eec6-4a68-8f32-8687efeff8d8","arxiv_id":"2506.23006","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A variational hydrodynamic model with a step-function pair distribution gives accurate longitudinal dispersion relations for the Yukawa one-component plasma across wide coupling and screening ranges.","lead":"This paper extends a variational hydrodynamic framework, originally developed for the one-component plasma, to the screened Yukawa one-component plasma. It derives dispersion relations that match molecular dynamics simulations well at short wavelengths, even where ordinary hydrodynamics fails.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing risk is the equilibrium-g closure in Eq. (3), not the step-function shape; the strong-screening sound-speed error and the paper's own stated relaxation caveat point there, but the claimed MD check is undocumented.","rationale":"The reader's weakest assumption was the step-function approximation for g. The paper, however, reports its own MD-based check indicating that the step-function and equation-of-state inputs affect the strong-screening sound speed by only about 1%. If that check is sound, the step-function form is not the main liability. The more fundamental and self-identified weak point is the closure that g is the equilibrium pair distribution function evaluated at local thermodynamic variables. This closure enters the derivation at Eq. (3) and propagates through the longitudinal dispersion relation (16) via Eqs. (18)-(20); it is exactly the kind of assumption that can fail at strong screening, where relaxation and nonlocal memory are expected to matter. The paper acknowledges this possibility in the conclusion but provides no quantitative test. Because the central claim is specifically that the variational dispersion reproduces MD across a wide parameter range, an untested closure at the heart of the variational principle is the most load-bearing concern. The existing conditional verdict already captures this uncertainty, so I do not propose changing the verdict; rather, I would sharpen the condition: the finite-wavelength agreement and the strong-screening discrepancy should be re-examined under an independently tested out-of-equilibrium g closure. The proposed MD check would settle whether the closure or some other ingredient is responsible for the 17% sound-speed error and would also clarify the moderate-coupling overestimation visible in Figs. 3-4.","tokens_in":21151,"tokens_out":6556,"duration_ms":77386,"concrete_test":"Perform nonequilibrium MD at the worst case κ0=3, Γ0=40 (and one moderate finite-q case, e.g., κ0=2, Γ0=80): launch a long-wavelength longitudinal density wave, phase-average the instantaneous pair distribution function g(r,t) conditioned on local n and T, and compare it to the equilibrium g(n,T,r) used in Eq. (3). Quantify the weighted difference in the integrals (18)-(20). If the difference is small (< a few percent), the closure is supported and the 17% sound-speed error must be traced elsewhere; if large, the equilibrium-g ansatz, not the step-function form, is the load-bearing failure. A minimal first step is to document the authors' MD check: evaluate Eq. (16) with MD-computed static g and compare to the step-function result at the same state points.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The variational Lagrangian (3) closes the nonlocal interaction term by taking g to be the equilibrium pair distribution function evaluated at the local reference density and temperature, g(nref, T(n,s), |a-a'|). All q-dependent corrections in the longitudinal dispersion (16), specifically the integrals (18)-(20) and their Γ-derivatives, inherit this equilibrium closure. This is the model's central physical assumption, and the paper itself identifies it as the likely source of the strong-screening sound-speed failure: the conclusion states that 'our model assumes a particular form for the out-of-equilibrium behavior of the pair distribution function' and that 'at higher screening, relaxation behaves differently.' At κ0=3, Γ0=40 the sound speed is overestimated by 17%; the authors state that their MD-based check removes the step-function and equation-of-state approximations as the cause (effect ~1%), but no details or data for this check are given. If that check is correct, the residual error must come from the equilibrium-g closure, not from the step-function shape. The finite-wavelength comparisons (Figs. 3-4) do not isolate this, and at moderate coupling (e.g., κ0=1, Γ0=40) the model visibly overestimates the MD peak positions while QLCA agrees better. Until the closure is tested independently, the central claim that the variational Lagrangian quantitatively reproduces MD across the full parameter range is not settled.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends a previously developed variational Lagrangian formulation for the classical one-component plasma to the Yukawa one-component plasma. The authors derive equations of motion in reference and laboratory coordinates, momentum and energy conservation laws, and longitudinal and transverse dispersion relations. The pair distribution function is approximated by a step function whose cutoff radius is fixed by matching the internal energy to a simulation-based fit. The longitudinal sound speed and finite-wavelength dispersion are compared with molecular dynamics data and with QLCA, extended QLCA, and Euler-plus-mean-field theory. The model performs well at weak and moderate screening and at finite wavelengths, while at strong screening (kappa = 2, 3) the sound speed is overestimated by up to 17%.","tokens_in":21406,"tokens_out":9247,"duration_ms":95579,"significance":"If the claims hold, this is a useful extension of a variational hydrodynamic framework to a model relevant for dusty plasmas, ultracold neutral plasmas, and related systems. The paper provides a fairly complete derivation, checks consistency with thermodynamics, includes the authors' own molecular dynamics spectra, and compares with several established theories. A notable strength is that, once the step-function cutoff is fixed by the internal energy fit, the dispersion predictions contain no adjustable parameters for the comparisons shown. The main weakness is that the central closure assumption, namely that the out-of-equilibrium pair distribution function equals the equilibrium pair distribution evaluated at local density and temperature, is not tested independently, and the claimed MD check that isolates the source of the strong-screening sound-speed error is undocumented.","major_comments":[{"comment":"The statement that 'our own MD simulations, which remove both approximations, together with the dispersion law given by Eq. (16), indicate that their effect is around 1%' is load-bearing: it is used to rule out the equation of state and the step-function approximation as causes of the 17% sound-speed error at κ0=3, Γ0=40. No details of these simulations or of the calculation are provided. Please specify how the exact pair distribution function and its derivatives with respect to temperature were obtained from MD and inserted into Eq. (16), how the pressure/equation-of-state input was replaced, what the resulting sound speeds are for the affected states, and what the statistical uncertainties are. A table or figure presenting these results is needed; without it, the attribution of the discrepancy to the closure in Eq. (3) is not verifiable.","section":"§VI.A and Conclusions"},{"comment":"The equilibrium-g closure is the central physical assumption of the model: the nonlocal term in Eq. (3) uses the equilibrium pair distribution function evaluated at the local reference density and temperature, and all q-dependent terms in the longitudinal dispersion, specifically the integrals j, ℓ, and b in Eqs. (18)-(20) and their Γ-derivatives, inherit this closure. The finite-wavelength comparisons in Figs. 3 and 4 do not isolate this assumption because the same closure is used in every model curve; moreover, at κ0=1, Γ0=40 and κ0=2, Γ0=80 the variational model visibly overestimates the MD peak positions. The authors themselves acknowledge in the Conclusions that 'our model assumes a particular form for the out-of-equilibrium behavior' and that 'at higher screening, relaxation behaves differently.' To support the claim of quantitative accuracy across the full parameter range, the closure should be tested independently, for example by computing Eq. (16) with the exact equilibrium g from MD at the problematic states and comparing the resulting sound speed and dispersion with both the step-function version and the MD data.","section":"§II, Eq. (3), and §IV, Eqs. (16)-(20)"},{"comment":"The paper claims that, 'considering the full coupling range, our approach provides the best overall agreement with the MD data compared to the other theoretical models tested.' This claim is based on visual inspection of Figs. 3 and 4, where the finite-wavelength comparison is restricted to |q0| ≤ 3 and dissipative effects are excluded because peaks become broad at larger wave vectors. Since the strong-screening sound-speed error is sizable and the closure is untested, a quantitative error metric (e.g., root-mean-square deviation between predicted and MD peak positions over the displayed range) would strengthen the claim. As written, the conclusion that the model is the most consistent across the full range is plausible but not quantitatively established.","section":"§VI.B and Conclusions"}],"minor_comments":[{"comment":"The text twice writes 'EQCLA' where 'EQLCA' is meant; please correct the typographical error.","section":"§VI.B"},{"comment":"Please define the equilibrium Wigner-Seitz radius a0 explicitly and clarify the notation c_L^2/(ω_p a0/κ0)^2; the current expression is not immediately transparent dimensionally.","section":"§VI.A, Eq. (28)"},{"comment":"The caption states that error bars reflect an uncertainty of 5% for panels (a) and (b) and 10% for panels (c) and (d), but no source or justification for these values is given. If these are taken from Ref. [37], cite the relevant values; if they are chosen by the authors, state the criterion.","section":"Fig. 2 caption"},{"comment":"The sentence 'A possible explanation for the disagreement of all theories with the simulations at κ0=0.5 and large Γ0 could be related to difficulties in extracting the sound speed from finite wave number MD data' is speculative. Either support it with a quantitative estimate, for example by recomputing the sound speed from the MD dispersion at different smallest wave vectors, or remove it.","section":"§VI.A"},{"comment":"The derivation of the adiabatic derivative leading to Eq. (15) is stated in one sentence. A brief derivation or a reference to the corresponding OCP calculation in Ref. [22] would help the reader verify the definition of f(Γ,κ).","section":"§V, Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a natural extension of the authors' earlier OCP work and the novelty is clear. The main concern for the editor is that the undocumented MD check in §VI.A is used to attribute the strong-screening discrepancy to the nonequilibrium closure; if that check cannot be provided or does not reproduce the claimed 1% effect, the conclusions about the source of error would need substantial revision. The paper otherwise appears sound and within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a legitimate extension of the variational hydrodynamics framework from the OCP to the Yukawa OCP, and the finite-wavelength dispersion results are the real payoff. The derivation is careful, the equilibrium energy check passes, and the transverse mode coinciding exactly with QLCA is a nice consistency result. The comparison against MD for wavelengths near the interparticle spacing is genuinely useful, and the new longitudinal dispersion contains thermal terms QLCA misses, which explains the improved performance at weak coupling.\n\nThe main soft spot is not the step-function g, though that is crude. It is the closure in Eq. (3): the nonlocal term uses the equilibrium pair distribution function evaluated at local density and temperature, which is an assumption about out-of-equilibrium behavior. The authors know this—they say in the conclusion that at higher screening relaxation might behave differently. The 17% sound-speed error at κ=3, Γ=40 is consistent with that closure being wrong, not with the step function, and the paper's own MD-based check (which would remove the step-function and EOS effects) is mentioned but not documented. That is the weakest link. For a referee, I would ask for details of that check and for a direct test of the closure, e.g., computing the dispersion integrals with the MD pair distribution function.\n\nAlso worth noting: at moderate coupling the model visibly overestimates the MD peaks at κ=2, and QLCA fits better around Γ/Γm ≈ 0.17. The claim of \"best overall agreement across the full coupling range\" is fair only in an averaged sense. The sound-speed agreement at κ ≤ 1 is genuinely excellent, though.\n\nNo code or data are shipped, which limits reproducibility, but the method is described well enough to reimplement. The citation pattern looks fine; Ref. 22 is the source of the framework and the prior OCP validation, so leaning on it is appropriate.\n\nBottom line: this deserves peer review. It is a solid, honest paper with a clear new result and an acknowledged limitation. The equilibrium-g closure is the thing to scrutinize. I would send it to a serious referee and ask for the undocumented MD check to be reported; if that check holds up, the paper is in good shape.","headline":"Solid variational extension to Yukawa OCP with real finite-wavelength payoff; the main open question is the equilibrium-g closure, which the authors themselves flag.","tokens_in":21954,"tokens_out":1978,"would_cite":true,"duration_ms":22495,"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":"A variational Lagrangian with a step-function pair distribution reproduces the longitudinal dispersion of the Yukawa one-component plasma from weak coupling to near melting, including wavelengths comparable to the interparticle spacing.","keywords":["Yukawa one-component plasma","variational hydrodynamics","pair distribution function","step-function approximation","longitudinal dispersion relation","speed of sound","strong coupling","molecular dynamics"],"falsifier":"Replace the step-function $g$ in Eqs. (18)-(20) with a numerically exact $g(\\Gamma,\\kappa,x)$ obtained from molecular dynamics or an integral-equation theory, recompute the sound speed at $\\kappa_0=2,3$ for $\\Gamma$ near a few percent of $\\Gamma_m$, and compare with the molecular-dynamics data: if the 9-17% overestimate persists, the step function is not the cause; if it disappears, the closure is.","tokens_in":20921,"feed_emoji":"⚛️","tokens_out":8427,"duration_ms":84737,"temperature":0.7,"pith_summary":"The paper sets out to show that a single variational Lagrangian—built from kinetic energy, local thermodynamics, and a nonlocal interaction term carrying the pair distribution function—can serve as a closed hydrodynamic description of the Yukawa one-component plasma at strong coupling, down to wavelengths comparable to the interparticle spacing. It extends a framework devised for the Coulomb one-component plasma to the screened potential $\\varphi(r)=q^2 e^{-r/\\lambda}/(4\\pi\\varepsilon_0 r)$, derives the equations of motion, the momentum and energy conservation laws, and the linearized dispersion relations, and then evaluates the longitudinal branch using a step-function pair distribution whose cutoff radius is matched to a simulation-based internal-energy fit. The computed speed of sound agrees with molecular-dynamics data within a few percent for $\\kappa_0 \\le 1$, while the finite-wavelength dispersion reproduces the simulated mode peaks even for $|\\vec q_0|$ of order unity. A sympathetic reader would conclude that variational hydrodynamics is a quantitatively useful, computationally cheap alternative to particle simulations for short-wavelength collective modes in strongly coupled Yukawa plasmas.","feed_headline":"One Lagrangian predicts Yukawa plasma waves at atomic spacing","feed_subtitle":"Variational hydrodynamics with a step-function pair distribution matches simulated sound speeds and short-wavelength modes.","key_machinery":"The load-bearing object is the averaged Lagrangian of Eq. (3), whose nonlocal third term couples the displacement field to the screened Coulomb potential through the pair distribution function $g$; this is the only term that carries correlation physics beyond ideal-gas thermodynamics. Evaluation of the dispersion integrals then uses the step-function ansatz $g=0$ for $r<r_c(\\Gamma,\\kappa)$ and $g=1$ otherwise, with the dimensionless cutoff radius $x_c = r_c/a$ fixed by requiring the step-function excess energy of Eq. (24) to match the simulation-based fit of Eq. (21). That ansatz reduces the angular integrals $j$, $\\ell$, and $b$ of Eqs. (18)-(20) to the closed elementary expressions of Eqs. (25)-(27), which is what makes the dispersion law fast to evaluate and transparent enough to expose which thermodynamic inputs control the sound speed.","core_discovery":"On its own terms, the paper's central discovery is that the longitudinal dispersion relation $\\omega_L^2/\\omega_p^2$ obtained from the variational Lagrangian, evaluated with the step-function pair distribution, matches molecular-dynamics results across $0.5 \\le \\kappa_0 \\le 3$ and from weakly coupled to near-melting states. In the long-wavelength limit, the speed of sound is accurate to about 3% for $\\kappa_0 = 0.5$ and $\\kappa_0 = 1$ over the entire coupling range considered, while at $\\kappa_0 = 2$ and $\\kappa_0 = 3$ it overestimates the simulations by up to roughly 9% and 17%, respectively; the authors' own simulations indicate those discrepancies are not caused by the equation of state or by the step-function approximation. At finite wavenumbers, including $|\\vec q_0| \\sim 1$ where the wavelength is comparable to the interparticle spacing, the predicted dispersion tracks the peaks of the longitudinal current fluctuation spectrum over the full coupling range for $\\kappa_0 = 1$ and $\\kappa_0 = 2$. The paper claims this gives the best overall agreement with molecular dynamics among the analytic approaches it tests, while noting that other methods remain better in isolated regimes.","pith_inferences":["A direct test of the closure is to feed a numerically exact $g(\\Gamma,\\kappa,x)$ from molecular dynamics or an integral-equation theory into Eqs. (18)-(20) and see whether the overestimate at $\\kappa_0=2,3$ shrinks; if it does not, the step function is exonerated and the error lives in the out-of-equilibrium closure or in missing dissipation.","The pattern of accurate finite-wavelength dispersion alongside a biased long-wavelength sound speed at strong screening suggests the missing ingredient is wavenumber-selective, such as viscous damping or a relaxation-time correction that affects small $|\\vec q_0|$ most.","A pragmatic hybrid is plausible: pin the long-wavelength sound speed to the Euler+MF thermodynamic value, which is accurate at strong screening, and keep the variational finite-wavenumber terms; this combination has not been tested in the paper."],"forward_implications":["Because the derivation keeps the interaction potential general, the same Lagrangian can be applied to other strongly coupled fluids once an energy fit and a pair-distribution ansatz are supplied.","The transverse dispersion relation coincides with QLCA, so the variational framework inherits the established transverse-wave predictions for Yukawa systems while adding thermal and correlation corrections to the longitudinal branch.","At wavelengths comparable to the interparticle spacing, the theory tracks the molecular-dynamics current-fluctuation peaks, so it can be used to interpret collective-mode spectra in dusty-plasma and ultracold-neutral-plasma experiments without running particle simulations.","The known failure at strong screening is specific to the long-wavelength sound speed and low coupling; fixing it likely requires a revised out-of-equilibrium pair-distribution ansatz or the isothermal limit, both of which the paper identifies as next steps."],"supporting_citations":[{"why":"Supplies the variational Lagrangian, the equations of motion, and the step-function strategy that this paper extends to the screened Yukawa potential.","marker":"[22]"},{"why":"Provides the molecular-dynamics sound-speed data across coupling and screening regimes that the long-wavelength comparison is measured against.","marker":"[37]"},{"why":"Provides the practical excess internal energy and pressure fits used for the thermodynamic input of the dispersion law.","marker":"[34]"},{"why":"Gives the Yukawa equation-of-state and sound-velocity relations used in thermodynamic consistency and in the Euler+MF comparison.","marker":"[35]"},{"why":"Introduces the step-function pair-distribution approximation for Yukawa long-wavelength dispersion and shows that energy-matched and pressure-matched cutoffs give similar curves.","marker":"[36]"},{"why":"Provides accurate thermodynamic properties and a pair-distribution-based pressure check used to justify the free-energy route to the pressure.","marker":"[31]"},{"why":"Supplies the melting-coupling parameterization $\\Gamma_m(\\kappa)$ used to normalize coupling strengths and to define the range of the energy fit.","marker":"[56, 57]"},{"why":"The molecular-dynamics code used for the authors' own longitudinal current fluctuation spectra at finite wavelengths.","marker":"[63]"}],"fun_headline_variants":["Variational hydrodynamics predicts Yukawa plasma waves at atomic scale","Step-function pair distribution nails Yukawa plasma sound speed","Variational Lagrangian matches Yukawa plasma waves at interparticle spacing","Variational Lagrangian matches Yukawa plasma dispersion across screening regimes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the true pair distribution function can be replaced by the step function $g=0$ for $r<r_c(\\Gamma,\\kappa)$ and $g=1$ otherwise, with $r_c$ fixed by matching the internal-energy fit, so that all correlation structure beyond a single cutoff is discarded from the dispersion integrals.","fun_headline_variants_meta":{"raw":{"variants":["Variational hydrodynamics predicts Yukawa plasma waves at atomic scale","Step-function pair distribution nails Yukawa plasma sound speed","Variational Lagrangian matches Yukawa plasma waves at interparticle spacing","Variational Lagrangian matches Yukawa plasma dispersion across screening regimes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00131,"raw_usage":{"total_tokens":5364,"prompt_tokens":992,"completion_tokens":4372,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":4303}},"tokens_in":608,"tokens_out":4372,"duration_ms":30800,"temperature":1.0,"reasoning_tokens":4303,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:52:35.654613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Replace the step-function $g$ in Eqs. (18)-(20) with a numerically exact $g(\\Gamma,\\kappa,x)$ obtained from molecular dynamics or an integral-equation theory, recompute the sound speed at $\\kappa_0=2,3$ for $\\Gamma$ near a few percent of $\\Gamma_m$, and compare with the molecular-dynamics data: if the 9-17% overestimate persists, the step function is not the cause; if it disappears, the closure is.","supporting_citations":[{"cited_title":"Krimans and S","cited_arxiv_id":null,"evidence_quote":"Supplies the variational Lagrangian, the equations of motion, and the step-function strategy that this paper extends to the screened Yukawa potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the molecular-dynamics sound-speed data across coupling and screening regimes that the long-wavelength comparison is measured against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the practical excess internal energy and pressure fits used for the thermodynamic input of the dispersion law."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Yukawa equation-of-state and sound-velocity relations used in thermodynamic consistency and in the Euler+MF comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the step-function pair-distribution approximation for Yukawa long-wavelength dispersion and shows that energy-matched and pressure-matched cutoffs give similar curves."},{"cited_title":"Tolias, S","cited_arxiv_id":null,"evidence_quote":"Provides accurate thermodynamic properties and a pair-distribution-based pressure check used to justify the free-energy route to the pressure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The molecular-dynamics code used for the authors' own longitudinal current fluctuation spectra at finite wavelengths."}],"review_version":1}