I observe what appears to be an artifact in the 2d contour plot (roughly at centre of the plot.)
The 3d plot does not show the weird crinkle that the 2d plot shows. At first I thought a floating point precision issue, but if that were the case, it would show in the 3d plot as well?
The glitch is due to numerical instability around a stationary point of the function. In Mathematica it can be eliminated, for example, by increasing the PlotPoints: