Analytic first and second derivatives of the energy are developed for the fragment molecular orbital method interfaced with molecular mechanics in the electrostatic embedding scheme at the level of Hartree‐Fock and density functional theory. The importance of the orbital response terms is demonstrated. The role of electrostatic embedding upon molecular vibrations is analyzed, comparing force field and quantum mechanical treatments for an ionic liquid and a solvated protein. The method is applied for 100 protein conformations sampled in molecular dynamics (MD) to take into account the complexity of a flexible protein structure in solution, and a good agreement with experimental data is obtained: Frequencies from an experimental infrared (IR) spectrum are reproduced within 17 cm−1.
Financed by the National Centre for Research and Development under grant No. SP/I/1/77065/10 by the strategic scientific research and experimental development program:
SYNAT - “Interdisciplinary System for Interactive Scientific and Scientific-Technical Information”.