# 19.7: Appendix - Aspects of Multivariate Calculus

- Page ID
- 32462

Multivariate calculus provides the framework for handling systems having many variables associated with each of several bodies. It is assumed that the reader has studied linear differential equations plus multivariate calculus and thus has been exposed to the calculus used in classical mechanics. Chapter \(5\) of this book introduced variational calculus which covers several important aspects of multivariate calculus such as Euler’s variational calculus and Lagrange multipliers. This appendix provides a brief review of a selection of other aspects of multivariate calculus that feature prominently in classical mechanics.

## Partial Differentiation

The extension of the derivative to multivariate calculus involves use of partial derivatives. The partial derivative with respect to the variable \(x_i\) of a multivariate function \(f(x_1, x_2,...., x_N )\) involves taking the normal one-variable derivative with respect to \(x_i\) assuming that the other \(N − 1\) variables are held constant. That is,

\[ \dfrac{\partial f\left(x_{1}, x_{2}, \ldots x_{N}\right)}{\partial x_{i}}=\lim _{h_{i} \rightarrow 0}\left[\dfrac{f\left(x_{1}, x_{2}, \ldots x_{i-1},\left(x_{i}+h_{i}\right), \ldots x_{N}\right)-f\left(x_{1}, x_{2}, \ldots, x_{N}\right)}{h_{i}}\right] \label{F.1}\]

where it will be assumed that the function \(f(x)\) is a continuously-differentiable function to \(n^{th}\) order, then all partial derivatives of that order or less are independent of the order in which they are performed. That is,

\[\dfrac{\partial^2 f(x)}{\partial x_i \partial x_j} = \dfrac{\partial^2 f(x)}{\partial x_j \partial x_i} \label{F.2}\]

The chain rule for partial differentiation gives that

\[\dfrac{\partial f (y_1, y_2, ...., y_N )}{ \partial y_j} = \sum^N_{k=1} \dfrac{\partial f(x)}{ \partial x_k} \dfrac{\partial x_k (y)}{\partial y_j} \label{F.3}\]

The total differential of a multivariate function \(f(x)\) is

\[df = \sum^N_{k=1} \dfrac{\partial f(x)}{ \partial x_k} dx_k \label{F.4}\]

This can be extended to higher-order derivatives using the operator formalism

\[ d^{n} f(x)=\left(d x_{1} \dfrac{\partial}{\partial x_{1}}+\ldots+d x_{N} \dfrac{\partial}{\partial x_{N}}\right)^{n} f(x)=\sum d x_{j_{1}} \ldots d x_{j_{n}} \dfrac{\partial^{n} f(x)}{\partial x_{j_{1}} \ldots \partial x_{j_{n}}} \label{F.5}\]

## Linear Operators

The linear operator notation provides a powerful, elegant, and compact way to express, and apply, the equations of multivariate calculus; it is used extensively in mathematics and physics. The linear operators typically comprise partial derivatives that act on scalar, vector, or tensor fields. Table \(\PageIndex{1}\) lists a few elementary examples of the use of linear operators in this textbook. The first four linear operators involve the widely used del operator \(\boldsymbol{\nabla}\) to generate the gradient, divergence and curl as described in appendices \(19.7\) and \(19.8\). The fifth and sixth linear operators act on the Lagrangian in Lagrangian mechanics applications. The final two linear operators act on the wavefunction for wave mechanics.

Name | Partial derivative | Field | Action |
---|---|---|---|

Gradient | \( \boldsymbol{\nabla} \equiv \hat{i} \dfrac{\partial}{ \partial x} + \hat{j} \dfrac{\partial}{ \partial y} + \mathbf{\hat{k}} \dfrac{\partial}{ \partial z}\) | Scalar potential \(V\) | \(\mathbf{E} = \boldsymbol{\nabla}V\) |

Divergence | \( \boldsymbol{\nabla} \cdot \equiv \left( \hat{i} \dfrac{\partial}{ \partial x} + \hat{j} \dfrac{\partial}{ \partial y} + \mathbf{\hat{k}} \dfrac{\partial}{ \partial z} \right) \cdot \) | Vector field \(\mathbf{E}\) | \(\boldsymbol{\nabla} \cdot \mathbf{E}\) |

Curl | \( \boldsymbol{\nabla} \times \equiv \left( \hat{i} \dfrac{\partial}{ \partial x} + \hat{j} \dfrac{\partial}{ \partial y} + \mathbf{\hat{k}} \dfrac{\partial}{ \partial z} \right) \times \) | Vector field \(\mathbf{E}\) | \(\boldsymbol{\nabla} \times \mathbf{E}\) |

