Polyspherical Approach: Methods & Applications

• The polyspherical approach is a strategy that decomposes geometric spaces into hierarchical radial-angular layers to simplify linearization, separation of variables, and operator factorization.
• It is applied across disciplines, from conformal geometry and harmonic analysis to computational molecular dynamics and quantum algorithms.
• The method underpins both classical coordinate systems and modern multicomplex topologies, facilitating advanced analytic and numerical techniques.

Searching arXiv for recent and foundational papers on “polyspherical” to ground the article in the literature. The polyspherical approach is a family of coordinate and representation methods in which a geometric, analytic, or dynamical problem is reformulated by decomposing variables into radial–angular hierarchies or by embedding the original space into a higher-dimensional homogeneous coordinate space subject to a quadratic constraint. In Euclidean geometry and conformal theory, Darboux’s polyspherical coordinates linearize inversion and the full conformal group on an ambient cone; in harmonic analysis they organize separation of variables on spheres through Vilenkin’s tree-based coordinate systems; in potential theory they provide higher-order radial normalizers; and in contemporary computational science they furnish coordinate systems, operator factorizations, and product-manifold models adapted to high-dimensional inference and quantum dynamics (Kastrup, 2008). The term also appears in modern work on multicomplex mappings between spheres, where polyspherical coordinates are used to factor angle variables and construct a nearly invertible dimensional reduction from $S^{2^n-1}$ to $S^n$ (Livschitz et al., 2022).

1. Historical definition and canonical constructions

In the classical geometric literature, “polyspherical” and “hyperspherical” coordinates refer to the same class of angular–radial parameterizations, while in Darboux’s ambient formalism they also denote homogeneous coordinates such as tetracyclic, pentaspherical, and hexaspherical coordinates constrained by a quadratic form (Kastrup, 2008). In the Euclidean setting $\mathbb{R}^d$, standard hyperspherical coordinates are

$$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$

with

$$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$

where $0 \le r < \infty$, $0 \le \theta_i \le \pi$ for $i=1,\ldots,d-2$, and $0 \le \theta_{d-1} < 2\pi$ (Kastrup, 2008).

The same radial–angular logic reappears in Vilenkin’s polyspherical coordinates, where the angular variables are organized by a binary rooted tree. In that construction, the leaf nodes correspond to Cartesian components, and each traversal from root to leaf multiplies the radius by a sine or cosine factor depending on branch choice. The paper distinguishes four node types: type $a$ with azimuth $S^n$0, type $S^n$1 with $S^n$2, type $S^n$3 with $S^n$4, and type $S^n$5 with $S^n$6 (Cohl, 2012). Standard spherical coordinates arise as the chain type $S^n$7, whereas generalized Hopf coordinates arise on $S^n$8 through a recursive type-$S^n$9 construction (Cohl, 2012).

A distinct but related modern meaning appears in the LG fibration. There, polyspherical coordinates on $\mathbb{R}^d$0 are paired with multicomplex rotation groups to define a contraction $\mathbb{R}^d$1 and a projection $\mathbb{R}^d$2 (Livschitz et al., 2022). The general coordinate formula used there is

$\mathbb{R}^d$3

with one full-circle angle and the others taken on half-circles (Livschitz et al., 2022).

These constructions share a common structural feature: coordinates are introduced sequentially, and each new coordinate is weighted by products of sines and cosines inherited from earlier choices. This suggests that the “polyspherical approach” is less a single formalism than a transferable design principle for separating geometry into hierarchical angular layers.

2. Linearization of conformal geometry and compactification

In dimensions $\mathbb{R}^d$4, Liouville’s theorem implies that generic conformal maps are generated by translations, rotations, dilations, and inversion, forming a finite-dimensional Lie group of dimension $\mathbb{R}^d$5 (Kastrup, 2008). The nonlinearity enters through inversion

$\mathbb{R}^d$6

and through the resulting special conformal transformations obtained as inversion–translation–inversion:

$\mathbb{R}^d$7

with finite form

$\mathbb{R}^d$8

In the Minkowski-signature convention used in the cited paper, the denominator appears as $\mathbb{R}^d$9 (Kastrup, 2008).

Darboux’s polyspherical construction linearizes these fractional-linear actions. One introduces homogeneous coordinates $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$0 with

$$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$1

The quadratic cone $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$2 parametrizes the conformal compactification, and the conformal group acts linearly on the ambient variables. In particular, translations act as shears,

$$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$3

dilations rescale $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$4 and $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$5 oppositely, and special conformal transformations act as

$$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$6

(Kastrup, 2008). Projecting back by $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$7 recovers the nonlinear action in physical coordinates.

