1
|
Zhong K, Nguyen HL, Do TN, Tan HS, Knoester J, Jansen TLC. Coarse-Grained Approach to Simulate Signatures of Excitation Energy Transfer in Two-Dimensional Electronic Spectroscopy of Large Molecular Systems. J Chem Theory Comput 2024; 20:6111-6124. [PMID: 38996082 PMCID: PMC11270824 DOI: 10.1021/acs.jctc.4c00413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Revised: 06/25/2024] [Accepted: 06/27/2024] [Indexed: 07/14/2024]
Abstract
Two-dimensional electronic spectroscopy (2DES) has proven to be a highly effective technique in studying the properties of excited states and the process of excitation energy transfer in complex molecular assemblies, particularly in biological light-harvesting systems. However, the accurate simulation of 2DES for large systems still poses a challenge because of the heavy computational demands it entails. In an effort to overcome this limitation, we devised a coarse-grained 2DES method. This method encompasses the treatment of the entire system by dividing it into distinct weakly coupled segments, which are assumed to communicate predominantly through incoherent exciton transfer. We first demonstrate the efficiency of this method through simulation on a model dimer system, which demonstrates a marked improvement in calculation efficiency, with results that exhibit good concordance with reference spectra calculated with less approximate methods. Additionally, the application of this method to the light-harvesting antenna 2 (LH2) complex of purple bacteria showcases its advantages, accuracy, and limitations. Furthermore, simulating the anisotropy decay in LH2 induced by energy transfer and its comparison with experiments confirm that the method is capable of accurately describing dynamical processes in a biologically relevant system. This method presented lends itself to an extension that accounts for the effect of intrasegment relaxation processes on the 2DES spectra, which for computational efficiency are ignored in the implementation reported here. It is envisioned that the method will be employed in the future to accurately and efficiently calculate 2D spectra of more extensive systems, such as photosynthetic supercomplexes.
Collapse
Affiliation(s)
- Kai Zhong
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 3, 9747 AG Groningen, The Netherlands
- School
of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, 21 Nanyang Link, 637371 Singapore
| | - Hoang Long Nguyen
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 3, 9747 AG Groningen, The Netherlands
- School
of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, 21 Nanyang Link, 637371 Singapore
| | - Thanh Nhut Do
- School
of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, 21 Nanyang Link, 637371 Singapore
| | - Howe-Siang Tan
- School
of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, 21 Nanyang Link, 637371 Singapore
| | - Jasper Knoester
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 3, 9747 AG Groningen, The Netherlands
- Faculty
of Science, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
| | - Thomas L. C. Jansen
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 3, 9747 AG Groningen, The Netherlands
| |
Collapse
|
2
|
Bose A. Incorporation of Empirical Gain and Loss Mechanisms in Open Quantum Systems through Path Integral Lindblad Dynamics. J Phys Chem Lett 2024:3363-3368. [PMID: 38498435 DOI: 10.1021/acs.jpclett.4c00489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Path integrals offer a robust approach for simulating open quantum dynamics with advancements transcending initial system size limitations. However, accurately modeling systems governed by mechanisms that do not conserve the number of quantum particles, such as lossy cavity modes, remains a challenge. We present a method to incorporate such empirical source and drain mechanisms within a path integral framework using quantum master equations. This technique facilitates rigorous inclusion of bath degrees of freedom while accommodating empirical time scales via Lindbladian dynamics. Computational costs are primarily driven by the path integral method with minimal overhead from Lindbladian terms. We use it to study exciton transport in a four-site Fenna-Matthews-Olson model, examining the potential loss of the exciton to the reaction center. This path integral Lindblad method promises an enhanced ability to simulate dynamics and will be fundamental to simulation of spectra in diverse quantum processes in open systems.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Mumbai 400005, India
| |
Collapse
|
3
|
Bose A. Quantum correlation functions through tensor network path integral. J Chem Phys 2023; 159:214110. [PMID: 38051096 DOI: 10.1063/5.0174338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Accepted: 11/07/2023] [Indexed: 12/07/2023] Open
Abstract
Tensor networks have historically proven to be of great utility in providing compressed representations of wave functions that can be used for the calculation of eigenstates. Recently, it has been shown that a variety of these networks can be leveraged to make real time non-equilibrium simulations of dynamics involving the Feynman-Vernon influence functional more efficient. In this work, a tensor network is developed for non-perturbatively calculating the equilibrium correlation function for open quantum systems using the path integral methodology. These correlation functions are of fundamental importance in calculations of rates of reactions, simulations of response functions and susceptibilities, spectra of systems, etc. The influence of the solvent on the quantum system is incorporated through an influence functional, whose unconventional structure motivates the design of a new optimal matrix product-like operator that can be applied to the so-called path amplitude matrix product state. This complex-time tensor network path integral approach provides an exceptionally efficient representation of the path integral, enabling simulations for larger systems strongly interacting with baths and at lower temperatures out to longer time. The derivation, design, and implementation of this method are discussed along with a wide range of illustrations ranging from rate theory and symmetrized spin correlation functions to simulation of response of the Fenna-Matthews-Olson complex to light.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Mumbai 400005, India
| |
Collapse
|
4
|
Wang G, Cai Z. Real-Time Simulation of Open Quantum Spin Chains with the Inchworm Method. J Chem Theory Comput 2023. [PMID: 38039552 DOI: 10.1021/acs.jctc.3c00751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2023]
Abstract
We study the real-time simulation of open quantum systems, where the system is modeled by a spin chain with each spin associated with its own harmonic bath. Our method couples the inchworm method for the spin-boson model and the modular path integral methodology for spin systems. In particular, the introduction of the inchworm method can significantly suppress the numerical sign problem. Both methods are tweaked to make them work seamlessly with each other. We represent our approach in the language of diagrammatic methods and analyze the asymptotic behavior of the computational cost. Extensive numerical experiments are performed to validate our method.
Collapse
Affiliation(s)
- Geshuo Wang
- Department of Mathematics, National University of Singapore, Block S17, 10 Lower Kent Ridge Road, Singapore 119076, Singapore
| | - Zhenning Cai
- Department of Mathematics, National University of Singapore, Block S17, 10 Lower Kent Ridge Road, Singapore 119076, Singapore
| |
Collapse
|
5
|
Bose A, Walters PL. Impact of Spatial Inhomogeneity on Excitation Energy Transport in the Fenna-Matthews-Olson Complex. J Phys Chem B 2023; 127:7663-7673. [PMID: 37647510 DOI: 10.1021/acs.jpcb.3c03062] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
The dynamics of the excitation energy transfer (EET) in photosynthetic complexes is an interesting question both from the perspective of fundamental understanding and the research in artificial photosynthesis. Over the past decade, very accurate spectral densities have been developed to capture spatial inhomogeneities in the Fenna-Matthews-Olson (FMO) complex. However, challenges persist in numerically simulating these systems, both in terms of parameterizing them and following their dynamics over long periods of time because of long non-Markovian memories. We investigate the dynamics of FMO with the exact treatment of various theoretical spectral densities using the new tensor network path integral-based methods, which are uniquely capable of addressing the difficulty of long memory length and incoherent Förster theory. It is also important to be able to analyze the pathway of EET flow, which can be difficult to identify given the non-trivial structure of connections between bacteriochlorophyll molecules in FMO. We use the recently introduced ideas of relating coherence to population derivatives to analyze the transport process and reveal some new routes of transport. The combination of exact and approximate methods sheds light on the role of coherences in affecting the fine details of the transport and promises to be a powerful toolbox for future exploration of other open systems with quantum transport.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Mumbai 400005, India
| | - Peter L Walters
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Miller Institute for Basic Research in Science, University of California Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
6
|
Bose A, Walters PL. Impact of Solvent on State-to-State Population Transport in Multistate Systems Using Coherences. J Chem Theory Comput 2023. [PMID: 37466459 DOI: 10.1021/acs.jctc.3c00200] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Understanding the pathways taken by a quantum particle during a transport process is an enormous challenge. There are broadly two different aspects of the problem that affect the route taken. First is obviously the couplings between the various sites, which translates into the intrinsic "strength" of a state-to-state channel. Apart from these inter-state couplings, the relative coupling strengths and timescales of the solvent modes form the second factor. This impact of the dissipative environment is significantly more difficult to analyze. Building on the recently derived relations between coherences and population derivatives, we present an analysis of the transport that allows us to account for both the effects in a rigorous manner. We demonstrate the richness hidden behind the transport even for a relatively simple system, a 4-site coarse-grained model of the Fenna-Matthews-Olson complex. The effect of the local dissipative media is highly nontrivial. We show that while the impact on the total site population may be small, there are noticeable changes to the pathway taken by the transport process. We also demonstrate how an analysis in a similar spirit can be done using the Förster approximation. The ability to untangle the dynamics at a greater granularity opens up possibilities in terms of design of novel systems with an eye toward quantum control.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Mumbai 400005, India
| | - Peter L Walters
- Department of Chemistry and Biochemistry, George Mason University, Fairfax, Virginia 22030, United States
| |
Collapse
|
7
|
Maity S, Kleinekathöfer U. Recent progress in atomistic modeling of light-harvesting complexes: a mini review. PHOTOSYNTHESIS RESEARCH 2023; 156:147-162. [PMID: 36207489 PMCID: PMC10070314 DOI: 10.1007/s11120-022-00969-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 09/22/2022] [Indexed: 06/16/2023]
Abstract
In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.
Collapse
Affiliation(s)
- Sayan Maity
- 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
|
8
|
Bose A, Walters PL. Effect of temperature gradient on quantum transport. Phys Chem Chem Phys 2022; 24:22431-22436. [PMID: 36102915 DOI: 10.1039/d2cp03030f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The recently introduced multisite tensor network path integral (MS-TNPI) method [Bose and Walters, J. Chem. Phys., 2022, 156, 24101] for simulating quantum dynamics of extended systems has been shown to be effective in studying one-dimensional systems coupled with local baths. Quantum transport in these systems is typically studied at a constant temperature. However, temperature seems to be a very obvious parameter that can be spatially changed to control this transport. Here, MS-TNPI is used to study the "non-equilibrium" effects of an externally imposed temperature profile on the excitonic transport in one-dimensional Frenkel chains coupled with local vibrations. We show that in addition to being important for incorporating heating effects of excitation by lasers, temperature can also be an interesting parameter for quantum control.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA.
| | - Peter L Walters
- Department of Chemistry, University of California, Berkeley, California 94720, USA.,Miller Institute for Basic Research in Science, University of California Berkeley, Berkeley, California 94720, USA.
| |
Collapse
|
9
|
Bose A. Zero-Cost Corrections to Influence Functional Coefficients from Bath Response Functions. J Chem Phys 2022; 157:054107. [DOI: 10.1063/5.0101396] [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/15/2022] Open
Abstract
Recent work has shown that it is possible to circumvent the calculation of the spectral density and directly calculate the coefficients of the discretized influence functionals using data from classical trajectory simulations. However, the accuracy of this procedure depends on the validity of the high temperature approximation. In this work, an alternative derivation based on the Kubo formalism is provided. This enables the calculation of additional correction terms that increases the range of applicability of the procedure to lower temperatures. Because it is based on the Kubo-transformed correlation function, this approach enables the direct use of correlation functions obtained from methods like ring-polymer molecular dynamics and centroid molecular dynamics in determining the influence functional coefficients for subsequent system-solvent simulations. The accuracy of the original procedure and the corrected procedure is investigated across a range of parameters. It is interesting that the correction term comes at zero additional cost. Furthermore, it is possible to improve upon the correction using zero-cost physical intuition and heuristics making the method even more accurate.
Collapse
Affiliation(s)
- Amartya Bose
- Chemistry, Princeton University Department of Chemistry, United States of America
| |
Collapse
|