This thesis focuses on simulating large molecular systems within and beyond the Born-Oppenheimer framework from first principles. Two approaches have been developed for very different but important applications.
The first one is a hybrid method based on classical force fields that predicts the high-energy ensemble of three-dimensional structures of a class of proteins critical in human physiology: the G protein-coupled receptors (GPCRs). GPCRs' functions rely on their activation marked by a series of conformational changes related to binding of certain ligands, but the short of experimental structures has hampered the study of their activation mechanism and drug discovery. Our method, combining homology modeling, hierarchical sampling, and nanosecond-scale molecular dynamics, is one of the very few computational methods that can predict their active-state conformations and is one of the most computationally inexpensive. It enables the conformational landscape and the first quantitative energy landscape of GPCR activation to be efficiently mapped out.
This method, named ActiveGEnSeMBLE, allows the inactive- and active-state conformations of GPCRs without an experimental structure to be systematically predicted. We have validated the method with one of the most well-studied GPCRs, human β2 adrenergic receptor (hβ2AR), and applied the method on a GPCR without an experimental structure, human somatostatin receptor 5 (hSSTR5). Insights on GPCR activation as well as structure prediction methods are discussed.
The second one is a semiclassical approach for large-scale nonadiabatic dynamics of condensed systems in extreme conditions, termed Gaussian Hartree Approximated Quantum Mechanics (GHA-QM). Many nonadiabatic processes related to important applications (e.g. renewable energy) happen in large systems, but existing excited state dynamics methods are too computationally demanding for their long timescale simulations. GHA-QM is based on the electron force field (eFF) framework where we model electrons as Gaussian wavepackets and nuclei as classical point charges, and obtain a simplified solution to the time-dependent Schrödinger equation as the equation of motion. We employ a force field philosophy approximating the total energy as a sum of electronic kinetic energies, electrostatic energies and a Pauli correction, which corrects for the lack of explicit antisymmetry in the wavefunctions. New designs of the Pauli potential and preliminary results on hydrogen systems are discussed. With the new development, we hope to improve the accuracy and range of applications of eFF to simulate the nonadiabatic dynamics of hundreds of thousands of electrons on nanosecond timescale.</p
