Legendre transformation

From Wikipedia, the free encyclopedia

Diagram illustrating the Legendre transformation of the function f(x) . The function is shown in red, and the tangent line at x0  is shown in blue. The tangent line intersects the vertical axis at (0, −f*) and f*  is the value of the Legendre transform f*(p)  where . Note that for any other point on the red curve, a line drawn through that point with the same slope as the blue line will have a y-intercept above the point (0, −f*), showing that f* is indeed a maximum.
Enlarge
Diagram illustrating the Legendre transformation of the function f(x) . The function is shown in red, and the tangent line at x0  is shown in blue. The tangent line intersects the vertical axis at (0, −f*) and f*  is the value of the Legendre transform f*(p)  where p=\dot{f}(x_0). Note that for any other point on the red curve, a line drawn through that point with the same slope as the blue line will have a y-intercept above the point (0, −f*), showing that f* is indeed a maximum.

In mathematics, it is often desirable to express a functional relationship f(x)\, as a different function, whose argument is the derivative of f , rather than x . If we let p = df/dx  be the argument of this new function, then this new function is written f^\star(p)\, and is called the Legendre transform of the original function, after Adrien-Marie Legendre.

The Legendre transform f^\star of a function f\, is defined as follows:

f^\star(p) = \mathrm{max}_x(px-f(x)).

The notation maxx indicates the maximization of the expression with respect to the variable x while p is held constant. The Legendre transform is its own inverse. Like the familiar Fourier transform, the Legendre transform takes a function f(x) and produces a function of a different variable p. However, while the Fourier transform consists of an integration with a kernel, the Legendre transform uses maximization as the transformation procedure. The transform is well behaved only if f(x) is a convex function:

\frac{d^2f}{dx^2}> 0

The Legendre transformation is an application of the duality relationship between points and lines. The functional relationship specified by f(x) can be represented equally well as a set of (x, y) points, or as a set of tangent lines specified by their slope and intercept values.

The Legendre transformation can be generalized to the Legendre-Fenchel transformation. It is commonly used in thermodynamics and in the hamiltonian formulation of classical mechanics.

Contents

[edit] Definitions

The definition of the Legendre transform can be made more explicit. To maximize f * with respect to x, we set its derivative equal to zero:

\frac{d}{dx} \left(f^* \right) = 0 \,
\frac{d}{dx} \left(xp-f(x) \right) = p-{df(x) \over dx} = 0. \quad \quad (1)\,

Thus, the expression is maximized when

p = {df(x) \over dx} \quad \quad \quad \quad \quad \quad (2).

This is a maximum because the second derivative is negative:

{d^2 \over dx^2}(xp-f(x)) = -{d^2f(x) \over dx^2} < 0,

since f was assumed convex. Next we invert (2) to obtain x as a function of p and plug this into (1) , which gives the more useful form,

f^\star(p) = p \,\, x(p) - f(x(p)).

This definition gives the conventional procedure for calculating the Legendre transform of f(x): find p = {df \over dx}, invert for x and substitute into the expression xpf(x). This definition makes clear the following interpretation: the Legendre transform produces a new function, in which the independent variable x is replaced by p = {df \over dx}, which is the derivative of the original function with respect to x.

[edit] Another definition

There is a third definition of the Legendre transform: f\, and f^\star are said to be Legendre transforms of each other if their first derivatives are inverse functions of each other:

Df = \left( Df^\star \right)^{-1}.

We can see this by taking derivative of f^\star:

{df^\star(p) \over dp} = {d \over dp}(xp-f(x)) = x.

Combining this equation with the maximization condition results in the following pair of reciprocal equations:

p = {df \over dx}(x),
x = {df^\star \over dp}(p).

We see that Df and Df^\star are inverses, as promised. They are unique up to an additive constant which is fixed by the additional requirement that

f(x) + f^\star(y) = x\,y.

Although in some cases (e.g. thermodynamic potentials) a non-standard requirement is used:

f(x) - f^\star(y) = x\,y.

The standard constraint will be considered in this article unless otherwise noted. The Legendre transformation is its own inverse, and is related to integration by parts.

[edit] Applications

[edit] Thermodynamics

The strategy behind the use of Legendre transforms is to shift the dependence of a function from one independent variable to another (the derivative of the original function with regard to this independent variable) by taking the difference between the original function and their product. They are used to transform among the various thermodynamic potentials. For example, while the internal energy is an explicit function of the extensive variables, entropy, volume (and chemical composition)

