Evidence for Dominant Role of Tunneling in Condensed Phases and at High Temperatures: Double Hydrogen Transfer in Porphycenes

Nuclear Quantum Effects in Proton or Hydrogen Transfer

Proton or hydrogen transfers, basic chemical reactions, proceed either by thermally activated barrier crossing or via tunneling. Studies of molecules undergoing single or double proton or hydrogen transfer in the ground or excited electronic state reveal that tunneling can dominate under conditions usually considered to favor the thermal process. Moreover, the tunneling probability strongly varies for excitation of certain vibrational modes, which changes the effective barrier and/or proton transfer distance. When the reaction is fast compared to vibrational relaxation, the mode selectivity can still be maintained for molecules in solutions at 293 K. These observations point to dangers of relating the calculated minimum energy paths and the associated barriers to the experimentally obtained activation energies. The multidimensional character of the reaction coordinate is obvious; it can dramatically change for slowly and rapidly relaxing environments. We postulate that the hydrogen bond definition should be extended by specifically including the role of molecular vibrations.

Coupling between tautomerism and radiationless deactivation in porphycenes

Quantum yields of fluorescence of porphycenes – porphyrin isomers – can vary by orders of magnitude, even for very similar derivatives, such as meso-dimethyl- vs. meso-tetramethylporphycene. In weakly emitting porphycenes the fluorescence intensity strongly depends on viscosity and can be recovered by placing a molecule in a rigid environment. We postulate that the efficient nonradiative deactivation is due to the quantum effect, delocalization of the inner protons. The delocalization, which increases with the strength of intramolecular hydrogen bonds may induce structural changes that lead to distortion from planarity and, as a result, efficient S0 ← S1 internal conversion. The effect seems to be general, as indicated by good correlation between the quantum yield of fluorescence and the distance between H-bonded nitrogen atoms, the latter being a reliable measure of hydrogen bonding strength. Based on the available photophysical and X-ray data, such correlation was found so far for over 20 differently substituted porphycenes.

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.

Computational Evidence for Tunneling and a Hidden Intermediate in the Biosynthesis of Tetrahydrocannabinol

Quantum tunneling is computed for a reaction sequence that models the conversion of the ortho-quinone methide of cannabigerolic acid 1 to the decarboxylated product (-)-trans-Δ⁹-tetrahydrocannabinol (THC, 3). This calculation is the first to evaluate multidimensional tunneling in this sequence. Computations were carried out with POLYRATE and GAUSSRATE using B3LYP/6-31G(d,p) to examine the mechanism of THC 3 formation. The pentyl chain on THC 3 and its precursors were replaced with a methyl group to compute tunneling contributions to the rates of four separate steps: (i) initial Diels-Alder reaction of the quinone methide with the trisubstituted alkene end-group of the geranyl 1Z-CH₃ to give 2Z-CH₃, (ii) acid-catalyzed keto-enol tautomerization, which converts 2rZ-CH₃ to 4rZ-CH₃, (iii) carboxyl rotamerization converting 4rZ-CH₃ to 4E-CH₃, and (iv) decarboxylation that converts 4E-CH₃ to 3-CH₃. Tunneling contributions to the rate constants of steps (i)-(iv) are between 19 and 76% at 293 K. In step (ii), nonuniform changes in the zero-point vibrational energy along the reaction path created a shallow minimum in the 0 K free energy. It is a hidden intermediate because it is not a minimum on the potential energy surface and is detectable only when zero-point energy is taken into account along the reaction path. Predicted kinetic isotope effects would be experimentally observable at temperatures that are convenient to use. This is particularly relevant in the decarboxylation stage of the reaction sequence and has important implications because of its role in THC 3 formation.

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.

Theoretical Study of Proton Tunneling in the Imidazole-Imidazolium Complex

