Spherical coordinate system

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a given point in space is specified by three numbers, (r, θ, φ): the radial distance of the radial line r connecting the point to the fixed point of origin (which is located on a fixed polar axis, or zenith direction axis, or z-axis); the polar angle θ of the radial line r; and the azimuthal angle φ of the radial line r.

The polar angle θ is measured between the z-axis and the radial line r. The azimuthal angle φ is measured between the

Once the radius is fixed, the three coordinates (r, θ, φ), known as a 3-tuple, provide a coordinate system on a sphere, typically called the spherical polar coordinates. Nota bene: the physics convention is followed in this article; (See both graphics re "physics convention" and re "mathematics convention").

The radial distance from the fixed point of origin is also called the radius, or radial line, or radial coordinate. The polar angle may be called

Both the use of symbols and the naming order of tuple coordinates differ among the several sources and disciplines. This article will use the ISO convention[1] frequently encountered in physics, where the naming tuple gives the order as: radial distance, polar angle, azimuthal angle, or ${\displaystyle (r,\theta ,\varphi )}$. (See graphic re the "physics convention".) In contrast, the conventions in many mathematics books and texts give the naming order differently as: radial distance, "azimuthal angle", "polar angle", and ${\displaystyle (\rho ,\theta ,\varphi )}$ or ${\displaystyle (r,\theta ,\varphi )}$—which switches the uses and meanings of symbols θ and φ. Other conventions may also be used, such as r for a radius from the z-axis that is not from the point of origin. Particular care must be taken to check the meaning of the symbols.

According to the conventions of

The spherical coordinate system of the physics convention can be seen as a generalization of the polar coordinate system in three-dimensional space. It can be further extended to higher-dimensional spaces, and is then referred to as a hyperspherical coordinate system.

Definition

To define a spherical coordinate system, one must designate an origin point in space, O, and two orthogonal directions: the zenith reference direction and the azimuth reference direction. These choices determine a reference plane that is typically defined as containing the point of origin and the x– and y–axes, either of which may be designated as the azimuth reference direction. The reference plane is perpendicular (orthogonal) to the zenith direction, and typically is desiginated "horizontal" to the zenith direction's "vertical". The spherical coordinates of a point P then are defined as follows:

• The radius or radial distance is the Euclidean distance from the origin O to P.
• The inclination (or polar angle) is the signed angle from the zenith reference direction to the line segment OP. (Elevation may be used as the polar angle instead of inclination; see below.)
• The azimuth (or azimuthal angle) is the signed angle measured from the azimuth reference direction to the orthogonal projection of the radial line segment OP on the reference plane.

The sign of the azimuth is determined by designating the rotation that is the positive sense of turning about the zenith. This choice is arbitrary, and is part of the coordinate system definition. (If the inclination is either zero or 180 degrees (= π radians), the azimuth is arbitrary. If the radius is zero, both azimuth and inclination are arbitrary.)

The elevation is the signed angle from the x-y reference plane to the radial line segment OP, where positive angles are designated as upward, towards the zenith reference. Elevation is 90 degrees (= π/2 radians) minus inclination. Thus, if the inclination is 60 degrees (= π/3 radians), then the elevation is 30 degrees (= π/6 radians).

In

Conventions

Several different conventions exist for representing spherical coordinates and prescribing the naming order of their symbols. The 3-tuple number set ${\displaystyle (r,\theta ,\varphi )}$ denotes radial distance, the polar angle—"inclination", or as the alternative, "elevation"—and the azimuthal angle. It is the common practice within the physics convention, as specified by

However, some authors (including mathematicians) use the symbol ρ (rho) for radius, or radial distance, φ for inclination (or elevation) and θ for azimuth—while others keep the use of r for the radius; all which "provides a logical extension of the usual polar coordinates notation".[3] As to order, some authors list the azimuth before the inclination (or the elevation) angle. Some combinations of these choices result in a left-handed coordinate system. The standard "physics convention" 3-tuple set ${\displaystyle (r,\theta ,\varphi )}$ conflicts with the usual notation for two-dimensional polar coordinates and three-dimensional cylindrical coordinates, where θ is often used for the azimuth.^[3]

