1
|
Zhang ZT, Vaníček JJL. Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics. J Chem Phys 2024; 160:084103. [PMID: 38385512 DOI: 10.1063/5.0187823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Accepted: 12/22/2023] [Indexed: 02/23/2024] Open
Abstract
We present a numerically exact approach for evaluating vibrationally resolved electronic spectra at finite temperatures using the coherence thermofield dynamics. In this method, which avoids implementing an algorithm for solving the von Neumann equation for coherence, the thermal vibrational ensemble is first mapped to a pure-state wavepacket in an augmented space, and this wavepacket is then propagated by solving the standard, zero-temperature Schrödinger equation with the split-operator Fourier method. We show that the finite-temperature spectra obtained with the coherence thermofield dynamics in a Morse potential agree exactly with those computed by Boltzmann-averaging the spectra of individual vibrational levels. Because the split-operator thermofield dynamics on a full tensor-product grid is restricted to low-dimensional systems, we briefly discuss how the accessible dimensionality can be increased by various techniques developed for the zero-temperature split-operator Fourier method.
Collapse
Affiliation(s)
- Zhan Tong Zhang
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří J L Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
2
|
Moghaddasi Fereidani R, Vaníček JJL. High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics. J Chem Phys 2024; 160:044113. [PMID: 38284658 DOI: 10.1063/5.0180070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Accepted: 12/22/2023] [Indexed: 01/30/2024] Open
Abstract
Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.
Collapse
Affiliation(s)
- Roya Moghaddasi Fereidani
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří J L Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
3
|
Wenzel M, Mitric R. Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation. J Chem Phys 2023; 159:234113. [PMID: 38108487 DOI: 10.1063/5.0178106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 11/26/2023] [Indexed: 12/19/2023] Open
Abstract
Spontaneous emission and internal conversion rates are calculated within harmonic approximations and compared to the results obtained within the semi-classical extended thawed Gaussian approximation (ETGA). This is the first application of the ETGA in the calculation of internal conversion and emission rates for real molecular systems, namely, formaldehyde, fluorobenzene, azulene, and a dicyano-squaraine dye. The viability of the models as black-box tools for prediction of spontaneous emission and internal conversion rates is assessed. All calculations were done using a consistent protocol in order to investigate how different methods perform without previous experimental knowledge using density functional theory (DFT) and time-dependent DFT (TD-DFT) with B3LYP, PBE0, ωB97XD, and CAM-B3LYP functionals. Contrasting the results with experimental data shows that there are further improvements required before theoretical predictions of emission and internal conversion rates can be used as reliable indicators for the photo-luminescence properties of molecules. We find that the ETGA performs rather similar to the vertical harmonical model. Including anharmonicities in the calculation of internal conversion rates has a moderate effect on the quantitative results in the studied systems. The emission rates are fairly stable with respect to computational parameters, but the internal conversion rate reveals itself to be highly dependent on the choice of the spectral line shape function, particularly the width of the Lorentzian function, associated with homogeneous broadening.
Collapse
Affiliation(s)
- Michael Wenzel
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| | - Roland Mitric
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
4
|
Franzke Y, Holzer C, Andersen JH, Begušić T, Bruder F, Coriani S, Della Sala F, Fabiano E, Fedotov DA, Fürst S, Gillhuber S, Grotjahn R, Kaupp M, Kehry M, Krstić M, Mack F, Majumdar S, Nguyen BD, Parker SM, Pauly F, Pausch A, Perlt E, Phun GS, Rajabi A, Rappoport D, Samal B, Schrader T, Sharma M, Tapavicza E, Treß RS, Voora V, Wodyński A, Yu JM, Zerulla B, Furche F, Hättig C, Sierka M, Tew DP, Weigend F. TURBOMOLE: Today and Tomorrow. J Chem Theory Comput 2023; 19:6859-6890. [PMID: 37382508 PMCID: PMC10601488 DOI: 10.1021/acs.jctc.3c00347] [Citation(s) in RCA: 37] [Impact Index Per Article: 37.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Indexed: 06/30/2023]
Abstract
TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.
Collapse
Affiliation(s)
- Yannick
J. Franzke
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Christof Holzer
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Josefine H. Andersen
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Tomislav Begušić
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Florian Bruder
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Sonia Coriani
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Fabio Della Sala
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Eduardo Fabiano
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Daniil A. Fedotov
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Susanne Fürst
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Sebastian Gillhuber
- Institute
of Inorganic Chemistry, Karlsruhe Institute
of Technology (KIT), Engesserstr. 15, 76131 Karlsruhe, Germany
| | - Robin Grotjahn
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Max Kehry
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Marjan Krstić
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Sourav Majumdar
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Brian D. Nguyen
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Shane M. Parker
- Department
of Chemistry, Case Western Reserve University, 10900 Euclid Ave, Cleveland, Ohio 44106 United States
| | - Fabian Pauly
- Institute
of Physics, University of Augsburg, Universitätsstr. 1, 86159 Augsburg, Germany
| | - Ansgar Pausch
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Eva Perlt
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Gabriel S. Phun
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Ahmadreza Rajabi
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Dmitrij Rappoport
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Bibek Samal
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Tim Schrader
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Manas Sharma
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Enrico Tapavicza
- Department
of Chemistry and Biochemistry, California
State University, Long Beach, 1250 Bellflower Boulevard, Long
Beach, California 90840-9507, United States
| | - Robert S. Treß
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Vamsee Voora
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Artur Wodyński
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Jason M. Yu
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Benedikt Zerulla
- Institute
of Nanotechnology, Karlsruhe Institute of
Technology (KIT), Hermann-von-Helmholtz-Platz
1, 76344 Eggenstein-Leopoldshafen Germany
| | - Filipp Furche
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Christof Hättig
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Marek Sierka
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - David P. Tew
- Physical
and Theoretical Chemistry Laboratory, University
of Oxford, South Parks
Road, Oxford OX1 3QZ, United Kingdom
| | - Florian Weigend
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| |
Collapse
|
5
|
Klētnieks Ē, Alonso YC, Vaníček JJL. Isotope Effects on the Electronic Spectra of Ammonia from Ab Initio Semiclassical Dynamics. J Phys Chem A 2023; 127:8117-8125. [PMID: 37738161 PMCID: PMC10561269 DOI: 10.1021/acs.jpca.3c04607] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2023] [Revised: 08/22/2023] [Indexed: 09/24/2023]
Abstract
Despite its simplicity, the single-trajectory thawed Gaussian approximation has proven useful for calculating the vibrationally resolved electronic spectra of molecules with weakly anharmonic potential energy surfaces. Here, we show that the thawed Gaussian approximation can capture surprisingly well even more subtle observables, such as the isotope effects in the absorption spectra, and we demonstrate it on the four isotopologues of ammonia (NH3, NDH2, ND2H, and ND3). The differences in their computed spectra are due to the differences in the semiclassical trajectories followed by the four isotopologues, and the isotope effects─narrowing of the transition band and reduction of the peak spacing─are accurately described by this semiclassical method. In contrast, the adiabatic harmonic model shows a double progression instead of the single progression seen in the experimental spectra. The vertical harmonic model correctly shows only a single progression but fails to describe the anharmonic peak spacing. Analysis of the normal-mode activation upon excitation provides insight into the elusiveness of the symmetric stretching progression in the spectra.
Collapse
Affiliation(s)
- Ēriks Klētnieks
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| | - Yannick Calvino Alonso
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| | - Jiří J. L. Vaníček
- Laboratory of Theoretical Physical
Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne
(EPFL), Lausanne CH-1015, Switzerland
| |
Collapse
|
6
|
Szychta K, Koszarna B, Banasiewicz M, Sobolewski A, O’Mari O, Clark JA, Vullev VI, Barboza CA, Gryko DT. Conformation of the Ester Group Governs the Photophysics of Highly Polarized Benzo[ g]coumarins. JACS AU 2023; 3:1918-1930. [PMID: 37502148 PMCID: PMC10369411 DOI: 10.1021/jacsau.3c00169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 06/10/2023] [Accepted: 06/13/2023] [Indexed: 07/29/2023]
Abstract
Photosensitizers that display "unusual" emission from upper electronically excited states offer possibilities for initiating higher-energy processes than what the governing Kasha's rule postulates. Achieving conditions for dual fluorescence from multiple states of the same species requires molecular design and conditions that favorably tune the excited-state dynamics. Herein, we switch the position of the electron-donating NMe2 group around the core of benzo[g]coumarins (BgCoum) and tune the electronic coupling and the charge-transfer character of the fluorescent excited states. For solvents with intermediate polarity, three of the four regioisomers exhibit fluorescence from two different excited states with bands that are well separated in the visible and the near-infrared spectral regions. Computational analysis, employing ab initio methods, reveals that the orientation of an ester on the pyrone ring produces two conformers responsible for the observed dual fluorescence. Studies with solid solvating media, which restricts the conformational degrees of freedom, concur with the computational findings. These results demonstrate how "seemingly inconsequential" auxiliary substituents, such as the esters on the pyrone coumarin rings, can have profound effects leading to "anti-Kasha" photophysical behavior important for molecular photonics, materials engineering, and solar-energy science.
Collapse
Affiliation(s)
- Kamil Szychta
- Institute
of Organic Chemistry of Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Beata Koszarna
- Institute
of Organic Chemistry of Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Marzena Banasiewicz
- Institute
of Physics of Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Andrzej Sobolewski
- Institute
of Physics of Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Omar O’Mari
- Department
of Bioengineering, University of California, Riverside, California 92521, United States
| | - John A. Clark
- Department
of Bioengineering, University of California, Riverside, California 92521, United States
| | - Valentine I. Vullev
- Department
of Bioengineering, University of California, Riverside, California 92521, United States
- Department
of Chemistry, Department of Biochemistry, and Materials Science and
Engineering Program, University of California, Riverside, California 92521, United States
| | - Cristina A. Barboza
- Institute
of Physics of Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
- Department
of Physical and Quantum Chemistry, Faculty of Chemistry, Wroclaw University of Science and Technology, Wrocław 50-370, Poland
| | - Daniel T. Gryko
- Institute
of Organic Chemistry of Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| |
Collapse
|
7
|
Du J, Tao X, Begušić T, Wei L. Computational Design of Molecular Probes for Electronic Preresonance Raman Scattering Microscopy. J Phys Chem B 2023; 127:4979-4988. [PMID: 37226966 PMCID: PMC10676804 DOI: 10.1021/acs.jpcb.3c00699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Recently developed electronic preresonance stimulated Raman scattering (epr-SRS) microscopy, in which the Raman signal of a dye is significantly boosted by setting the incident laser frequency near the electronic excitation energy, has pushed the sensitivity of SRS microscopy close to that offered by confocal fluorescence microscopy. Prominently, the maintained narrow line-width of epr-SRS also offers high multiplexity that breaks the "color barrier" in optical microscopy. However, detailed understanding of the fundamental mechanism in these epr-SRS dyes still remains elusive. Here, we combine experiments with theoretical modeling to investigate the structure-function relationship, aiming to facilitate the design of new probes and expanding epr-SRS palettes. Our ab initio approach employing the displaced harmonic oscillator (DHO) model provides a consistent agreement between simulated and experimental SRS intensities of various triple-bond bearing epr-SRS probes with distinct scaffolds. We further review two popular approximate expressions for epr-SRS, namely the short-time and Albrecht A-term equations, and compare them to the DHO model. Overall, the theory allows us to illustrate how the observed intensity differences between molecular scaffolds stem from the coupling strength between the electronic excitation and the targeted vibrational mode, leading to a general design strategy for highly sensitive next-generation vibrational imaging probes.
Collapse
Affiliation(s)
- Jiajun Du
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Lu Wei
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
8
|
Wenzel M, Mitric R. Internal conversion rates from the extended thawed Gaussian approximation: Theory and validation. J Chem Phys 2023; 158:034105. [PMID: 36681643 DOI: 10.1063/5.0130340] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The theoretical prediction of the rates of nonradiative processes in molecules is fundamental in assessing their emissive properties. In this context, global harmonic models have been widely used to simulate vibronic spectra as well as internal conversion rates and to predict photoluminescence quantum yields. However, these simplified models suffer from the limitations that are inherent to the harmonic approximation and can have a severe effect on the calculated internal conversion rates. Therefore, the development of more accurate semiclassical methods is highly desirable. Here, we introduce a procedure for the calculation of nonradiative rates in the framework of the time-dependent semi-classical Extended Thawed Gaussian Approximation (ETGA). We systematically investigate the performance of the ETGA method by comparing it to the adiabatic and vertical harmonic methods, which belong to the class of widely used global harmonic models. Its performance is tested in potentials that cannot be treated adequately by global harmonic models, beginning with Morse potentials of varying anharmonicity followed by a double well potential. The calculated radiative and nonradiative internal conversion rates are compared to reference values based on exact quantum dynamics. We find that the ETGA has the capability to predict internal conversion rates in anharmonic systems with an appreciable energy gap, whereas the global harmonic models prove to be insufficient.
Collapse
Affiliation(s)
- Michael Wenzel
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| | - Roland Mitric
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
9
|
Solowan HP, Malý P, Brixner T. Direct comparison of molecular-beam versus liquid-phase pump-probe and two-dimensional spectroscopy on the example of azulene. J Chem Phys 2022; 157:044201. [DOI: 10.1063/5.0088365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Although azulene's anomalous fluorescence originating from S2 rather than from S1 is the textbook example for the violation of Kasha's rule, the understanding of the underlying processes is still a subject of investigation. Here, we use action-based coherent two-dimensional electronic spectroscopy (2DES) to measure a single Liouville-space response pathway from S0 via S1 to the S2 state of azulene. We directly compare this sequential excitation in liquid phase detecting S2 fluorescence and in a molecular beam detecting photoionized cations, using the S2 anomalous emission to our advantage. We complement the 2DES study by pump-probe measurements of S1 excitation dynamics, including vibrational relaxation and passage through a conical intersection. The direct comparison of liquid and gas phase allows us to assess the effect of the solvent and the interplay of intra- and inter-molecular energy relaxation.
Collapse
Affiliation(s)
| | - Pavel Malý
- Institute of Physical and Theoretical Chemistry, Julius-Maximilians-Universität Würzburg, Germany
| | - Tobias Brixner
- Institut fuer Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Germany
| |
Collapse
|
10
|
Begušić T, Tapavicza E, Vaníček J. Applicability of the Thawed Gaussian Wavepacket Dynamics to the Calculation of Vibronic Spectra of Molecules with Double-Well Potential Energy Surfaces. J Chem Theory Comput 2022; 18:3065-3074. [PMID: 35420803 DOI: 10.1021/acs.jctc.2c00030] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Simulating vibrationally resolved electronic spectra of anharmonic systems, especially those involving double-well potential energy surfaces, often requires expensive quantum dynamics methods. Here, we explore the applicability and limitations of the recently proposed single-Hessian thawed Gaussian approximation for the simulation of spectra of systems with double-well potentials, including 1,2,4,5-tetrafluorobenzene, ammonia, phosphine, and arsine. This semiclassical wavepacket approach is shown to be more robust and to provide more accurate spectra than the conventional harmonic approximation. Specifically, we identify two cases in which the Gaussian wavepacket method is especially useful due to the breakdown of the harmonic approximation: (i) when the nuclear wavepacket is initially at the top of the potential barrier but delocalized over both wells, e.g., along a low-frequency mode, and (ii) when the wavepacket has enough energy to classically go over the low potential energy barrier connecting the two wells. The method is efficient and requires only a single classical ab initio molecular dynamics trajectory, in addition to the data required to compute the harmonic spectra. We also present an improved algorithm for computing the wavepacket autocorrelation function, which guarantees that the evaluated correlation function is continuous for arbitrary size of the time step.
Collapse
Affiliation(s)
- Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
11
|
Manian A, Shaw RA, Lyskov I, Russo SP. Exciton Dynamics of a Diketo-Pyrrolopyrrole Core for All Low-Lying Electronic Excited States Using Density Functional Theory-Based Methods. J Chem Theory Comput 2022; 18:1838-1848. [PMID: 35196857 DOI: 10.1021/acs.jctc.2c00070] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ab initio treatments of interexcited state internal conversion (IC) are more often than not missing from exciton dynamic descriptions, because of their inherent complexity. Here, we define "interexcited state IC" as a same-spin nonradiative transition between states i and j, where i ≠ j ≠ 0. Competing directly with multiexciton processes such as singlet fission or triplet photoupconversion, inclusion of this mechanism in the narrative of molecular photophysics would allow for strategic synthesis of chromophores for more efficient photon-harvesting applications. Herein, we present a robust formalism which can model these rates using density functional theory (DFT)-based methods within the Franck-Condon and Herzberg-Teller regime. Using an unsubstituted diketo-pyrrolopyrrole (DPP) core as a case study, we illustrate the exciton dynamics along the first four excited states for both singlet and triplet manifolds, showing ultrafast same-spin transfer mechanisms due to all excited states, excluding the first triplet level, being in close energetic proximity (within 0.8 eV of each other). The resulting electron same-spin rates outcompete the electron spin-flipping intersystem crossing (ISC) rates, with excitons firmly obeying Kasha's rule as they cascade down from the high-lying excited states toward the lower states. Furthermore, we calculated that only the first singlet excited state displayed a reasonable probability of triplet exciton generation, of ∼40%, with a near-zero chance of the exciton reverting to the singlet manifold once the electron-hole pair are of parallel spin.
Collapse
Affiliation(s)
- Anjay Manian
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, 3000, Australia
| | - Robert A Shaw
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| | - Igor Lyskov
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, 3000, Australia
| | - Salvy P Russo
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, 3000, Australia
| |
Collapse
|
12
|
Awuku S, Bradley SJ, Ghiggino KP, Steer RP, Stevens AL, White JM, Yeow C. Photophysics and spectroscopy of 1,2-Benzazulene. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.139114] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
13
|
Roulet J, Vaníček J. An implicit split-operator algorithm for the nonlinear time-dependent Schrödinger equation. J Chem Phys 2021; 155:204109. [PMID: 34852494 DOI: 10.1063/5.0071153] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The explicit split-operator algorithm is often used for solving the linear and nonlinear time-dependent Schrödinger equations. However, when applied to certain nonlinear time-dependent Schrödinger equations, this algorithm loses time reversibility and second-order accuracy, which makes it very inefficient. Here, we propose to overcome the limitations of the explicit split-operator algorithm by abandoning its explicit nature. We describe a family of high-order implicit split-operator algorithms that are norm-conserving, time-reversible, and very efficient. The geometric properties of the integrators are proven analytically and demonstrated numerically on the local control of a two-dimensional model of retinal. Although they are only applicable to separable Hamiltonians, the implicit split-operator algorithms are, in this setting, more efficient than the recently proposed integrators based on the implicit midpoint method.
Collapse
Affiliation(s)
- Julien Roulet
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénieries Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénieries Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
14
|
Prompt and Long-Lived Anti-Kasha Emission from Organic Dyes. Molecules 2021; 26:molecules26226999. [PMID: 34834093 PMCID: PMC8623836 DOI: 10.3390/molecules26226999] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 11/12/2021] [Accepted: 11/15/2021] [Indexed: 11/17/2022] Open
Abstract
Anti-Kasha behavior has been the subject of intense debate in the last few years, as demonstrated by the high number of papers appearing in the literature on this topic, dealing with both mechanistic and applicative aspects of this phenomenon. Examples of anomalous emitters reported in the last 10 years are collected in the present review, which is focused on strictly anti-Kasha organic molecules displaying radiative deactivation from Sn and/or Tn, with n greater than 1.
Collapse
|
15
|
Behera SK, Park SY, Gierschner J. Duale Emission: Klassen, Mechanismen und Bedingungen. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202009789] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Santosh Kumar Behera
- Madrid Institute for Advanced Studies IMDEA Nanociencia Ciudad Universitaria de Cantoblanco C/ Faraday 9 28049 Madrid Spanien
| | - Soo Young Park
- Laboratory of Supramolecular Optoelectronic Materials and Research Institute of Advanced Materials (RIAM) Department of Materials Science and Engineering Seoul National University ENG 445 Seoul 08826 Korea
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies IMDEA Nanociencia Ciudad Universitaria de Cantoblanco C/ Faraday 9 28049 Madrid Spanien
| |
Collapse
|
16
|
Begušić T, Vaníček J. Finite-Temperature, Anharmonicity, and Duschinsky Effects on the Two-Dimensional Electronic Spectra from Ab Initio Thermo-Field Gaussian Wavepacket Dynamics. J Phys Chem Lett 2021; 12:2997-3005. [PMID: 33733773 PMCID: PMC8006135 DOI: 10.1021/acs.jpclett.1c00123] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Accepted: 03/08/2021] [Indexed: 05/28/2023]
Abstract
Accurate description of finite-temperature vibrational dynamics is indispensable in the computation of two-dimensional electronic spectra. Such simulations are often based on the density matrix evolution, statistical averaging of initial vibrational states, or approximate classical or semiclassical limits. While many practical approaches exist, they are often of limited accuracy and difficult to interpret. Here, we use the concept of thermo-field dynamics to derive an exact finite-temperature expression that lends itself to an intuitive wavepacket-based interpretation. Furthermore, an efficient method for computing finite-temperature two-dimensional spectra is obtained by combining the exact thermo-field dynamics approach with the thawed Gaussian approximation for the wavepacket dynamics, which is exact for any displaced, distorted, and Duschinsky-rotated harmonic potential but also accounts partially for anharmonicity effects in general potentials. Using this new method, we directly relate a symmetry breaking of the two-dimensional signal to the deviation from the conventional Brownian oscillator picture.
Collapse
|
17
|
Begušić T, Vaníček J. On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy. J Chem Phys 2020; 153:184110. [DOI: 10.1063/5.0031216] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
18
|
Golubev NV, Begušić T, Vaníček J. On-the-Fly Ab Initio Semiclassical Evaluation of Electronic Coherences in Polyatomic Molecules Reveals a Simple Mechanism of Decoherence. PHYSICAL REVIEW LETTERS 2020; 125:083001. [PMID: 32909765 DOI: 10.1103/physrevlett.125.083001] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Accepted: 06/30/2020] [Indexed: 06/11/2023]
Abstract
Irradiation of a molecular system by an intense laser field can trigger dynamics of both electronic and nuclear subsystems. The lighter electrons usually move on much faster, attosecond timescale but the slow nuclear rearrangement damps ultrafast electronic oscillations, leading to the decoherence of the electronic dynamics within a few femtoseconds. We show that a simple, single-trajectory semiclassical scheme can evaluate the electronic coherence time in polyatomic molecules accurately by demonstrating an excellent agreement with full-dimensional quantum calculations. In contrast to numerical quantum methods, the semiclassical one reveals the physical mechanism of decoherence beyond the general blame on nuclear motion. In the propiolic acid, the rate of decoherence and the large deviation from the static frequency of electronic oscillations are quantitatively described with just two semiclassical parameters-the phase space distance and signed area between the trajectories moving on two electronic surfaces. Because it evaluates the electronic structure on the fly, the semiclassical technique avoids the "curse of dimensionality" and should be useful for preselecting molecules for experimental studies.
Collapse
Affiliation(s)
- Nikolay V Golubev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
19
|
Behera SK, Park SY, Gierschner J. Dual Emission: Classes, Mechanisms, and Conditions. Angew Chem Int Ed Engl 2020; 60:22624-22638. [PMID: 32783293 DOI: 10.1002/anie.202009789] [Citation(s) in RCA: 81] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2020] [Revised: 08/10/2020] [Indexed: 12/28/2022]
Abstract
There has been much interest in dual-emission materials in the past few years for materials and life science applications; however, a systematic overview of the underlying processes is so-far missing. We resolve this issue herein by classifying dual-emission (DE) phenomena as relying on one emitter with two emitting states (DE1), two independent emitters (DE2), or two correlated emitters (DE3). Relevant DE mechanisms for materials science are then briefly described together with the electronic and/or geometrical conditions under which they occur. For further reading, references are given that offer detailed insight into the complex mechanistic aspects of the various DE processes or provide overviews on materials families or their applications. By avoiding ambiguities and misinterpretations, this systematic, insightful Review might inspire future targeted designs of DE materials.
Collapse
Affiliation(s)
- Santosh Kumar Behera
- Madrid Institute for Advanced Studies, IMDEA Nanociencia, Ciudad Universitaria de Cantoblanco, C/ Faraday 9, 28049, Madrid, Spain
| | - Soo Young Park
- Laboratory of Supramolecular Optoelectronic Materials and Research Institute of Advanced Materials (RIAM), Department of Materials Science and Engineering, Seoul National University, ENG 445, Seoul, 08826, Korea
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies, IMDEA Nanociencia, Ciudad Universitaria de Cantoblanco, C/ Faraday 9, 28049, Madrid, Spain
| |
Collapse
|
20
|
Begušić T, Vaníček J. On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature. J Chem Phys 2020; 153:024105. [PMID: 32668922 DOI: 10.1063/5.0013677] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
To compute and analyze vibrationally resolved electronic spectra at zero temperature, we have recently implemented the on-the-fly ab initio extended thawed Gaussian approximation [A. Patoz et al., J. Phys. Chem. Lett. 9, 2367 (2018)], which accounts for anharmonicity, mode-mode coupling, and Herzberg-Teller effects. Here, we generalize this method in order to evaluate spectra at non-zero temperature. In line with thermo-field dynamics, we transform the von Neumann evolution of the coherence component of the density matrix to the Schrödinger evolution of a wavefunction in an augmented space with twice as many degrees of freedom. Due to the efficiency of the extended thawed Gaussian approximation, this increase in the number of coordinates results in nearly no additional computational cost. More specifically, compared to the original, zero-temperature approach, the finite-temperature method requires no additional ab initio electronic structure calculations. At the same time, the new approach allows for a clear distinction among finite-temperature, anharmonicity, and Herzberg-Teller effects on spectra. We show, on a model Morse system, the advantages of the finite-temperature thawed Gaussian approximation over the commonly used global harmonic methods and apply it to evaluate the symmetry-forbidden absorption spectrum of benzene, where all of the aforementioned effects contribute.
Collapse
Affiliation(s)
- Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
21
|
Veys K, Escudero D. Computational Protocol To Predict Anti-Kasha Emissions: The Case of Azulene Derivatives. J Phys Chem A 2020; 124:7228-7237. [DOI: 10.1021/acs.jpca.0c05205] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Koen Veys
- Quantum Chemistry and Physical Chemistry Section, Department of Chemistry, KU Leuven, Celestijnenlaan 200f, 3001 Leuven, Belgium
| | - Daniel Escudero
- Quantum Chemistry and Physical Chemistry Section, Department of Chemistry, KU Leuven, Celestijnenlaan 200f, 3001 Leuven, Belgium
| |
Collapse
|