Non-Markovian modeling of Fermi-Bose systems coupled to one or several Fermi-Bose thermal baths

A method is proposed to describe Fermi or Bose systems coupled to one or several heat baths composed of fermions and/or bosons. The method, called the coupled equations of motion method, properly includes non-Markovian effects. The approach is exact in the full-coupling approximation when only bosonic particles are present in the system and baths. The approach provides an approximate treatment when fermions are present either in the system and/or in one or several environments. Our approach has the advantage of properly respecting the Pauli exclusion principle for fermions during the evolution. We illustrate the approach for the single fermionic or bosonic oscillator coupled to one or two heat baths assuming different types of quantum statistics (fermion or boson) for them. The cases of a Fermi system coupled to fermion or boson heat baths or a mixture of both are analyzed in detail. With the future goal of treating Fermi systems formed of an increasing number of two-level systems (qubits), we discuss possible simplifications that could be made in the equations of motion and their limits of validity in terms of the system-bath coupling or of the initial heat bath temperatures.