The point is that the difference between the input and output is not the same as the gain error. It combines the phase shift and amplitude error.
Did you bother to actually look at the various plots and previous analysis?? The |Vout-Vin| error is the same regardless of op-amp phase!! This is verified in the AC Linear Plots, Time Domain Plots and running RMS computations, the |Vout-Vin| Error is almost independent of the op-amp phase. Must we remind that
Unity Gain means Unity Gain, that's exactly 1 and nothing else!! The Unity Gain Error is the difference between the output and unity with an input of 1, and since the input in AC analysis is Unity it can also be expressed as Unity Gain Error = Vout - 1, or Vout - Vin since Vin is 1 for AC analysis.
With the simple OP-amp and not very low frequencies (e.g. 100 Hz) the difference between the input and output is essentially out of phase from the main signal. Adding such an out of phase voltage to the input gives mainly a phase shift, but little amplitude error. The simulation also shows the small amplitude error, if one directly looks at the amplitude.
So the simple buffer is already much better in the amplitude than just looking at the input to output difference suggests. Most of the input to output difference is from the phase shift.
The simple buffer is
NOT much better in magnitude error and has been shown by simple analysis to suffer from finite op-amp loop gain at frequency. Look at the Error plots in both Linear AC analysis, Time Domain and running RMS computational results, these all agree with the analysis we provided.
The simple answer, if you chose to accept it, is that the op-amp gain magnitude at frequency limits the Unity Gain Follower Precision, regardless of op-amp phase, and Unity Gain is not anywhere near as good as you portray. We've shown this with analysis, and multiple simulations in AC, Time Domain and RMS computations.
Do you agree that the Unity Gain Follower behaves as Gain Transfer Function = A(f)/((A(f)+1), where A(f) is the op-amp gain at frequency f, or do you have another better model? Do you agree that the error in a Unity Gain Follower is anything that deviates away from a gain of exactly 1? How do you explain the simulation results we've provided which all agree with the previous analysis we also provided??
Look at the OP07 Follower, it has a A(f) of 78.6dB at 100Hz and in the Follower shows a Vout-Vin Magnitude Error of -78.36dB, that's 120ppm, which agrees with the Time Domain sinewave error of 1.2mV/10V (Vin is 10V peak, or 120ppm)!!
This is a similar result for the ideal VCVS Follower which has A(f) of 78dB everywhere with no phase shift and a Vout-Vin error of -78dB (125ppm), the Time Domain simulation also shows ~1.2mV/10V, or ~120ppm for the VCVS Follower. Can you explain this result with a op-amp VCVS model that has no inherent phase shift anywhere, yet exhibits almost identical magnitude behavior to the OP07 Follower? If you examine the Time Domain result carefully, you will discover that the VCVS Follower has a voltage minimum of -1.2mV exactly at the Input Sinewave peak 10V (120ppm) which indicates that the loop gain isn't high enough to drive that towards zero. Just for fun we increased the VCVS gain to 120dB and the result was exactly as expected a minimum of 10uV at the input 10V peak (1ppm), the LTspice simulation is cleaner since the VCVS is ideal and the "numerical noise" less.
The Dual Amp approach alleviates the dependance upon op-amp gain at lower frequencies and provides a much closer Unity Gain response. The Dual Amp Follower Transfer Function is (A(f)^2 + A(f))/(A(f)^2 + A(f) + 1) for matched op-amps, and (A(f)B(f) + A(f))/(A(f)B(f) + B(f) + 1) where B(f) is the feedback op-amp gain at frequency. It should be clear to anyone that these produce a much better result than the Standard Unity Gain Follower above, mainly due to the A(f)^2 or A(f)B(f) terms in the numerator and denominator. Again simulation results with Linear AC, Time Domain and RMS computations prove this out!!!
Best