Abstract
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.
Introduction
The transfer of a proton or hydrogen, often referred to as “the simplest chemical reaction”, is, in fact, a very complicated process. The reason is the quantum nature of proton and its propensity to tunnel. Since tunneling is extremely sensitive to the potential energy barrier “seen” by the proton, and this barrier changes as the atoms vibrate, accurate description of tautomerization is a formidable task. Our work on systems engaged in intra- or intermolecular hydrogen bonds focuses on the understanding of tautomerization mechanisms, which is mandatory for the precise description of the reaction path. I will present several examples of “unusual” proton transfers, which should no longer be called unusual once the mechanisms responsible for tautomerization are elucidated. The second part will focus on the outstanding features of tautomerism in porphycenes. Finally, studies of tautomerization in single molecules will be presented.
Example 1. Mode-selective excited state proton transfer in 2-(2′-pyridyl)pyrrole. This molecule exhibits excited state intramolecular proton transfer (ESIPT), manifested by the appearance of dual fluorescence: F1, corresponding to the initially excited species, and F2, originating from the tautomer formed in S1 after proton translocation from the pyrrole to pyridyl nitrogen (Fig. 1). Dual emission is observed not only in solutions [1], but also for molecules isolated in supersonic jets [2], where the reaction becomes vibrational-mode specific: Tautomerization occurs after selective excitation into the vibronic levels corresponding to two in-plane modes, those that lead to the decrease of the N···N distance (and thus, to the increase of hydrogen bond strength). Excitation of other modes does not lead to the reaction, even though the energy deposited in the molecule may be larger than in the previous case. This finding immediately shows the danger involved in using 1D reaction models, which predict that when more energy is put into the system, it gets closer to the barrier, and the reaction is accelerated. In reality, vibrational excitation does not always help. Depending on the form of a particular vibrational mode, tautomerization may be enhanced, slowed down, or not influenced by vibrational excitation.

Vibrational mode-specific excited state intramolecular proton transfer in 2-(2′-pyridyl)pyrrole isolated in supersonic jet. Depending on which vibration is excited the molecule can either (a) remain in the initially excited form and emit “normal” fluorescence (F1) or (b) tunnel to the tautomeric form, identified by the red-shifted fluorescence (F2).
Kinetic isotope effect of 30–60, proving the tunneling mechanism, were obtained from the comparison of F1/F2 intensity ratios in non-deuterated and N-deuterated 2-(2′-pyridyl)pyrrole [2].
Example 2. Coupling between excited state proton transfer and a large-amplitude twisting motion. Photoexcited syn rotamer of 7-(2′-pyridyl)indole decays very rapidly: The S1 lifetime in solution at 293 K is ca 1 ps. Remarkably, in the supersonic jet regime the decay becomes faster than in room temperature solution (280 fs and 390 fs in N-deuterated molecule). This was interpreted [3] as evidence of coupling of proton transfer with twisting of the pyridyl group with respect to the indole moiety (Fig. 2). Theory predicts the proton transfer to be barrierless and coupled to the indole-pyridine torsional coordinate [4]. The apparently unusual increase of the reaction rate at lower temperatures has been explained by the presence of viscosity barrier in solution and its absence in jet-isolated molecules.

Excited state intramolecular proton transfer in 7-(2′-pyridyl)indole is accompanied by twisting. Proper description of the reaction path requires using at least two coordinates, denoted here as rNH and φ.
Example 3. Double hydrogen transfer is faster than a single hydrogen transfer reaction. Coupling of tautomerization with torsional motions of methyl groups. Meso-tetramethyl-substituted porphycene (Fig. 3) exhibits ultrafast (<100 fs) trans-trans double hydrogen transfer, whereas the translocation of a single hydrogen atom trans-cis-trans is much slower [5]. This has been explained by coupling of single hydrogen transfer with the torsional motion of methyl substituents, and the lack of such coupling in the case of the double hydrogen transfer process. In other words, the reaction coordinates are very different for single and double hydrogen transfer processes.