Laplacian | \(\nabla^2 = \boldsymbol{\nabla} \cdot \boldsymbol{\nabla} \equiv \dfrac{\partial^2}{ \partial x^2} + \dfrac{\partial^2}{ \partial y^2} + \dfrac{\partial^2}{ \partial z^2}\) | Scalar potential \(V\) | \(\nabla^2V\) |

Euler-Lagrange | \(\Lambda_j \equiv \dfrac{d}{dt} \dfrac{\partial}{ \partial \dot{q}_j} − \dfrac{\partial}{ \partial q_j}\) | Scalar Lagrangian \(L\) | \(\Lambda L = 0\) |

Canonical momentum | \(p_j \equiv \dfrac{\partial}{ \partial \dot{q}_j}\) | Scalar Lagrangian \(L\) | \(p_j \equiv \dfrac{\partial L}{ \partial \dot{q}_j}\) |

Canonical momentum | \(p_j \equiv \dfrac{\hbar}{ i} \dfrac{\partial}{ \partial \dot{q}_j}\) | Wavefunction \(\Psi\) | \(p_j\Psi \equiv \dfrac{\hbar}{ i} \dfrac{\partial \Psi}{ \partial \dot{q}_j}\) |

Hamiltonian | \(H = i\hbar \dfrac{ \partial }{ \partial t}\) | Wavefunction \(\Psi\) | \(H\Psi = i\hbar \dfrac{ \partial \Psi}{ \partial t} = E\Psi\) |

There are three ways of expressing operations such as addition, multiplication, transposition or inversion of operations that are completely equivalent because they all are based on the same principles of linear algebra. For example, a transformation \(\mathbf{O}\) acting on a vector \(\mathbf{A}\) can produce the vector \(\mathbf{B}\). The simplest way to express this transformation is in terms of components

\[B_i = \sum^3_{j=1} O_{ij}A_j \label{F.6}\]

Another way is to use matrix mechanics where the \(3 \times 3\) matrix \((\mathbf{O})\) transforms the column vector \((\mathbf{A})\) to the column vector \((\mathbf{B})\), that is,

\[(\mathbf{B})=(\mathbf{O}) (\mathbf{A}) \label{F.7}\]

The third approach is to assume an operator \(\mathbf{O}\) acts on the vector \(\mathbf{A}\)

\[\mathbf{B} = \mathbf{OA} \label{F.8}\]

In classical mechanics, and quantum mechanics, these three equivalent approaches are used and exploited extensively and interchangeably. In particular the rules of matrix manipulation, that are given in appendix \(19.1\), are synonymous, and equivalent to, those that apply for operator manipulation. If the operator is complex then the operator properties are summarized as follows.

The generalization of the transpose for complex operators is the *Hermitian conjugate* \(O^{\dagger}\)

\[O^{\dagger}_{ij} = O^*_{ji} \label{F.9}\]

Note also that

\[\mathbf{O}^{\dagger} = (O^*)^T = (O^T )^* \label{F.10}\]

The generalization of a symmetric matrix is *Hermitian*, that is, \(O\) is equal to its Hermitian conjugate

\[O^{\dagger}_{ij} = O^*_{ji} = O_{ij} \label{F.11}\]

For a real matrix the complex conjugation has no effect so the matrix is real and symmetric.

The generalization of orthogonal is *unitary* for which the operator is unitary if it is non-singular and

\[O^{−1} = O^{\dagger} \label{F.12}\]

which implies

\[OO^{\dagger} = U = O^{\dagger}O \label{F.13}\]

## Transformation Jacobian

The Jacobian determinant, which is usually called the Jacobian, is used extensively in mechanics for both rotational and translational coordinate transformations. The Jacobian determinant is defined as being the ratio of the \(n\)-dimensional volume element \(dx_1dx_2...dx_n\) in one coordinate system, to the volume element \(dy_1dy_2...dy_n\) in the second coordinate system. That is

\[J\left(y_{1} y_{2} \ldots y_{n}\right) \equiv \dfrac{\partial x_{1} \partial x_{2} \ldots \partial x_{n}}{\partial y_{1} \partial y_{2} \ldots \partial y_{n}}=\begin{vmatrix}

\dfrac{\partial x_{1}}{\partial y_{1}} & \dfrac{\partial x_{1}}{\partial y_{2}} & \ldots & \dfrac{\partial x_{1}}{\partial y_{n}} \\ \dfrac{\partial x_{2}}{\partial y_{1}} & \dfrac{\partial x_{2}}{\partial y_{2}} & \cdots & \dfrac{\partial x_{2}}{\partial y_{n}} \\ \vdots & \vdots & \vdots & \vdots \\ \dfrac{\partial x_{n}}{\partial y_{1}} & \dfrac{\partial x_{n}}{\partial y_{2}} & \ldots & \dfrac{\partial x_{n}}{\partial y_{n}} \end{vmatrix} \label{F.14}\]

