Contact Hamilton–Jacobi Equation Overview

• The contact Hamilton–Jacobi equation is a first-order PDE defined by u-dependence that captures nonconservative and dissipative dynamics.
• It replaces the classical symplectic framework with a contact structure, enabling geometric formulations and accommodating effects like discounting.
• The theory advances analysis via viscosity solutions, semigroup representations, and applications in thermodynamics, cosmology, and discrete/stochastic systems.

The contact Hamilton–Jacobi equation is a first-order partial differential equation in which the Hamiltonian depends not only on the base point and momentum variables, but also on the unknown function itself. In its evolutionary form it is written as

$u_t(x,t)+H(x,Du(x,t),u(x,t))=0,$

while the stationary form is

$H(x,Du(x),u(x))=c$

or, in the critical case studied in weak KAM theory,

$H(x,Du(x),u(x))=0.$

This $u$-dependence is the defining contact feature: it replaces the conservative symplectic Hamilton–Jacobi framework by a contact one on $T^*M\times \mathbb R$, introduces dissipative or nonconservative effects, and alters both the geometric and analytic structure of the theory (Ni et al., 2021). When $H$ is independent of $u$, one recovers the classical autonomous Hamilton–Jacobi equation; when the $u$-dependence is linear, one obtains discounted or conformally symplectic models (Wang et al., 2020).

1. Geometric formulation and contact dynamics

A standard contact phase space is $T^*Q\times \mathbb R$ with Darboux coordinates $(q^i,p_i,z)$ and contact form

$H(x,Du(x),u(x))=c$0

The associated Reeb field is

$H(x,Du(x),u(x))=c$1

and the contact Hamiltonian vector field of a Hamiltonian $H(x,Du(x),u(x))=c$2 has local form

$H(x,Du(x),u(x))=c$3

Its integral curves satisfy

$H(x,Du(x),u(x))=c$4

(León et al., 2021). In the notation used for manifolds $H(x,Du(x),u(x))=c$5, this same system is written as

$H(x,Du(x),u(x))=c$6

(Wang et al., 2020).

The geometric Hamilton–Jacobi problem asks for a section $H(x,Du(x),u(x))=c$7 whose image is compatible with the contact flow. For a Legendrian first jet $H(x,Du(x),u(x))=c$8, one formulation gives

$H(x,Du(x),u(x))=c$9

while in the broader contact Hamiltonian formulation one obtains modified equations involving the Reeb derivative of $H(x,Du(x),u(x))=0.$0 (Esen et al., 2022, León et al., 2021). In particular, the paper on contact Hamiltonian systems distinguishes between the contact Hamiltonian vector field $H(x,Du(x),u(x))=0.$1 and the evolution vector field

$H(x,Du(x),u(x))=0.$2

and these yield different Hamilton–Jacobi equations: in the Legendrian $H(x,Du(x),u(x))=0.$3-section formulation, $H(x,Du(x),u(x))=0.$4 leads to $H(x,Du(x),u(x))=0.$5, whereas $H(x,Du(x),u(x))=0.$6 leads to $H(x,Du(x),u(x))=0.$7 (León et al., 2021). This separation is specific to contact geometry and has no direct symplectic analogue.

The contact viewpoint is also compatible with generalized and implicit dynamics. One paper models implicit contact Hamiltonian systems as Legendrian submanifolds of an extended tangent contact manifold and derives an implicit contact Hamilton–Jacobi equation of the form

$H(x,Du(x),u(x))=0.$8

with additional Morse-family constraints

$H(x,Du(x),u(x))=0.$9

(Esen et al., 2021). This places singular Herglotz-type systems inside the same contact Hamilton–Jacobi framework.

2. PDE structure, viscosity solutions, and semigroup representation

For PDE purposes, a contact Hamilton–Jacobi equation is typically treated under continuity, convexity, coercivity in the momentum variable, and Lipschitz control in $u$0. A basic low-regularity setting assumes: $u$1 together with a uniform Lipschitz condition in $u$2 (Ni et al., 2021). Under such assumptions, the associated Lagrangian is the Legendre transform