Angles are typically measured in degrees (°) or in radians (rad), where 360° = 2π rad. The use of degrees is most common in geography, astronomy, and engineering, where radians are commonly used in mathematics and theoretical physics. The unit for radial distance is usually determined by the context, as occurs in applications of the 'unit sphere', see applications.

When the system is used to designate physical three-space, it is customary to assign positive to azimuth angles measured in the counterclockwise sense from the reference direction on the reference plane—as seen from the "zenith" side of the plane. This convention is used in particular for geographical coordinates, where the "zenith" direction is north and the positive azimuth (longitude) angles are measured eastwards from some prime meridian.

Note:

Unique coordinates

Any spherical coordinate triplet (or tuple) ${\displaystyle (r,\theta ,\varphi )}$ specifies a single point of three-dimensional space. On the reverse view, any single point has infinitely many equivalent spherical coordinates. That is, the user can add or subtract any number of full turns to the angular measures without changing the angles themselves, and therefore without changing the point. It is convenient in many contexts to use negative radial distances, the convention being ${\displaystyle (-r,-\theta ,\varphi )}$, which is equivalent to ${\displaystyle (r,\theta ,\varphi )}$ for any r, θ, and φ. Moreover, ${\displaystyle (r,-\theta ,\varphi )}$ is equivalent to ${\displaystyle (r,\theta ,\varphi {+}180^{\circ })}$.

When necessary to define a unique set of spherical coordinates for each point, the user must restrict the range, aka interval, of each coordinate. A common choice is:

But instead of the interval [0°, 360°), the azimuth φ is typically restricted to the half-open interval (−180°, +180°], or (−π, +π ] radians, which is the standard convention for geographic longitude.

For the polar angle θ, the range (interval) for inclination is [0°, 180°], which is equivalent to elevation range (interval) [−90°, +90°]. In geography, the latitude is the elevation.

Even with these restrictions, if the polar angle (inclination) is 0° or 180°—elevation is −90° or +90°—then the azimuth angle is arbitrary; and if r is zero, both azimuth and polar angles are arbitrary. To define the coordinates as unique, the user can assert the convention that (in these cases) the arbitrary coordinates are set to zero.

Plotting

To plot any dot from its spherical coordinates (r, θ, φ), where θ is inclination, the user would: move r units from the origin in the zenith reference direction (z-axis); then rotate by the amount of the azimuth angle (φ) about the origin from the designated azimuth reference direction, (i.e., either the x– or y–axis, see Definition, above); and then rotate from the z-axis by the amount of the θ angle.

Applications

Just as the two-dimensional Cartesian coordinate system is useful—has a wide set of applications—on a planar surface, a two-dimensional spherical coordinate system is useful on the surface of a sphere. For example, one sphere that is described in Cartesian coordinates with the equation x² + y² + z² = c² can be described in spherical coordinates by the simple equation r = c. (In this system—shown here in the mathematics convention—the sphere is adapted as a unit sphere, where the radius is set to unity and then can generally be ignored, see graphic.)

This (unit sphere) simplification is also useful when dealing with objects such as rotational matrices. Spherical coordinates are also useful in analyzing systems that have some degree of symmetry about a point, including: volume integrals inside a sphere; the potential energy field surrounding a concentrated mass or charge; or global weather simulation in a planet's atmosphere.

Three dimensional modeling of loudspeaker output patterns can be used to predict their performance. A number of polar plots are required, taken at a wide selection of frequencies, as the pattern changes greatly with frequency. Polar plots help to show that many loudspeakers tend toward omnidirectionality at lower frequencies.

An important application of spherical coordinates provides for the

Instead of inclination, the

