Optimally chosen Nakajima-Zwanzig master equation for mean field approximation

► Inhomogeneous non-Markovian master equation. ► Complete positivity preserved and correctly equilibrates. ► Calculations for NV center in diamond in good agreement with exact results.