Improved convergence of scattering calculations in the oscillator representation
This work provides a more efficient computational method for scattering calculations in nuclear physics, though it is an incremental improvement over existing techniques.
The authors present a hybrid representation for solving scattering states in the Schrödinger equation, combining harmonic oscillator eigenstates with a finite difference grid, and demonstrate significant convergence improvement over the JM-ECS method.
The Schrödinger equation for two and tree-body problems is solved for scattering states in a hybrid representation where solutions are expanded in the eigenstates of the harmonic oscillator in the interaction region and on a finite difference grid in the near-- and far--field. The two representations are coupled through a high--order asymptotic formula that takes into account the function values and the third derivative in the classical turning points. For various examples the convergence is analyzed for various physics problems that use an expansion in a large number of oscillator states. The results show significant improvement over the JM-ECS method [Bidasyuk et al, Phys. Rev. C 82, 064603 (2010)].