(both non-relativistic and relativistic variants of the method)
Last update: 29/01/2007
The J-matrix method is one of the algebraic methods in quantum scattering theory. It allows to find an approximate solution of the scattering problem on the radial potential V=V(r), vanishing faster than the Coulomb one. The method is based on fact, that the radial kinetic energy operator is tridiagonal in some suitable bases - forms so called Jacobi matrix. See documentation for full theoretical description of the method.
The main goal of one of the scientific projects realized in our Department, is further developing the J-matrix method as well as the computer code JMATRIX, implementing the method and allowing for calculations of both relativistic and non-relativistic phase shifts. The project started in 1999. After developing the relativistic variant of the method by Pawel Horodecki, it became necessary to check validity of the method in real calculations. The author of the first version of the JMATRIX code (prepared as a part of his MSc thesis), as well as the next versions, was Pawel Syty. Test calculations, performed on some model potentials, has proved correctness of the method, however, it became evident, that programming the method is inherently connected with many numerical problems. After these first results were published (see publications), the project was suspended, and the authors were engaged to work on the other projects.
In the current year the project has been revived. Computer code has been significantly improved, time of calculations greatly reduced, and some new useful functions introduced. Moreover, there is a work in progress on improving the method itself, to achieve a faster convergence.
The main idea of this page is to release effects of our work on developing the J-matrix method, in particular the code JMATRIX. In addition, Wolfram Research's Mathematica notebooks are available, allowing for visual tracing how the method works, and some other related stuff.
Author gives permission to use and modify provided codes and notebooks freely, under the following conditions.
Fortran 95 sources
You need a Fortran 95 compiler to build executables from provided sources. See links below for free Fortran 95 compilers. The package also includes basic documentation and sample input file.
Previous versions (obsolete):
We provide some precompiled binary executables for the most popular operating systems, they are built from the current sources. However, to achieve the best performance, consider compiling program on the destination machine, from the provided sources. Moreover, in case of using provided binaries, it is not possible to specify user-potential, since it requires recompilation of the code. You also need to prepare the input file to start the calculations, sample file you can find here.
Mathematica 5 notebooks
Provided notebooks may be useful in understanding the method, however, in real calculations Fortran 95 code is much quicker, therefore highly recommended. If you don't have access to the Mathematica package, you may use free MathReader, but only for reading.
The list is rather short, because the project was suspended in the years 2000-2006.
P. Syty, Programming the J-matrix method, MSc thesis, Gdansk University of Technology (1999)
Example calculations and results (graphs in the Adobe PDF format)
1. Scattering from square-well potential: V(r)=V0 for a<=r<=b, 0 for r<a or r>b. Non-relativistic and relativistic calculations for V0=-1, a=0.8, b=1. Projectile parameters: l=1, kappa=1, E=3 a.u.
2. Scattering from truncated Coulomb potential: V(r)=-Z/r for r<=r0, 0 for r>r0. Non-relativistic calculations for Z=30, r0=1. Projectile parameters: l=1, kappa=1, E=3 a.u.
3. Scattering from truncated Coulomb potential: V(r)=-Z/r for r<=r0, 0 for r>r0. Relativistic calculations for Z=30, r0=1. Projectile parameters: various l and kappa, E=0.4 a.u.
4. Scattering from Yukawa potential: V(r)=-g2Exp(-mr)/r for r<=r0, 0 for r>r0. Non-relativistic calculations for g=1, m=1, r0=1. Projectile parameters: l=1, kappa=1, E=3 a.u.
Some utilities related to project (i.e. applications to create input files) will be published soon.
The G95 project
- free Fortran 95 compiler, available for multiple cpu architectures and operating systems
All comments and suggestions concerning the method, code or this WWW-page are welcome. Please send them to code developer's e-mail address: email@example.com