U = U(S,V,\{N_i\})\,

the enthalpy, the (non standard) Legendre transform of U with respect to −PV

H = U + PV \, = H(S,P,\{N_i\})\,
P=\, -\left( \frac{\partial U}{\partial V}\right)_S\,

becomes a function of the entropy and the intensive quantity, pressure, as natural variables, and is useful when the (external) P is constant. The free energies (Helmholtz and Gibbs), are obtained through further Legendre transforms, by subtracting TS (from U and H respectively), shift dependence from the entropy S to its conjugate intensive variable temperature T, and are useful when it is constant.

[edit] Hamilton-Lagrange mechanics

A Legendre transform is used in classical mechanics to derive the Hamiltonian formulation from the Lagrangian one, and conversely. While the Lagrangian is an explicit function of the positional coordinates qj and generalized velocities dqj /dt (and time), the Hamiltonian shifts the functional dependence to the positions and momenta,defined as p_j=\frac{\partial L}{\partial \dot{q}_j}

H\left(q_j,p_j,t\right) = \sum_i \dot{q}_i p_i - L(q_j,\dot{q}_j,t) \,.

Each of the two formulations has its own applicability, both in the theoretical foundations of the subject, and in practice, depending on the ease of calculation for a particular problem. The coordinates are not necessarily rectilinear, but can also be angles, etc. An optimum choice takes advantage of the actual physical symmetries.

[edit] An example - variable capacitor

As another example from physics, consider a parallel-plate capacitor whose plates can approach or recede from one another, exchanging work with external mechanical forces which maintain the plate separation — analogous to a gas in a cylinder with a piston. We want the attractive force f between the plates as a function of the variable separation x. (The two vectors point in opposite directions.) If the charges on the plates remain constant as they move, the force is the negative gradient of the electrostatic energy

U (Q, \mathbf{x} ) = \begin{matrix} \frac{1}{2} \end{matrix} QV \,.

However, if the voltage between the plates V is maintained constant by connection to a battery, which is a reservoir for charge at constant potential difference, the force now becomes the negative gradient of the Legendre transform

U - QV = -\begin{matrix} \frac{1}{2} \end{matrix} QV \,.

The two functions happen to be negatives only because of the linearity of the capacitance. Of course, for given charge, voltage and distance, the static force must be the same by either calculation since the plates cannot "know" what might be held constant as they move.


[edit] Examples

The exponential function ex has  x ln x − x  as a Legendre transform since the respective first derivatives ex and ln x are inverse to each other. This example shows that the respective domains of a function and its Legendre transform need not agree.

Similarly, the quadratic form

f(x) = \begin{matrix} \frac{1}{2} \end{matrix} \, x^t \, A \, x

with A a symmetric invertible n-by-n-matrix has

f^\star(y) = \begin{matrix} \frac{1}{2} \end{matrix} \, y^t \, A^{-1} \, y

as a Legendre transform.

[edit] Legendre transformation in one dimension

In one dimension, a Legendre transform to a function f : R → R with an invertible first derivative may be found using the formula

f^\star(y) = y \, x - f(x), \, x = \dot{f}^{-1}(y)

This can be seen by integrating both sides of the defining condition restricted to one-dimension

\dot{f}(x) = \dot{f}^{\star-1}(x)

from x0 to x1, making use of the fundamental theorem of calculus on the left hand side and substituting

y = \dot{f}^{\star-1}(x)

on the right hand side to find

f(x_1) - f(x_0) = \int_{y_0}^{y_1} y \, \ddot{f}^\star(y) \, dy

with f*′(y0) = x0, f*′(y1) = x1. Using integration by parts the last integral simplifies to

y_1 \, \dot{f}^\star(y_1) - y_0 \, \dot{f}^\star(y_0) - \int_{y_0}^{y_1} \dot{f}^\star(y) \, dy = y_1 \, x_1 - y_0 \, x_0 - f^\star(y_1) + f^\star(y_0).

Therefore,

f(x_1) + f^\star(y_1) - y_1 \, x_1 = f(x_0) + f^\star(y_0) - y_0 \, x_0.

Since the left hand side of this equation does only depend on x1 and the right hand side only on x0, they have to evaluate to the same constant.

f(x) + f^\star(y) - y \, x = C,\, x = \dot{f}^\star(y) = \dot{f}^{-1}(y).

Solving for f* and choosing C to be zero results in the above-mentioned formula.

[edit] Geometric interpretation

