2020/02/13 by Sheng, Changtao, Ma, Suna, Li, Huiyuan +2
#33C45 #35Q40 #65M70 #65N25 #65N35 #FOS: Mathematics #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.2002.05334
In this paper, we introduce two new families of generalised Hermite polynomials/functions (GHPs/GHFs) in arbitrary dimensions, and develop efficient and accurate generalised Hermite spectral algorithms for PDEs with integral fractional Laplacian (IFL) and/or Schrödinger operators in \mathbb Rd. As a generalisation of the G. Szegö's family in 1D (1939), the first family of GHPs (resp. GHFs) are orthogonal with respect to |\bx|2μ \e-|\bx|2 (resp. |\bx |2μ) in \mathbb Rd. We further define adjoint generalised Hermite functions (A-GHFs) which have an interwoven connection with the corresponding GHFs through the Fourier transform, and which are orthogonal with respect to the inner product [u,v]Hs(\mathbb Rd)=((-Δ)s/ 2u, (-Δ)s/2 v )\mathbb Rd associated with the IFL of order s>0. Thus, the spectral-Galerkin method using A-GHFs as basis functions leads to a diagonal stiffness matrix for the IFL (which is known to be notoriously difficult and expensive to discretise). The new basis also finds efficient and accurate in solving PDEs with the fractional Schrödinger operator: (-Δ)s +|\bs x|2μ with s∈ (0,1] and μ>-1/2. Following the same spirit, we construct the second family of GHFs, dubbed as Müntz-type generalised Hermite functions (M-GHFs), which are orthogonal with respect to an inner product associated with the underlying Schrödinger operator, and are tailored to the singularity of the solution at the origin. We demonstrate that the Müntz-type GHF spectral method leads to sparse matrices and spectrally accurate to some Schrödinger eigenvalue problems.