The azimuth angle (or

.

Instead of the radial distance r geographers commonly use

radius of Earth

, which is approximately 6,360 ± 11 km (3,952 ± 7 miles).

However, modern geographical coordinate systems are quite complex, and the positions implied by these simple formulae may be inaccurate by several kilometers. The precise standard meanings of latitude, longitude and altitude are currently defined by the World Geodetic System (WGS), and take into account the flattening of the Earth at the poles (about 21 km or 13 miles) and many other details.

Planetary coordinate systems use formulations analogous to the geographic coordinate system.

In astronomy

A series of astronomical coordinate systems are used to measure the elevation angle from several fundamental planes. These reference planes include: the observer's horizon, the galactic equator (defined by the rotation of the Milky Way), the celestial equator (defined by Earth's rotation), the plane of the ecliptic (defined by Earth's orbit around the Sun), and the plane of the earth terminator (normal to the instantaneous direction to the Sun).

Coordinate system conversions

See also: List of common coordinate transformations § To spherical coordinates

As the spherical coordinate system is only one of many three-dimensional coordinate systems, there exist equations for converting coordinates between the spherical coordinate system and others.

Cartesian coordinates

The spherical coordinates of a point in the ISO convention (i.e. for physics: radius r, inclination θ, azimuth φ) can be obtained from its Cartesian coordinates (x, y, z) by the formulae

$${\displaystyle {\begin{aligned}r&={\sqrt {x^{2}+y^{2}+z^{2}}}\\\theta &=\arccos {\frac {z}{\sqrt {x^{2}+y^{2}+z^{2}}}}=\arccos {\frac {z}{r}}={\begin{cases}\arctan {\frac {\sqrt {x^{2}+y^{2}}}{z}}&{\text{if }}z>0\\\pi +\arctan {\frac {\sqrt {x^{2}+y^{2}}}{z}}&{\text{if }}z<0\\+{\frac {\pi }{2}}&{\text{if }}z=0{\text{ and }}xy\neq 0\\{\text{undefined}}&{\text{if }}x=y=z=0\\\end{cases}}\\\varphi &=\operatorname {sgn}(y)\arccos {\frac {x}{\sqrt {x^{2}+y^{2}}}}={\begin{cases}\arctan({\frac {y}{x}})&{\text{if }}x>0,\\\arctan({\frac {y}{x}})+\pi &{\text{if }}x<0{\text{ and }}y\geq 0,\\\arctan({\frac {y}{x}})-\pi &{\text{if }}x<0{\text{ and }}y<0,\\+{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y>0,\\-{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y<0,\\{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}\end{aligned}}}$$

The

inverse tangent denoted in φ = arctan y/x must be suitably defined, taking into account the correct quadrant of (x, y). See the article on atan2

.

Alternatively, the conversion can be considered as two sequential rectangular to polar conversions: the first in the Cartesian xy plane from (x, y) to (R, φ), where R is the projection of r onto the xy-plane, and the second in the Cartesian zR-plane from (z, R) to (r, θ). The correct quadrants for φ and θ are implied by the correctness of the planar rectangular to polar conversions.

These formulae assume that the two systems have the same origin, that the spherical reference plane is the Cartesian xy plane, that θ is inclination from the z direction, and that the azimuth angles are measured from the Cartesian x axis (so that the y axis has φ = +90°). If θ measures elevation from the reference plane instead of inclination from the zenith the arccos above becomes an arcsin, and the cos θ and sin θ below become switched.

Conversely, the Cartesian coordinates may be retrieved from the spherical coordinates (radius r, inclination θ, azimuth φ), where r ∈ [0, ∞), θ ∈ [0, π], φ ∈ [0, 2π), by

$${\displaystyle {\begin{aligned}x&=r\sin \theta \,\cos \varphi ,\\y&=r\sin \theta \,\sin \varphi ,\\z&=r\cos \theta .\end{aligned}}}$$

Cylindrical coordinates

Main article: Cylindrical coordinate system

Cylindrical coordinates (axial radius ρ, azimuth φ, elevation z) may be converted into spherical coordinates (central radius r, inclination θ, azimuth φ), by the formulas

$${\displaystyle {\begin{aligned}r&={\sqrt {\rho ^{2}+z^{2}}},\\\theta &=\arctan {\frac {\rho }{z}}=\arccos {\frac {z}{\sqrt {\rho ^{2}+z^{2}}}},\\\varphi &=\varphi .\end{aligned}}}$$

Conversely, the spherical coordinates may be converted into cylindrical coordinates by the formulae

$${\displaystyle {\begin{aligned}\rho &=r\sin \theta ,\\\varphi &=\varphi ,\\z&=r\cos \theta .\end{aligned}}}$$

These formulae assume that the two systems have the same origin and same reference plane, measure the azimuth angle φ in the same senses from the same axis, and that the spherical angle θ is inclination from the cylindrical z axis.

Generalization

See also: Ellipsoidal coordinates

It is also possible to deal with ellipsoids in Cartesian coordinates by using a modified version of the spherical coordinates.

Let P be an ellipsoid specified by the level set

$${\displaystyle ax^{2}+by^{2}+cz^{2}=d.}$$

The modified spherical coordinates of a point in P in the ISO convention (i.e. for physics: radius r, inclination θ, azimuth φ) can be obtained from its Cartesian coordinates (x, y, z) by the formulae

$${\displaystyle {\begin{aligned}x&={\frac {1}{\sqrt {a}}}r\sin \theta \,\cos \varphi ,\\y&={\frac {1}{\sqrt {b}}}r\sin \theta \,\sin \varphi ,\\z&={\frac {1}{\sqrt {c}}}r\cos \theta ,\\r^{2}&=ax^{2}+by^{2}+cz^{2}.\end{aligned}}}$$

An infinitesimal volume element is given by

$${\displaystyle \mathrm {d} V=\left|{\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}\right|\,dr\,d\theta \,d\varphi ={\frac {1}{\sqrt {abc}}}r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi ={\frac {1}{\sqrt {abc}}}r^{2}\,\mathrm {d} r\,\mathrm {d} \Omega .}$$

The square-root factor comes from the property of the determinant that allows a constant to be pulled out from a column:

$${\displaystyle {\begin{vmatrix}ka&b&c\\kd&e&f\\kg&h&i\end{vmatrix}}=k{\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}}.}$$

Integration and differentiation in spherical coordinates

The following equations (Iyanaga 1977) assume that the colatitude θ is the inclination from the positive z axis, as in the physics convention discussed.

The line element for an infinitesimal displacement from (r, θ, φ) to (r + dr, θ + dθ, φ + dφ) is

$${\displaystyle \mathrm {d} \mathbf {r} =\mathrm {d} r\,{\hat {\mathbf {r} }}+r\,\mathrm {d} \theta \,{\hat {\boldsymbol {\theta }}}+r\sin {\theta }\,\mathrm {d} \varphi \,\mathbf {\hat {\boldsymbol {\varphi }}} ,}$$

where

$${\displaystyle {\begin{aligned}{\hat {\mathbf {r} }}&=\sin \theta \cos \varphi \,{\hat {\mathbf {x} }}+\sin \theta \sin \varphi \,{\hat {\mathbf {y} }}+\cos \theta \,{\hat {\mathbf {z} }},\\{\hat {\boldsymbol {\theta }}}&=\cos \theta \cos \varphi \,{\hat {\mathbf {x} }}+\cos \theta \sin \varphi \,{\hat {\mathbf {y} }}-\sin \theta \,{\hat {\mathbf {z} }},\\{\hat {\boldsymbol {\varphi }}}&=-\sin \varphi \,{\hat {\mathbf {x} }}+\cos \varphi \,{\hat {\mathbf {y} }}\end{aligned}}}$$

are the local orthogonal

unit vectors

in the directions of increasing r, θ, and φ, respectively, and x̂, ŷ, and ẑ are the unit vectors in Cartesian coordinates. The linear transformation to this right-handed coordinate triplet is a rotation matrix,

$${\displaystyle R={\begin{pmatrix}\sin \theta \cos \varphi &\sin \theta \sin \varphi &{\hphantom {-}}\cos \theta \\\cos \theta \cos \varphi &\cos \theta \sin \varphi &-\sin \theta \\-\sin \varphi &\cos \varphi &{\hphantom {-}}0\end{pmatrix}}.}$$

This gives the transformation from the spherical to the cartesian, the other way around is given by its inverse. Note: the matrix is an orthogonal matrix, that is, its inverse is simply its transpose.

The Cartesian unit vectors are thus related to the spherical unit vectors by:

$${\displaystyle {\begin{bmatrix}\mathbf {\hat {x}} \\\mathbf {\hat {y}} \\\mathbf {\hat {z}} \end{bmatrix}}={\begin{bmatrix}\sin \theta \cos \varphi &\cos \theta \cos \varphi &-\sin \varphi \\\sin \theta \sin \varphi &\cos \theta \sin \varphi &{\hphantom {-}}\cos \varphi \\\cos \theta &-\sin \theta &{\hphantom {-}}0\end{bmatrix}}{\begin{bmatrix}{\boldsymbol {\hat {r}}}\\{\boldsymbol {\hat {\theta }}}\\{\boldsymbol {\hat {\varphi }}}\end{bmatrix}}}$$

The general form of the formula to prove the differential line element, is^[5]

$${\displaystyle \mathrm {d} \mathbf {r} =\sum _{i}{\frac {\partial \mathbf {r} }{\partial x_{i}}}\,\mathrm {d} x_{i}=\sum _{i}\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|{\frac {\frac {\partial \mathbf {r} }{\partial x_{i}}}{\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|}}\,\mathrm {d} x_{i}=\sum _{i}\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|\,\mathrm {d} x_{i}\,{\hat {\boldsymbol {x}}}_{i},}$$

that is, the change in ${\displaystyle \mathbf {r} }$ is decomposed into individual changes corresponding to changes in the individual coordinates.

To apply this to the present case, one needs to calculate how ${\displaystyle \mathbf {r} }$ changes with each of the coordinates. In the conventions used,

$${\displaystyle \mathbf {r} ={\begin{bmatrix}r\sin \theta \,\cos \varphi \\r\sin \theta \,\sin \varphi \\r\cos \theta \end{bmatrix}}.}$$

Thus,

$${\displaystyle {\frac {\partial \mathbf {r} }{\partial r}}={\begin{bmatrix}\sin \theta \,\cos \varphi \\\sin \theta \,\sin \varphi \\\cos \theta \end{bmatrix}}=\mathbf {\hat {r}} ,\quad {\frac {\partial \mathbf {r} }{\partial \theta }}={\begin{bmatrix}r\cos \theta \,\cos \varphi \\r\cos \theta \,\sin \varphi \\-r\sin \theta \end{bmatrix}}=r\,{\hat {\boldsymbol {\theta }}},\quad {\frac {\partial \mathbf {r} }{\partial \varphi }}={\begin{bmatrix}-r\sin \theta \,\sin \varphi \\{\hphantom {-}}r\sin \theta \,\cos \varphi \\0\end{bmatrix}}=r\sin \theta \,\mathbf {\hat {\boldsymbol {\varphi }}} .}$$

The desired coefficients are the magnitudes of these vectors:^[5]

$${\displaystyle \left|{\frac {\partial \mathbf {r} }{\partial r}}\right|=1,\quad \left|{\frac {\partial \mathbf {r} }{\partial \theta }}\right|=r,\quad \left|{\frac {\partial \mathbf {r} }{\partial \varphi }}\right|=r\sin \theta .}$$

The surface element spanning from θ to θ + dθ and φ to φ + dφ on a spherical surface at (constant) radius r is then

$${\displaystyle \mathrm {d} S_{r}=\left\|{\frac {\partial {\mathbf {r} }}{\partial \theta }}\times {\frac {\partial {\mathbf {r} }}{\partial \varphi }}\right\|\mathrm {d} \theta \,\mathrm {d} \varphi =\left|r{\hat {\boldsymbol {\theta }}}\times r\sin \theta {\boldsymbol {\hat {\varphi }}}\right|=r^{2}\sin \theta \,\mathrm {d} \theta \,\mathrm {d} \varphi ~.}$$

Thus the differential solid angle is

$${\displaystyle \mathrm {d} \Omega ={\frac {\mathrm {d} S_{r}}{r^{2}}}=\sin \theta \,\mathrm {d} \theta \,\mathrm {d} \varphi .}$$

The surface element in a surface of polar angle θ constant (a cone with vertex at the origin) is

$${\displaystyle \mathrm {d} S_{\theta }=r\sin \theta \,\mathrm {d} \varphi \,\mathrm {d} r.}$$

The surface element in a surface of azimuth φ constant (a vertical half-plane) is

$${\displaystyle \mathrm {d} S_{\varphi }=r\,\mathrm {d} r\,\mathrm {d} \theta .}$$

The

Jacobian matrix of partial derivatives

,

$${\displaystyle J={\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}={\begin{pmatrix}\sin \theta \cos \varphi &r\cos \theta \cos \varphi &-r\sin \theta \sin \varphi \\\sin \theta \sin \varphi &r\cos \theta \sin \varphi &{\hphantom {-}}r\sin \theta \cos \varphi \\\cos \theta &-r\sin \theta &{\hphantom {-}}0\end{pmatrix}},}$$

namely

$${\displaystyle \mathrm {d} V=\left|{\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}\right|\,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi =r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi =r^{2}\,\mathrm {d} r\,\mathrm {d} \Omega ~.}$$

Thus, for example, a function f(r, θ, φ) can be integrated over every point in R³ by the triple integral

$${\displaystyle \int \limits _{0}^{2\pi }\int \limits _{0}^{\pi }\int \limits _{0}^{\infty }f(r,\theta ,\varphi )r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi ~.}$$

The

Laplacian

for scalar fields,

$${\displaystyle {\begin{aligned}\nabla f&={\partial f \over \partial r}{\hat {\mathbf {r} }}+{1 \over r}{\partial f \over \partial \theta }{\hat {\boldsymbol {\theta }}}+{1 \over r\sin \theta }{\partial f \over \partial \varphi }{\hat {\boldsymbol {\varphi }}},\\[8pt]\nabla ^{2}f&={1 \over r^{2}}{\partial \over \partial r}\left(r^{2}{\partial f \over \partial r}\right)+{1 \over r^{2}\sin \theta }{\partial \over \partial \theta }\left(\sin \theta {\partial f \over \partial \theta }\right)+{1 \over r^{2}\sin ^{2}\theta }{\partial ^{2}f \over \partial \varphi ^{2}}\\[8pt]&=\left({\frac {\partial ^{2}}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial }{\partial r}}\right)f+{1 \over r^{2}\sin \theta }{\partial \over \partial \theta }\left(\sin \theta {\frac {\partial }{\partial \theta }}\right)f+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}}{\partial \varphi ^{2}}}f~,\\[8pt]\end{aligned}}}$$

