---
title: Hermite Spectral Collocation Method
url: https://www.emergentmind.com/topics/hermite-spectral-collocation-method
type: topic
---

# Hermite Spectral Collocation Method

Hermite spectral collocation method denotes a class of global or mapped pseudospectral discretizations in which the approximation is built from Hermite polynomials or Hermite functions and the governing equation is enforced at Hermite nodes, typically the zeros of a Hermite polynomial. Its canonical setting is an unbounded domain, most often \(\mathbb{R}\), where the Gaussian structure of the Hermite basis encodes decay at infinity and makes the method especially suitable for Schrödinger-type operators, Fokker–Planck-type equations, and other problems with rapidly decaying solutions [2507.15350].

## 1. Basis functions, approximation spaces, and node sets

The algebraic foundation is the family of Hermite polynomials
\[
H_n(x)=(-1)^n e^{x^2}\frac{d^n}{dx^n}e^{-x^2},
\]
which satisfy the weighted orthogonality
\[
\int_{-\infty}^{\infty} H_n(x)H_m(x)e^{-x^2}\,dx=\gamma_n\delta_{n,m},\qquad \gamma_n=2^n n!\sqrt{\pi}.
\]
A closely related orthonormal basis on \(L^2(\mathbb{R})\) is furnished by Hermite functions, obtained by absorbing the Gaussian weight into the basis itself [2507.15350].

For collocation and interpolation on \(\mathbb{R}\), one works in finite-dimensional Hermite spaces such as
\[
\mathbb{H}_n=\operatorname{span}\{\psi_0,\psi_1,\dots,\psi_n\},
\]
or, in alternate notation,
\[
\mathcal{H}_N=\operatorname{span}\{\widetilde H_0,\widetilde H_1,\dots,\widetilde H_N\}.
\]
The nodes are chosen as the real zeros of \(H_{n+1}\), equivalently of \(\psi_{n+1}\), and inherit the symmetry \(x_j=-x_{n-j}\) [2507.15350]. In the collocation interpretation used for kinetic problems, the velocity variable is represented by Lagrange cardinal functions built on the zeros of \(H_M\), so that differentiation and quadrature are both tied to Gauss–Hermite structure [1807.02418].

Two Hermite-type bases are prominent in unbounded-domain collocation for nonlocal operators: the normalized Hermite functions \(\widehat H_n\) and the over-scaled functions \(\widetilde H_n=e^{-x^2/2}\widehat H_n\). Both are paired with Gauss–Hermite collocation points, but they lead to markedly different conditioning properties [1801.09073].

## 2. Collocation formulation and differentiation structure

In its simplest global form, Hermite spectral interpolation seeks \(h_n\in\mathbb{H}_n\) such that
\[
h_n(x_j)=f(x_j),\qquad j=0,\dots,n,
\]
at the zeros of \(H_{n+1}\). The same nodal set underlies spectral collocation for differential equations: one approximates the solution by a Hermite expansion and enforces the differential operator pointwise at the Hermite nodes [2507.15350].

