Hex  1.0
Hydrogen-electron collision solver
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Hex Documentation
Author
Jakub Benda, MFF UK
Date
1. 2. 2014

Hex-dwba

Hex-dwba is a high-energy computational module of the Hex package. It implements the distorted wave Born approximation of the first order and is expected to give trustworthy results for ehergies above 1 keV.

Theory

The distorted wave Born approximation of the first order (DWBA-1) is similar to the well-known plane wave Born approximation. However, to speed up the convergence of the Born series (of which we are using just a single term in the first-order computation), part of the potential is excluded from the perurbative description and solved precisely. The spliting is the following

\[ V = W + U \, \]

(1)

where \( V \) is the full potential felt by the projectile, \( U \) is the so called "distorting potential" that will be solved precisely as mentioned above and explained below. This splitting has two consequences:

Implementation

The program iterates over partial waves of the asymptotic final state. For every iterations is will loop over all initial states that can contribute to the requested scattering process and still result in the current asymptotic final state.

The radial functions \( \chi_{k\ell}(r) \) are computed using numerical solution of the equation (2) by the adaptive Cash-Karp Runge-Kutta method provided by GSL (see below).

Requirements

You will need a modern C++ compiler to compile the program as for example

DWBA uses some free external libraries

The first is used for numerical integration and special and utitlity functions. The second is used for large-precision rational numbers that are necessary when working with analytic formulas. Both libraries are open-source and maintained. They can be downloaded free of charge and compiled on any Posix-compatible system. The compilation in Windows system can be done e.g. in the Code::Blocks IDE (that is not at all straightforward, though).