Fact-checked by Grok 2 weeks ago

Ellipsoidal coordinates

Ellipsoidal coordinates, also known as confocal ellipsoidal coordinates, constitute a three-dimensional orthogonal curvilinear coordinate system defined by the intersection of three families of confocal quadrics: ellipsoids, hyperboloids of one sheet, and hyperboloids of two sheets, all sharing the same foci.^[1] In this system, a point in space is specified by coordinates (\lambda, \mu, \nu) satisfying \lambda < c^2 < \mu < b^2 < \nu < a^2 (assuming a > b > c > 0), where the Cartesian coordinates (x, y, z) relate via the equations \frac{x^2}{a^2 - \lambda} + \frac{y^2}{b^2 - \lambda} + \frac{z^2}{c^2 - \lambda} = 1 and analogous forms for \mu and \nu.^[1] Introduced by Carl Gustav Jacob Jacobi in 1839 as part of his solution to the problem of geodesics on a triaxial ellipsoid, these coordinates extend the two-dimensional elliptical system to three dimensions and enable the separation of variables in key partial differential equations.^[2] The system's orthogonality arises from the confocal nature of the quadrics, ensuring that the coordinate surfaces intersect at right angles, with scale factors given by h_\lambda = \sqrt{\frac{(\mu - \lambda)(\nu - \lambda)}{4(a^2 - \lambda)(b^2 - \lambda)(c^2 - \lambda)}}, and similar expressions for h_\mu and h_\nu.^[1] This property makes ellipsoidal coordinates particularly valuable in mathematical physics for solving boundary value problems in domains bounded by ellipsoids, as the Laplacian and Helmholtz operators separate completely in these variables.^[1] Historically termed "elliptic coordinates" by some authors like Hilbert and Cohn-Vossen, the system has been formalized in texts such as Morse and Feshbach's Methods of Theoretical Physics (1953), where it is simply called "ellipsoidal coordinates."^[1] Notable applications include the analytical solution of the Schrödinger equation for two-center quantum systems, such as the hydrogen molecular ion, where the non-spherical symmetry aligns with the coordinate surfaces.^[3] In geophysics and geodesy, ellipsoidal coordinates facilitate computations of geodesics and gravitational potentials on triaxial ellipsoids modeling planetary shapes.^[2] Additionally, they underpin the theory of ellipsoidal harmonics, used in solving elliptic boundary value problems and in fields like electroencephalography for modeling bioelectric fields within ellipsoidal approximations of the human head.^[4]^[5] Despite their utility, the complexity of the metric and associated special functions, such as Lamé functions, limits their routine use compared to spherical or cylindrical systems, though numerical advancements continue to expand their practical scope.^[6]

Overview

Definition and Properties

Ellipsoidal coordinates, also known as confocal ellipsoidal coordinates, form a three-dimensional orthogonal curvilinear coordinate system in Euclidean space, typically parameterized by variables (λ, μ, ν). This system generalizes the two-dimensional elliptic coordinates and features coordinate surfaces belonging to three families of confocal quadrics: ellipsoids (for constant λ), hyperboloids of one sheet (for constant μ), and hyperboloids of two sheets (for constant ν).^[1]^[4] A defining property of ellipsoidal coordinates is their orthogonality, whereby the gradients of the coordinate functions are mutually perpendicular at every point, ensuring that the coordinate surfaces intersect at right angles throughout space. The confocal property means that all quadrics share a common set of six foci, usually positioned along the principal axes at (±h₁, 0, 0), (0, ±h₂, 0), and (0, 0, ±h₃), where the h_i represent semi-focal distances derived from the differences in the squares of the semi-axes. Geometrically, these coordinates can be intuited as arising from the intersections of surfaces defined by ratios of distances from a point to focal conics—such as ellipses and hyperbolas—in the coordinate planes, providing a natural partitioning of space aligned with ellipsoidal symmetry.^[1]^[4] This orthogonal and confocal structure renders ellipsoidal coordinates especially suitable for analytical solutions to problems exhibiting ellipsoidal boundaries or symmetry, including boundary value problems in electrostatics, gravitation, and quantum mechanics where partial differential equations like Laplace's or the Helmholtz equation separate in these variables. In contrast to non-confocal systems like geodetic coordinates, which are tailored to a fixed reference ellipsoid for Earth-surface modeling in surveying and are not based on intersecting families of quadrics, ellipsoidal coordinates offer a versatile, space-filling framework independent of any specific physical reference.^[1]^[7]^[4]

Historical Development

