Acceleration of Correlation-corrected Vibrational Self-consistent Field Calculation Times for Large Polyatomic Molecules

Acceleration of the correlation-corrected Vibrational self-consistent field (CC-VSCF) method for anharmonic calculations of vibrational states of polyatomic molecules is described. The acceleration assumes pairwise additive interactions between different normal modes, and employs orthogonality of the single-mode vibrational wavefunctions. This greatly reduces the effort in computing correlation effects between different vibrational modes, which is treated by second order perturbation theory in CC-VSCF. The acceleration can improve the scaling of the overall computational effort from N6 to N4, where N is the number of vibrational modes. Sample calculation times, using semi-empirical potential surfaces (PM3), are given for a series of glycine peptides. Large computational acceleration, and significant reduction of the scaling of the effort with system size, is found and discussed.