Pith. sign in

REVIEW 5 major objections 5 minor 26 references

Exact Stark analytical function for H{\alpha} line based on the FFM Model related with plasma parameters

T0 review · 5 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Fifteen Lorentzian components with ten parameters reproduce all 2,992 simulated H-alpha Stark profiles in the paper's grid.

desk verdict 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. read the letter →

arxiv 2501.11113 v1 pith:XNG2VZ2Z submitted 2025-01-19 physics.plasm-ph physics.atom-ph

A. Sarsa , A. Jiménez-Solano , M.S. Dimitrijević , C. Yubero This is my paper · ORCID
classification physics.plasm-phphysics.atom-ph
keywords StarkbroadeningH-alphalineplasmaspectroscopyOpticalEmissionFrequencyFluctuationModelLorentzianprofilesgeneticalgorithmfittingdiagnostics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 5 minor

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.

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 (5)
  1. [Sec. 2.2, Fig. 1, Eq. (7)] 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.
  2. [Sec. 3, Fig. 3, Eq. (6)] 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.
  3. [Title, Secs. 1 and 2.1] 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.
  4. [Sec. 3, Figs. 5-7] 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.
  5. [Sec. 2.1, Sec. 4] 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.
minor comments (5)
  1. [Sec. 2.2] 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.
  2. [Sec. 2.1, Table 2] 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.
  3. [Eq. (2)] The text calls epsilon_0 the 'vacuum permeability'; it should be the vacuum permittivity.
  4. [Fig. 1] 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.
  5. [References] 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').

Circularity Check

2 steps flagged · score 4.0 of 10

The core 15-Lorentzian parameterization is a fit to external CS data, so it is not circular, but the paper's confirmation of the linear displacement law d_k = s k is built in by construction, and the extracted shielding coefficient C is a back-calculation from the fitted slope.

  1. self definitional [Section 2.1, Eq. (6) and Section 3, Figure 3]
    "because of Equation (3), a linear dependence of the displacements, dk, with k is assumed, dk = s k (6)... In Figure 3 we show the results for the slope, s... The results obtained confirm the hypothesis of the linear dependence of the displacements, Eq. (6)."

    The linear dependence is not a result of the fit; it is imposed by Eq. (6), where s is one of the ten fitted parameters. Since no independent estimates of the individual displacements d_k are reported, plotting the fitted s cannot confirm Eq. (6). Any successfully fitted profile built with d_k = s k will, by construction, display exactly the linear ordering the paper claims to verify. The confirmation is therefore equivalent to the model's own constraint.

  2. fitted input called prediction [Section 3, paragraph following Figure 3]
    "Starting from Eq. (3) and using the values of the slope, the shielding parameter of the electric field can be estimated, obtaining values in the range 1.65 ≤ C ≤ 2.45."

    Eq. (3) gives d_k ≈ 10^{-16} C k n_e^{2/3}, and Eq. (6) sets d_k = s k with s fitted at each plasma condition. Therefore C = s / (10^{-16} n_e^{2/3}) is just a rescaling of the fitted slope, not an independently determined shielding coefficient. Presenting the resulting range as an estimate adds no new information and does not test Eq. (3), since the form d_k ∝ k n_e^{2/3} was already assumed in the fit.

full rationale

The core of the paper is an explicit parameterization: ten parameters of a 15-Lorentzian model are fitted, condition by condition, to 2992 computer-simulated Stark profiles of Gigosos et al. In that respect the central claim, that the analytical expression reproduces the CS benchmark across the grid, is not circular: the benchmark is external, and the fit is a legitimate interpolation. The cited prior works [21-24] are the same group's earlier reduced-Lorentzian fits, but they are used as motivation rather than as the sole justification of the new model. The circularity concerns are limited to two secondary interpretive steps: Fig. 3 is presented as confirming the linear displacement law although Eq. (6) builds that law into the model, and the shielding coefficient C is back-computed from the fitted slope rather than independently extracted. Because these steps are interpretive rather than the core benchmark-fitting claim, overall circularity is partial and moderate. The absence of a quantitative fit-quality metric is a correctness or validation concern, not a circularity finding.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central result is a fit, not a derivation: all ten parameters per plasma condition are obtained by genetic algorithm fitting to external CS data. The model rests on the FFM physical picture and the Gigosos benchmark, but no new entity is introduced.

free parameters (4)
  • s (slope of dk = s k) = Varies over 2992 conditions; values in supplementary material
    Fitted by genetic algorithm for each (rho, ne, mu_r); sets the displacement of the Stark components.
  • omega (common FWHM) = Varies over 2992 conditions; values in supplementary material
    Fitted for each condition; assumed identical for all fifteen Lorentzians.
  • a0 (central component area) = Varies over 2992 conditions; values in supplementary material
    Fitted area for k=0 Lorentzian.
  • a1, a2, a3, a4, a5, a6, a8 (component areas) = Varies over 2992 conditions; values in supplementary material
    Fitted areas for the symmetric k-components; seven distinct values because symmetric pairs share the same area.
