A structure-preserving relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Liu, Huini, Yi, Nianyu, Yin, Peimeng
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909798838042624
author Liu, Huini
Yi, Nianyu
Yin, Peimeng
author_facet Liu, Huini
Yi, Nianyu
Yin, Peimeng
contents In this paper, we propose a mass- and modified energy-conservative relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation. Utilizing only a single auxiliary variable, we simultaneously reformulate the distinct nonlinear terms present in both the Schrödinger equation and the Poisson equation into their equivalent expressions, constructing a system equivalent to the original Schrödinger-Poisson equation. Our proposed scheme, derived from this equivalent system, is implemented linearly, avoiding the need for iterative techniques to solve the nonlinear equation. Additionally, it is executed sequentially, eliminating the need to solve a coupled large linear system. We in turn rigorously derive the optimal error estimates for the proposed scheme, demonstrating second order accuracy in time and $(k+1)$th order accuracy in space when employing polynomials of degree up to $k$. Numerical experiments validate the accuracy and effectiveness of our method and emphasize its conservation properties over long-time simulations.
format Preprint
id arxiv_https___arxiv_org_abs_2405_12848
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A structure-preserving relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation
Liu, Huini
Yi, Nianyu
Yin, Peimeng
Numerical Analysis
35Q55, 65M15, 65M60
In this paper, we propose a mass- and modified energy-conservative relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation. Utilizing only a single auxiliary variable, we simultaneously reformulate the distinct nonlinear terms present in both the Schrödinger equation and the Poisson equation into their equivalent expressions, constructing a system equivalent to the original Schrödinger-Poisson equation. Our proposed scheme, derived from this equivalent system, is implemented linearly, avoiding the need for iterative techniques to solve the nonlinear equation. Additionally, it is executed sequentially, eliminating the need to solve a coupled large linear system. We in turn rigorously derive the optimal error estimates for the proposed scheme, demonstrating second order accuracy in time and $(k+1)$th order accuracy in space when employing polynomials of degree up to $k$. Numerical experiments validate the accuracy and effectiveness of our method and emphasize its conservation properties over long-time simulations.
title A structure-preserving relaxation Crank-Nicolson finite element method for the Schrödinger-Poisson equation
topic Numerical Analysis
35Q55, 65M15, 65M60
url https://arxiv.org/abs/2405.12848