And it leads to the following expressions for the divergence and curl of vector fields,

$${\displaystyle \nabla \cdot \mathbf {A} ={\frac {1}{r^{2}}}{\partial \over \partial r}\left(r^{2}A_{r}\right)+{\frac {1}{r\sin \theta }}{\partial \over \partial \theta }\left(\sin \theta A_{\theta }\right)+{\frac {1}{r\sin \theta }}{\partial A_{\varphi } \over \partial \varphi },}$$

$${\displaystyle {\begin{aligned}\nabla \times \mathbf {A} ={}&{\frac {1}{r\sin \theta }}\left[{\partial \over \partial \theta }\left(A_{\varphi }\sin \theta \right)-{\partial A_{\theta } \over \partial \varphi }\right]{\hat {\mathbf {r} }}\\[4pt]&{}+{\frac {1}{r}}\left[{1 \over \sin \theta }{\partial A_{r} \over \partial \varphi }-{\partial \over \partial r}\left(rA_{\varphi }\right)\right]{\hat {\boldsymbol {\theta }}}\\[4pt]&{}+{\frac {1}{r}}\left[{\partial \over \partial r}\left(rA_{\theta }\right)-{\partial A_{r} \over \partial \theta }\right]{\hat {\boldsymbol {\varphi }}},\end{aligned}}}$$

