Search

Rovibrational energies, wave functions, and Raman transition moments are reported for the lowest-energy states of the H$_3^+$ molecular ion including the magnetic couplings of the proton spins and molecular rotation in the presence of a weak external magnetic field. The rovibrational-hyperfine-Zeeman Hamiltonian matrix is constructed and diagonalized using the rovibrational eigenstates and the proton spin functions. The developed methodology can be used to compute hyperfine-Zeeman effects also for higher-energy rovibrational excitations of H$_3^+$ and other polyatomic molecules. These developments will guide future experiments extending quantum logic spectroscopy to polyatomic systems.

Full-dimensional (12D) vibrational states of the methanol molecule (CH$_3$OH) have been computed using the GENIUSH-Smolyak approach and the potential energy surface from Qu and Bowman (2013). All vibrational energies are converged better than 0.5 cm$^{-1}$ with respect to the basis and grid size up to the first overtone of the CO stretch, ca. 2000 cm$^{-1}$ beyond the zero-point vibrational energy. About seventy torsion-vibration states are reported and assigned. The computed vibrational energies agree with the available experimental data within less than a few cm$^{-1}$ in most cases, which confirms the good accuracy of the potential energy surface. The computations are carried out using curvilinear normal coordinates with the option of path-following coefficients which minimize the coupling of the small- and large-amplitude motions. It is important to ensure tight numerical fulfilment of the $C_{3\mathrm{v}}$(M) molecular symmetry for every geometry and coefficient set used to define the curvilinear normal coordinates along the torsional coordinate to obtain a faithful description of degeneracy in this floppy system. The reported values may provide a computational reference for fundamental spectroscopy, astrochemistry, and for the search of the proton-to-electron mass ratio variation using the methanol molecule.

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible $\Delta K\neq 0$ transition moments arise in this symmetric top system due to strong rovibrational couplings.

All vibrational energies of the (trans-, cis-, delocalized-) formic acid molecule are converged up to 4500 cm$^{-1}$ beyond the zero-point vibrational energy with the GENIUSH-Smolyak variational approach and using an ab initio potential energy surface [D. P. Tew and W. Mizukami, J. Phys. Chem. A, 120, 9815-9828 (2016)]. Full-dimensional dipole and polarizability surfaces are fitted to points computed at the CCSD/aug-cc-pVTZ level of theory. Then, body-fixed vibrational dipole and polarizability transition moments are evaluated and used to simulate jet-cooled infrared and Raman spectra of HCOOH. The benchmark-quality vibrational energy, transition moment, and wave function list will be used in further work in comparison with vibrational experiments, and in further rovibrational computations.

Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schr\"odinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are `delocalized' over multiple potential energy wells. These `large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schr\"odinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.

The automated development of a new ab initio full-dimensional potential energy surface (PES) is reported for the CH$_4\cdot$F$^-$ complex using the ROBOSURFER program package. The new potential provides a near-spectroscopic quality description over a broad configuration range including the methane-ion dissociation, as well as isolated methane vibrations. In particular, it improves upon the earlier [Czak\'o, Braams, Bowman (2008)] PES over intermediate methane-fluoride distances. Full-dimensional (12D) variational vibrational computations using the new PES and the GENIUSH-Smolyak algorithm show that tunneling splittings larger than 0.1 cm$^{-1}$ appear below the top of the interconversion barrier of the four equivalent minima of the complex.

Vibrational states of the formic acid molecule are converged using the GENIUSH-Smolyak approach and the potential energy surface taken from [D. Tew and W. Mizukami, J. Phys. Chem. A 120, 9815 (2016)]. The quantum nuclear motion is described by using the $cis$-$trans$ torsional coordinate and eight curvilinear normal coordinates defined with respect to an instantaneous reference configuration changing as a function of the torsional degree of freedom. Harmonic oscillator basis functions are used for the curvilinear normal coordinates, a Fourier basis for the torsional coordinate, and a simple basis pruning condition is combined with a Smolyak integration grid. $Trans$, $cis$, and $delocalized$ vibrational states are reported up to and slightly beyond the isomerization barrier.