Quantum master equations for non-linear optical response of molecular systems

► We consider non-linear response functions for a simple two-level molecular system. ► It is shown that it cannot be constructed using a single reduced density matrix master equation. ► A set of master equations is derived which together reproduce the two-level response function. ► Parametric projection operators are introduced to derive proper equations in general cases.