Ellipsoidal coordinates were introduced by Carl Gustav Jacob Jacobi in 1839, in his foundational work on geodesics on a triaxial ellipsoid, where he employed these coordinates to achieve separation of variables and reduce the problem to elliptic integrals.^[8] Independently in the same year, Gabriel Lamé utilized ellipsoidal coordinates to separate variables in the Laplace equation, thereby establishing the basis for Lamé functions as solutions to partial differential equations in ellipsoidal geometries.^[9] This three-dimensional system evolved from two-dimensional elliptic coordinates, which had earlier been applied to solve boundary value problems involving confocal ellipses and hyperbolas in the plane.^[10] The influence of these coordinates extended to the development of Lamé functions, which facilitated the separation of variables in Helmholtz's equation and other PDEs within ellipsoidal domains.^[11] In the early 20th century, David Hilbert and Richard Courant referenced ellipsoidal coordinates in their systematic approach to boundary value problems, detailing their properties in Methods of Mathematical Physics (1924).^[12] Their adoption in quantum mechanics became prominent with Philip M. Morse and Herman Feshbach's Methods of Theoretical Physics (1953), where the coordinates were used to address separable potentials in atomic and molecular problems.^[13] Notable later texts include Lev D. Landau and Evgeny M. Lifshitz's Electrodynamics of Continuous Media (1984), which employed an alternative notation for ellipsoidal coordinates in discussions of dielectric responses and field problems in anisotropic media.^[14] These developments underscore the coordinates' enduring role in mathematical physics, particularly for applications in potential theory.

Geometry

Confocal Quadrics

Confocal quadrics form the geometric foundation of the ellipsoidal coordinate system, consisting of a one-parameter family of quadratic surfaces that share the same pair of focal conics: an ellipse in the xy-plane and a hyperbola in the xz-plane. These surfaces encompass ellipsoids, hyperboloids of one sheet, and hyperboloids of two sheets, depending on the value of the parameter labeling each surface.^[15] The mathematical representation of these confocal quadrics is given by the equation

\frac{x^2}{a^2 + \sigma} + \frac{y^2}{b^2 + \sigma} + \frac{z^2}{c^2 + \sigma} = 1,

where a > b > c > 0 are fixed parameters representing the semi-axes lengths of a reference ellipsoid, and \sigma varies over appropriate intervals to generate the family.^[4] The focal conics arise as degenerate cases of this equation when \sigma takes specific values related to -a^2, -b^2, and -c^2, with their locations and shapes derived from the inter-focal distances \sqrt{a^2 - b^2}, \sqrt{a^2 - c^2}, and \sqrt{b^2 - c^2}.^[15] In the ellipsoidal coordinate system, each coordinate value specifies a particular quadric surface from one of three orthogonal families generated by this equation, enabling the parametrization of space through intersections of these surfaces.^[4]

Coordinate Surfaces

In ellipsoidal coordinates (\lambda, \mu, \nu), the coordinate surfaces consist of three families of confocal quadrics that intersect orthogonally to define the coordinate system. These surfaces share the same foci, a property that underpins their geometric utility in three-dimensional space. With -c^2 < \lambda < \infty, -b^2 < \mu < -c^2, and -a^2 < \nu < -b^2, the surfaces are as follows.^[1] Surfaces of constant \lambda are ellipsoids described by the equation

\frac{x^2}{a^2 + \lambda} + \frac{y^2}{b^2 + \lambda} + \frac{z^2}{c^2 + \lambda} = 1,

where a > b > c > 0 are fixed semi-axes parameters. These ellipsoids are nested, with increasing \lambda producing successively larger surfaces that enclose one another, filling the exterior region beyond the focal conic.^[1]^[4] Surfaces of constant \mu form hyperboloids of one sheet given by

\frac{x^2}{a^2 + \mu} + \frac{y^2}{b^2 + \mu} - \frac{z^2}{c^2 + \mu} = 1.

These hyperboloids connect opposite branches along the principal axes and intersect the ellipsoidal surfaces transversely.^[1] Surfaces of constant \nu are hyperboloids of two sheets expressed as

\frac{x^2}{a^2 + \nu} - \frac{y^2}{b^2 + \nu} - \frac{z^2}{c^2 + \nu} = 1.

These disconnected sheets align primarily along the longest axis and complement the one-sheet hyperboloids in partitioning space.^[1] The orthogonal intersections among these families—ellipsoids with each type of hyperboloid—generate a curvilinear grid that conforms to the ellipsoidal symmetry, enabling efficient representation of fields varying along these natural geometric boundaries.^[1]^[4]

Formulation

Transformation Equations

Ellipsoidal coordinates (\lambda, \mu, \nu) are defined with respect to a reference ellipsoid characterized by semi-axes lengths a > b > c > 0 along the x, y, and z directions, respectively. These parameters represent the confocal family of ellipsoids, hyperboloids of one sheet, and hyperboloids of two sheets, respectively.^[16] The transformation from ellipsoidal coordinates to Cartesian coordinates is given by the following equations:

