{"id":"6bd9c605-4146-49be-aa20-bb967e607a03","arxiv_id":"2501.11113","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A 15-Lorentzian analytical Stark profile for H-alpha was fitted to 2992 simulated spectra, producing parameter tables for plasma diagnosis.","lead":"This paper builds an analytical formula for the shape of the hydrogen-alpha emission line in plasmas, using fifteen overlapping Lorentzian curves whose parameters are fitted to 2992 simulated spectra. If the fit holds up, plasma scientists could use it to extract electron density and temperature from simple measurements of this line.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No quantitative fit-quality metric is reported: the claim that Eq. (5) reproduces 2992 CS profiles rests on an unscaled FoM (Fig. 1) and a single visual comparison (Fig. 2).","rationale":"The reader's weakest assumption concerned the external validity of the CS benchmark data. My concern is complementary: even taking the CS profiles as truth, the paper does not demonstrate quantitatively that its model actually reproduces them. The FoM in Eq. (7) is a relative chi-square-like quantity, but no acceptance criterion, no error bars, and no residual analysis are provided. Figure 1 shows FoM values for all fits, but without a reference scale those numbers are uninterpretable; Figure 2 is a single eyeballed example. Since the model forces a common width and linear shifts, deviations from the CS profiles are to be expected, and the reported parameter tables inherit whatever bias results from those constraints. This is a load-bearing issue because the diagnostic use of the tables requires the fits to be accurate everywhere, not just on average. The proposed test would settle whether the fits are actually good and whether the linear-shift constraint is a meaningful source of bias. I do not think the paper's central idea is wrong, but the current evidence is insufficient to support the 'exact' label and the implied diagnostic reliability. The reader's CONDITIONAL verdict remains appropriate, so I recommend no change to the verdict, but the condition should include a quantitative validation of fit quality and a test of the linear-shift assumption.","tokens_in":9749,"tokens_out":3652,"duration_ms":29510,"concrete_test":"Select 100 CS profiles randomly from the 2992 conditions. For each, reconstruct Eq. (5) using the supplementary parameters and compute the normalized L2 residual ||P_model - P_CS||_2 / max(P_CS) over the spectral grid; report the maximum, median, and 95th percentile. If the maximum exceeds 5%, the 'exact' claim fails. In addition, refit the same 100 profiles with an unconstrained version where each d_k is an independent parameter (dropping Eq. 6) and test whether the best-fit d_k deviate systematically from s*k beyond the grid resolution; if they do, the linear-shift assumption biases the fitted s and areas.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the 15-Lorentzian model of Eq. (5) with the ten parameters fitted via Eq. (7) reproduces the CS Stark profiles across 2992 conditions—is not quantitatively established in the manuscript. The only aggregate measure presented is the FoM of Eq. (7), but no scale, threshold, or distribution of FoM values is given; Figure 1 plots values without a reference, and the only visual check (Figure 2) is a single case at µr = 6, Te = 24984 K, ne = 2.14×10^22 m^-3. The model imposes two strong constraints: a common width ω for all fifteen components and a strictly linear displacement d_k = s k (Eq. 6). These are physically motivated only approximately; in the low-density/low-µr regime the authors themselves report oscillations of s and ω (Figs. 3–5), suggesting the constraints are strained exactly there. Without residuals or an error metric (e.g., normalized L2 error over the spectral grid), the supplementary parameter tables may contain biased values whose diagnostic inversion would propagate unknown systematic errors into ne, Te, and Tg. The paper's own 'exact function' label is therefore unsupported by the evidence shown.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an analytical model for the Stark-broadened H-alpha line profile as a sum of fifteen Lorentzian components (Eq. 5), with a common width omega and a linearly k-dependent displacement d_k = s k (Eq. 6). The ten parameters {s, omega, a0, a1, ..., a6, a8} are fitted, using a genetic algorithm, to 2992 computer-simulated Stark profiles from Gigosos et al. The fitted parameters are tabulated in the supplementary material as functions of ne, Te, and the reduced fictitious mass mu_r. The authors analyze the parameter trends, interpret them in terms of ion and electron collision regimes, and propose the model as a diagnostic tool for extracting ne, Te, and Tg from measured H-alpha profiles.","tokens_in":10026,"tokens_out":3959,"duration_ms":41033,"significance":"If the model is as accurate as claimed and the parameter tables are reliable, this would be a practically useful resource for optical emission spectroscopy diagnostics, because it provides a closed analytical form for the H-alpha Stark profile over a wide range of plasma conditions. The use of the full fifteen-component splitting is physically motivated, and fitting all 2992 CS profiles is a substantial computational effort. The paper also makes a useful step by analyzing the parameter dependencies and relating them to plasma regimes. However, the central quantitative claim is not demonstrated in the manuscript: no scale-invariant fit metric, residual analysis, or independent validation is reported, and the 'exact' in the title overstates the status of what is a fitted approximation to a simulation model. The practical diagnostic value of the supplementary tables depends on validation that is currently missing.","major_comments":[{"comment":"The central claim that Eq. (5) reproduces the 2992 CS profiles is not quantitatively established. The only aggregate measure is the FoM of Eq. (7), which is an unscaled chi-square-like quantity, and Fig. 1 plots these values without a comparison scale, a threshold, or a statement of what constitutes an acceptable fit. The paper should report a normalized error metric (for example, normalized L2 error or mean absolute error over the spectral grid) with its distribution across all 2992 conditions, and give at least a few representative residuals. Without this, the fitted parameters in the supplementary tables cannot be assessed for bias, and the diagnostic inversion would propagate unknown systematic errors.","section":"Sec. 2.2, Fig. 1, Eq. (7)"},{"comment":"The statement that 'the results obtained confirm the hypothesis of the linear dependence of the displacements' is circular, because Eq. (6) imposes d_k = s k by construction; fitting a single slope s cannot test the linearity in k. The paper should either reframe this as a consistency check of the model assumption or test it by allowing the displacements d_k to vary independently at a few representative conditions and comparing the fitted values with a linear trend. Similarly, the shielding coefficient C quoted from Eq. (3) is back-calculated from the fitted slope and is not an independent prediction; the text should say so explicitly.","section":"Sec. 3, Fig. 3, Eq. (6)"},{"comment":"The title's phrase 'Exact Stark analytical function' overstates the status of the result. Equation (5) is a parametric model fitted to computer-simulated profiles, and the CS data themselves rest on approximations (e.g., the treatment of ion dynamics and microfield distributions). Please replace 'exact' with a more accurate qualifier such as 'analytical fit' or 'parametric model' throughout the manuscript, including the label 'Our model (exact function)' in Fig. 2.","section":"Title, Secs. 1 and 2.1"},{"comment":"The claim that 'the onset of stability is omega = 2 s = d2' is asserted but not quantitatively supported. Stability is not defined, and the relation is presented as an observation from a few plotted cases. The paper should define a stability criterion (for example, a tolerance on the change of parameters with mu_r) and demonstrate the relation statistically, or remove the claim. The oscillations in s and omega at low ne and low mu_r visible in Figs. 3-5 suggest that the model's constraints are strained in exactly that regime, which is not addressed by the current analysis.","section":"Sec. 3, Figs. 5-7"},{"comment":"The paper provides no validation against independent experimental Stark profiles or against an independent line-shape calculation. Since all fitted parameters are calibrated to the Gigosos et al. CS profiles [19,20], the diagnostic output of the model inherits any systematic errors in those simulations. At minimum, the authors should compare the model's predicted profiles with a few measured H-alpha profiles from the literature, or with the results of an independent Stark-broadening code, to demonstrate that the parameter tables are usable for experimental diagnosis.","section":"Sec. 2.1, Sec. 4"}],"minor_comments":[{"comment":"The genetic algorithm description is internally inconsistent: the text first states a maximum of 500 generations and later states a stopping criterion of 1000 generations. Please correct this.","section":"Sec. 2.2"},{"comment":"Table 2 gives relative intensities for the k-components in a static field, but the model treats the areas a_k as free fitting parameters and does not appear to use these values. Please clarify whether the table values are used as initial guesses, as constraints, or only as a physical reference.","section":"Sec. 2.1, Table 2"},{"comment":"The text calls epsilon_0 the 'vacuum permeability'; it should be the vacuum permittivity.","section":"Eq. (2)"},{"comment":"Figure 1 would be much more informative if it included labeled axes and perhaps a histogram of the FoM values rather than a point plot, since the current figure does not allow the reader to assess the distribution of fit quality.","section":"Fig. 1"},{"comment":"A few reference format issues appear: reference [17] is missing page numbers, and the title of reference [25] contains a typo ('ad excessive' should be 'and excessive').","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the underlying fitting effort is substantial. My main concern is evidentiary rather than conceptual: the central accuracy claim needs quantitative support and independent validation. The title and the 'confirmation' of linear displacements should be moderated. I therefore recommend major revision rather than rejection, since the required additions are within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: this is an incremental but careful extension of the same group's earlier one-, three-, and five-Lorentzian FFM fits to the full fifteen-component splitting of H-alpha, with fitted parameter tables for 2992 computer-simulated plasma conditions. If you do OES diagnostics of hydrogen plasmas, the supplementary tables give a ready-made closed-form profile to invert for ne, Te, and Tg. That is genuinely useful. The physical model is clearly laid out: common width, linear displacement slope, areas tied to line intensities, and calibration against Gigosos et al.'s simulations, which is the standard benchmark. The soft spots are real but not fatal. First, no quantitative fit-quality metric is reported. Figure 1 plots an unscaled FoM for all fits, with no threshold, distribution, or residuals; the only visual check is one profile. Saying the fits are satisfactory on that evidence is thin. A normalized error or worst-case residual plot would settle it. Second, the title's 'exact' is marketing. It is a fitted analytic approximation. Third, the paper says Figure 3 confirms the linear displacement hypothesis, but Eq. (6) imposes d_k = s k, so the fit cannot confirm the linearity it assumes. The data show that a linear displacement with a fitted slope reproduces the simulations well, which is useful but not confirmation. The extracted shielding coefficient C is back-calculated from the fitted slope and inherits the assumption. The stress-test note is right that there are no parameter uncertainties and no experimental validation. That limits confidence in inversion, especially at low density and low reduced mass, where the authors themselves report oscillations; those regions are exactly where the model is strained. None of this sinks the paper. The core claim—fifteen Lorentzians with ten fitted parameters reproduce the CS Stark profiles across a wide range—is plausible and supported by the evidence shown, even if not quantified as tightly as it should be. This is honest, competent work. I would send it to peer review with a request for residuals, error bars, an experimental cross-check, and a title change. If those come back, it is publishable and useful to the OES community.","headline":"Useful incremental 15-Lorentzian fit of H-alpha Stark profiles with a large parameter table, but the 'exact' label and the claimed confirmation of linear displacement overstate what is a calibrated fit.","tokens_in":780,"tokens_out":1631,"would_cite":false,"duration_ms":32333,"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":"Fifteen Lorentzian components with ten parameters reproduce all 2,992 simulated H-alpha Stark profiles in the paper's grid.","keywords":["Stark broadening","H-alpha line","plasma spectroscopy","Optical Emission Spectroscopy","Frequency Fluctuation Model","Lorentzian profiles","genetic algorithm fitting","plasma diagnostics"],"falsifier":"Take a measured H-alpha line from a plasma whose electron density and temperatures are known independently (by Thomson scattering or Langmuir probes), invert the line with the published parameter tables, and check whether the recovered ne, Te, and Tg match the independent values to within the model's reported accuracy; a systematic mismatch would show the parameterization does not transfer from simulated benchmarks to real plasmas.","tokens_in":9553,"feed_emoji":"⚛️","tokens_out":10275,"duration_ms":94324,"temperature":0.7,"pith_summary":"This paper tries to establish that the Stark-broadened shape of the hydrogen-alpha line can be represented exactly by a sum of fifteen Lorentzian profiles—one central component plus seven symmetric pairs—with only ten free parameters. Those parameters are fitted once, by a genetic algorithm, to the 2,992 computer-simulated Stark profiles that the paper takes as benchmark data over a wide grid of electron density, electron temperature, and gas temperature. If the representation holds, a measured H-alpha line can be inverted through the published parameter tables to recover ne, Te, and Tg, which is exactly what optical emission spectroscopy of hydrogen plasmas needs. The model is built on the Frequency Fluctuation Model and on the physical decomposition of the total width into electron- and ion-collision contributions, so the fitted parameters carry meaning about the plasma rather than being pure fitting coefficients.","feed_headline":"Fifteen Lorentzians fit 2,992 H-alpha Stark profiles","feed_subtitle":"A ten-parameter formula turns a measured H-alpha line into estimates of electron density and temperature.","key_machinery":"The machinery is the analytical profile Eq. (5): $$P_S = \\sum_{k=-8, \\, |k|\\neq 7}^{8} \\frac{2 a_{|k|}}{\\pi} \\frac{\\omega}{4(\\$\\lambda$ - d_k)^2 + \\$omega^{2}$},$$ with the linear shift $d_k = s k$ and the width $\\omega = \\omega_e(n_e, T_e) + \\omega_i(n_e, T_g)$. In the Frequency Fluctuation Model, the microfield around the emitting hydrogen atom fluctuates due to electron and ion collisions, producing Stark-split components that overlap constructively, destructively, or separately; the full set of fifteen components ($k = 0, \\pm1, \\pm2, \\pm3, \\pm4, \\pm5, \\pm6, \\pm8$) accounts for both $\\sigma$ and pi polarizations with relative intensities fixed by atomic wavefunctions. The ten fitted parameters are optimized with a genetic algorithm, and the reduced fictitious mass $\\mu_r = \\mu\\, T_e/T_g$ encodes ion mobility, which is what connects the fit to gas temperature.","core_discovery":"The central claim is that Eq. (5), the complete splitting of the H-alpha transition into fifteen Lorentzian components with a common width omega, linear shifts d_k = s k, and area coefficients a_|k|, reproduces the computer-simulated Stark profiles across all 2,992 plasma conditions studied. The ten temperature- and density-dependent parameters {s, omega, a0, a1, a2, a3, a4, a5, a6, a8} are fitted independently for each condition, and the paper reports them in supplementary material as a diagnostic library. From the fitted slope s the shielding coefficient C is estimated in the range 1.65 to 2.45, and the behavior of the parameters identifies regimes where ions are quasi-static versus dynamically mobile. The paper therefore claims to provide an analytical Stark function, based on the Frequency Fluctuation Model, that converts measured H-alpha profiles into plasma parameters without deconvolution of overlapping broadening mechanisms.","pith_inferences":["A natural next test, not performed in the paper, is to apply the published tables to experimental H-alpha profiles from plasmas with independently known ne, Te, and Tg (e.g., by Thomson scattering) to see whether the recovered parameters agree.","The same fifteen-Lorentzian template could be fitted to computer-simulated H-beta and H-gamma profiles to build comparable diagnostic tables for other Balmer lines, but the paper does not report such fits.","Because the fitting is done independently per plasma condition, the published tables could be interpolated or smoothed to give a continuous mapping from profile to parameters; the paper does not specify an interpolation scheme."],"forward_implications":["In any plasma within the fitted (ne, Te, Tg) range, a measured H-alpha line can be analyzed with the published tables to recover electron density, electron temperature, and gas temperature.","The parameter tables replace the need to run or store 2,992 simulated profiles, making the diagnostic cheap enough for real-time optical emission spectroscopy.","The slope-derived shielding coefficient C, between 1.65 and 2.45, gives a direct probe of electric-field shielding in the plasma.","The stability domains in the parameters mark where ion dynamics are negligible, telling a user when a simpler model would suffice."],"supporting_citations":[{"why":"supplies the computer-simulated H-alpha Stark profiles, including ion dynamics and non-equilibrium effects, that are the fitting targets for all 2,992 conditions.","marker":"[19,20]"},{"why":"introduces the Frequency Fluctuation Model on which the Lorentzian decomposition is based.","marker":"[14,15,16]"},{"why":"gives the collision-regime picture used to justify constructive and destructive superposition of split components.","marker":"[17]"},{"why":"provides the earlier reduced-Lorentzian analytical model and the linear shift relation extended in this work.","marker":"[24]"},{"why":"supplies the Stark shift and shielding coefficient framework used to interpret the slope parameter.","marker":"[25]"},{"why":"provides the genetic algorithm used for the parameter optimization.","marker":"[26]"}],"fun_headline_variants":["Ten-parameter formula fits 2,992 H-alpha Stark profiles","H-alpha line gives electron density and temperature","Genetic algorithm tunes H-alpha Stark function","Exact Stark analytical function: 15 Lorentzians","H-alpha broadening: 10 parameters, all plasmas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model's diagnostic value rests entirely on the accuracy of the computer-simulated Stark profiles used as fitting targets; if those simulations contain systematic errors, every fitted parameter inherits them, since no independent experimental validation is provided.","fun_headline_variants_meta":{"raw":{"variants":["Ten-parameter formula fits 2,992 H-alpha Stark profiles","H-alpha line gives electron density and temperature","Genetic algorithm tunes H-alpha Stark function","Exact Stark analytical function: 15 Lorentzians","H-alpha broadening: 10 parameters, all plasmas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002256,"raw_usage":{"total_tokens":8714,"prompt_tokens":940,"completion_tokens":7774,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":7699}},"tokens_in":556,"tokens_out":7774,"duration_ms":53481,"temperature":1.0,"reasoning_tokens":7699,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:36:57.060026+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a measured H-alpha line from a plasma whose electron density and temperatures are known independently (by Thomson scattering or Langmuir probes), invert the line with the published parameter tables, and check whether the recovered ne, Te, and Tg match the independent values to within the model's reported accuracy; a systematic mismatch would show the parameterization does not transfer from simulated benchmarks to real plasmas.","supporting_citations":[{"cited_title":"Griem, A.C","cited_arxiv_id":null,"evidence_quote":"gives the collision-regime picture used to justify constructive and destructive superposition of split components."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the earlier reduced-Lorentzian analytical model and the linear shift relation extended in this work."},{"cited_title":"Cvetanovic, S.S","cited_arxiv_id":null,"evidence_quote":"supplies the Stark shift and shielding coefficient framework used to interpret the slope parameter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the genetic algorithm used for the parameter optimization."}],"review_version":1}