assumptions (5)
  • domain assumption The CS model of Gigosos et al. [19,20] provides accurate Stark-only benchmark profiles for the H-alpha line.
    Invoked in Section 2.1 as the data source for all 2992 fits; the paper does not validate this benchmark against independent experiments.
  • domain assumption All Stark components share the same FWHM omega = omega_e + omega_i (Eq. 4).
    Used in Eqs. (4)-(5); no test of whether component widths differ.
  • domain assumption The component displacements are strictly linear in k, dk = s k (Eq. 6).
    This linear form is imposed in the model and then 'confirmed' from the fitted slope s, which is circular.
  • domain assumption The reduced fictitious mass parametrization mu_r = mu Te/Tg (Eq. 1) captures ion dynamics effects.
    Used to span the simulation grid; assumes ion mobility is fully described by this single parameter.
  • domain assumption Fine structure effects are negligible for H-alpha at electron densities above 10^20 m^-3.
    Stated in Section 2.1 following reference [13]; used to justify the 15-component model.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exact Stark analytical function for H{\alpha} line based on the FFM Model related with plasma parameters." pith.science (2026). https://pith.science/paper/XNG2VZ2Z

@misc{pith2026250111113,
  author       = {Pith},
  title        = {Pith review of: Exact Stark analytical function for H\alpha line based on the FFM Model related with plasma parameters},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XNG2VZ2Z}},
  note         = {Machine review of arXiv:2501.11113}
}
read the original abstract

Optical Emission Spectroscopy is a widely used technique for plasma diagnosis, with particular interest in hydrogen atomic emission due to its prevalence in plasmas. However, accurately determining plasma parameters like electron density, electron temperature, and gas temperature starting from the experimental profiles remains a challenge. This paper introduces a comprehensive model for Stark broadening of the H{\alpha} line in a wide range of plasma conditions, addressing the limitations of existing analytical expressions for line shapes. The proposed model encompasses the full splitting of the transition into fifteen Lorentzian profiles and electric micro-field fluctuations surrounding the emitting atoms due to collisions with charged particles. Starting from accurate spectral data obtained from realistic computer simulations, fitting parameters of the model, have been obtained by using an optimization method based on a genetic algorithm. The set of parameters of the model are reported for a wide range of plasma conditions. The behavior of these parameters is analyzed to understand their dependence in terms of the electron density and temperature and gas density of the plasma. The model parameters here obtained constitute a useful tool in plasma diagnosis to obtaining the values of the physical parameters of the plasma starting from the experimental profiles.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Philip, B

    N. Philip, B. Saoudi, M. Crevier, M. Moisan, J. Barbeau, J. Pelletier, The respective roles of UV photons and oxygen atoms in plasma sterilization at reduced gas pressure: the case of N2–O2 mixtures, IEEE Trans. Plasma Sci. 30 (2002) 1429–36

    work page 2002

  2. [2]

    Moisan, J

    M. Moisan, J. Barbeau , M.C. Crevier, J. Pelletier , N. Philip , B. Saoudi, Plasma sterilization. Methods and mechanisms, Pure Appl. Chem. 74 (2003) 349–58

    work page 2003

  3. [3]

    Moisan, J

    M. Moisan, J. Barbeau, J. Pelletier, Plasma sterilization - methods and mechanisms, Le vide, Sci. Tech. Appl. 299 (2001) 15–28. 21

    work page 2001

  4. [4]

    Feugeas, B

    J. Feugeas, B. Gomez, A. Craievich, Ion nitriding of stainless steels. Real time surface characterization by synchrotron X-ray diffraction, Surface & Coatings Technology. 154 (2-3) (2002) 167-75

    work page 2002

  5. [5]

    Kabouzi , M

    Y. Kabouzi , M. Moisan , J.C. Rostaing , C. Trassy , C. Guérin , D. Kéoroak , Z. Zakrzewski, Abatement of perfluorinated compounds using microwave plasmas at atmospheric pressure, J. Appl. Phys. 93 (2003) 9483–96

    work page 2003

  6. [6]

    Rostaing, F

    J.C. Rostaing, F. Bryselbout , M. Moisan , J.C. Parent, An efficient method for purifying noble gases using high frequency electric field discharges, C.R. Acad. Sci. Paris Sér. IV Phys. Astrophys. 1 (1) (2001) 99–105

    work page 2001

  7. [7]

    Parigger and E

    C.G. Parigger and E. Oks, Hydrogen Balmer Series Spectroscopy in laser-induced breakdown plasmas, International Review of Atomic and Molecular Physics. 1 (1) (2010) 13-23

    work page 2010

  8. [8]

    Borkowska-Burnecka, W

    J.B. Borkowska-Burnecka, W. Zyrnichi, M.Welna ann P. Jamroz , Electron density from Balmer Series Hydrogen Lines and Ionization Temperatures in Inductively Coupled Argon Plasmas Sup plied by Aerosol and Volatile Species, International Journal of Spectroscopy. Article ID 7521050 (2016) 1-9

    work page 2016

