Analytic second derivatives for the efficient electrostatic embedding in the fragment molecular orbital method.

Area of Science:

• Computational Chemistry
• Theoretical Chemistry
• Spectroscopy

Background:

• Accurate calculation of molecular properties is crucial for understanding chemical processes.
• Fragment Molecular Orbital (FMO) methods offer a way to study large molecular systems.
• Simulating vibrational spectra provides insights into molecular structure and dynamics.

Purpose of the Study:

• To develop and validate an efficient computational method for calculating analytic second energy derivatives.
• To apply the method to simulate infrared spectra of proteins.
• To analyze the nature of localized vibrations within protein structures.

Main Methods:

• Developed analytic second derivatives for restricted Hartree-Fock and density functional theory using the two-body FMO method.
• Incorporated an electrostatic embedding potential with self-consistent determination using point charges and electron densities.
• Validated accuracy against exact embedding potentials and full calculations; assessed computational efficiency via parallelization.

Main Results:

• The developed FMO-based method accurately reproduces vibrational spectra.
• Simulations of Trp-cage and crambin proteins successfully reproduced infrared spectra.
• Analysis revealed localized vibrational modes within the Amide I peak of crambin and Tyr symmetric stretch in Trp-cage.

Conclusions:

• The FMO method with electrostatic embedding provides an accurate and computationally efficient approach for vibrational spectroscopy of large molecules.
• This method enables detailed analysis of localized vibrations in complex biological systems.
• The findings facilitate a deeper understanding of protein dynamics and function through vibrational analysis.