Vibrational predissociation of H5+

Probing the Relationship Between Large-Amplitude Motions in H5(+) and Proton Exchange Between H3(+) and H2

Understanding the spectroscopy and dynamics of H5(+) is central in gaining insights into the H3(+) + H2 → H5(+) → H2 + H3(+) proton transfer reaction. This molecular ion exhibits large-amplitude vibrations, which allow for the transfer of a proton between H3(+) and H2 even in its ground vibrational state. With vibrational excitation, the number of open channels for permutations of protons increase. In this work, the minimized energy path variant of diffusion Monte Carlo is used to investigate how the energetically accessible proton permutations evolve as H5(+) is dissociated into H3(+) + H2. Two mechanisms for proton permutation are investigated. The first is the proton hop, which correlates to large-amplitude vibrations of the central proton in H5(+). The second is the exchange of a pair of hydrogen atoms between H3(+) and H2. This mechanism requires several proton hops along with a 120° rotation of H3(+) within the H5(+) molecular ion. This analysis shows that while there is a narrow region of configuration space over which both isomerization processes are energetically accessible, full permutation of the five protons in H5(+) more likely occurs through a stepwise mechanism. Such full permutation of the protons becomes accessible when the shared proton stretch is excited to the vpt = 2 or 3 excited state. The effects of deuteration and rotational excitation of the H2 and H3(+) products are also investigated. Deuteration inhibits permutation of protons, while rotational excitation has only a small impact on these processes.

Vibrational dynamics of the H5(+) and its isotopologues from multiconfiguration time-dependent Hartree calculations

Full-dimensional multiconfiguration time-dependent Hartree (MCTDH) computations are reported for the vibrational states of the H(5)(+) and its H(4)D(+), H(3)D(2)(+), H(2)D(3)(+), HD(4)(+), D(5)(+) isotopologues employing two recent analytical potential energy surfaces of Xie et al. [J. Chem. Phys. 122, 224307 (2005)] and Aguado et al. [J. Chem. Phys. 133, 024306 (2010)]. The potential energy operators are constructed using the n-mode representation adapted to a four-combined mode cluster expansion, including up to seven-dimensional grids, chosen adequately to take advantage in representing the MCTDH wavefunction. An error analysis is performed to quantify the convergence of the potential expansion to reproduce the reference surfaces at the energies of interest. An extensive analysis of the vibrational ground state properties of these isotopes and comparison with the reference diffusion Monte Carlo results by Acioli et al. [J. Chem. Phys. 128, 104318 (2008)] are presented. It is found that these systems are highly delocalized, interconverting between equivalent minima through rotation and internal proton transfer motions even at their vibrational ground state. Isotopic substitution affects the zero-point energy and structure, showing preference in the arrangements of the H and D within the mixed clusters, and the most stable conformers of each isotopomer are the ones with the H in the central position. Vibrational excited states are also computed and by comparing the energies and structures predicted from the two surfaces, the effect of the potential topology on them is discussed.

IR photodissociation spectroscopy of H7(+), H9(+), and their deuterated analogues

Cluster ions of H7(+)/D7(+) and H9(+)/D9(+) produced in a supersonic molecular beam with a pulsed discharge source are mass selected and studied with infrared laser photodissociation spectroscopy. Photodissociation occurs by the loss of H2 (D2) from each cluster, producing resonances in the 2000-4500 cm(-1) region. Vibrational patterns indicate that these ions consist of an H3(+) (D3(+)) core ion solvated by H2 (D2) molecules. There is no evidence for the shared proton structure seen previously for H5(+). The H3(+) ion core vibrational bands are weakened and broadened significantly, presumably by enhanced rates of intramolecular vibrational relaxation. Computational studies at the DFT/B3LYP or MP2 levels of theory (including scaling) are adequate to reproduce qualitative details of the vibrational spectra, but neither provides quantitative agreement with vibrational frequencies.

Mid- and Far-IR Spectra of H5(+) and D5(+) Compared to the Predictions of Anharmonic Theory

H5(+) is the smallest proton-bound dimer. As such, its potential energy surface and spectroscopy are highly complex, with extreme anharmonicity and vibrational state mixing; this system provides an important benchmark for modern theoretical methods. Unfortunately, previous measurements covered only the higher-frequency region of the infrared spectrum. Here, spectra for H5(+) and D5(+) are extended to the mid- and far-IR, where the fundamental of the proton stretch and its combinations with other low-frequency vibrations are expected. Ions in a supersonic molecular beam are mass-selected and studied with multiple-photon dissociation spectroscopy using the FELIX free electron laser. A transition at 379 cm(-1) is assigned tentatively to the fundamental of the proton stretch of H5(+), and bands throughout the 300-2200 cm(-1) region are assigned to combinations of this mode with bending and torsional vibrations. Coupled vibrational calculations, using ab initio potential and dipole moment surfaces, account for the highly anharmonic nature of these complexes.

Toward a realistic density functional theory potential energy surface for the H5+ cluster

The potential energy surface of H(5)(+) is characterized using density functional theory. The hypersurface is evaluated at selected configurations employing different functionals, and compared with results obtained from ab initio CCSD(T) calculations. The lowest ten stationary points (minima and saddle-points) on the surface are located, and the features of the short-, intermediate-, and long-range intermolecular interactions are also investigated. A detailed analysis of the surface's topology, and comparisons with extensive CCSD(T) results, as well as a recent ab initio analytical surface, shows that density functional theory calculations using the B3(H) functional represent very well all aspects studied on the H(5)(+) potential. These include the tiny energy difference between the minimum at 1-C(2v) configuration and the 2-D(2d) one corresponding to the transition state for the proton transfer between the two equivalent C(2v) minima, and also the correct asymptotic behavior of the long-range interactions. The calculated binding energy and dissociation enthalpies compare very well with previous benchmark coupled-cluster ab initio data, and with experimental data available. Based on these results the use of such approach to perform first-principles molecular dynamics simulations could provide reliable information regarding the dynamics of protonated hydrogen clusters.

Submillimetre-wave lines of H2D+ and D2H+ as probes into chemistry in cold dark clouds

We discuss past and recent progress of our continuing project of submillimetre-wave spectroscopic investigations of H2D+ and D2H+. Three new lines of H2D+ in the 2.5-3.5 THz range are measured with a tunable far-infrared laser system. Since these molecules are very light asymmetric molecules, analysis based on a conventional effective Hamiltonian is not very useful in predicting the transition frequencies to the accuracy of the order of several MHz or better. In this respect, any addition of new accurate measurements of transition frequencies is important. In this paper, some discussions will be made on H5+ and its deuterated species as probable interstellar species in cold dark clouds. In particular, D3+, which is predicted to be abundant in cold dark clouds, can be (indirectly) detected by observing D3+ x H2.