15min:

CHARLOTTE E. HINKLE, ANNE B. MCCOY, Department of Chemistry, The Ohio State University, Columbus, OH 43210.

Protonated methane has long been of interest to astrochemists due to its presumed importance as a reaction intermediate in the reaction involving CH₃⁺ + HD within the interstellar medium. Within the interstellar medium there is a nonstatistical H/D isotopic abundance observed for isotopologues of CH₃⁺. While classical trajectory calculations have been performed dissociating CH₅⁺ and CH₄D⁺ into the fragments, CH₃⁺ + H₂, CH₂D⁺ + H₂ and CH₃⁺ + HD, these calculations do not account for a large portion of the available energy being tied up in the zero point energy of the reactants and products. Earlier work in our group on CH₂D₃⁺ showed the deuterium atoms were localized to the CH₃⁺ group, rather than the H₂ moiety. Classical calculations fail to account for this observed localization, instead showing full delocalization of D between both CH₃⁺ and H₂. With a quantum mechanical treatment, the energetics and wave functions will depend on which asymptotic channel is chosen, while in the classical treatment, these channels will be energetically equivalent. By performing Diffusion Monte Carlo simulations in Jacobi coordinates, we can constrain the distance between the CH₃⁺ and H₂ subunits. Using this technique we have evaluated a one-dimensional reaction potential that includes the full anharmonic zero point energy in the remaining degrees of freedom and can determine how energetics of this reaction change upon partial deuteration of CH₃⁺ or H₂. We have also evaluated the probability amplitude associated with the wave functions that are obtained in the DMC simulations at various values of the reaction coordinate.