Volume element

In mathematics, a volume element provides a means for integrating a function with respect to volume in various coordinate systems such as spherical coordinates and cylindrical coordinates. Thus a volume element is an expression of the form

dV = \rho(u_1,u_2,u_3)\,du_1\,du_2\,du_3

where the $u_i$ are the coordinates, so that the volume of any set $B$ can be computed by

\operatorname{Volume}(B) = \int_B \rho(u_1,u_2,u_3)\,du_1\,du_2\,du_3.

For example, in spherical coordinates $dV = u_1^2\sin u_2\,du_1\,du_2\,du_3$ , and so $\rho = u_1^2\sin u_2$ .

The notion of a volume element is not limited to three dimensions: in two dimensions it is often known as the area element, and in this setting it is useful for doing surface integrals. Under changes of coordinates, the volume element changes by the absolute value of the Jacobian determinant of the coordinate transformation (by the change of variables formula). This fact allows volume elements to be defined as a kind of measure on a manifold. On an orientable differentiable manifold, a volume element typically arises from a volume form: a top degree differential form. On a non-orientable manifold, the volume element is typically the absolute value of a (locally defined) volume form: it defines a 1-density.

Volume element in Euclidean space

In Euclidean space, the volume element is given by the product of the differentials of the Cartesian coordinates

dV = dx\,dy\,dz.

In different coordinate systems of the form $x=x(u_1,u_2,u_3), y=y(u_1,u_2,u_3), z=z(u_1,u_2,u_3)$ , the volume element changes by the Jacobian of the coordinate change:

dV = \left|\frac{\partial (x,y,z)}{\partial (u_1,u_2,u_3)}\right|\,du_1\,du_2\,du_3.

For example, in spherical coordinates

\begin{align} x&=\rho\cos\theta\sin\phi\\ y&=\rho\sin\theta\sin\phi\\ z&=\rho\cos\phi \end{align}

the Jacobian is

\left |\frac{\partial(x,y,z)}{\partial (\rho,\theta,\phi)}\right| = \rho^2\sin\phi

so that

dV = \rho^2\sin\phi\,d\rho\,d\theta\,d\phi.

This can be seen as a special case of the fact that differential forms transform through a pullback $F^*$ as

F^*(u \; dy^1 \wedge \cdots \wedge dy^n) = (u \circ F) \det \left(\frac{\partial F^j}{\partial x^i}\right) dx^1 \wedge \cdots \wedge dx^n

Volume element of a linear subspace

Consider the linear subspace of the n-dimensional Euclidean space Rⁿ that is spanned by a collection of linearly independent vectors

X_1,\dots,X_k.

To find the volume element of the subspace, it is useful to know the fact from linear algebra that the volume of the parallelepiped spanned by the $X_i$ is the square root of the determinant of the Gramian matrix of the $X_i$ :

\sqrt{\det(X_i\cdot X_j)_{i,j=1\dots k}}.

Any point p in the subspace can be given coordinates $(u_1,u_2,\dots,u_k)$ such that

p = u_1X_1 + \cdots + u_kX_k.

At a point p, if we form a small parallelepiped with sides $du_i$ , then the volume of that parallelepiped is the square root of the determinant of the Grammian matrix

\sqrt{\det\left((du_i X_i)\cdot (du_j X_j)\right)_{i,j=1\dots k}} = \sqrt{\det(X_i\cdot X_j)_{i,j=1\dots k}}\; du_1\,du_2\,\cdots\,du_k.

This therefore defines the volume form in the linear subspace.

Volume element of manifolds

On an oriented Riemannian manifold of dimension n, the volume element is a volume form equal to the Hodge dual of the unit constant function, $f(x) = 1$ :

\omega = \star 1

Equivalently, the volume element is precisely the Levi-Civita tensor $\epsilon$ .^[1] In coordinates,

\omega = \sqrt{|\det g|}\, dx^1 \wedge \cdots \wedge dx^n

where $\det g$ is the determinant of the metric tensor g written in the coordinate system.

