Fixed-Node, Importance-Sampling Diffusion Monte Carlo for Vibrational Structure with Accurate and Compact Trial States

Vibrational second-order perturbation theory based on curvilinear coordinates: Thermochemical applications

This work improves and extends a general and robust workflow for the computation of anharmonic vibrational frequencies to thermodynamic functions, paving the way toward the study of large flexible molecules. The key new feature is the extension of closed-form expressions for both zero-point vibrational energies and partition functions to second-order vibrational perturbation theory based on curvilinear internal coordinates. The use of curvilinear coordinates enables the reduction of couplings between different degrees of freedom, enriching the arsenal of existing vibrational approaches, and can lead to effective, low-dimensional linear-scaling models. The accuracy of the results obtained for some prototypical systems paves the way toward the systematic use of this new implementation in the study of molecules containing a few dozen atoms, as exemplified by the test cases of a molecular motor, a nucleoside, and two hormones.

Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order

The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.