\begin{align} x^2 &= \frac{(a^2 + \lambda)(a^2 + \mu)(a^2 + \nu)}{(a^2 - b^2)(a^2 - c^2)}, \\ y^2 &= \frac{(b^2 + \lambda)(b^2 + \mu)(b^2 + \nu)}{(b^2 - a^2)(b^2 - c^2)}, \\ z^2 &= \frac{(c^2 + \lambda)(c^2 + \mu)(c^2 + \nu)}{(c^2 - b^2)(c^2 - a^2)}. \end{align}

^[16] Since the equations express squared Cartesian coordinates, the signs of x, y, and z are chosen independently as positive or negative to specify the position within one of the eight octants of space. This allows the coordinate system to cover the entire three-dimensional space without singularities except along the degenerate surfaces where the denominators vanish.

Range of Coordinates

In ellipsoidal coordinates, the parameters satisfy \lambda > \mu > \nu, with \lambda in the ellipsoids family, \mu in the hyperboloids of one sheet family, and \nu in the hyperboloids of two sheets family. The specific ranges are -c^2 < \lambda < \infty, -b^2 < \mu < -c^2, and -a^2 < \nu < -b^2, where a > b > c > 0 are the semi-axes of the reference confocal ellipsoid. These disjoint intervals ensure that each coordinate corresponds uniquely to one family of confocal quadrics without overlap.^[16] The full three-dimensional Euclidean space is covered by these ranges, with the sign choices in the Cartesian transformation allowing coverage of all octants. Uniqueness holds such that, except at the focal conics themselves, every point in space (x, y, z) \neq (0, 0, 0) admits a unique (\lambda, \mu, \nu) in the specified ranges, modulo the octant symmetry; this one-to-one correspondence (up to signs) facilitates the orthogonal decomposition of space and separation of variables in the Hamilton-Jacobi equation.^[16] Boundary conditions occur at the limits of these ranges: for instance, \lambda = -c^2 defines the degenerate case where the \lambda-coordinate surface (ellipsoid) collapses to the focal ellipse in the xy-plane; \mu = -b^2 and \nu = -a^2 correspond to degenerate hyperboloids collapsing to the focal conics in the xz- and yz-planes, respectively. These singular curves serve as the common foci for all confocal quadrics in the system.

Orthogonal Properties

Scale Factors

In ellipsoidal coordinates (\lambda, \mu, \nu), the scale factors h_\lambda, h_\mu, and h_\nu are the diagonal elements of the metric tensor, determining the infinitesimal arc length along each coordinate direction. These factors are crucial for expressing the geometry of the coordinate system and are derived from the position vector \mathbf{r} = (x, y, z) expressed in terms of the coordinates. The line element takes the orthogonal form

ds^2 = h_\lambda^2 \, d\lambda^2 + h_\mu^2 \, d\mu^2 + h_\nu^2 \, d\nu^2,

where the scale factors are given by

h_\lambda = \frac{1}{2} \sqrt{ \frac{(\lambda - \mu)(\lambda - \nu)}{S(\lambda)} },

with the auxiliary function S(\sigma) = (\sigma + a^2)(\sigma + b^2)(\sigma + c^2) and a > b > c > 0 the semi-axes lengths of the confocal reference ellipsoid. The expressions for the other scale factors follow by appropriate permutation to ensure positivity given the ordering \lambda \geq \mu \geq \nu:

h_\mu = \frac{1}{2} \sqrt{ \frac{(\mu - \nu)(\lambda - \mu)}{S(\mu)} }, \quad h_\nu = \frac{1}{2} \sqrt{ \frac{(\lambda - \nu)(\mu - \nu)}{S(\nu)} }.