Proton tunneling in the hydrogen-bonded imidazole–imidazolium complex ion has been studied theoretically. Ab initio CASSCF/6-311++G(d,p) calculations concerning geometry optimization and vibrational frequencies have been carried out for equilibrium and transition state structures of the system. Two-dimensional double-well model potentials were constructed on the basis of ab initio results and used to analyze the proton dynamics in the hydrogen bond and the influence of the excitation of low-frequency hydrogen-bond vibrations on the proton tunneling splittings. The energy of tunneling-split vibrational sublevels of the high-frequency tunneling mode have been calculated for its ground and first excited vibrational state for the series of excitations of the coupled low-frequency intramolecular hydrogen-bond modes. The promoting and suppressing effect of the low-frequency modes on the proton splittings was shown in the ground and first excited vibrational state of the tunneling mode. The vibrational sublevels form the two separate semicontinuous bands between which the absorption transitions may occur. This mechanism explains the experimentally observed splitting and doublet-component broadening of the high-frequency N–H stretching infrared (IR) absorption band.

Controlling Emissive Properties by Intramolecular Hydrogen Bonds: Alkyl and Aryl meso-Substituted Porphycenes

Porphycene, a porphyrin isomer, is an efficient fluorophore. However, four-fold meso substitution with alkyl groups decreases the fluorescence quantum yield by orders of magnitude. For aryl substituents, this effect is small. To explain this difference, we have synthesized and studied a mixed aryl-alkyl-substituted compound, 9,20-diphenyl-10,19-dimethylporphycene, as well as the 9,20-diphenyl and 9,20-dimethyl derivatives. Analysis of the structural, spectroscopic, and photophysical data of the six porphycenes, combined with quantum chemical calculations, shows a clear correlation between the strength of the intramolecular NH⋅⋅⋅N hydrogen bonds and the efficiency of the radiationless depopulation of the lowest-excited singlet state. This result led us to propose a model in which the delocalization of the inner protons in the cavity of the macrocycle is responsible for the nonradiative deactivation channel. The applicability of the model is confirmed by the literature data for other alkyl- or aryl-substituted porphycenes. The finding of a correlation between structural and emissive characteristics enables a rational design of porphycenes with desired photophysical properties.

2 + 2 Can Make Nearly a Thousand! Comparison of Di- and Tetra-Meso-Alkyl-Substituted Porphycenes

Two porphycenes, substituted at the meso positions with two and four methyl groups, respectively, reveal similar absorption spectra, but their photophysical properties are completely different. 9,20-dimethylporphycene emits fluorescence with about 20% quantum yield, independent of the solvent. In contrast, fluorescence of 9,10,19,20-tetramethylporphycene is extremely weak in nonviscous solvents, but it can be recovered by placing the chromophore in a rigid environment. We propose a model that explains these differences, based on calculations and structural analogies with other extremely weakly emitting derivatives, dibenzo[cde,mno]porphycenes. The efficient S₁ deactivation involves delocalization of two inner cavity protons coupled with proton translocation toward a high-energy cis tautomer. The latter process leads to distortion from planarity. The probability of deactivation increases with the strength of the intramolecular NH···N hydrogen bonds. The model also explains the observation of biexponential fluorescence decay in weakly emitting porphycenes. It can be extended to other derivatives, in particular, the asymmetrically substituted ones. We also point to the possibility of using specific porphycenes as viscosity sensors, in particular, when working in single molecule regime.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

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.

Influence of local microenvironment on the double hydrogen transfer in porphycene

We performed time-resolved transient absorption and fluorescence anisotropy measurements in order to study tautomerization of porphycene in rigid polymer matrices at cryogenic temperatures. Studies were carried out in poly(methyl methacrylate) (PMMA), poly(vinyl butyral) (PVB), and poly(vinyl alcohol) (PVA). The results prove that in all studied media hydrogen tunnelling plays a significant role in the double hydrogen transfer which becomes very sensitive to properties of the environment below approx. 150 K. We also demonstrate that there exist two populations of porphycene molecules in rigid media: "hydrogen-transferring" molecules, in which tautomerization occurs on time scales below 1 ns and "frozen" molecules in which double hydrogen transfer is too slow to be monitored with nanosecond techniques. The number of "frozen" molecules increases when the sample is cooled. We explain this effect by interactions of guest molecules with a rigid host matrix which disturbs symmetry of porphycene and hinders tunnelling. Temperature dependence of the number of hydrogen-transferring molecules suggests that the factor which restores the symmetry of the double-minimum potential well in porphycene are intermolecular vibrations localized in separated regions of the amorphous polymer.

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.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