For second-order model problems on \(\mathbb{R}\), two standard collocation formulations are
\[
u_n''(x_j)+(\alpha-x_j^2)u_n(x_j)=f(x_j),\qquad j=0,\dots,n,
\]
and
\[
-u_n''(x_j)+\alpha u_n(x_j)=f(x_j),\qquad j=0,\dots,n.
\]
Using cardinal functions
\[
\sigma_j(x)=\frac{\psi_{n+1}(x)}{\psi_{n+1}'(x_j)(x-x_j)},
\]
the collocation ansatz yields matrix systems of the form
\[
(D+S)\hat u_n=\hat f_n
\]
or
\[
(-D+\alpha I)\hat u_n=\hat f_n,
\]
where \(D=[\sigma_j''(x_i)]_{i,j=0}^n\) is the Hermite second-derivative matrix and \(S\) is the diagonal contribution from the potential term [2507.15350].

For fractional PDEs on \(\mathbb{R}^d\), the principal technical task is the construction of differentiation matrices for \((-\Delta)^{\alpha/2}\). In that setting, the method consists of expanding the solution with proper global basis functions and imposing collocation conditions on the Gauss-Hermite points; the corresponding differentiation matrices are computed recursively, using explicit formulas involving confluent hypergeometric functions and scaling relations for the fractional Laplacian [1801.09073].

A distinct but related collocation realization appears in semi-Lagrangian Vlasov–Poisson discretization. There, the velocity dependence is represented by Lagrange cardinal functions at Gauss–Hermite nodes, derivatives in \(v\) are produced by Hermite differentiation matrices, and velocity integrals are approximated by Gauss–Hermite quadrature [1807.02418]. This retains the collocational character but applies it to phase-space transport rather than directly to a spatial boundary-value problem.

## 3. Convergence, error bounds, and superconvergence

For analytic data on \(\mathbb{R}\), Hermite spectral interpolation exhibits exponential convergence. Under analyticity in the strip
\[
\mathcal S_\rho=\{z\in\mathbb C:\Im z\in[-\rho,\rho]\},
\]
Wang and Zhang derive a global max-norm error estimate with decay factor \(e^{-\rho\sqrt{2n}}\); derivative errors carry polynomial prefactors, roughly \(n^{m/2}\) for the \(m\)-th derivative [2507.15350].

The central refinement is superconvergence. For Hermite interpolation in \(\mathbb{H}_n\), the first derivative superconverges at the zeros \(\{\tau_j\}\) of \(\psi_{n+1}'\), while the second derivative superconverges at the zeros \(\{\eta_j\}\) of \(\psi_{n+1}''\). Globally, the first-derivative error at \(\tau_j\) improves from \(O(n^{1/2}e^{-\rho\sqrt{2n}})\) to at least \(O(n^{1/6}e^{-\rho\sqrt{2n}})\), and in the central region \(|x|\le \varrho\sqrt{2n+1}\) to \(O(e^{-\rho\sqrt{2n}})\). For the second derivative, the error at \(\eta_j\) improves from \(O(ne^{-\rho\sqrt{2n}})\) to \(O(n^{1/2}e^{-\rho\sqrt{2n}})\) [2507.15350].

The same pattern carries into collocation for differential equations. If the exact solution belongs to \(\mathbb{H}_{n+1}\), then the collocation error has the exact form
\[
u(x)-u_n(x)=a_{n+1}\psi_{n+1}(x),
\]
so the function error vanishes at the collocation nodes, the derivative error superconverges at the zeros of \(\psi_{n+1}'\), and the second-derivative error superconverges at the zeros of \(\psi_{n+1}''\) [2507.15350]. An important contrast with Chebyshev interpolation is that, for Hermite interpolation, the number of superconvergence points increases with derivative order [2507.15350].

For unbounded-domain fractional PDEs, the dominant practical issue is not only approximation error but conditioning. Numerical evidence shows that the differentiation matrices associated with the normalized Hermite basis have algebraically growing condition numbers, whereas the over-scaled basis leads to much worse conditioning; scaling factors strongly influence accuracy and convergence speed [1801.09073].

## 4. Unbounded and semi-infinite domains

The native domain of the Hermite basis is \(\mathbb{R}\), and several formulations work directly there. In Wang and Zhang’s analysis, no mapping from a finite interval is used; the Gaussian factor in the Hermite functions provides the decay needed to approximate solutions vanishing at infinity [2507.15350]. This direct formulation is particularly natural when the operator itself contains \(x^2\), as in harmonic-oscillator or Schrödinger-type problems [2507.15350].

For semi-infinite domains \((0,\infty)\), two mapping strategies recur. In the Thomas–Fermi setting, Hermite functions on \(\mathbb{R}\) are composed with
\[
\phi(z)=\frac{1}{k}\ln(z),\qquad z\in(0,\infty),
\]
so that transformed basis functions \(\hat H_n(z)=\widetilde H_n(\phi(z))\) can be collocated at mapped Hermite nodes \(x_j=\phi^{-1}(w_j)=e^{kw_j}\) [1604.01454]. In Lane–Emden and Volterra-type problems, a logarithmic–hyperbolic map
\[
\phi(z)=\ln(\sinh(kz))
\]
transfers \((0,\infty)\) to \(\mathbb{R}\), after which transformed Hermite functions are used in a collocation framework with boundary behavior built into the ansatz [1008.2063, 1008.2337].

These mapped formulations address two recurrent obstacles of semi-infinite problems: singular behavior at the origin and boundary conditions at infinity. In the Lane–Emden formulation, the approximate solution is modified to incorporate \(y(0)=A\) and \(y'(0)=B\) directly through
\[
\widehat\xi_N y(x)=A+Bx+x\,\widehat\xi_N y(x),
\]
so that the transformed Hermite part does not have to enforce derivative data at the singular point itself [1008.2063]. In the Thomas–Fermi formulation, a lifting function
\[
p(x)=1+\lambda x+x^2
\]
is combined with the transformed Hermite expansion, and the residual is set to zero at transformed Hermite collocation nodes, producing a nonlinear algebraic system solved by Newton’s method [1604.01454].

## 5. High-dimensional, fractional, and sparse formulations

In high dimensions, the principal obstacle is the tensor-product growth of the Hermite basis. For generalized Hermite functions
\[
\Phi_{\boldsymbol n}^{\boldsymbol\alpha,\boldsymbol\beta}(\boldsymbol x)=\prod_{j=1}^d \phi_{n_j}^{\alpha_j,\beta_j}(x_j),
\]
full grids lead to \((N+1)^d\) degrees of freedom, whereas regular hyperbolic cross and optimized hyperbolic cross index sets reduce the approximation space dramatically [1306.3207]. The regular hyperbolic cross has cardinality \(\mathcal O(N(\log N)^{d-1})\), while for optimized hyperbolic cross with \(\gamma\in(0,1)\) the cardinality is \(\mathcal O(N)\) [1306.3207].

Although the high-dimensional parabolic theory is presented in Galerkin form, the same generalized Hermite basis, hyperbolic cross index sets, and sparse Smolyak Gauss–Hermite grids directly inform collocation design. The paper’s approximation estimates show that, in weighted Korobov spaces, optimized hyperbolic cross Hermite approximations can retain dimension-independent convergence exponents, with dimensional dependence entering only through constants [1306.3207]. This suggests that sparse-grid Hermite collocation is the natural high-dimensional analogue of the one-dimensional Gauss–Hermite collocation method.

For space-fractional PDEs on \(\mathbb{R}^d\), Hermite collocation replaces local derivative matrices by differentiation matrices for \((-\Delta)^{\alpha/2}\). Two Hermite-type bases are employed, the collocation points are Gauss–Hermite points, and the matrices are computed recursively; multi-term fractional PDEs are handled by summing the matrices for the individual fractional orders [1801.09073].

A further extension appears in kinetic equations, where Hermite collocation is applied in velocity space rather than in physical space. In the semi-Lagrangian Vlasov–Poisson method, the unbounded velocity domain is treated with Gauss–Hermite nodes and Hermite differentiation matrices, and the numerical behavior depends strongly on a scaling parameter in the Gaussian weight, indicating that adaptive scaling is a central unresolved issue for Hermite-based solvers on unbounded velocity domains [1807.02418].

## 6. Applications, terminology, and limitations

Hermite spectral collocation is used for unbounded-domain ODEs and PDEs, including Schrödinger-type operators, fractional PDEs, Lane–Emden equations, Thomas–Fermi equations, Volterra population models, and kinetic equations in velocity space [2507.15350, 1801.09073, 1008.2063, 1604.01454, 1008.2337, 1807.02418]. Across these settings, three implementation themes recur: collocation at Hermite or transformed Hermite nodes, basis scaling or mapping to match decay, and the reduction of the continuous problem to linear or nonlinear algebraic systems.

A common misconception is that every “Hermite collocation” method is a global spectral method on \(\mathbb{R}^d\). In the Maxwell literature, however, “Hermite methods” can denote local tensor-product Hermite–Birkhoff schemes on Cartesian grids. In that formulation there is no global Hermite basis, no Hermite transform, and no Hermite quadrature; instead, one stores derivative data at grid vertices, reconstructs a local polynomial on each cell, and advances the solution explicitly with an energy-conserving update [2401.12043]. The relation to Hermite spectral collocation is structural rather than literal: both use Hermite interpolation data, but one is global and modal or nodal on Gauss–Hermite points, whereas the other is local and cell-based [2401.12043].

The main limitations are also consistent across the literature. Conditioning can deteriorate rapidly, especially for over-scaled bases in fractional problems [1801.09073]. On semi-infinite domains, performance depends on mapping and scaling parameters, which are usually tuned empirically [1604.01454, 1008.2063]. In high dimensions, full tensor Hermite spaces remain prohibitive without hyperbolic-cross or sparse-grid compression [1306.3207]. For collocation in velocity space, a poor choice of Gaussian scaling can substantially degrade accuracy or even destabilize long-time computations, so adaptive scaling remains an important practical problem [1807.02418].

Taken together, these developments define Hermite spectral collocation method as a family of Gaussian-adapted collocation techniques for unbounded or mapped-unbounded geometries, distinguished by Hermite node sets, sparse derivative structure in Hermite coordinates, and, in the most recent work, a detailed superconvergence theory that identifies where collocation and interpolation are most accurate [2507.15350].

Source: https://www.emergentmind.com/topics/hermite-spectral-collocation-method