1
|
Le Dé B, Jaouadi A, Mangaud E, Chin AW, Desouter-Lecomte M. Managing temperature in open quantum systems strongly coupled with structured environments. J Chem Phys 2024; 160:244102. [PMID: 38913841 DOI: 10.1063/5.0214051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 06/06/2024] [Indexed: 06/26/2024] Open
Abstract
In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.
Collapse
Affiliation(s)
- Brieuc Le Dé
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | - Amine Jaouadi
- LyRIDS, ECE Paris, Graduate School of Engineering, Paris F-75015, France
| | - Etienne Mangaud
- MSME, Université Gustave Eiffel, UPEC, CNRS, F-77454 Marne-La-Vallée, France
| | - Alex W Chin
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | | |
Collapse
|
2
|
Ke Y, Richardson JO. Insights into the mechanisms of optical cavity-modified ground-state chemical reactions. J Chem Phys 2024; 160:224704. [PMID: 38856061 DOI: 10.1063/5.0200410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Accepted: 05/27/2024] [Indexed: 06/11/2024] Open
Abstract
In this work, we systematically investigate the mechanisms underlying the rate modification of ground-state chemical reactions in an optical cavity under vibrational strong-coupling conditions. We employ a symmetric double-well description of the molecular potential energy surface and a numerically exact open quantum system approach-the hierarchical equations of motion in twin space with a matrix product state solver. Our results predict the existence of multiple peaks in the photon frequency-dependent rate profile for a strongly anharmonic molecular system with multiple vibrational transition energies. The emergence of a new peak in the rate profile is attributed to the opening of an intramolecular reaction pathway, energetically fueled by the cavity photon bath through a resonant cavity mode. The peak intensity is determined jointly by kinetic factors. Going beyond the single-molecule limit, we examine the effects of the collective coupling of two molecules to the cavity. We find that when two identical molecules are simultaneously coupled to the same resonant cavity mode, the reaction rate is further increased. This additional increase is associated with the activation of a cavity-induced intermolecular reaction channel. Furthermore, the rate modification due to these cavity-promoted reaction pathways remains unaffected, regardless of whether the molecular dipole moments are aligned in the same or opposite direction as the light polarization.
Collapse
Affiliation(s)
- Yaling Ke
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Jeremy O Richardson
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
3
|
Yang H, Li W, Ren J, Shuai Z. Time-Dependent Density Matrix Renormalization Group Method for Quantum Transport with Phonon Coupling in Molecular Junction. J Chem Theory Comput 2023; 19:6070-6081. [PMID: 37669099 DOI: 10.1021/acs.jctc.3c00340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.
Collapse
Affiliation(s)
- Hengrui Yang
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, People's Republic of China
| |
Collapse
|
4
|
Dan X, Shi Q. Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method. J Chem Phys 2023; 159:044101. [PMID: 37486050 DOI: 10.1063/5.0155172] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Accepted: 06/30/2023] [Indexed: 07/25/2023] Open
Abstract
Hydrogen atom scattering on metal surfaces is investigated based on a simplified Newns-Anderson model. Both the nuclear and electronic degrees of freedom are treated quantum mechanically. By partitioning all the surface electronic states as the bath, the hierarchical equations of motion method for the fermionic bath is employed to simulate the scattering dynamics. It is found that, with a reasonable set of parameters, the main features of the recent experimental studies of hydrogen atom scattering on metal surfaces can be reproduced. Vibrational states on the chemisorption state whose energies are close to the incident energy are found to play an important role, and the scattering process is dominated by a single-pass electronic transition forth and back between the diabatic physisorption and chemisorption states. Further study on the effects of the atom-surface coupling strength reveals that, upon increasing the atom-surface coupling strength, the scattering mechanism changes from typical nonadiabatic transitions to dynamics in the electronic friction regime.
Collapse
Affiliation(s)
- Xiaohan Dan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China and University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
5
|
Cainelli M, Borrelli R, Tanimura Y. Effect of mixed Frenkel and charge transfer states in time-gated fluorescence spectra of perylene bisimides H-aggregates: Hierarchical equations of motion approach. J Chem Phys 2022; 157:084103. [DOI: 10.1063/5.0102000] [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
We theoretically investigated the effect of mixed Frenkel (F) and charge transfer (CT) states on the spectral properties of perylene bisimide (PBI) derivatives, focusing on the role of strong electron-phonon interactions. The model consists of a four-level system described by the Holstein Hamiltonian coupled to independent local heat-baths on each site, described by Brownian spectral distribution functions. We employ the reduced hierarchical equations of motion (HEOM) approach to calculate the time evolution of the system and compare it to the pure F exciton cases. We compute the absorption and time-gated fluorescence (TGF) spectra for different exciton transfer integrals and F-CT band gap conditions. The coherence length of excitons ($N_{coh}$) is evaluated employing two different definitions. We observe the presence of an excited hot state peak whose intensity is associated with the delocalization of the excited species and ultrafast dynamics that are solely dependent on the frequency of the local bath. The results indicate that the inclusion of CT states promotes localization of the excitons which is manifested in a decrease in the intensity of the hot state peak and the 0--1 peak, and an increase in the intensity of the 0--0 emission peak in TGF spectrum, leading to a decrease of $N_{coh}$.
Collapse
Affiliation(s)
| | - Raffaele Borrelli
- Department of Agricoltural Science, Università degli Studi di Torino, Italy
| | | |
Collapse
|
6
|
Fay TP. A simple improved low temperature correction for the hierarchical equations of motion. J Chem Phys 2022; 157:054108. [DOI: 10.1063/5.0100365] [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/14/2022] Open
Abstract
The study of open system quantum dynamics has been transformed by the hierarchical equations of motion (HEOM) method, which gives the exact dynamics for a system coupled to a harmonic bath at arbitrary temperature and system-bath coupling strength. However in its standard form the method is only consistent with the weak-coupling quantum master equation at all temperatures when many auxiliary density operators are included in the hierarchy, even when low temperature corrections are included. Here we propose a new low temperature correction scheme for the termination of the hierarchy based on Zwanzig projection which alleviates this problem, and restores consistency with the weak-coupling master equation with a minimal hierarchy. The utility of the new correction scheme is demonstrated on a range of model systems, including the Fenna-Metthews-Olson complex. The new closure is found to improve convergence of the HEOM even beyond the weak-coupling limit and is very straightforward to implement in existing HEOM codes.
Collapse
Affiliation(s)
- Thomas Patrick Fay
- Department of Chemistry, University of California Berkeley Department of Chemistry, United States of America
| |
Collapse
|
7
|
Ke Y, Kaspar C, Erpenbeck A, Peskin U, Thoss M. Nonequilibrium reaction rate theory: Formulation and implementation within the hierarchical equations of motion approach. J Chem Phys 2022; 157:034103. [DOI: 10.1063/5.0098545] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The study of chemical reactions in environments under nonequilibrium conditions has been of interest recently in a variety of contexts, including current-induced reactions in molecular junctions and scanning tunneling microscopy experiments. In this work, we outline a fully quantum mechanical, numerically exact approach to describe chemical reaction rates in such nonequilibrium situations. The approach is based on an extension of the flux correlation function formalism to nonequilibrium conditions and uses a mixed real and imaginary time hierarchical equations of motion approach for the calculation of rate constants. As a specific example, we investigate current-induced intramolecular proton transfer reactions in a molecular junction for different applied bias voltages and molecule-lead coupling strengths.
Collapse
Affiliation(s)
- Yaling Ke
- Institute of Physics, Albert-Ludwigs-Universität Freiburg, Germany
| | | | | | - Uri Peskin
- Chemistry, Technion Israel Institute of Technology, Israel
| | - Michael Thoss
- University of Freiburg Institute of Physics, Germany
| |
Collapse
|