These formulas arise from the standard transformation between Cartesian and ellipsoidal coordinates, where x, y, z are expressed using products involving ( \sigma + \alpha^2 ) for \sigma = \lambda, \mu, \nu and \alpha = a, b, c. To derive the scale factors, one computes the partial derivatives \partial \mathbf{r}/\partial \lambda, \partial \mathbf{r}/\partial \mu, and \partial \mathbf{r}/\partial \nu from the coordinate transformation equations. Each scale factor is then h_i = \left| \partial \mathbf{r} / \partial q_i \right|, where q_i denotes the respective coordinate. Squaring these magnitudes yields h_i^2 = \left( \partial x / \partial q_i \right)^2 + \left( \partial y / \partial q_i \right)^2 + \left( \partial z / \partial q_i \right)^2. Substituting the explicit forms of x, y, z (which involve square roots of products like (\lambda + a^2)(\mu + a^2)(\nu + a^2) normalized by differences of the semi-axes squared) and simplifying using the definitions of the confocal quadrics results in the expressions above after algebraic manipulation and taking square roots. The factor of $1/2 emerges from the differentiation of the nested square roots in the transformation. Orthogonality of the system ensures the off-diagonal metric terms vanish. In asymptotic regimes near the foci of the confocal quadrics, the scale factors simplify significantly. For instance, as coordinates approach the degenerate cases corresponding to the focal ellipse or hyperbola (e.g., when \lambda, \mu, \nu tend toward -a^2, -b^2, -c^2), the expressions reduce to those of limiting coordinate systems such as elliptic cylindrical or spherical coordinates, where one or more scale factors become constant or follow simpler hyperbolic/trigonometric forms. This behavior facilitates approximations in regions close to the origin or along the principal axes.

Line Element

In ellipsoidal coordinates (\lambda, \mu, \nu), the infinitesimal line element ds expresses the squared distance between neighboring points and takes the form

ds^2 = h_\lambda^2 \, d\lambda^2 + h_\mu^2 \, d\mu^2 + h_\nu^2 \, d\nu^2,

where h_\lambda, h_\mu, and h_\nu are the scale factors associated with each coordinate direction.^[17] The corresponding infinitesimal volume element dV is the product of the scale factors times the coordinate differentials:

dV = h_\lambda h_\mu h_\nu \, d\lambda \, d\mu \, d\nu.

Substituting the explicit expressions for the scale factors yields

dV = \frac{(\lambda - \mu)(\lambda - \nu)(\mu - \nu)}{8 \sqrt{S(\lambda) S(\mu) S(\nu)}} \, d\lambda \, d\mu \, d\nu,

where S(\xi) = (\xi + a^2)(\xi + b^2)(\xi + c^2) for semi-axes lengths a > b > c > 0, and the ordering \lambda > \mu > \nu ensures dV > 0.^[18] The infinitesimal surface elements on the coordinate surfaces are products of the two relevant scale factors and differentials. On the ellipsoid surface of constant \lambda,

dS_\lambda = h_\mu h_\nu \, d\mu \, d\nu.

On the hyperboloid surface of constant \mu,

dS_\mu = h_\lambda h_\nu \, d\lambda \, d\nu,

and on the hyperboloid surface of constant \nu,

dS_\nu = h_\lambda h_\mu \, d\lambda \, d\mu.

These expressions arise directly from the orthogonality of the coordinate system.^[17] The Jacobian determinant of the coordinate transformation from (\lambda, \mu, \nu) to Cartesian (x, y, z) coordinates is |h_\lambda h_\mu h_\nu|, providing the scaling factor for volume integrals in the ellipsoidal system.^[18]

Differential Operators

Gradient and Divergence

In ellipsoidal coordinates (\lambda, \mu, \nu), the gradient of a scalar function f is expressed using the scale factors h_\lambda, h_\mu, and h_\nu as

\nabla f = \frac{1}{h_\lambda} \frac{\partial f}{\partial \lambda} \hat{e}_\lambda + \frac{1}{h_\mu} \frac{\partial f}{\partial \mu} \hat{e}_\mu + \frac{1}{h_\nu} \frac{\partial f}{\partial \nu} \hat{e}_\nu,

where \hat{e}_\lambda, \hat{e}_\mu, and \hat{e}_\nu are the unit vectors along the respective coordinate directions.^[17] The divergence of a vector field \mathbf{A} = A_\lambda \hat{e}_\lambda + A_\mu \hat{e}_\mu + A_\nu \hat{e}_\nu takes the form

\nabla \cdot \mathbf{A} = \frac{1}{h_\lambda h_\mu h_\nu} \left[ \frac{\partial (h_\mu h_\nu A_\lambda)}{\partial \lambda} + \frac{\partial (h_\lambda h_\nu A_\mu)}{\partial \mu} + \frac{\partial (h_\lambda h_\mu A_\nu)}{\partial \nu} \right].

This expression accounts for the variation in the scale factors across the coordinate surfaces.^[17] The curl of \mathbf{A} is more involved due to the orthogonal but curvilinear nature of the system and is given by

\nabla \times \mathbf{A} = \frac{1}{h_\lambda h_\mu h_\nu} \begin{vmatrix} h_\lambda \hat{e}_\lambda & h_\mu \hat{e}_\mu & h_\nu \hat{e}_\nu \\ \frac{\partial}{\partial \lambda} & \frac{\partial}{\partial \mu} & \frac{\partial}{\partial \nu} \\ h_\lambda A_\lambda & h_\mu A_\mu & h_\nu A_\nu \end{vmatrix},