Further, the inverse Jacobian in Cartesian coordinates is

$${\displaystyle J^{-1}={\begin{pmatrix}{\dfrac {x}{r}}&{\dfrac {y}{r}}&{\dfrac {z}{r}}\\\\{\dfrac {xz}{r^{2}{\sqrt {x^{2}+y^{2}}}}}&{\dfrac {yz}{r^{2}{\sqrt {x^{2}+y^{2}}}}}&{\dfrac {-\left(x^{2}+y^{2}\right)}{r^{2}{\sqrt {x^{2}+y^{2}}}}}\\\\{\dfrac {-y}{x^{2}+y^{2}}}&{\dfrac {x}{x^{2}+y^{2}}}&0\end{pmatrix}}.}$$

The metric tensor in the spherical coordinate system is ${\displaystyle g=J^{T}J}$.

Distance in spherical coordinates

In spherical coordinates, given two points with φ being the azimuthal coordinate

$${\displaystyle {\begin{aligned}{\mathbf {r} }&=(r,\theta ,\varphi ),\\{\mathbf {r} '}&=(r',\theta ',\varphi ')\end{aligned}}}$$

The distance between the two points can be expressed as^[6]

$${\displaystyle {\begin{aligned}{\mathbf {D} }&={\sqrt {r^{2}+r'^{2}-2rr'(\sin {\theta }\sin {\theta '}\cos {(\varphi -\varphi ')}+\cos {\theta }\cos {\theta '})}}\end{aligned}}}$$

Kinematics

In spherical coordinates, the position of a point or particle (although better written as a triple${\displaystyle (r,\theta ,\varphi )}$) can be written as^[7]

$${\displaystyle \mathbf {r} =r\mathbf {\hat {r}} .}$$

Its velocity is then^[7]

$${\displaystyle \mathbf {v} ={\frac {\mathrm {d} \mathbf {r} }{\mathrm {d} t}}={\dot {r}}\mathbf {\hat {r}} +r\,{\dot {\theta }}\,{\hat {\boldsymbol {\theta }}}+r\,{\dot {\varphi }}\sin \theta \,\mathbf {\hat {\boldsymbol {\varphi }}} }$$

and its acceleration is^[7]

$${\displaystyle {\begin{aligned}\mathbf {a} ={}&{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}\\[1ex]={}&{\hphantom {+}}\;\left({\ddot {r}}-r\,{\dot {\theta }}^{2}-r\,{\dot {\varphi }}^{2}\sin ^{2}\theta \right)\mathbf {\hat {r}} \\&{}+\left(r\,{\ddot {\theta }}+2{\dot {r}}\,{\dot {\theta }}-r\,{\dot {\varphi }}^{2}\sin \theta \cos \theta \right){\hat {\boldsymbol {\theta }}}\\&{}+\left(r{\ddot {\varphi }}\,\sin \theta +2{\dot {r}}\,{\dot {\varphi }}\,\sin \theta +2r\,{\dot {\theta }}\,{\dot {\varphi }}\,\cos \theta \right){\hat {\boldsymbol {\varphi }}}\end{aligned}}}$$

The angular momentum is

$${\displaystyle \mathbf {L} =\mathbf {r} \times \mathbf {p} =\mathbf {r} \times m\mathbf {v} =mr^{2}\left(-{\dot {\varphi }}\sin \theta \,\mathbf {\hat {\boldsymbol {\theta }}} +{\dot {\theta }}\,{\hat {\boldsymbol {\varphi }}}\right)}$$

Where ${\displaystyle m}$ is mass. In the case of a constant φ or else θ = π/2, this reduces to vector calculus in polar coordinates.

The corresponding angular momentum operator then follows from the phase-space reformulation of the above,

$${\displaystyle \mathbf {L} =-i\hbar ~\mathbf {r} \times \nabla =i\hbar \left({\frac {\hat {\boldsymbol {\theta }}}{\sin(\theta )}}{\frac {\partial }{\partial \phi }}-{\hat {\boldsymbol {\phi }}}{\frac {\partial }{\partial \theta }}\right).}$$

The torque is given as^[7]

$${\displaystyle \mathbf {\tau } ={\frac {\mathrm {d} \mathbf {L} }{\mathrm {d} t}}=\mathbf {r} \times \mathbf {F} =-m\left(2r{\dot {r}}{\dot {\varphi }}\sin \theta +r^{2}{\ddot {\varphi }}\sin {\theta }+2r^{2}{\dot {\theta }}{\dot {\varphi }}\cos {\theta }\right){\hat {\boldsymbol {\theta }}}+m\left(r^{2}{\ddot {\theta }}+2r{\dot {r}}{\dot {\theta }}-r^{2}{\dot {\varphi }}^{2}\sin \theta \cos \theta \right){\hat {\boldsymbol {\varphi }}}}$$

The kinetic energy is given as^[7]

$${\displaystyle E_{k}={\frac {1}{2}}m\left[\left({\dot {r}}\right)^{2}+\left(r{\dot {\theta }}\right)^{2}+\left(r{\dot {\varphi }}\sin \theta \right)^{2}\right]}$$

See also

• Celestial coordinate system
   – System for specifying positions of celestial objectsPages displaying short descriptions of redirect targets
• Coordinate system – Method for specifying point positions
• Del in cylindrical and spherical coordinates – Mathematical gradient operator in certain coordinate systems
• Double Fourier sphere method
• Elevation (ballistics) – Angle in ballistics
• Euler angles – Description of the orientation of a rigid body
• Gimbal lock – Loss of one degree of freedom in a three-dimensional, three-gimbal mechanism
• Hypersphere
   – Generalized sphere of dimension n (mathematics)Pages displaying short descriptions of redirect targets
• Jacobian matrix and determinant – Matrix of all first-order partial derivatives of a vector-valued function
• List of canonical coordinate transformations
• Sphere – Set of points equidistant from a center
• Spherical harmonic
   – Special mathematical functions defined on the surface of a spherePages displaying short descriptions of redirect targets
• Theodolite – Optical surveying instrument
• Vector fields in cylindrical and spherical coordinates – Vector field representation in 3D curvilinear coordinate systems
• Yaw, pitch, and roll
   – Principal directions in aviationPages displaying short descriptions of redirect targets

Notes

1. ^ "ISO 80000-2:2019 Quantities and units – Part 2: Mathematics". ISO. 19 May 2020. pp. 20–21. Item no. 2-17.3. Retrieved 2020-08-12.
2. ^ Duffett-Smith, P and Zwart, J, p. 34.
3. ^ ^{a b} Eric W. Weisstein (2005-10-26). "Spherical Coordinates". MathWorld. Retrieved 2010-01-15.
4. ^ "Video Game Math: Polar and Spherical Notation". Academy of Interactive Entertainment (AIE). Retrieved 2022-02-16.
5. ^ ^{a b} "Line element (dl) in spherical coordinates derivation/diagram". Stack Exchange. October 21, 2011.
6. ^ https://math.stackexchange.com/questions/833002/distance-between-two-points-in-spherical-coordinates
7. ^
   OCLC 1104053368.{{cite book}}: CS1 maint: location (link) CS1 maint: location missing publisher (link
   )

Bibliography

• Iyanaga, Shōkichi; Kawada, Yukiyosi (1977).
  ISBN 978-0262090162
  .
• LCCN 52011515
  .
• LCCN 55010911
  .
• Korn GA,
  LCCN 59014456
  . ASIN B0000CKZX7.
• Sauer R, Szabó I (1967). Mathematische Hilfsmittel des Ingenieurs. New York: Springer Verlag. pp. 95–96.
  LCCN 67025285
  .
• Moon P, Spencer DE (1988). "Spherical Coordinates (r, θ, ψ)". Field Theory Handbook, Including Coordinate Systems, Differential Equations, and Their Solutions (corrected 2nd ed., 3rd print ed.). New York: Springer-Verlag. pp. 24–27 (Table 1.05).
  ISBN 978-0-387-18430-2
  .
• Duffett-Smith P, Zwart J (2011). Practical Astronomy with your Calculator or Spreadsheet, 4th Edition. New York: Cambridge University Press. p. 34.
  ISBN 978-0521146548
  .

External links

• "Spherical coordinates", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• MathWorld description of spherical coordinates
• Coordinate Converter – converts between polar, Cartesian and spherical coordinates