---
title: Hybrid Spherical-Plane Wave Model (HSPWM)
url: https://www.emergentmind.com/topics/hybrid-spherical-plane-wave-model-hspwm
type: topic
---

# Hybrid Spherical-Plane Wave Model (HSPWM)

Hybrid Spherical-Plane Wave Model (HSPWM) denotes a class of hybrid representations that combine spherical-wave and plane-wave descriptions when neither a purely spherical-wave model (SWM) nor a purely plane-wave model (PWM) is adequate over the full domain of interest. In extremely large-scale array communications, HSPWM is constructed from explicit PWM/SWM demarcations based on channel gain and effective rank, and it selects SWM or PWM according to geometry-dependent thresholds [2301.06036]. In Terahertz ultra-massive MIMO, a closely related formulation, often written as HSPM, uses PWM within subarrays and SWM among subarrays [2106.05491]. The same label, or a structurally analogous hybridization, also appears in integrated UM-MIMO–IRS channels, modular XL-array localization, distributed multi-UAV near-field communications, active-source wavefield modeling, Gaussian-basis scattering theory, and hybrid-basis galaxy clustering inference [2205.13113], [2504.13455], [2509.06967], [2402.00573], [2605.18564], [2007.14962].

## 1. Terminology and scope

The literature does not use HSPWM as a single standardized construction. Instead, the recurring idea is to combine a spherical description where curvature, range variation, or global geometry matters, with a plane-wave description where local structure or computational tractability makes planarization valid. This suggests that HSPWM is best understood as a modeling pattern rather than a unique formula.

| Domain | Hybridization principle | Representative paper |
|---|---|---|
| XL-array communications | SWM inside demarcated near-field region, PWM outside | [2301.06036] |
| THz UM-MIMO | PWM within subarray, SWM among subarrays | [2106.05491] |
| UM-MIMO–IRS cascaded channels | Subarray-wise planar steering with spherical inter-subarray phases | [2205.13113] |
| Modular XL-array localization | Planar steering within each SA, spherical variation across SAs | [2504.13455] |
| Multi-UAV near-field communications | Kronecker form: spherical across UAVs, planar within each UPA | [2509.06967] |
| Layered half-space wavefields | Cylindrical Hankel source field with plane-wave layer solution | [2402.00573] |
| Quantum scattering / chemistry | SGTOs combined with plane-wave-modulated SGTOs | [2605.18564] |
| Galaxy clustering | Spherical Fourier modes at low $k$, plane-wave FFT block at high $k$ | [2007.14962] |

Within this broad family, the most explicit decision-theoretic formulation is the XL-array communication model derived from PWM/SWM applicability boundaries. Other papers instead implement hybridization through subarray partitioning, basis design, or asymptotic matching.

## 2. Demarcation-based HSPWM in extremely large-scale arrays

For XL-array communications, the underlying question is not whether SWM or PWM is universally preferable, but where each is applicable. The baseline PWM assumes that the transmitter–receiver distance is sufficiently large that the wavefront is locally planar across the array, so path gain depends on a common distance $r$ and phase slope is determined by an incident angle. For a ULA, the far-field steering vector is
$$
\bar h_n=\frac{\lambda}{4\pi r}\exp\!\Bigl(-j\frac{2\pi}{\lambda}\bigl(r-d(n-\tfrac{N-1}{2})\sin\theta\bigr)\Bigr).
$$
The SWM removes the planar approximation and assigns a distinct propagation distance $r_n$ to each element:
$$
h_n=\frac{\lambda}{4\pi r_n}\exp\!\bigl(-j\frac{2\pi}{\lambda}r_n\bigr),
$$
with
$$
r_n=\sqrt{r^2-2rd\,(n-\tfrac{N-1}{2})\sin\theta+\bigl(d(n-\tfrac{N-1}{2})\bigr)^2 }.
$$
These definitions are treated for four single-LoS setups: point-to-ULA, point-to-UPA, ULA-to-ULA, and ULA-to-UPA [2301.06036].

For point-to-ULA, the normalized received-power ratio under SWM versus PWM with MRC is
$$
\mu(r,\theta)=\frac{G_{\rm near}(r,\theta)}{G_{\rm far}(r,\theta)}
=\frac{r^2}{N}\sum_{n=1}^N\frac{1}{r_n^2}.
$$
A closed form is
$$
\mu(r,\theta)
=\frac{r}{N d\cos\theta}\Biggl[
\arctan\!\Bigl(\frac{N d}{2r\cos\theta}+\tan\theta\Bigr)
+\arctan\!\Bigl(\frac{N d}{2r\cos\theta}-\tan\theta\Bigr)
\Biggr].
$$
Analysis of $\partial^2\mu/\partial r^2$ yields an inflection point at $\pm \pi/6$. For $|\theta|<\pi/6$, $\mu(r,\theta)<1$ always, so PWM uniformly over-estimates SWM power; for $|\theta|>\pi/6$, $\mu(r,\theta)$ first exceeds $1$ and then returns to $1$ [2301.06036].