which expands to components involving cross-derivatives scaled by products of the h_i.^[17] These operators facilitate the analysis of vector and scalar fields exhibiting ellipsoidal symmetry, such as gravitational or electrostatic potentials in triaxial configurations, by aligning with the natural geometry of confocal quadrics.^[17]

Laplacian

In orthogonal curvilinear coordinates such as the ellipsoidal system denoted by (λ, μ, ν), the Laplacian operator ∇² acting on a scalar function φ is given by the general expression

\nabla^2 \phi = \frac{1}{h_\lambda h_\mu h_\nu} \left[ \frac{\partial}{\partial \lambda} \left( \frac{h_\mu h_\nu}{h_\lambda} \frac{\partial \phi}{\partial \lambda} \right) + \frac{\partial}{\partial \mu} \left( \frac{h_\lambda h_\nu}{h_\mu} \frac{\partial \phi}{\partial \mu} \right) + \frac{\partial}{\partial \nu} \left( \frac{h_\lambda h_\mu}{h_\nu} \frac{\partial \phi}{\partial \nu} \right) \right],

where h_λ, h_μ, and h_ν are the scale factors associated with each coordinate direction.^[19] The scale factors for ellipsoidal coordinates are

h_\lambda = \frac{1}{2} \sqrt{\frac{(\mu - \lambda)(\nu - \lambda)}{S(\lambda)}}, \quad h_\mu = \frac{1}{2} \sqrt{\frac{(\lambda - \mu)(\nu - \mu)}{S(\mu)}}, \quad h_\nu = \frac{1}{2} \sqrt{\frac{(\lambda - \nu)(\mu - \nu)}{S(\nu)}},

with S(\sigma) = (\sigma - a^2)(\sigma - b^2)(\sigma - c^2) defined using the squared interfocal distances a^2 > b^2 > c^2 > 0.^[20] Substituting these scale factors into the general Laplacian yields an expanded form involving

\nabla^2 \phi = \frac{4 \sqrt{ -S(\lambda) S(\mu) S(\nu) } }{ (\lambda - \mu)(\mu - \nu)(\nu - \lambda) } \left[ \frac{\partial}{\partial \lambda} \left( \sqrt{ \frac{ (\lambda - \mu)(\lambda - \nu) }{ S(\lambda) } } \frac{\partial \phi}{\partial \lambda} \right) + \frac{\partial}{\partial \mu} \left( \sqrt{ \frac{ (\mu - \lambda)(\mu - \nu) }{ S(\mu) } } \frac{\partial \phi}{\partial \mu} \right) + \frac{\partial}{\partial \nu} \left( \sqrt{ \frac{ (\nu - \lambda)(\nu - \mu) }{ S(\nu) } } \frac{\partial \phi}{\partial \nu} \right) \right].

This expression accounts for the coordinate ordering λ > μ > ν and the signs within the square roots, ensuring positive scale factors.^[20] The form of the Laplacian facilitates the additive separation of variables for the Helmholtz equation ∇²φ + k²φ = 0 in ellipsoidal coordinates, resulting in ordinary differential equations solvable by Lamé functions.^[21] In the limiting case where the interfocal distances vanish (a → b → c), ellipsoidal coordinates degenerate to spherical coordinates, and the Laplacian reduces to the familiar spherical form involving radial and angular derivatives.^[20]

Alternative Representations

Angular Parametrization

In the angular parametrization of ellipsoidal coordinates, a non-orthogonal system is employed to describe points within or on the surface of an ellipsoid, utilizing parameters (s, \theta, \phi) where s \in [0,1] serves as a radial-like parameter scaling from the origin to the ellipsoidal boundary, while \theta \in [0, \pi] and \phi \in [0, 2\pi) act as polar and azimuthal angular variables, respectively, analogous to those in spherical coordinates. The Cartesian coordinates are related by the transformation equations \begin{align*} x &= a s \sin\theta \cos\phi, \ y &= b s \sin\theta \sin\phi, \ z &= c s \cos\theta, \end{align*} where a, b, and c denote the lengths of the semi-axes along the x-, y-, and z-directions, respectively. This mapping stretches the unit ball in (s, \theta, \phi)-space affinely to fill the interior of the ellipsoid \frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2} \leq 1.^[22] Unlike the standard orthogonal ellipsoidal coordinates, which rely on confocal quadric surfaces, this parametrization centers the system at the origin and directly modifies spherical coordinates via axis-aligned scaling, resulting in non-orthogonal coordinate surfaces due to the differing semi-axes. The line element ds^2 reflects this non-orthogonality through off-diagonal metric components:

