Self-consistent calculation of nuclear photoabsorption cross section: Finite amplitude method with Skyrme functionals in the three-dimensional real space

The finite amplitude method (FAM), which we have recently proposed (T. Nakatsukasa, T. Inakura, and K. Yabana, Phys. Rev. C 76, 024318 (2007)), simplifies significantly the fully self-consistent RPA calculation. Employing the FAM, we are conducting systematic, fully self-consistent response calculations for a wide mass region. This paper is intended to present a computational scheme to be used in the systematic investigation and to show the performance of the FAM for a realistic Skyrme energy functional. We implemented the method in the mixed representation in which the forward and backward RPA amplitudes are represented by indices of single-particle orbitals for occupied states and the spatial grid points for unoccupied states. We solve the linear response equation for a given frequency. The equation is a linear algebraic problem with a sparse non-hermitian matrix, which is solved with an iterative method. We show results of the dipole response for selected spherical and deformed nuclei. The peak energies of the giant dipole resonance agree well with measurements for heavy nuclei, while they are systematically underestimated for light nuclei. We also discuss the width of the giant dipole resonance in the fully self-consistent RPA calculation.Comment: 11 pages, 10 figure

Linear Responses in Time-dependent Hartree-Fock-Bogoliubov Method with Gogny Interaction

A numerical method to integrate the time-dependent Hartree-Fock Bogoliubov (TDHFB) equations with Gogny interaction is proposed. The feasibility of the TDHFB code is illustrated by the conservation of the energy, particle numbers, and center-of-mass in the small amplitude vibrations of oxygen 20. The TDHFB code is applied to the isoscalar quadrupole and/or isovector dipole vibrations in the linear (small amplitude) region in oxygen isotopes (masses A = 18,20,22 and 24), titanium isotopes (A = 44,50,52 and 54), neon isotope (A = 26), and magnesium isotopes (A = 24 and 34). The isoscalar quadrupole and isovector dipole strength functions are calculated from the expectation values of the isoscalar quadrupole and isovector dipole moments.Comment: 10 pages, 13 figure

Beta-gamma spectroscopy of the neutron-rich 150Ba

International audience; Excited states in the neutron-rich nucleus 150Ba have been observed via β–γ spectroscopy at the Radioactive Isotope Beam Factory, RIKEN Nishina Center. The 150Ba ions were produced by the in-flight fission of a 238U beam with an energy of 345 MeV/nucleon. The E(2+) energy of 150Ba was identified at 100 keV, which is the lowest known in the neutron-rich Ba isotopes. A γ -ray peak was also observed at 597 keV. A mean-field calculation with a fully 3D real space was performed and a static octupole deformation was obtained for the Ba isotopes. Kπ = 0− and 1− excited states with significant octupole collectivity were newly predicted at around or lower than 1 MeV on the ground state of 150Ba by a random-phase approximation calculation. The 597 keV γ ray can be interpreted as a negative-parity state, showing that 150Ba may possess octupole collectivity

Evaluations of uncertainties in simulations of propagation of ultrahigh-energy cosmic-ray nuclei derived from microscopic nuclear models

Photodisintegration is a main energy loss process of ultrahigh-energy cosmic-ray (UHECR) nuclei in intergalactic space at the highest energies. Therefore, it is crucial to understand photodisintegration's systematic uncertainty to simulate the propagation of UHECR nuclei. In this work, we calculated the cross sections using the random phase approximation (RPA) of density functional theory (DFT), a microscopic nuclear model. We calculated the $E1$ strength of 29 nuclei using three different density functionals. We obtained the cross sections of photonuclear reactions, including photodisintegration, with the $E1$ strength. Then, we implemented the cross sections in a cosmic ray propagation code CRPropa. We found that the difference between the RPA calculations and TALYS in CRPropa in the energy spectrum can be more than the statistical uncertainty of UHECR energy spectrum assuming some astrophysical parameters. We also found that the difference of some astrophysical parameters obtained by a combined fit of UHECR energy spectrum and composition data can be more than the uncertainty of the data between the RPA calculations and TALYS assuming a phenomenological model of UHECR sources.Comment: 36 pages, 18 figure