Functional derivative

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In the calculus of variations, a field of mathematical analysis, the functional derivative (or variational derivative)[1] relates a change in a functional to a change in a function that the functional depends on.

In the calculus of variations, functionals are usually expressed in terms of an integral of functions, their arguments, and their derivatives. In an integrand L of a functional, if a function f is varied by adding to it another function δf that is arbitrarily small, and the resulting L is expanded in powers of δf, the coefficient of δf in the first order term is called the functional derivative.

For example, consider the functional

 J[f] = \int \limits_a^b L[ \, x, f(x), f \, '(x) \, ] \, dx \ ,

where f ′(x) ≡ df/dx. If f is varied by adding to it a function δf, and the resulting integrand L(x, f +δf, f '+δf ′) is expanded in powers of δf, then the change in the value of J to first order in δf can be expressed as follows:[1][Note 1]

  \delta J = \int_a^b  \frac{\delta J}{\delta f(x)} {\delta f(x)} dx \, .

The coefficient of δf(x), denoted as δJ/δf(x), is called the functional derivative of J with respect to f at the point x.[2] For this example functional, the functional derivative is the left hand side of the Euler-Lagrange equation,[3]

 \frac{\delta J}{\delta f(x)} = \frac{\part L}{\part f} -\frac{d}{dx} \frac{\part L}{\part f'} \, .


The functional derivative is defined. Then the functional differential is defined in terms of the functional derivative.

Functional derivative[edit]

Given a manifold M representing (continuous/smooth/with certain boundary conditions/etc.) functions ρ and a functional F defined as

F\colon M \rightarrow \mathbb{R} \quad \mbox{or} \quad F\colon M \rightarrow \mathbb{C} \, ,

the functional derivative of F[ρ], denoted δF/δρ, is defined by[4]

  \int \frac{\delta F}{\delta\rho(x)} \phi(x) dx 
&= \lim_{\varepsilon\to 0}\frac{F[\rho+\varepsilon \phi]-F[\rho]}{\varepsilon} \\
&= \left [ \frac{d}{d\epsilon}F[\rho+\epsilon \phi]\right ]_{\epsilon=0},

where \epsilon \phi is the variation of ρ, and \phi is an arbitrary function.

Functional differential[edit]

The differential (or variation or first variation) of the functional F[ρ] is,[5] [Note 2]

  \delta F = \int   \frac {\delta F} {\delta \rho(x)} \ \delta \rho(x)   \ dx \ ,

where δρ(x) = εϕ(x) is the variation of ρ(x). This is similar in form to the total differential of a function F(ρ1, ρ2, ..., ρn),

 dF =   \sum_{i=1} ^n  \frac {\partial F} {\partial \rho_i} \ d\rho_i  \ ,

where ρ1, ρ2, ... , ρn are independent variables. Comparing the last two equations, the functional derivative δF/δρ(x) has a role similar to that of the partial derivative ∂F/∂ρi , where the variable of integration x is like a continuous version of the summation index i.[2]

Formal description[edit]

The definition of a functional derivative may be made more mathematically precise and formal by defining the space of functions more carefully. For example, when the space of functions is a Banach space, the functional derivative becomes known as the Fréchet derivative, while one uses the Gâteaux derivative on more general locally convex spaces. Note that the well-known Hilbert spaces are special cases of Banach spaces. The more formal treatment allows many theorems from ordinary calculus and analysis to be generalized to corresponding theorems in functional analysis, as well as numerous new theorems to be stated.


Like the derivative of a function, the functional derivative satisfies the following properties, where F[ρ] and G[ρ] are functionals:

\frac{\delta(\lambda F + \mu G)}{\delta \rho(x)} = \lambda \frac{\delta F}{\delta \rho(x)} + \mu \frac{\delta G}{\delta \rho(x)},\ \qquad  \lambda,\mu   constant,
  • Product rule:[7]
\frac{\delta(FG)}{\delta \rho(x)} = \frac{\delta F}{\delta \rho(x)} G + F \frac{\delta G}{\delta \rho(x)} \, ,
  • Chain rules:
If f is a differentiable function, then
\displaystyle\frac{\delta F[f(\rho)] }{\delta\rho(x)}  = \frac{\delta F[f(\rho)]}{\delta f(\rho(x) )} \ \frac {df(\rho(x))} {d\rho(x)} \ , [8]
\frac{\delta f( F[\rho])}{\delta\rho(x)} = \frac {df(F[\rho])} {dF[\rho]} \ \frac{\delta F[\rho]}{\delta\rho(x)} \, . [9]

Determining functional derivatives[edit]

We give a formula to determine functional derivatives for a common class of functionals that can be written as the integral of a function and its derivatives. This is a generalization of the Euler–Lagrange equation: indeed, the functional derivative was introduced in physics within the derivation of the Lagrange equation of the second kind from the principle of least action in Lagrangian mechanics (18th century). The first three examples below are taken from density functional theory (20th century), the fourth from statistical mechanics (19th century).


Given a functional

F[\rho] = \int f( \boldsymbol{r}, \rho(\boldsymbol{r}), \nabla\rho(\boldsymbol{r}) )\, d\boldsymbol{r},

and a function ϕ(r) that vanishes on the boundary of the region of integration, from a previous section Definition,

\int \frac{\delta F}{\delta\rho(\boldsymbol{r})} \, \phi(\boldsymbol{r}) \, d\boldsymbol{r}  
&  = \left [ \frac{d}{d\varepsilon} \int f( \boldsymbol{r}, \rho + \varepsilon \phi, \nabla\rho+\varepsilon\nabla\phi )\, d\boldsymbol{r} \right ]_{\varepsilon=0} \\
&  = \int \left( \frac{\partial f}{\partial\rho} \, \phi + \frac{\partial f}{\partial\nabla\rho} \cdot \nabla\phi \right) d\boldsymbol{r} \\
&  = \int \left[ \frac{\partial f}{\partial\rho} \, \phi + \nabla \cdot \left( \frac{\partial f}{\partial\nabla\rho} \, \phi \right) - \left( \nabla \cdot \frac{\partial f}{\partial\nabla\rho} \right) \phi \right] d\boldsymbol{r} \\
&  = \int \left[ \frac{\partial f}{\partial\rho} \, \phi - \left( \nabla \cdot \frac{\partial f}{\partial\nabla\rho} \right) \phi \right] d\boldsymbol{r} \\
&  = \int \left( \frac{\partial f}{\partial\rho} -  \nabla \cdot \frac{\partial f}{\partial\nabla\rho} \right) \phi(\boldsymbol{r}) \ d\boldsymbol{r} \, .

The second line is obtained using the total derivative, where ∂f /∂∇ρ is a derivative of a scalar with respect to a vector.[Note 3] The third line was obtained by use of a product rule for divergence. The fourth line was obtained using the divergence theorem and the condition that ϕ=0 on the boundary of the region of integration. Since ϕ is also an arbitrary function, applying the fundamental lemma of calculus of variations to the last line, the functional derivative is

\frac{\delta F}{\delta\rho(\boldsymbol{r})} = \frac{\partial f}{\partial\rho} - \nabla \cdot \frac{\partial f}{\partial\nabla\rho}

where ρ = ρ(r) and f = f (r, ρ, ∇ρ). This formula is for the case of the functional form given by F[ρ] at the beginning of this section. For other functional forms, the definition of the functional derivative can be used as the starting point for its determination. (See the example Coulomb potential energy functional.)

The above equation for the functional derivative can be generalized to the case that includes higher dimensions and higher order derivatives. The functional would be,

F[\rho(\boldsymbol{r})] = \int f( \boldsymbol{r}, \rho(\boldsymbol{r}), \nabla\rho(\boldsymbol{r}), \nabla^{(2)}\rho(\boldsymbol{r}), \dots, \nabla^{(N)}\rho(\boldsymbol{r}))\, d\boldsymbol{r},

where the vector r ∈ ℝn, and (i) is a tensor whose ni components are partial derivative operators of order i,

 \left [ \nabla^{(i)} \right ]_{\alpha_1 \alpha_2 \cdots \alpha_i} = \frac {\partial^{\, i}} {\partial r_{\alpha_1}  \partial r_{\alpha_2} \cdots \partial r_{\alpha_i} } \qquad \qquad \text{where} \quad  \alpha_1, \alpha_2, \cdots, \alpha_i = 1, 2, \cdots , n \ . [Note 4]

An analogous application of the definition of the functional derivative yields

\frac{\delta F[\rho]}{\delta \rho} &{} = \frac{\partial f}{\partial\rho} - \nabla \cdot \frac{\partial f}{\partial(\nabla\rho)} + \nabla^{(2)} \cdot \frac{\partial f}{\partial\left(\nabla^{(2)}\rho\right)} + \dots + (-1)^N \nabla^{(N)} \cdot \frac{\partial f}{\partial\left(\nabla^{(N)}\rho\right)} \\
&{} =   \frac{\partial f}{\partial\rho} + \sum_{i=1}^N (-1)^{i}\nabla^{(i)} \cdot \frac{\partial f}{\partial\left(\nabla^{(i)}\rho\right)} \ .

In the last two equations, the ni components of the tensor  \frac{\partial f}{\partial\left(\nabla^{(i)}\rho\right)} are partial derivatives of f with respect to partial derivatives of ρ,

 \left [ \frac {\partial f} {\partial \left (\nabla^{(i)}\rho \right ) } \right ]_{\alpha_1 \alpha_2 \cdots \alpha_i} = \frac {\partial f} {\partial \rho_{\alpha_1 \alpha_2 \cdots \alpha_i} } \qquad \qquad \text{where} \quad \rho_{\alpha_1 \alpha_2 \cdots \alpha_i} \equiv \frac {\partial^{\, i}\rho} {\partial r_{\alpha_1} \, \partial r_{\alpha_2} \cdots \partial r_{\alpha_i} }   \ ,

and the tensor scalar product is,

 \nabla^{(i)} \cdot \frac{\partial f}{\partial\left(\nabla^{(i)}\rho\right)} = \sum_{\alpha_1, \alpha_2, \cdots, \alpha_i = 1}^n \ \frac {\partial^{\, i} } {\partial r_{\alpha_1} \, \partial r_{\alpha_2} \cdots \partial r_{\alpha_i} }  \ \frac {\partial f} {\partial \rho_{\alpha_1 \alpha_2 \cdots \alpha_i} }   \  .  [Note 5]


Thomas–Fermi kinetic energy functional[edit]

The Thomas–Fermi model of 1927 used a kinetic energy functional for a noninteracting uniform electron gas in a first attempt of density-functional theory of electronic structure:

T_\mathrm{TF}[\rho] = C_\mathrm{F} \int \rho^{5/3}(\mathbf{r}) \, d\mathbf{r} \, .

Since the integrand of TTF[ρ] does not involve derivatives of ρ(r), the functional derivative of TTF[ρ] is,[10]

\frac{\delta T_{\mathrm{TF}}}{\delta \rho (\boldsymbol{r}) } 
& = C_\mathrm{F} \frac{\partial \rho^{5/3}(\mathbf{r})}{\partial \rho(\mathbf{r})}  \\
& = \frac{5}{3} C_\mathrm{F}  \rho^{2/3}(\mathbf{r}) \, .

Coulomb potential energy functional[edit]

For the electron-nucleus potential, Thomas and Fermi employed the Coulomb potential energy functional

V[\rho] =  \int \frac{\rho(\boldsymbol{r})}{|\boldsymbol{r}|} \ d\boldsymbol{r}.

Applying the definition of functional derivative,

\int \frac{\delta V}{\delta \rho(\boldsymbol{r})} \ \phi(\boldsymbol{r}) \ d\boldsymbol{r} 
& {} = \left [ \frac{d}{d\varepsilon}  \int \frac{\rho(\boldsymbol{r}) + \varepsilon \phi(\boldsymbol{r})}{|\boldsymbol{r}|} \ d\boldsymbol{r} \right ]_{\varepsilon=0} \\
& {} =  \int  \frac {1} {|\boldsymbol{r}|} \, \phi(\boldsymbol{r}) \ d\boldsymbol{r} \, .


 \frac{\delta V}{\delta \rho(\boldsymbol{r})} = \frac{1}{|\boldsymbol{r}|} \ .

For the classical part of the electron-electron interaction, Thomas and Fermi employed the Coulomb potential energy functional

J[\rho] = \frac{1}{2}\iint \frac{\rho(\mathbf{r}) \rho(\mathbf{r}')}{\vert \mathbf{r}-\mathbf{r}' \vert}\, d\mathbf{r} d\mathbf{r}' \, .

From the definition of the functional derivative,

\int \frac{\delta J}{\delta\rho(\boldsymbol{r})} \phi(\boldsymbol{r})d\boldsymbol{r}  
& {} = \left [ \frac {d \ }{d\epsilon} \, J[\rho + \epsilon\phi] \right ]_{\epsilon = 0} \\
& {} = \left [ \frac {d \ }{d\epsilon} \, \left ( \frac{1}{2}\iint \frac {[\rho(\boldsymbol{r}) + \epsilon \phi(\boldsymbol{r})] \,  [\rho(\boldsymbol{r}') + \epsilon \phi(\boldsymbol{r}')]  }{\vert \boldsymbol{r}-\boldsymbol{r}' \vert}\, d\boldsymbol{r} d\boldsymbol{r}'  \right ) \right ]_{\epsilon = 0} \\
& {} =  \frac{1}{2}\iint \frac {\rho(\boldsymbol{r}') \phi(\boldsymbol{r})  }{\vert \boldsymbol{r}-\boldsymbol{r}' \vert}\, d\boldsymbol{r} d\boldsymbol{r}' +
            \frac{1}{2}\iint \frac {\rho(\boldsymbol{r}) \phi(\boldsymbol{r}')  }{\vert \boldsymbol{r}-\boldsymbol{r}' \vert}\, d\boldsymbol{r} d\boldsymbol{r}'   \\

The first and second terms on the right hand side of the last equation are equal, since r and r′ in the second term can be interchanged without changing the value of the integral. Therefore,

 \int \frac{\delta J}{\delta\rho(\boldsymbol{r})} \phi(\boldsymbol{r})d\boldsymbol{r} = \int \left ( \int \frac {\rho(\boldsymbol{r}') }{\vert \boldsymbol{r}-\boldsymbol{r}' \vert} d\boldsymbol{r}' \right )  \phi(\boldsymbol{r}) d\boldsymbol{r}

and the functional derivative of the electron-electron coulomb potential energy functional J[ρ] is,[11]

 \frac{\delta J}{\delta\rho(\boldsymbol{r})} = \int \frac {\rho(\boldsymbol{r}') }{\vert \boldsymbol{r}-\boldsymbol{r}' \vert} d\boldsymbol{r}' \, .

The second functional derivative is

\frac{\delta^2 J[\rho]}{\delta \rho(\mathbf{r}')\delta\rho(\mathbf{r})}  = \frac{\partial}{\partial \rho(\mathbf{r}')} \left ( \frac{\rho(\mathbf{r}')}{\vert \mathbf{r}-\mathbf{r}' \vert} \right ) = \frac{1}{\vert \mathbf{r}-\mathbf{r}' \vert}.

Weizsäcker kinetic energy functional[edit]

In 1935 von Weizsäcker proposed to add a gradient correction to the Thomas-Fermi kinetic energy functional to make it suit better a molecular electron cloud:

T_\mathrm{W}[\rho] = \frac{1}{8} \int \frac{\nabla\rho(\mathbf{r}) \cdot \nabla\rho(\mathbf{r})}{ \rho(\mathbf{r}) } d\mathbf{r} =  \int t_\mathrm{W} \ d\mathbf{r} \, ,


 t_\mathrm{W} \equiv  \frac{1}{8}  \frac{\nabla\rho \cdot \nabla\rho}{ \rho } \qquad \text{and} \ \ \rho = \rho(\boldsymbol{r}) \ .

Using a previously derived formula for the functional derivative,

\frac{\delta T_\mathrm{W}}{\delta \rho(\boldsymbol{r})} 
& = \frac{\partial t_\mathrm{W}}{\partial \rho} - \nabla\cdot\frac{\partial t_\mathrm{W}}{\partial \nabla \rho} \\
& = -\frac{1}{8}\frac{\nabla\rho \cdot \nabla\rho}{\rho^2} - \left ( \frac {1}{4} \frac {\nabla^2\rho} {\rho} -  \frac {1}{4} \frac {\nabla\rho \cdot \nabla\rho} {\rho^2} \right ) \qquad \text{where} \ \ \nabla^2 = \nabla \cdot \nabla \ ,

and the result is,[12]

 \frac{\delta T_\mathrm{W}}{\delta \rho(\boldsymbol{r})} =  \ \ \, \frac{1}{8}\frac{\nabla\rho \cdot \nabla\rho}{\rho^2} - \frac{1}{4}\frac{\nabla^2\rho}{\rho} \ .


The entropy of a discrete random variable is a functional of the probability mass function.

H[p(x)] = -\sum_x p(x) \log p(x)


\sum_x \frac{\delta H}{\delta p(x)} \, \phi(x) 
& {} = \left[ \frac{d}{d\epsilon} H[p(x) + \epsilon\phi(x)] \right]_{\epsilon=0}\\
& {} = \left [- \, \frac{d}{d\varepsilon}  \sum_x \, [p(x) + \varepsilon\phi(x)] \ \log [p(x) + \varepsilon\phi(x)] \right]_{\varepsilon=0} \\
& {} = \displaystyle -\sum_x \, [1+\log p(x)] \ \phi(x) \, .


\frac{\delta H}{\delta p(x)} = -1-\log p(x).



 F[\varphi(x)]= e^{\int \varphi(x) g(x)dx}.

Using the delta function as a test function,

\frac{\delta F[\varphi(x)]}{\delta \varphi(y)} 
& {} = \lim_{\varepsilon\to 0}\frac{F[\varphi(x)+\varepsilon\delta(x-y)]-F[\varphi(x)]}{\varepsilon}\\
& {} = \lim_{\varepsilon\to 0}\frac{e^{\int (\varphi(x)+\varepsilon\delta(x-y)) g(x)dx}-e^{\int \varphi(x) g(x)dx}}{\varepsilon}\\
& {} = e^{\int \varphi(x) g(x)dx}\lim_{\varepsilon\to 0}\frac{e^{\varepsilon \int \delta(x-y) g(x)dx}-1}{\varepsilon}\\
& {} = e^{\int \varphi(x) g(x)dx}\lim_{\varepsilon\to 0}\frac{e^{\varepsilon g(y)}-1}{\varepsilon}\\
& {} = e^{\int \varphi(x) g(x)dx}g(y).


 \frac{\delta F[\varphi(x)]}{\delta \varphi(y)} = g(y) F[\varphi(x)].

This is particularly useful in calculating the correlation functions from the partition function in quantum field theory.

Functional derivative of a function[edit]

A function can be written in the form of an integral like a functional. For example,

\rho(\boldsymbol{r}) = F[\rho] = \int \rho(\boldsymbol{r}') \delta(\boldsymbol{r}-\boldsymbol{r}')\, d\boldsymbol{r}'.

Since the integrand does not depend on derivatives of ρ, the functional derivative of ρ(r) is,

\frac {\delta \rho(\boldsymbol{r})} {\delta\rho(\boldsymbol{r}')} \equiv \frac {\delta F} {\delta\rho(\boldsymbol{r}')} 
& = \frac{\partial \ \ }{\partial \rho(\boldsymbol{r}')} \, [\rho(\boldsymbol{r}') \delta(\boldsymbol{r}-\boldsymbol{r}')] \\
& = \delta(\boldsymbol{r}-\boldsymbol{r}').

Using the delta function as a test function[edit]

In physics, it's common to use the Dirac delta function \delta(x-y) in place of a generic test function \phi(x), for yielding the functional derivative at the point y (this is a point of the whole functional derivative as a partial derivative is a component of the gradient):

\frac{\delta F[\rho(x)]}{\delta \rho(y)}=\lim_{\varepsilon\to 0}\frac{F[\rho(x)+\varepsilon\delta(x-y)]-F[\rho(x)]}{\varepsilon}.

This works in cases when F[\rho(x)+\varepsilon f(x)] formally can be expanded as a series (or at least up to first order) in \varepsilon. The formula is however not mathematically rigorous, since F[\rho(x)+\varepsilon\delta(x-y)] is usually not even defined.

The definition given in a previous section is based on a relationship that holds for all test functions ϕ, so one might think that it should hold also when ϕ is chosen to be a specific function such as the delta function. However, the latter is not a valid test function.

In the definition, the functional derivative describes how the functional F[\varphi(x)] changes as a result of a small change in the entire function \varphi(x). The particular form of the change in \varphi(x) is not specified, but it should stretch over the whole interval on which x is defined. Employing the particular form of the perturbation given by the delta function has the meaning that \varphi(x) is varied only in the point y. Except for this point, there is no variation in \varphi(x).


  1. ^ According to Giaquinta & Hildebrandt (1996, p. 18), this notation is customary in physical literature.
  2. ^ Called differential in (Parr & Yang 1989, p. 246), variation or first variation in (Courant & Hilbert 1953, p. 186), and variation or differential in (Gelfand & Fomin 2000, p. 11, § 3.2).
  3. ^ For a three-dimensional cartesian coordinate system,
\frac{\partial f}{\partial\nabla\rho} =   \frac{\partial f}{\partial\rho_x}  \mathbf{\hat{i}}   + \frac{\partial f}{\partial\rho_y}  \mathbf{\hat{j}} + \frac{\partial f}{\partial\rho_z}  \mathbf{\hat{k}}\, ,  \qquad 
& \text{where} \ \rho_x = \frac{\partial \rho}{\partial x}\, ,  \ \rho_y = \frac{\partial \rho}{\partial y}\, , \ \rho_z = \frac{\partial \rho}{\partial z}\, \\
& \text{and} \ \ \mathbf{\hat{i}}, \ \mathbf{\hat{j}}, \ \mathbf{\hat{k}} \ \ \text {are unit vectors along the x, y, z axes.}
  4. ^ For example, for the case of three dimensions (n = 3) and second order derivatives (i = 2), the tensor (2) has components,
     \left [ \nabla^{(2)} \right ]_{\alpha \beta} = \frac {\partial^{\,2}} {\partial r_{\alpha} \, \partial r_{\beta}} \qquad \qquad \text{where} \quad \alpha, \beta = 1, 2, 3 \, .
  5. ^ For example, for the case n = 3 and i = 2, the tensor scalar product is,
     \nabla^{(2)} \cdot \frac{\partial f}{\partial\left(\nabla^{(2)}\rho\right)} = \sum_{\alpha, \beta = 1}^3 \ \frac {\partial^{\, 2} } {\partial r_{\alpha} \, \partial r_{\beta} }  \ \frac {\partial f} {\partial \rho_{\alpha \beta} }    \qquad \text{where} \ \ \rho_{\alpha \beta} \equiv \frac {\partial^{\, 2}\rho} {\partial r_{\alpha} \, \partial r_{\beta} } \ .


  1. ^ a b (Giaquinta & Hildebrandt 1996, p. 18)
  2. ^ a b (Parr & Yang 1989, p. 246).
  3. ^ (Gelfand & Fomin 2000, p. 28)
  4. ^ (Parr & Yang 1989, p. 246, Eq. A.2).
  5. ^ (Parr & Yang 1989, p. 246, Eq. A.1).
  6. ^ (Parr & Yang 1989, p. 247, Eq. A.3).
  7. ^ (Parr & Yang 1989, p. 247, Eq. A.4).
  8. ^ (Greiner & Reinhardt 1996, p. 38, Eq. 7).
  9. ^ (Parr & Yang 1989, p. 251, Eq. A.34).
  10. ^ (Parr & Yang 1989, p. 247, Eq. A.6).
  11. ^ (Parr & Yang 1989, p. 248, Eq. A.11).
  12. ^ (Parr & Yang 1989, p. 247, Eq. A.9).