{"id":"f5da5301-8c96-4c7c-bf89-7625a9fb1549","arxiv_id":"2411.18358","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A six-term local heat flux closure, with coefficients linked to wave growth or damping and to linear kinetic theory, explains over 95% of the heat flux variance in 1D electrostatic electron plasma simulations.","lead":"Sparse regression on particle-in-cell simulations of 1D electrostatic plasmas identifies a compact six-term formula that captures more than 95% of the heat flux variation in Landau-damped Langmuir waves and two-stream instabilities. A reader might care because interpretable data-driven closures could let fluid models capture kinetic wave-particle effects at much lower cost than full kinetic simulations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The six-term closure may be an artifact of the hand-selected local library; since §2.1 calls the restriction arbitrary and linear theory (§A) indicates k-dependent, effectively nonlocal terms, the library's sufficiency is a load-bearing, untested premise.","rationale":"After reading the full paper, I find the reader's weakest assumption to be the correct point of entry. The term library is described in §2.1 as an 'arbitrary choice' made for computational convenience, and the paper's own linear-theory analysis shows that matching the kinetic heat flux perturbation requires coefficients that depend on wavenumber and growth rate (Eq. 3.3). This means the discovered six-term local closure is at best a local approximation for a narrow band of modes, not a demonstrated general closure. The low FVU values and the recovery of the momentum equation in Appendix B show the SR pipeline works on these data, but they do not show that the restricted library is sufficient; FVU is measured within the same library, and the momentum recovery uses a different target variable. I therefore propose one concrete test: enlarge the library with E, second derivatives, and nonlocal representations, and check whether the six terms survive. This is exactly the sensitivity test the paper declines to run. The reader's CONDITIONAL verdict already anticipates this class of concern, so I do not change it. Credit is due for the linear-theory consistency checks and the careful data-alignment validation, which support the pipeline's correctness even if the library sufficiency remains open.","tokens_in":25756,"tokens_out":8328,"duration_ms":82259,"concrete_test":"Rerun the SR pipeline on the two-stream nb/ne=0.1 and the vth=0.01c Langmuir datasets with an augmented library adding E and ∂xE, second-order derivative terms (e.g., vth³ ∂x²n, n vth² ∂x²V, n vth ∂x²vth), and a nonlocal (Fourier-space) representation of the first-derivative terms. Record selected terms and FVU on held-out folds. If the six original terms remain and FVU changes by <1 percentage point, the library restriction is not the source of the found closure; if new terms are selected or FVU drops by >1 pp, the six-term closure is library-dependent and the central claim must be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central discovery claim rests on the assumption that the candidate library (Eq. 2.5 plus a few first-derivative terms, chosen 'arbitrarily' in §2.1) contains the true closure. This is tested only implicitly: low FVU and linear-theory checks are necessary but not sufficient, because both use the same restricted library. Linear collisionless theory (Appendix A, Eq. A8) gives q̃ ∝ r n̄ v̄th² ω/k, i.e. a 1/k nonlocality. To emulate this with local first-derivative terms, coefficients A2, A3, A6 must absorb k and γ, and indeed the constraints (Eq. 3.3) make A1 and kA6 explicit functions of k, ωr, γ. A local closure with mode-dependent coefficients is not a transferable closure. If the library were enlarged to include E (which in 1D is the inverse-gradient of n), second derivatives, or nonlocal operators, the six terms might change or FVU might drop substantially. The paper acknowledges this possibility (§4) but performs no such test. Thus the claim that these six terms are 'physically relevant' is not yet established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses sparse regression (SINDy/PDE-FIND) on OSIRIS particle-in-cell data of one-dimensional electrostatic electron dynamics, covering Landau-damped Langmuir waves and two-stream instabilities, to discover a local heat-flux closure. It reports a six-term model q = A1 n vth^2 V + A2 vth^3 ∂x n + A3 n vth^2 ∂x vth + A4 + A5 n vth^3 + A6 n vth^2 ∂x V, with fraction-of-variance-unexplained (FVU) values between 0.5% and 7% depending on case. The authors quantify term importance, compare coefficients with constraints derived from linear collisionless theory, identify a possible seventh term ∝ n vth V^2, and validate the regression workflow by recovering the electron momentum equation in Appendix B.","tokens_in":26079,"tokens_out":4173,"duration_ms":39285,"significance":"If the discovered closure is robust, it offers a parsimonious, interpretable local closure that could inform collisionless fluid models and complement nonlocal closures such as Hammett-Perkins. The paper has clear strengths: the linear-theory constraints in Sec. 3.4/Appendix A are derived rather than fitted, the recovery of the momentum equation in Appendix B demonstrates that the SR workflow correctly identifies known relations when the library contains them, and the treatment of PIC field staggering and cross-validation is careful. However, the significance is tempered by two issues: the candidate term library is explicitly acknowledged in Sec. 2.1 as an arbitrary restriction, and the central claim of 'discovering' a closure is based on goodness of fit within that library rather than on predictive validation against an enlarged function space or truly out-of-sample regimes.","major_comments":[{"comment":"The central claim that the six terms in Eq. (3.1) are 'physically relevant' rests on the assumption that the hand-selected candidate library (dimensionally consistent products n vth^α V^(3−α) plus a small set of first-order derivative terms) is sufficiently expressive to contain the true closure. Section 2.1 explicitly calls this restriction 'an arbitrary choice,' and the library excludes E and B, second and higher spatial derivatives, and nonlocal operators. The reported FVU values are computed within the same restricted library, so they demonstrate only that the six terms are the best local first-derivative fit available, not that they capture the dominant physics. A concrete robustness test is needed: repeat the regression with an enlarged library that includes E, ∂x^2 n, ∂x^2 vth, ∂x(V) products, and perhaps a few nonlocal/bilinear terms, and show whether the six-term structure persists and how much the FVU changes. Without such a test, the discovery claim is not fully established.","section":"§2.1, Eq. (2.5)"},{"comment":"Equation (3.3) and the surrounding discussion show that the coefficients A1–A6 are explicit functions of k, ωr, γ, and vth. This means the discovered 'closure' is not a parameter-free local closure: applying it in a fluid code would require knowledge of the mode wavenumber and growth/damping rate, which are typically not known a priori. The paper acknowledges this in §4 but does not explain how the closure could be used predictively. Moreover, the linear-theory comparisons in Figs. 7–8 insert SR-fitted coefficients (e.g., A5 or A2, A3) into the theoretical expressions to predict other SR-fitted coefficients (A1 or A6); this is a consistency check between two quantities derived from the same data, not an independent out-of-sample validation. The authors should either demonstrate predictive use of the closure (e.g., in a fluid simulation) or clearly state that the result is a data-driven parametrization of q for diagnosed mode parameters, not a transferable closure.","section":"§3.4, Eq. (3.3) and §4"},{"comment":"The closure coefficients are reported without uncertainty bounds. Given the particle noise in PIC data and the fact that 10-fold cross-validation is used only for term selection, the paper should quantify the sampling variability of A1–A6 and C. For example, Fig. 2b and Fig. 3a show coefficient trends versus vth and nb/ne with no error bars, and statements such as 'some coefficients are close to zero' or 'near-zero coefficients are generally not found consistently' lack a quantitative threshold. Without confidence intervals, the reader cannot assess whether the observed parameter dependences are significant or whether the claimed six-term form is stable under repeated sampling.","section":"§3.1–§3.2, Figs. 2–3"}],"minor_comments":[{"comment":"The sentence 'we then collect the volume-integrated y and θj quantities are over these small volumes' contains a grammatical error ('are' should be removed).","section":"§2.1, Eq. (2.1)"},{"comment":"The phrase '10^5 (10^4) electrons (ions) per cell' is ambiguous; it should read '10^5 electrons and 10^4 ions per cell'.","section":"§2.2.1"},{"comment":"The shorthand Φ± is used in Eq. (3.3) but not defined until Eq. (3.4); the definition should appear before first use.","section":"§3.4, Eq. (3.3)"},{"comment":"The axis label 'γ × 10ω−1 pe' is confusing; specify the units of γ explicitly (e.g., 'γ/ωpe (multiplied by 10)').","section":"Fig. 4"},{"comment":"In the discussion of Cheng et al. (2023), the phrase 'the the local approximation' contains a duplicated article; it should read 'the local approximation'.","section":"§1, paragraph 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of J. Plasma Phys. and the authors are transparent about the arbitrary library choice and the parametric dependence of the coefficients. The main weakness is that the discovery claim is not yet supported by a robustness test against an enlarged library or by a genuine out-of-sample prediction. I recommend major revision with a request for such a test; the manuscript is otherwise well written and the analytical linear-theory constraints are a valuable contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a genuine advance over prior SR closure work: it covers both Landau damping and two-stream instability, identifies a consistent six-term local closure, and derives linear-theory constraints (Appendix A) that the coefficients satisfy to first order. Second, the biggest caveat is the one the authors state themselves in Section 2.1: the candidate term library is an arbitrary restriction. The discovery claim is therefore conditioned on that library containing the true closure.\n\nWhat's good: the method validation is real. Appendix B recovers the electron momentum equation to sub-percent coefficient accuracy, and shows the E-field staggering correction. That gives confidence the pipeline works. The FVU values (0.5-7%) and the deltaFVU term-importance ranking are useful. The linear-theory constraints (Eq. 3.3) are derived, not fitted, and the check in Figs. 7-8 is a reasonable consistency test, though it uses SR-fitted coefficients for some terms to predict others, so it is not independent validation.\n\nWhere it's soft: the term library issue. Section 2.1 literally calls the restriction 'an arbitrary choice.' Linear theory (Appendix A, Eq. A8) gives q~ proportional to omega/k, i.e., a 1/k nonlocality. To mimic that with local first-derivative terms, coefficients A2, A3, A6 must absorb k and gamma; indeed Eq. (3.3) makes A1 and kA6 explicit functions of k, omega_r, gamma. So the discovered closure is not a transferable closure—it is a local fit with mode-dependent coefficients. The paper acknowledges this in Section 4, but does not test a larger library (e.g., including E, second derivatives, or nonlocal operators). That is a real gap, not a fatal one. The authors are transparent about it, and for the stated purpose—an interpretable starting point for sub-grid closure—the six terms are a reasonable candidate.\n\nOther minor issues: coefficients are reported without uncertainty bounds; no code or data shipped; the two-stream case combines counter-streaming populations into a single species, which limits the growth-phase closure's practical use (they note this). These are addressable.\n\nBottom line: this deserves a serious referee. It is a credible, honest contribution to reduced plasma modeling, with a central claim that is plausible but not fully proven. I would send it out, and I would want the referee to push on library sensitivity and coefficient uncertainties.","headline":"A credible, honest sparse-regression closure paper whose six-term discovery is real but conditioned on an arbitrary term library that the authors flag but never test; worth sending to a serious referee.","tokens_in":26571,"tokens_out":2211,"would_cite":true,"duration_ms":20501,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.35.Fp","52.25.Dg","52.65.Rr"],"model":"deepseek-v4-flash","headline":"The paper claims that sparse regression consistently discovers a six-term local heat flux closure that accounts for more than 95% of the heat flux variation in one-dimensional electrostatic plasma phenomena.","keywords":["sparse regression","heat flux closure","Landau damping","two-stream instability","collisionless plasma","particle-in-cell simulation","Langmuir waves","fluid closure"],"falsifier":"Train the same sparse-regression pipeline on the same simulations with a candidate library that adds second-order derivative terms such as $n v_{\\rm th}\\partial_x^2 V$ and $v_{\\rm th}^2\\partial_x^2 n$, and electric-field terms such as $nE$ and $n v_{\\rm th}E$; if any of these are consistently selected with a material drop in FVU, the six-term local closure would be an artifact of the hand-picked library. A second concrete check is to compare the closure's predictions against OSIRIS output in a two- or three-dimensional electrostatic simulation, where the $k$-odd terms would have to acquire directional structure.","tokens_in":25601,"feed_emoji":"⚡","tokens_out":11472,"duration_ms":86931,"temperature":0.7,"pith_summary":"This paper sets out to show that in one-dimensional electrostatic plasmas, the heat flux—the quantity that a collisionless fluid model must add by hand to close its moment equations—can be replaced by a short, local formula built from the density, flow velocity, and thermal speed and their first spatial derivatives. Applying sparse regression to particle-in-cell data from Landau-damped Langmuir waves and two-stream instabilities, the authors find the same six terms recurring across cases, together accounting for more than 95% of the variation in the heat flux. The result is relevant because collisionless fluid models are far cheaper than kinetic simulations but need closures, and most analytic closures are nonlocal or rely on assumptions such as linearity or adiabaticity that break down in turbulence and reconnection. A compact, interpretable local closure is directly usable in fluid solvers and sub-grid models.","feed_headline":"Six local terms capture 95% of plasma heat flux variation","feed_subtitle":"Compact local formula reproduces kinetic heat flux in Landau-damped waves and two-stream instabilities.","key_machinery":"The machinery is sparse regression in the SINDy/PDE-FIND style, applied to small space-time volume integrals of the fluid moments so that particle noise is reduced. The candidate library is deliberately restricted to dimensionally consistent products $n v_{\\rm th}^{\\alpha} V^{3-\\alpha}$ plus selected first-order derivative combinations such as $v_{\\rm th}^3\\partial_x n$, $n v_{\\rm th}^2\\partial_x v_{\\rm th}$, and $n v_{\\rm th}^2\\partial_x V$, along with a constant; the paper states that this restriction is an arbitrary choice made to keep the library tractable. Ten-fold cross-validation is used to keep only terms that appear consistently across folds, and the resulting formula is checked against two analytic constraints derived from linear collisionless Vlasov theory that relate the coefficients to the wave frequency, growth rate, and wavenumber.","core_discovery":"On the paper's own terms, the central discovery is that the electron heat flux $q$ in these 1D electrostatic settings is captured by $q = A_1 n v_{\\rm th}^2 V + A_2 v_{\\rm th}^3 \\partial_x n + A_3 n v_{\\rm th}^2 \\partial_x v_{\\rm th} + A_4 + A_5 n v_{\\rm th}^3 + A_6 n v_{\\rm th}^2 \\partial_x V$, where $n$, $V$, and $v_{\\rm th}$ are the electron density, flow velocity, and thermal speed. The six coefficients split into $k$-even and $k$-odd groups: the even group is independent of propagation direction, while the odd group switches sign when the wave direction is reversed and vanishes for standing waves or symmetric counter-streaming setups. Sparse regression returns this six-term form with an unexplained-variance fraction of 2–7% for Landau-damped Langmuir waves and 0.5–5% for two-stream instabilities, i.e. regular capture of more than 95% of the heat-flux variation. The coefficients $A_2$, $A_3$, and $A_6$ are approximately proportional to the instantaneous growth or damping rate, and the recovered values satisfy two constraint relations derived from linear collisionless Vlasov theory, Eq. (3.3). In the saturated phase of the two-stream instability, where net growth vanishes, a three-term form with $A_1$, $A_4$, and $A_5$ suffices; since $A_4$ does not affect $\\partial_x q$, only two terms remain important for pressure evolution.","pith_inferences":["The paper leaves the parameterization of the coefficients as future work; a natural extension is a quasi-local closure in which $A_2$, $A_3$, and $A_6$ are explicit functions of the instantaneous growth rate and characteristic wavenumber, which would make the closure portable across phases of an instability.","The same regression could be run on the beam and core electron populations separately; the paper reports preliminary similarity in the recovered terms, and separate closures may avoid the destabilizing sign of $A_3$ that appears when a single combined electron species is modelled in the two-stream-unstable case.","The six-term expression is a natural physics-informed ansatz for machine-learned closures: using it as the functional form shrinks the parameter space compared with unrestricted neural networks while retaining interpretability and a built-in guarantee of localness.","A direct fluid-code test would settle practical value: implement Eq. (3.1) with growth-phase coefficients in a 10-moment solver, run past saturation, and compare long-time behaviour with OSIRIS; the paper notes that small coefficient errors can seed numerical instabilities, so stability, not just FVU, is the decisive test."],"forward_implications":["A 10-moment fluid solver for 1D electrostatic problems can use Eq. (3.1) directly as a closure and inherit more than 95% of the kinetic heat-flux variance without solving the Vlasov equation.","During saturated phases with no net wave growth, the closure collapses to $q = A_1 n v_{\\rm th}^2 V + A_4 + A_5 n v_{\\rm th}^3$; because $A_4$ does not appear in $\\partial_x q$, only two terms drive the pressure evolution in that regime.","Because $A_2$, $A_3$, and $A_6$ track the instantaneous growth or damping rate, a fixed-coefficient closure is phase-specific; making these coefficients functions of $\\gamma$ and a characteristic wavenumber would extend the closure across the growth-to-saturation transition.","The $k$-odd terms imply that any three-dimensional generalization must involve a direction vector, so isotropic or standing-wave settings retain only the $k$-even terms.","The seventh term $\\propto n v_{\\rm th} V^2$ acts as a second-order braking term that becomes important when $|\\gamma|$ is decreasing, i.e. near saturation of an instability rather than in deep linear growth."],"supporting_citations":[{"why":"Introduces the SINDy algorithm used for sparse regression.","marker":"Brunton et al. 2016"},{"why":"Presents the PDE-FIND variant that the implementation builds on for discovering PDEs from data.","marker":"Rudy et al. 2017"},{"why":"Supplies the modified sparse-regression workflow and its earlier application to reduced plasma models.","marker":"Alves & Fiuza 2022"},{"why":"Defines the collisionless Landau-damping closure whose local approximation is the paper's main theoretical reference.","marker":"Hammett & Perkins 1990"},{"why":"Gives the local approximant of Hammett-Perkins with a characteristic wavenumber, motivating the gradient terms in the candidate library.","marker":"Sharma et al. 2006"},{"why":"Provides an improved ten-moment closure used as the local Hammett-Perkins comparison point.","marker":"Ng et al. 2020"},{"why":"Closest prior work, recovering the local Hammett-Perkins closure in a Landau-damped Langmuir standing-wave setup.","marker":"Cheng et al. 2023"},{"why":"Earlier sparse-regression closure discovery for reconnection, which avoided the linear regime that this paper focuses on.","marker":"Donaghy & Germaschewski 2023"},{"why":"Documents the OSIRIS particle-in-cell code used to generate the kinetic data.","marker":"Fonseca et al. 2002, 2008"}],"fun_headline_variants":["Data-driven heat flux closure: six terms do the job","Sparse regression finds six-term plasma heat flux closure","Six local terms explain 95% of plasma heat flux","Heat flux closure from sparse regression: six terms suffice","Six-term closure captures plasma heat flux variation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the hand-selected candidate library—dimensionally consistent products of $n$, $v_{\\rm th}$, and $V$ plus a few first-order derivative terms—is rich enough to contain the true heat-flux dependence; the paper itself calls that restriction an arbitrary choice, so if the real closure needs higher-order derivatives, nonlocal operators, or explicit electric- and magnetic-field terms, the six recovered terms would be an artifact of the preselected library rather than the dominant physics.","fun_headline_variants_meta":{"raw":{"variants":["Data-driven heat flux closure: six terms do the job","Sparse regression finds six-term plasma heat flux closure","Six local terms explain 95% of plasma heat flux","Heat flux closure from sparse regression: six terms suffice","Six-term closure captures plasma heat flux variation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000194,"raw_usage":{"total_tokens":1420,"prompt_tokens":1075,"completion_tokens":345,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":691,"completion_tokens_details":{"reasoning_tokens":283}},"tokens_in":691,"tokens_out":345,"duration_ms":3794,"temperature":1.0,"reasoning_tokens":283,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:16:48.642153+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Train the same sparse-regression pipeline on the same simulations with a candidate library that adds second-order derivative terms such as $n v_{\\rm th}\\partial_x^2 V$ and $v_{\\rm th}^2\\partial_x^2 n$, and electric-field terms such as $nE$ and $n v_{\\rm th}E$; if any of these are consistently selected with a material drop in FVU, the six-term local closure would be an artifact of the hand-picked library. A second concrete check is to compare the closure's predictions against OSIRIS output in a two- or three-dimensional electrostatic simulation, where the $k$-odd terms would have to acquire directional structure.","supporting_citations":[{"cited_title":", Proctor, Joshua L","cited_arxiv_id":null,"evidence_quote":"Introduces the SINDy algorithm used for sparse regression."},{"cited_title":", Brunton, Steven L","cited_arxiv_id":null,"evidence_quote":"Presents the PDE-FIND variant that the implementation builds on for discovering PDEs from data."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the modified sparse-regression workflow and its earlier application to reduced plasma models."},{"cited_title":"& Perkins, Francis W","cited_arxiv_id":null,"evidence_quote":"Defines the collisionless Landau-damping closure whose local approximation is the paper's main theoretical reference."},{"cited_title":", Quataert, Eliot & Stone, James M","cited_arxiv_id":null,"evidence_quote":"Gives the local approximant of Hammett-Perkins with a characteristic wavenumber, motivating the gradient terms in the candidate library."},{"cited_title":"Journal of Plasma Physics 89 (1), 895890105","cited_arxiv_id":null,"evidence_quote":"Earlier sparse-regression closure discovery for reconnection, which avoided the linear regime that this paper focuses on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the OSIRIS particle-in-cell code used to generate the kinetic data."}],"review_version":1}