Simulation of Vibrational Circular Dichroism Spectra Using Second-Order Møller-Plesset Perturbation Theory and Configuration Interaction Doubles

Brendan M. Shumberger,T. Daniel Crawford
2024-06-08
Abstract:We present the first single-reference calculations of the atomic axial tensors (AATs) using wave-function-based methods including dynamic electron correlation effects using second-order Møller-Plesset perturbation theory (MP2) and configuration interaction doubles (CID). Our implementation involves computing the overlap of numerical derivatives of the correlated wave functions with respect to both nuclear displacement coordinates and the external magnetic field. Out test set included three small molecules, including the axially chiral hydrogen molecule dimer and (P)-hydrogen peroxide, and the achiral H2O. For our molecular test set, we observed deviations of the AATs for MP2 and CID from that of the Hartree-Fock (HF) method upwards of 49%, varying with the choice of basis set. For (P)-hydrogen peroxide, electron correlation effects on the VCD rotatory strengths and corresponding spectra were particularly significant, with maximum deviations of the rotatory strengths of 62% and 49% for MP2 and CID, respectively, using our largest basis set. The inclusion of dynamic electron correlation to the computation of the AATs can have a significant impact on the resulting rotatory strengths and VCD spectra.
Chemical Physics,Quantum Physics
What problem does this paper attempt to address?
This paper mainly discusses the use of second-order Møller-Plesset perturbation theory (MP2) and coupled interaction doubles (CID) method to simulate the vibrational circular dichroism (VCD) spectra. The authors implemented, for the first time, the calculation of the nuclear axis tensor (AAT) for VCD spectra based on wavefunction dynamics. They computed the wavefunction derivatives with respect to nuclear coordinates and external magnetic fields using numerical differentiation methods and combined this information to obtain AAT. The paper mentions that the rotational intensity of VCD spectra is zero under the Born-Oppenheimer approximation, which is a theoretical challenge. Previous studies proposed some heuristic models to predict VCD spectra, but there are limitations. Nafie and Freedman introduced nuclear perturbations to overcome this problem, while Stephens proposed a formula based on the ground-state wavefunction to calculate the magnetic dipole vibrational transition moment. The computational part of the paper involves the calculation of AAT rotational intensities for small molecules, including hydrogen dimer, water, and (P)-hydrogen peroxide, using the MP2 and CID methods. The results show that the effect of electron correlation on AAT is significant, with maximum deviations reaching 49%, and the impact on the VCD rotational intensities of (P)-hydrogen peroxide is particularly large, with a maximum change of 62%. Furthermore, the authors compare the results obtained with different basis sets and highlight the importance of dynamic electron correlation effects on the final rotational intensities and VCD spectra. The finite-difference method implemented with Python code MagPy is used to validate and compute these effects. The study shows that as the size of the basis set increases, the differences between the MP2 and CID methods also increase, and CID generally exhibits larger variations. In summary, this paper aims to improve the computational methods for VCD spectra by considering dynamic electron correlation effects, providing more accurate tools for understanding and predicting the stereochemical properties of biomolecules.