For point-to-UPA, the analogous quantity $\mu(r,\phi,\varphi)$ produces an equi-power surface. In the uniform circular PA case with $\beta\equiv \cos^2\phi\cos^2\varphi$, the dividing curve is
$$
\cos^2\phi\cos^2\varphi=\frac{1}{2}.
$$
When $\beta\ge \tfrac12$, PWM always over-estimates, whereas for $\beta<\tfrac12$, $\mu(r)$ first exceeds $1$ and then converges to $1$.

For MIMO, the paper adopts the effective-rank metric
$$
\operatorname{erank}(A)=\exp\!\left[-\sum_i \frac{\sigma_i}{\|A\|_*}\ln\!\left(\frac{\sigma_i}{\|A\|_*}\right)\right]\in[1,\min(p,q)].
$$
Under pure PWM LoS MIMO, $\operatorname{erank}=1$; under SWM, $\operatorname{erank}>1$ in the near field. For ULA-to-ULA, the threshold surface is approximated by
$$
r(\theta,\phi)\approx r_1\Bigl|
\cos^2\bigl(\theta+\tfrac{\phi}{2}\bigr)-\sin^2\bigl(\tfrac{\phi}{2}\bigr)
\Bigr|,
$$
and for ULA-to-UPA an upper bound is
$$
r(\phi,\varphi)\le
r_1-r_1\bigl(1-|\sin\phi|\bigr)\bigl(1-\cos^2\varphi\bigr).
$$
The same study also notes that scatterers cause $\mu(r)$ to fluctuate and extend the near-field region when scatterer power rises [2301.06036].

These results lead directly to a regime-switching HSPWM. For point-to-ULA, one uses SWM if $r<r_p(\theta)$, where $r_p(\theta)$ solves $\mu(r,\theta)=\delta_\Delta\approx1$, and PWM otherwise. In indicator form,
$$
h_n^{\rm HSPW}
=
I_{\rm SW}(r,\theta)\,\frac{\lambda}{4\pi r_n}e^{-j2\pi r_n/\lambda}
+
\bigl(1-I_{\rm SW}(r,\theta)\bigr)\,
\frac{\lambda}{4\pi r}e^{-j2\pi(r-d(n-\frac{N-1}{2})\sin\theta)/\lambda}.
$$
For other geometries, $r_n$ is replaced by the appropriate elementwise distance.

## 3. Subarray-based channel modeling in THz UM-MIMO and IRS systems

A second major interpretation of HSPWM appears in THz ultra-massive MIMO, where the hybridization is spatially hierarchical rather than boundary-based. The key assumption is that a subarray aperture is small enough for a plane-wave approximation to hold within the subarray, while the distances and angles among subarray reference antennas require spherical modeling. In this formulation, the $k_r\times k_t$ block of path $p$ is
$$
\mathbf H_p^{\,k_r k_t}
=
\alpha_p^{\,k_r k_t}\,
\mathbf a_{r,p}^{\,k_r k_t}
\bigl(\mathbf a_{t,p}^{\,k_r k_t}\bigr)^H,
$$
where $\alpha_p^{\,k_r k_t}$ retains exact distance-dependent phase and amplitude, and $\mathbf a_{r,p}^{\,k_r k_t}$ and $\mathbf a_{t,p}^{\,k_r k_t}$ are planar steering vectors within the corresponding subarrays [2106.05491].

This model was evaluated against full SWM by the normalized Frobenius errors
$$
\varepsilon_P=\frac{\|\mathbf H_P-\mathbf H_S\|_F}{\|\mathbf H_S\|_F},
\qquad
\varepsilon_{HSPM}=\frac{\|\mathbf H_{\rm HSPM}-\mathbf H_S\|_F}{\|\mathbf H_S\|_F}.
$$
For $f=0.4$ THz, subarray spacing $32\lambda$, and $N_t=N_r=1024$, $\varepsilon_{HSPM}$ is up to $14$ dB lower than $\varepsilon_P$ at $20$ m. The parameter count is also intermediate: SWM requires $2N_pN_tN_r$ real parameters, PWM requires $6N_p$, and HSPM requires
$$
N_p\bigl[1+5K_tK_r\bigr].
$$
The same paper develops a two-phase channel-estimation mechanism using a DCNN for reference subarrays and geometric extrapolation to the remaining subarray pairs [2106.05491].

