You have a model of the opamp, including the gain_bandwidth plot that has high gain at DC, and then perhaps at 10Hz? begins to roll off, with a 90 degree phaseshift.
Now insert a Rout, perhaps 100 ohms or 20 ohms, in series with the voltage_controlled voltage source that is the opamp high_gain model. Ohhh I recall a opamp with 10,000Hz UGBW and the resistor Rout of 80,000 ohms, using only 1microAmp.
With the Rout, you now have an adequate model to derive the Zout with the phase shifts and frequency variation.
Now --- alter the model to become a CLOSED LOOP, gain of +1.
Now --- drive the Closed Loop model with a Current Source.
Now --- work thru the equations, and find Vout/Iout as the frequency varies. And you must bring along all the phase information.
==========================
Now for the surprise.
The falling gain of the opamp, and the 90 degree phase shift, and the Rout lumped component in the model --------- result in an INDUCTOR result.
================================
Let us use the excellent answer (almost a 1_line derivation) of sarthak. And make the open_loop output impedance be just a resistor Ro. Now we can see:
- Zo = Ro / (1 + A * B ) where B is ratio of 1, and A has -90 degree phase shift and A becomes very small as frequency approaches UGBW.
which becomes INDUCTIVE.
As frequency increases, the Zo increases just like an inductive reactance increases with frequency, while keeping +90 phaseshift.
======================================
For fun, evaluate the maths, or run a .AC simulation, with a external Capacitor shunting to Ground.
Notice for frequencies well below the opamp's UGBW, the Zo is very close to a pure inductor, thus there is NO DAMPENING for the ringing generated by the resonance of an external capacitor.
Its is YOUR RESPONSIBLITY to pick a DISCRETE RESISTOR, external to the silicon opamp, that will dampen.
I suggest you consider Rdampen = sqrt(L / C).
Can you now compute L?