{"id":"54f720f9-783d-436a-aa15-0620c72340e3","arxiv_id":"2502.04851","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Higher gas flow suppresses the ozone-to-nitrogen-oxide mode transition in a grid-type surface dielectric barrier discharge, an effect the model attributes to reduced vibrationally excited nitrogen.","lead":"The paper measures ozone, nitrogen dioxide, and dinitrogen pentoxide in a surface dielectric barrier discharge at gas flows from 0 to 10 liters per minute, and fits a chemistry model to the data. It finds that higher gas flow suppresses the switch from ozone-producing to nitrogen-oxide-producing chemistry, which matters for tuning plasma reactors used in water treatment and agriculture.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central N2(v) mechanism is an input, not a tested result; fitted tau_v trend also contradicts a convective-loss picture.","rationale":"The reader's weakest assumption—that Eqs. 4-5 introduce T0_v and tau_v as unconstrained fit parameters and that the flow dependence of N2(v) is an input rather than a measurement—is exactly the load-bearing concern. My stress-test sharpens it in two ways: (1) a thermal mechanism via gas temperature is explicitly shown by the authors' own sensitivity analysis to cause earlier transition, yet Tg/Tp are held essentially constant, so the model excludes a plausible alternative without justification; (2) the fitted tau_v increases with flow, which is opposite to the behavior expected if convective removal of N2(v) were the physical cause. Neither point is ad hominem; both are internal to the model and its stated sensitivities. The paper is transparent about non-unique fits and about the model's limitations, and the experimental measurements appear solid, so the conditional verdict is appropriate. The proposed re-fit would settle whether the N2(v) attribution is necessary or merely one of several acceptable parameterizations. No change to the reader's CONDITIONAL verdict is warranted.","tokens_in":17591,"tokens_out":5186,"duration_ms":53715,"concrete_test":"Refit the model to the measured O3, NO2, and N2O5 time series for all flow rates with T0_v and tau_v fixed at their 0-slm values, allowing instead Tg, Tp, and nO to vary with flow using a physically motivated relation (e.g., reduced specific power deposition and gas heating as flow increases). If this alternative reproduces the flow-dependent suppression of the mode transition with comparable or better agreement, the claim that N2(v) density is the key factor is unsupported. Also test whether a single tau_v set by the shortest residence time (about 1.5 s at 10 slm) with T0_v as the only N2(v) parameter can match the data; if it cannot, the prescribed tau_v trend is unjustified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that higher gas flow suppresses the O3-to-NOx mode transition by reducing the density of N2(v>12)—is not independently supported. nN2(v>12) is never measured and is not a transported species in the model; it is prescribed by Eq. 4-5 with T0_v and tau_v freely adjusted per flow rate in Table III. Figure 9 therefore shows an input, not a derived or measured quantity. The authors acknowledge the fits are non-unique and that the Maxwellian vibrational distribution is unlikely. A thermal alternative is not excluded: the model's own sensitivity analysis states that increasing Tg promotes earlier transition, yet Tg is held at 310 K and Tp only changes 335→340 K across all flow rates; gas heating at low flow could explain the observed mode transition without invoking N2(v). Moreover, Table III shows tau_v increases from 18 s at 0 slm to 35 s at 6-10 slm. If flow removed vibrationally excited N2 by convection, the effective rise time of Tv should shorten, not lengthen. The longer tau_v at high flow is a fit knob that delays the transition but has no physical support in the model. Because the key quantity is manually imposed, the agreement between measured and simulated species densities cannot validate the mechanistic claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports optical absorption spectroscopy measurements of O3, NO2, and N2O5 densities in a grid-type surface dielectric barrier discharge operated in dry synthetic air at flow rates from 0 to 10 slm. At low flow rates the system transitions from an O3-dominated to an NOx-dominated mode over tens of seconds, while at high flow rates no clear transition is observed. A two-zone zero-dimensional chemical kinetics model with eight free parameters is fitted to the measured density profiles. The authors attribute the flow-rate dependence of the mode transition to a flow-induced decrease in the density of vibrationally excited N2(v>12), represented in the model by a decreasing steady-state vibrational temperature T0_v and an increasing time constant tau_v with flow. The model reproduces O3 dynamics well at all flow rates, while NO2 and N2O5 agreement degrades at higher flow.","tokens_in":17873,"tokens_out":5120,"duration_ms":47019,"significance":"The experimental data set, covering absolute densities of key reactive species over a range of flow rates and operating times, is a useful contribution and the optical absorption methodology is well established. The two-zone model description is transparent, the reaction set is clearly documented, and the authors are commendably explicit about the non-uniqueness of their fitted parameters and the questionable validity of the Maxwellian vibrational distribution. If the central mechanistic claim were independently supported, the result would be significant for controlling SDBD mode transitions in applications. However, the claim as presented is not sufficiently supported: the key quantity nN2(v>12) is prescribed by hand-fitted parameters rather than derived from or measured against the flow physics, and alternative mechanisms such as gas heating are not excluded. The paper's strength is the honest and detailed presentation of a model that can be tuned to reproduce the data, not yet a validated mechanism.","major_comments":[{"comment":"The central mechanistic claim, that increasing gas flow suppresses the O3-to-NOx mode transition by reducing nN2(v>12), is not supported by the model as presented. nN2(v>12) is prescribed by Eqs. (4)–(5) with T0_v and tau_v as free parameters, and Table III shows T0_v declining from 6500 K at 0 slm to 4500 K at 8–10 slm while tau_v increases from 18 s to 35 s. Figure 9 therefore displays an input of the fit, not a derived or measured quantity. Because the flow dependence of nN2(v>12) is entirely imposed by these hand-set parameters, the agreement between measured and simulated O3 densities in Figs. 10–12 cannot validate the proposed mechanism; it only demonstrates that the model can be tuned to reproduce the data. The authors' own statement that the parameter combinations are not unique further weakens the inference. To support the claim, the model would need to derive the temporal evolution of nN2(v>=12) from a transport or kinetics equation that couples gas flow to vibrational excitation/loss, or the vibrational density would need to be measured independently.","section":"§IV (Eqs. 4–5), Table III, Fig. 9"},{"comment":"The sensitivity discussion states that an increase in T_g leads to an earlier mode transition, yet in Table III T_g is held at 310 K for all flow rates and T_p changes only from 335 K to 340 K. Longer residence times at low flow would be expected to increase gas heating, providing an alternative explanation for the early transition at low flow that the model does not test. Because the fitted gas temperatures are nearly flow-independent by construction, the model cannot discriminate between the proposed N2(v) mechanism and a thermal mechanism. A test that varies T_g or T_p while holding T0_v fixed would be required to show that vibrational excitation, rather than gas temperature, is the controlling factor.","section":"§V.B, Table III"},{"comment":"The fitted time constant tau_v increases from 18 s at 0 slm to 35 s at 6–10 slm. If higher gas flow reduced the density of vibrationally excited N2 by convective removal, one would expect the effective rise time of T_v to shorten, not lengthen, because the loss rate increases with flow. The monotonic increase of tau_v with flow appears to be an ad hoc adjustment that delays the transition at high flow rather than a physically motivated parameter; no rationale is given for why convection should slow the vibrational excitation or relaxation process. This inconsistency further undermines the credibility of the fitted parameter trends and, with them, the central claim.","section":"Table III, tau_v trend"},{"comment":"The measured density profiles are presented without any error bars, confidence intervals, or uncertainty propagation from the Lambert–Beer analysis (Eq. 3). Given that the mode transition is identified from the temporal behavior of the measured curves and that the model is fit to these curves, the absence of uncertainty estimates leaves unquantified how clearly the transition is distinct at 2 slm versus 4 slm, and whether the simulated fits are within experimental error. The experimental foundation of the flow-dependent disappearance of the transition would be far more convincing if measurement uncertainties were reported.","section":"§III.A, Figs. 6–7"}],"minor_comments":[{"comment":"The conclusion uses both \"N2(v) with v ≥ 12\" and \"N2(v >12)\" inconsistently; pick one notation and use it throughout the manuscript.","section":"§VI"},{"comment":"The notation \"kflow;p,g\" uses a semicolon that is puzzling; also the subscript \"p,g\" could be confused with \"pg\". Please clarify or rename this coefficient.","section":"§IV, Eq. (8)"},{"comment":"Units are typeset inconsistently: \"0.0 ms-1\" and \"0.1 m -1\" should be \"0.0 m s^-1\" and \"0.1 m s^-1\" respectively, with proper spacing and exponents.","section":"§II.A, §V.A"},{"comment":"The data availability statement appears with garbled spacing and repeated words (\"A V AILABILITY OF DAT A ST A TEMENT...\"); this should be cleaned up before publication.","section":"Data Availability Statement"}],"recommendation":"major_revision","confidential_remarks":"The experimental data are likely reproducible and the model is described in sufficient detail, but the conclusions overstate the support for the vibrational mechanism. I would encourage the authors to either substantially reframe the central claim as a hypothesis consistent with the model, or add a discriminating test (e.g., varying gas temperature independently or deriving nN2(v>12) from a transport model) before publication. The paper may be better suited to a specialty plasma chemistry journal after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a read for the dataset. The authors measured O3, NO2, and N2O5 densities over time in a grid-type SDBD at flows from 0 to 10 slm, and they clearly document the known O3-to-NOx transition at low flow and its suppression at high flow. The two-zone zero-dimensional model is a sensible extension of the Shimizu and Park approaches, and the authors are unusually candid that the parameter fits are non-unique, that the Maxwellian vibrational distribution is unlikely to be physically exact, and that the model really only tracks O3 well.\n\nThe soft spot is exactly where the stress-test lands. The central claim—that higher flow suppresses the transition by reducing nN2(v>12)—is not a prediction; it is an input. T0_v is hand-set from 6500 K at 0 slm to 4500 K at 8-10 slm, and tau_v from 18 to 35 s. Figure 9 therefore shows the assumption, not a derived or measured quantity. The authors admit this in the text, but the abstract and conclusion still frame it as a finding, which overstates the evidence.\n\nThe tau_v trend is an additional real problem that the paper does not address. If flow removes vibrationally excited N2 by convection, the effective rise time of Tv should get shorter with increasing flow, not longer. The fitted tau_v increases by roughly a factor of two, which points in the opposite direction. That inconsistency deserves an explanation or the mechanism loses force. Also, the sensitivity analysis says raising Tg promotes earlier transition, yet Tg is fixed at 310 K at all flows; gas heating at low flow is a plausible alternative that is not excluded by the measurements.\n\nThat said, the experimental result itself is probably solid. The flow-rate dependence of the mode transition is clearly shown, and the reaction-balance analysis gives useful insight into which channels matter for O3 and NOx. Missing error bars and the manual fitting procedure are minor-to-moderate issues for a measurements-first plasma chemistry paper.\n\nWho benefits: applied groups working on SDBDs for ozone or nitrogen oxide production, VOC treatment, and anyone modeling these reactors. They will use the flow data and the reaction scheme even if they treat the N2(v) mechanism as a hypothesis.\n\nRecommendation: engage with it. It deserves serious review, but the referee should push for a clean separation between fitted inputs and tested predictions, for a physically constrained treatment of tau_v, and for error bars on the absorption measurements. A revised version that presents the N2(v) dependence as an assumption consistent with the data, rather than a demonstrated cause, would be honest and still useful.","headline":"Good new flow-rate data for a grid SDBD, but the paper's central N2(v) mechanism is fitted input, not a tested result.","tokens_in":18450,"tokens_out":2739,"would_cite":true,"duration_ms":29326,"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 gas flow rate through a grid-type surface dielectric barrier discharge controls whether the plasma's chemistry settles into an ozone-dominated or a nitrogen-oxide-dominated mode, via the density of vibrationally…","keywords":["dielectric barrier discharge","reactive oxygen and nitrogen species","plasma chemistry","optical absorption spectroscopy","chemical kinetics simulation","mode transition","vibrational excitation","gas flow rate"],"falsifier":"Measure the density of $N_2(v \\ge 12)$ directly in the reactor at 0 slm and 10 slm using state-selective spectroscopy (e.g., vacuum-ultraviolet absorption or two-photon laser-induced fluorescence). The claimed mechanism predicts a significantly lower density at 10 slm at the time the transition would be expected; if the density is not lower, or if injecting pre-vibrationally excited $N_2$ into a high-flow gas stream does not restore the transition, the proposed flow-rate–vibrational-excitation link is falsified.","tokens_in":17404,"feed_emoji":"⚡","tokens_out":11300,"duration_ms":95768,"temperature":0.7,"pith_summary":"This paper tries to establish that the gas flow rate in a grid-type surface dielectric barrier discharge operated in dry synthetic air determines which of two chemical modes the plasma occupies: an ozone-dominated mode at high flow, and a nitrogen-oxide-dominated mode that emerges after a transition time at low flow. The authors measure $O_3$, $NO_2$, and $N_2O_5$ densities by optical absorption spectroscopy at flow rates from 0 to 10 slm and find that the $O_3$-to-$NO_x$ transition occurs at 0-2 slm but is absent at 4-10 slm. A two-zone zero-dimensional chemical kinetics model reproduces the measured $O_3$ dynamics and attributes the flow-rate effect to a lower density of highly vibrationally excited nitrogen, $N_2(v \\ge 12)$, represented in the model by a reduced vibrational temperature at higher flow. If the mechanism is correct, flow rate provides a simple external control for plasma chemistry relevant to ozone generation, pollution remediation, and biomedical applications.","feed_headline":"Higher gas flow suppresses the ozone-to-NOx mode switch","feed_subtitle":"Low flow lets vibrationally excited nitrogen flip the chemistry; high flow holds the ozone mode.","key_machinery":"The central object is the time-dependent vibrational temperature $T_v(t)$ of nitrogen, defined by $T_v = T_g + T_v^0(1 - \\exp(-t/\\tau_v))$, whose Maxwellian tail gives the density $n_{N_2(v \\ge 12)} = n_{N_2}\\exp(-12\\Delta\\varepsilon_v / (k_B T_v))$. This density drives the key reaction R9, $O + N_2(v) \\rightarrow NO + N$, and the NO produced then consumes $O_3$ via R18, $O_3 + NO \\rightarrow NO_2 + O_2$, switching the system from the ozone mode to the $NO_x$ mode. The two-zone model couples a plasma-near region and a gas-far region through diffusion, plasma-induced drift, and a flow-removal term, and its eight free parameters ($n_O$, $T_v^0$, $\\tau_v$, $T_g$, $T_p$, $n_{N_2(A^3\\Sigma)}$, $n_{O_2(a^1\\Delta)}$, $r_N$) are fitted to the measured densities at each flow rate.","core_discovery":"The paper's central discovery is that the mode transition from the $O_3$ to the $NO_x$ state in this surface dielectric barrier discharge is governed by the gas flow rate through the density of $N_2(v \\ge 12)$, the source of NO via $O + N_2(v) \\rightarrow NO + N$. Experimentally, the transition is observed at low flow rates and not at high flow rates. The model captures this by assigning each flow rate a steady-state vibrational temperature $T_v^0$ and a rise time $\\tau_v$: $T_v^0$ decreases from 6500 K at 0 slm to 4500 K at 5-10 slm, and $\\tau_v$ increases from 18 s to 35 s, so that the computed $n_{N_2(v \\ge 12)}$ is lower at high flow. The simulations follow the $O_3$ rise, peak, and decay at 0 slm within about 50% agreement, while the $NO_2$ and $N_2O_5$ dynamics deviate more strongly at high flow, where the role of vibrationally excited nitrogen is diminished. The authors conclude that this mechanism is the key factor for the absence of the mode transition at higher flow rates, and note that N-atom production becomes a potentially significant pathway for $NO_2$ formation in that regime.","pith_inferences":["If the Maxwellian vibrational distribution assumption in equations (4)-(5) is replaced by a more realistic non-equilibrium distribution, the fitted absolute values of $T_v^0$ would change, but the qualitative flow-rate dependence of $n_{N_2(v \\ge 12)}$ could persist; the paper's mechanism is therefore testable at the level of the vibrational distribution shape itself.","The model's underprediction of $NO_2$ at high flow suggests that atomic nitrogen, not vibrationally excited $N_2$, may dominate $NO_x$ production in the high-flow regime; a direct measurement of N-atom densities (for example by two-photon absorption laser-induced fluorescence) would determine whether the proposed mechanism is incomplete.","The same residence-time argument may carry over to other plasma-chemical converters: any means of shortening the effective gas residence time—pulsed flow, oscillating gas velocity, or reactor geometry—might suppress the mode transition in the same way as raising the flow rate."],"forward_implications":["At high gas flow rates the discharge remains in the ozone-dominated mode, so raising the flow is a direct way to favor $O_3$ production over $NO_x$ production in applications such as water treatment and sterilization.","The mode transition is set by the slow build-up of vibrationally excited nitrogen (time constant $\\tau_v$ of tens of seconds), so operating the discharge in short bursts—or with residence times shorter than $\\tau_v$—should hold the ozone mode even at low flow.","Reactions R9 and R18 form the nonlinear switch pair: any measure that lowers the rate of $O + N_2(v) \\rightarrow NO + N$ or that accelerates $O_3$ consumption will shift the transition time.","The fitted parameters show that suppressing the steady-state vibrational temperature from 6500 K to 4500 K is sufficient to prevent the transition, giving a quantitative target for flow-based mode control."],"supporting_citations":[{"why":"Supplies the original semi-empirical model linking the ozone-to-NOx mode transition to the build-up of vibrationally excited N2, which this work extends to a two-zone geometry.","marker":"[16]"},{"why":"Supplies the extension of the mode-transition model with fixed densities of N2(A3Σ) and O2(a1Δ) and the rapid-switch chemistry this work adopts.","marker":"[18]"},{"why":"Supplies the two-zone transport treatment (diffusion between plasma-near and gas-far regions) used in the model equations.","marker":"[31]"},{"why":"Supplies the optical absorption spectroscopy methodology, reactor setup, and species-specific filter choices used for the density measurements.","marker":"[26]"},{"why":"Provides the basis for interpreting 193 nm absorption as N2O5 rather than N2O4, which fixes how the measured N2O5 densities are derived.","marker":"[22]"},{"why":"Supplies evaluated rate coefficients for neutral species reactions in the chemical kinetics scheme.","marker":"[35]"},{"why":"Supplies the temperature-dependent rate coefficients for ozone and nitrogen oxide reactions, including the key O3 + NO and O + NO2 reactions.","marker":"[36]"}],"fun_headline_variants":["Flow rate flips discharge chemistry from ozone to NOx","Excited nitrogen powers flow-sensitive plasma modes","Low flow stokes nitrogen vibrations, favoring NOx","Surface plasma mode switch hinges on gas velocity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that higher flow suppresses the ozone-to-NOx transition by reducing the density of vibrationally excited nitrogen is an input to the model, not a measured fact: the vibrational temperature and its rise time are fitted separately for each flow rate, and the assumed Maxwellian shape of the vibrational distribution is acknowledged in the paper to be experimentally unlikely.","fun_headline_variants_meta":{"raw":{"variants":["Flow rate flips discharge chemistry from ozone to NOx","Excited nitrogen powers flow-sensitive plasma modes","Low flow stokes nitrogen vibrations, favoring NOx","Surface plasma mode switch hinges on gas velocity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":3004,"prompt_tokens":1035,"completion_tokens":1969,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":651,"completion_tokens_details":{"reasoning_tokens":1909}},"tokens_in":651,"tokens_out":1969,"duration_ms":15940,"temperature":1.0,"reasoning_tokens":1909,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T21:12:00.917527+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the density of $N_2(v \\ge 12)$ directly in the reactor at 0 slm and 10 slm using state-selective spectroscopy (e.g., vacuum-ultraviolet absorption or two-photon laser-induced fluorescence). The claimed mechanism predicts a significantly lower density at 10 slm at the time the transition would be expected; if the density is not lower, or if injecting pre-vibrationally excited $N_2$ into a high-flow gas stream does not restore the transition, the proposed flow-rate–vibrational-excitation link is falsified.","supporting_citations":[{"cited_title":"Shimizu , author Y","cited_arxiv_id":null,"evidence_quote":"Supplies the original semi-empirical model linking the ozone-to-NOx mode transition to the build-up of vibrationally excited N2, which this work extends to a two-zone geometry."},{"cited_title":"Park , author W","cited_arxiv_id":null,"evidence_quote":"Supplies the extension of the mode-transition model with fixed densities of N2(A3Σ) and O2(a1Δ) and the rapid-switch chemistry this work adopts."},{"cited_title":"Sakiyama , author D","cited_arxiv_id":null,"evidence_quote":"Supplies the two-zone transport treatment (diffusion between plasma-near and gas-far regions) used in the model equations."},{"cited_title":"Schücke , author A","cited_arxiv_id":null,"evidence_quote":"Supplies the optical absorption spectroscopy methodology, reactor setup, and species-specific filter choices used for the density measurements."},{"cited_title":"\\ Huh , author J","cited_arxiv_id":null,"evidence_quote":"Provides the basis for interpreting 193 nm absorption as N2O5 rather than N2O4, which fixes how the measured N2O5 densities are derived."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies evaluated rate coefficients for neutral species reactions in the chemical kinetics scheme."}],"review_version":1}