Neutron densities from a global analysis of medium-energy proton-nucleus elastic scattering

A new method for extracting neutron densities from intermediate-energy elastic proton-nucleus scattering observables uses a global Dirac phenomenological approach based on the relativistic impulse approximation. Datasets for ^40Ca, ^48Ca, and ^208Pb in the energy range from 500 MeV to 1040 MeV are considered. The global fits are successful in reproducing the data and in predicting datasets not included in the analysis. Using this global approach, energy-independent neutron densities are obtained. The vector point proton density distribution ₕ^p is determined from the empirical charge density after unfolding the proton form factor. The other densities, ₕ^n, ₒ^p, ₒ^n, are parametrized. This work provides energy-independent values for the rms neutron radius R₍ and the neutron skin thickness S₍, in contrast to the energy-dependent values obtained by previous studies. In addition, the results presented in this paper show that the expected rms neutron radius and the skin thickness for ^40Ca are accurately reproduced. The values of R₍ and S₍ obtained from the global fits that we consider to be the most reliable are given as follows: for ^40Ca, 3. 314>R₍>3. 310fm and -0. 063>S₍ >-0. 067fm; for ^48Ca, 3. 459>R₍>3. 413fm and 0. 102>S₍>0. 056fm; and for ^208Pb, 5. 550>R₍>5. 522fm and 0. 111>S₍>0. 083fm. These values are in reasonable agreement with nonrelativistic Skyrme-Hartree-Fock models and with relativistic Hartree-Bogoliubov models with density-dependent meson-nucleon couplings. The results from the global fits for ^48Ca and ^208Pb are generally not in agreement with the usual relativistic mean-field models.