The same ambient perspective is tied to stereographic compactification. If $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$8 has coordinates $$x_1 = r \cos \theta_1,\quad x_2 = r \sin \theta_1 \cos \theta_2,\quad x_3 = r \sin \theta_1 \sin \theta_2 \cos \theta_3,\ \ldots$$9, stereographic projection from the north pole gives

$$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$0

with inverse

$$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$1

This compactifies $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$2 to $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$3 by adjoining a point at infinity (Kastrup, 2008). In Euclidean signature the conformal group is $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$4, and in Minkowski signature it is $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$5, both acting linearly on ambient coordinates (Kastrup, 2008).

In this setting, the polyspherical approach is an algebraic device: it replaces fractional-linear conformal action in $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$6-space by linear orthogonal action in a constrained $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$7-dimensional space. A plausible implication is that its enduring usefulness in field theory and AdS/CFT stems from this exact linearization property.

3. Separation of variables, harmonic analysis, and addition theorems

In $$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$8, hyperspherical coordinates induce the standard Jacobian

$$x_d = r \sin \theta_1 \cdots \sin \theta_{d-2} \sin \theta_{d-1},$$9

and the Laplacian decomposes as

$0 \le r < \infty$0

with spherical harmonics as eigenfunctions of $0 \le r < \infty$1 and eigenvalues $0 \le r < \infty$2 (Kastrup, 2008). Vilenkin’s tree-based polyspherical coordinates extend this separation mechanism beyond the standard chain by encoding hierarchical angular quantum numbers associated with the node types of the tree (Cohl, 2012).

This structure is central to the analysis of polyharmonic kernels. For the polyharmonic operator $0 \le r < \infty$3 on $0 \le r < \infty$4, a fundamental solution $0 \le r < \infty$5 satisfies

$0 \le r < \infty$6

and in the power-law regime one has kernels of the form $0 \le r < \infty$7, whereas in even dimensions with $0 \le r < \infty$8 logarithmic factors appear [(Cohl, 2012); (Cohl et al., 2022)]. The cited work develops Jacobi, Gegenbauer, and Chebyshev expansions for these kernels and compares Fourier expansions in rotationally invariant systems with Gegenbauer expansions in polyspherical coordinates to derive addition theorems (Cohl, 2012).

For $0 \le r < \infty$9 and $0 \le \theta_i \le \pi$0, the canonical Gegenbauer expansion is

$0 \le \theta_i \le \pi$1

where $0 \le \theta_i \le \pi$2 is the inter-point angle (Cohl, 2012). The associated Fourier expansion in a rotationally invariant coordinate system uses the toroidal parameter $0 \le \theta_i \le \pi$3 and Chebyshev polynomials or Fourier phases in the azimuthal difference (Cohl, 2012).

In even-dimensional polyharmonic theory, parameter differentiation of the master Gegenbauer expansion yields binomial and logarithmic kernels in Vilenkin polyspherical coordinates and in generalized Hopf coordinates (Cohl et al., 2022). The resulting addition theorems are stated in terms of normalized hyperspherical harmonics, products of angular factors $0 \le \theta_i \le \pi$4, and in the Hopf case meridional Jacobi blocks $0 \le \theta_i \le \pi$5 (Cohl et al., 2022). These formulas are explicitly designed for separation-of-variables solvers, multipole algorithms, and boundary integral methods (Cohl, 2012).

The polyspherical approach here is analytical rather than merely coordinatizing: it provides the angular hierarchy that makes the orthogonal polynomial structure visible. This suggests that its importance in higher-dimensional potential theory lies in furnishing coordinate systems whose separated factors are already aligned with the representation theory of $0 \le \theta_i \le \pi$6 and related subgroup chains.

4. Higher-order potential theory and polyspherical normalizers

On the Poincaré disk $0 \le \theta_i \le \pi$7 with hyperbolic metric

$0 \le \theta_i \le \pi$8

the hyperbolic Laplace–Beltrami operator is

$0 \le \theta_i \le \pi$9

and a $i=1,\ldots,d-2$0-polyharmonic function of order $i=1,\ldots,d-2$1 satisfies

$i=1,\ldots,d-2$2

(Picardello et al., 2023). The basic $i=1,\ldots,d-2$3-Poisson kernel is defined by

$i=1,\ldots,d-2$4

where $i=1,\ldots,d-2$5 is specified by $i=1,\ldots,d-2$6, and the classical Poisson kernel is

$i=1,\ldots,d-2$7

(Picardello et al., 2023).

The higher-order $i=1,\ldots,d-2$8-Poisson kernels are

$i=1,\ldots,d-2$9

