A new method of treating open systems is presented. The normal treatment using the generalized master equation with the projection of Argyres and Kelley is meaningful only if the state of the reservoir never deviates appreciably from the reference state which appears in the projection. Otherwise, one must make at least a partial resummation of the perturbative expansion of the kernel of the generalized master equation. The present method avoids the introduction of a projection operator and allows us to overcome such resummation difficulties. It is based on an integrodifferential equation for the statistical operator of the composite system, which naturally provides a hierarchy of equations involving the statistical operator ϱ(t) of the open system and suitable quantities describing higher and higher order bath-system correlations. Treating the deviations of the bath from its initial equilibrium or stationary state as expansion parameters, one gets an approximation scheme, each step of which gives a closed system of equations for ϱ(t) and a suitable set of correlation quantities.
