【文章內(nèi)容簡介】 ******///cc F ion impact collisions if (nflag_col_spec(3).) then call mc_ion_neut(3,0,dif_engy_fm,n_reac_fm, * num_col_ij_4,extra_col_ij_4,rate_col_ij_4) end ifcc********************************************************///cc CF3 ion impact collisions if (nflag_col_spec(4).) then call mc_ion_neut(4,0,dif_engy_cf3m,n_reac_cf3m, * num_col_ij_5,extra_col_ij_5,rate_col_ij_5) end ifcc********************************************************/// end if end ifcc********************************************************/// return endcc*********************************************************///cc end of MCMCElectron_neutral collisions1）Elastic collisione+Ar→e+Are+CF4→e+CF4e+ N2→e+ N22)Inelastic collisiona) Excitationi. The mian reactionTotalelectronic excitation e+Ar→e+Ar* Vibrational excitation e+CF4→e+CF4 Vibrational excitation e+CF4→e+CF4 Vibrational excitation e+CF4→e+CF4 Electronic excitation e+CF4→e+CF4* Ratational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 Vibrational excitation e+ N2→e+N2 e+ N2→e+N2 * ii. The main formula () ()Where Eth is the threshold energy of the inelastic collision, excitation in this case . Taking into consideration that M+m≈M and g≈v we obtain ( ) ()Where E=mv2/2 is the electron energy before collision. The excitation process is treated as if it were an elastic collision with pre_collision velocities and V. The post_collision velocities are given by Eqs.()() in which all v’s are replaced by . b) Ionization(electron_A)Where A denotes the neutral particle with a mass , which can be assumed equal to that of the ion M .The process is represented as ()Where A+ denotes an ion,e1 is the incident electron , e2 is the ejected electron , and the symbols in the parentheses denote the velocities .The energy balance equation is ()Where Eth is the threshold energy of the ionization . Because of the large ion_to_electron mass ratio ,we can assume that the momentum of the incident electron is much lower than the momentum of the neutral particle,. the incident electron removes an electron from the neutral ,and the neutral bees an ion,continuing on its trajectron undisturbed . On one hand ,this assumption means that the created ion takes the velocity and direction of the neutral particle before collision ,. V’=V .On the other hand , the energy balance can be rewritten ()The left_hand side is known as the excess energyΔE after ionization and we need to find an algorithm how to divide it into the scattered and ejected electrons. In general, when there is no published work on the division of the excess energy , it is divided randomly into two () ()For the electron_Ar ionzation ,however,we use the expression for sampling the Eej by a random number R,()Where Einc =mv2/2 is the energy of the incident electron. The units of α0 andα1 are electronvolts .Once the excess energy is divided, v’ is calculated from Eq.() () ()The post_collision velocity of the scattered electron v’ is obtained by Eq() in which v is replaced by .The pre_collision velocity for the ejected electron , which does not exist in reality, is assumed ()and the post_collision velocity is obtained again from Eq.() by replacing v with and v’ with v’ej . For many types of collision no data on differential cross_section is assumed not to depend to the deflection angle x .The total cross_section calculated from Eq.()gives σT=4πσ（g）.Therefore , the probabilityDefined by ()bees sinxdxdφ/ means that the post_collision velocities take random direction. The scattering with this property is called there are no sufficient data for the differential cross_sections of electron_CF4 and electron_N2 collisions we assume that the scattering is isotropic .The approximation of the interaction potential with the screened Coulomb potential made for collisions with atomic Ar is not valid for molecular case of isotropic scattering the deflection angle , which is in the interval [0, π],is randomly sampled by () The angle φ is randomly sampled from () and the post_collision velocities for electron_CF4 and electron_N2 elastic collisions can be found from Eqs.()_().Similar to electron_Ar elastic collisions,these equations can be somewhat simplified by use of M+m≈M and g≈v. Electron_CF4 and electron_N2 inelastic collisions are treated similar to electron_Ar inelastic collisions,taking into account the isotropic scattering after case of ionization,the excess energy is divided by a ra