The integrated UM-MIMO–IRS case extends the same idea to cascaded channels. The near-/far-field partition is expressed through the Rayleigh distance
$$
D_{\rm Rayleigh}=\frac{2S^2}{\lambda}.
$$
For the HSPWM channel matrix,
$$
\min\{K_rN_p,\;K_tN_p,\;N_r,\;N_t\}
\le
\rank(H_{\rm HSPM})
\le
\min\{K_rK_tN_p,\;N_r,\;N_t\}.
$$
For the cascaded channel $H^{\rm cas}=A\,\diag(\mathbf p)\,B$, the rank obeys
$$
\rank(H^{\rm cas})\le \min\{\rank(A),\rank(B)\}.
$$
Numerically, at $D=40$ m and $N=256$, the capacity of HSPWM deviates from full SWM by only $5\times10^{-4}$ bits/s/Hz, whereas PWM is approximately $37$ bits/s/Hz worse. The same work introduces a subarray-based sparse representation and two compressive-sensing estimation algorithms, separate-side estimation (SSE) and dictionary-shrinkage estimation (DSE); DSE is reported to be about $0.8$ dB better than SSE at low SNR [2205.13113].

## 4. Localization, compressed sensing, and tensor formulations

In modular XL-array localization, HSPWM is tied to spatial non-stationarity (SNS). A base station with $K$ subarrays uses planar steering within each subarray and spherical variation of $(\theta_{k,\ell,p},\phi_{k,\ell,p},d_{k,\ell,p})$ across subarrays. Visibility regions are modeled by binary indicators $\chi_{k,\ell,p}\in\{0,1\}$. The resulting localization pipeline has three stages: visible-SA selection and AoA estimation via SOMP, coarse 3-D position estimation via weighted least squares (WLS), and reduced-dictionary CS followed by a final WLS refinement. The paper reports that SA interval $D$ must balance angular diversity and link quality, with a simulation optimum of about $1$ m in a $5\times5$ grid at THz $320$ GHz; with total $M=576$ elements, the reported optimum is $K_x=K_z=4$, $M_x=M_z=6$ [2504.13455].

In distributed multi-UAV near-field communications, HSPWM is also called the cross-field model. The array response is approximated by a Kronecker product
$$
\mathbf g_{\rm HSPWM}(\theta,\varphi,r)
\approx
\mathbf b(\theta,\varphi,r)\otimes \mathbf a(\theta,\varphi),
$$
where $\mathbf a$ captures the planar manifold within each UAV subarray and $\mathbf b$ captures spherical variation across UAVs. Under MRC, the SNR becomes
$$
\gamma_{\rm HSPWM}
=
\bar P\,\|\mathbf a\|^2\,\|\mathbf b\|^2,
$$
which reduces to a single double sum over UAV indices rather than the four-fold sum of full SWM. The paper states that $\gamma_{\rm HSPWM}$ nearly coincides with $\gamma_{\rm SWM}$ over all tested angles and spacings, while retaining lower analytical complexity. This structure yields a rank-1 tensor representation for each path and motivates tensor-OMP; simulation results show tensor-OMP achieves NMSE comparable to spherical-domain OMP (SD-OMP), with reduced computational complexity and improved scalability [2509.06967].

These algorithmic developments show that HSPWM is not only a forward model. It also defines the geometry of the inference problem: support sets in CS, visible-region selection, reduced dictionaries, and low-rank tensor atoms all inherit the specific way spherical and planar components are separated.

## 5. Extensions beyond wireless communications

The same hybrid logic appears in several other fields, although the mathematical objects differ substantially.

In active-source wavefield modeling for a layered half-space, the model replaces the usual planar-wave assumption of free-vibration forward models with a cylindrically spreading source field described by Hankel functions, while retaining a plane-wave-based layered eigenproblem. The vertical and radial surface responses are written as modal sums of $H_0^{(2)}(k_m r)$ and $H_1^{(2)}(k_m r)$ terms, including both propagating and decaying modes. Reported runtimes are at least two orders of magnitude faster than numerical methods; for example, Profile I gives HSPWM $\sim 5$ s versus SGFD $\sim 2200$ s, and Profile II gives HSPWM $\sim 0.8$ s versus SGFD $\sim 7360$ s. The method also captures modal osculation and leaky waves [2402.00573].