ds^2 = (a^2 \sin^2\theta \cos^2\phi + b^2 \sin^2\theta \sin^2\phi + c^2 \cos^2\theta) \, ds^2 + s^2 (a^2 \cos^2\theta \cos^2\phi + b^2 \cos^2\theta \sin^2\phi + c^2 \sin^2\theta) \, d\theta^2 + s^2 \sin^2\theta (a^2 \sin^2\phi + b^2 \cos^2\phi) \, d\phi^2

+ 2 s \sin\theta \cos\theta (a^2 \cos^2\phi + b^2 \sin^2\phi - c^2) \, ds \, d\theta + 2 s \sin^2\theta \sin\phi \cos\phi (b^2 - a^2) \, ds \, d\phi + 2 s^2 \sin\theta \cos\theta \sin\phi \cos\phi (b^2 - a^2) \, d\theta \, d\phi.

The cross terms vanish only in the degenerate spherical case where a = b = c.^[22] The volume element in these coordinates is dV = a b c \, s^2 \sin\theta \, ds \, d\theta \, d\phi, which integrates to the ellipsoid volume \frac{4}{3} \pi a b c over the full domain.^[22] This parametrization offers advantages in applications involving bounded ellipsoidal domains, such as cavities in acoustics or electromagnetics, where the finite radial range s \in [0,1] simplifies boundary specifications and numerical discretizations compared to unbounded or confocal systems, enabling efficient finite-domain simulations.^[23]

Notation Variants

Ellipsoidal coordinates are most commonly denoted by the triplet (\lambda, \mu, \nu), where the parameters satisfy the ordering -a^2 < \nu < -b^2 < \mu < -c^2 < \lambda < \infty (assuming a > b > c > 0) in some conventions, differing from the primary article range c^2 < \lambda < b^2 < \mu < a^2 < \nu < \infty. This convention aligns the coordinate surfaces with confocal ellipsoids (\lambda = constant), one-sheeted hyperboloids (\mu = constant), and two-sheeted hyperboloids (\nu = constant).^[24] An alternative notation appears in the work of Landau and Lifshitz, who employ (\xi, \eta, \zeta) with ranges -a^2 < \zeta < -b^2 < \eta < -c^2 < \xi < \infty, having units of distance squared; these are related to certain standard parameters by shifts such as \xi = \lambda + c^2, \eta = \mu + b^2, \zeta = \nu + a^2 to facilitate analysis. In Arfken's treatment, a similar variant uses (\xi_1, \xi_2, \xi_3) with ranges a^2 > \xi_3 > b^2 > \xi_2 > c^2 > \xi_1, reflecting a reversed ordering for the confocal parameters.^[17] Further variants include (u, \beta, \lambda), where u \in [0, +\infty) represents the ellipsoidal parameter, \beta \in [-\pi/2, +\pi/2] the ellipsoidal latitude, and \lambda \in (-\pi, +\pi] the ellipsoidal longitude, allowing a natural transition to oblate spheroidal systems. Another form is (\sigma_1, \sigma_2, \sigma_3) emphasizing positive parameters. Conversions between notations, such as shifts from the standard to the Landau variant, preserve the underlying confocal geometry while adapting to analytical needs.^[25]

Applications

Separation of Variables

Ellipsoidal coordinates facilitate the separation of variables for partial differential equations such as Laplace's equation \nabla^2 \phi = 0, which is particularly useful for problems involving ellipsoidal symmetry. In these coordinates (\lambda, \mu, \nu), satisfying c^2 < \lambda < b^2 < \mu < a^2 < \nu < \infty (assuming $0 < c < b < a), the Laplacian operator separates into independent ordinary differential equations (ODEs), one for each coordinate.^[26]^[27] The separated ODEs are Lamé equations, elliptic analogs of the Legendre differential equation encountered in spherical coordinates. The solutions to these ODEs are Lamé functions, denoted typically as Ec_n^m(\xi, k) or Es_n^m(\xi, k) for even and odd types, serving as the basis functions for expansions in ellipsoidal harmonics.^[26]^[27] The separation constants introduce an eigenvalue problem, where the constants ensure compatibility across the three equations, leading to a spectrum of eigenvalues for the Lamé operator. This eigenvalue structure arises from the need for single-valued, bounded solutions over the respective domains of \lambda, \mu, \nu, with the degree \sigma playing a role similar to that in spherical harmonics. The Lamé equation thus governs the spectral decomposition, with solutions classified by degree n and order m, where for each n, there are $2n+1 independent Lamé functions.^[26]^[28] Boundary conditions in ellipsoidal coordinates align naturally with the geometry of constant-coordinate surfaces, which are confocal ellipsoids (\lambda = constant), hyperboloids (\mu = constant), and another family of hyperboloids (\nu = constant). For instance, Dirichlet conditions specifying \phi = 0 on an ellipsoidal boundary \lambda = \lambda_0 can be imposed directly on the corresponding Lamé function, simplifying the solution for interior or exterior problems within ellipsoidal domains. This geometric compatibility makes ellipsoidal coordinates ideal for such boundaries, unlike Cartesian or spherical systems.^[26]^[27] Beyond boundary value problems, ellipsoidal coordinates find mathematical applications in conformal mappings and integral transforms. The coordinate system's relation to confocal quadrics enables representations of analytic functions via ellipsoidal harmonics, aiding in the study of complex potentials and mappings in three dimensions. Additionally, Lamé functions underpin integral transforms analogous to Fourier-Legendre expansions, used for inverting certain integral equations in potential theory and solving integral equations over ellipsoidal regions.^[27]^[28]

