Abstract
We propose a simple polynomial basis set that is easily extendable to any desired higher-order accuracy. This method is based on the Constrained Interpolation Profile (CIP) method and the profile is chosen so that the subgrid scale solution approaches the real solution by the constraints from the spatial derivative of the original equation. By adopting the higher-order derivatives of the master equations as constraints to generate a self-consistent subgrid profile, this solution quickly converges. 3rd and 5th order polynomials are tested on the one-dimensional Schrödinger equation and are proved to give solutions a few orders of magnitude higher in accuracy than conventional methods for lower-lying eigenstates.
Original language | English |
---|---|
Pages (from-to) | 1-7 |
Number of pages | 7 |
Journal | Journal of Computational Physics |
Volume | 196 |
Issue number | 1 |
DOIs | |
State | Published - 1 May 2004 |
Externally published | Yes |
Keywords
- Basis sets
- CIP method
- CIP-BS method
- Dirichlet1 boundary conditions
- Generalized eigenvalue equation
- Neuman boundary conditions
- Time-dependent Schrödinger equation