{"id":"55dfed06-cab5-42c7-a159-44d53fb400b6","arxiv_id":"2608.07721","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Introduces PASCHEN-1D, a 1D drift-diffusion-Poisson plasma solver with multi-mechanism emission and circuit coupling, validated against published benchmarks, though the Paschen curve match requires pd-dependent fit parameters.","lead":"PASCHEN-1D is a new one-dimensional fluid plasma solver that couples drift-diffusion transport, electrode emission models, and flexible external circuit networks. The paper validates it against nanosecond pulsed discharges, DC breakdown, and photoemission-driven transients, but the Paschen curve match requires pressure-dependent fitting of emission and ionization parameters.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Paschen-curve benchmark is a parameter fit, not a predictive validation; the solver's central breakdown claim is therefore overstated.","rationale":"The paper is a software-description manuscript, and the core numerical machinery—drift-diffusion transport, Kurganov-Tadmor fluxes, RK4 integration, circuit coupling, and emission boundary conditions—is presented in sufficient detail to be plausible. The validation set includes several qualitatively different benchmarks, and the paper is transparent about many limitations in Section 5. The most load-bearing concern is not the fluid-model idealization itself (which the authors acknowledge and which is a known modeling choice) but the fact that one of the headline validations, the Paschen curve, is achieved by varying gamma and B with pd. This turns that benchmark from a test of the solver into a demonstration of sensitivity to fitted parameters, weakening the abstract's claim that results are consistent with published studies in a way that implies predictive capability. The reader's verdict of CONDITIONAL is appropriate: the solver is promising but not yet independently verified, and the Paschen benchmark should be reframed or replaced with a fixed-parameter test. I do not see an internal inconsistency that would justify rejection; the concern is about the weight of the validation evidence, not the soundness of the numerical method. The concrete test of running with fixed parameters would settle whether the full-curve agreement is a genuine prediction or a fit. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":19995,"tokens_out":7096,"duration_ms":74607,"concrete_test":"Request the code or run scripts from the authors, then rerun the Paschen benchmark with a single fixed gamma and the A,B values used at pd=1 Torr cm across the full pd range (0.24 to 5 Torr cm). Record the breakdown voltages obtained by the stated negative-differential-slope criterion. If the low-pd branch yields no breakdown or deviates strongly from the experimental curve, then the full-curve agreement in Fig. 6 should be explicitly labeled as a sensitivity study with pd-dependent fitted parameters, not as predictive validation. Additionally, report the exact gamma(pd) and B(pd) sequences used for the discrete markers in Fig. 7(b) so the fit can be audited against independent measurements of secondary emission yields and Townsend coefficients.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract states that PASCHEN-1D is validated against benchmarks including Paschen curve construction for argon and nitrogen, with results consistent with published studies. Section 4.4, however, shows that no single constant secondary-electron yield gamma can reproduce the experimental curve: agreement across the full pd range is recovered only by allowing gamma to increase on the low-pd branch and by reducing B on the high-pd branch. These are free parameters adjusted to match the data, so the reported agreement is partly constructed rather than a test of the solver's physics. The paper itself acknowledges this at the end of Section 4.4 and in Section 5, where it notes that reproducing breakdown over a wide pd range requires careful parameter selection and that no single coefficient set is guaranteed to work. The low-pd branch is especially fragile because the paper reports only 2–3 ionization events across the gap, making breakdown dominated by the fitted surface emission rather than by volumetric transport. Consequently, the Paschen benchmark cannot support the strongest claim of predictive validation; it is a sensitivity demonstration with pd-dependent effective parameters. The remaining benchmarks—nanosecond DBD, DC glow transition, and photoemission-driven discharge—are more credible but are qualitative and, in the photoemission case, compare against work by the same group. The central claim of a validated, self-consistent breakdown tool therefore rests on weaker evidence than the abstract implies.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"PASCHEN-1D is a one-dimensional drift-diffusion-Poisson fluid solver for gas discharges, with multi-mechanism electrode emission (ion-induced secondary, field, thermionic, and quantum photoemission) and flexible lumped-RLC external circuit coupling. The paper presents the governing equations, numerical methods (Kurganov-Tadmor fluxes, explicit diffusion, fourth-order Runge-Kutta, Picard iteration for density-Poisson coupling), and validation benchmarks: a nanosecond pulsed dielectric-barrier discharge, a numerical convergence study, DC Townsend-to-glow transitions, Paschen curves for argon and nitrogen, and an ultrafast photoemission-driven transient discharge. The central claim is that the solver is validated against published results with consistent agreement. However, the Paschen-curve agreement is obtained only after post hoc adjustment of the secondary-electron yield gamma and the Townsend coefficient B as functions of pd, which the paper acknowledges in Section 4.4 and Section 5. This weakens the predictive-validation claim for the Paschen benchmark, which is one of the headline validations.","tokens_in":20199,"tokens_out":8126,"duration_ms":71939,"significance":"If the solver is reliable, it addresses a real gap: many fluid plasma solvers lack self-consistent external circuit coupling and sophisticated emission physics. The modular architecture, the multiple emission mechanisms, and the unified circuit formulation are valuable. A clear strength is the convergence study in Section 4.2, which provides internal numerical consistency. The nanosecond-DBD and DC glow benchmarks are credible, though qualitative, and the photoemission case is compared against prior work by the same group and a WarpX PIC simulation. However, the most quantitative-looking validation, the Paschen curve, is effectively a parameter fit rather than a predictive test. The paper is honest about this limitation in its text, but the abstract and conclusions present the Paschen comparison as a validation benchmark without that qualification. With the validation claims recalibrated, the solver could be a useful community tool for breakdown and transient-discharge studies.","major_comments":[{"comment":"The Paschen-curve benchmark is not a predictive validation. The paper states (Section 4.4, after Fig. 7) that no single constant gamma can reproduce the experimental curve; agreement across the full pd range is recovered only by allowing the effective gamma to increase on the low-pd branch and reducing B on the high-pd branch. Since Eq. (34) shows the breakdown voltage is exponentially sensitive to gamma, adjusting these inputs to match the target data constitutes a parameter fit rather than a test of the solver's physics. The abstract and Section 6 list 'Paschen curve construction/reconstruction' among validation benchmarks without this qualification. Please reframe Section 4.4 explicitly as a sensitivity study of Paschen's law (or a demonstration of the solver's flexibility in exploring parameter dependence), and qualify the corresponding claims in the abstract and conclusions. This is load-bearing because the paper's central claim of being 'validated ... with results consistent with published studies' is overstated.","section":"Section 4.4 / Abstract / Section 6"},{"comment":"The convergence statement is stronger than the data support. Figure 3 shows that at fixed dt/dx, the quasi-steady electron density depends on the spatial resolution: N=500 and N=1000 overpredict the density relative to N=2000, with the text attributing this to numerical diffusion. The claim that 'all benchmark results presented in this study are obtained using discretization within the converged regime' requires reporting the specific spatial and temporal resolutions used in each benchmark (for example, the grid and time step in Section 4.1) and showing that these resolutions lie in the converged plateau for the observable of interest, or quantifying the discretization error. Without this, the blanket convergence claim is not substantiated for every benchmark.","section":"Section 4.2"},{"comment":"The photoemission-driven discharge validation is qualitative. The text states that PASCHEN-1D results 'are in agreement' with the experiments of Ref. [21] and with WarpX PIC simulations and that 'trends recovered' agree, but no quantitative error metrics are given: peak current values, voltage-sag magnitudes, or recovery times are not compared numerically. The same holds for the nanosecond-DBD benchmark in Section 4.1, where 'close agreement' is asserted without quantitative comparison. Please provide quantitative measures of agreement (e.g., relative errors in peak current, voltage, and timing) for these cases, or explicitly label them as qualitative comparisons that support code behavior but not quantitative predictive accuracy.","section":"Section 4.5 / Section 4.1"}],"minor_comments":[{"comment":"The sentence 'All benchmark cases are reproduced using parameters reported in the original references, without tuning or empirical fitting beyond those explicitly stated' is misleading, because Section 4.4 explicitly tunes gamma and B to match the experimental Paschen curve. Please revise this sentence to clarify that the Paschen comparison involves parameter adjustment.","section":"Section 4, opening paragraph"},{"comment":"Equation (17) is garbled in the typeset text (the expression for beta_2 appears as 'beta_2=2el epsilon, epsilon_3 L' with missing symbols). Please correct the equation and define all variables.","section":"Eq. (17)"},{"comment":"The paper reports that Eq. (34) gives a negative breakdown voltage of -147.98 V at pd=0.2 Torr cm for A=11.5, B=176, and gamma=0.05. A direct evaluation of the standard Paschen formula with these constants gives approximately -16 V. Please check the arithmetic and the parameter values used.","section":"Eq. (34) and Section 4.4"},{"comment":"Section 3.3 states that first-order (piecewise constant) reconstruction is used in the Kurganov-Tadmor scheme, which involves no slope limiter, yet Section 4.1 mentions setting 'a slope limiter theta=1.01 in the Kurganov-Tadmor Scheme.' Please reconcile this contradiction, either by removing the slope-limiter parameter from the benchmark description or by clarifying that the implementation supports higher-order reconstruction with slope limiting but first-order was used for the reported results.","section":"Section 3.3 vs. Section 4.1"},{"comment":"The Einstein-relation closure is typeset as D_i = mu_i k(T_i/e, missing the closing parenthesis. Please correct this to D_i = mu_i k T_i / e.","section":"Section 2.4.1"},{"comment":"Table 1 lists PASCHEN-1D as 'Open-source: Yes' and the text cites it as open source, but no repository URL or code-availability statement is provided anywhere in the manuscript. Please add the distribution information.","section":"Table 1"},{"comment":"The caption says 'resolution-independent plateau,' but the plateau is visibly dependent on N: N=500 and N=1000 give higher densities than N=2000. Please change the wording to 'N-dependent plateau' or 'resolution-dependent plateau' and discuss the implications.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The paper is an honest code-description manuscript, and the authors explicitly acknowledge the Paschen-curve limitation in Section 4.4 and Section 5. The main issue is that the abstract and conclusions do not carry this nuance, presenting the Paschen benchmark as a validation when it is in fact a sensitivity study with pd-dependent fitted parameters. I believe this is fixable within the scope of a revision by reframing the Paschen section and qualifying the validation language. The paper would also benefit from an explicit code-availability statement, since it advertises the solver as open source. The fit with the journal's scope is acceptable for a numerical-methods and validation paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid software-description for a 1D drift-diffusion-Poisson fluid solver with unusually complete external-circuit coupling and a flexible menu of surface emission models. The integration is genuinely useful: most public fluid codes treat the circuit crudely or bolt emission on as a fixed yield, and here the lumped RLC network, dielectric layers, and several emission mechanisms are handled in one modular framework. The equations are standard, the Kurganov-Tadmor/RK4 numerics are sensible, and the convergence study in Section 3 gives me reasonable confidence the code is internally consistent.\n\nWhat is new is the artifact, not the physics, and that is fine. But the central validation claim needs reining in. The abstract says the solver is validated against Paschen curve construction for argon and nitrogen. Section 4.4 shows something weaker: with a fixed parameter set the solver fails on both branches, and agreement is recovered only by letting gamma and B vary with pd. The paper is admirably transparent about this, and the discussion of why low pd makes breakdown emission-dominated (2–3 ionization events across the gap) is physically sensible. But that means the Paschen benchmark is a sensitivity demonstration with fitted effective parameters, not a predictive validation. The other benchmarks are more credible yet still qualitative: the nanosecond DBD reproduces trends from Adamovich et al., and the photoemission case compares against the same group's earlier work. None include quantitative error metrics.\n\nThe bigger practical problem is the open-source claim. Table 1 marks it \"Yes,\" and the paper is clearly intended as a community tool, but no repository, license, or download link appears anywhere. A software paper that ships code and reproducible scripts deserves credit; this one asks the reader to trust that the code exists. That should be fixed before publication.\n\nI would not call the Paschen section circular in a damaging sense, because the authors themselves flag it as a limitation and explain the sensitivity analysis they performed. The flaw is in the framing, not in the science. I would send this to peer review with a clear request: soften the abstract's validation language, reframe Section 4.4 as a sensitivity study rather than a validation, and make the code publicly available. With those changes it becomes a genuinely useful resource for low-temperature plasma modelers.\n\nFor a reading group, I would bring it up as an example of an honest software paper with a textbook overclaim in the abstract. I would not cite it yet, but I would revisit once the code is available.","headline":"Useful, honest software paper whose core validation claim is overstated: the Paschen benchmark is a pd-dependent parameter fit, and the open-source claim lacks a repository.","tokens_in":20791,"tokens_out":1262,"would_cite":false,"duration_ms":14853,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.80.-s","52.65.Kj","52.25.Fi"],"model":"deepseek-v4-flash","headline":"A one-dimensional fluid solver, PASCHEN-1D, claims to reproduce gas-discharge benchmarks by self-consistently coupling drift-diffusion transport, multi-mechanism electrode emission, and external RLC circuits.","keywords":["fluid plasma solver","drift-diffusion-Poisson","gas discharge breakdown","surface electron emission","external circuit coupling","Paschen curve","dielectric-barrier discharge","ultrafast photoemission"],"falsifier":"A direct test would be to measure the low-pd breakdown voltage of argon in a well-characterized cell with clean electrodes of known secondary-emission yield, and compare it with a PASCHEN-1D run that uses the measured yield rather than a fitted one; if the solver cannot reproduce the left branch of the Paschen curve with independently measured surface parameters, the claim that it captures breakdown physics rather than fitted coefficients would fail. Alternatively, a kinetic particle-in-cell simulation with identical cross-sections and surface yields would serve the same purpose.","tokens_in":19754,"feed_emoji":"⚡","tokens_out":8005,"duration_ms":73193,"temperature":0.7,"pith_summary":"PASCHEN-1D is a one-dimensional time-dependent fluid plasma solver that solves drift-diffusion equations for electrons and ions together with Poisson's equation, while dynamically coupling to lumped RLC circuits and to multiple electrode emission mechanisms. The paper's central claim is that this self-consistent coupling lets a single code reproduce published behavior across nanosecond pulsed dielectric-barrier discharges, DC Townsend-to-glow transitions, argon and nitrogen Paschen curves, and ultrafast photoemission-driven discharges. The reason a reader should care is that breakdown and transient discharge behavior is often controlled precisely by the interplay of surface emission, space charge, and circuit impedance, effects that most fluid codes treat separately or crudely. If the claim holds, PASCHEN-1D provides a low-cost, physically transparent platform for studying and designing pulsed and DC discharge systems.","feed_headline":"1D plasma solver ties emission and circuits to breakdown behavior","feed_subtitle":"Self-consistent drift-diffusion model matches benchmarks from pulsed discharges to Paschen curves.","key_machinery":"The load-bearing object is the coupled drift-diffusion-Poisson system itself, advanced by a finite-volume discretization: continuity equations $\\partial n_s/\\partial t + \\partial \\Gamma_s/\\partial x = S_s$ with $\\Gamma_e = -D_e \\partial n_e/\\partial x - \\mu_e n_e E$ and $\\Gamma_i = -D_i \\partial n_i/\\partial x + \\mu_i n_i E$, a Poisson equation for the potential, and source terms combining Townsend ionization $\\nu_i = \\alpha \\mu_e |E|$ with electron-ion recombination. The argument is carried by the two feedback loops built on this system: plasma fluxes produce a circuit current, the circuit returns the gap voltage imposed as a Dirichlet boundary condition, and the electrode emission module converts surface physics into an electron flux boundary condition, with a Picard iteration keeping the density boundary conditions consistent with the electric field. Transport and ionization coefficients come from user closures or from local-field interpolation of swarm-data tables, so the solver is modular across gases.","core_discovery":"The discovery is a numerical one: a drift-diffusion-Poisson solver with Kurganov-Tadmor fluxes, explicit fourth-order Runge-Kutta time stepping, and a Picard iteration between density boundary conditions and the electrostatic field can hold together a much wider loop than usual. The loop includes ion-induced secondary emission, Fowler-Nordheim and Murphy-Good field emission, Richardson-Dushman thermionic emission, and quantum-mechanical photoemission at the electrodes, plus an external circuit whose instantaneous gap voltage is fed back into the Poisson boundary condition and whose current is computed from the space-averaged charged-particle fluxes via the Ramo-Shockley relation. On the strength of the benchmark cases, the paper claims that this unified framework reproduces the reference results without tuning beyond the parameters reported in the original studies, and that the numerical convergence study shows the reported agreement reflects the modeled physics rather than discretization artifacts.","pith_inferences":["A natural extension the authors do not pursue is replacing pd-dependent effective $\\gamma$ with a single energy-dependent secondary emission model; if such a model could reproduce the left Paschen branch without tuning, it would strengthen the fluid approach in regimes where only a few ionization events occur per gap.","Because the circuit module accepts arbitrary RLC combinations, the solver could be used to co-design external drive waveforms and discharge geometry, for example sweeping shunt inductors or series capacitors to find operating points for pulsed-power or thruster applications.","The precomputed quantum photoemission pulse at the cathode suggests a route to modeling multi-pulse laser trains or RF-modulated photoemission by concatenating emission pulses, which would connect the solver to laser-triggered breakdown and high-repetition-rate experiments.","The modular transport interface means adding an energy-equation closure or multiple ion species could be done without altering the validated plasma-circuit coupling; such extensions would likely be required before the solver could reliably model molecular-gas chemistry or capacitively coupled discharges."],"forward_implications":["The solver reproduces nanosecond pulsed dielectric-barrier discharges, including avalanche onset, space-charge screening, and dielectric-mediated field redistribution, using parameters taken from the reference benchmark.","For DC discharges in argon at pd = 1 Torr cm, PASCHEN-1D produces pre-breakdown, Townsend, and glow regimes and identifies breakdown at the onset of negative differential gap-voltage slope, with microsecond relaxation to quasi-steady plateaus.","Argon and nitrogen Paschen curves are reproduced across the tested pd range, but only when the ion-induced secondary emission yield is allowed to vary with pd; the low-pd branch is governed by surface emission and the high-pd branch by the Townsend $\\alpha$-process.","The ultrafast photoemission benchmark shows the code bridging picosecond emission physics and microsecond plasma-circuit recovery, matching measured voltage and current transients along with independent particle-in-cell simulation results.","A convergence study on the quasi-steady mid-gap electron density shows resolution-independent behavior for sufficiently small $\\Delta t/\\Delta x$, supporting the claim that benchmark agreement is physical rather than numerical."],"supporting_citations":[{"why":"Supplies the nanosecond pulsed dielectric-barrier discharge benchmark, including geometry, dielectric layers, and nitrogen mobility coefficients, and provides the dielectric-mapping formulation for circuit coupling.","marker":"[33]"},{"why":"Provides the experimental argon and nitrogen Paschen curves and the Townsend/rate coefficients used in the breakdown and glow benchmarks.","marker":"[2]"},{"why":"Supplies the classical Townsend ionization coefficients and the argon transport coefficients used for DC breakdown, glow transition, and Paschen curve runs.","marker":"[1]"},{"why":"Provides the ultrafast photoemission-driven discharge experiment, including its R-C circuit parameters and measured voltage/current transients used as a benchmark.","marker":"[21]"},{"why":"Supplies the Kurganov-Tadmor central-upwind flux scheme that stabilizes the drift-dominated electron and ion transport.","marker":"[53]"},{"why":"Gives the Ramo-Shockley theorem used to convert space-averaged charged-particle flux profiles into an instantaneous circuit current.","marker":"[34–37]"},{"why":"Provides the quantum-mechanical photoemission model used for ultrafast cathode emission injection in the photoemission benchmark.","marker":"[44–48]"},{"why":"Defines the swarm-data table format the solver interpolates for local-field electron transport coefficients.","marker":"[30]"}],"fun_headline_variants":["PASCHEN-1D unifies emission, circuits, and plasma breakdown","1D plasma solver couples emission physics to external circuits","Breakdown sim with multi-mechanism emission and circuit feedback","New fluid solver matches pulsed discharges and Paschen curves","Drift-diffusion code links electrodes, circuits, and gas dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the drift-diffusion approximation with local-field transport coefficients and a single effective ion species adequately describes the electrons and ions in all benchmark regimes, a premise the paper itself flags as fragile at low pd, where the ionization mean free path becomes comparable to the gap and only two to three ionization events occur per gap.","fun_headline_variants_meta":{"raw":{"variants":["PASCHEN-1D unifies emission, circuits, and plasma breakdown","1D plasma solver couples emission physics to external circuits","Breakdown sim with multi-mechanism emission and circuit feedback","New fluid solver matches pulsed discharges and Paschen curves","Drift-diffusion code links electrodes, circuits, and gas dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1481,"prompt_tokens":979,"completion_tokens":502,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":414}},"tokens_in":595,"tokens_out":502,"duration_ms":5952,"temperature":1.0,"reasoning_tokens":414,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:21:41.664807+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to measure the low-pd breakdown voltage of argon in a well-characterized cell with clean electrodes of known secondary-emission yield, and compare it with a PASCHEN-1D run that uses the measured yield rather than a fitted one; if the solver cannot reproduce the left branch of the Paschen curve with independently measured surface parameters, the claim that it captures breakdown physics rather than fitted coefficients would fail. Alternatively, a kinetic particle-in-cell simulation with identical cross-sections and surface yields would serve the same purpose.","supporting_citations":[{"cited_title":"=−D\"∂n\"∂x+µ","cited_arxiv_id":null,"evidence_quote":"Supplies the classical Townsend ionization coefficients and the argon transport coefficients used for DC breakdown, glow transition, and Paschen curve runs."}],"review_version":1}