The Born–Oppenheimer approximation underlies much of chemical simulation and provides the framework defining the potential energy surfaces that are used for much of our pictorial understanding of ...chemical phenomena. However, this approximation breaks down when the dynamics of molecules in excited electronic states are considered. Describing dynamics when the Born–Oppenheimer approximation breaks down requires a quantum mechanical description of the nuclei. Chemical reaction dynamics on excited electronic states is critical for many applications in renewable energy, chemical synthesis, and bioimaging. Furthermore, it is necessary in order to connect with many ultrafast pump–probe spectroscopic experiments. In this review, we provide an overview of methods that can describe nonadiabatic dynamics, with emphasis on those that are able to simultaneously address the quantum mechanics of both electrons and nuclei. Such ab initio quantum molecular dynamics methods solve the electronic Schrödinger equation alongside the nuclear dynamics and thereby avoid the need for precalculation of potential energy surfaces and nonadiabatic coupling matrix elements. Two main families of methods are commonly employed to simulate nonadiabatic dynamics in molecules: full quantum dynamics, such as the multiconfigurational time-dependent Hartree method, and classical trajectory-based approaches, such as trajectory surface hopping. In this review, we describe a third class of methods that is intermediate between the two: Gaussian basis set expansions built around trajectories.
The increasing availability of high-quality experimental data and first-principles calculations creates opportunities for developing more accurate empirical force fields for simulation of proteins. ...We developed the AMBER-FB15 protein force field by building a high-quality quantum chemical data set consisting of comprehensive potential energy scans and employing the ForceBalance software package for parameter optimization. The optimized potential surface allows for more significant thermodynamic fluctuations away from local minima. In validation studies where simulation results are compared to experimental measurements, AMBER-FB15 in combination with the updated TIP3P-FB water model predicts equilibrium properties with equivalent accuracy, and temperature dependent properties with significantly improved accuracy, in comparison with published models. We also discuss the effect of changing the protein force field and water model on the simulation results.
We develop an extension of the variational quantum eigensolver (VQE) algorithm-multistate contracted VQE (MC-VQE)-that allows for the efficient computation of the transition energies between the ...ground state and several low-lying excited states of a molecule, as well as the oscillator strengths associated with these transitions. We numerically simulate MC-VQE by computing the absorption spectrum of an ab initio exciton model of an 18-chromophore light-harvesting complex from purple photosynthetic bacteria.
Few would dispute that theoretical chemistry tools can now provide keen insights into chemical phenomena. Yet the holy grail of efficient and reliable prediction of complex reactivity has remained ...elusive. Fortunately, recent advances in electronic structure theory based on the concepts of both element- and rank-sparsity, coupled with the emergence of new highly parallel computer architectures, have led to a significant increase in the time and length scales which can be simulated using first principles molecular dynamics. This opens the possibility of new discovery-based approaches to chemical reactivity, such as the recently proposed ab initio nanoreactor. We argue that due to these and other recent advances, the holy grail of computational discovery for complex chemical reactivity is rapidly coming within our reach.
Hybrid quantum mechanical–molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of ...key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500–600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11–16 residues or 200–300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.
The development of accurate molecular mechanics force fields is a significant challenge that must be addressed for the continued success of molecular simulation. We developed the ForceBalance method ...to automatically derive accurate force field parameters using flexible combinations of experimental and theoretical reference data. The method is demonstrated in the parametrization of two rigid water models, yielding new parameter sets (TIP3P-FB and TIP4P-FB) that accurately describe many physical properties of water.
Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and ...experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ∼2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.
Biological systems rely on recyclable materials resources such as amino acids, carbohydrates and nucleic acids. When biomaterials are damaged as a result of aging or stress, tissues undergo repair by ...a depolymerization-repolymerization sequence of remodelling. Integration of this concept into synthetic materials systems may lead to devices with extended lifetimes. Here, we show that a metastable polymer, end-capped poly(o-phthalaldehyde), undergoes mechanically initiated depolymerization to revert the material to monomers. Trapping experiments and steered molecular dynamics simulations are consistent with a heterolytic scission mechanism. The obtained monomer was repolymerized by a chemical initiator, effectively completing a depolymerization-repolymerization cycle. By emulating remodelling of biomaterials, this model system suggests the possibility of smart materials where aging or mechanical damage triggers depolymerization, and orthogonal conditions regenerate the polymer when and where necessary.