🌱 Tony's Memex Tony's Intertwingled Memex
Digital Garden | OER Curriculum
Home / Digital Garden / Second-Order Partial Derivatives

Second-Order Partial Derivatives

Core Definition

For a function of two variables z=f(x,y)z = f(x, y), the second-order partial derivatives are obtained by differentiating its first-order partial derivatives (fxf_x and fyf_y). There are four standard second-order partial derivatives:

  1. Pure Second Partial with respect to xx:
    fxx=βˆ‚2fβˆ‚x2=βˆ‚βˆ‚x(βˆ‚fβˆ‚x)f_{xx} = \frac{\partial^2 f}{\partial x^2} = \frac{\partial}{\partial x}\left(\frac{\partial f}{\partial x}\right)
  2. Pure Second Partial with respect to yy:
    fyy=βˆ‚2fβˆ‚y2=βˆ‚βˆ‚y(βˆ‚fβˆ‚y)f_{yy} = \frac{\partial^2 f}{\partial y^2} = \frac{\partial}{\partial y}\left(\frac{\partial f}{\partial y}\right)
  3. Mixed Second Partial (xx first, then yy):
    fxy=βˆ‚2fβˆ‚yβˆ‚x=βˆ‚βˆ‚y(βˆ‚fβˆ‚x)f_{xy} = \frac{\partial^2 f}{\partial y \partial x} = \frac{\partial}{\partial y}\left(\frac{\partial f}{\partial x}\right)
  4. Mixed Second Partial (yy first, then xx):
    fyx=βˆ‚2fβˆ‚xβˆ‚y=βˆ‚βˆ‚x(βˆ‚fβˆ‚y)f_{yx} = \frac{\partial^2 f}{\partial x \partial y} = \frac{\partial}{\partial x}\left(\frac{\partial f}{\partial y}\right)

By Clairaut's Theorem, if fxyf_{xy} and fyxf_{yx} are continuous on an open region, then fxy=fyxf_{xy} = f_{yx}.

The symmetric square matrix containing all second-order partials is the Hessian Matrix:

Hf(x,y)=[fxxfxyfyxfyy]\mathbf{H}_f(x, y) = \begin{bmatrix} f_{xx} & f_{xy} \\ f_{yx} & f_{yy} \end{bmatrix}

Concept Explanation

Geometric Meaning of Pure vs. Mixed Partials

  • Pure Partials (fxx,fyyf_{xx}, f_{yy}): Measure the concavity (rate of change of slope) of the surface z=f(x,y)z = f(x, y) along vertical slices parallel to the coordinate axes:
    • If fxx>0f_{xx} > 0, the surface is concave upward along the xx-direction (curving like a bowl or trough).
    • If fxx<0f_{xx} < 0, the surface is concave downward along the xx-direction (curving like a dome or ridge).
  • Mixed Partials (fxy,fyxf_{xy}, f_{yx}): Measure the twisting or cross-coupling shear of the surface:
    • fxyf_{xy} describes how the slope in the xx-direction tilts as you step in the yy-direction.
    • Non-zero mixed partials indicate that changes in one variable directly influence the system's sensitivity to another variable.

Notation Order Comparison

Notice the difference between subscript (Lagrange) notation and differential operator (Leibniz) notation:

  • Subscript Notation reads Left-to-Right:
    fxy=(fx)yβ€…β€ŠβŸΉβ€…β€ŠdifferentiateΒ withΒ respectΒ toΒ xΒ first,Β thenΒ yf_{xy} = (f_x)_y \implies \text{differentiate with respect to } x \text{ first, then } y
  • Leibniz Notation reads Right-to-Left (from inside out):
    βˆ‚2fβˆ‚yβˆ‚x=βˆ‚βˆ‚y(βˆ‚fβˆ‚x)β€…β€ŠβŸΉβ€…β€ŠdifferentiateΒ withΒ respectΒ toΒ xΒ first,Β thenΒ y\frac{\partial^2 f}{\partial y \partial x} = \frac{\partial}{\partial y}\left(\frac{\partial f}{\partial x}\right) \implies \text{differentiate with respect to } x \text{ first, then } y
   Subscript Order:    f_{ x  y }   ===>  (1) x first,   (2) y second
                          |  |
                          1  2

   Leibniz Order:      rac{\partial^2 f}{\partial y \partial x}  ===>  (1) x first,   (2) y second
                                                     |        |
                                                     2        1

The Hessian Matrix and Curvature

The Hessian matrix Hf\mathbf{H}_f is the multivariable generalization of the single-variable second derivative fβ€²β€²(x)f''(x):

Hf(x,y)=[fxxfxyfyxfyy]\mathbf{H}_f(x, y) = \begin{bmatrix} f_{xx} & f_{xy} \\ f_{yx} & f_{yy} \end{bmatrix}
The Hessian determinant (or discriminant) is:
D=det⁑(Hf)=fxxfyyβˆ’(fxy)2D = \det(\mathbf{H}_f) = f_{xx} f_{yy} - (f_{xy})^2
In multivariable optimization (Second Derivative Test):

  • If D>0D > 0 and fxx>0f_{xx} > 0, the surface has positive Gaussian curvature in all directions, corresponding to a local minimum.
  • If D>0D > 0 and fxx<0f_{xx} < 0, the surface has negative curvature in all directions, corresponding to a local maximum.
  • If D<0D < 0, the surface bends upward in one direction and downward in another, forming a saddle point.
  • If D=0D = 0, the test is inconclusive.