Show all 26 references
  1. [9]

    García, M.S

    Yubero, M.C. García, M.S. Dimitrijevi ć, A. Sola and A. Gamero, Measuring the electron density in plasmas from the difference of Lorentzian part of the widths of two Balmer series hydrogen lines. Spectrochim. Acta Part B. 107 (2015) 164-169

    2015

  2. [10]

    Laux, T.G

    C.O. Laux, T.G. Spence, C.H. Kruger and R.N. Zare, Optical diagnostics of atmospheric pressure air plasmas, Plasma Sources Sci. Technol. 12 (2003) 125-138

    2003

  3. [11]

    Nikiforo v, Ch

    A.Y. Nikiforo v, Ch. Leys, M.A. González and J.L. Walsh, Electron density measurement in atmospheric pressure plasma jets: Stark broadening of hydrogenated and non-hydrogenated lines. Plasma Sources. Sci. Technol. 24 (2005) 034001. 22

    2005

  4. [12]

    Yubero, M.C

    C. Yubero, M.C. Garcia and M.D. Calzada, On the use of the H alpha spectral line to determine the electron density in a microwave (2.45 GHz) plasma torch at atmospheric pressure, Spectrochim. Acta Part B. 61(5) (2006) 540-544

    2006

  5. [13]

    Konjevic, M

    N. Konjevic, M. Ivkovic and N. Sakan, Hydrogen Balmer lines for low electron number density plasma diagnosis. Spectrochim. Acta Part B. 76 (2012) 16-26

    2012

  6. [14]

    Calisti, L

    A. Calisti, L. Godbert, R. Stamm and B. Talin. Fast numerical method for line shape studies in hot and dense plasmas. J.Quant. Spectrosc. Radiat. Transfer. 51 (1994) 59-64

    1994

  7. [15]

    Calisti, C

    A. Calisti, C. Mosse, M. Koubiti, R. Stamm and B. Talin. Study of some plasma radiative properties based on the frequency fluctuation model approach. J.Quant. Spectrosc. Radiat. Transfer. 54 (1995) 89-95

    1995

  8. [16]

    Calisti, F

    A. Calisti, F. Khelfaoui, R. Stamm and B. Talin. Spectral line Shapes. Ed. L. Frommhold, AIP (New York) (1990) volume 6

    1990

  9. [17]

    Griem, A.C

    H.R. Griem, A.C. Kolb and K.Y. Shen , Stark Broadening of Hydrogen Lines in a Plasma. Phys.Rev. 116 (4)

  10. [18]

    Gigosos, Stark broadening models for plasma diagnostics

    M.A. Gigosos, Stark broadening models for plasma diagnostics. J. Phys. D: Appl. Phys. 47 (2014) 343001 (33pp)

    2014

  11. [19]

    Gigosos and V

    M.A. Gigosos and V. Cardeñoso, New plasma diagnosis tables of hydrogen Stark broadening including ion dynamics, J. Phys. B: At. Mol. Opt. Phys. 29 (1996) 4795-828

    1996

  12. [20]

    Gigosos, M.A

    M.A. Gigosos, M.A. González and V. Cardeñoso, Computer simulated Balmer - alpha, -beta and –gamma Stark line profiles for non -equilibrium plasmas diagnostics, Spectrochim. Acta Part B. 58 (2003) 1489-1504

    2003

  13. [21]

    Díaz-Soriano, J.M

    A. Díaz-Soriano, J.M. Alcaraz-Pelegrina, A. Sarsa, M.S. Dimitrijevic, C. Yubero, A simple and accurate analytical model of the Stark profile and its application to plasma characterization, J. Quant. Spectrosc. Radiat. Transfer 207 (2018) 89-94. 23

    2018

  14. [22]

    Díaz-Soriano, M.S

    A. Díaz-Soriano, M.S. Dimitrijevic, J.M. Alcaraz -Pelegrina, A. Sarsa, C. Yubero, Simple and analytical function for the Stark profile of the H α line and its application to plasma characterization, J. Quant. Spectrosc. Radiat. Transfer. 217 (2018) 111-115

    2018

  15. [23]

    Ortiz -Mora, A

    A. Ortiz -Mora, A. Diaz -Soriano, A. Sarsa, M.S. Dimitrijevic and C. Yubero, A practical method for plasma diagnosis with Balmer series hydrogen lines, Spectrochim. Acta Part B.163 (2020) 15728

    2020

  16. [24]

    J. M. Alcaraz-Pelegrina, A. Sarsa, M.S. Dimitrijevic and C. Yubero, Analysis of the ion collisional contribution over the Stark profile in H α line. Spectrochim. Acta Part B: Atomic Spectroscopy. 194 (2022) 106455

    2022

  17. [25]

    Cvetanovic, S.S

    N. Cvetanovic, S.S. Ivkovic, B.M. Obradovic and M. M. Kuraica, Simultaneous influence of Stark effect ad excessive line broadening on the H α line. Eur. Phys. J. D 71 (2017) 317-327

    2017

  18. [26]

    H. Holland . Adaptation in natural and artificial systems. MIT Press (Cambridge, MA) (1992)

    1992

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.