Multivariable calculus

Critical Points of the Gaussian Bump e^(-x^2-y^2)

The Gaussian bump f(x,y) = e^(-(x^2+y^2)) has exactly one critical point, the origin. There f_xx = f_yy = -2 and f_xy = 0, so D = 4 and the point is a local maximum with value 1. It is the global maximum too, since the exponent -(x^2+y^2) is largest at the origin and e^t is increasing.

f(x,y)=e(x2+y2)f(x,y) = e^{-(x^{2}+y^{2})}
  • (0,0)\left(0, 0\right) local maximumdiscriminant D = 4

The gradient points straight at the origin

Both partials come from the chain rule on the exponent u=(x2+y2)u = -(x^{2}+y^{2}), whose derivatives are 2x-2x and 2y-2y.

fx=2xe(x2+y2),fy=2ye(x2+y2)\frac{\partial f}{\partial x} = -2x\,e^{-(x^{2}+y^{2})}, \qquad \frac{\partial f}{\partial y} = -2y\,e^{-(x^{2}+y^{2})}

Written as a vector, the gradient is 2e(x2+y2)(x,y)-2e^{-(x^{2}+y^{2})}\,(x, y), a negative multiple of the position vector, so at every point other than the origin it aims straight back at the origin. That is the picture of a surface which climbs toward the centre from every direction.

Since the exponential never vanishes, the gradient is zero only when (x,y)=(0,0)(x, y) = (0,0). One critical point, and the list is complete.

For the second partials, differentiate the products above. Each one needs the product rule, since the polynomial factor and the exponential both depend on the variable.

fxx=(4x22)e(x2+y2),fyy=(4y22)e(x2+y2),fxy=4xye(x2+y2)f_{xx} = (4x^{2}-2)e^{-(x^{2}+y^{2})}, \qquad f_{yy} = (4y^{2}-2)e^{-(x^{2}+y^{2})}, \qquad f_{xy} = 4xy\,e^{-(x^{2}+y^{2})}

At the origin the exponential is 1, so fxx=fyy=2f_{xx} = f_{yy} = -2 and fxy=0f_{xy} = 0. Then D=4>0D = 4 > 0 with a negative diagonal entry, a local maximum of value f(0,0)=1f(0,0) = 1.

Two independent confirmations

The discriminant is one route to the answer. Here is a second one that needs no derivatives at all. The exponent (x2+y2)-(x^{2}+y^{2}) is at most 0, with equality only at the origin, and tett \mapsto e^{t} is strictly increasing. So fe0=1f \le e^{0} = 1 everywhere, with equality only at the origin.

(x2+y2)0    e(x2+y2)1,equality only at (0,0)-(x^{2}+y^{2}) \le 0 \implies e^{-(x^{2}+y^{2})} \le 1, \qquad \text{equality only at } (0,0)

That is stronger than what the second derivative test gives. The test says local maximum; this argument says global maximum, and says the maximiser is unique. Whenever a function is a monotone transformation of something obviously extremal, that comparison is worth doing, because it settles global questions the Hessian cannot see.

A third confirmation comes from the quadratic approximation near the origin. Expanding the exponential gives f1(x2+y2)f \approx 1 - (x^{2}+y^{2}), a downward paraboloid whose Hessian is 2-2 on the diagonal and 0 off it. That matches the second partials computed above, which is the check worth running whenever a Hessian looks surprising.

The mistake students make

This surface is constant along every circle centred at the origin, and that leads to a tempting wrong conclusion: since ff does not change as you walk around a circle of radius 1, every point on that circle must be critical. It is not. Being constant along a curve makes the directional derivative along the curve zero, which is one direction out of infinitely many.

The gradient measures all directions at once. At a point on that circle, the gradient is 2e1(x,y)-2e^{-1}(x, y), which has length 2e10.7362e^{-1} \approx 0.736 and points radially inward. The surface is flat around the circle and steep across it, so no point on the circle is critical.

  • A level curve tells you where the directional derivative is zero along that curve.
  • A critical point needs every directional derivative to be zero, which means the whole gradient vector is zero.
  • On this surface only the origin qualifies, because the gradient is a nonzero multiple of (x,y)(x,y) everywhere else.

The other slip is a chain rule error: writing fx=x2e(x2+y2)f_x = -x^{2}e^{-(x^{2}+y^{2})} instead of 2xe(x2+y2)-2x\,e^{-(x^{2}+y^{2})}, by multiplying by the exponent rather than by its derivative. Paired with the matching fy=y2e(x2+y2)f_y = -y^{2}e^{-(x^{2}+y^{2})}, that system still vanishes only at the origin, so the critical point survives and nothing looks wrong yet. The Hessian is where it shows. Differentiating x2e(x2+y2)-x^{2}e^{-(x^{2}+y^{2})} gives (2x32x)e(x2+y2)\left(2x^{3}-2x\right)e^{-(x^{2}+y^{2})} in xx and 2x2ye(x2+y2)2x^{2}y\,e^{-(x^{2}+y^{2})} in yy, and both are zero at the origin, as is the matching fyyf_{yy}. Every entry vanishes, so D=0D = 0 and the test reports degenerate on a surface with an obvious peak. Differentiate the exponent, then multiply, and check the power of xx that comes out.

Frequently asked questions

Why is D exactly 4 with no factor of e?

Because the exponential is evaluated at the critical point itself, and e0=1e^{0} = 1 there. The general second partials are (4x22)e(x2+y2)(4x^{2}-2)e^{-(x^{2}+y^{2})} and so on, but at the origin the exponential factor contributes a clean 1 and only the polynomial part survives. On surfaces whose critical points sit away from the origin, that factor does show up in the discriminant.

How would this change for e^(-(x^2 + 4y^2))?

The single critical point stays at the origin, since the gradient is 2eu(x,4y)-2e^{u}(x, 4y), which still vanishes only there. The second partials become fxx=2f_{xx} = -2 and fyy=8f_{yy} = -8 with fxy=0f_{xy} = 0, so D=16D = 16 rather than 4. It is still a local maximum of value 1, just a narrower bump in the yy direction, and a larger discriminant is the algebraic signature of that steeper curvature.