Trans-trans interconversion in 9,10,19,20-tetramethylporphycene is ultrafast, but the cis-trans and trans-cis processes are much slower.
Example 4. Tautomerization rate in a single molecule becomes temporarily slower by many orders of magnitude. Trans-trans conversion has been studied for 2,7,12,17-tetraphenylporphycene on a single molecule level using fluorescence microscopy [6]. The analysis of the spatial fluorescence patterns of more than 850 single molecules, revealed, for a small fraction (about 5%), an unusual behavior: change ot tautomerization rate with time. The average tautomerization time, determined from measurements on bulk samples, is about 1 ps. In contrast, in some of the single molecules, this time was extended to minutes. The huge decrease of the tautomerization rate was explained by the inability of the molecule to assume a conformation corresponding to a symmetric double minimum potential (Fig. 4). Such behavior may be due to the low frequency torsional motions of phenyl substituents, coupled to the motions of the polymer matrix.

Bottom, confocal fluorescence microscopy images of fluorescence of single molecules of 2,7,12,17-tetraphenylporphycene embedded in a polymer film at 293 K. Images (a) and (b) correspond to two consecutive 20 min scans. The molecule indicated by the circle switches from a “locked” state (a) to a regime of fast tautomerization (b). Top, schematic representation of the potential energy profiles for the two cases.
Outstanding features of tautomerism in porphycenes
Due to strong intramolecular hydrogen bonds and well defined geometry of the inner cavity, porphycenes are well-suited for studies of non-classical effects in tautomerization. They are also interesting for investigations in other areas, such as photodynamic therapy, photostability, or possible application as dyes.
Electronic spectra recorded for cold porphycene, isolated in supersonic jets or helium nanodroplets reveal tunneling splittings [7]. The magnitude of the splittings is much larger in S0 than in S1: 4.4 vs<0.2 cm−1, respectively, for the 0–0 transition. Single vibronic level fluorescence spectra (Fig. 5) show that the tunneling splittings are vibrational mode-specific [8], [9]. In other words, the barrier to tautomerization may change with vibrational excitation. Both an increase (12 cm−1 for the 2Ag mode) and decrease (to a value <1 cm−1 for 1Ag) have been observed. Most bands reveal splittings similar to that observed for the vibrationless level. However, a “neutral” mode may change its character when the molecule is isotopically substituted, even at positions remote from the inner cavity. Such a case was observed for the 4Ag mode: in the nondeuterated porphycene it exhibits practically the same splitting as the zero point level, but in the molecule with all 12 peripheral protons replaced by deuterons, this value increases to 9 cm−1 [10]. Such a “reversed isotope effect” clearly demonstrates the multidimensional character of the potential responsible for hydrogen transfer. It is possible to reproduce the increased/decreased splittings by using the models of 2D symmetric or antisymmetric coupling between the high frequency NH stretching vibration and a low frequency mode [11]. In order to simulate the case when both tautomerization enhancing and hindering modes are simultaneously excited, symmetric and antisymmetric mode coupling can be combined. These results demonstrate that the correct potential should be multidimensional.

Single vibronic level fluorescence spectra of porphycene isolated in a supersonic jet. Spectra (a), (b), and (d) correspond to the emission obtained by exciting the 2Ag vibrational mode in the S1 state (v′=1, 2, and 3, respectively), whereas the spectrum (c) was obtained for 4Ag excitation. Ground state vibrational levels are denoted by v″.
In a condensed phase (solutions, glasses, polymers), the two inner protons are no longer coherently delocalized. Tautomerization now becomes a rate process. Since the trans-trans conversion is a self-exchange process, i.e. the substrate and product are chemically equivalent, special techniques are required for rate determination. They are based on pump-probe spectroscopy with polarized light [12], and the reaction rate can be determined from the analysis of transient absorption anisotropy (r(t)) signals:
ΔA(t)par and ΔA(t)perp denote changes in absorbance measured for probe polarization parallel and perpendicular to the polarization of pump pulse, respectively. Examples of transient anisotropy curves are provided in Fig. 6. The measured signal may contain contributions from (a) ground state depopulation by the pump pulse (bleaching); (b) stimulated emission; (c) transient absorption originating from the excited state. A proper choice of pump and probe wavelengths enables the separation of these contributions and the simultaneous determination of ground

