Scaling of non-Markovian Monte Carlo wave-function methods

We demonstrate a scaling method for non-Markovian Monte Carlo wave-function simulations used to study open quantum systems weakly coupled to their environments. We derive a scaling equation, from which the result for the expectation values of arbitrary operators of interest can be calculated, all the quantities in the equation being easily obtainable from the scaled Monte Carlo simulations. In the optimal case, the scaling method can be used, within the weak coupling approximation, to reduce the size of the generated Monte Carlo ensemble by several orders of magnitude. Thus, the developed method allows faster simulations and makes it possible to solve the dynamics of the certain class of non-Markovian systems whose simulation would be otherwise too tedious because of the requirement for large computational resources.