### Transformation of integrals

Consider a coordinate transformation for the integral of the function \(f(x_1, x_2, ..x_n)\) to the integral of a function \(g(y_1, y_2, ...y_n)\) where \(y_i = h (x_1, x_2, ...x_n)\). The coordinate transformation of the integral equation can be expressed in terms of the Jacobian \(J(y_1y_2...y_n)\)

\[\begin{align} \label{F.15} \int f\left(x_{1}, x_{2}, \ldots x_{n}\right) d x_{1} d x_{2} \ldots d x_{n} &=\int g\left(y_{1}, y_{2}, \ldots y_{n}\right) d y_{1} d y_{2} \ldots d y_{n}=\\ \int f\left(x_{1}, x_{2}, \ldots x_{n}\right) \dfrac{\partial x_{1} \partial x_{2} \ldots \partial x_{n}}{\partial y_{1} \partial y_{2} \ldots \partial y_{n}} d y_{1} d y_{2} \ldots d y_{n} &=\int f\left(y_{1}, y_{2}, . . y_{n}\right) J\left(y_{1}, y_{2}, \ldots y_{n}\right) d y_{1} d y_{2} \ldots d y_{n} \nonumber \end{align} \]

### Transformation of differential equations

The differential cross sections for scattering can be defined either by the number of a definite kind of particle/per event, going into the volume element in momentum space \(dp_1dp_2dp_3\), or by the number going into the solid angle element having momentum between \(p\) and \(p + dp\). That is, the first definition can be written as a differential equation

\[\dfrac{\partial^3S(p_1, p_2, p_3)}{ \partial p_1\partial p_2\partial p_3 } dp_1dp_2dp_3 = \dfrac{\partial^3 S (p_1(p\theta \phi ), p_2(p\theta \phi ), p_3(p\theta \phi )) }{\partial p_1\partial p_2\partial p_3 } \dfrac{\partial (p_1, p_2, p_3) }{\partial (p, \theta , \phi )} dpd\theta d\phi \label{F.16}\]

As shown in table \(19.3.4\), \(dp_1dp_2dp_3 = p^2 \sin \theta dpd\theta d\phi \), that is, the Jacobian equals \(p^2 \sin \theta \). Thus Equation \ref{F.16} can be written as

\[\dfrac{\partial^3S(p_1, p_2, p_3)}{ \partial p_1\partial p_2\partial p_3} dp_1dp_2dp_3 = \left[\dfrac{ \partial^3S }{\partial p_1\partial p_2\partial p_3} p^2 \right] (\sin \theta dpd\theta d\phi ) = \dfrac{\partial^2 \sigma (p, \theta , \phi )}{ \partial p\partial \Omega} dpd\Omega \label{F.17}\]

The differential cross section is defined by

\[\dfrac{\partial^2\sigma (p, \theta , \phi )}{ \partial p\partial \Omega} \equiv \dfrac{\partial^3S}{ \partial p_1\partial p_2\partial p_3} p^2 \label{F.18}\]

where the \(p^2\) factor is absorbed into the cross section and the solid angle term is factored out

### Properties of the Jacobian

In classical mechanics the Jacobian often is extended from 3 dimensions to \(n\)-dimensional transformations. The Jacobian is unity for unitary transformations such as rotations and linear translations which implies that the volume element is preserved. It will be shown that this also is true for a certain class of transformations in classical mechanics that are called canonical transformations. The Jacobian transforms the local density to be correct for any scale transformations such as transforming linear dimensions from centimeters to inches.

Example \(\PageIndex{1}\): Jacobian for transform from cartesian to spherical coordinates

Consider the transform in the three-dimensional integral \(\int (x_1, x_2, x_3)dx_1dx_2dx_3\) under transformation from cartesian coordinates \((x_1, x_2, x_3)\) to spherical coordinates \((r, \theta , \phi )\). The transformation is governed by the geometric relations \(x_1 = r \sin \theta \cos \phi , x_2 = r \sin \theta \sin \phi , x_3 = r \cos \theta \). For this transformation the Jacobian determinant equals

