1
|
Le Dé B, Jaouadi A, Mangaud E, Chin AW, Desouter-Lecomte M. Managing temperature in open quantum systems strongly coupled with structured environments. J Chem Phys 2024; 160:244102. [PMID: 38913841 DOI: 10.1063/5.0214051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 06/06/2024] [Indexed: 06/26/2024] Open
Abstract
In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.
Collapse
Affiliation(s)
- Brieuc Le Dé
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | - Amine Jaouadi
- LyRIDS, ECE Paris, Graduate School of Engineering, Paris F-75015, France
| | - Etienne Mangaud
- MSME, Université Gustave Eiffel, UPEC, CNRS, F-77454 Marne-La-Vallée, France
| | - Alex W Chin
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, F-75005 Paris, France
| | | |
Collapse
|
2
|
Hino K, Kurashige Y. Encoding a Many-Body Potential Energy Surface into a Grid-Based Matrix Product Operator. J Chem Theory Comput 2024; 20:3839-3849. [PMID: 38647101 DOI: 10.1021/acs.jctc.4c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
An efficient algorithm for compressing a given many-body potential energy surface (PES) of molecular systems into a grid-based matrix product operator (MPO) is proposed. The PES is once represented by a full-dimensional or truncated many-body expansion form, which is obtained by ab initio calculations at each grid mesh point, and then all terms in the expansion are compressed and merged into a single MPO while maintaining the bond dimension of the MPO as small as possible. It was shown that the ab initio PES of the H2CO was compressed by more than 2 orders of magnitude in the size of the site operators without loss of accuracy. By the use of grid basis, the tensor rank of the site operators of the MPO is reduced from four to three due to the diagonal nature of the position-dependent operators on grid basis, which significantly reduces the computational cost of the tensor contractions required in the real and imaginary time evolution of the matrix product state (MPS) wave functions with the grid-based MPO (Grid-MPO) Hamiltonian. Similar to other grid-based methods, Grid-MPO is easily applicable to any kinds of potentials of molecular systems, such as analytical empirical model potentials expressed by position operators and ab initio potentials, if the values at the grid points are available. Using the Grid-MPO combined with the MPS, we calculated the time correlation function of the Eigen cation H 3 O + ( H 2 O ) 3 to predict the infrared spectrum and compared with the experimental and the previous theoretical studies. The actual scaling with the size of systems was examined for the multidimensional Henon-Heiles Hamiltonian. It was shown that the method is considerably accelerated by the graphic processing unit (GPU) because the sizes of site operators were kept small and all tensors were able to be stored on the VRAM of a GPU.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
- FOREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
- CREST, JST, Honcho 4-1-8, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
3
|
Xiang C, Jia W, Fang WH, Li Z. Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase. J Chem Theory Comput 2024; 20:775-786. [PMID: 38198503 DOI: 10.1021/acs.jctc.3c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).
Collapse
Affiliation(s)
- Chunyang Xiang
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Weile Jia
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
4
|
Li G, Shi Z, Guo X, Wang L. What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching. J Phys Chem Lett 2023; 14:9855-9863. [PMID: 37890155 DOI: 10.1021/acs.jpclett.3c02690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/29/2023]
Abstract
When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.
Collapse
Affiliation(s)
- Guijie Li
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Zhecun Shi
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
5
|
Yang H, Li W, Ren J, Shuai Z. Time-Dependent Density Matrix Renormalization Group Method for Quantum Transport with Phonon Coupling in Molecular Junction. J Chem Theory Comput 2023; 19:6070-6081. [PMID: 37669099 DOI: 10.1021/acs.jctc.3c00340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.
Collapse
Affiliation(s)
- Hengrui Yang
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, Guangdong 518172, People's Republic of China
| |
Collapse
|
6
|
Xu Y, Cheng Y, Song Y, Ma H. New Density Matrix Renormalization Group Approaches for Strongly Correlated Systems Coupled with Large Environments. J Chem Theory Comput 2023. [PMID: 37471519 DOI: 10.1021/acs.jctc.2c01316] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
Thanks to the high compression of the matrix product state (MPS) form of the wave function and the efficient site-by-site iterative sweeping optimization algorithm, the density matrix normalization group (DMRG) and its time-dependent variant (TD-DMRG) have been established as powerful computational tools in accurately simulating the electronic structure and quantum dynamics of strongly correlated molecules with a large number (101-2) of quantum degrees of freedom (active orbitals or vibrational modes). However, the quantitative characterization of the quantum many-body behaviors of realistic strongly correlated systems requires a further consideration of the interaction between the embedded active subsystem and the remaining correlated environment, e.g., a larger number (102-3) of external orbitals in electronic structure or infinite condensed-phase phononic modes in nucleus dynamics. To this end, we introduced three new post-DMRG and TD-DMRG approaches, namely (1) DMRG2sCI-MRCI and DMRG2sCI-ENPT by the reconstruction of selected configuration interaction (sCI) type of compact reference function from DMRG coefficients and the use of externally contracted MRCI (multireference configuration interaction) and Epstein-Nesbet perturbation theory (ENPT), without recourse to the expensive high order n-electron reduced density matrices (n-RDMs). (2) DMRG combined with RR-MRCI (renormalized residue-based MRCI), which improves the computational accuracy and efficiency of internally contracted (ic) MRCI by renormalizing the contracted bases with small-sized buffer environment(s) of a few external orbitals as probes based on quantum information theory. (3) HM (hierarchical mapping)-TD-DMRG in which a large environment is reduced to a small number of renormalized environmental modes (which accounts for the most vital system-environment interactions) through stepwise mapping transformation. These advances extend the efficacy of highly accurate DMRG/TD-DMRG computations to the quantitative characterization of the electronic structure and quantum dynamics in realistic strongly correlated systems coupled with large environments and are reviewed in this paper.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| |
Collapse
|
7
|
Lyu N, Mulvihill E, Soley MB, Geva E, Batista VS. Tensor-Train Thermo-Field Memory Kernels for Generalized Quantum Master Equations. J Chem Theory Comput 2023; 19:1111-1129. [PMID: 36719350 DOI: 10.1021/acs.jctc.2c00892] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The generalized quantum master equation (GQME) approach provides a rigorous framework for deriving the exact equation of motion for any subset of electronic reduced density matrix elements (e.g., the diagonal elements). In the context of electronic dynamics, the memory kernel and inhomogeneous term of the GQME introduce the implicit coupling to nuclear motion and dynamics of electronic density matrix elements that are projected out (e.g., the off-diagonal elements), allowing for efficient quantum dynamics simulations. Here, we focus on benchmark quantum simulations of electronic dynamics in a spin-boson model system described by various types of GQMEs. Exact memory kernels and inhomogeneous terms are obtained from short-time quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) simulations and are compared with those obtained from an approximate linearized semiclassical method, allowing for assessment of the accuracy of these approximate memory kernels and inhomogeneous terms. Moreover, we have analyzed the computational cost of the full and reduced-dimensionality GQMEs. The scaling of the computational cost is dependent on several factors, sometimes with opposite scaling trends. The TT-TFD memory kernels can provide insights on the main sources of inaccuracies of GQME approaches when combined with approximate input methods and pave the road for the development of quantum circuits that implement GQMEs on digital quantum computers.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Ellen Mulvihill
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States.,Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
8
|
Komarova K. Density Matrix via Few Dominant Observables for the Ultrafast Non-Radiative Decay in Pyrazine. J Chem Theory Comput 2023; 19:746-757. [PMID: 36657738 PMCID: PMC11137821 DOI: 10.1021/acs.jctc.2c01211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Indexed: 01/21/2023]
Abstract
Unraveling the density matrix of a non-stationary quantum state as an explicit function of a few observables provides a complementary view of quantum dynamics. We have recently developed a practical way to identify the minimal set of the dominant observables that govern the quantal dynamics even in the case of strong non-adiabatic effects and large anharmonicity [Komarova et al., J. Chem. Phys. 155, 204110 (2021)]. Fast convergence in the number of the dominant contributions is achieved when instead of the density matrix we describe the time-evolution of the surprisal, the logarithm of the density operator. In the present work, we illustrate the efficiency of the proposed approach using an example of the early time dynamics in pyrazine in a Hilbert space accounting for up to four vibrational normal modes, {Q10a, Q6a, Q1, and Q9a}, and two coupled electronic states, the optically dark B 1 3 u ( n π * ) and the bright B 1 2 u ( π π * ) states. Dynamics in four-dimensional (4D) configurational space involve 19,600 vibronic eigenstates. Our results reveal that the rate of the ultrafast population decay as well as the shape of the nuclear wave packets in 2D, accounting only for {Q10a,Q6a} normal modes, are accurately captured with only six dominant time-independent observables in the surprisal. Extension of the dynamics to 3D and 4D vibrational subspace requires only five additional constraints. The time-evolution of a quantum state in 4D vibrational space on two electronic states is thus compacted to only 11 time-dependent coefficients of these observables.
Collapse
Affiliation(s)
- Ksenia Komarova
- The Fritz Haber Center for Molecular
Dynamics and Institute of Chemistry, The
Hebrew University of Jerusalem, Jerusalem91904, Israel
| |
Collapse
|
9
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis. J Chem Theory Comput 2022; 18:5837-5855. [DOI: 10.1021/acs.jctc.2c00702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| |
Collapse
|
10
|
Ge Y, Li W, Ren J, Shuai Z. Computational Method for Evaluating the Thermoelectric Power Factor for Organic Materials Modeled by the Holstein Model: A Time-Dependent Density Matrix Renormalization Group Formalism. J Chem Theory Comput 2022; 18:6437-6446. [PMID: 36174220 DOI: 10.1021/acs.jctc.2c00651] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining the Kubo formula with the finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF = α2 σ, where α is the Seebeck coefficient and σ is the electrical conductivity, and apply the algorithm to Holstein Hamiltonian with electron-phonon coupling to model organic materials. Our algorithm can provide a unified description covering the weak coupling limit described by the bandlike Boltzmann transport to the strong coupling hopping limit.
Collapse
Affiliation(s)
- Yufei Ge
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China
| | - Jiajun Ren
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, 100875Beijing, P. R. China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084Beijing, P. R. China.,School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen518172, Guangdong, P. R. China
| |
Collapse
|
11
|
Wang Y, Ren J, Li W, Shuai Z. Hybrid Quantum-Classical Boson Sampling Algorithm for Molecular Vibrationally Resolved Electronic Spectroscopy with Duschinsky Rotation and Anharmonicity. J Phys Chem Lett 2022; 13:6391-6399. [PMID: 35802770 DOI: 10.1021/acs.jpclett.2c01475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Using a photonic quantum computer for boson sampling has demonstrated a tremendous advantage over classical supercomputers. It is highly desirable to develop boson sampling algorithms for realistic scientific problems. In this work, we propose a hybrid quantum-classical sampling (HQCS) algorithm to calculate the optical spectrum for complex molecules considering Duschinsky rotation effects and anharmonicity. The classical sum-over-states method for this problem has a computational complexity that exponentially increases with system size. The HQCS algorithm creates an intermediate harmonic potential energy surface (PES) to bridge the initial and final PESs. The magnitude and sign of the overlap between the initial and the intermediate state are estimated by boson sampling and classical algorithms, respectively. The overlap between the intermediate and the final state is efficiently evaluated by classical algorithms. The feasibility of HQCS is demonstrated in calculations of the emission spectrum of a Morse model as well as the pyridine molecule.
Collapse
Affiliation(s)
- Yuanheng Wang
- 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
| | - Weitang Li
- 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
|
12
|
Lyu N, Soley MB, Batista VS. Tensor-Train Split-Operator KSL (TT-SOKSL) Method for Quantum Dynamics Simulations. J Chem Theory Comput 2022; 18:3327-3346. [PMID: 35649210 DOI: 10.1021/acs.jctc.2c00209] [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/30/2022]
Abstract
Numerically exact simulations of quantum reaction dynamics, including nonadiabatic effects in excited electronic states, are essential to gain fundamental insights into ultrafast chemical reactivity and rigorous interpretations of molecular spectroscopy. Here, we introduce the tensor-train split-operator KSL (TT-SOKSL) method for quantum simulations in tensor-train (TT)/matrix product state (MPS) representations. TT-SOKSL propagates the quantum state as a tensor train using the Trotter expansion of the time-evolution operator, as in the tensor-train split-operator Fourier transform (TT-SOFT) method. However, the exponential operators of the Trotter expansion are applied using a rank-adaptive TT-KSL scheme instead of using the scaling and squaring approach as in TT-SOFT. We demonstrate the accuracy and efficiency of TT-SOKSL as applied to simulations of the photoisomerization of the retinal chromophore in rhodopsin, including nonadiabatic dynamics at a conical intersection of potential energy surfaces. The quantum evolution is described in full dimensionality by a time-dependent wavepacket evolving according to a two-state 25-dimensional model Hamiltonian. We find that TT-SOKSL converges faster than TT-SOFT with respect to the maximally allowed memory requirement of the tensor-train representation and better preserves the norm of the time-evolving state. When compared to the corresponding simulations based on the TT-KSL method, TT-SOKSL has the advantage of avoiding the need to construct the matrix product state Laplacian by exploiting the linear scaling of multidimensional tensor-train Fourier transforms.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, United States.,Yale Quantum Institute, Yale University, P.O. Box 208334, New Haven, Connecticut 06520-8263, United States
| |
Collapse
|
13
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
14
|
Li W, Ren J, Yang H, Shuai Z. On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:254003. [PMID: 35378514 DOI: 10.1088/1361-648x/ac640e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2022] [Accepted: 04/04/2022] [Indexed: 06/14/2023]
Abstract
Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N2molecule, and the S1/S2internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Hengrui Yang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| |
Collapse
|
15
|
Ren J, Li W, Jiang T, Wang Y, Shuai Z. Time‐dependent density matrix renormalization group method for quantum dynamics in complex systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| |
Collapse
|
16
|
Xu Y, Xie Z, Xie X, Schollwöck U, Ma H. Stochastic Adaptive Single-Site Time-Dependent Variational Principle. JACS AU 2022; 2:335-340. [PMID: 35252984 PMCID: PMC8889605 DOI: 10.1021/jacsau.1c00474] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Indexed: 06/14/2023]
Abstract
In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the traditional advantages of both the high efficiency of the single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random-matrix theory, additional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.
Collapse
Affiliation(s)
- Yihe Xu
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| | - Xiaoyu Xie
- Department
of Chemistry, University of Liverpool, Liverpool L69 3BX, United Kingdom
| | - Ulrich Schollwöck
- Arnold
Sommerfeld Center of Theoretical Physics, Department of Physics, University of Munich, Theresienstrasse 37, 80333 Munich, Germany
- Munich
Center for Quantum Science and Technology (MCQST), 80799 Munich, Germany
| | - Haibo Ma
- School
of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle
Emissions Control, Nanjing University, Nanjing 210023, China
| |
Collapse
|
17
|
Baiardi A, Grimmel SA, Steiner M, Türtscher PL, Unsleber JP, Weymuth T, Reiher M. Expansive Quantum Mechanical Exploration of Chemical Reaction Paths. Acc Chem Res 2022; 55:35-43. [PMID: 34918903 DOI: 10.1021/acs.accounts.1c00472] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Quantum mechanical methods have been well-established for the elucidation of reaction paths of chemical processes and for the explicit dynamics of molecular systems. While they are usually deployed in routine manual calculations on reactions for which some insights are already available (typically from experiment), new algorithms and continuously increasing capabilities of modern computer hardware allow for exploratory open-ended computational campaigns that are unbiased and therefore enable unexpected discoveries. Highly efficient and even automated procedures facilitate systematic approaches toward the exploration of uncharted territory in molecular transformations and dynamics. In this work, we elaborate on such explorative approaches that range from reaction network explorations with (stationary) quantum chemical methods to explorative molecular dynamics and migrant wave packet dynamics. The focus is on recent developments that cover the following strategies. (i) Pruning search options for elementary reaction steps by heuristic rules based on the first-principles of quantum mechanics: Rules are required for reducing the combinatorial explosion of potentially reactive atom pairings, and rooting them in concepts derived from the electronic wave function makes them applicable to any molecular system. (ii) Enforcing reactive events by external biases: Inducing a reaction requires constraints that steer and direct elementary-step searches, which can be formulated in terms of forces, velocities, or supplementary potentials. (iii) Manual steering facilitated by interactive quantum mechanics: As ultrafast quantum chemical methods allow for real-time manual interactions with molecular systems, human-intuition-guided paths can be easily explored with suitable human-machine interfaces. (iv) New approaches for transition-state optimization with continuous curve representations can provide stable schemes to be driven in an automated way by allowing for an efficient tuning of the curve's parameters (instead of a manipulation of a collection of structures along the path), and (v) reactive molecular dynamics and direct wave packet propagation exploit the equations of motion of an underlying mechanical theory (usually, classical Newtonian mechanics or Schrödinger quantum mechanics). Explorative approaches are likely to replace the current state of the art in computational chemistry, because they reduce the human effort to be invested in reaction path elucidations, they are less prone to errors and bias-free, and they cover more extensive regions of the relevant configuration space. As a result, computational investigations that rely on these techniques are more likely to deliver surprising discoveries.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A. Grimmel
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L. Türtscher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P. Unsleber
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
18
|
Komarova K, Remacle F, Levine RD. Compacting the density matrix in quantum dynamics: Singular value decomposition of the surprisal and the dominant constraints for anharmonic systems. J Chem Phys 2021; 155:204110. [PMID: 34852474 DOI: 10.1063/5.0072351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a practical method for compacting the time evolution of the quantum state of a closed physical system. The density matrix is specified as a function of a few time-independent observables where their coefficients are time-dependent. The key mathematical step is the vectorization of the surprisal, the logarithm of the density matrix, at each time point of interest. The time span used depends on the required spectral resolution. The entire course of the system evolution is represented as a matrix where each column is the vectorized surprisal at the given time point. Using the singular value decomposition (SVD) of this matrix, we generate realistic approximations for the time-independent observables and their respective time-dependent coefficients. This allows for a simplification of the algebraic procedure for determining the dominant constraints (the time-independent observables) in the sense of the maximal entropy approach. A non-stationary coherent initial state of a Morse oscillator is used to introduce the approach. We derive the analytical exact expression for the surprisal as a function of time, and this offers a benchmark for comparison with the accurate but approximate SVD results. We discuss two examples of a Morse potential of different anharmonicities, H2 and I2 molecules. We further demonstrate the approach for a two-coupled electronic state problem, the well-studied non-radiative decay of pyrazine from its bright state. Five constraints are found to be enough to capture the ultrafast electronic population exchange and to recover the dynamics of the wave packet in both electronic states.
Collapse
Affiliation(s)
- K Komarova
- The Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Francoise Remacle
- The Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - R D Levine
- The Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| |
Collapse
|
19
|
Ren JJ, Wang YH, Li WT, Jiang T, Shuai ZG. Time-dependent density matrix renormalization group coupled with n-mode representation potentials for the excited state radiationless decay rate: Formalism and application to azulene. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2108138] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jia-jun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Yuan-heng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Wei-tang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Zhi-gang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| |
Collapse
|
20
|
Jiang T, Ren J, Shuai Z. Chebyshev Matrix Product States with Canonical Orthogonalization for Spectral Functions of Many-Body Systems. J Phys Chem Lett 2021; 12:9344-9352. [PMID: 34549961 DOI: 10.1021/acs.jpclett.1c02688] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We propose a method to calculate the spectral functions of many-body systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigenenergies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. For the first time we demonstrate that Chebyshev MPS can be used to deal with ab initio electronic Hamiltonian effectively. We emphasize the strength of coCheMPS to calculate the low excited states of systems with sparse discrete spectrum. We also caution the application for electron-phonon systems with dense density of states.
Collapse
Affiliation(s)
- Tong Jiang
- 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
|
21
|
Li W, Ren J, Shuai Z. A general charge transport picture for organic semiconductors with nonlocal electron-phonon couplings. Nat Commun 2021; 12:4260. [PMID: 34253724 PMCID: PMC8275621 DOI: 10.1038/s41467-021-24520-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Accepted: 06/22/2021] [Indexed: 11/16/2022] Open
Abstract
The nonlocal electron-phonon couplings in organic semiconductors responsible for the fluctuation of intermolecular transfer integrals has been the center of interest recently. Several irreconcilable scenarios coexist for the description of the nonlocal electron-phonon coupling, such as phonon-assisted transport, transient localization, and band-like transport. Through a nearly exact numerical study for the carrier mobility of the Holstein-Peierls model using the matrix product states approach, we locate the phonon-assisted transport, transient localization and band-like regimes as a function of the transfer integral (V) and the nonlocal electron-phonon couplings (ΔV), and their distinct transport behaviors are analyzed by carrier mobility, mean free path, optical conductivity and one-particle spectral function. We also identify an "intermediate regime" where none of the established pictures applies, and the generally perceived hopping regime is found to be at a very limited end in the proposed regime paradigm.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, China.
| |
Collapse
|
22
|
Zhai H, Chan GKL. Low communication high performance ab initio density matrix renormalization group algorithms. J Chem Phys 2021; 154:224116. [DOI: 10.1063/5.0050902] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
23
|
Wang Y, Ren J, Shuai Z. Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG. J Chem Phys 2021; 154:214109. [PMID: 34240969 DOI: 10.1063/5.0052804] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S1 to S0 on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%-40% higher than the rates under HA.
Collapse
Affiliation(s)
- Yuanheng Wang
- 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
|
24
|
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
|
25
|
Peng J, Xie Y, Hu D, Lan Z. Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis. J Chem Phys 2021; 154:094122. [PMID: 33685149 DOI: 10.1063/5.0039743] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.
Collapse
Affiliation(s)
- Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
26
|
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
|
27
|
Shuai Z, Li W, Ren J, Jiang Y, Geng H. Applying Marcus theory to describe the carrier transports in organic semiconductors: Limitations and beyond. J Chem Phys 2020; 153:080902. [PMID: 32872875 DOI: 10.1063/5.0018312] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.
Collapse
Affiliation(s)
- Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, People's Republic of China
| | - Yuqian Jiang
- Laboratory for Nanosystem and Hierarchy Fabrication, National Center for Nanoscience and Technology, Chinese Academy of Sciences, 100084 Beijing, People's Republic of China
| | - Hua Geng
- Department of Chemistry, Capital Normal University, 100048 Beijing, People's Republic of China
| |
Collapse
|
28
|
Ren J, Li W, Jiang T, Shuai Z. A general automatic method for optimal construction of matrix product operators using bipartite graph theory. J Chem Phys 2020; 153:084118. [PMID: 32872857 DOI: 10.1063/5.0018149] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Constructing matrix product operators (MPOs) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For the DMRG to be conveniently used in different problems described by different Hamiltonians, in this work, we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) it is automatic in that only the definition of the operator is required; (ii) it is symbolic thus free of any numerical error; (iii) the complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; and (iv) the symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and the Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with the sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.
Collapse
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Tong Jiang
- 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
|
29
|
Shuai Z. Thermal Vibration Correlation Function Formalism for Molecular Excited State Decay Rates. CHINESE J CHEM 2020. [DOI: 10.1002/cjoc.202000226] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Zhigang Shuai
- MOE Key Laboratory for Organic OptoElectronics and Molecular Enguneering, Department of Chemistry, Tsinghua University 100084 Beijing China
| |
Collapse
|
30
|
Li W, Ren J, Shuai Z. Finite-Temperature TD-DMRG for the Carrier Mobility of Organic Semiconductors. J Phys Chem Lett 2020; 11:4930-4936. [PMID: 32492339 DOI: 10.1021/acs.jpclett.0c01072] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
A large number of nonadiabatic dynamical studies have been applied to reveal the nature of carrier transport in organic semiconductors with different approximations. We present here a "nearly exact" graphical-process-unit-based finite-temperature time-dependent density matrix renormalization group (TD-DMRG) method to evaluate the carrier mobility in organic semiconductors, as described by the electron-phonon model, in particular, in rubrene crystal, one of the prototypical organic semiconductors, with parameters derived from first-principles. We find that (i) TD-DMRG is a general and robust method that can bridge the gap between hopping and band pictures, covering a wide range of electronic coupling strengths and (ii) with realistic parameters, TD-DMRG is able to account for the experimentally observed "band-like" transport behavior (∂μ/∂T < 0) in rubrene. We further study the long-standing puzzle of the isotope effect for charge transport and unambiguously demonstrate that the negative isotope effect (∂μ/∂m < 0 where m is the atomic mass) should be universal.
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
|
Jiang T, Li W, Ren J, Shuai Z. Finite Temperature Dynamical Density Matrix Renormalization Group for Spectroscopy in Frequency Domain. J Phys Chem Lett 2020; 11:3761-3768. [PMID: 32316732 DOI: 10.1021/acs.jpclett.0c00905] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present a novel methodology through casting the dynamical density matrix renormalization group (DDMRG) into the matrix product state (MPS) formulation to calculate the spectroscopy at finite temperature for molecular aggregates. The frequency domain algorithm can avoid the time evolution accumulation of error and is naturally suitable for parallelization, in addition to facile graphic processing unit (GPU) acceleration. The high accuracy is demonstrated by simulating the optical spectra of vibronic model systems ranging from an exactly solvable dimer model to a more complex real-world perylene bisimide (PBI) J-aggregate. The relationship between the 0-0 emission strength and the exciton thermal coherent length is discussed for linearly stacked aggregates. The computing performance largely boosted by GPU demonstrates that DDMRG emerges as a promising method to study dynamical properties for complex systems.
Collapse
Affiliation(s)
- Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of 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
|
32
|
Zheng J, Peng J, Xie Y, Long Y, Ning X, Lan Z. Study of the exciton dynamics in perylene bisimide (PBI) aggregates with symmetrical quasiclassical dynamics based on the Meyer–Miller mapping Hamiltonian. Phys Chem Chem Phys 2020; 22:18192-18204. [DOI: 10.1039/d0cp00648c] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The exciton dynamics in one-dimensional stacked PBI (Perylene Bisimide) aggregates was studied with SQC-MM dynamics (Symmetrical Quasiclassical Dynamics based on the Meyer–Miller mapping Hamiltonian).
Collapse
Affiliation(s)
- Jie Zheng
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - Jiawei Peng
- SCNU Environmental Research Institute
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment
- South China Normal University
- Guangzhou 510006
- China
| | - Yu Xie
- SCNU Environmental Research Institute
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment
- South China Normal University
- Guangzhou 510006
- China
| | - Yunze Long
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - Xin Ning
- Industrial Research Institute of Nonwovens & Technical Textiles
- Shandong Center for Engineered Nonwovens (SCEN)
- College of Textiles Clothing
- Qingdao University
- Qingdao 266071
| | - Zhenggang Lan
- SCNU Environmental Research Institute
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Theoretical Chemistry of Environment
- South China Normal University
- Guangzhou 510006
- China
| |
Collapse
|