Physical Problems

Ellipsoidal coordinates find significant application in electrostatics, particularly for solving the Laplace equation around charged ellipsoidal conductors or dielectrics. The potential outside a uniformly charged or conducting ellipsoid can be expressed using ellipsoidal harmonics derived from Lamé functions, enabling efficient computation of the electric field and surface charge density. This approach, rooted in Chandrasekhar's theory, provides an exact solution for the electrostatic potential of homogeneous ellipsoids and extends to heterogeneous cases via perturbation methods, avoiding numerical approximations for symmetric geometries. Multipole expansions in ellipsoidal coordinates further facilitate the representation of the potential for non-spherical charge distributions, with Lamé functions serving as the basis for higher-order terms in the expansion. In quantum mechanics, ellipsoidal coordinates are employed to solve the Schrödinger equation for particles in triaxial potentials, such as those modeling deformed atomic nuclei. For triaxial ellipsoidal quantum billiards, which approximate the confinement in irregularly shaped nuclei, the wavefunctions separate into products of one-dimensional solutions along the confocal coordinates, revealing semiclassical features like gross-shell structures. Additionally, the three-dimensional anisotropic harmonic oscillator separates in ellipsoidal coordinates when the frequency ratios satisfy specific commensurability conditions (e.g., \omega_x : \omega_y : \omega_z = 1 : \sqrt{2} : \sqrt{3}), yielding exact energy levels and eigenfunctions that capture the effects of deformation on quantum states. In acoustics and elasticity, ellipsoidal coordinates enable the analysis of vibrations and scattering from ellipsoidal bodies using ellipsoidal harmonics. Low-frequency scattering problems for sound waves or elastic waves incident on ellipsoidal obstacles are solved by expanding the wavefield in terms of Lamé functions, providing insights into resonance modes and radiation patterns. For free vibrations of elastic ellipsoids, the method yields complete displacement solutions, accounting for the body's triaxial shape and material properties, which is crucial for modeling seismic responses or structural dynamics in non-spherical media. Although distinct from geodetic coordinates used in surveying, ellipsoidal coordinates relate to geodesy through the modeling of gravitational potentials around the Earth's reference ellipsoid, where harmonic expansions approximate terrain effects on the geoid. In numerical methods, ellipsoidal coordinates form the basis for finite element meshes in simulations involving ellipsoidal domains, such as potential theory problems, by leveraging open-source implementations of ellipsoidal harmonics for efficient discretization and boundary element solutions.

References

