Evidence for Two Forms, Double Hydrogen Tunneling, and Proximity of Excited States in Bridge-Substituted Porphycenes: Supersonic Jet Studies

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.

An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm^-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

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.

Hydrogen Delocalization in an Asymmetric Biomolecule: The Curious Case of Alpha-Fenchol

Rotational microwave jet spectroscopy studies of the monoterpenol α-fenchol have so far failed to identify its second most stable torsional conformer, despite computational predictions that it is only very slightly higher in energy than the global minimum. Vibrational FTIR and Raman jet spectroscopy investigations reveal unusually complex OH and OD stretching spectra compared to other alcohols. Via modeling of the torsional states, observed spectral splittings are explained by delocalization of the hydroxy hydrogen atom through quantum tunneling between the two non-equivalent but accidentally near-degenerate conformers separated by a low and narrow barrier. The energy differences between the torsional states are determined to be only 16(1) and 7(1) cm-1hc for the protiated and deuterated alcohol, respectively, which further shrink to 9(1) and 3(1) cm-1hc upon OH or OD stretch excitation. Comparisons are made with the more strongly asymmetric monoterpenols borneol and isopinocampheol as well as with the symmetric, rapidly tunneling propargyl alcohol. In addition, the third-in contrast localized-torsional conformer and the most stable dimer are assigned for α-fenchol, as well as the two most stable dimers for propargyl alcohol.

Simple models for the quick estimation of ground state hydrogen tunneling splittings in alcohols and other compounds

Models for the quick estimation of energy splittings caused by coherent tunneling of hydrogen atoms are evaluated with available experimental data for alcohols and improvements are proposed. The discussed models are mathematically simple and require only results from routine quantum chemical computations, i.e. hybrid DFT calculation of the equilibrium geometry and the transition state within the harmonic approximation. A benchmark of experimental splittings spanning four orders of magnitude for 27 alcohol species is captured by three evaluated models with a mean symmetric deviation factor of 1.7, 1.5 and 1.4, respectively, i.e. the calculated values deviate on average by this factor in either direction. Limitations of the models are explored with alcohols featuring uncommon properties, such as an inverted conformational energy sequence, a very light molecular frame, an elevated torsional frequency, or a coupling with a second internal degree of freedom. If the splitting of either the protiated or deuterated form of an alcohol is already experimentally determined, the one of the second isotopolog can be estimated by three additional models with a mean symmetric deviation factor of 1.14, 1.19 and 1.15, respectively. It is shown that this can be achieved with a novel approach without any quantum chemical calculation by directly correlating experimental splittings of isotopologs across related species. This is also demonstrated for other classes of compounds with hydrogen tunneling, such as amines, thiols, and phenols. Furthermore, it is found that the isotope effect can even be anticipated without any further knowledge about the system solely from the size of either splitting with a mean symmetric deviation factor of 1.3. This is based on an extensive sample of 77 pairs of splittings spanning eight orders of magnitude for isotopologs of chemically diverse compounds.

Double intramolecular hydrogen transfer assisted dual emission in a carbazole-embedded porphyrin-like macrocycle

The introduction of a pyrrole ring at one of the meso positions of carbazole-based porphyrins lowers the structural symmetry and results in dual emission, which strongly depends on the excitation wavelength and temperature. The origin of dual emission induced by NH-tautomerism is confirmed via photophysical and DFT calculations.

Efficient Quantum Mechanical Calculations of Mode-Specific Tunneling Splittings upon Fundamental Excitation in the Dimer of Formic Acid

The formic acid dimer (FAD) is an important benchmark system for understanding the double hydrogen transfer process. Most recently, Zhang et al. measured a few tunneling splittings upon fundamental excitation of FAD precisely (Zhang, Y. et al. J. Chem. Phys. 2017, 146, 244306); however, relevant theoretical studies are very limited. Here, we present a multidimensional quantum dynamics study on mode-specific tunneling splittings upon fundamental excitation in FAD with an efficient theoretical scheme developed by our group in which the process-oriented basis function customization strategy is combined with the preconditioned inexact spectral transform method. Various mode-specific tunneling splittings upon fundamental excitation are systematically calculated, and interesting mode-specific excitation effects on tunneling rate are identified. In particular, the calculated tunneling splittings for the ν₂₂ and ν₂₁ states are in good agreement with experiment, and the remarkable mode-specific suppression effects upon excitation should result from that the antisymmetric vibrational modes hinder the concerted double H-transfer. The present work is helpful to acquire a better understanding of the mode-specific excitation effects on tunneling processes.

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.

Direct Visualization of Hydrogen-Transfer Intermediate States by Scanning Tunneling Microscopy

