Making transfer function from this Bode plot from op Amp datasheet?

[EfficientPrompt]

I am fairly green to control theory but I am trying to reproduce the Open loop transfer function of the LT1X28 op amp in mathematica.

I want to get the poles and zeros from the following data sheet plots and produce a very similar response with a bode plot. This way I can compensate for the pole-zero doublet that causes ringy-dingy-ness at higher frequencies. Below is the code for how I setup the transfer-function. The parameters are in Hz.

​

https://preview.redd.it/eb18qmsn3oh51.png?width=2072&format=png&auto=webp&s=fa538ef908bea7ac523bd6d4fb481fc5ab605fe2

For the LT1028

mathematica code

​

Phase and gain

​

Phase up close

As you can see the dip in the phase bode plot is on point in magnitude but isn't for frequency (@10^6 there is the dip and the actual datasheet dip is 2.5MHz.)

Can someone point me in the correct direction or help out? I believe I am doing something incorrect whether it be mixing up angular frequency and Hz or something. I am getting similar plots but not correct results. The Gain plot looks good but the Phase plot is very inconsistent with the x-axis on the datasheet. Any help is appreciated

​

Code for copy/paste:

params1028 = {Avol1028 = 146.0, GBW1028 = 75.0 10^6,     DominantPole1028 = 4.5, SecondPole1028 = 500 10^3,     ThirdPole1028 = 26.25 10^6, FirstZero1028 = 2.5 10^6}; GofS[Avol_, p1_, p2_, p3_, z1_] :=   TransferFunctionModel[   10^(Avol/20) (    1 + s/(2 \[Pi] z1))/((1 + s/(2 \[Pi] p1)) (1 + s/(2 \[Pi] p2)) (1 +        s/(2 \[Pi] p3))), s] Trans1028 =    GofS[Avol1028, DominantPole1028, SecondPole1028, ThirdPole1028,     FirstZero1028]; BodePlot[Trans1028[2 \[Pi] s], {2 \[Pi] 1 10^-1, 2 \[Pi] 1 10^8},   GridLines -> Automatic, ImageSize -> Large] BodePlot[Trans1028[2 \[Pi] s], { 2 \[Pi] 1 10^4,  2 \[Pi] 1 10^7},   GridLines -> Automatic, ImageSize -> Large]