Area element of a surface

A simple example of a volume element can be explored by considering a two-dimensional surface embedded in n-dimensional Euclidean space. Such a volume element is sometimes called an area element. Consider a subset $U \subset \mathbf{R}^2$ and a mapping function

\varphi:U\to \mathbf{R}^n

thus defining a surface embedded in $\mathbf{R}^n$ . In two dimensions, volume is just area, and a volume element gives a way to determine the area of parts of the surface. Thus a volume element is an expression of the form

f(u_1,u_2)\,du_1\,du_2

that allows one to compute the area of a set B lying on the surface by computing the integral

\operatorname{Area}(B) = \int_B f(u_1,u_2)\,du_1\,du_2.

Here we will find the volume element on the surface that defines area in the usual sense. The Jacobian matrix of the mapping is

\lambda_{ij}=\frac{\partial \varphi_i} {\partial u_j}

with index i running from 1 to n, and j running from 1 to 2. The Euclidean metric in the n-dimensional space induces a metric $g=\lambda^T\lambda$ on the set U, with matrix elements

g_{ij}=\sum_{k=1}^n \lambda_{ki} \lambda_{kj} = \sum_{k=1}^n \frac{\partial \varphi_k} {\partial u_i} \frac{\partial \varphi_k} {\partial u_j}.

The determinant of the metric is given by

\det g = \left| \frac{\partial \varphi} {\partial u_1} \wedge \frac{\partial \varphi} {\partial u_2} \right|^2 = \det (\lambda^T \lambda)

For a regular surface, this determinant is non-vanishing; equivalently, the Jacobian matrix has rank 2.

Now consider a change of coordinates on U, given by a diffeomorphism

f \colon U\to U , \,\!

so that the coordinates $(u_1,u_2)$ are given in terms of $(v_1,v_2)$ by $(u_1,u_2)= f(v_1,v_2)$ . The Jacobian matrix of this transformation is given by

F_{ij}= \frac{\partial f_i} {\partial v_j}.

In the new coordinates, we have

\frac{\partial \varphi_i} {\partial v_j} = \sum_{k=1}^2 \frac{\partial \varphi_i} {\partial u_k} \frac{\partial f_k} {\partial v_j}

and so the metric transforms as

\tilde{g} = F^T g F

where $\tilde{g}$ is the pullback metric in the v coordinate system. The determinant is

\det \tilde{g} = \det g (\det F)^2.

Given the above construction, it should now be straightforward to understand how the volume element is invariant under an orientation-preserving change of coordinates.

In two dimensions, the volume is just the area. The area of a subset $B\subset U$ is given by the integral

\begin{align} \mbox{Area}(B) &= \iint_B \sqrt{\det g}\; du_1\; du_2 \\ &= \iint_B \sqrt{\det g} \;|\det F| \;dv_1 \;dv_2 \\ &= \iint_B \sqrt{\det \tilde{g}} \;dv_1 \;dv_2. \end{align}

Thus, in either coordinate system, the volume element takes the same expression: the expression of the volume element is invariant under a change of coordinates.

Note that there was nothing particular to two dimensions in the above presentation; the above trivially generalizes to arbitrary dimensions.

Example: Sphere

For example, consider the sphere with radius r centered at the origin in R³. This can be parametrized using spherical coordinates with the map

\phi(u_1,u_2) = (r\cos u_1\sin u_2,r\sin u_1\sin u_2,r\cos u_2).

Then

g = \begin{pmatrix}r^2\sin^2u_2 & 0 \\ 0 & r^2\end{pmatrix},

and the area element is

\omega = \sqrt{\det g}\; du_1 du_2 = r^2\sin u_2\, du_1 du_2.

References

Besse, Arthur L. (1987), Einstein manifolds, Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)], vol. 10, Berlin, New York: Springer-Verlag, pp. xii+510, ISBN 978-3-540-15279-8

↑ Carroll, Sean. Spacetime and Geometry. Addison Wesley, 2004, p. 90