Eliminating imaginary vibrational frequencies in quantum-chemical cluster models of enzymatic active sites

John Herbert,Paige Bowling,Saswata Dasgupta
DOI: https://doi.org/10.26434/chemrxiv-2024-fvvms
2024-02-09
Abstract:In constructing finite models of enzyme active sites for use in quantum-chemical calculations, atoms at the periphery of the model system are often constrained to prevent structural collapse during geometry relaxation. A simple fixed-atom or ``coordinate lock'' approach is commonly employed but leads to undesirable artifacts including the appearance of small imaginary frequencies. These preclude the evaluation of finite-temperature free energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials, and here we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While that approach does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy, in addition to the artificial rigidity already introduced by fixed-atom constraints. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct models that can be used in unconstrained geometry relaxations.
Chemistry
What problem does this paper attempt to address?
This paper addresses the issue of imaginary vibrational frequencies that arise when studying enzyme active sites using a quantum chemical cluster model. These imaginary frequencies hinder the evaluation of finite temperature free energy corrections and limit thermal chemistry calculations to enthalpy-based only. The study compares two methods: fixing atoms and using a harmonic constraint potential to eliminate imaginary frequencies. The fixed atom method leads to an underestimation of zero-point energy and vibrational entropy as well as artificial rigidity of the system, while the harmonic constraint potential eliminates imaginary frequencies and allows for unconstrained geometry optimization.