# Determination of System poles through open/short circuit time constant technique

I was reading about the open and short circuit time constant technique for estimating the dominant and non-dominant poles of the system respectively. For the dominant pole (open circuit time constant) we short all the independent voltage sources and replace all but one of the capacitors with open circuit. We then calculate the time constant of the resultant circuit and add the time constants for all the capacitors and the dominant pole frequency is given by: $$p_{\text{dom}} = 1/\Sigma_i(\tau_i)$$ where $\tau_i$ is the is one of the time constants of the circuit.

Similarly for the short circuit time constant technique to estimate the non-dominant pole we replace each capacitor by a short circuit and the pole is given by:

$$p_{\text{non_dom}} = \Sigma_i(1/\tau_i)$$ (For details: http://paginas.fe.up.pt/~fff/eBook/MDA/Met_con_tem.html)

Consider the following circuit: For this circuit the non dominant poles are calculated in Gray Meyer using the above procedure:

$\tau_1 = C_fR_L$ and $\tau_2 = C_{gs}[R_S||1/g_m||R_L]$ and $p_{\text{non_dom}} = 1/\tau_1 +1/\tau_2$ But suppose I add a load capacitor to the output as shown: Now when I try to calculate the non-dominant pole from the above method, every other capacitor seems to be shorted making the time delay for each capacitor to be zero.
How should we calculate the non-dominant pole of this new circuit?