advanced30 min studyLesson 444 in Mathematics

Symplectomorphisms and generating functions

Verify canonical maps by pullback and derive a local transformation from a generating function with a regular projection.

On this page

What you’ll understand

  • Test a linear map using the symplectic matrix identity.
  • Derive a type-two generating function transformation with correct signs.
  • Identify the local invertibility and topology limits of a chosen generating representation.

A canonical transformation preserves the pairing that turns energies into Hamiltonian motion. It may shear, rotate or stretch phase space, provided conjugate directions compensate appropriately. Generating functions encode certain such maps through a scalar derivative rule, making preservation automatic once the required local equations can be solved. They are coordinate representations with regularity assumptions, not a promise that every symplectomorphism admits one preferred global formula. We will verify both the pullback condition and the solvability condition rather than treating a generating expression as a complete transformation by itself.

Preserve the two-form, not merely a determinant

A diffeomorphism Φ between symplectic manifolds is symplectic when Φ*ω_target=ω_source. In canonical coordinates, a linear transformation S must satisfy SᵀJS=J, with J the matrix representing ω on tangent vectors. Preserving the top wedge gives det S=1 for real symplectic matrices, but the converse fails in phase dimension four and above. For example, scaling both q₁ and p₁ by two and both q₂ and p₂ by one half has determinant one while multiplying the two canonical area blocks by four and one quarter.

In a single canonical plane the condition reduces to determinant one, so every orientation-preserving area-preserving linear map is symplectic. Hamiltonian time maps give many nonlinear examples, and composing symplectic maps preserves the property. A coordinate change can preserve ω while changing the expression of H; transformed trajectories are generated by the correspondingly transformed Hamiltonian. It does not follow that the numerical values of every observable stay unchanged. The form identifies an invariant structure, not a prohibition on changing coordinates or redistributing Euclidean widths.

Φ∗ω=ω,S⊤JS=J,J=(0I−I0)in the order (q,p)\Phi^*\omega=\omega,\qquad S^\top JS=J,\qquad J=\begin{pmatrix}0&I\\-I&0\end{pmatrix}\quad\text{in the order }(q,p)
The matrix identity preserves every canonical pairing. Determinant one follows, but in dimensions at least four it does not imply the full identity.↔ Scroll wide equations horizontally.

Differentiate a type-two generating function

Let F₂(q,P) depend on old positions and new momenta, and define p=∂F₂/∂q and Q=∂F₂/∂P. To solve for P as a function of (q,p), require the mixed Hessian F_{qP} to be invertible at the point. The implicit-function theorem then gives a local transformation. Taking the differential of F₂ yields p·dq+Q·dP. Rearranging shows P·dQ−p·dq is the exact differential of P·Q−F₂. Taking another exterior derivative proves the symplectic pullback identity.

The invertibility condition is not a decorative assumption. Without it, the equations may define a relation with several branches or a lower-dimensional constraint instead of a local map. Other types of generating functions use different choices of old and new independent coordinates. For a type-one function F(q,Q), the compatible signs are p=F_q and P=−F_Q. Those minus signs arise from the same tautological-form comparison, not from arbitrary convention. A map whose graph projects badly onto the chosen independent variables can require another generating type.

p=∂qF2(q,P),Q=∂PF2(q,P),det⁡FqP≠0,P⋅dQ−p⋅dq=d(P⋅Q−F2)p=\partial_qF_2(q,P),\quad Q=\partial_PF_2(q,P),\quad\det F_{qP}\ne0,\qquad P\cdot dQ-p\cdot dq=d(P\cdot Q-F_2)
The determinant condition gives a local map from (q,p) to (Q,P). The exact differential proves symplecticity in our ω=−dθ convention.↔ Scroll wide equations horizontally.

Local generating charts can fail globally

A quarter-turn Q=p,P=−q is symplectic, but old q and new P satisfy P=−q and cannot serve as independent variables for a type-two representation on an open set. A different generating chart handles it. Even where a local representation exists, global topology can obstruct a single-valued scalar potential for the tautological-form difference. The closed nonexact momentum translation on a circle from the previous lesson is a concrete example. A global exact symplectomorphism is therefore a stronger specification than a merely symplectic one.

Generating functions also organize numerical integrators. Solving derivative relations step by step can produce maps preserving the symplectic form exactly even when the original energy is not exactly preserved. This is a geometric design principle, not an unconditional accuracy or stability theorem. A coarse time step can still yield poor trajectories, implicit equations can have unsuitable branches, and a preserved form does not impose a desired phase error. Verification should include the intended limiting Hamiltonian, solvability and error behavior alongside the canonical-map identity.

Compose a kick and a drift through one function

Use F₂(q,P)=qP+q²+P²/4. Derive the map, verify its determinant and evaluate it at (q,p)=(1,3).

  1. Differentiation gives p=P+2q and Q=q+P/2. Thus P=p−2q and Q=p/2.

  2. In (q,p) order the map has matrix S=(01/2−21)S=\begin{pmatrix}0&1/2\\-2&1\end{pmatrix}, with determinant one. In this single phase plane that verifies symplecticity. The mixed derivative F_{qP}=1 also verifies local solvability everywhere.

  3. At (1,3), the new momentum is P=1 and the new position is Q=3/2. Substituting these values back into p=P+2q and Q=q+P/2 checks both defining equations.

The result

The map is a global linear symplectomorphism here. Its exact generating structure does not mean it preserves every chosen energy function.

Explore the model

Interactive lab

Interactive lab · matrix

A canonical shear in one phase plane

Pause and Reset, then set a=1,b=1. Predict (q,p)↦(q+p,p), determinant one and unchanged oriented area. This shear is generated by F₂=qP+P²/2. For a fresh comparison Reset and set a=0.5,b=0; predict area scale 0.5, so this map is not symplectic on a single fixed canonical plane. Readouts describe the selected final matrix during interpolation.

A canonical shear in one phase planeThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. 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 forThe same plotted family contains both canonical and noncanonical transformations; the form-preservation test decides which is which.

Check your understanding

Why is det F_{qP}≠0 required for this type-two generating construction?

Try it yourself

For F₂(q,P)=qP+P²/4, derive the transformation, verify the canonical form directly, and find the Hamiltonian whose time-one flow gives this shear.

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

Give me a hint

Differentiate first. Then compare the map with q̇=H_p,ṗ=0 for a momentum-only quadratic.

Compare with the solution

The equations give p=P and Q=q+P/2=q+p/2. Therefore dQ∧dP=(dq+dp/2)∧dp=dq∧dp. The Hamiltonian H=p²/4 has q̇=p/2 and ṗ=0, whose time-one map is exactly this shear. This particular generating map is an exact flow of that H; a general generating-function discretization need not be an exact flow of the originally specified Hamiltonian.

Take it with you

  • Symplectomorphisms preserve the full two-form; volume alone is weaker in higher dimensions.
  • A generating function requires a regular projection to define a local map.
  • Different generating types and global exactness address different representation issues.

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.