Quantum tunneling in real space: Tautomerization of single porphycene molecules on the (111) surface of Cu, Ag, and Au

Anisotropic coupling of individual vibrational modes to a Cu(110) substrate

We present a study on the excitation of individual vibrational modes with ballistic charge carriers propagating along the Cu(110) surface. By means of the molecular nanoprobe technique, where the reversible switching of a molecule-in this case tautomerization of porphycene-is utilized to detect excitation events, we reveal anisotropic coupling of two distinct vibrational modes to the substrate. The N-H bending mode, excited below |E| ≈ 376 meV, exhibits maxima perpendicular to the rows of the Cu(110) substrate and minima along the rows. In contrast, the N-H stretching mode, excited above |E| ≈ 376 meV, displays maxima along the rows and is constant otherwise. This inversion of the anisotropy reflects the orthogonality between the N-H bending and stretching mode. Additionally, we observe an energy-dependent asymmetry in the propagation direction of charge carriers injected into the Cu(110) surface state. Hereby, the anisotropic band structure results in a correlation between the group velocity and the tunneling probability into electronic states of the substrate.

First-Principles Simulations of Tip Enhanced Raman Scattering Reveal Active Role of Substrate on High-Resolution Images

Tip-enhanced Raman scattering (TERS) has emerged as a powerful tool to obtain subnanometer spatial resolution fingerprints of atomic motion. Theoretical calculations that can simulate the Raman scattering process and provide an unambiguous interpretation of TERS images often rely on crude approximations of the local electric field. In this work, we present a novel and first-principles-based method to compute TERS images by combining Time Dependent Density Functional Theory (TD-DFT) and Density Functional Perturbation Theory (DFPT) to calculate Raman cross sections with realistic local fields. We present TERS results on free-standing benzene and C60 molecules, and on the TCNE molecule adsorbed on Ag(100). We demonstrate that chemical effects on chemisorbed molecules, often ignored in TERS simulations of larger systems, dramatically change the TERS images. This observation calls for the inclusion of chemical effects for predictive theory-experiment comparisons and an understanding of molecular motion at the nanoscale.

Tautomerization of single asymmetric oxahemiporphycene molecules on Cu(111)

We have studied 22-oxahemiporphycene molecules by a combination of scanning tunneling microscopy at low temperatures and density functional theory calculations. In contrast to other molecular switches with typically two switching states, these molecules can in principle exist in three different tautomers, due to their asymmetry and three inequivalent binding positions of a hydrogen atom in their macrocycle. Different tautomers are identified from the typical appearance on the surface and tunneling electrons can be used to tautomerize single molecules in a controllable way with the highest rates if the STM tip is placed close to the hydrogen binding positions in the cavity. Characteristic switching processes are explained by the different energy pathways upon adsorption on the surface. Upon applying higher bias voltages, deprotonation occurs instead of tautomerization, which becomes evident in the molecular appearance.

Spectroscopic investigation of photophysics and tautomerism of amino- and nitroporphycenes

Parent, unsubstituted porphycene and its two derivatives: 2,7,12,17-tetra-n-propylporphycene and 2,7,12,17-tetra-t-butylporphycene were substituted at the meso position with amino and nitro groups. These two families of porphycenes were characterized in detail with respect to their spectral, photophysical, and tautomeric properties. Two trans tautomers of similar energies coexist in the ground electronic state, but only one form dominates in the lowest excited singlet state. Absorption, magnetic circular dichroism (MCD), and emission anisotropy combined with quantum-chemical calculations led to the assignment of S₁ and S₂ transitions in both tautomers. Compared with the parent porphycene, the S₁-S₂ energy gap significantly increases; for one tautomeric form, the effect is twice as large as for the other. Both amino- and nitroporphycenes emit single fluorescence; previously reported dual emission of aminoporphycenes is attributed to a degradation product. Introduction of bulky t-butyl groups leads to a huge decrease in fluorescence intensity; this effect, arising from the interaction of the meso substituent with the adjacent t-butyl moiety, is particularly strong in the nitro derivative.

Determination of concerted or stepwise mechanism of hydrogen tunneling from isotope effects: Departure between experiment and theory

