A spherical electron cloud hopping model for studying product branching ratios of dissociative recombination

Knowledge-based probabilistic representations of branching ratios in chemical networks: the case of dissociative recombinations

Experimental data about branching ratios for the products of dissociative recombination of polyatomic ions are presently the unique information source available to modelers of natural or laboratory chemical plasmas. Yet, because of limitations in the measurement techniques, data for many ions are incomplete. In particular, the repartition of hydrogen atoms among the fragments of hydrocarbons ions is often not available. A consequence is that proper implementation of dissociative recombination processes in chemical models is difficult, and many models ignore invaluable data. We propose a novel probabilistic approach based on Dirichlet-type distributions, enabling modelers to fully account for the available information. As an application, we consider the production rate of radicals through dissociative recombination in an ionospheric chemistry model of Titan, the largest moon of Saturn. We show how the complete scheme of dissociative recombination products derived with our method dramatically affects these rates in comparison with the simplistic H-loss mechanism implemented by default in all recent models.

Spherical electron cloud hopping molecular dynamics simulation on dissociative recombination of protonated water

Dissociative recombination (DR) of H(3)O(+) with electrons at zero collision energy has been studied by a direct ab initio molecular dynamics method on four low-lying electronic states of the system. Initial conditions for trajectories are determined by a spherical electron cloud hopping (SECH) model, while nonadiabatic effects are considered through a surface hopping scheme. The energies, forces, and nonadiabatic coupling strengths (NACS) used in trajectory propagations are calculated on-the-fly via state-average complete active self-consistent field (CASSCF) theory with full valence electrons. Dynamics results show that the H(3)O(+) DR is ultrafast and yields diversity of products. Product branching fractions are predicted to be 0.660 for (OH + 2H), 0.230 for (H(2)O + H), 0.108 for (OH + H(2)), and 0.002 for (O + H + H(2)), which are in excellent agreement with the heavy-ion storage ring experimental results. Kinetic energies of the eliminated hydrogen atoms are large and show a bimodal distribution.