A novel, structure-preserving, second-order-in-time relaxation scheme for Schrödinger-Poisson systems

Agissilaos Athanassoulis (Lead / Corresponding author), Theodoros Katsaounis, Irene Kyza (Lead / Corresponding author), Stephen Metcalfe

Research output: Contribution to journalArticlepeer-review

1 Citation (Scopus)
20 Downloads (Pure)

Abstract

We introduce a new structure preserving, second order in time relaxation-type scheme for approximating solutions of the Schrödinger-Poisson system. More specifically, we use the Crank-Nicolson scheme as a time stepping mechanism, whilst the nonlinearity is handled by means of a relaxation approach in the spirit of [10], [11], [34] for the nonlinear Schrödinger equation. For the spatial discretisation we use the standard conforming finite element scheme. The resulting scheme is explicit with respect to the nonlinearity, i.e. it requires the solution of a linear system for each time-step, and satisfies discrete versions of the system's mass conservation and energy balance laws for constant meshes. The scheme is seen to be second order in time. We conclude by presenting some numerical experiments, including an example from cosmology and an example with variable time-steps which demonstrate the effectiveness and robustness of the new scheme.
Original languageEnglish
Article number112307
Number of pages18
JournalJournal of Computational Physics
Volume490
Early online date22 Jun 2023
DOIs
Publication statusPublished - 1 Oct 2023

Keywords

  • Schrödinger-Poisson system
  • energy preserving scheme
  • relaxation scheme in time
  • finite element method
  • Finite element method
  • Energy preserving scheme
  • Relaxation scheme in time

ASJC Scopus subject areas

  • Computational Mathematics
  • General Physics and Astronomy
  • Applied Mathematics
  • Numerical Analysis
  • Computer Science Applications
  • Modelling and Simulation
  • Physics and Astronomy (miscellaneous)

Fingerprint

Dive into the research topics of 'A novel, structure-preserving, second-order-in-time relaxation scheme for Schrödinger-Poisson systems'. Together they form a unique fingerprint.

Cite this