For the case that ωA and ωB , on the other hand, are complex, we allow the replacement of Eq. ( ) and the response function becomes complex, i.e., damped. It follows that the implementation of the damped quadratic response function of Eq. ( ) is a (conceptually) simple generalization of the implementation for real frequencies.