Kinetic profiles of the anisotropy r(t) measured, after excitation to the S1 state, for four different porphycenes in paraffin oil solution at room temperature.
where f is the relative contribution from S0,
This procedure has now been applied for more than 20 differently substituted porphycenes [13], [14], [15]. The rates span more than three orders of magnitude and correlate well with the parameters that characterize the strength of intramolecular hydrogen bonds: the NH stretching frequencies, N-N distances, proton chemical shifts. The large rate variations are indicative of tunneling. In order to determine its role in tautomerization, parent porphycene and 2,7,12,17-tetra-tert-butyl porphycene, together with their –d2 isotopologues (inner protons substituted with deuterons) have been studied in a wide temperature range, 20–500 K [16]. It was found that the rate can be described assuming contributions from three different channels, of which two correspond to (i) tunneling from a vibrational ground state; (ii) vibrationally-activated tunneling (via the 2Ag mode). Even at room temperature, the role of these two channels is dominant, showing that tunneling may be the major reaction path even under “normal” conditions.
It is interesting to check if the data obtained for the two regimes – isolated molecules and condensed phases – are compatible. For incoherent tunneling, the rate is proportional to Δ2. Using the experimental data obtained for porphycene in the vibrationless state and for the v=1 level of the 2Ag mode, one obtains: (12/4.4)2=7.4. This is in perfect agreement with the corresponding ratio of the rates determined for the two tunneling channels: k1/k0=7.0±0.9.
Tautomerism in single molecules
It has now become possible to observe tautomerism in single molecules. For porphycene, this has been achieved using three different techniques. The first was fluorescence microscopy [17], [18] (Fig. 7). Spatial images of emission from single chromophores are recorded while using polarized light for excitation (the so-called radial and azimuthal polarization modes). Certain image shapes, such as a ring pattern, observed for molecule Σ in Fig. 7, cannot be obtained when the emission is due to a single oscillating dipole, but they can be simulated assuming two transition moments forming a certain angle (ca. 72° for porphycene). This corresponds to two trans tautomers rapidly interconverting on the time scale of the experiment (usually seconds). Investigation of single molecules enables observing events that are not detected while studying ensembles, and thus averaging over large number of molecules. Example 4, presented above (Fig. 4), describes such a situation: a small fraction of single molecules of 2,7,12,17-tetraphenylporphycene was found to exhibit an unusual behavior: switching between the regime of slow and fast tautomerization over a period of minutes.
![Fig. 7:
(a) Confocal image of two porphycene molecules showing two typical fluorescence patterns: a ring pattern for molecule Σ and a double-lobe pattern for molecule Π. (b) 3D plot of the same area as in (a). (c), (d) Simulated fluorescence patterns, assuming two differently oriented molecules with two transition dipole moments forming an angle of 72°. Scan area: 2.5×2.5 μm2). Reprinted with permission from Ref. [17]. Copyright (2005) American Chemical Society.](/document/doi/10.1515/pac-2016-0706/asset/graphic/j_pac-2016-0706_fig_017.jpg)
(a) Confocal image of two porphycene molecules showing two typical fluorescence patterns: a ring pattern for molecule Σ and a double-lobe pattern for molecule Π. (b) 3D plot of the same area as in (a). (c), (d) Simulated fluorescence patterns, assuming two differently oriented molecules with two transition dipole moments forming an angle of 72°. Scan area: 2.5×2.5 μm2). Reprinted with permission from Ref. [17]. Copyright (2005) American Chemical Society.
Not only fluorescence, but also Raman spectra of single molecules can now be obtained. Surface-enhanced resonance Raman spectra were recorded for single porphycene molecules located on gold nanostructures [19], [20]. Thousands of spectra could be registered from the same molecule at a rate of several spectra per second. Analysis of the consecutive spectra revealed infrequent, reversible changes: disappearance of some bands and the appearance of new ones. Such behavior was interpreted as due to trans-cis-trans tautomerization in single porphycene molecules. This conclusion was supported by calculations of vibrational structure.
Unlike optical methods that monitor the reaction via fluorescence or Raman spectra, scanning probe microscopies enable looking “directly” at protons in single, tautomerizing molecules (Fig. 8). Scanning tunneling microscopy experiments were performed for porphycene located on Cu(110) surface, under high vacuum at 5 K [21], [22], [23], [24], [25]. Under such conditions, the molecule was found to exist in the cis form, a structure not observed previously for parent porphycene. Cis-cis reaction rates could be determined from the analysis of time intervals observed between tautomerization events. Tautomerization could be induced by heating, but also by inelastic electron scattering, leading to vibrational excitation of porphycene. The reaction activated in this way was found to be vibrational mode-specific [22]. Further STM/AFM experiments demonstrated that tautomerization in single porphycene molecules on copper surface can also be induced by two other stimuli: (a) light [24]; (b) mechanical force [25]. An impressive result was the possibility to switch between two cis tautomers by moving a single copper atom in the vicinity of a single porphycene molecule [21].
![Fig. 8:
Top, STM images of a single porphycene molecule on a Cu(110) surface at 5 K, demonstrating the existence of two cis tautomeric forms. Bottom, current traces obtained for a fixed tip position, illustrating switching between the two tautomers. Adapted from Ref. [21].](/document/doi/10.1515/pac-2016-0706/asset/graphic/j_pac-2016-0706_fig_018.jpg)
Top, STM images of a single porphycene molecule on a Cu(110) surface at 5 K, demonstrating the existence of two cis tautomeric forms. Bottom, current traces obtained for a fixed tip position, illustrating switching between the two tautomers. Adapted from Ref. [21].
Conclusions
The experiments described above leave no doubt that the proper explanation of proton/hydrogen transfer reactions cannot neglect tunneling and the contributions from different vibrational modes to the reaction path. Thus, accurate description of tautomerism requires a multidimensional approach with explicit consideration of tunneling.
Studies of meso-substituted porphycenes reveal that proton transfer can be coupled to conformational changes, such as methyl group rotation. Another example is provided by photoinduced proton transfer in 7-(2′-pyridyl)indole.
Tunneling splittings due to coherent delocalization of two internal protons in porphycene are vibrational mode-specific. In other words, the barrier to hydrogen transfer may strongly depend on the vibrational state of the molecule. This was also demonstrated for the electronically excited 2-(2′-pyridyl)pyrrole.
A general conclusion that emerges from extensive studies of porphycenes, carried out under diverse experimental regimes, such as condensed phases, isolated or single molecules is that the quantum effects, tunneling in particular, are important, and perhaps even dominant, even for “normal”, routinely used molecular environments: room temperature solutions. Non-classical effects in proton/hydrogen transfer should be always taken into account while trying to understand the intricacies of tautomerism.
Article note
A collection of invited papers based on presentations at the XXVIth IUPAC Symposium on Photochemistry, Osaka, Japan, April 3–8 2016.
Acknowledgments
This work was supported by the Polish National Science Centre (NCN) grants no. DEC-2011/02/A/ST5/00443, DEC-2013/10/M/ST4/00069, by the PL-Grid Infrastructure grant, as well as the computing grant from the Interdisciplinary Centre for Mathematical and Computational Modeling.
References
[1] M. Kijak, A. Zielińska, C. Chamchoumis, J. Herbich, R. P. Thummel, J. Waluk. Chem. Phys. Lett.400, 279 (2004).10.1016/j.cplett.2004.10.121Search in Google Scholar
[2] M. Kijak, E. Nosenko, A. Singh, R. P. Thummel, J. Waluk. J. Am. Chem. Soc. 129, 2738 (2007).10.1021/ja068109fSearch in Google Scholar
[3] E. Nosenko, G. Wiosna-Sałyga, M. Kunitski, I. Petkova, A. Singh, W. J. Buma, R. P. Thummel, J. Waluk. Angew. Chem. Int. Ed. 47, 6037 (2008).10.1002/anie.200801350Search in Google Scholar
[4] A. L. Sobolewski, W. Domcke. J. Phys. Chem. A111, 11725 (2007).10.1021/jp075803oSearch in Google Scholar
[5] M. Gil, J. Dobkowski, G. Wiosna-Sałyga, N. Urbańska, P. Fita, C. Radzewicz, M. Pietraszkiewicz, P. Borowicz, D. Marks, M. Glasbeek, J. Waluk. J. Am. Chem. Soc.132, 13472 (2010).10.1021/ja105353mSearch in Google Scholar
[6] H. Piwoński, A. Sokołowski, M. Kijak, S. Nonell, J. Waluk. J. Phys. Chem. Lett.4, 3967 (2013).10.1021/jz4022602Search in Google Scholar
[7] J. Sepioł, Y. Stepanenko, A. Vdovin, A. Mordziński, E. Vogel, J. Waluk. Chem. Phys. Lett. 296, 549 (1998).10.1016/S0009-2614(98)01085-9Search in Google Scholar
[8] A. Vdovin, J. Waluk, B. Dick, A. Slenczka. ChemPhysChem10, 761 (2009).10.1002/cphc.200900022Search in Google Scholar
[9] E. T. Mengesha, J. Sepioł, P. Borowicz, J. Waluk. J. Chem. Phys.138, 174201 (2013).10.1063/1.4802769Search in Google Scholar
[10] E. T. Mengesha, A. Zehnacker-Rentien, J. Sepioł, M. Kijak, J. Waluk. J. Phys. Chem. B119, 2193 (2015).10.1021/jp505553zSearch in Google Scholar
[11] S. Takada, H. Nakamura. J. Chem Phys. 102, 3977 (1995).10.1063/1.468526Search in Google Scholar
[12] P. Fita, N. Urbańska, C. Radzewicz, J. Waluk. Chem. Eur. J. 15, 4851 (2009).10.1002/chem.200802428Search in Google Scholar
[13] P. Ciąćka, P. Fita, A. Listkowski, M. Kijak, S. Nonell, D. Kuzuhara, H. Yamada, C. Radzewicz, J. Waluk. J. Phys. Chem. B119, 2292 (2015).10.1021/jp506150rSearch in Google Scholar
[14] P. Fita, P. Garbacz, M. Nejbauer, C. Radzewicz, J. Waluk. Chem. Eur. J. 17, 3672 (2011).Search in Google Scholar
[15] P. Fita, M. Pszona, G. Orzanowska, D. Sánchez-García, S. Nonell, E. Vauthey, J. Waluk. Chem.- Eur. J.18, 13160 (2012).Search in Google Scholar
[16] P. Ciąćka, P. Fita, A. Listkowski, C. Radzewicz, J. Waluk. J. Phys. Chem. Lett.7, 283 (2016).10.1021/acs.jpclett.5b02482Search in Google Scholar
[17] H. Piwoński, C. Stupperich, A. Hartschuh, J. Sepioł, A. Meixner, J. Waluk. J. Am. Chem. Soc.127, 5302 (2005).10.1021/ja043265cSearch in Google Scholar
[18] A. M. Chizhik, R. Jäger, A. I. Chizhik, S. Bär, H. G. Mack, M. Sackrow, C. Stanciu, A. Lyubimtsev, M. Hanack, A. J. Meixner. Phys. Chem. Chem. Phys.13, 1722 (2011).10.1039/c0cp02228dSearch in Google Scholar
[19] S. Gawinkowski, M. Pszona, A. Gorski, J. Niedziółka-Jönsson, I. Kamińska, W. Nogala, J. Waluk. Nanoscale8, 3337 (2016).10.1039/C5NR08627BSearch in Google Scholar
[20] W. Nogala, P. Kannan, S. Gawinkowski, M. Jönsson-Niedziolka, M. Kominiak, J. Waluk, M. Opallo. Nanoscale7, 10767 (2015).10.1039/C5NR02077HSearch in Google Scholar
[21] T. Kumagai, F. Hanke, S. Gawinkowski, J. Sharp, K. Kotsis, J. Waluk, M. Persson, L. Grill. Nature Chem.6, 41 (2014).10.1038/nchem.1804Search in Google Scholar
[22] T. Kumagai, F. Hanke, S. Gawinkowski, J. Sharp, K. Kotsis, J. Waluk, M. Persson, L. Grill. Phys. Rev. Lett.111, 246101 (2013).10.1103/PhysRevLett.111.246101Search in Google Scholar
[23] J. N. Ladenthin, L. Grill, S. Gawinkowski, S. Liu, J. Waluk, T. Kumagai. ACS Nano9, 7287 (2015).10.1021/acsnano.5b02147Search in Google Scholar
[24] H. Böckmann, S. Liu, J. Mielke, S. Gawinkowski, J. Waluk, L. Grill, M. Wolf, T. Kumagai. Nano Lett. 16, 1034 (2016).10.1021/acs.nanolett.5b04092Search in Google Scholar
[25] J. Ladenthin, T. Frederiksen, M. Persson, J. Sharp, S. Gawinkowski, J. Waluk, T. Kumagai. Nature Chem. DOI: 10.1038/NCHEM.2552 (2016).Search in Google Scholar
©2016 IUPAC & De Gruyter. This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License. For more information, please visit: http://creativecommons.org/licenses/by-nc-nd/4.0/
Articles in the same Issue
- Frontmatter
- In this issue
- Editorial
- The 2016 Nobel Prize in Chemistry
- Conference papers
- Cyclotriphosphazene, an old compound applied to the synthesis of smart dendrimers with tailored properties
- Milestones in microwave-assisted organophosphorus chemistry
- Synthetic routes to polyphosphoesters as solid polymer electrolytes for lithium ion batteries
- Soluble polymer supports for homogeneous catalysis in flow reactions
- Box-Behnken experimental design for the production of precipitated calcium carbonate
- Synthesis and characterisation of lignin-like oligomers as a bio-inspired consolidant for waterlogged archaeological wood
- Synthesis, characterizations and Pb(II) sorption properties of cobalt phosphonate materials
- Performance of poly(styrene-co-divinylbenzene) functionalized with different aminophosphonate pendant groups, in the removal of phenolic compounds from aqueous solutions
- Synthesis of new dithia[3.3]parapara- and metapara-cyclophane based tectons: toward an universal surface-confined 2D/3D molecular binding motif
- Effects of surface modification and ultrasonic agitation on the properties of PHBV/ZnO nanocomposites
- Modified macromolecules in the prevention of silica scale
- Ionizing radiation: a versatile tool for nanostructuring of polymers
- Non-classical effects in proton or hydrogen transfer
- IUPAC Recommendations
- Source-based nomenclature for single-strand homopolymers and copolymers (IUPAC Recommendations 2016)
Articles in the same Issue
- Frontmatter
- In this issue
- Editorial
- The 2016 Nobel Prize in Chemistry
- Conference papers
- Cyclotriphosphazene, an old compound applied to the synthesis of smart dendrimers with tailored properties
- Milestones in microwave-assisted organophosphorus chemistry
- Synthetic routes to polyphosphoesters as solid polymer electrolytes for lithium ion batteries
- Soluble polymer supports for homogeneous catalysis in flow reactions
- Box-Behnken experimental design for the production of precipitated calcium carbonate
- Synthesis and characterisation of lignin-like oligomers as a bio-inspired consolidant for waterlogged archaeological wood
- Synthesis, characterizations and Pb(II) sorption properties of cobalt phosphonate materials
- Performance of poly(styrene-co-divinylbenzene) functionalized with different aminophosphonate pendant groups, in the removal of phenolic compounds from aqueous solutions
- Synthesis of new dithia[3.3]parapara- and metapara-cyclophane based tectons: toward an universal surface-confined 2D/3D molecular binding motif
- Effects of surface modification and ultrasonic agitation on the properties of PHBV/ZnO nanocomposites
- Modified macromolecules in the prevention of silica scale
- Ionizing radiation: a versatile tool for nanostructuring of polymers
- Non-classical effects in proton or hydrogen transfer
- IUPAC Recommendations
- Source-based nomenclature for single-strand homopolymers and copolymers (IUPAC Recommendations 2016)