Chemistry in laser‑induced plasmas at local thermodynamic equilibrium

The equation of state for plasmas containing negative and positive ions of elements and molecules formed by these elements is modeled under the assumption that all ionization processes and chemical reactions are at local thermal equilibrium and the Coulomb interaction in the plasma is described by the Debye–Huckel theory. The hierarchy problem for constants of molecular reactions is resolved by using three different algorithms for high, medium, and low temperatures: the contraction principle, the Newton–Raphson method, and a scaled Newton–Raphson method, respectively. These algorithms are shown to have overlapping temperature ranges in which they are stable. The latter allows one to use the developed method for calculating the equation of state in combination with numerical solvers of Navier–Stokes equations to simulate laser-induced plasmas initiated in an atmosphere and to study formation of molecules and their ions in such plasmas. The method is applicable to a general chemical network. It is illustrated with examples of Ca–Cl and C–Si–N laser-induced plasmas.