Hydrogen atoms bonded within molecular cavities often undergo tunneling or thermal-transfer processes that play major roles in diverse physical phenomena. Such transfers may or may not entail intermediate states. The existence of such fleeting states is typically determined by indirect means, while their direct visualization has not been achieved, largely because their concentrations under equilibrium conditions are negligible. Here we use density-functional-theory calculations and scanning-tunneling-microscopy (STM) image simulations to predict that, under specially designed nonequilibrium conditions of voltage-enhanced high transfer rates, the cis-intermediate of the two-hydrogen transfer process in metal-free naphthalocyanine molecules adsorbed on Ag(111) surfaces would be visualizable in a composite image of double-C morphology. As guided by the theoretical predictions, at adjusted scanning temperature and bias, STM experiments achieve a direct visualization of the cis-intermediate. This work demonstrates a practical way to directly visualize elusive intermediates, which enhances understanding of the quantum dynamics of hydrogen atoms.

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.

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.

Expanded, Contracted, and Isomeric Porphyrins: Theoretical Aspects

The use of cyclic polyene perimeter-model approaches, such as Gouterman's four-orbital model and Michl's perimeter model, to analyze trends in the electronic structures and optical properties of expanded, contracted, and isomeric porphyrins is described with an emphasis on the use of magnetic circular dichroism (MCD) spectroscopy to validate the results of TD-DFT calculations. Trends in the electronic structures and optical properties of isomeric porphyrins are examined by comparing the properties of porphycenes, corrphycenes, hemiporphycenes, isoporphycenes, N-confused and neoconfused porphyrins, and norroles, whereas those of ring-contracted porphyrins are examined by comparing the properties of subporphyrins, triphyrins, and vacataporphyrins. The ring-expanded compounds that are examined include cyclo[n]pyrroles, [22]pentaphyrins(1.1.1.1.1), sapphyrins, smaragdyrins, isosmaragdyrins, orangarins, ozaphyrins, [26]hexaphyrins(1.1.1.1.1.1), rubyrins, rosarins, amethyrins, isoamethyrins, bronzaphyrins, and doubly N-confused hexaphyrins.

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.

Single molecule Raman spectra of porphycene isotopologues

Single molecule surface-enhanced resonance Raman scattering (SERRS) spectra have been obtained for the parent porphycene (Pc-d0) and its deuterated isotopologue (Pc-d12), located on gold and silver nanoparticles. Equal populations of "hot spots" by the two isotopologues are observed for 1 : 1 mixtures in a higher concentration range of the single molecule regime (5 × 10(-9) M). For decreasing concentrations, hot spots are preferentially populated by undeuterated molecules. This is interpreted as an indication of a lower surface diffusion coefficient of Pc-d12. The photostability of single Pc molecules placed on nanoparticles is strongly increased in comparison with polymer environments. Trans tautomeric species dominate the spectra, but the analysis of time traces reveals transient intermediates, possibly due to rare cis tautomeric forms.

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

Investigation of the double hydrogen transfer in porphycene, its 2,7,12,17-tetra-tert-butyl derivative, and their N-deuterated isotopologues revealed the dominant role of tunneling, even at room temperature in condensed phase. Ultrafast optical spectroscopy with polarized light employed in a wide range of temperatures allowed the identification and evaluation of contributions of two tunneling modes: vibrational ground-state tunneling, occurring from the zero vibrational level, and vibrationally activated, via a large amplitude, low-frequency mode. Good correspondence was found between the rates of incoherent tunneling occurring in condensed phase and the values estimated on the basis of tunneling splittings observed in molecules isolated in supersonic jets or helium nanodroplets. The results provide solid experimental insight into widely proposed quantum facets of ubiquitous hydrogen-transfer phenomena.

Calculations of Mode-Specific Tunneling of Double-Hydrogen Transfer in Porphycene Agree with and Illuminate Experiment

We report a theoretical study of mode-specific tunneling splittings in double-hydrogen transfer in trans-porphycene. We use a novel, mode-specific "Qim path method", in which the reaction coordinate is the imaginary-frequency normal mode of the saddle point separating the equivalent minima. The model considers all 108 normal modes and uses no adjustable parameters. The method gives the ground vibrational-state tunneling splitting, as well the increase in the splitting upon excitation of certain modes, in good agreement with experiment. Interpretation of these results is also transparent with this method. In addition, predictions are made for mode excitations not investigated experimentally. Results for d1 and d2 isotopolgues are also in agreement with experiment.

Double Hydrogen Transfer in Low Symmetry Porphycenes

The rate constants of intramolecular double hydrogen transfer have been determined for of tert-butyl-substituted porphycenes: 2-tert-butyl-, 2,7-di-tert-butyl-, 2,7,12-tri-tert-butyl-, and 2,7,12,17-tetra-tert-butylporphycene using femtosecond pump-probe polarization spectroscopy. The rates increase monotonically with the number of substituents. As is the case for other porphycenes studied so far, the tautomerization becomes slower after excitation to the lowest excited singlet state. The potential for the trans-trans self-exchange reaction is symmetrical in di- and tetra-tert-butyl substituted derivatives, but not for the singly and triply substituted ones. Our studies enabled determination of the relative populations of nonequivalent tautomers and thus of the equilibrium constants in S0 and S1 states, as well as estimation of ground and excited state activation energies for the tautomerization process.