[1]
Confocal Ellipsoidal Coordinates -- from Wolfram MathWorld
These coordinates correspond to three confocal quadrics all sharing the same pair of foci. Surfaces of constant xi are confocal ellipsoids.
[2]
[PDF] Jacobi's solution for geodesics on a triaxial ellipsoid - arXiv
Nov 3, 2025 · On Boxing Day, 1838, Jacobi found a solution to the problem of geodesics on a triaxial ellipsoid, with the course of the geodesic and the ...
[3]
Analytical Solution of Schrödinger Equation in Ellipsoidal Coordinates
Research has shown that ellipsoidal coordinates are very useful for studying two-center problems in quantum mechanics, especially for systems with non-spherical ...
[4]
[PDF] The ellipsoidal system and its geometry
The ellipsoidal coor- dinate system is characterized by the fact that it specifies a point by using solely non-degenerate second-degree surfaces.Missing: mathematics | Show results with:mathematics
[5]
[PDF] Electric field and potential calculation for a bioelectric current dipole ...
Aug 31, 2005 · Abstract. Computer-based modelling in ellipsoidal geometry can be useful in electro- encephalography and gastrography due to the resemblance ...
[6]
Semi-analytical approach for computation of ellipsoidal harmonics
In particular, it is the most suitable system for solving problems of boundary, limit and initial values in three-dimensional domains with ellipsoidal boundary ...
[7]
[PDF] Geometric Reference Systems in Geodesy - Knowledge Bank
type of ellipsoidal coordinates for regional and global applications, and Cartesian for local ... orthogonal, meaning that surfaces of constant coordinates ...
[8]
[PDF] Moduli spaces for Lamé functions Speaker - Purdue Math
Jan 18, 2022 · It was discovered by Gabriel Lamé in 18391 when separating variables in the Laplace equation in the ellipsoidal coordinates in R3. Lamé ...
[9]
Elliptic coordinates - Encyclopedia of Mathematics
Jun 5, 2020 · In elliptic coordinates the Laplace equation allows separation of variables. Degenerate elliptic coordinates are two ... ds=(cosh2˜σ−cos2˜τ)d˜σd ...<|control11|><|separator|>
[10]
[PDF] Jacobi elliptic coordinates, functions of Heun and Lame type ... - CORE
Lamé and Heun functions arise via separation of the Laplace equation in general Jacobi ellipsoidal or conical coordinates.Missing: 2D | Show results with:2D
[11]
[PDF] methods of mathematical physics
an English "Courant-Hilbert" has become irresistible. Even so, it is hoped ... 3 Ellipsoidal coordinates. 222. 224. 226. Page 10. CONTENTS. §9 ...
[12]
[PDF] methods of theoretical physics - astrosen.unam.mx
Morse and Feshbach Methods of Theoretical Physics. Muskat Physical ... Degenerate Forms of Ellipsoidal Coordinates. Con- finence of Singularities ...
[13]
[PDF] arXiv:0902.1421v2 [math.DG] 10 Feb 2010
Feb 10, 2010 · a line tangent to two confocal quadrics and of its linear element in elliptic coordinates plays a fundamental rôle. On one hand the linear ...
[14]
Methods of theoretical physics : Morse, Philip M ... - Internet Archive
Nov 11, 2022 · Methods of theoretical physics. by: Morse, Philip M. (Philip McCord), 1903-1985. Publication date: 1953. Topics: Mathematical physics. Publisher ...
[15]
[PDF] COORDINATE SYSTEMS , ox ox - LSU Math
In Section 2.8 two sets of orthogonal confocal parabolas were described. ... (nonintcrsecting) conducting plane. 2.15 Confocal Ellipsoidal Coordinates (~,'~', ~3).
[16]
[PDF] methods of theoretical physics - astrosen.unam.mx
Morse and Feshbach Methods of Theoretical Physics. Muskat Physical Principles ... Ellipsoidal Coordinates. Problems. Trigonometric and Hyperbolic ...
[17]
5.10 Three-Dimensional Solutions to Laplace's Equation - MIT
The first step in exploiting these new systems is to write the Laplacian and other differential operators in terms of those coordinates. This is also ...
[18]
23.21 Physical Applications
§23.21(iii) Ellipsoidal Coordinates ; 1 ·, ; 2 ·, ; 3 · are constants. The Laplacian operator ∇ 2 (§1.5(ii)) is given by ...
[19]
29.18 Mathematical Applications
Contents ; §29.18(i) Sphero-Conal Coordinates ; §29.18(ii) Ellipsoidal Coordinates ; §29.18(iii) Spherical and Ellipsoidal Harmonics ; §29.18(iv) Other Applications ...
[20]
Ellipsoid -- from Wolfram MathWorld
The general ellipsoid, also called a triaxial ellipsoid, is a quadratic surface which is given in Cartesian coordinates by (x^2)/(a^2)+(y^2)/(b^2)+(z^2)/(c
[21]
[PDF] UCLA Electronic Theses and Dissertations - eScholarship
Notably, the angular parametrization of ellipsoidal coordinates parallels that of the spherical system, albeit with an additional factor of a, b and c to ...
[22]
(PDF) A Study on Geodetic Boundary Value Problems in Ellipsoidal ...
Jul 5, 2018 · The ellipsoidal coordinate system is described in some detail. For a one-to-one correspondence between ellipsoidal and Cartesian coordinates two ...Missing: 326211145 | Show results with:326211145
[23]
https://escholarship.org/content/qt73g562zg/qt73g562zg.pdf
[24]
the theory of spherical and ellipsoidal harmonics - Internet Archive
Sep 8, 2022 · the theory of spherical and ellipsoidal harmonics. by: E.W. Hobson, Sc.D., LL.D., F.R.S. Publication date: 1931. Collection ...
[25]
[PDF] methods of theoretical physics
Morse and Feshbach Methods of Theoretical Physics. Muskat Physical ... Degenerate Forms of Ellipsoidal Coordinates. Con- fluence of Singularities ...