Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published April 2023 | public
Journal Article

Quantum harmonic free energies for biomolecules and nanomaterials

Abstract

Obtaining the free energy of large molecules from quantum mechanical energy functions is a long-standing challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of large size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocrystals, and show that the error per atom decreases with system size in the limit of large systems. We further show that we can obtain the vibrational contributions to the binding free energies of prototypical protein–ligand complexes where exact computation is too expensive to be practical. Our work raises the possibility of routine quantum mechanical estimates of thermodynamic quantities in complex systems.

Additional Information

© The Author(s), under exclusive licence to Springer Nature America, Inc. 2023. We thank S. Murlidaran for help with protein preparation. A.F.W. (theoretical analysis) was supported by the US Department of Energy (grant no. DE-SC0018140). C.L. (computational studies) was supported by the US National Science Foundation (grant no. 1931328). X.Z. (improved DFT implementation for proteins) was supported by the US Department of Energy (grant no. DE-SC0019330). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. Contributions. A.F.W. and G.K.C. conceived the project. A.F.W., C.L. and X.Z. performed the work. All authors contributed to the writing of the paper. Data availability. Data for all systems can be regenerated by running our stochastic Lanczos code on the geometries available at https://doi.org/10.6084/m9.figshare.22258447.v1 (ref. 35). Source Data are provided with this paper. Code availability. The stochastic Lanczos code is available in Supplementary Software 1. The Gaussian and plane waves implementation is available through PySCF at www.pyscf.org. Competing interests. Garnet K. Chan is a co-owner of QSimulate, Inc. The other authors declare no competing interests.

Additional details

Created:
August 22, 2023
Modified:
October 20, 2023