State-to-State Kinetics of Molecular and Atomic Hydrogen Plasmas MARIO CAPITELLI Department of Chemistry, University of Bari (Italy) CNR Institute of Inorganic Methodologies and Plasmas Bari (Italy)
MOLECULAR DYNAMICS KINETIC MODELS state-resolved elementary process probability non-equilibrium distributions on internal degrees of freeedom FLUID DYNAMICS
MOLECULAR DYNAMICS electron-impact induced processes atom-molecule collision processes gas-surface interaction processes KINETIC MODELS state-resolved elementary process probability non-equilibrium distributions on internal degrees of freeedom FLUID DYNAMICS
ELECTRON-IMPACT INDUCED PROCESSES PROCESSES vibronic excitation direct dissociation THEORETICAL APPROACHES semiclassical impact parameter method Adamson et al., Chem. Phys. Lett., 436 (2007) similarity function approach INPUT DATA molecular potential of involved electronic states transition moment Born high-energy cross section value DYNAMICAL INFORMATION energy and vibrational dependence of cross sections (complete set) isotopic effect investigation state-resolved and macroscopic rate coefficients radiative-decay vibrational excitation cross sections (EV) FEATURES satisfactory accuracy (agreement with experimental data within the error of other methods) low computational load (except in the ab-initio step) suitable only for dipole-allowed transitions no treatment of resonances
INTERFACING SCHEME of DIFFERENT CODES AB-INITIO CALCULATION potentials, transition dipole moments and GOS (Multi Reference Configuration Interaction) GAMESS (General Atomic and Molecular Electronic Structure System ) DYNAMIC CALCULATION solution of the radial Schroedinger equation for quantum vibrational levels and estimation of structural and dynamical factors (evaluation of integrals) IMPACT Cross Section maxwellian eedf distribution Rate Coefficients FAUST BORN Cross Section Einstein coeffs of spontaneous emission Radiative-Decay Vibrational Excitation
INTERFACING SCHEME of DIFFERENT CODES AB-INITIO CALCULATION GAMESS (General Atomic and Molecular Electronic Structure System ) DYNAMIC CALCULATION rapid estimation of cross sections from oscillator strength and transition energy FAUST BORN Cross Section SIMILARITY APPROACH Cross Section IMPACT Cross Section maxwellian eedf distribution Rate Coefficients
electron-impact induced VIBRONIC EXCITATION CROSS SECTIONS for H 2 comparison of different approaches X 1 Σ g (ν i =0) B 1 Σ u X 1 Σ g (ν i =0) C 1 Π u Impact Parameter Method Celiberto 2001 Impact Parameter Method Celiberto 2001 recommended values Itikawa 2008 recommended values Itikawa 2008 similarity approach similarity approach D Ammando, in preparation Itikawa et al., J. Phys. Chem. Ref. Data, 37 (2008) Celiberto et al., ADNDT, 77 (2001)
VIBRATIONALLY-RESOLVED CROSS SECTIONS for H 2 X 1 Σ g (ν i ) B 1 Σ u Impact Parameter Method Celiberto 2001 similarity approach
VIBRATIONALLY-RESOLVED CROSS SECTIONS for H 2 X 1 Σ g (ν i ) C 1 Π u Impact Parameter Method Celiberto 2001 similarity approach
RESONANT PROCESSES in e- H 2 collisions PROCESSES dissociative attachment resonant vibrational excitation (ev processes) THEORETICAL APPROACH local resonance theory INPUT DATA molecular potential of involved electronic states resonance characterization (width and energy dependence) DYNAMICAL INFORMATION angular, energy and vibrational dependence of cross sections (complete set) isotopic effect investigation FEATURES excellent agreement with experimental data resonable computational load
ATOM-MOLECULE COLLISIONS: from PES to complete sets of rovibrational CROSS SECTIONS and RATES PES from literature LIVX Calculation of complete rovibrational ladders relative to all the possible reactants/products, including quasibound states by WKB method DINX Quasiclassical dynamics of atom-diatom collision processes. Distribution of computations on a wide computational grid with periodic check of results, automatic re-calculations for uncompleted jobs and collection of good ones
QCTX Trajectory analysis and cross section results for reactive and non-reactive processes as well as dissociation. The quasibound states are here considered as bound states. Various reactant and product weight functions can be used to analyze the trajectory results Recombination can be obtained in this module by two methods: detailed balance applied to dissociation rates and orbiting resonance applied to rovibrational state-to-state rates XRATE From translational energy dependent rovibrational cross sections to rate coe fficient. This module includes the possibility of selecting and mixing rate coe fficients on the base of lifetime of initial/final states and of the level of detail requested (e.g. total thermal rate, sum of rates on final rotation, dissociation summed to state-to-state rate coefficients with low lifetime, etc.)
TERMOLECULAR HYDROGEN RECOMBINATION Two possible mechanisms: orbiting resonance theory and direct three body recombination (by detailed balance of dissociation data) the two mechanisms can be used in a complementary way, by choosing accurately the states involved in the calculations a detailed database of rate coefficients including initial and final rotation and state selected dissociation for H+H 2 collision process is necessary
TERMOLECULAR HYDROGEN RECOMBINATION normalized recombination vs final vibrational quantum number recombination rates vs final rotational quantum number Esposito&Capitelli J. Phys. Chem. A 133 (2009)
ION-MOLECULE COLLISION PROCESSES PES used: L.P. Viegas, A. Alijah & A.J.C. Varandas, J. Chem. Phys. 126 (2007) couplings calculated from the potential matrix, with particular care for sign consistency of eigenvectors Esposito, in preparation a vibrationally detailed database of rate coefficients should be calculated, using trajectory surface hopping, including initial/final rotation and dissociation Esposito, in preparation
DECONVOLUTION of temperature-dependent RATE COEFFICIENTS to energy-dependent CROSS SECTIONS This method is useful when, for a given chemical process, there are not cross sections data, but there are information about rate coefficients in a given range of temperatures (Minelli, to be published). The only data required are a good functional form for cross section and a temperature dependent rate. The reference rate coefficients are input to a standard nonlinear optimization algorithm (Downhill Simplex Method, DSM (Nelder and Mead, Computer Journal 7 (1965)) to give reconstructed cross sections that satisfy a convergence criterion. Experimental or Theoretical Rate Coefficients deconvolution Cross Sections Case study to test method dissociation of the hydrogen molecule by atom impact comparison between cross sections obtained by DSM (Downhill Simplex Method - full lines) and by QCT (Quasi Classical Trajectory - symbols)
MOLECULAR DYNAMICS for GAS/SURFACE INTERACTION PROCESSES
ISOTOPIC EFFECT: H/D RECOMBINATION on GRAPHITE
PLASMA NUMERICAL EXPERIMENTS ion sources plasma-wall transition electrostatic probe laser-induced plasma hypersonic plasma flow Electric thrusters Negative ion sources Plasma shock tube Divertor region Sheath physics: - see - instability - electronegativity Ion collection Electron collection Dusty LIBS Under-water cavitation bubble dynamics MHD bow shock Space charge region kinetic models (Vlasov, PIC-MCC) fluid models (CFD,MHD,SPH)
ATOMIC HYDROGEN under SHOCK WAVE 50000 x=0 shock front Gas Temperature (K) 40000 30000 20000 upstream T = 5000 K 1 P = 10-2 atm 1 u = 5 1 downstream T (x=0) > T 2 1 P (x=0)> P 2 1 u (x=0)< u 2 1 Euler Equations 10000 x = 1.6 m x = 3.1 m 0-3 -2-1 0 1 2 3 4 5 10 1 10 0 spatial coordinate x (m) thick! = 0 H thin! = 1 CR Master & Electron Boltzmann Equations Radiative Transfer Equation Molar Fractions 10-1 10-2 10-3 H + - e - H 2 kinetics 10-4 -3-2 -1 0 1 2 3 4 5 spatial coordinate x(m)
x(m) x = 1.6 = 1.6 m x(m) x = 1.6 = 1.6 m eedf (ev -3/2 ) 10-1 10-3 10-5 10-7 10-9 10-11 10-13 10-15!=0 N = 1.30 10 15 cm -3 e N = 3.32 10 17 cm -3 H N = 2.02 10 11 cm -3 n=2! = 1 N = 1.20 10 13 cm -3 e N = 5.44 10 16 cm -3 H N = 3.32 10 10 cm -3 n=2 N i /(N H g i ) 10-1 10-3 10-5 10-7 10-9 10-11! = 1 N = 1.20 10 13 cm -3 e N = 5.44 10 16 cm -3 H N = 3.32 10 10 cm -3 n=2!=0 N = 1.30 10 15 cm -3 e N = 3.32 10 17 cm -3 H N = 2.02 10 11 cm -3 n=2 10-17 0 5 10 15 20 25 30 electron energy (ev) 10-13 0 2 4 6 8 10 12 14 level energy (ev) eedf (ev -3/2 ) 10 0 10-2 10-4 10-6 10-8 10-10 10-12 10-14 10-16 x(m) x = 3.1 = 3.1 m! = 1 N = 1.88 10 15 cm -3 e N = 2.54 10 17 cm -3 H N = 7.52 10 10 cm -3 n=2! = 0 N = 7.64 10 14 cm -3 e N = 3.43 10 17 cm -3 H N = 9.99 10 10 cm -3 n=2 N i /(N H g i ) 10 0 10-2 10-4 10-6 10-8 x(m) x = 3.1 = 3.1 m! = 1 N = 1.88 10 15 cm -3 e N = 2.54 10 17 cm -3 H N = 7.52 10 10 cm -3 n=2! = 0 N = 7.64 10 14 cm -3 e N = 3.43 10 17 cm -3 H N = 9.99 10 10 cm -3 n=2 10-18 0 5 10 15 20 25 30 10-10 0 2 4 6 8 10 12 14 electron energy (ev) level energy (ev)
LOCAL SOURCE FUNCTION optically thin plasma (λ=1) B(T) and j e! /k! [J m-2 s -1 sr -1 Hz -1 ] 10-6 10-7 10-8 10-9 10-10 upstream condition P 1 =10-2 atm T 1 =1000K u=20 (5) x = 3.28m Tgas = 12775 K T=12775 K Temperature (K) S n =j e n /k n B n (T) 0 2 4 6 8 10 12 14 16 100000 10000 1000 B(T) and j e! /k! [J m-2 s -1 sr -1 Hz -1 ] T GAS (K) T H [0-1] -2-1 0 1 2 3 4 5 10-5 10-6 10-7 10-8 10-9 (4) x [m] x = 3.29m Tgas = 17817 K x=0 shock front Tgas (K) Telectrons (K) TH[1-2] (K) T=17817 K 0 2 4 6 8 10 12 14 16 T e- S! =j e! /"! B! (T)!! B(T) and j e! /k! [J m-2 s -1 sr -1 Hz -1 ] 10-6 10-7 10-8 10-9 10-10 Source function is the ratio of emissivity over absorption coefficient S " = j " /# " at equilibrium equals the blackbody intensity S " = B " (T ) non-equilibrium distributions of spectral quantities (6) x = 3.30m Tgas = 11461 K T=11461 K S n =j e n /k n B n (T) 0 2 4 6 8 10 12 14 16 E [ev] Capitelli et al., J. Phys. B, 43 (2010) E [ev] E [ev]
MOLECULAR DYNAMICS electron-impact induced processes atom-molecule collision processes gas-surface interaction processes collision integrals THERMODYNAMICS TRANSPORT THEORY transport properties KINETIC MODELS FLUID DYNAMICS state-resolved elementary process probability non-equilibrium distributions on internal degrees of freeedom
THERMODYNAMIC PROPERTIES THERMODYNAMIC PROPERTIES of multicomponent plasmas including Debye-Hückel correction SIMPLIFIED TWO and THREE-LEVEL PARTITION FUNCTION ATOMIC HYDROGEN LEVELS in a BOX in the presence of Coulomb and Debye-Hückel potentials
SIMPLIFIED TWO-LEVEL PARTITION FUNCTION the method is based on the idea of combining a multitude of atomic energy levels into two or three grouped levels. The partition function for atomic hydrogen can be approximated as the method has been applied to many atomic and ionic systems n max =10 n max =100 D Ammando et al., Spectrochim. Acta B, 65 (2010)
ENERGY LEVELS for ATOMIC HYDROGEN in a SPHERICAL BOX box radius = 10 3 a o particle in a spherical box Bohr atom Coulomb potential screened Coulomb potential (Debye length = 10 2 a o ) Capitelli&Giordano, Physical Review A, 80 (2009) n
TRANSPORT PROPERTIES DATABASE OF TRANSPORT CROSS SECTION and COLLISION INTEGRALS for He-H 2 interactions in a wide temperature range [100-50000 K] Bruno et al., Physics of Plasmas (2010) in press EFFECT OF MAGNETIC FIELD (tensorial reformulation) EFFECT OF ELECTRONICALLY EXCITED STATES (EES)
INTERFACING SCHEME of DIFFERENT CODES DYNAMICAL CALCULATION classical collision integrals for elastic collision between species from (accurate, model or phenomenological) interaction potentials OMEGA calculation of the speciesconcentration as a function of temperature HIERARCHICAL CODE TRANSPORT PROPERTY CALCULATION solution of a system of linear equations algebraic equations of order n*k, n being the number of species and K the order of approximation of the expansion in Sonine polynomials CHAPMAN-ESKOG CODE heat flux & diffusion
INTERNAL CONDUCTIVITY for HYDROGEN PLASMA including EES ground-state collision integrals including EES collision integrals Bruno et al., Physics of Plasmas 14 (2007)