Volume 55, Number 5, September-October 2021
|Page(s)||2141 - 2168|
|Published online||13 October 2021|
Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
School of Mathematics, Shanghai University of Finance and Economics, Shanghai
2 College of Science, Nanjing University of Posts and Telecommunications, Nanjing 210023, China
3 State Key Laboratory of Computer Science/Laboratory of Parallel Computing, Institute of Software, Chinese Academy of Sciences, Beijing 100190, China
4 Division of Mathematical Sciences, School of Physical and Mathematical Sciences, Nanyang Technological University, Nanyang 637371, Singapore
5 Beijing Computational Science Research Center, Beijing, 100193, China
* Corresponding author: firstname.lastname@example.org
Accepted: 20 August 2021
In this paper, we introduce two families of nontensorial generalised Hermite polynomials/functions (GHPs/GHFs) in arbitrary dimensions, and develop efficient and accurate spectral methods for solving PDEs with integral fractional Laplacian (IFL) and/or Schrödinger operators in ℝ d . As a generalisation of the G. Szegö’s family in 1D (1939), the first family of multivariate GHPs (resp. GHFs) are orthogonal with respect to the weight function |x|2μe-|x|2 (resp. |x|2μ ) in ℝd. We further construct the adjoint generalised Hermite functions (A-GHFs), which have an interwoven connection with the corresponding GHFs through the Fourier transform, and are orthogonal with respect to the inner product [u,v]Hs(ℝd) = ((-Δ)s/2u, (-Δ)s/2v)ℝd associated with the IFL of order s > 0. As an immediate consequence, 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 remarkably efficient in solving PDEs with the fractional Schrödinger operator: (-Δ) s + | x |2μ with s ∈ (0,1] and μ > −1/2 in ℝ d We construct the second family of multivariate nontensorial Müntz-type 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 solution to some Schrödinger eigenvalue problems.
Mathematics Subject Classification: 65N35 / 65N25 / 35Q40 / 33C45 / 65M70
Key words: Generalised Hermite polynomials/functions / integral fractional Laplacian / Schrödinger operators with fractional spower potential / Müntz-type generalised Hermite functions
© The authors. Published by EDP Sciences, SMAI 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.