Abstract We present a code for solving the single-particle, time-independent Schrodinger equation in two dimensions. Our program utilizes the imaginary time propagation (ITP) algorithm, and it includes the most recent developments in the ITP method: the arbitrary order operator factorization and the exact inclusion of a (possibly very strong) magnetic field. Our program is able to solve thousands of eigenstates of a two-dimensional quantum system in reasonable time with commonly available hardware. The main motivation behind our work is to allow the study of highly excited states and energy spectra of two-dimensional quantum dots and billiard systems with a single versatile code, e.g., in quantum chaos research. In our implementation we emphasize a modern and easily extensible design, simple and user-friendly interfaces, and an open-source development philosophy. Program summary Program title: itp2d Catalogue identifier: AENR_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AENR_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU General Public License version 3 No. of lines in distributed program, including test data, etc.: 11310 No. of bytes in distributed program, including test data, etc.: 97720 Distribution format: tar.gz Programming language: C++ and Python. Computer: Tested on x86 and x86-64 architectures. Operating system: Tested under Linux with the g++ compiler. Any POSIX-compliant OS with a C++ compiler and the required external routines should suffice. Has the code been vectorised or parallelized?: Yes, with OpenMP. RAM: 1 MB or more, depending on system size. Classification: 7.3. External routines: FFTW3 ( http://www.fftw.org ), CBLAS ( http://netlib.org/blas ), LAPACK ( http://www.netlib.org/lapack ), HDF5 ( http://www.hdfgroup.org/HDF5 ), OpenMP ( http://openmp.org ), TCLAP ( http://tclap.sourceforge.net ), Python ( http://python.org ), Google Test ( http://code.google.com/p/googletest/ ) Nature of problem: Numerical calculation of the lowest energy solutions (up to a few thousand, depending on available memory), of a single-particle, time-independent Schrdinger equation in two dimensions with or without a homogeneous magnetic field. Solution method: Imaginary time propagation (also known as the diffusion algorithm), with arbitrary even order factorization of the imaginary time evolution operator Additional comments: Please see the README file distributed with the program for more information. The source code of our program is also available at https: //bitbucket.org/luukko/itp2d. Running time: Seconds to hours
[1]
Andrew Hassell,et al.
Boundary Quasi-Orthogonality and Sharp Inclusion Bounds for Large Dirichlet Eigenvalues
,
2010,
SIAM J. Numer. Anal..
[2]
Andreas Stathopoulos,et al.
PRIMME: preconditioned iterative multimethod eigensolver—methods and software description
,
2010,
TOMS.
[3]
E. Krotscheck,et al.
An arbitrary order diffusion algorithm for solving Schrödinger equations
,
2009,
Comput. Phys. Commun..
[4]
P. Löwdin.
On the Non‐Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals
,
1950
.
[5]
Á. Rubio,et al.
octopus: a first-principles tool for excited electron-ion dynamics.
,
2003
.
[6]
Š. Janeček,et al.
Billiards in magnetic fields: a molecular dynamics approach.
,
2009,
Physical review. E, Statistical, nonlinear, and soft matter physics.
[7]
H. Stöckmann,et al.
Quantum Chaos: An Introduction
,
1999
.
[8]
Š. Janeček,et al.
Gauge-invariant real-space method for density functional calculations in an external magnetic field
,
2008
.
[9]
Vicente Hernández,et al.
SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems
,
2005,
TOMS.
[10]
Square billiard with a magnetic flux
,
2000,
Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.
[11]
Aguiar.
Eigenvalues and eigenfunctions of billiards in a constant magnetic field.
,
1996,
Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.
[12]
J. Toivanen,et al.
Solution of time-independent Schrödinger equation by the imaginary time propagation method
,
2007,
J. Comput. Phys..
[13]
B. C. Carlson,et al.
Orthogonalization Procedures and the Localization of Wannier Functions
,
1957
.
[14]
E. Krotscheck,et al.
A fast and simple program for solving local Schrödinger equations in two and three dimensions
,
2008,
Comput. Phys. Commun..
[15]
B. M. Fulk.
MATH
,
1992
.
[16]
Nils Berglund,et al.
Integrability and ergodicity of classical billiards in a magnetic field
,
1995,
chao-dyn/9501009.
[17]
Steven G. Johnson,et al.
The Design and Implementation of FFTW3
,
2005,
Proceedings of the IEEE.
[18]
S. Chin.
Multi-product splitting and Runge-Kutta-Nyström integrators
,
2008,
0809.0914.
[19]
M. Aichinger,et al.
Fourth-order algorithms for solving local Schrödinger equations in a strong magnetic field
,
2005,
Comput. Phys. Commun..