
(both nonrelativistic and relativistic variants of the method) Last update: 29/01/2007
Overview The Jmatrix 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 Jmatrix method as well as the computer code JMATRIX, implementing the method and allowing for calculations of both relativistic and nonrelativistic 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 Jmatrix 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. Program features
Fortran 95 sources Current version: 2.01 / 28.01.2007 [gzipped tar archive] [rar archive] [changelog] 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): Binaries 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 userpotential, 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.
Publications The list is rather short, because the project was suspended in the years 20002006.
P. Syty, Programming the Jmatrix method, MSc thesis, Gdansk University of Technology (1999) Example calculations and results (graphs in the Adobe PDF format) 1. Scattering from squarewell potential: V(r)=V_{0} for a<=r<=b, 0 for r<a or r>b. Nonrelativistic and relativistic calculations for V_{0}=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<=r_{0}, 0 for r>r_{0}. Nonrelativistic calculations for Z=30, r_{0}=1. Projectile parameters: l=1, kappa=1, E=3 a.u.
3. Scattering from truncated Coulomb potential: V(r)=Z/r for r<=r_{0}, 0 for r>r_{0}. Relativistic calculations for Z=30, r_{0}=1. Projectile parameters: various l and kappa, E=0.4 a.u. In preparation 4. Scattering from Yukawa potential: V(r)=g^{2}Exp(mr)/r for r<=r_{0}, 0 for r>r_{0}. Nonrelativistic calculations for g=1, m=1, r_{0}=1. Projectile parameters: l=1, kappa=1, E=3 a.u. In preparation Utilities 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 Contact All comments and suggestions concerning the method, code or this WWWpage are welcome. Please send them to code developer's email address: sylas@mif.pg.gda.pl 