Engineering Applications

  • Structural Engineering: Strain energy functions U(x,y)U(x, y) in steel columns evaluate stability under axial loading xx and thermal expansion yy using Hessian determinants to verify structural resistance against buckling (see Structural Column Strain Energy Stability (Civil & Mechanical Engineering)).
  • Electrical & Computer Engineering: Closed-loop automatic gain control circuits in audio amplifiers use Hessian matrices of the amplification function G(v,r)G(v, r) to locate optimal bias voltages and feedback resistances (see Automated Audio Amplifier Loop Gain (Electrical & Computer Engineering)).
  • Robotics & Kinematics: Joint torque equations utilize second partials to compute centrifugal and Coriolis coupling terms in dynamic multi-link robot arms.
  • Geotechnical Engineering: Subsurface settlement basins model soil compaction where pure partials quantify settlement curvature under foundation footings.

Example Problems

Example 1: Full Second-Order Partial Evaluation and Clairaut Verification

Compute all four second-order partial derivatives for f(x,y)=x3y4βˆ’5x2y+2y2f(x, y) = x^3 y^4 - 5x^2 y + 2y^2 and verify Clairaut's Theorem.

Solution
  1. Compute First-Order Partials:

    fx=βˆ‚βˆ‚x[x3y4βˆ’5x2y+2y2]=3x2y4βˆ’10xyf_x = \frac{\partial}{\partial x}[x^3 y^4 - 5x^2 y + 2y^2] = 3x^2 y^4 - 10xy
    fy=βˆ‚βˆ‚y[x3y4βˆ’5x2y+2y2]=4x3y3βˆ’5x2+4yf_y = \frac{\partial}{\partial y}[x^3 y^4 - 5x^2 y + 2y^2] = 4x^3 y^3 - 5x^2 + 4y

  2. Compute Pure Second Partials:

    fxx=βˆ‚βˆ‚x[3x2y4βˆ’10xy]=6xy4βˆ’10yf_{xx} = \frac{\partial}{\partial x}[3x^2 y^4 - 10xy] = 6x y^4 - 10y
    fyy=βˆ‚βˆ‚y[4x3y3βˆ’5x2+4y]=12x3y2+4f_{yy} = \frac{\partial}{\partial y}[4x^3 y^3 - 5x^2 + 4y] = 12x^3 y^2 + 4

  3. Compute Mixed Partials:

    fxy=βˆ‚βˆ‚y[3x2y4βˆ’10xy]=12x2y3βˆ’10xf_{xy} = \frac{\partial}{\partial y}[3x^2 y^4 - 10xy] = 12x^2 y^3 - 10x
    fyx=βˆ‚βˆ‚x[4x3y3βˆ’5x2+4y]=12x2y3βˆ’10xf_{yx} = \frac{\partial}{\partial x}[4x^3 y^3 - 5x^2 + 4y] = 12x^2 y^3 - 10x

  4. Clairaut Verification: Since fxy=12x2y3βˆ’10x=fyxf_{xy} = 12x^2 y^3 - 10x = f_{yx}, Clairaut's Theorem is verified.

Example 2: Hessian Matrix Assembly and Determinant Evaluation

For the potential energy function V(x,y)=x3+y3βˆ’3xyV(x, y) = x^3 + y^3 - 3xy, construct the Hessian matrix HV(x,y)\mathbf{H}_V(x, y) and evaluate its determinant at the critical points (0,0)(0, 0) and (1,1)(1, 1).

Solution
  1. First-Order Partials:
    Vx=3x2βˆ’3y,Vy=3y2βˆ’3xV_x = 3x^2 - 3y, \qquad V_y = 3y^2 - 3x
  2. Second-Order Partials:
    Vxx=6x,Vyy=6y,Vxy=βˆ’3V_{xx} = 6x, \qquad V_{yy} = 6y, \qquad V_{xy} = -3
  3. Hessian Matrix:
    HV(x,y)=[6xβˆ’3βˆ’36y]\mathbf{H}_V(x, y) = \begin{bmatrix} 6x & -3 \\ -3 & 6y \end{bmatrix}
    D(x,y)=det⁑(HV)=(6x)(6y)βˆ’(βˆ’3)2=36xyβˆ’9D(x, y) = \det(\mathbf{H}_V) = (6x)(6y) - (-3)^2 = 36xy - 9
  4. Evaluate at Points:
    • At (0,0)(0, 0): D(0,0)=36(0)(0)βˆ’9=βˆ’9<0β€…β€ŠβŸΉβ€…β€ŠD(0, 0) = 36(0)(0) - 9 = -9 < 0 \implies Saddle Point.
    • At (1,1)(1, 1): D(1,1)=36(1)(1)βˆ’9=27>0D(1, 1) = 36(1)(1) - 9 = 27 > 0. Since Vxx(1,1)=6>0β€…β€ŠβŸΉβ€…β€ŠV_{xx}(1, 1) = 6 > 0 \implies Local Minimum.

Connections & References