Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.
Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.
Abstract The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond‐dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab‐initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many‐body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree‐Fock self‐consistent field (HF‐SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation‐based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both and symmetries; (2) an efficient second‐order DMRG‐self‐consistent field (SCF) implementation; (3) an externally contracted multi‐reference CI (MRCI) and Epstein‐Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.
Accurate evaluation of electron correlations is essential for the reliable quantitative description of electronic structures in strongly correlated systems, including bond-dissociating molecules, polyradicals, large conjugated molecules, and transition metal complexes. To provide a user-friendly tool for studying such challenging systems, our team developed Kylin 1.0 [J. Comput. Chem. 44, 1316 (2023)], an ab initio quantum chemistry program designed for efficient density matrix renormalization group (DMRG) and post-DMRG methods, enabling high-accuracy calculations with large active spaces. We have now further advanced the software with the release of Kylin 1.3, featuring optimized DMRG algorithms and an improved tensor contraction scheme in the diagonalization step. Benchmark calculations on the Mn4CaO5 cluster demonstrate a remarkable speed-up of up to 16 fater than Kylin 1.0. Moreover, a more user-friendly and efficient algorithm [J. Chem. Theory Comput. 17, 3414 (2021)] for sampling configurations from DMRG wavefunction is implemented as well. Additionally, we have also implemented a spin-adapted version of the externally contracted multi-reference configuration interaction (EC-MRCI) method [J. Phys. Chem. A 128, 958 (2024)], further enhancing the program’s efficiency and accuracy for electron correlation calculations.
The accurate theoretical interpretation of ultrafast time-resolved spectroscopy experiments relies on full quantum dynamics simulations for the investigated system, which is nevertheless computationally prohibitive for realistic molecular systems with a large number of electronic and/or vibrational degrees of freedom. In this work, we propose a unitary transformation approach for realistic vibronic Hamiltonians, which can be coped with using the adaptive time-dependent density matrix renormalization group (t-DMRG) method to efficiently evolve the nonadiabatic dynamics of a large molecular system. We demonstrate the accuracy and efficiency of this approach with an example of simulating the exciton dissociation process within an oligothiophene/fullerene heterojunction, indicating that t-DMRG can be a promising method for full quantum dynamics simulation in large chemical systems. Moreover, it is also shown that the proper vibronic features in the ultrafast electronic process can be obtained by simulating the two-dimensional (2D) electronic spectrum by virtue of the high computational efficiency of the t-DMRG method.
Electronic and/or vibronic coherence has been found by recent ultrafast spectroscopy experiments in many chemical, biological, and material systems. This indicates that there are strong and complicated interactions between electronic states and vibration modes in realistic chemical systems. Therefore, simulations of quantum dynamics with a large number of electronic and vibrational degrees of freedom are highly desirable. Due to the efficient compression and localized representation of quantum states in the matrix-product state (MPS) formulation, time-evolution methods based on the MPS framework, which we summarily refer to as tDMRG (time-dependent density-matrix renormalization group) methods, are considered to be promising candidates to study the quantum dynamics of realistic chemical systems. In this work, we benchmark the performances of four different tDMRG methods, including global Taylor, global Krylov, and local one-site and two-site time-dependent variational principles (1TDVP and 2TDVP), with a comparison to multiconfiguration time-dependent Hartree and experimental results. Two typical chemical systems of internal conversion and singlet fission are investigated: one containing strong and high-order local and nonlocal electron-vibration couplings and the other exhibiting a continuous phonon bath. The comparison shows that the tDMRG methods (particularly, the 2TDVP method) can describe the full quantum dynamics in large chemical systems accurately and efficiently. Several key parameters in the tDMRG calculation including the truncation error threshold, time interval, and ordering of local sites were also investigated to strike the balance between efficiency and accuracy of results.
Modern computational techniques used to simulate quantum chemistry are on the boundary of tractability due to the exponential growth of the molecular wavefunction, requiring a careful balance between molecule size and simulation accuracy. In recent years, quantum computing has risen in popularity as a potential alternative to conventional (classical) techniques; however, most methods rely on access to "digital" quantum computers composed of qubits and quantum gates, which at present are severely limited by noise. We have developed a real-time, analog approach to simulate vibronic chemical dynamics with existing quantum technology. Our approach uses an intuitive mapping of molecular electronic and vibrational degrees of freedom onto quantum resonators and qudit (d-level system) states, with controllable couplings between degrees of freedom. The measurement output can be mapped onto different time-dependent observables, including the time-domain simulation of vibronic spectra. Our approach can also incorporate controlled sources of noise to simulate system-bath interactions and dissipative dynamics at a minimal cost. We present experimental results using a trapped-ion device, thus showing the potential for near-term simulation of chemical dynamics in complex environments beyond the abilities of classical computers.
High Resolution Image Download MS PowerPoint Slide Theoretical predictions of photochemical processes are essential for interpreting and understanding spectral features. Reliable quantum dynamics calculations of vibronic systems require precise modeling of anharmonic effects in the potential energy surfaces and off-diagonal nonadiabatic coupling terms. In this work, we present the n -mode quantization of all vibronic Hamiltonian terms comprised of general high-dimensional model representations. We expand the existing vibrational DMRG formalism by applying the n -mode quantization to all potential energy surfaces entering the Hamiltonian as well as to all off-diagonal coupling terms. This is accompanied by the introduction of a novel matrix product state architecture employing tailored local site operators, which allow an effective encoding of the vibronic wave function in a tensor-train format. This results in a second-quantized framework for accurate vibronic calculations employing the density matrix renormalization group algorithm. We demonstrate the accuracy and reliability of this approach by calculating the excited-state quantum dynamics of maleimide. We analyze convergence and the choice of parameters of the underlying time-dependent density matrix renormalization group algorithm for the n -mode vibronic Hamiltonian, demonstrating that it enables accurate calculations of complex photochemical dynamics.