Search papers, labs, and topics across Lattice.
This paper introduces a novel computational method called NEO-sDFT, which integrates Nuclear-Electronic Orbital and Subsystem Density Functional Theory to facilitate quantum chemical calculations for large molecular systems containing thousands of atoms while treating selected protons quantum mechanically. The method's flexibility allows for various embedding schemes within a unified framework, significantly enhancing its applicability. Accuracy assessments reveal that NEO-sDFT achieves interaction energy errors of only a few meV compared to traditional Nuclear-Electronic Orbital Density Functional Theory, while also demonstrating favorable computational scaling for large systems.
NEO-sDFT achieves quantum mechanical treatment of protons in massive molecular systems with unprecedented accuracy and efficiency, opening new avenues for simulating complex biochemical processes.
We present a new computational approach based on Nuclear-Electronic Orbital and Subsystem Density Functional Theory. The resulting NEO-sDFT methodology enables quantum chemical calculations of large molecular systems composed of thousands of atoms while also describing selected protons quantum mechanically. A key feature of the associated program implementation is its flexibility in assigning quantum protons to subsystems, allowing different embedding schemes to be employed within a unified framework. The accuracy of NEO-sDFT is assessed by computing water dimer interaction energies, which exhibit errors of only a few meV compared to reference Nuclear-Electronic Orbital Density Functional Theory results. The favorable computational scaling of NEO-sDFT is further showcased through calculations on water clusters of increasing size, with up to 303 protons being treated quantum mechanically. Our results demonstrate that NEO-sDFT is an accurate and efficient approach for incorporating nuclear quantum effects into large-scale simulations, providing a foundation for future applications to proton transfer processes in biochemical systems.