In quantum physics and chemistry, HSPWM denotes an analytic framework for free-particle Green’s-function matrix elements over spherical Gaussian-type orbitals (SGTOs) and plane-wave-modulated SGTOs. Plane-wave modulation shifts Gaussian centers into the complex domain, with $\mathbf C^\dagger=\mathbf C+\tfrac{i\mathbf k}{2\alpha}$ for the $s$-type case, and the general PW-SGTO matrix elements reduce to SGTO core integrals evaluated at the complex separation vector
$$
\mathbf R^\dagger
=
(\mathbf A-\mathbf B)-i\Bigl(\frac{\mathbf k_1}{2\alpha}+\frac{\mathbf k_2}{2\beta}\Bigr).
$$
The framework provides recurrence relations for two-center radial integrals and asymptotic expansions for large $R/\sqrt\eta$ and large $\sqrt\eta\,k_0$, with the stated purpose of stable continuum-electron calculations in scattering and autoionization studies [2605.18564].

In galaxy clustering inference, the hybridization occurs in spectral space rather than physical space. The density contrast is decomposed into a spherical Fourier–Bessel block for $k\le k_{\rm hyb}$ and a Cartesian plane-wave FFT block for $k>k_{\rm hyb}$. The implementation adopts a sharp cutoff at $k_{\rm hyb}\approx0.04\,h\,{\rm Mpc}^{-1}$ with $\ell_{\max}=15$ for $R\approx 500\,h^{-1}\,{\rm Mpc}$, giving $N_d=456$ SFB modes. The total likelihood is taken as
$$
L_{\rm hyb}(\theta)=L_{\rm sph}(\theta)\,L_{\rm pw}(\theta),
$$
under the approximation that the two blocks are weakly correlated; mock-based cross-correlations are reported to be smaller than $0.06$ [2007.14962].

Related mathematical constructions make the same spherical/plane-wave duality explicit. One work shows how exact spherical electromagnetic and Robinson–Trautman gravitational waves approach plane-fronted limits and then constructs a hybrid metric
$$
g_{\mu\nu}^{\rm HSPWM}
=
H(r)\,g_{\mu\nu}^{\rm spherical}
+
[1-H(r)]\,g_{\mu\nu}^{\rm plane},
$$
with matching in a buffer zone [2109.04158]. Another proves that any Helmholtz solution in a ball can be represented as a continuous superposition of evanescent plane waves, and numerical tests report machine-precision accuracy up to $\ell\lesssim 4\kappa$ with bounded coefficients for EPW approximation sets [2305.02175]. A further mathematical line expresses reproducing kernels of spherical, complex, and symplectic harmonics as plane-wave integrals over Stiefel manifolds [1702.00611].

## 6. Computational profile, limitations, and recurrent misconceptions

In the XL-array formulation derived from PWM/SWM demarcations, the computational contrast is explicit. SWM-based channel generation costs $O(\#\text{elements})$ distance-square-root and exponential evaluations per user, whereas PWM uses a precomputed steering vector and a common factor $1/r$. The stated implementation strategy is to precompute equi-power or equi-rank boundary surfaces $r_{\rm threshold}(\text{angles})$ offline, then, at runtime, compute $(r,\theta)$ or $(r,\phi,\varphi)$, look up $r_{\rm th}$, select SWM or PWM, and form the channel vector or matrix accordingly. The threshold parameter is chosen in practice as $\delta_\Delta\approx 0.99$–$1.05$ [2301.06036].

Across fields, the computational motivation is similar but not identical. In layered half-space modeling, the hybrid model is reported to be hundreds to thousands of times faster than SGFD or DSGFD on representative profiles [2402.00573]. In the multi-UAV setting, tensor-OMP exploits the Kronecker structure induced by HSPWM to reduce projection cost relative to SD-OMP [2509.06967]. In the cosmology hybrid basis, a single evaluation of the compressed spherical likelihood takes about $1$ min, while the FFT-based plane-wave step is sub-second per likelihood call [2007.14962].

Several misconceptions recur. First, HSPWM is not a universally fixed acronym with a single canonical formula; the literature uses the label for non-equivalent constructions in communications, wave physics, chemistry, cosmology, and mathematical analysis. Second, HSPWM is not always a literal interpolation between a spherical formula and a planar formula. In wireless channels it is often a regime selector or a subarray decomposition; in galaxy clustering it is a split in $k$-space; in Gaussian-basis scattering it is a hybrid basis with complex-shifted centers. Third, the hybrid model is not a claim that full SWM is unnecessary. The XL-array results show that applicability depends on angle, array structure, effective rank, and scatterers, while the THz and distributed-array papers tie hybridization to subarray aperture and spacing rather than to a single universal near-/far-field boundary [2301.06036], [2106.05491], [2509.06967].

Taken together, these works define HSPWM as a family of technically distinct but structurally related methods: each preserves spherical behavior where curvature, range dependence, or global geometry is decisive, and each retains plane-wave structure where local approximation or algorithmic efficiency is decisive.

Source: https://www.emergentmind.com/topics/hybrid-spherical-plane-wave-model-hspwm