$u$3

which may take the value $u$4 when only coercivity, rather than superlinearity, is assumed (Ni et al., 2021).

A central analytic object is the implicit backward Lax–Oleinik semigroup

$u$5

This is the contact analogue of the classical explicit Lax–Oleinik operator, but it is implicit because the running cost depends on the evolving value $u$6 itself (Ni et al., 2021). The corresponding forward semigroup $u$7 is defined dually. Under the standard assumptions, $u$8 is the unique continuous viscosity solution of

$u$9

(Ni et al., 2021).

The semigroup satisfies monotonicity and an exponential Lipschitz estimate: $T^*M\times \mathbb R$0 with the same type of bound for $T^*M\times \mathbb R$1 (Ni et al., 2021). Fixed points of $T^*M\times \mathbb R$2 are precisely stationary viscosity solutions: $T^*M\times \mathbb R$3 (Ni et al., 2021). In the weak KAM language, backward weak KAM solutions, viscosity solutions, and fixed points of $T^*M\times \mathbb R$4 coincide; similarly, forward weak KAM solutions correspond to fixed points of $T^*M\times \mathbb R$5 (Ni et al., 2021, Ni et al., 27 Apr 2026).

For the Cauchy problem on $T^*M\times \mathbb R$6,

$T^*M\times \mathbb R$7

one also obtains a solution semigroup $T^*M\times \mathbb R$8 with a variational formula

$T^*M\times \mathbb R$9

where $H$0 is defined by an implicit variational principle (Jin, 20 Mar 2026). This Euclidean representation supports semigroup methods, commutation results, and multi-time compatibility.

3. Weak KAM theory, critical values, and solution selection

The stationary contact Hamilton–Jacobi equation is closely linked to weak KAM theory on compact manifolds. In the generalized stationary equation

$H$1

one studies the set of constants $H$2 for which viscosity solutions exist. For $H$3-independent Hamiltonians $H$4, the critical value is characterized by

$H$5

(Ni et al., 27 Apr 2026). In the contact case, especially when $H$6 is merely non-decreasing in $H$7, the solvability set

$H$8

is a connected interval rather than a singleton (Liu et al., 22 Sep 2025).

For Hamiltonians that are convex and superlinear in $H$9 and strictly increasing in $u$0, a distinguished vanishing-contact limit was established for the family

$u$1

where $u$2 is the Mañé critical value of the conservative Hamiltonian $u$3. The viscosity solution $u$4 converges uniformly, as $u$5, to a uniquely selected critical solution $u$6 of

$u$7

(Wang et al., 2020). This limit is characterized in two equivalent ways.

The first is via weighted Mather inequalities. Defining

$u$8

the selected limit is

$u$9

(Wang et al., 2020). The second uses the Peierls barrier $u$0 of the conservative Hamiltonian: $u$1 In the discounted case $u$2, hence $u$3, this reduces to the known formula

$u$4

(Wang et al., 2020). This shows that the contact dependence contributes to the limit through the first $u$5-jet $u$6.

When the $u$7-dependence is only non-decreasing, classical comparison can fail. A generalized comparison principle replaces pointwise monotonicity by integral inequalities over special Mather measures. If $u$8 solve

$u$9

and

$T^*Q\times \mathbb R$0

then $T^*Q\times \mathbb R$1 on $T^*Q\times \mathbb R$2 (Liu et al., 22 Sep 2025). Here $T^*Q\times \mathbb R$3 denotes the “ordinal” Mather measures satisfying

$T^*Q\times \mathbb R$4

This identifies the degeneracy of $T^*Q\times \mathbb R$5-monotonicity, rather than convexity in $T^*Q\times \mathbb R$6 alone, as the source of nonuniqueness.

4. Stability, long-time behavior, and perturbed equations

Lyapunov stability of stationary solutions has become a central topic for contact-type Hamilton–Jacobi equations. For the evolutionary equation

$T^*Q\times \mathbb R$7

