Search

The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. The method is rooted in an exponential expansion of the bath correlation function, which in essence strategically reshapes a continuous environment into a set of effective bath modes that allow for more efficient cutoff at finite temperatures. Based on this understanding, one can map the HEOM method into a Schrödinger-like equation, with a non-Hermitian super-Hamiltonian for an extended wave function being the tensor product of the central system wave function and the Fock state of these effective bath modes. In this work, we explore the possibility of representing the extended wave function as a tree tensor network state (TTNS) and the super-Hamiltonian as a tree tensor network operator of the same structure as the TTNS, as well as the application of a time propagation algorithm using the time-dependent variational principle. Our benchmark calculations based on the spin-boson model with a slow-relaxing bath show that the proposed HEOM+TTNS approach yields consistent results with those of the conventional HEOM method, while the computation is considerably sped up. In addition, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme. [ABSTRACT FROM AUTHOR]

We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation [R. Borrelli, J. Chem. Phys. 150, 234102 (2019)] to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation. [ABSTRACT FROM AUTHOR]

Within the well-established optical response function formalism, a new strategy with the central idea of employing the forward-backward stochastic Schrödinger equations in a segmented way to accurately obtain the two-dimensional (2D) electronic spectrum is presented in this paper. Based on the simple excitonically coupled dimer model system, the validity and efficiency of the proposed schemes are demonstrated in detail, along with the comparison against the deterministic hierarchy equations of motion and perturbative second-order time-convolutionless quantum master equations. In addition, an important insight is provided in this paper that the characteristic frequency of the overdamped environment is an extremely crucial factor to regulate the lifetimes of the oscillating signals in 2D electronic spectra and of quantum coherence features of system dynamics. It is worth noting that the proposed scheme benefiting from its stochastic nature and wavefunction framework and many other advantages of substantially reducing the numerical cost has a great potential to systematically investigate various quantum effects observed in realistic large-scale natural and artificial photosynthetic systems. [ABSTRACT FROM AUTHOR]