1
|
Takahashi H, Borrelli R. Tensor-Train Format Hierarchical Equations of Motion Formalism: Charge Transfer in Organic Semiconductors via Dissipative Holstein Models. J Chem Theory Comput 2024. [PMID: 39152908 DOI: 10.1021/acs.jctc.4c00711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2024]
Abstract
Hierarchical Equations of Motion (HEOM) in the Tensor-Train (TT) representation is applied to study the charge-transfer dynamics in organic semiconductors (OSCs). The theoretical formulation as well as the basic computational aspects of HEOM-TT are discussed in detail. Charge transfer in OSCs is modeled using dissipative polaronic models that incorporate the effects of both high- and low-frequency molecular vibrations, and it is simulated in a fully quantum and nonperturbative manner, which has not been studied intensively. The capability of treating complex electron-vibrational systems is examined by analyzing and comparing the numerical behavior of the time-dependent variational approach and the time-Alternating Minimal Energy methods and by calculating the current autocorrelation function and diffusivity across various models. Our results indicate that the HEOM-TT framework offers a robust tool for the detailed analysis of complex polaronic systems, suggesting its potential for broader applications.
Collapse
|
2
|
Yao Q, Ji Q, Li X, Zhang Y, Chen X, Ju MG, Liu J, Wang J. Machine Learning Accelerates Precise Excited-State Potential Energy Surface Calculations on a Quantum Computer. J Phys Chem Lett 2024; 15:7061-7068. [PMID: 38950102 DOI: 10.1021/acs.jpclett.4c01445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2024]
Abstract
Electronically excited-state problems represent a crucial research field in quantum chemistry, closely related to numerous practical applications in photophysics and photochemistry. The emerging of quantum computing provides a promising computational paradigm to solve the Schrödinger equation for predicting potential energy surfaces (PESs). Here, we present a deep neural network model to predict parameters of the quantum circuits within the framework of variational quantum deflation and subspace search variational quantum eigensolver, which are two popular excited-state algorithms to implement on a quantum computer. The new machine learning-assisted algorithm is employed to study the excited-state PESs of small molecules, achieving highly accurate predictions. We then apply this algorithm to study the excited-state properties of the ArF system, which is essential to a gas laser. Through this study, we believe that with future advancements in hardware capabilities, quantum computing could be harnessed to solve excited-state problems for a broad range of systems.
Collapse
Affiliation(s)
- Qianjun Yao
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Qun Ji
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Xiaopeng Li
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei 230026, China
| | - Yehui Zhang
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Xinyu Chen
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Ming-Gang Ju
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
| | - Jie Liu
- Hefei National Laboratory, Hefei 230088, China
| | - Jinlan Wang
- Key Laboratory of Quantum Materials and Devices of Ministry of Education, School of Physics, Southeast University, Nanjing 211189, China
- Suzhou Laboratory, Suzhou 215009, China
| |
Collapse
|
3
|
Delgado F, Enríquez M. Quantum Entanglement and State-Transference in Fenna-Matthews-Olson Complexes: A Post-Experimental Simulation Analysis in the Computational Biology Domain. Int J Mol Sci 2023; 24:10862. [PMID: 37446061 DOI: 10.3390/ijms241310862] [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: 05/18/2023] [Revised: 06/13/2023] [Accepted: 06/27/2023] [Indexed: 07/15/2023] Open
Abstract
Fenna-Mathews-Olson complexes participate in the photosynthetic process of Sulfur Green Bacteria. These biological subsystems exhibit quantum features which possibly are responsible for their high efficiency; the latter may comprise multipartite entanglement and the apparent tunnelling of the initial quantum state. At first, to study these aspects, a multidisciplinary approach including experimental biology, spectroscopy, physics, and math modelling is required. Then, a global computer modelling analysis is achieved in the computational biology domain. The current work implements the Hierarchical Equations of Motion to numerically solve the open quantum system problem regarding this complex. The time-evolved states obtained with this method are then analysed under several measures of entanglement, some of them already proposed in the literature. However, for the first time, the maximum overlap with respect to the closest separable state is employed. This authentic multipartite entanglement measure provides information on the correlations, not only based on the system bipartitions as in the usual analysis. Our study has led us to note a different view of FMO multipartite entanglement as tiny contributions to the global entanglement suggested by other more basic measurements. Additionally, in another related trend, the initial state, considered as a Förster Resonance Energy Transfer, is tracked using a novel approach, considering how it could be followed under the fidelity measure on all possible permutations of the FMO subsystems through its dynamical evolution by observing the tunnelling in the most probable locations. Both analyses demanded significant computational work, making for a clear example of the complexity required in computational biology.
Collapse
Affiliation(s)
- Francisco Delgado
- School of Engineering and Sciences, Tecnologico de Monterrey, Atizapan 52926, Mexico
| | - Marco Enríquez
- School of Engineering and Sciences, Tecnologico de Monterrey, Santa Fe 01389, Mexico
| |
Collapse
|
4
|
Atsango AO, Montoya-Castillo A, Markland TE. An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra. J Chem Phys 2023; 158:074107. [PMID: 36813724 DOI: 10.1063/5.0138671] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump-probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.
Collapse
Affiliation(s)
- Austin O Atsango
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
5
|
Chen L, Bennett DIG, Eisfeld A. Calculating non-linear response functions for multi-dimensional electronic spectroscopy using dyadic non-Markovian quantum state diffusion. J Chem Phys 2022; 157:114104. [DOI: 10.1063/5.0107925] [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
We present a methodology for simulating multi-dimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space, but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.
Collapse
Affiliation(s)
- Lipeng Chen
- Department of Chemistry, Max-Planck-Institute for the Physics of Complex Systems, Germany
| | - Doran I. G Bennett
- Chemistry, Southern Methodist University Department of Chemistry, United States of America
| | | |
Collapse
|
6
|
Fang J, Chen ZH, Su Y, Zhu ZF, Wang Y, Xu RX, Yan Y. Coherent excitation energy transfer in model photosynthetic reaction center: Effects of non-Markovian quantum environment. J Chem Phys 2022; 157:084119. [DOI: 10.1063/5.0104641] [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
Excitation energy transfer (EET) and electron transfer (ET) are crucially involved in photosynthetic processes. In reality, the photosynthetic reaction center constitutes an open quantum system of EET and ET, which manifests an interplay of pigments, solar light and phonon baths. So far theoretical studies have been mainly based on master equation approaches in the Markovian condition. The non-Markovian environmental effect, which may play a crucial role, has not been sufficiently considered. In this work, we propose a mixed dynamic approach to investigate this open system. The influence of phonon bath is treated via the exact dissipaton equation of motion (DEOM) while that of photon bath is via the Lindblad master equation. Specifically, we explore the effect of non-Markovian quantum phonon bath on the coherent transfer dynamics and its manipulation on the current--voltage behavior. Distinguished from the results of completely Markovian Lindblad equation and those adopting classical environment description, the mixed DEOM--Lindblad simulations exhibittransfer coherence up to a few hundreds femtosecondsand the related environmental manipulation effect on current.These non-Markovian quantum coherent effects may be extended tomore complex and realistic systems and be helpful to thedesign of organic photovoltaic devices.
Collapse
Affiliation(s)
- Jie Fang
- University of Science and Technology of China, China
| | - Zi-Hao Chen
- University of Science and Technology of China, China
| | - Yu Su
- Department of Chemical Physics, University of Science and Technology of China, China
| | - Zi-Fan Zhu
- University of Science and Technology of China, China
| | - Yao Wang
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, China
| | - Rui-Xue Xu
- University of Science and Technology of China, China
| | - YiJing Yan
- Department of Chemical Physics, USTC, China
| |
Collapse
|
7
|
Ko L, Cook RL, Whaley KB. Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states. J Chem Phys 2022; 156:244108. [DOI: 10.1063/5.0082822] [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
We develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.
Collapse
Affiliation(s)
- Liwen Ko
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy Nanoscience Institute at Berkeley, Berkeley, California 94720, USA
| | - Robert L. Cook
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy Nanoscience Institute at Berkeley, Berkeley, California 94720, USA
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy Nanoscience Institute at Berkeley, Berkeley, California 94720, USA
| |
Collapse
|
8
|
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
|
9
|
Cho KH, Rhee YM. Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes. Phys Chem Chem Phys 2021; 23:26623-26639. [PMID: 34842245 DOI: 10.1039/d1cp04615b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
10
|
Abstract
Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.
Collapse
Affiliation(s)
- Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
11
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 166] [Impact Index Per Article: 55.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
12
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
13
|
Kaspar C, Thoss M. Efficient Steady-State Solver for the Hierarchical Equations of Motion Approach: Formulation and Application to Charge Transport through Nanosystems. J Phys Chem A 2021; 125:5190-5200. [PMID: 34102055 DOI: 10.1021/acs.jpca.1c02863] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.
Collapse
Affiliation(s)
- C Kaspar
- Institute of Physics, University of Freiburg, Hermann-Herder-Straße 3, D-79104 Freiburg, Germany
| | - M Thoss
- Institute of Physics, University of Freiburg, Hermann-Herder-Straße 3, D-79104 Freiburg, Germany.,EUCOR Centre for Quantum Science and Quantum Computing, University of Freiburg, Hermann-Herder-Straße 3, D-79104 Freiburg, Germany
| |
Collapse
|
14
|
Yan Y, Xu M, Li T, Shi Q. Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors. J Chem Phys 2021; 154:194104. [PMID: 34240893 DOI: 10.1063/5.0050720] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 104 effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.
Collapse
Affiliation(s)
- Yaming Yan
- 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Meng Xu
- 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Tianchu Li
- 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, 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; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| |
Collapse
|
15
|
Borrelli R, Dolgov S. Expanding the Range of Hierarchical Equations of Motion by Tensor-Train Implementation. J Phys Chem B 2021; 125:5397-5407. [DOI: 10.1021/acs.jpcb.1c02724] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Raffaele Borrelli
- DISAFA, University of Torino, Largo Paolo Braccini 2, Grugliasco 10095, Italy
| | - Sergey Dolgov
- University of Bath, Claverton Down, BA2 7AY Bath, United Kingdom
| |
Collapse
|
16
|
Borrelli R, Gelin MF. Finite temperature quantum dynamics of complex systems: Integrating
thermo‐field
theories and
tensor‐train
methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1539] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
| | - Maxim F. Gelin
- School of Sciences Hangzhou Dianzi University Hangzhou China
| |
Collapse
|
17
|
Maity S, Daskalakis V, Elstner M, Kleinekathöfer U. Multiscale QM/MM molecular dynamics simulations of the trimeric major light-harvesting complex II. Phys Chem Chem Phys 2021; 23:7407-7417. [PMID: 33876100 DOI: 10.1039/d1cp01011e] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Photosynthetic processes are driven by sunlight. Too little of it and the photosynthetic machinery cannot produce the reductive power to drive the anabolic pathways. Too much sunlight and the machinery can get damaged. In higher plants, the major Light-Harvesting Complex (LHCII) efficiently absorbs the light energy, but can also dissipate it when in excess (quenching). In order to study the dynamics related to the quenching process but also the exciton dynamics in general, one needs to accurately determine the so-called spectral density which describes the coupling between the relevant pigment modes and the environmental degrees of freedom. To this end, Born-Oppenheimer molecular dynamics simulations in a quantum mechanics/molecular mechanics (QM/MM) fashion utilizing the density functional based tight binding (DFTB) method have been performed for the ground state dynamics. Subsequently, the time-dependent extension of the long-range-corrected DFTB scheme has been employed for the excited state calculations of the individual chlorophyll-a molecules in the LHCII complex. The analysis of this data resulted in spectral densities showing an astonishing agreement with the experimental counterpart in this rather large system. This consistency with an experimental observable also supports the accuracy, robustness, and reliability of the present multi-scale scheme. To the best of our knowledge, this is the first theoretical attempt on this large complex system is ever made to accurately simulate the spectral density. In addition, the resulting spectral densities and site energies were used to determine the exciton transfer rate within a special pigment pair consisting of a chlorophyll-a and a carotenoid molecule which is assumed to play a role in the balance between the light harvesting and quenching modes.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany.
| | | | | | | |
Collapse
|
18
|
Tanimura Y. Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM). J Chem Phys 2021; 153:020901. [PMID: 32668942 DOI: 10.1063/5.0011599] [Citation(s) in RCA: 140] [Impact Index Per Article: 46.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.
Collapse
Affiliation(s)
- Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
19
|
Janković V, Mančal T. Nonequilibrium steady-state picture of incoherent light-induced excitation harvesting. J Chem Phys 2020; 153:244110. [PMID: 33380098 DOI: 10.1063/5.0029918] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We formulate a comprehensive theoretical description of excitation harvesting in molecular aggregates photoexcited by weak incoherent radiation. An efficient numerical scheme that respects the continuity equation for excitation fluxes is developed to compute the nonequilibrium steady state (NESS) arising from the interplay between excitation generation, excitation relaxation, dephasing, trapping at the load, and recombination. The NESS is most conveniently described in the so-called preferred basis in which the steady-state excitonic density matrix is diagonal. The NESS properties are examined by relating the preferred-basis description to the descriptions in the site or excitonic bases. Focusing on a model photosynthetic dimer, we find that the NESS in the limit of long trapping time is quite similar to the excited-state equilibrium in which the stationary coherences originate from the excitation-environment entanglement. For shorter trapping times, we demonstrate how the properties of the NESS can be extracted from the time-dependent description of an incoherently driven but unloaded dimer. This relation between stationary and time-dependent pictures is valid, provided that the trapping time is longer than the decay time of dynamic coherences accessible in femtosecond spectroscopy experiments.
Collapse
Affiliation(s)
- Veljko Janković
- Scientific Computing Laboratory, Center for the Study of Complex Systems, Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Tomáš Mančal
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 121 16 Prague 2, Czech Republic
| |
Collapse
|
20
|
Janković V, Mančal T. Exact description of excitonic dynamics in molecular aggregates weakly driven by light. J Chem Phys 2020; 153:244122. [PMID: 33380075 DOI: 10.1063/5.0029914] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a rigorous theoretical description of excitonic dynamics in molecular light-harvesting aggregates photoexcited by weak-intensity radiation of arbitrary properties. While the interaction with light is included up to the second order, the treatment of the excitation-environment coupling is exact and results in an exact expression for the reduced excitonic density matrix that is manifestly related to the spectroscopic picture of the photoexcitation process. This expression takes fully into account the environmental reorganization processes triggered by the two interactions with light. This is particularly important for slow environments and/or strong excitation-environment coupling. Within the exponential decomposition scheme, we demonstrate how our result can be recast as the hierarchy of equations of motion (HEOM) that explicitly and consistently includes the photoexcitation step. We analytically describe the environmental reorganization dynamics triggered by a delta-like excitation of a single chromophore and demonstrate how our HEOM, in appropriate limits, reduces to the Redfield equations comprising a pulsed photoexcitation and the nonequilibrium Förster theory. We also discuss the relation of our formalism to the combined Born-Markov-HEOM approaches in the case of excitation by thermal light.
Collapse
Affiliation(s)
- Veljko Janković
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 121 16 Prague 2, Czech Republic
| | - Tomáš Mančal
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 121 16 Prague 2, Czech Republic
| |
Collapse
|
21
|
Sahu A, Kurian JS, Tiwari V. Vibronic resonance is inadequately described by one-particle basis sets. J Chem Phys 2020; 153:224114. [DOI: 10.1063/5.0029027] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Amitav Sahu
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | - Jo Sony Kurian
- Department of Chemistry, Indian Institute of Science Education and Research, Bhopal, Madhya Pradesh 462066, India
| | - Vivek Tiwari
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| |
Collapse
|
22
|
Wang YC, Zhao Y. The hierarchical stochastic schrödinger equations: Theory and applications. CHINESE J CHEM PHYS 2020. [DOI: 10.1063/1674-0068/cjcp2009165] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
23
|
Yan Y, Xing T, Shi Q. A new method to improve the numerical stability of the hierarchical equations of motion for discrete harmonic oscillator modes. J Chem Phys 2020; 153:204109. [DOI: 10.1063/5.0027962] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Tao Xing
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| |
Collapse
|
24
|
Suzuki Y, Watanabe H, Okiyama Y, Ebina K, Tanaka S. Comparative study on model parameter evaluations for the energy transfer dynamics in Fenna–Matthews–Olson complex. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110903] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
25
|
Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10196821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.
Collapse
|
26
|
Yan Y, Liu Y, Xing T, Shi Q. Theoretical study of excitation energy transfer and nonlinear spectroscopy of photosynthetic light‐harvesting complexes using the nonperturbative reduced dynamics method. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1498] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Yanying Liu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Tao Xing
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| |
Collapse
|
27
|
Ullah A, Han L, Yan YA, Zheng X, Yan Y, Chernyak V. Stochastic equation of motion approach to fermionic dissipative dynamics. II. Numerical implementation. J Chem Phys 2020; 152:204106. [DOI: 10.1063/1.5142166] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Arif Ullah
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lu Han
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yun-An Yan
- School of Physics and Optoelectronic Engineering, Ludong University, Shandong 264025, China
| | - Xiao Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Laboratory for Physical Sciences at the Microscale and iChEM, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Vladimir Chernyak
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemistry, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, USA
| |
Collapse
|
28
|
Han L, Ullah A, Yan YA, Zheng X, Yan Y, Chernyak V. Stochastic equation of motion approach to fermionic dissipative dynamics. I. Formalism. J Chem Phys 2020; 152:204105. [DOI: 10.1063/1.5142164] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Lu Han
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Arif Ullah
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yun-An Yan
- School of Physics and Optoelectronic Engineering, Ludong University, Shandong 264025, China
| | - Xiao Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Laboratory for Physical Sciences at the Microscale & iChEM, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Vladimir Chernyak
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemistry, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, USA
| |
Collapse
|
29
|
Kramer T, Rodríguez M. Effect of disorder and polarization sequences on two-dimensional spectra of light-harvesting complexes. PHOTOSYNTHESIS RESEARCH 2020; 144:147-154. [PMID: 31872335 DOI: 10.1007/s11120-019-00699-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Accepted: 12/10/2019] [Indexed: 06/10/2023]
Abstract
Two-dimensional electronic spectra (2DES) provide unique ways to track the energy transfer dynamics in light-harvesting complexes. The interpretation of the peaks and structures found in experimentally recorded 2DES is often not straightforward, since several processes are imaged simultaneously. The choice of specific pulse polarization sequences helps to disentangle the sometimes convoluted spectra, but brings along other disturbances. We show by detailed theoretical calculations how 2DES of the Fenna-Matthews-Olson complex are affected by rotational and conformational disorder of the chromophores.
Collapse
|
30
|
Li W, Ren J, Shuai Z. Numerical assessment for accuracy and GPU acceleration of TD-DMRG time evolution schemes. J Chem Phys 2020; 152:024127. [PMID: 31941314 DOI: 10.1063/1.5135363] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.
Collapse
Affiliation(s)
- 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 Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, 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
| |
Collapse
|
31
|
Marcus M, Knee GC, Datta A. Towards a spectroscopic protocol for unambiguous detection of quantum coherence in excitonic energy transport. Faraday Discuss 2020; 221:110-132. [DOI: 10.1039/c9fd00068b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We propose a witness for quantum coherence in EET that can be extracted directly from two-pulse pump–probe spectroscopy experimental data.
Collapse
Affiliation(s)
- Max Marcus
- Department of Physics
- University of Warwick
- Coventry
- UK
| | | | - Animesh Datta
- Department of Physics
- University of Warwick
- Coventry
- UK
| |
Collapse
|
32
|
Rajabi S, Mančal T. Energy transfer times in fluorographene-based biomimetic light harvesting antennae. Chem Phys 2019. [DOI: 10.1016/j.chemphys.2019.110477] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
33
|
Provazza J, Coker DF. Multi-level description of the vibronic dynamics of open quantum systems. J Chem Phys 2019; 151:154114. [DOI: 10.1063/1.5120253] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Justin Provazza
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|
34
|
Maguire H, Iles-Smith J, Nazir A. Environmental Nonadditivity and Franck-Condon physics in Nonequilibrium Quantum Systems. PHYSICAL REVIEW LETTERS 2019; 123:093601. [PMID: 31524488 DOI: 10.1103/physrevlett.123.093601] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Indexed: 06/10/2023]
Abstract
We show that for a quantum system coupled to both vibrational and electromagnetic environments, enforcing additivity of their combined influences results in nonequilibrium dynamics that does not respect the Franck-Condon principle. We overcome this shortcoming by employing a collective coordinate representation of the vibrational environment, which permits the derivation of a nonadditive master equation. When applied to a two-level emitter our treatment predicts decreasing photon emission rates with increasing vibrational coupling, consistent with Franck-Condon physics. In contrast, the additive approximation predicts the emission rate to be completely insensitive to vibrations. We find that nonadditivity also plays a key role in the stationary nonequilibrium model behavior, enabling two-level population inversion under incoherent electromagnetic excitation.
Collapse
Affiliation(s)
- Henry Maguire
- Photon Science Institute & School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| | - Jake Iles-Smith
- Photon Science Institute & School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
- Department of Physics and Astronomy, University of Sheffield, Hounsfield Road, Sheffield, S3 7RH, United Kingdom
| | - Ahsan Nazir
- Photon Science Institute & School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
35
|
Rahman H, Kleinekathöfer U. Chebyshev hierarchical equations of motion for systems with arbitrary spectral densities and temperatures. J Chem Phys 2019; 150:244104. [PMID: 31255062 DOI: 10.1063/1.5100102] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The time evolution in open quantum systems, such as a molecular aggregate in contact with a thermal bath, still poses a complex and challenging problem. The influence of the thermal noise can be treated using a plethora of schemes, several of which decompose the corresponding correlation functions in terms of weighted sums of exponential functions. One such scheme is based on the hierarchical equations of motion (HEOM), which is built using only certain forms of bath correlation functions. In the case where the environment is described by a complex spectral density or is at a very low temperature, approaches utilizing the exponential decomposition become very inefficient. Here, we utilize an alternative decomposition scheme for the bath correlation function based on Chebyshev polynomials and Bessel functions to derive a HEOM approach up to an arbitrary order in the environmental coupling. These hierarchical equations are similar in structure to the popular exponential HEOM scheme, but are formulated using the derivatives of the Bessel functions. The proposed scheme is tested up to the fourth order in perturbation theory for a two-level system and compared to benchmark calculations for the case of zero-temperature quantum Ohmic and super-Ohmic noise. Furthermore, the benefits and shortcomings of the present Chebyshev-based hierarchical equations are discussed.
Collapse
Affiliation(s)
- Hasan Rahman
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
36
|
Dunn IS, Tempelaar R, Reichman DR. Removing instabilities in the hierarchical equations of motion: Exact and approximate projection approaches. J Chem Phys 2019; 150:184109. [DOI: 10.1063/1.5092616] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ian S. Dunn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Roel Tempelaar
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
37
|
Kato A, Tanimura Y. Hierarchical Equations of Motion Approach to Quantum Thermodynamics. ACTA ACUST UNITED AC 2019. [DOI: 10.1007/978-3-319-99046-0_24] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023]
|
38
|
|
39
|
|
40
|
Wang Y, Ke Y, Zhao Y. The hierarchical and perturbative forms of stochastic Schrödinger equations and their applications to carrier dynamics in organic materials. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1375] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Yu‐Chen Wang
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yaling Ke
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| | - Yi Zhao
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University Xiamen China
| |
Collapse
|
41
|
Han L, Zhang HD, Zheng X, Yan Y. On the exact truncation tier of fermionic hierarchical equations of motion. J Chem Phys 2018; 148:234108. [PMID: 29935503 DOI: 10.1063/1.5034776] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The hierarchical equations of motion (HEOM) theory is in principle exact for describing the dissipative dynamics of quantum systems linearly coupled to Gaussian environments. In practice, the hierarchy needs to be truncated at a finite tier. We demonstrate that, for general systems described by the fermionic HEOM, the (n+L̃)th-tier truncation with L̃=2NσNν yields the exact density operators up to the nth tier. Here, Nσ = 2 for fermionic systems and Nν is the system degrees of freedom. For noninteracting systems, L̃ is further reduced by half. Such an exact termination pattern originates from the Pauli exclusion principle for fermions, and it holds true regardless of the system-environment coupling strength, the number of coupling reservoirs, or the specific scheme employed to unravel the environment memory contents. The relatively small L̃ emphasizes the nonperturbative nature of the HEOM theory. We also propose a simplified HEOM approach to further reduce the memory cost for practical calculations.
Collapse
Affiliation(s)
- Lu Han
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hou-Dao Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiao Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Laboratory for Physical Sciences at the Microscale and iChEM, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
42
|
Kramer T, Noack M, Reinefeld A, Rodríguez M, Zelinskyy Y. Efficient calculation of open quantum system dynamics and time-resolved spectroscopy with distributed memory HEOM (DM-HEOM). J Comput Chem 2018; 39:1779-1794. [DOI: 10.1002/jcc.25354] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2018] [Revised: 04/20/2018] [Accepted: 04/22/2018] [Indexed: 12/17/2022]
Affiliation(s)
- Tobias Kramer
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
- Department of Physics; Harvard University, 17 Oxford Street; Cambridge Massachusetts 02138
| | - Matthias Noack
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
| | | | - Mirta Rodríguez
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
| | | |
Collapse
|
43
|
Richter M, Fingerhut BP. Coarse-grained representation of the quasi adiabatic propagator path integral for the treatment of non-Markovian long-time bath memory. J Chem Phys 2018; 146:214101. [PMID: 28576089 DOI: 10.1063/1.4984075] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The description of non-Markovian effects imposed by low frequency bath modes poses a persistent challenge for path integral based approaches like the iterative quasi-adiabatic propagator path integral (iQUAPI) method. We present a novel approximate method, termed mask assisted coarse graining of influence coefficients (MACGIC)-iQUAPI, that offers appealing computational savings due to substantial reduction of considered path segments for propagation. The method relies on an efficient path segment merging procedure via an intermediate coarse grained representation of Feynman-Vernon influence coefficients that exploits physical properties of system decoherence. The MACGIC-iQUAPI method allows us to access the regime of biological significant long-time bath memory on the order of hundred propagation time steps while retaining convergence to iQUAPI results. Numerical performance is demonstrated for a set of benchmark problems that cover bath assisted long range electron transfer, the transition from coherent to incoherent dynamics in a prototypical molecular dimer and excitation energy transfer in a 24-state model of the Fenna-Matthews-Olson trimer complex where in all cases excellent agreement with numerically exact reference data is obtained.
Collapse
Affiliation(s)
- Martin Richter
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| |
Collapse
|
44
|
Megier N, Strunz WT, Viviescas C, Luoma K. Parametrization and Optimization of Gaussian Non-Markovian Unravelings for Open Quantum Dynamics. PHYSICAL REVIEW LETTERS 2018; 120:150402. [PMID: 29756890 DOI: 10.1103/physrevlett.120.150402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/24/2017] [Indexed: 06/08/2023]
Abstract
We derive a family of Gaussian non-Markovian stochastic Schrödinger equations for the dynamics of open quantum systems. The different unravelings correspond to different choices of squeezed coherent states, reflecting different measurement schemes on the environment. Consequently, we are able to give a single shot measurement interpretation for the stochastic states and microscopic expressions for the noise correlations of the Gaussian process. By construction, the reduced dynamics of the open system does not depend on the squeezing parameters. They determine the non-Hermitian Gaussian correlation, a wide range of which are compatible with the Markov limit. We demonstrate the versatility of our results for quantum information tasks in the non-Markovian regime. In particular, by optimizing the squeezing parameters, we can tailor unravelings for improving entanglement bounds or for environment-assisted entanglement protection.
Collapse
Affiliation(s)
- Nina Megier
- Institut für Theoretische Physik, Technische Universität Dresden, D-01062, Dresden, Germany
| | - Walter T Strunz
- Institut für Theoretische Physik, Technische Universität Dresden, D-01062, Dresden, Germany
| | - Carlos Viviescas
- Departamento de Física, Universidad Nacional de Colombia, Carrera 30 No. 45-03, Bogota D.C., Colombia
| | - Kimmo Luoma
- Institut für Theoretische Physik, Technische Universität Dresden, D-01062, Dresden, Germany
| |
Collapse
|
45
|
Blau SM, Bennett DIG, Kreisbeck C, Scholes GD, Aspuru-Guzik A. Local protein solvation drives direct down-conversion in phycobiliprotein PC645 via incoherent vibronic transport. Proc Natl Acad Sci U S A 2018; 115:E3342-E3350. [PMID: 29588417 PMCID: PMC5899487 DOI: 10.1073/pnas.1800370115] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The mechanisms controlling excitation energy transport (EET) in light-harvesting complexes remain controversial. Following the observation of long-lived beats in 2D electronic spectroscopy of PC645, vibronic coherence, the delocalization of excited states between pigments supported by a resonant vibration, has been proposed to enable direct excitation transport from the highest-energy to the lowest-energy pigments, bypassing a collection of intermediate states. Here, we instead show that for phycobiliprotein PC645 an incoherent vibronic transport mechanism is at play. We quantify the solvation dynamics of individual pigments using ab initio quantum mechanics/molecular mechanics (QM/MM) nuclear dynamics. Our atomistic spectral densities reproduce experimental observations ranging from absorption and fluorescence spectra to the timescales and selectivity of down-conversion observed in transient absorption measurements. We construct a general model for vibronic dimers and establish the parameter regimes of coherent and incoherent vibronic transport. We demonstrate that direct down-conversion in PC645 proceeds incoherently, enhanced by large reorganization energies and a broad collection of high-frequency vibrations. We suggest that a similar incoherent mechanism is appropriate across phycobiliproteins and represents a potential design principle for nanoscale control of EET.
Collapse
Affiliation(s)
- Samuel M Blau
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
| | - Doran I G Bennett
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
- Bio-Inspired Solar Energy Program, Canadian Institute for Advanced Research, Toronto, ON M5G 1Z8, Canada
| | - Christoph Kreisbeck
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
| | - Gregory D Scholes
- Bio-Inspired Solar Energy Program, Canadian Institute for Advanced Research, Toronto, ON M5G 1Z8, Canada
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138;
- Bio-Inspired Solar Energy Program, Canadian Institute for Advanced Research, Toronto, ON M5G 1Z8, Canada
| |
Collapse
|
46
|
Zhang PP, Bentley CDB, Eisfeld A. Flexible scheme to truncate the hierarchy of pure states. J Chem Phys 2018; 148:134103. [DOI: 10.1063/1.5022225] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P.-P. Zhang
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - C. D. B. Bentley
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - A. Eisfeld
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| |
Collapse
|
47
|
Affiliation(s)
- Hou-Dao Zhang
- Department of Chemical Physics and Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Rui-Xue Xu
- Department of Chemical Physics and Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiao Zheng
- Department of Chemical Physics and Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Department of Chemical Physics and Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
- iChEM (Collaborative Innovation Center of Chemistry for Energy Materials), University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
48
|
Jansen TLC. Simple Quantum Dynamics with Thermalization. J Phys Chem A 2018; 122:172-183. [PMID: 29199829 PMCID: PMC5770886 DOI: 10.1021/acs.jpca.7b10380] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2017] [Revised: 12/04/2017] [Indexed: 02/05/2023]
Abstract
In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.
Collapse
Affiliation(s)
- Thomas L. C. Jansen
- Zernike Institute for Advanced
Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
49
|
Provazza J, Segatta F, Garavelli M, Coker DF. Semiclassical Path Integral Calculation of Nonlinear Optical Spectroscopy. J Chem Theory Comput 2018; 14:856-866. [DOI: 10.1021/acs.jctc.7b01063] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Justin Provazza
- Department
of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Francesco Segatta
- European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK), 38123 Trento, Italy
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - David F. Coker
- Department
of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
50
|
Häse F, Kreisbeck C, Aspuru-Guzik A. Machine learning for quantum dynamics: deep learning of excitation energy transfer properties. Chem Sci 2017; 8:8419-8426. [PMID: 29619189 PMCID: PMC5863613 DOI: 10.1039/c7sc03542j] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2017] [Accepted: 10/23/2017] [Indexed: 12/20/2022] Open
Abstract
Understanding the relationship between the structure of light-harvesting systems and their excitation energy transfer properties is of fundamental importance in many applications including the development of next generation photovoltaics. Natural light harvesting in photosynthesis shows remarkable excitation energy transfer properties, which suggests that pigment-protein complexes could serve as blueprints for the design of nature inspired devices. Mechanistic insights into energy transport dynamics can be gained by leveraging numerically involved propagation schemes such as the hierarchical equations of motion (HEOM). Solving these equations, however, is computationally costly due to the adverse scaling with the number of pigments. Therefore virtual high-throughput screening, which has become a powerful tool in material discovery, is less readily applicable for the search of novel excitonic devices. We propose the use of artificial neural networks to bypass the computational limitations of established techniques for exploring the structure-dynamics relation in excitonic systems. Once trained, our neural networks reduce computational costs by several orders of magnitudes. Our predicted transfer times and transfer efficiencies exhibit similar or even higher accuracies than frequently used approximate methods such as secular Redfield theory.
Collapse
Affiliation(s)
- Florian Häse
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; ; Tel: +1-617-384-8188
| | - Christoph Kreisbeck
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; ; Tel: +1-617-384-8188
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; ; Tel: +1-617-384-8188
| |
Collapse
|