Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators

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...

Full description

Saved in:
Bibliographic Details
Main Authors: Sheng, Changtao, Ma, Suna, Li, Huiyuan, Wang, Li-Lian, Jia, Lueling
Other Authors: School of Physical and Mathematical Sciences
Format: Article
Language:English
Published: 2022
Subjects:
Online Access:https://hdl.handle.net/10356/159739
Tags: Add Tag
No Tags, Be the first to tag this record!
Institution: Nanyang Technological University
Language: English
id sg-ntu-dr.10356-159739
record_format dspace
spelling sg-ntu-dr.10356-1597392023-02-28T20:02:20Z Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators Sheng, Changtao Ma, Suna Li, Huiyuan Wang, Li-Lian Jia, Lueling School of Physical and Mathematical Sciences Science::Mathematics Generalised Hermite Polynomials Integral Fractional Laplacian 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. Ministry of Education (MOE) Published version The work of the first author is partially supported by Shanghai Pujiang Program 21PJ1403500. The work of the second author is partially supported by the National Natural Science Foundation of China (No. 12101325). The work of the third author is partially supported by the National Natural Science Foundation of China (No. 11871455 and 11971016). The research of the fourth author is partially supported by Singapore MOE AcRF Tier 2 Grant: MOE2018- T2-1-059 and Tier 1 Grant: RG15/21. The research of the fifth author is supported in part by the National Natural Science Foundation of China (No. 11871092 and 12101035). 2022-06-30T07:04:23Z 2022-06-30T07:04:23Z 2021 Journal Article Sheng, C., Ma, S., Li, H., Wang, L. & Jia, L. (2021). Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators. ESAIM: Mathematical Modelling and Numerical Analysis, 55(5), 2141-2168. https://dx.doi.org/10.1051/m2an/2021049 0764-583X https://hdl.handle.net/10356/159739 10.1051/m2an/2021049 2-s2.0-85117264195 5 55 2141 2168 en MOE2018-T2-1-059 RG15/21 ESAIM: Mathematical Modelling and Numerical Analysis © 2021 The authors. Published by EDP Sciences, SMAI. 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. application/pdf
institution Nanyang Technological University
building NTU Library
continent Asia
country Singapore
Singapore
content_provider NTU Library
collection DR-NTU
language English
topic Science::Mathematics
Generalised Hermite Polynomials
Integral Fractional Laplacian
spellingShingle Science::Mathematics
Generalised Hermite Polynomials
Integral Fractional Laplacian
Sheng, Changtao
Ma, Suna
Li, Huiyuan
Wang, Li-Lian
Jia, Lueling
Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
description 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.
author2 School of Physical and Mathematical Sciences
author_facet School of Physical and Mathematical Sciences
Sheng, Changtao
Ma, Suna
Li, Huiyuan
Wang, Li-Lian
Jia, Lueling
format Article
author Sheng, Changtao
Ma, Suna
Li, Huiyuan
Wang, Li-Lian
Jia, Lueling
author_sort Sheng, Changtao
title Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
title_short Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
title_full Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
title_fullStr Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
title_full_unstemmed Nontensorial generalised hermite spectral methods for PDEs with fractional Laplacian and Schrödinger operators
title_sort nontensorial generalised hermite spectral methods for pdes with fractional laplacian and schrödinger operators
publishDate 2022
url https://hdl.handle.net/10356/159739
_version_ 1759854869408120832