Abstract
Open quantum systems are ubiquitous in nature and central to quantum technologies. A common description of their dynamics is given by the celebrated Lindblad master equation, which can be generalized to the non-Markovian scenario. In this work, we introduce the Non-Markovian Ensemble Propagation (NMEP) method, which extends the Monte Carlo Wave-Function (MCWF) method to the non-Markovian case in a simple and general manner. We demonstrate its accuracy and effectiveness in a selection of examples, and compare the results with either analytic expressions or direct numerical integration of the master equation.