Reaction of cobalt porphycene with hydride reagents: spectroscopic detection of Co–H porphycene species and formation of Co–SnR3 porphycene species

The reaction of tetrapropylporphycenatocobalt(III) with tributyltin hydride generates a cobalt(III)–hydride porphycene detectable by UV-vis spectroscopy under diluted conditions, whereas it is impossible to characterize hydride species of cobalt porphyrins. One of the reasons for the stability of the porphycene hydride species is that the porphycene ring has a lower LUMO energy level due to the decrease in the symmetry of the ligand character. However, the hydride species in a highly concentrated solution of the complex is easily converted into the cobalt(II) species via dimerization or reaction of the hydride with excess tributyltin hydride through hemolysis of the Sn–H/Co–H bonds. When the Co(III) porphycene is reacted with LiBHEt3 , the final product is the cobalt(III)–ethyl complex formed by β-rearrangement during the reaction of the hydride species and diethylborane in a solvent cage. In the reaction of tetrakistrifluoromethylporphycenatocobalt(III) with tributyltin hydride, the dominant reaction pathway includes one-electron reduction of the porphycene ring together with radical coupling of the tin reagent rather than the net hydride transfer. This finding suggests that the delicate control of the LUMO energy level influences the stability of the hydride species. The tetrapropylporphycenatocobalt(III) complex with tributyltin or triphenyltin hydride in the presence of AIBN produces the corresponding Co(III)– trialkyltin complex. This complex was characterized by 1H NMR spectroscopy.

The long and winding road to new porphycenes

We describe attempts — not always successful — made over the years to improve the efficiency of porphycene synthesis and to produce novel compounds, custom-designed for specific purposes. New porphycenes are reported, some of them obtained rather unexpectedly as by-products of the planned reactions. Structure and energy computations of possible tautomeric forms in porphycenes substituted by one, two, three, and four tert-butyl groups lead to predictions regarding the kinetics and mechanisms of intramolecular double hydrogen transfer. The occurrence of tautomerization in single molecules of tert-butylsubstituted porphycenes is demonstrated by using fluorescence polarization techniques.

Efficient synthesis of porphycene

Unsubstituted porphycene was prepared in 65% yield from 2,7,12,17-tetra-tert-butylporphycene. Taking into account the yield of the substrate, this represents more than a five-fold improvement compared to the methodology used to date. Enhanced availability of the parent porphycene may be exploited in synthetic procedures based on derivatization of the unsubstituted compound.

Unusually Slow Intermolecular Proton-Deuteron Exchange in Porphycene

The rates of the intermolecular exchange of two internal protons by deuterons have been determined for alcohol solutions of porphycene and its 2,7,12,17-tetra-t-butyl derivative using femtosecond pump-probe polarization spectroscopy. The process is very slow, requiring approximately twenty hours to substitute both protons by deuterons in more than 90% of the molecules dissolved in bulk EtOD or BuOD at 293 K. Equal rates were found for the exchange of the first and the second proton.

Fluorescence Studies of Terrylene in a Supersonic Jet: Indication of A Dark Electronic State Below the Allowed Transition

Jet-cooled terrylene has been studied in helium buffer gas using a pulsed nozzle by means of laser-induced fluorescence. Fluorescence excitation and two-color depletion experiments (resulting in hole burning spectra) are presented. Analysis of the spectra leads to the conclusion that another excited electronic state is present in the vicinity of the allowed 1B1u state. Assuming (according to previous literature suggestions Karabunarliev, S.; Baumgarten, M.; Müllen, K. J. Phys. Chem. A 1998, 102, 7029) that this dark state is the 21Ag state, we discuss the vibrational structure of the fluorescence excitation spectrum in terms of two manifolds of vibronic states belonging to Sd(21Ag) and S1(1B1u) states. The anomalous shift between excitation and dispersed fluorescence spectra observed earlier for terrylene in a neon matrix is discussed as a consequence of terrylene electronic relaxation to the low-energy dark state.

Vibrational Gating of Double Hydrogen Tunneling in Porphycene

A procedure that enables determining the reaction rate from the analysis of fluorescence anisotropy is described and applied to the investigation of double hydrogen transfer between inner-cavity nitrogen atoms in electronically excited porphycene. Tautomerization proceeds as a thermally activated synchronous double hydrogen tunneling. The barrier to the reaction is dynamically modulated by a vibration that simultaneously changes the strength of two intramolecular hydrogen bonds. Different mechanisms of tautomerization in porphycene and its parent isomer, porphyrin, can be understood by analyzing the potentials for hydrogen transfer.