Diffusion Monte Carlo Approaches for Evaluating Rotationally Excited States of Symmetric Top Molecules: Application to H₃O⁺ and D₃O⁺

Abstract

An approach is described for evaluating energies and wave functions for rotationally excited states of symmetric top molecules using diffusion Monte Carlo methods. The approach is based on the fact that, for many systems, the rotation/vibration Hamiltonian can be modeled by terms that depend on the vibrational coordinates and powers of the components of the rotational angular momentum vector, J. In the case of symmetric top molecules with M = 0, the rotational part of the wave function is given by the tesseral harmonics. We construct rotationally excited states within the diffusion Monte Carlo approach by imposing nodal surfaces that are obtained from the roots of the tesseral harmonics. Results are presented for H(3)O(+) and D(3)O(+) with J <or= 10. Where comparisons to previous calculations can be made, the agreement is excellent.

Guided Diffusion Monte Carlo: A Method for Studying Molecules and Ions That Display Large Amplitude Vibrational Motions

“…For hydronium, both the guided and unguided approaches yield results that are in good agreement with the previously reported DMC ground state energy of 7453 cm −1 , which was obtained using an unguided simulation with 20 000 walkers and a 1 au time step. 38 For the unguided calculations, ensembles larger than 10 000 walkers yield energies that are well-converged, while for the guided calculations, 1000 walkers are needed to achieve similar accuracy. These results are consistent with our earlier studies of water clusters.…”

GPU-Accelerated Neural Network Potential Energy Surfaces for Diffusion Monte Carlo

“…The zero-point energy is evaluated by taking the time-averaged value of V ref (τ) in eq . In the expression for V ref , V̅ is the average potential energy of the ensemble, α is 0.5/Δτ, , and N w (τ) is the number of walkers at imaginary time τ. Once equilibrated, the density of the walkers near the molecular configuration represented by the coordinates, x , is proportional to the value of the ground state wave function at that geometry.…”

Rotation/Torsion Coupling in H5+, D5+, H4D+, and HD4+ Using Diffusion Monte Carlo

“…A more complete description of the theory and our implementation can be found in earlier publications. 18,22,24,25,27 Briefly, the algorithm that we use is based on the studies of Anderson. 23, 27 The approach derives from the parallel structure of the diffusion equation, which includes a coordinate-dependent rate process…”