2,7,12,17-Tetra(2,5-thienylene)-substituted porphycenes

We report syntheses of thiophene and dithiophene-substituted porphycenes (ThPc and DThPc) at 2,7,12,17-positions by McMurry coupling. The crystal structure of ThPc revealed that the porphycene plane shows a highly planar structure, and the dihedral angles between the porphycene core and thiophene are relatively small at 21[Formula: see text] and 18[Formula: see text]. ThPc and DThPc exhibit red-shifted and broadened absorption because of the extension of [Formula: see text] conjugations through porphycene to the substituted thiophenes. We found that introduction of thiophene units onto porphycene results in decreasing the HOMO–LUMO differences effectively.

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.

Supersonic jet spectroscopy of parent hemiporphycene: Structural assignment and vibrational analysis for S0 and S1 electronic states

Hemiporphycene (HPc), a constitutional isomer of porphyrin, is studied under supersonic expansion conditions by means of laser-induced fluorescence, visible-visible hole-burning experiments, single vibronic level fluorescence techniques, and quantum chemical calculations. Only one trans form of jet-cooled HPc is observed, in contrast to solution studies that evidence a mixture of two trans tautomeric forms separated in energy by ∼1 kcal/mol. Reliable structural assignment is provided by simulating absorption and emission patterns at the density functional theory and time-dependent density functional theory levels of theory. The vibronic spectra are nicely reproduced for both electronic ground and lowest excited singlet states for the most stable trans form. In contrast to another porphyrin isomer, porphycene (Pc), no tunneling or photo-induced hydrogen transfer is detected. The lower symmetry of HPc compared with Pc and the concomitant non-equivalent positions of the inner-cavity nitrogen atoms result in a non-symmetric double minimum potential for tautomerization, larger energy barrier, and a longer tunneling distance, with the average intramolecular hydrogen bond length larger in HPc than in Pc. HPc readily forms hydrates that show red-shifted absorption relative to the bare molecule.

Nature of Large Temporal Fluctuations of Hydrogen Transfer Rates in Single Molecules

Double hydrogen transfer was monitored in single molecules of parent porphycene and its tetra- t-butyl derivative using confocal fluorescence microscopy. The molecules have been embedded in a polymer matrix. Under such conditions, a significant fraction of the population reveals a huge decrease of the tautomerization rate with respect to the value obtained from ensemble studies in solution. This effect is explained by a model that assumes that the rate is determined by the reorganization coordinate that involves slow relaxation of the polymer matrix. The model provides indirect evidence for the dominant role of tunneling. It is proposed that tautomerization in single molecules of the porphycene family can be used to probe polymer relaxation dynamics on the time scale ranging from picoseconds to minutes.

Direct observation of single-molecule hydrogen-bond dynamics with single-bond resolution

The hydrogen bond represents a fundamental interaction widely existing in nature, which plays a key role in chemical, physical and biochemical processes. However, hydrogen bond dynamics at the molecular level are extremely difficult to directly investigate. Here, in this work we address direct electrical measurements of hydrogen bond dynamics at the single-molecule and single-event level on the basis of the platform of molecular nanocircuits, where a quadrupolar hydrogen bonding system is covalently incorporated into graphene point contacts to build stable supramolecule-assembled single-molecule junctions. The dynamics of individual hydrogen bonds in different solvents at different temperatures are studied in combination with density functional theory. Both experimental and theoretical results consistently show a multimodal distribution, stemming from the stochastic rearrangement of the hydrogen bond structure mainly through intermolecular proton transfer and lactam-lactim tautomerism. This work demonstrates an approach of probing hydrogen bond dynamics with single-bond resolution, making an important contribution to broad fields beyond supramolecular chemistry.

Unusual effects in single molecule tautomerization: hemiporphycene

Temporal evolution of the fluorescence spectra of individual hemiporphycene molecules reveals unusual double hydrogen transfer reaction dynamics.

Proton in a Confined Space: Structural Studies of H+ ⊂Crypt-111 Iodide and Some Halogen-Bonded Derivatives

