On the computation of ground state and dynamics of Schrödinger-Poisson-Slater system

  • Authors:

  • Yong Zhang;Xuanchun Dong

  • Affiliations:

  • Department of Mathematical Sciences, Tsinghua University, Beijing 100084, PR China;Department of Mathematics and Center for Computational Science and Engineering, National University of Singapore, 119076 Singapore, Singapore

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.46

Visualization

Abstract

In this paper, we deal with the computation of ground state and dynamics of the Schrodinger-Poisson-Slater (SPS) system. To this end, backward Euler and time-splitting pseudospectral methods are proposed for the nonlinear Schrodinger equation with the nonlocal Hartree potential approximated by solving a Poisson equation. The approximation approaches for the Hartree potential include fast convolution algorithms, which are accelerated by using FFT in 1D and fast multipole method (FMM) in 2D and 3D, and sine/Fourier pseudospectral methods. The inconsistency in 0-mode in Fourier pseudospectral approach is pointed out, which results in a significant loss of high-order of accuracy as expected for spectral methods. Numerical comparisons show that in 1D the fast convolution and sine pseudospectral approaches are compatible. While, in 3D the fast convolution approach based on FMM is second-order accurate and the Fourier pseudospectral approach is better than it from both efficiency and accuracy point of view. Among all these approaches, the sine pseudospectral one is the best candidate in the numerics of the SPS system. Finally, we apply the backward Euler and time-splitting sine pseudospectral methods to study the ground state and dynamics of 3D SPS system in different setups.