Hierarchical equations of motion in the Libra software package

Abstract

We report the implementation of a hierarchical equations of motion (HEOM) module within the open-source Libra software. It includes the standard and scaled HEOM algorithms for computing the dynamics of open quantum systems interacting with a harmonic bath. The module allows the computing of the evolution of the reduced density matrix, as well as spectral lineshapes. The truncation, filtering, and "update list" schemes, as well as OpenMP parallelization, allow for further computational saving. The package is written in a mix of C++ and Python languages, delivering the best compromise between user friendliness and efficiency. The Python layer of the package takes advantage of standard Python libraries, such as h5py, which allows efficient storage and retrieval of the generated results. The package can be seamlessly used within Jupyter notebooks; its careful design shall provide the maximal convenience and intuitiveness to its users.

Improving Lattice Rigidity and Charge Carrier Lifetime by Engineering Spacer Cation of Ruddlesden–Popper Perovskites: A Time-Domain Ab Initio Study

“…The last snapshot was used to generate 6 ps adiabatic MD trajectories for all systems in the microcanonical ensemble with the atomic time step of 1 fs at the Γ-point. Finally, the NA-MD simulations were ran using the last 3 ps geometries as the initial conditions for study of nonradiative charge recombination using the Libra package. − …”

Pedagogical Overview of the Fewest Switches Surface Hopping Method

“…The couplings g i are described using spectral density function, chosen here to be Drude spectral density: J ( ω ) = π 2 ∑ i g i 2 m ω i δ ( ω − ω i ) = 2 λ γ ω γ 2 + ω 2 . The dynamics of this particular Hamiltonian can be computed numerically exactly using a method called the hierarchy equation of motion (HEOM). − …”

“…The nonexponential dynamics and incorrect time scale of decay without decoherence remain even when a thermal bath is included. If a thermal bath is added to the Hamiltonian of eq to simulate the condensed phase, under certain regimes, exact quantum dynamics can be done (using the hierarchical equation of motion − in a strong system bath coupling regime), and Marcus theory can be used under certain regimes (weak diabatic coupling).…”

“…The dynamics of this particular Hamiltonian can be computed numerically exactly using a method called the hierarchy equation of motion (HEOM). 26 − 29 …”

Coherence and Efficient Energy Transfer in Molecular Wires: Insights from Surface Hopping Simulations

“…This method and the algorithm are well described in Refs. [52][53][54][55][56][57]. Our open-source codes to do these calculations can be found in Ref.…”