Numerical approach of some three-body problems

Abstract On taking the two-proton decay in nuclear physics as a sample case we formulate a numerical method for three-body problems which are mathematically described by systems of coupled 2d Schrodinger equations in polar coordinates r , ϕ . With some minimal adaptations the method becomes applicable on other cases. The specific feature of the procedure consists in a shuttle propagation along r : two quantities are propagated, the log-derivative matrix and the solution itself, with a backwards propagation of the former, followed by a forwards propagation of the latter. The results show remarkable stability. Numerical illustrations from a simple test model are reported and we also explain how the data obtained in this way can be exploited to obtain useful physical information.

[1]  L.Gr. Ixaru LILIX—A package for the solution of the coupled channel Schrödinger equation☆ , 2002 .

[2]  H. De Meyer,et al.  CP methods for the Schro¨dinger equation revisited , 1998 .

[3]  B. A. Brown,et al.  Two-proton radioactivity of 45Fe. , 2002 .

[4]  B. M. Fulk MATH , 1992 .

[5]  Liviu Gr. Ixaru New numerical method for the eigenvalue problem of the 2D Schrödinger equation , 2010, Comput. Phys. Commun..

[6]  M. Zhukov,et al.  Two-proton radioactivity and three-body decay. II. Exploratory studies of lifetimes and correlations , 2003 .

[7]  Guido Vanden Berghe,et al.  CP methods of higher order for Sturm-Liouville and Schrödinger equations , 2004, Comput. Phys. Commun..

[8]  P. Salamon,et al.  New simple form for a phenomenological nuclear potential , 2008 .

[9]  V. Goldansky On neutron-deficient isotopes of light nuclei and the phenomena of proton and two-proton radioactivity , 1960 .

[10]  R. C. Johnson,et al.  Two proton radioactivity and three-body decay: General problems and theoretical approach , 2001 .

[11]  Guido Vanden Berghe,et al.  MATSLISE: A MATLAB package for the numerical solution of Sturm-Liouville and Schrödinger equations , 2005, TOMS.

[12]  Liviu Gr Ixaru,et al.  Numerical methods for differential equations and applications , 1984 .

[13]  A. Stolz,et al.  First observation of two-proton radioactivity in Ni48 , 2011 .

[14]  H. De Meyer,et al.  SLCPM12 - A program for solving regular Sturm-Liouville problems. , 1999 .

[15]  Marnix Van Daele,et al.  Automatic computation of quantum-mechanical bound states and wavefunctions , 2013, Comput. Phys. Commun..

[16]  D. Delion,et al.  Simple approach to two-proton emission , 2013 .