on a compact manifold, a stationary solution $T^*Q\times \mathbb R$8 is Lyapunov stable if small $T^*Q\times \mathbb R$9-perturbations remain small under the backward semigroup $(q^i,p_i,z)$0, and asymptotically stable if they converge back to $(q^i,p_i,z)$1 as $(q^i,p_i,z)$2 (Ni et al., 27 Apr 2026).

A recent PDE formulation treats separated Hamiltonians

$(q^i,p_i,z)$3

with $(q^i,p_i,z)$4 continuous, convex, and coercive in $(q^i,p_i,z)$5, and $(q^i,p_i,z)$6 continuous and bounded (Ni et al., 27 Apr 2026). The key object is the shifted critical value

$(q^i,p_i,z)$7

Its one-sided derivatives satisfy

$(q^i,p_i,z)$8

$(q^i,p_i,z)$9

From this, PDE-type criteria are obtained. If one of the conditions

$H(x,Du(x),u(x))=c$00

or

$H(x,Du(x),u(x))=c$01

or

$H(x,Du(x),u(x))=c$02

holds, then $H(x,Du(x),u(x))=c$03 is locally asymptotically stable, with exponential decay rate in the sense

$H(x,Du(x),u(x))=c$04

for some $H(x,Du(x),u(x))=c$05 (Ni et al., 27 Apr 2026). Dual sign conditions imply Lyapunov instability.

A one-dimensional non-monotone theory on the circle identifies a scalar averaged quantity

$H(x,Du(x),u(x))=c$06

associated with the unique periodic Aubry orbit (Xu et al., 2024). If $H(x,Du(x),u(x))=c$07, the stationary solution $H(x,Du(x),u(x))=c$08 is asymptotically stable on a local order interval and

$H(x,Du(x),u(x))=c$09

If $H(x,Du(x),u(x))=c$10, $H(x,Du(x),u(x))=c$11 is Lyapunov unstable, and the evolutionary equation has infinitely many nontrivial time-periodic viscosity solutions (Xu et al., 2024). This replaces a global sign assumption on $H(x,Du(x),u(x))=c$12 by an orbit-average criterion.

Perturbation theory fits the same pattern. If $H(x,Du(x),u(x))=c$13 is a locally Lyapunov asymptotically stable viscosity solution of

$H(x,Du(x),u(x))=c$14

then for the perturbed stationary equation

$H(x,Du(x),u(x))=c$15

there exist viscosity solutions $H(x,Du(x),u(x))=c$16 with

$H(x,Du(x),u(x))=c$17

(Wu et al., 9 Jan 2025). Under the stronger positivity condition

$H(x,Du(x),u(x))=c$18

the nearby perturbed solution is locally unique and remains locally asymptotically stable (Wu et al., 9 Jan 2025). A key semigroup estimate is

$H(x,Du(x),u(x))=c$19

which quantitatively controls the perturbation (Wu et al., 9 Jan 2025).

5. Commutation, discrete formulations, and stochastic extensions

Contact Hamilton–Jacobi theory also supports algebraic and geometric extensions beyond single-time deterministic flows. For two contact Hamiltonians $H(x,Du(x),u(x))=c$20 on $H(x,Du(x),u(x))=c$21, the correct compatibility condition is not Poisson commutation but the vanishing of the contact or Jacobi bracket: $H(x,Du(x),u(x))=c$22 Under the structural assumptions of convexity, coercivity in $H(x,Du(x),u(x))=c$23, and uniform Lipschitz dependence in $H(x,Du(x),u(x))=c$24, this implies semigroup commutation: $H(x,Du(x),u(x))=c$25 for all $H(x,Du(x),u(x))=c$26 (Jin, 20 Mar 2026). Pairwise Jacobi commutation of several Hamiltonians yields existence and uniqueness for multi-time systems

$H(x,Du(x),u(x))=c$27

(Jin, 20 Mar 2026).

A discrete Hamilton–Jacobi theory for contact dynamics has also been developed. Starting from a discrete contact Lagrangian $H(x,Du(x),u(x))=c$28, one defines discrete Legendre transforms and right discrete contact Hamilton equations