Isotope substitution is an important experimental technique that offers deep insight into reaction mechanisms, as the measured kinetic isotope effects (KIEs) can be directly compared with theory. For multiple proton transfer processes, there are two types of mechanisms: stepwise transfer and concerted transfer. The Bell-Limbach model provides a simple theory to determine whether the proton transfer mechanism is stepwise or concerted from KIEs. Recent STM experiments have studied the proton switching process in water tetramers on NaCl(001). Theoretical studies predict that this process occurs via a concerted mechanism, however, the experimental KIEs resemble the Bell-Limbach model for stepwise tunneling, raising question on the underlying mechanism or the validity of the model. We study this system using ab initio instanton theory, and in addition to thermal rates, we also considered microcanonical rates, as well as tunneling splittings. Instanton theory predicts a concerted mechanism, and the KIEs for tunneling rates (both thermal and microcanonical) upon deuteration are consistent with the Bell-Limbach model for concerted tunneling, but could not explain the experiments. For tunneling splittings, partial and full deuteration changes the size of it in a similar fashion to how it changes the rates. We further examined the Bell-Limbach model in another system, porphycene, which has both stepwise and concerted tunneling pathways. The KIEs predicted by instanton theory are again consistent with the Bell-Limbach model. This study highlights differences between KIEs in stepwise and concerted tunneling, and the discrepancy between theory and recent STM experiments. New theory/experiments are desired to settle this problem.

Influence of bulky substituents on single-molecule SERS sensitivity

The surface-enhanced Raman spectroscopy (SERS) detection limit strongly depends on the molecular structure, which we demonstrate for a family of tert-butyl-substituted porphycenes. Even though the investigated species present very similar photophysical properties, the ratio between the SERS signal and fluorescence background depends on the number of bulky tert-butyl groups. Moreover, the probability of single molecule detection systematically drops with the number of the moieties attached to the pyrrole ring. As steric hindrance is the only significantly changing feature among the studied chromophores, we attribute the observed phenomena to the spatial structure. We also show that the sensitivity of the SERS technique can be improved by lowering the temperature. We managed to observe single-molecule spectra for derivatives for which this was unattainable at room temperature.

Multidimensional Hydrogen Tunneling in Supported Molecular Switches: The Role of Surface Interactions

The nuclear tunneling crossover temperature (T_{c}) of hydrogen transfer reactions in supported molecular-switch architectures can lie close to room temperature. This calls for the inclusion of nuclear quantum effects (NQEs) in the calculation of reaction rates even at high temperatures. However, computations of NQEs relying on standard parametrized dimensionality-reduced models quickly become inadequate in these environments. In this Letter, we study the paradigmatic molecular switch based on porphycene molecules adsorbed on metallic surfaces with full-dimensional calculations that combine density-functional theory for the electrons with the semiclassical ring-polymer instanton approximation for the nuclei. We show that the double intramolecular hydrogen transfer (DHT) rate can be enhanced by orders of magnitude due to surface fluctuations in the deep-tunneling regime. We also explain the origin of an Arrhenius temperature dependence of the rate below T_{c} and why this dependence differs at different surfaces. We propose a simple model to rationalize the temperature dependence of DHT rates spanning diverse fcc [110] surfaces.

Real-Space Imaging of a Single-Molecule Monoradical Reaction

Organic radicals consisting of light elements exhibit a low spin-orbit coupling and weak hyperfine interactions with a long spin coherence length, which are crucial for future applications in molecular spintronics. However, the synthesis and characterization of these organic radicals have been a formidable challenge due to their chemical instability arising from unpaired electrons. Here, we report a direct imaging of the surface chemical transformation of an organic monoradical synthesized via the monodehydrogenation of a chemically designed precursor. Bond-resolved scanning tunneling microscopy unambiguously resolves various products formed through a complex structural dissociation and rearrangement of organic monoradicals. Density functional theory calculations reveal detailed reaction pathways from the monoradical to different cyclized products. Our study provides unprecedented insights into complex surface reaction mechanisms of organic radical reactions at the single molecule level, which may guide the design of stable organic radicals for future quantum technology applications.

Force-Activated Isomerization of a Single Molecule

Understanding and controlling isomerization at the single molecular level should provide new insight into the molecular dynamics and design guidelines of functional devices. Scanning tunneling microscopy (STM) has been demonstrated to be a powerful tool to study isomerization of single molecules on a substrate, by either electric field or inelastic electron tunneling mechanisms. A similar molecular isomerization process can in principle be induced by mechanical force; however, relevant study has remained elusive. Here, we demonstrate that isomerization of a N,N-dimethylamino-dianthryl-benzene molecule on Ag(100) can be mechanically driven by the STM tip. The existence of an out-of-plane dimethylamino group in the molecule is found to play a pivotal role in the isomerization process by providing a steric hindrance effect for asymmetric interaction between the STM tip and the molecule. This underlying mechanism is further confirmed by performing molecular dynamics simulations, which show agreement with experimental results. Our work opens the opportunity to manipulate the molecular configuration on the basis of mechanical force.

