advanced32 min studyLesson 347 in Mathematics

Jacobi fields and conjugate points

Linearize the stationary equation and use endpoint zeros of its variational equation to detect a loss of second-variation positivity.

On this page

What you’ll understand

  • Derive the scalar Jacobi equation from the second variation.
  • Define conjugacy relative to a fixed initial endpoint.
  • Analyze interval-length thresholds using an exactly solvable quadratic model.

Two nearby stationary curves can start at the same point with slightly different slopes and later refocus. Their infinitesimal separation satisfies a linear equation called the Jacobi equation. A second zero of that separation is a conjugate point. The idea connects the geometry of neighboring solutions with the quadratic form that tests local stability. In a scalar quadratic example, the first conjugate point marks the exact interval length at which a previously positive energy acquires a zero direction; extending the interval further creates a negative direction.

Linearize the stationary equation

Take a scalar smooth stationary curve for a one-dimensional fixed-endpoint functional. Assume enough regularity, for example a three-times continuously differentiable integrand and a smooth curve, so the coefficients below are defined and the linearization is justified. In the second variation, integrate the mixed term by parts using zero endpoint perturbations. The quadratic form becomes the integral of P(h′)2+Qh2P(h')^2+Qh^2, where P=FvvP=F_{vv} and Q=Fuu−(Fuv)′Q=F_{uu}-(F_{uv})', all evaluated along the stationary curve.

Applying the Euler–Lagrange procedure to this quadratic form gives the linear Jacobi equation. The same equation arises by differentiating a smoothly parameterized family of stationary curves with respect to its family parameter. Thus its solution is an infinitesimal displacement between neighboring extremals, not necessarily a finite difference that satisfies the original nonlinear equation exactly. The word field is also used in higher-dimensional geometric settings; here it is simply a scalar function along the chosen curve.

δ2J[u](h,h)=∫ab[P(h′)2+Qh2]dx,−(Pj′)′+Qj=0,P=Fvv,Q=Fuu−ddxFuv\delta^2J[u](h,h)=\int_a^b[P(h')^2+Qh^2]dx,\qquad -(Pj')'+Qj=0,\quad P=F_{vv},\quad Q=F_{uu}-\frac d{dx}F_{uv}
The mixed-term boundary contribution vanishes for fixed endpoints. The Jacobi equation is the stationary equation of the quadratic second variation.↔ Scroll wide equations horizontally.

Conjugate endpoints identify zero directions

Assume PP is continuously differentiable and strictly positive on the interval, and QQ is continuous. Solve the Jacobi equation with j(a)=0j(a)=0 and j′(a)=1j'(a)=1. Any later point where this nonzero solution vanishes is conjugate to the initial endpoint. Changing its nonzero initial derivative merely rescales the solution, so the zero locations do not depend on that normalization. At a conjugate terminal endpoint, the Jacobi field itself is an admissible zero direction of the quadratic form, as follows by multiplying its equation by the field and integrating.

Before the first conjugate point, a factorization explains positivity. Where the Jacobi solution does not vanish, write the quadratic form as a squared expression involving h′−(j′/j)hh'-(j'/j)h plus an endpoint term. For smooth homogeneous perturbations the apparent singular term at the initial zero has a vanishing boundary limit. If the terminal value of the Jacobi solution is nonzero, equality in the square forces a proportionality incompatible with the terminal perturbation being zero. This is a positivity argument for the second variation; nonlinear minimality still needs the appropriate uniform estimates and remainder control.

P(h′)2+Qh2=P(h′−j′jh)2+ddx(Pj′jh2)(j≠0)P(h')^2+Qh^2=P\left(h'-\frac{j'}j h\right)^2+\frac d{dx}\left(P\frac{j'}j h^2\right)\quad(j\ne0)
This identity uses the Jacobi equation. Boundary limits and endpoint zeros must be handled rather than dividing by a Jacobi field at its zeros.↔ Scroll wide equations horizontally.

Length matters, even for a fixed local coefficient

For the quadratic action with density [(u′)2−ω2u2]/2[(u')^2-\omega^2u^2]/2, the zero curve is stationary and the Jacobi equation is j′′+ω2j=0j''+\omega^2j=0. With the stated normalization its solution is sin⁡(ωx)/ω\sin(\omega x)/\omega when the initial endpoint is zero. The first positive conjugate point is π/ω\pi/\omega. The derivative coefficient remains one for every interval length, so a local Legendre condition cannot detect this global-in-the-interval change by itself.

Poincaré's sharp interval inequality makes the threshold exact: the second variation is bounded below by [1−(ωT/π)2]∫(h′)2[1-(\omega T/\pi)^2]\int(h')^2 on an interval of length TT. Below the threshold it is coercive; at the threshold a sine mode has zero energy; above it the first sine mode has negative energy. These are complete statements for this quadratic functional. In a nonlinear problem, a zero second variation alone cannot decide higher-order stability, and a geometric conjugate-point argument should not silently replace a full sufficient-condition theorem.

A stable interval before refocusing

Consider J[u]=12∫01[(u′)2−4u2]dxJ[u]=\tfrac12\int_0^1[(u')^2-4u^2]dx on H01(0,1)H_0^1(0,1). Find the first conjugate point to zero and prove the zero curve is the unique minimizer on this interval.

  1. The Jacobi equation is j′′+4j=0j''+4j=0. With j(0)=0,j′(0)=1j(0)=0,j'(0)=1, its solution is j(x)=12sin⁡(2x)j(x)=\tfrac12\sin(2x), whose first positive zero is π/2\pi/2.

  2. The terminal coordinate one lies before π/2\pi/2. More quantitatively, Poincaré gives ∫01h2≤π−2∫01(h′)2\int_0^1h^2\le\pi^{-2}\int_0^1(h')^2.

  3. Therefore δ2J[0](h,h)≥(1−4/π2)∫(h′)2>0\delta^2J[0](h,h)\ge(1-4/\pi^2)\int(h')^2>0 for every nonzero homogeneous perturbation. The functional is exactly half this quadratic form, so zero is its unique global minimizer.

The result

The first conjugate length is π/2, and the shorter unit interval has a positive stability margin 1 − 4/π² in the derivative norm.

Explore the model

Interactive lab

Interactive lab · oscillator

A zero crossing is not automatically a conjugate point

Reset and set both ζ=0,ω0=2\zeta=0,\omega_0=2. Predict the instrument's displacement cos⁡(2t)\cos(2t) reaches its first zero at t=π/4t=\pi/4, while the Jacobi solution sin⁡(2t)/2\sin(2t)/2 first returns to zero at t=π/2t=\pi/2. Use Play or Step and compare displacement and velocity: the instrument has velocity −2sin⁡(2t)-2\sin(2t), so the Jacobi field equals −v/4-v/4. Its initial data differ from the displacement's. This is the same linear equation with different initial conditions, not a boundary-value optimizer.

A zero crossing is not automatically a conjugate pointUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forConjugacy depends on the Jacobi equation and its specified endpoint-zero initial data.

Check your understanding

For j″ + 4j = 0 with j(0) = 0 and j′(0) = 1, where is the first positive conjugate point?

Try it yourself

For the same quadratic density, compare terminal lengths T=π/2T=\pi/2 and T=πT=\pi. Exhibit a zero direction at the first length and a negative direction at the second.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

At the first length use the Jacobi field. At the longer length use the first Dirichlet sine mode adapted to that interval.

Compare with the solution

At T=π/2T=\pi/2, h=sin⁡(2x)h=\sin(2x) has zero endpoints and ∫[(h′)2−4h2]=0\int[(h')^2-4h^2]=0. Poincaré shows the whole quadratic form is nonnegative there, so zero is a nonunique global minimum with all multiples of this sine also minimizing. At T=πT=\pi, take h=sin⁡xh=\sin x: the quadratic form is −3π/2<0-3\pi/2<0, so zero is not a local minimum.

Take it with you

  • The Jacobi equation linearizes the stationary equation along a candidate.
  • A conjugate point is a second zero of the endpoint-normalized Jacobi field.
  • For the quadratic model, the first conjugate length exactly separates positive, zero, and negative directions.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.