Experimental observations and modeling data are reported on the solid-state structural features of crypt- 111⋅HI (1) and the three-component co-crystals that 1 forms with α,ω-diiodoperfluoroalkanes 2 a-d. X-ray analyses indicate that, in all five systems and at low temperature, the caged proton is covalently bonded to a single nitrogen atom and is involved in a network of intramolecular hydrogen bonds. In contrast, room-temperature, solid-state ¹⁵ N NMR spectroscopy suggests magnetic equivalency of the two N atoms of crypt-111 in both 1 and co-crystals of 1 with diiodoperfluoroalkanes. Computational modelling confirms that the acidic hydrogen inside the cavity preferentially sits along the internitrogen axis and is covalently bonded to one nitrogen. The computed energy barriers suggest that the hopping of the encapsulated proton between the two N atoms of the cage can occur in the halogen-bonded co-crystals of 1⋅2, but it is hardly possible in the pure H⁺ ⊂crypt-111 iodide 1. These different pictures of the proton position and dynamics obtained by using different techniques and conditions confirm the unique characteristics of the confined space within the cavity of crypr-111 and the distinctive features of processes occurring therein.

Non-typical fluorescence studies of excited and ground state proton and hydrogen transfer

Fluorescence studies of tautomerization have been carried out for various systems that exhibit single and double proton or hydrogen translocation in various environments, such as liquid and solid condensed phases, ultracold supersonic jets, and finally, polymer matrices with single emitters. We focus on less explored areas of application of fluorescence for tautomerization studies, using porphycene, a porphyrin isomer, as an example. Fluorescence anisotropy techniques allow investigations of self-exchange reactions, where the reactant and product are formally identical. Excitation with polarized light makes it possible to monitor tautomerization in single molecules and to detect their three-dimensional orientation. Analysis of fluorescence from single vibronic levels of jet-isolated porphycene not only demonstrates coherent tunneling of two internal protons, but also indicates that the process is vibrational mode-specific. Next, we present bifunctional proton donor-acceptor systems, molecules that are able, depending on the environment, to undergo excited state single intramolecular or double intermolecular proton transfer. For molecules that have donor and acceptor groups located in separate moieties linked by a single bond, excited state tautomerization can be coupled to mutual twisting of the two subunits.

Spectroscopic and microscopic investigations of tautomerization in porphycenes: condensed phases, supersonic jets, and single molecule studies

Tautomerization of porphycene, coherent in supersonic jets and a rate process in solutions, can be controlled for single molecules on surfaces.

Non-classical effects in proton or hydrogen transfer

Results of thorough investigations of tautomerism in the ground and excited electronic states performed for various organic molecules under different regimes: condensed phases, isolated and single molecules, with a special emphasis on porphycene, a porphyrin isomer, demonstrate that, in order to precisely describe the proton/hydrogen transfer path, one has to explicitly consider tunneling. Tautomerization is a multidimensional process, controlled by excitation of specific vibrational modes. Vibrational excitation can both enhance or hinder the reaction. The role of specific vibrational modes can now be assessed, even at the level of single molecules.

Spectroscopy and Tautomerization Studies of Porphycenes

Tautomerization in porphycenes, constitutional isomers of porphyrins, is strongly entangled with spectral and photophysical parameters. The intramolecular double hydrogen transfer occurring in the ground and electronically excited states leads to uncommon spectroscopic characteristics, such as depolarized emission, viscosity-dependent radiationless depopulation, and vibrational-mode-specific tunneling splittings. This review starts with documentation of the electronic spectra of porphycenes: Absorption and magnetic circular dichroism are discussed, together with their analysis based on the perimeter model. Next, photophysical characteristics are presented, setting the stage for the final part, which discusses the developments in research on tautomerism. Porphycenes have been studied in different experimental regimes: molecules in condensed phases, isolated in supersonic jets and helium nanodroplets, and, recently also on the level of single molecules investigated by optical and scanning probe microscopies. Because of the rich and detailed information obtained from these diverse investigations, porphycenes emerge as very good models for studying the complex, multidimensional phenomena involved in the process of intramolecular double hydrogen transfer.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).