IIP2 using Monte Carlo

For a zero-IF double-balanced mixer, is transient analysis with Monte Carlo the only valid way to simulate for IIP2? Can QPSS/Monte Carlo be used in this situation? IF the IIP2 is around 50dBm, what's the appropriate tolarance control(reltol, vabstol, iabstol)? Thank you.