\[J(r, \theta, \phi)= \begin{vmatrix} \sin \theta \cos \phi & r \cos \theta \cos \phi & -r \sin \theta \sin \phi \\ \sin \theta \sin \phi & r \cos \theta \sin \phi & r \sin \theta \cos \phi \\ \cos \theta & -r \sin \theta & 0 \end{vmatrix} = r^{2} \sin \theta \nonumber\]

Thus the three-dimensional volume integral transforms to

\[\int f(x_1, x_2, x_3)dx_1dx_2dx_3 = \int f(r, \theta , \phi ) J (r, \theta , \phi ) drd\theta d\phi = \int f(r, \theta , \phi )r^2 \sin \theta drd\theta d\phi \nonumber\]

which is the well-known volume integral in spherical coordinates.

## Legendre transformation

Hamiltonian mechanics can be derived directly from Lagrange mechanics by considering the Legendre transformation between the conjugate variables \((\mathbf{q}, \mathbf{\dot{q}}, t)\) and \((\mathbf{q}, \mathbf{p}, t)\). Such a derivation is of considerable importance in that it shows that Hamiltonian mechanics is based on the same variational principles as those used to derive Lagrangian mechanics; that is d’Alembert’s Principle or Hamilton’s Principle. The general problem of converting Lagrange’s equations into the Hamiltonian form hinges on the inversion of equation \((8.1.3)\) that defines the generalized momentum \(\mathbf{p}\). This inversion is simplified by the fact that \((8.1.3)\) is the first partial derivative of the Lagrangian \(L(\mathbf{q}, \mathbf{\dot{q}}, t)\) which is a scalar function.

Consider transformations between two functions \(F(\mathbf{u}, \mathbf{w})\) and \(G(\mathbf{v}, \mathbf{w})\) where \(\mathbf{u}\) and \(\mathbf{v}\) are the active variables related by the functional form

\[\mathbf{v} = \boldsymbol{\nabla}_{\mathbf{u}} F(\mathbf{u}, \mathbf{w}) \label{F.19}\]

and where \(\mathbf{w}\) designates passive variables and \(\boldsymbol{\nabla}_{\mathbf{u}}F(\mathbf{u}, \mathbf{w})\) is the first-order derivative of \(F(\mathbf{u}, \mathbf{w})\), i.e. the gradient, with respect to the components of the vector \(\mathbf{u}\). The Legendre transform states that the inverse formula can always be written in the form

\[\mathbf{u} = \boldsymbol{\nabla}_{\mathbf{v}}G(\mathbf{v}, \mathbf{w}) \label{F.20}\]

where the function \(G(\mathbf{v}, \mathbf{w})\) is related to \(F(\mathbf{u}, \mathbf{w})\) by the symmetric relation

\[G(\mathbf{v}, \mathbf{w}) + F(\mathbf{u}, \mathbf{w}) = \mathbf{u} \cdot \mathbf{v} \label{F.21}\]

and where the scalar product \(\mathbf{u} \cdot \mathbf{ v} = \sum^N_{i = 1} u_iv_i\).

Furthermore the derivatives with respect to all the passive variables \(\{w_i\}\) are related by

\[\boldsymbol{\nabla}_{\mathbf{w}}F(\mathbf{u}, \mathbf{w}) = −\boldsymbol{\nabla}_{\mathbf{w}} G(\mathbf{v}, \mathbf{w}) \label{F.22}\]

The relationship between the functions \(F(\mathbf{u}, \mathbf{w})\) and \(G(\mathbf{v}, \mathbf{w})\) is symmetrical and each is said to be the Legendre transform of the other.

## Exercises

1. Below you will find a set of integrals. Your teaching assistant will divide you into groups and each group will be assigned one integral to work on. Once your group has solved the integral, write the solution on the board in the space provided by the teaching assistant.

(a) \(\int^{2\pi}_0 \int^{\pi/4}_{0} \int^{\cos \theta}_0 r^2 \sin \theta dr d\theta d \phi\)

(b) \(\int (\dfrac{\mathbf{\dot{r}}}{r} - \dfrac{\mathbf{r}\dot{r}}{r^2}) dt\)

(c) \(\int_S \mathbf{A} \cdot d\mathbf{a}\) where \(\mathbf{A} = x\hat{i} + y\hat{j} + z\mathbf{\hat{k}}\) and \(S\) is the sphere \(x^2 + y^2 + z^2 = 9\).

(d) \(\int_S ( \boldsymbol{\nabla} \times \mathbf{A}) \cdot d\mathbf{a}\) where \(\mathbf{A} = y\hat{i} + z\hat{j} + x\mathbf{\hat{k}}\) and \(S\) is the surface defined by the paraboloid \(z = 1−x^2 − y^2\), where \(z \geq 0\).