while for $0 \le \theta_{d-1} < 2\pi$0 one has

$0 \le \theta_{d-1} < 2\pi$1

with $0 \le \theta_{d-1} < 2\pi$2 the Busemann function (Picardello et al., 2023). The $0 \le \theta_{d-1} < 2\pi$3-polyspherical functions are then defined by boundary averaging,

$0 \le \theta_{d-1} < 2\pi$4

with $0 \le \theta_{d-1} < 2\pi$5 equal to the classical $0 \le \theta_{d-1} < 2\pi$6-spherical function (Picardello et al., 2023).

Their role is asymptotic normalization. For $0 \le \theta_{d-1} < 2\pi$7 and $0 \le \theta_{d-1} < 2\pi$8,

$0 \le \theta_{d-1} < 2\pi$9

whereas for the critical value $a$0,

$a$1

(Picardello et al., 2023). For $a$2 not in the interior of the $a$3-spectrum, the zeros of these functions do not accumulate at the boundary, allowing normalized kernels

$a$4

to behave as approximate identities (Picardello et al., 2023).

An analogous normalization mechanism appears on regular trees. For the simple random walk operator $a$5 on the $a$6-regular tree, with spectral radius $a$7, $a$8-polyharmonicity of order $a$9 means

$S^n$00

For $S^n$01, the Martin kernel is

$S^n$02

where $S^n$03 is the small root of

$S^n$04

(Sava-Huss et al., 2019). The tree polyspherical functions are

$S^n$05

with $S^n$06 the spherical function and asymptotics

$S^n$07

(Sava-Huss et al., 2019).

In both settings, polyspherical functions are not coordinate charts but higher-order radial normalizers. They remove the exact growth generated by repeated application of the resolvent or Poisson kernel hierarchy. The recurrence of this pattern across the disk and the tree suggests a broader principle: higher-order boundary problems often require not only spherical functions but a full polyspherical tower adapted to the polyharmonic order.

5. Multicomplex and topological variants

In the LG fibration, the polyspherical approach is integrated with multicomplex algebra. The paper introduces commuting imaginary units

$S^n$08

their power-set basis, and the commutative ring

$S^n$09

(Livschitz et al., 2022). For $S^n$10, this gives the bicomplex algebra

$S^n$11

and the basic rotation group is

$S^n$12

More generally, the simple multicomplex rotation group is

$S^n$13

(Livschitz et al., 2022).

A core identity is the product decomposition

$S^n$14

proved inductively in the paper (Livschitz et al., 2022). This identity matches the polyspherical coordinate decomposition of $S^n$15.

The LG fibration is the composite $S^n$16, where the contraction map groups the $S^n$17 source angles of $S^n$18 into $S^n$19 multicomplex rotation parameters, and the projection $S^n$20 maps the multicomplex product to $S^n$21 using parity factors derived from modular reduction by $S^n$22 (Livschitz et al., 2022). The paper defines an exceptional set

$S^n$23

on which invertibility fails, and calls the map “almost bijective” because the kernel has measure zero (Livschitz et al., 2022).

The same work defines the multicomplex inner product

$S^n$24

for

$S^n$25

and introduces a distance difference function

$S^n$26

to determine when inner products are invariant under the dimensional reduction (Livschitz et al., 2022).

This is a nonclassical extension of the polyspherical idea. Instead of linearizing conformal action or separating the Laplacian, it organizes a hierarchy of commuting phase variables on $S^n$27 and projects them to $S^n$28. A plausible implication is that the common feature is again structural factorization: polyspherical layers provide a way to expose latent product geometry on spheres and sphere-like manifolds.

6. Contemporary computational and applied uses

In molecular reaction dynamics, the polyspherical approach is used to build curvilinear internal coordinates and a sum-of-products kinetic energy operator compatible with ML-MCTDH. The framework begins with a hierarchy of coordinate frames, from the space-fixed frame to internal and body-fixed fragment frames, with relative orientations specified by Euler angles $S^n$29 and fragment internals represented by Jacobi or Radau vectors (Zhang et al., 8 Jul 2025). In mass-weighted curvilinear coordinates $S^n$30, the kinetic energy operator is written in Laplace–Beltrami form,

$S^n$31

and in deepest-layer body-fixed frames it is decomposed into vibrational, Coriolis, and rotational contributions using tensors $S^n$32, $S^n$33, and $S^n$34 (Zhang et al., 8 Jul 2025). The same paper emphasizes that this naturally yields a sum-of-products Hamiltonian

$S^n$35

the precise structure required by ML-MCTDH (Zhang et al., 8 Jul 2025).

