Abstract:The "hierarchical equations of motion" (HEOM) method is a powerful exact numerical approach to solve the dynamics and find the steady-state of a quantum system coupled to a non-Markovian and nonperturbative environment. Originally developed in the context of physical chemistry, it has also been extended and applied to problems in solid-state physics, optics, single-molecule electronics, and biological physics. Here we present a numerical library in Python, integrated with the powerful QuTiP platform, which implements the HEOM for both bosonic and fermionic environments. We demonstrate its utility with a series of examples consisting of benchmarks against important known results and examples demonstrating insights gained with this library for this article. For the bosonic case, our results include demonstrations of how to fit arbitrary spectral densities with different approaches, and a study of the dynamics of energy transfer in the Fenna-Matthews-Olson photosynthetic complex. For the latter, we both clarify how a suitable non-Markovian environment can protect against pure dephasing, and model recent experimental results demonstrating the suppression of electronic coherence. Importantly, we show that by combining the HEOM method with the reaction coordinate method we can observe nontrivial system-environment entanglement on timescales substantially longer than electronic coherence alone. We also demonstrate results showing how the HEOM can be used to benchmark different strategies for dynamical decoupling of a system from its environment, and show that the Uhrig pulse-spacing scheme is less optimal than equally spaced pulses when the environment's spectral density is very broad. For the fermionic case, we present an integrable single-impurity example, used as a benchmark of the code, and a more complex example of an impurity strongly coupled to a single vibronic mode, with applications to single-molecule electronics. https://doi.org/10.1103/PhysRevResearch.5.013181 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Published by the American Physical Society

SmoQyDQMC.jl: A flexible implementation of determinant quantum Monte Carlo for Hubbard and electron-phonon interactions

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

QMCPACK : An open source ab initio Quantum Monte Carlo package for the electronic structure of atoms, molecules, and solids

QuantumDynamics.jl: A modular approach to simulations of dynamics of open quantum systems

PyQMC: An all-Python real-space quantum Monte Carlo module in PySCF

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

A Hybrid Quantum-Classical Method for Electron-Phonon Systems

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

DeepQMC: an open-source software suite for variational optimization of deep-learning molecular wave functions

MLQD: A package for machine learning-based quantum dissipative dynamics

HierarchicalEOM.jl: An efficient Julia framework for hierarchical equations of motion in open quantum systems

Quantum Package 2.0: An Open-Source Determinant-Driven Suite of Programs

qlbm -- A Quantum Lattice Boltzmann Software Framework

SmoQyDEAC.jl: A differential evolution package for the analytic continuation of imaginary time correlation functions

JuliaQCD: Portable lattice QCD package in Julia language

A coupled cluster framework for electrons and phonons

QUBO.jl: A Julia Ecosystem for Quadratic Unconstrained Binary Optimization

A comparative determinant quantum Monte Carlo study of the acoustic and optical variants of the Su-Schrieffer-Heeger model

QuTiP-BoFiN: A bosonic and fermionic numerical hierarchical-equations-of-motion library with applications in light-harvesting, quantum control, and single-molecule electronics

Pyqcm: An open-source Python library for quantum cluster methods

Pymablock: an algorithm and a package for quasi-degenerate perturbation theory