$H(x,Du(x),u(x))=c$29

The associated discrete contact Hamilton–Jacobi equation is

$H(x,Du(x),u(x))=c$30

(Esen et al., 2022). Its geometric meaning is the commutation relation

$H(x,Du(x),u(x))=c$31

between the discrete contact flow and the discrete section (Esen et al., 2022).

In the stochastic setting, a stochastic contact Hamilton–Jacobi equation is introduced for Stratonovich stochastic contact systems: $H(x,Du(x),u(x))=c$32 Together with

$H(x,Du(x),u(x))=c$33

this generates the stochastic contact Hamiltonian flow (Zhan et al., 2024). The same paper constructs contact structure-preserving numerical schemes from truncated expansions of the stochastic generating function and proves mean-square convergence order $H(x,Du(x),u(x))=c$34 for the resulting one-step approximation (Zhan et al., 2024).

6. Applications and specialized variants

One prominent application is thermodynamics. On a thermodynamic contact manifold with coordinates $H(x,Du(x),u(x))=c$35, contact Hamilton–Jacobi theory yields

$H(x,Du(x),u(x))=c$36

where $H(x,Du(x),u(x))=c$37 is a thermodynamic potential on the equilibrium-state space (Ghosh, 2022). In the energy representation, for an ideal gas, one has

$H(x,Du(x),u(x))=c$38

and the contact HJ equation reproduces characteristic equations for isochoric and isothermal thermodynamic transformations (Ghosh, 2022). A related geometric HJ formulation on $H(x,Du(x),u(x))=c$39 uses Legendrian sections $H(x,Du(x),u(x))=c$40 and the condition

$H(x,Du(x),u(x))=c$41

(Esen et al., 2022).

Another application arises in cosmology. For strong sustained rapid turn inflation in two-field models, after fixing the scalar potential $H(x,Du(x),u(x))=c$42 and the conformal class of the scalar field metric, the strong-SRRT consistency condition becomes a geometric contact Hamilton–Jacobi equation on the first jet bundle $H(x,Du(x),u(x))=c$43: $H(x,Du(x),u(x))=c$44 In local isothermal Liouville coordinates $H(x,Du(x),u(x))=c$45, this takes the form

$H(x,Du(x),u(x))=c$46

where $H(x,Du(x),u(x))=c$47 is the unknown metric in the chosen conformal class (Babalic et al., 2024). The paper analyzes this as a proper first-order nonlinear PDE, studies its characteristic flow, and extracts asymptotic conditions near nondegenerate critical points of the potential (Babalic et al., 2024).

A further specialized direction is the sign-changing discounted equation

$H(x,Du(x),u(x))=c$48

which is a contact Hamilton–Jacobi equation with affine $H(x,Du(x),u(x))=c$49-dependence. When $H(x,Du(x),u(x))=c$50 changes sign, the set of admissible constants is

$H(x,Du(x),u(x))=c$51

and the stationary problem admits maximal and minimal viscosity solutions (Ni et al., 2022). Large-time behavior of the evolutionary equation

$H(x,Du(x),u(x))=c$52

is then organized around extremal weak KAM solutions and semigroup monotonicity (Ni et al., 2022). This suggests that even the linear-in-$H(x,Du(x),u(x))=c$53 subclass already exhibits several distinctly contact phenomena: failure of global comparison, multiplicity of stationary solutions, and basin-dependent asymptotics.

Across these variants, the contact Hamilton–Jacobi equation retains a common structural role. It is simultaneously a geometric reduction principle, a viscosity-solution PDE with semigroup representation, a weak KAM selection mechanism, and a framework for nonconservative dynamics. What changes from problem to problem is the dominant structure: in Tonelli settings, Aubry–Mather theory and Peierls barriers are central; in low-regularity PDE settings, critical values and strict subsolutions dominate; in discrete, stochastic, and thermodynamic settings, the contact generating-function interpretation becomes primary (Wang et al., 2020, Ni et al., 27 Apr 2026, Esen et al., 2022).