The dynamics of water dimers was investigated at the single-molecule level by using a scanning tunneling microscope. The two molecules in a water dimer, bound on a Cu(110) surface at 6 K, were observed to exchange their roles as hydrogen-bond donor and acceptor via hydrogen-bond rearrangement. The interchange rate is $\sim 60$ times higher for $({\mathrm{H}}_2\mathrm{O}{)}_2$ than for $({\mathrm{D}}_2\mathrm{O}{)}_2$, suggesting that quantum tunneling is involved in the process. The interchange rate is enhanced upon excitation of the intermolecular mode that correlates with the reaction coordinate.

We address the double hydrogen transfer (DHT) dynamics of the porphycene molecule, a complex paradigmatic system in which the making and breaking of H-bonds in a highly anharmonic potential energy surface require 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 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 × 10¹¹ s^–1 at 150 K and 0.63 × 10¹¹ 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 thermostated ring-polymer molecular dynamics simulations. The position and line shape of this peak, centered at around 2600 cm^–1 and spanning 750 cm^–1, stem 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.

We report the direct observation of intramolecular hydrogen atom transfer reactions (tautomerization) within a single porphycene molecule on a Cu(110) surface by scanning tunneling microscopy. It is found that the tautomerization can be induced via inelastic electron tunneling at 5 K. By measuring the bias-dependent tautomerization rate of isotope-substituted molecules, we can assign the scanning tunneling microscopy-induced tautomerization to the excitation of specific molecular vibrations. Furthermore, these vibrations appear as characteristic features in the $dI/dV$ spectra measured over individual molecules. The vibrational modes that are associated with the tautomerization are identified by density functional theory calculations. At higher temperatures above $\sim 75\mathrm{K}$, tautomerization is induced thermally and an activation barrier of about 168 meV is determined from an Arrhenius plot.