Gaussian functions in Hylleraas‐configuration interaction calculations. V. An accurate ab initio H+3 potential‐energy surface

Ground state of the H3(+) molecular ion: physics behind

Five physics mechanisms of interaction leading to the binding of the H3(+) molecular ion are identified. They are realized in a form of variational trial functions, and their respective total energies are calculated. Each of them provides subsequently the most accurate approximation for the Born–Oppenheimer (BO) ground state energy among (two–three–seven)-parametric trial functions being, correspondingly, H2-molecule plus proton (two variational parameters), H2(+)-ion plus H-atom (three variational parameters), and generalized Guillemin–Zener (seven variational parameters). These trial functions are chosen following a criterion of physical adequacy. They include the electronic correlation in the exponential form, exp(γr12), where γ is a variational parameter. Superpositions of two different mechanisms of binding are investigated, and a particular one, which is a generalized Guillemin–Zener plus H2-molecule plus proton (ten variational parameters), provides the total energy at the equilibrium of E = −1.3432 au. The superposition of three mechanisms, generalized Guillemin–Zener plus (H2-molecule plus proton) plus (H2(+)-ion plus H) (14 parameters) leads to the total energy, which deviates from the best known BO energy to 0.0004 au, it reproduces two–three significant digits in exact, non-BO total energy. In general, our variational energy agrees in two–three–four significant digits with the most accurate results available at present as well as major expectation values.

Classical dynamics of laser-driven D₃⁺

A classical model of the triatomic D₃⁺ molecule subjected to an intense, few-cycle laser pulse is introduced. The model is capable of describing the laser-induced correlated motion of both electrons and nuclei in three dimensions, and allows us to follow the motion of the two electrons and three deuterons from the initial field-free state, during the pulse, and until the bond breaking into the final fragments. By averaging over many trajectories, we calculate the relative yields of the ionization and dissociation channels, as well as the kinetic energy release (KER) from the fragment ions. A comparison with recent experimental KER spectra shows good qualitative agreement. In addition, we find a pathway in which an emitted electron recombines into a high-lying Rydberg state, resulting in D + D⁺ + D⁺ fragments with the same KER as in the D⁺ + D⁺ + D⁺ channel.

Explicitly correlated potential energy surface of H3+, including relativistic and adiabatic corrections

After a short historical account of the theory of the H3+ ion, two ab initio methods are reviewed that allow the computation of the ground-state potential energy surface (PES) of H3+ in the Born-Oppenheimer (BO) approximation, with microhartree or even sub-microhartree accuracy, namely the R12 method and the method of explicitly correlated Gaussians. The BO-PES is improved by the inclusion of relativistic effects and adiabatic corrections. It is discussed how non-adiabatic effects on rotation and vibration can be simulated by corrections to the moving nuclear masses. The importance of the appropriate analytic fit to the computed points of the PES for the subsequent computation of the rovibronic spectrum is addressed. Some recent extensions of the computed PES in the energy region above the barrier to linearity are reviewed. This involves a large set of input geometries and the correct treatment of the dissociation asymptotics, including the coupling with the first excited singlet state. Some comments on this state as well as on the lowest triplet state of H3+ are made. The paper ends with a few remarks on the ion H5+.