We derive the general form of a master equation describing the reduced time evolution of a sequence of subsystems “propagating” in an environment which can be described as a sequence of subenvironments. The interaction between subsystems and subenvironments is described in terms of a collision model, with the irreversible dynamics of the subenvironments between collisions explicitly taken into account. In the weak coupling regime, we show that the collisional model produces a correlated Markovian evolution for the joint density matrix of the multipartite system. The associated Lindblad superoperator contains pairwise terms describing cross correlation between the different subsystems. Such a model can describe a broad range of physical situations, ranging from quantum channels with memory to photon propagation in concatenated quantum optical systems.

• Received 17 May 2011

DOI:https://doi.org/10.1103/PhysRevLett.108.040401