Multivariable calculus

Second partials of ln(e^x + e^y)

For f(x,y) = ln(e^x + e^y) write S = e^x + e^y. Then f_xx = f_yy = e^x e^y / S^2 and the mixed partial f_xy = f_yx = -e^x e^y / S^2. Every entry of the Hessian has the same magnitude, so the discriminant f_xx f_yy - f_xy^2 is exactly zero at every point and the function is convex but flat along one direction.

f(x,y)=ln⁡(ex+ey)f(x,y) = \ln\bigl(e^{x} + e^{y}\bigr)
∂f∂x=exex+ey\frac{\partial f}{\partial x} = \frac{e^{x}}{e^{x} + e^{y}}
∂f∂y=eyex+ey\frac{\partial f}{\partial y} = \frac{e^{y}}{e^{x} + e^{y}}

Second and mixed partials

∂2f∂x2=exey(ex+ey)2\frac{\partial^2 f}{\partial x^2} = \frac{e^{x}e^{y}}{\left(e^{x} + e^{y}\right)^{2}}
∂2f∂y ∂x=−exey(ex+ey)2\frac{\partial^2 f}{\partial y\,\partial x} = \frac{-e^{x}e^{y}}{\left(e^{x} + e^{y}\right)^{2}}
∂2f∂y2=exey(ex+ey)2\frac{\partial^2 f}{\partial y^2} = \frac{e^{x}e^{y}}{\left(e^{x} + e^{y}\right)^{2}}

Differentiating the log-sum-exp

Write S=ex+eyS = e^{x} + e^{y}, which is positive for all real inputs, so the logarithm is defined everywhere. The chain rule gives inner over outer, and only one term of SS survives each partial derivative.

∂f∂x=exS,∂f∂y=eyS\frac{\partial f}{\partial x} = \frac{e^{x}}{S}, \qquad \frac{\partial f}{\partial y} = \frac{e^{y}}{S}

Those two first partials are the softmax weights: each is between 00 and 11, and they add to exactly 11. Differentiate ex/Se^{x}/S in xx with the quotient rule, remembering that SS contains exe^{x}.

∂2f∂x2=exS−exexS2=ex(S−ex)S2=exeyS2\frac{\partial^2 f}{\partial x^2} = \frac{e^{x}S - e^{x}e^{x}}{S^{2}} = \frac{e^{x}\left(S - e^{x}\right)}{S^{2}} = \frac{e^{x}e^{y}}{S^{2}}

For the mixed partial the numerator exe^{x} is constant in yy, so only the denominator contributes, and it contributes with a minus sign.

∂2f∂y ∂x=−ex⋅eyS2\frac{\partial^2 f}{\partial y\,\partial x} = -\frac{e^{x}\cdot e^{y}}{S^{2}}

The mistake: treating the denominator as a constant

The tempting wrong move is to say that since fx=ex/Sf_x = e^{x}/S looks like an exponential, differentiating again just returns ex/Se^{x}/S. That ignores the exe^{x} hiding inside SS. The wrong answer climbs toward 11 as xx pulls ahead of yy, whereas the true curvature decays to 00 out there and never exceeds 1/41/4 anywhere.

A useful check is the identity fxx=fx(1−fx)f_{xx} = f_x(1 - f_x), which follows because fx=ex/Sf_x = e^x/S and 1−fx=ey/S1 - f_x = e^y/S. Since fxf_x lies strictly between 00 and 11, the curvature is always positive and peaks at 1/41/4 on the line x=yx = y.

The other slip is the sign of the mixed partial. It has to be negative: raising yy increases SS, which shrinks the weight ex/Se^{x}/S, so the slope in the xx direction goes down.

A Hessian with determinant zero

All three second partials have the same magnitude exey/S2e^{x}e^{y}/S^{2}, with a minus sign only on the mixed one, so the discriminant collapses.

D=fxxfyy−fxy 2=(exeyS2)2−(exeyS2)2=0D = f_{xx}f_{yy} - f_{xy}^{\,2} = \left(\frac{e^{x}e^{y}}{S^{2}}\right)^{2} - \left(\frac{e^{x}e^{y}}{S^{2}}\right)^{2} = 0

The reason is a symmetry you can see without any calculus. Adding the same constant cc to both inputs multiplies SS by ece^{c}, so f(x+c,y+c)=f(x,y)+cf(x+c, y+c) = f(x,y) + c. The function is exactly linear along the diagonal direction, which means zero curvature there and forces a zero eigenvalue in the Hessian.

The other eigenvalue is 2exey/S22e^{x}e^{y}/S^{2}, which is positive, so the surface curves upward across the diagonal. A Hessian that is positive semidefinite everywhere makes this function convex, which is why log-sum-exp is a standard smooth stand-in for the maximum of xx and yy.

Frequently asked questions

Why is the discriminant of ln(e^x + e^y) always zero?

Because f(x+c,y+c)=f(x,y)+cf(x+c,y+c) = f(x,y) + c, the function is linear along the diagonal, so the Hessian has a zero eigenvalue in that direction. Concretely all three second partials share the magnitude exey/S2e^{x}e^{y}/S^{2}, and fxxfyyf_{xx}f_{yy} cancels fxy2f_{xy}^2.

How large can f_xx get?

At most 1/41/4. Since fxx=fx(1−fx)f_{xx} = f_x(1-f_x) with fx=ex/(ex+ey)f_x = e^{x}/(e^{x}+e^{y}) strictly between 00 and 11, the product is largest when fx=1/2f_x = 1/2, which happens exactly on the line x=yx = y.