Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory

Proton transfer and hydrogen tunneling play key roles in many processes of chemical and biological importance. The generalized nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method was developed in order to capture hydrogen tunneling effects in systems involving the transfer and tunneling of one or more protons. The generalized NEO-MSDFT method treats the transferring protons quantum mechanically on the same level as the electrons and obtains the delocalized vibronic states associated with hydrogen tunneling by mixing localized NEO-DFT states in a nonorthogonal configuration interaction scheme. Herein, we present the derivation and implementation of analytical gradients for the generalized NEO-MSDFT vibronic state energies and the nonadiabatic coupling vectors between these vibronic states. We use this methodology to perform adiabatic and nonadiabatic dynamics simulations of the double proton transfer reactions in the formic acid dimer and the heterodimer of formamidine and formic acid. The generalized NEO-MSDFT method is shown to capture the strongly coupled synchronous or asynchronous tunneling of the two protons in these processes. Inclusion of vibronically nonadiabatic effects is found to significantly impact the double proton transfer dynamics. This work lays the foundation for a variety of nonadiabatic dynamics simulations of multiple proton transfer systems, such as proton relays and hydrogen-bonding networks.