search

njoy / NJOY2016

Nuclear data processing with legacy NJOY
https://www.njoy21.io/NJOY2016
Other
96 stars 86 forks source link

Coulomb phase shift when APE .ne. APT #182

Closed toshihikokawano closed 3 years ago

toshihikokawano commented 3 years ago

Line 3224 in samm.f90, phase-shift should be calculated by the effect…ive radius APE

whaeck commented 3 years ago

@toshihikokawano Interesting, I've stared at this part of the code a few times already in the last couple of years and never noticed.

Since the phase shift appears to be calculated using the wrong rho due to this mixup, this might explain the small differences with the new resonance reconstructionimplementation if the Ca40 resonance parameters do define two different radii (which I suspect they do). If that is true, that would also mean that your coulomb wave function implementation that I've been using could have been correct all along.

By the way: I like your avatar picture ;-)

whaeck commented 3 years ago

Just verified: Ca40 does indeed have 2 radii defined: 0.4993153 for calculating P and S 0.3828792 for the phase shift