Before I get to know how to decide the frequency dependance of an op amp's gain by simulation program, let me say something I've learned recently about the open-loop gain of an operational amplifier.
I haven't study the internal circuits of an op amp yet, but
if the open-loop gain Aol of an op amp can be described by an expression which is the same as that of a low-pass RC filter:
Aol(s) = A / (1+s/wb) -------- (1)
where 'A' is the op amp's dc gain and 'wb' is its 3-dB frequency.
For frequency w >> wb, the magnitude of Aol can be approximated by
| Aol(jw) | ≈ A*wb / w = wt/w ---- (2)
in which wt=(A*wb) is known as the unity-gain bandwidth of the op amp, and (2) can also be written as
Aol*w ≈ A*wb = wt
or
Aol*f ≈ A*fb = ft =GBW ---- (3)
which says that the product of an op amp's gain and the frequency at which the gain is measured is a constant equal to the op amp's unity-gain bandwidth for frequencies that are much larger than the op amp's 3-dB frequency. The product happens to be the definition of the gain bandwidth product (GBW). Therefore, if an op amp's open-loop gain can be described by (1), then for any frequency that is much larger than the op amp's break frequency, the op amp's GBW is equal to its unity-gain frequency.
As a test, let's see how an inverting opeational amplifier's gain relates to frequency and its GBW.
Assuming a finite open-loop op amp gain Aol, the closed-loop gain of the inverting amplifier is
Vout/Vi = (-R2/R1) / [1 + (1+R2/R1)/Aol] ----- (4)
inserting (1) for Aol into (4)
Vout/Vi = (-R2/R1) / [ 1+ (1+R2/R1)/A + (1+R2/R1)*s/A*wb ]
≈ (-R2/R1) / [ 1 + (1+R2/R1)*s/A*wb ], for A >> (1+R2/R1), A = dc gain of the op amp
Replacing s by ( jw)=(j*2*π*f) and wb by (wt/A)=(2*π*ft/A),
Vout/Vi = (-R2/R1) / [1 + (1+R2/R1)*(jf/ft)
= (-R2/R1) / [1+ (1+R2/R1)*(jf/GBW)
With R1=1k, R2=3k, GBW=1megHz, f=1meg
|Vout/Vi| = 3/√17 = -2.762 dB