For a strictly convex function the Legendre-transformation can be interpreted as a mapping between the graph of the function and the family of tangents of the graph. (For a function of one variable, the tangents are well-defined at all but at most countably many points since a convex function is differentiable at all but at most countably many points.)

The equation of a line with slope m and y-intercept b is given by

y = mx + b\,.

For this line to be tangent to the graph of a function f at the point (x0, f(x0)) requires

f\left(x_0\right) = m x_0 + b

and

m = \dot{f}\left(x_0\right)

f' is strictly monotone as the derivative of a strictly convex function, and the second equation can be solved for x0, allowing to eliminate x0 from the first giving the y-intercept b of the tangent as a function of its slope m:

b = f\left(\dot{f}^{-1}\left(m\right)\right) - m \cdot \dot{f}^{-1}\left(m\right) = -f^\star(m)

Here f* denotes the Legendre transform of f.

The family of tangents of the graph of f is therefore (parameterized by m) given by

y = mx - f^\star(m)

or, written implicitly, by the solutions of the equation

F(x,y,m) = y + f^\star(m) - mx = 0.

The graph of the original function can be reconstructed from this family of lines as the envelope of this family by demanding

{\partial F(x,y,m)\over\partial m} = \dot{f}^\star(m) - x = 0.

Eliminating m from these two equations gives

y = x \cdot \dot{f}^{\star-1}(x) - f^\star\left(\dot{f}^{\star-1}(x)\right).

Identifying y with f(x) and recognizing the right side of the preceding equation as the Legendre transform of f* we find

f(x) = f^{\star\star}(x).

[edit] Legendre transformation in more than one dimension

For a differentiable real-valued function on an open subset U of Rn the Legendre conjugate of the pair (U, f) is defined to be the pair (V, g), where V is the image of U under the gradient mapping Df, and g is the function on V given by the formula

g(y) = \left\langle y, x \right\rangle - f\left(x\right), \, x = \left(Df\right)^{-1}(y)

where

\left\langle u,v\right\rangle = \sum_{k=1}^{n}u_{k} \cdot v_{k}

is the scalar product on Rn.

Alternatively, if X is a real vector space and Y is its dual vector space, then for each point x of X and y of Y, there is a natural identification of the cotangent spaces T*Xx with Y and T*Yy with X. If f is a real differentiable function over X, then ∇f is a section of the cotangent bundle T*X and as such, we can construct a map from X to Y. Similarly, if g is a real differentiable function over Y, ∇g defines a map from Y to X. If both maps happen to be inverses of each other, we say we have a Legendre transform.

[edit] Further properties

In the following the Legendre transform of a function f is denoted as f*.

[edit] Scaling properties

The Legendre transformation has the following scaling properties:

f(x) = a \cdot g(x) \Rightarrow f^\star(p) = a \cdot g^\star\left(\frac{p}{a}\right)
f(x) = g(a \cdot x) \Rightarrow f^\star(p) = g^\star\left(\frac{p}{a}\right)

It follows that if a function is homogeneous of degree r then its image under the Legendre transformation is a homogeneous function of degree s, where 1/r + 1/s = 1.

[edit] Behavior under translation

f(x) = g(x) + b \Rightarrow f^\star(p) = g^\star(p) - b
f(x) = g(x + y) \Rightarrow f^\star(p) = g^\star(p) - p \cdot y

[edit] Behavior under inversion

f(x) = g^{-1}(x) \Rightarrow f^\star(p) = - p \cdot g^\star\left(\frac{1}{p}\right)

[edit] Behavior under linear transformations

Let A be a linear transformation from Rn to Rm. For any convex function f on Rn, one has

\left(A f\right)^\star = f^\star A^\star

where A* is the adjoint operator of A defined by

\left \langle Ax, y^\star \right \rangle = \left \langle x, A^\star y^\star \right \rangle

A closed convex function f is symmetric with respect to a given set G of orthogonal linear transformations,

f\left(A x\right) = f(x), \; \forall x, \; \forall A \in G

if and only if f* is symmetric with respect to G.

[edit] Infimal convolution

The infimal convolution of two functions f and g is defined as

\left(f \star_\inf  g\right)(x) = \inf \left \{ f(x-y) + g(y) \, | \, y \in \mathbb{R}^n \right \}

Let f1, …, fm be proper convex functions on Rn. Then

\left( f_1 \star_\inf \cdots \star_\inf f_m \right)^\star = f_1^\star + \cdots + f_m^\star

[edit] See also

[edit] References