No Arabic abstract
In this paper we present a new formalism to implement the nuclear particle-vibration coupling (PVC) model. The key issue is the proper treatment of the continuum, that is allowed by the coordinate space representation. Our formalism, based on the use of zero-range interactions like the Skyrme forces, is microscopic and fully self-consistent. We apply it to the case of neutron single-particle states in $^{40}$Ca, $^{208}$Pb and $^{24}$O. The first two cases are meant to illustrate the comparison with the usual (i.e., discrete) PVC model. However, we stress that the present approach allows to calculate properly the effect of PVC on resonant states. We compare our results with those from experiments in which the particle transfer in the continuum region has been attempted. The latter case, namely $^{24}$O, is chosen as an example of a weakly-bound system. Such a nucleus, being double-magic and not displaying collective low-lying vibrational excitations, is characterized by quite pure neutron single-particle states around the Fermi surface.
Nuclear $beta$-decay in magic nuclei is investigated, taking into account the coupling between particle and collective vibrations,on top of self-consistent random phase approximation calculations based on Skyrme density functionals. The low-lying Gamow-Teller strength is shifted downwards and at times becomes fragmented; as a consequence, the $beta$-decay half-lives are reduced due to the increase of the phase space available for the decay. In some cases, this leads to a very good agreement between theoretical and experimental lifetimes: this happens, in particular, in the case of the Skyrme force SkM*, that can also reproduce the line shape of the high energy Gamow-Teller resonance as it was previously shown.
The microscopic description of neutron scattering by $^{16}$O below 30 MeV is carried out by means of the continuum particle-vibration coupling (cPVC) method with the Skyrme nucleon-nucleon ($NN$) effective interaction. In the cPVC method, a proper boundary condition on a nucleon in continuum states is imposed, which enables one to evaluate the transition matrix in a straightforward manner. Experimental data of the total and total-elastic cross sections are reproduced quite well by the cPVC method. An important feature of the result is the fragmentation of the single-particle resonance into many peaks as well as the shift of its centroid energy. Thus, some part of the fine structure of the experimental cross sections at lower energies is well described by the cPVC framework. The cPVC method based on a real $NN$ effective interaction is found to successfully explain about 85% of the reaction cross section, through explicit channel-coupling effects.
It has been known that the time-dependent Hartree-Fock (TDHF) method, or the time-dependent density functional theory (TDDFT), fails to describe many-body quantum tunneling. We overcome this problem by superposing a few time-dependent Slater determinants with the time-dependent generator coordinate method (TDGCM). We apply this method to scattering of two $alpha$ particles in one dimension, and demonstrate that the TDGCM method yields a finite tunneling probability even at energies below the Coulomb barrier, at which the tunneling probability is exactly zero in the TDHF. This is the first case in which a many-particle tunneling is simulated with a microscopic real-time approach.
In order to study structure of proto-neutron stars and those in subsequent cooling stages, it is of great interest to calculate inhomogeneous hot and cold nuclear matter in a variety of phases. The finite-temperature Hartree-Fock-Bogoliubov (FT-HFB) theory is a primary choice for this purpose, however, its numerical calculation for superfluid (superconducting) many-fermion systems in three dimensions requires enormous computational costs. To study a variety of phases in the crust of hot and cold neutron stars, we propose an efficient method to perform the FT-HFB calculation with the three-dimensional (3D) coordinate-space representation. Recently, an efficient method based on the contour integral of Greens function with the shifted conjugate-orthogonal conjugate-gradient method has been proposed [Phys. Rev. C 95, 044302 (2017)]. We extend the method to the finite temperature, using the shifted conjugate-orthogonal conjugate-residual method. We benchmark the 3D coordinate-space solver of the FT-HFB calculation for hot isolated nuclei and fcc phase in the inner crust of neutron stars at finite temperature. The computational performance of the present method is demonstrated. Different critical temperatures of the quadrupole and the octupole deformations are confirmed for $^{146}$Ba. The robustness of the shape coexistence feature in $^{184}$Hg is examined. For the neutron-star crust, the deformed neutron-rich Se nuclei embedded in the sea of superfluid low-density neutrons appear in the fcc phase at the nucleon density of 0.045 fm$^{-3}$ and the temperature of $k_B T=200$ keV. The efficiency of the developed solver is demonstrated for nuclei and inhomogeneous nuclear matter at finite temperature. It may provide a standard tool for nuclear physics, especially for the structure of the hot and cold neutron-star matters.
The Energy Density Functional theory is one of the most used methods developed in nuclear structure. It is based on the assumption that the energy of the ground state is a functional only of the density profile. The method is extremely successful within the effective force approach, noticeably the Skyrme or Gogny forces, in reproducing the nuclear binding energies and other bulk properties along the whole mass table. Although the Density Functional is in this case represented formally as the Hartree-Fock mean field of an effective force, the corresponding single-particle states in general do not reproduce the phenomenology particularly well. To overcome this difficulty, a strategy has been developed where the effective force is adjusted to reproduce directly the single particle energies, trying to keep the ground state energy sufficiently well reproduced. An alternative route, that has been developed along several years, for solving this problem is to introduce the mean field fluctuations, as represented by the collective vibrations of the nuclear system, and their influence on the single particle dynamics and structure. This is the basis of the particle-vibration coupling model. In this paper we present a formal theory of the particle-vibration coupling model based on the Green s function method. The theory extends to realistic effective forces the macroscopic particle-vibration coupling models and the (microscopic) Nuclear Field Theory. It is formalized within the functional derivative approach to many-body theory. An expansion in diagrams is devised for the single particle self-energy and the phonon propagator. Critical aspects of the particle-vibration coupling model are analysed in general. Applications at the lowest order of the expansion are presented and discussed.