A Simple Quantum Integro-Differential Solver (SQuIDS)

作者：Delgado Carlos A Argueelles; Salvado Jordi^*; Weaver Christopher N

来源：Computer Physics Communications, 2015, 196: 569-591.

DOI：10.1016/j.cpc.2015.06.022

• Quantum mechanics
• Ordinary differential equations
• SU groups

Simple Quantum Integro-Differential Solver (SQUIDS) is a C++ code designed to solve semi-analytically the evolution of a set of density matrices and scalar functions. This is done efficiently by expressing all operators in an SU (N) basis. SQuIDS provides a base class from which users can derive new classes to include new non-trivial terms from the right hand sides of density matrix equations. The code was designed in the context of solving neutrino oscillation problems, but can be applied to any problem that involves solving the quantum evolution of a collection of particles with Hilbert space of dimension up to six.
Program summary
Program title: SQuIDS
Catalogue identifier: AEXG_v1_0
Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEXG_v1_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: GNU Lesser General Public License, version 3
No. of lines in distributed program, including test data, etc.: 18198
No. of bytes in distributed program, including test data, etc.: 137607
Distribution format: tar.gz
Programming language: C++11.
Computer: 32- and 64-bit x86.
Operating system: Linux, Mac OS X, FreeBSD.
RAM: Proportional to the number of nodes, the dimension of the Hilbert space, the number of scalar functions, and the number of density matrices used in the problem. Classification: 11.1.
External routines: GNU Scientific Library (http://www.gnu.org/software/gsl/).
Nature of problem: Solve the evolution of open quantum systems of Hilbert space dimension N with self interactions and interaction with classical fields.
Solution method: The SU (N) algebra is implemented as a C++ object and is embedded into the GSL ordinary differential equation solver.
Restrictions: The code is only implemented up to Hilbert spaces of dimension six, but a Mathematica notebook is provided in order to generate higher dimensional solutions. Furthermore, only ordinary differential equation solution methods that require only the first derivative can be used.
Running time: Proportional to the number of nodes, the dimension of the Hilbert space, the number of scalar functions, the number of density matrices, and the numerical precision used in the problem.

• 出版日期2015-11