1
|
Nanni L. Investigating proton tunneling dynamics in the time-dependent Schrödinger equation. J Comput Chem 2024; 45:1614-1623. [PMID: 38523084 DOI: 10.1002/jcc.27348] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 03/02/2024] [Accepted: 03/07/2024] [Indexed: 03/26/2024]
Abstract
Understanding the temporal evolution of the wave function in the proton tunneling reactions allows us to make theoretical predictions on the possible femtosecond spectroscopy patterns. However, the analytical solution of the time-dependent Schrödinger equation of complex molecules is challenging and requires a high computational cost. In this study, we solve the time-dependent Schrödinger equation using the Fourier grid Hamiltonian method, highlighting its simplicity of calculation, even for multidimensional tunneling reactions. The obtained model is applied to studying malonaldehyde's two-dimensional intramolecular proton tunneling, comparing the results with those obtained using other computational methods.
Collapse
Affiliation(s)
- Luca Nanni
- Department of Physics, University of Ferrara, Ferrara, Italy
| |
Collapse
|
2
|
Kasprzycki P, Kopycki P, Listkowski A, Gorski A, Radzewicz C, Birch DJS, Waluk J, Fita P. Influence of local microenvironment on the double hydrogen transfer in porphycene. Phys Chem Chem Phys 2020; 22:17117-17128. [PMID: 32687131 DOI: 10.1039/d0cp02687e] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Piotr Kasprzycki
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland. and Institute of Physical Chemistry, Polish Academy of Sciences, 44/52 Kasprzaka, Warsaw 01-224, Poland.
| | - Przemysław Kopycki
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland.
| | - Arkadiusz Listkowski
- Institute of Physical Chemistry, Polish Academy of Sciences, 44/52 Kasprzaka, Warsaw 01-224, Poland. and Faculty of Mathematics and Science, Cardinal Stefan Wyszyński University, Dewajtis 5, 01-815 Warsaw, Poland
| | - Aleksander Gorski
- Institute of Physical Chemistry, Polish Academy of Sciences, 44/52 Kasprzaka, Warsaw 01-224, Poland.
| | - Czesław Radzewicz
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland.
| | - David J S Birch
- Photophysics Group, Centre for Molecular Nanometrology, Department of Physics, Scottish Universities Physics Alliance, University of Strathclyde, 107 Rottenrow East, Glasgow G4 0NG, UK
| | - Jacek Waluk
- Institute of Physical Chemistry, Polish Academy of Sciences, 44/52 Kasprzaka, Warsaw 01-224, Poland. and Faculty of Mathematics and Science, Cardinal Stefan Wyszyński University, Dewajtis 5, 01-815 Warsaw, Poland
| | - Piotr Fita
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland.
| |
Collapse
|
3
|
Litman Y, Richardson JO, Kumagai T, Rossi M. Elucidating the Nuclear Quantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene. J Am Chem Soc 2019; 141:2526-2534. [PMID: 30648386 PMCID: PMC6728096 DOI: 10.1021/jacs.8b12471] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
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 × 1011 s–1 at 150 K and 0.63 × 1011 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.
Collapse
Affiliation(s)
- Yair Litman
- Theory Department , Fritz Haber Institute of the Max Planck Society , Faradayweg 4-6 , 14195 Berlin , Germany
| | | | - Takashi Kumagai
- Physical Chemistry Department , Fritz Haber Institute of the Max Planck Society , Faradayweg 4-6 , 14195 Berlin , Germany
| | - Mariana Rossi
- Theory Department , Fritz Haber Institute of the Max Planck Society , Faradayweg 4-6 , 14195 Berlin , Germany
| |
Collapse
|
4
|
Smedarchina Z, Siebrand W, Fernández-Ramos A. Entanglement and co-tunneling of two equivalent protons in hydrogen bond pairs. J Chem Phys 2018; 148:102307. [DOI: 10.1063/1.5000681] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zorka Smedarchina
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), e Departamento de Química Física, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
- National Research Council of Canada, Ottawa, Ontario K1A 0R6 Canada
| | - Willem Siebrand
- National Research Council of Canada, Ottawa, Ontario K1A 0R6 Canada
| | - Antonio Fernández-Ramos
- Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), e Departamento de Química Física, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| |
Collapse
|
5
|
Abstract
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.
Collapse
Affiliation(s)
- Jacek Waluk
- Institute of Physical Chemistry, Polish Academy of Sciences , 01-224 Warsaw, Kasprzaka 44/52, Poland.,Faculty of Mathematics and Science, Cardinal Stefan Wyszyński University , Dewajtis 5, 01-815 Warsaw, Poland
| |
Collapse
|
6
|
Meisner J, Kästner J. Atom Tunneling in Chemistry. Angew Chem Int Ed Engl 2016; 55:5400-13. [DOI: 10.1002/anie.201511028] [Citation(s) in RCA: 140] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2015] [Revised: 01/08/2016] [Indexed: 01/01/2023]
Affiliation(s)
- Jan Meisner
- Institut für Theoretische Chemie Universität Stuttgart Pfaffenwaldring 55 70569 Stuttgart Germany
| | - Johannes Kästner
- Institut für Theoretische Chemie Universität Stuttgart Pfaffenwaldring 55 70569 Stuttgart Germany
| |
Collapse
|
7
|
Affiliation(s)
- Jan Meisner
- Institut für Theoretische Chemie Universität Stuttgart Pfaffenwaldring 55 70569 Stuttgart Deutschland
| | - Johannes Kästner
- Institut für Theoretische Chemie Universität Stuttgart Pfaffenwaldring 55 70569 Stuttgart Deutschland
| |
Collapse
|
8
|
Alves TV, Simón-Carballido L, Ornellas FR, Fernández-Ramos A. Hindered rotor tunneling splittings: an application of the two-dimensional non-separable method to benzyl alcohol and two of its fluorine derivatives. Phys Chem Chem Phys 2016; 18:8945-53. [PMID: 26960818 DOI: 10.1039/c5cp05307b] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work we present a novel application of the two-dimensional non-separable (2D-NS) method to the calculation of torsional tunneling splittings in systems with two hindered internal rotors. This method could be considered an extension of one-dimensional methods for the case of compounds with two tops. The 2D-NS method includes coupling between torsions in the kinetic and potential energy. Specifically, it has been applied to benzyl alcohol (BA) and two of its fluorine derivatives: 3-fluorobenzyl alcohol (3FBA) and 4-fluorobenzyl alcohol (4FBA). These molecules present two torsions, i.e., about the -CH2OH (ϕ1) and -OH (ϕ2) groups. The electronic structure calculations to build the two-dimensional torsional potential energy surface were performed at the DF-LMP2-F12//DF-LMP2/cc-pVQZ level of theory. For BA and 4FBA the calculated ground-state vibrational level splittings are 429 and 453 MHz, respectively, in good agreement with the experimental values of 337.10 and 492.82 MHz, respectively. In these two cases there are four equivalent wells and the tunneling splitting is the result of transitions between the two closer minima along ϕ1. The analysis of the wavefunctions, as well as the previous experimental work on the system, supports this conclusion. For 3FBA the observed ground-state splitting is 0.82 MHz, whereas in this case the calculated value amounts only to 0.02 MHz. The 2D-NS method, through the analysis of the wavefunctions, shows that this tiny tunneling splitting occurs between the two most stable minima of the potential energy surface. Additionally, we predict that the first vibrationally excited tunneling splitting will also be small and exclusively due to the interconversion between the second lowest minima.
Collapse
Affiliation(s)
- Tiago Vinicius Alves
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, São Paulo, SP 05508-000, Brazil
| | | | | | | |
Collapse
|
9
|
Ciąćka P, Fita P, Listkowski A, Radzewicz C, Waluk J. Evidence for Dominant Role of Tunneling in Condensed Phases and at High Temperatures: Double Hydrogen Transfer in Porphycenes. J Phys Chem Lett 2016; 7:283-288. [PMID: 26727277 DOI: 10.1021/acs.jpclett.5b02482] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Piotr Ciąćka
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw , Pasteura 5, 02-093 Warsaw, Poland
| | - Piotr Fita
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw , Pasteura 5, 02-093 Warsaw, Poland
| | - Arkadiusz Listkowski
- Institute of Physical Chemistry, Polish Academy of Sciences , Kasprzaka 44/52, 01-224 Warsaw, Poland
- Faculty of Mathematics and Natural Sciences, College of Science, Cardinal Stefan Wyszyński University , Dewajtis 5, 01-815 Warsaw, Poland
| | - Czesław Radzewicz
- Institute of Experimental Physics, Faculty of Physics, University of Warsaw , Pasteura 5, 02-093 Warsaw, Poland
| | - Jacek Waluk
- Institute of Physical Chemistry, Polish Academy of Sciences , Kasprzaka 44/52, 01-224 Warsaw, Poland
- Faculty of Mathematics and Natural Sciences, College of Science, Cardinal Stefan Wyszyński University , Dewajtis 5, 01-815 Warsaw, Poland
| |
Collapse
|
10
|
Ren Y, Bian W. Mode-Specific Tunneling Splittings for a Sequential Double-Hydrogen Transfer Case: An Accurate Quantum Mechanical Scheme. J Phys Chem Lett 2015; 6:1824-1829. [PMID: 26263255 DOI: 10.1021/acs.jpclett.5b00672] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We present the first accurate quantum dynamics calculations of mode-specific tunneling splittings in a sequential double-hydrogen transfer process. This is achieved in the vinylidene-acetylene system, the simplest molecular system of this kind, and by large-scale parallel computations with an efficient theoretical scheme developed by us. In our scheme, basis functions are customized for the hydrogen transfer process; a 4-dimensional basis contraction strategy is combined with the preconditioned inexact spectral transform method; efficient parallel implementation is achieved. Mode-specific permutation tunneling splittings of vinylidene states are reported and tremendous mode-specific promotion effects are revealed; in particular, the CH2 rock mode enhances the ground-state splitting by a factor of 10(3). We find that the ground-state vinylidene has a reversible-isomerization time of 622 ps, much longer than all previous estimates. Our calculations also shed light on the importance of the deep intermediate well and vibrational excitation in the double-hydrogen transfer processes.
Collapse
Affiliation(s)
- Yinghui Ren
- †Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- ‡University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wensheng Bian
- †Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|