Elucidating the NuclearQuantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene

Abstract

We address the double hydrogen transfer (DHT) dynamics of the porphycene molecule: A complex paradigmatic system where the making and breaking of H-bonds in a highly anharmonic potential energy surface requires a quantum mechanical treatment not only of the electrons, but also of the nuclei. We combine density-functional theory calculations, employing hybrid functionals and van der Waals corrections, with recently proposed and optimized path-integral ring-polymer methods for the approximation of quantum vibrational spectra and reaction rates. Our full-dimensional ring-polymer instanton simulations show that below 100 K the concerted DHT tunneling pathway dominates, but between 100 K and 300 K there is a competition between concerted and stepwise pathways when nuclear quantum effects are included. We obtain ground-state reaction rates of $2.19\times 10^{11}{\mathrm{s}}^{-1}$ at 150 K and $0.63\times 10^{11}{\mathrm{s}}^{-1}$ at 100 K, in good agreement with experiment. We also reproduce the puzzling N-H stretching band of porphycene with very good accuracy from thermostatted ring-polymer molecular dynamics simulations. The position and lineshape of this peak, centered at around 2600 cm${}^{-1}$ and spanning 750 cm${}^{-1}$, stems from a combination of very strong H-bonds, the coupling to low-frequency modes, and the access to $cis$-like isomeric conformations, which cannot be appropriately captured with classical-nuclei dynamics. These results verify the appropriateness of our general theoretical approach and provide a framework for a deeper physical understanding of hydrogen transfer dynamics in complex systems.

References in corpus (8)

• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• How to remove the spurious resonances from ring polymer molecular dynamics
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Atom-Tunneling in Chemistry
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• Dual-Level Approach to Instanton Theory
• A diabatic state model for double proton transfer in hydrogen bonded complexes

Cited by in corpus (19)

• Machine Learning Force Fields
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Inexpensive modelling of quantum dynamics using path integral generalized Langevin equation thermostats
• Dynamical Strengthening of Covalent and Non-Covalent Molecular Interactions by Nuclear Quantum Effects at Finite Temperature
• Instanton formulation of Fermi's golden rule in the Marcus inverted regime
• Quantum dynamics using path integral coarse-graining
• Semiclassical instanton formulation of Marcus-Levich-Jortner theory
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Nonadiabatic instanton rate theory beyond the golden-rule limit
• Significance of nuclear quantum effects in hydrogen bonded molecular chains
• Multidimensional Hydrogen Tunneling in Supported Molecular Switches: The Role of Surface Interactions
• Dissipative Tunneling Rates through the Incorporation of First-Principles Electronic Friction in Instanton Rate Theory I: Theory
• Divide-and-Conquer Method for Instanton Rate Theory
• Exact tunneling splittings from symmetrized path integrals
• Classical and quantum dynamical manifestations of index-2 saddles: concerted versus sequential reaction mechanisms
• Aqueous Solution Chemistry In Silico and the Role of Data Driven Approaches
• Quantum rates in dissipative systems with spatially varying friction
• Third-order many-body expansion of OSV-MP2 wavefunction for low-order scaling analytical gradient computation