In fault-tolerant quantum computation for vibrational wave functions, the polyspherical approach denotes curvilinear internal coordinates built from bond lengths, bond angles, and dihedral angles, used particularly for large-amplitude motion such as torsion in $S^n$36 (Majland et al., 22 Aug 2025). There the general Laplace–Beltrami kinetic energy structure is assumed, but the practical operator is obtained from a low-order $S^n$37-mode expansion of the $S^n$38-matrix, specifically a $S^n$39-mode expansion producing an overall $S^n$40-body kinetic operator (Majland et al., 22 Aug 2025). The resulting Hamiltonian is encoded in a sum-over-products form, compressed by SVD or CP/HOOI, and then block-encoded for qubitization (Majland et al., 22 Aug 2025). The paper reports a hydrogen peroxide benchmark in polyspherical coordinates with Toffoli cost $S^n$41, qubit count $S^n$42, and $S^n$43 modals per mode (Majland et al., 22 Aug 2025).

In nonparametric statistics on product manifolds, the polyspherical approach refers to smoothing and inference on the polysphere

$S^n$44

For i.i.d. data $S^n$45 on this manifold, the kernel density estimator is

$S^n$46

with

$S^n$47

(García-Portugués et al., 2024). The theory derives pointwise bias and variance, asymptotic normality, AMISE-optimal bandwidths, product and spherically symmetric polyspherical kernels, and a Jensen–Shannon divergence $S^n$48-sample test (García-Portugués et al., 2024). An application embeds infant hippocampal morphology on $S^n$49 via skeletal representations and applies the smoothing methodology there (García-Portugués et al., 2024).

In geometric design, a related but specialized usage appears in polynomial approximation of $S^n$50 by spline patches on regular spherical triangulations. The approach uses canonical equilateral spherical triangles obtained by radial projection of tetrahedral, octahedral, or icosahedral faces and optimizes the simplified radial error

$S^n$51

for polynomial triangular patches (Vavpetič et al., 2021). Although this work does not use ambient conformal coordinates or harmonic-analysis trees, it still exploits symmetry reduction through elementary symmetric polynomials and spherical triangulations. This suggests that the label “polyspherical” can also mark a design methodology centered on decomposing spherical geometry into repeated canonical pieces.

7. Conceptual unification and scope

Across these literatures, several distinct meanings of “polyspherical approach” can be identified. In classical geometry and field theory, it is an ambient homogeneous-coordinate formalism that linearizes conformal transformations on a quadratic cone (Kastrup, 2008). In harmonic analysis and potential theory, it is a separation framework based on nested angular coordinates or higher-order spherical normalizers [(Cohl, 2012); (Picardello et al., 2023)]. In multicomplex topology, it is a layered angular factorization supporting a contraction–projection map between spheres of dimensions $S^n$52 and $S^n$53 (Livschitz et al., 2022). In computational molecular dynamics and quantum algorithms, it is a curvilinear internal-coordinate strategy that turns high-dimensional Hamiltonians into tractable sum-of-products forms (Zhang et al., 8 Jul 2025, Majland et al., 22 Aug 2025). In statistics, it denotes native inference on products of spheres rather than on a single sphere or in Euclidean ambient space (García-Portugués et al., 2024).

A concise comparison is useful.

Domain	Core object	Role of polyspherical structure
Conformal geometry	Ambient cone $S^n$54	Linearizes fractional-linear conformal action
Harmonic analysis	Tree-based angular coordinates	Enables separation, expansions, and addition theorems
Polyharmonic potential theory	Higher-order radial normalizers $S^n$55	Stabilizes boundary limits and Poisson representations
Multicomplex topology	$S^n$56	Factors sphere data into commuting rotation layers
ML-MCTDH / quantum algorithms	Curvilinear internal coordinates	Produces SOP kinetic and potential operators
Statistics on product spheres	$S^n$57	Supports geometry-aware KDE and inference

The principal misconception is that the term must refer to a single coordinate chart on $S^n$58. The literature shows at least three non-equivalent uses: standard hyperspherical coordinates, Darboux-type homogeneous ambient coordinates, and higher-order spherical normalizers. Another common simplification is to identify polyspherical with Hopf-type constructions only; in fact, the tree-based Vilenkin systems, the conformal ambient cone, and the product-sphere statistical framework are independent developments.

What unifies the topic is not a single formula but a recurrent strategy: represent geometry by iterated angular layers, by product sphere factors, or by an ambient homogeneous cone so that symmetry, separability, or operator factorization becomes explicit. This suggests why the polyspherical approach continues to reappear in settings as different as conformal symmetry, polyharmonic analysis, reaction dynamics, and high-dimensional manifold statistics.