Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born–Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.
The Hohenberg-Kohn theorem of density-functional theory establishes the existence of a bijection between the ground-state electron density and the external potential of a many-body system. This guarantees a one-to-one map from the electron density to all observables of interest including electronic excited-state energies. Time-Dependent Density-Functional Theory (TDDFT) provides one framework to resolve this map; however, the approximations inherent in practical TDDFT calculations, together with their computational expense, motivate finding a cheaper, more direct map for electronic excitations. Here, we show that determining density and energy functionals via machine learning allows the equations of TDDFT to be bypassed. The framework we introduce is used to perform the first excited-state molecular dynamics simulations with a machine-learned functional on malonaldehyde and correctly capture the kinetics of its excited-state intramolecular proton transfer, allowing insight into how mechanical constraints can be used to control the proton transfer reaction in this molecule. This development opens the door to using machine-learned functionals for highly efficient excited-state dynamics simulations.
more » « less- Award ID(s):
- 1955381
- NSF-PAR ID:
- 10380699
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 13
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract The description of π‐donor amido moieties as ‘weak‐field’ ligands can belie the influence of metal‐ligand covalency on the overall ligand field of coordination complexes, which can in turn influence properties including the magnetic ground state and those of their excited states. In this contribution, the ligand fields of pseudo‐octahedral Ni(II) complexes supported by diarylamido pincer‐type amido ligands – three previously reported examples supported by asymmetric (2‐R‐phenanthridin‐4‐yl)(8‐quinolinyl)amido ligands (R = Cl, CF3,
t Bu;R L1 ) along with a new congener bearing a symmetricbis (8‐quinolinyl)amido ligand (BQA;L2 ) – were investigated in two ways. First, high‐frequency and ‐field electron paramagnetic resonance spectroscopy (HFEPR), SQUID magnetometry, and electronic absorption spectroscopy were used to determine the ligand field parameters. Second, the ability to electrochemically address ligand‐based oxidations despite metal‐centered SOMOs in the parentS =1 paramagnets was investigated, supported by time‐dependent density functional theory (TDDFT) identification of strong intervalence charge‐transfer (IVCT) transitions attributed to electronic communication between two Namidomoieties mediated by a Ni(II) bridge. These findings are discussed in the broader context of 3d transition metal coordination complexes of weak‐field π‐donor ligands. -
more » « less
Abstract The theorems of density functional theory (DFT) establish bijective maps between the local external potential of a many-body system and its electron density, wavefunction and, therefore, one-particle reduced density matrix. Building on this foundation, we show that machine learning models based on the one-electron reduced density matrix can be used to generate surrogate electronic structure methods. We generate surrogates of local and hybrid DFT, Hartree-Fock and full configuration interaction theories for systems ranging from small molecules such as water to more complex compounds like benzene and propanol. The surrogate models use the one-electron reduced density matrix as the central quantity to be learned. From the predicted density matrices, we show that either standard quantum chemistry or a second machine-learning model can be used to compute molecular observables, energies, and atomic forces. The surrogate models can generate essentially anything that a standard electronic structure method can, ranging from band gaps and Kohn-Sham orbitals to energy-conserving ab-initio molecular dynamics simulations and infrared spectra, which account for anharmonicity and thermal effects, without the need to employ computationally expensive algorithms such as self-consistent field theory. The algorithms are packaged in an efficient and easy to use Python code, QMLearn, accessible on popular platforms.
-
more » « less
Hydrogen tunneling plays a critical role in many biologically and chemically important processes. The nuclear–electronic orbital multistate density functional theory (NEO-MSDFT) method was developed to describe hydrogen transfer systems. In this approach, the transferring proton is treated quantum mechanically on the same level as the electrons within multicomponent DFT, and a nonorthogonal configuration interaction scheme is used to produce delocalized vibronic states from localized vibronic states. The NEO-MSDFT method has been shown to provide accurate hydrogen tunneling splittings for fixed molecular systems. Herein, the NEO-MSDFT analytical gradients for both ground and excited vibronic states are derived and implemented. The analytical gradients and semi-numerical Hessians are used to optimize and characterize equilibrium and transition state geometries and to generate minimum energy paths (MEPs), for proton transfer in the deprotonated acetylene dimer and malonaldehyde. The barriers along the resulting MEPs are lower when the transferring proton is quantized because the NEO-MSDFT method inherently includes the zero-point energy of the transferring proton. Analysis of the proton densities along the MEPs illustrates that the proton density can exhibit symmetric or asymmetric bilobal character associated with symmetric or slightly asymmetric double-well potential energy surfaces and hydrogen tunneling. Analysis of the contributions to the intrinsic reaction coordinate reveals that changes in the C–O bond lengths drive proton transfer in malonaldehyde. This work provides the foundation for future reaction path studies and direct nonadiabatic dynamics simulations of a wide range of hydrogen transfer reactions.
-
more » « less
Abstract Density functional theory (DFT) is the standard formalism to study the electronic structure of matter at the atomic scale. In Kohn–Sham DFT simulations, the balance between accuracy and computational cost depends on the choice of exchange and correlation functional, which only exists in approximate form. Here, we propose a framework to create density functionals using supervised machine learning, termed NeuralXC. These machine-learned functionals are designed to lift the accuracy of baseline functionals towards that provided by more accurate methods while maintaining their efficiency. We show that the functionals learn a meaningful representation of the physical information contained in the training data, making them transferable across systems. A NeuralXC functional optimized for water outperforms other methods characterizing bond breaking and excels when comparing against experimental results. This work demonstrates that NeuralXC is a first step towards the design of a universal, highly accurate functional valid for both molecules and solids.