Temperature dependence of the vibrational spectrum of porphycene: a qualitative failure of classical-nuclei molecular dynamics

Approximate quantum dynamics succeed in predicting a temperature-dependent blue-shift of the high-frequency stretch bands that arise from vibrational coupling between low-frequency thermally activated modes and high-frequency quantized ones. Classical nuclei molecular dynamics fail and instead predict a red-shift.

Controlling Hydrogen-Transfer Rate in Molecules on Graphene by Tunable Molecular Orbital Levels

Molecular switches are building blocks of potential interest to store binary information, especially when they can be organized in periodic lattices. Among the variety of possible systems, switches based on hydrogen transfer are of special importance because they allow the switching operation to occur without severe conformational change that may interfere with neighboring molecular units. We have studied the excitation process of hydrogen transfer inside porphyrin molecules assembled on a graphene surface, using a low-temperature scanning tunneling microscope. We show that this hydrogen transfer is induced by an electronic resonant tunneling process through the molecular orbitals. Using nitrogen doping of graphene, we tune the rate of hydrogen transfer by shifting the molecular orbital energies owing to the charge transfer at nitrogen dopant sites in the graphene lattice. The control of the switching process allows the storage of information inside a molecular lattice, which is demonstrated by writing an artificial pattern inside a molecular island.

Deep Molecular Orbital Driven High-Temperature Hydrogen Tautomerization Switching

Hydrogen tautomerization molecular switches, a promising class of molecular components for the construction of complex nanocircuits, have been extensively studied using low-temperature scanning tunneling microscopy. However, these molecules are generally only reliably controllable in cryogenic environments, obstructing their utility in real devices. Here, we use dispersion-inclusive density functional theory and systematically investigate the adsorption and tautomerization behaviors of porphycene on six transition-metal surfaces. Among these surfaces, we found that hydrogen tautomerization on the Pt(110) surface corresponds to the largest switching barrier, allowing a controllable transition at high temperature. The switching behavior is closely related to the exceptional degree of charge transfer in the HOMO-2 orbital, illustrating the important role of deep orbital-surface interactions in porphycene molecular switching. Our work provides an in-depth understanding of the porphycene tautomerization mechanism and highlights new research avenues toward the practical application of molecular switches.

Real-Space Observation of Quantum Tunneling by a Carbon Atom: Flipping Reaction of Formaldehyde on Cu(110)

We present a direct observation of carbon-atom tunneling in the flipping reaction of formaldehyde between its two mirror-reflected states on a Cu(110) surface using low-temperature scanning tunneling microscopy (STM). The flipping reaction was monitored in real time, and the reaction rate was found to be temperature independent below 10 K. This indicates that this reaction is governed by quantum mechanical tunneling, albeit involving a substantial motion of the carbon atom (∼1 Å). In addition, deuteration of the formaldehyde molecule resulted in a significant kinetic isotope effect ( R_CH2O/ R_CD2O ≈ 10). The adsorption structure, reaction pathway, and tunneling probability were examined by density functional theory calculations, which corroborate the experimental observations.

Elucidating the Nuclear Quantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene

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.

Anharmonicity in a double hydrogen transfer reaction studied in a single porphycene molecule on a Cu(110) surface

Anharmonicity plays a crucial role in hydrogen transfer reactions in hydrogen-bonding systems, which leads to a peculiar spectral line shape of the hydrogen stretching mode as well as highly complex intra/intermolecular vibrational energy relaxation. Single-molecule study with a well-defined model is necessary to elucidate a fundamental mechanism. Recent low-temperature scanning tunnelling microscopy (STM) experiments revealed that the cis↔cis tautomerization in a single porphycene molecule on Cu(110) at 5 K can be induced by vibrational excitation via an inelastic electron tunnelling process and the N-H(D) stretching mode couples with the tautomerization coordinate [Kumagai et al. Phys. Rev. Lett. 2013, 111, 246101]. Here we discuss a pronounced anharmonicity of the N-H stretching mode observed in the STM action spectra and the conductance spectra. Density functional theory calculations find a strong intermode coupling of the N-H stretching with an in-plane bending mode within porphycene on Cu(110).