No Arabic abstract
To study the exotic odd nuclear systems, the self-consistent continuum Skyrme-Hartree-Fock-Bogoliubov theory formulated with Greens function technique is extended to include blocking effects with the equal filling approximation. Detailed formula are presented.To perform the integrals of the Greens function properly, the contour paths $C_{rm b}^{-}$ and $C_{rm b}^{+}$ introduced for the blocking effects should include the blocked quasi-particle state but can not intrude into the continuum area. By comparing with the box-discretized calculations, the great advantages of the Greens function method in describing the extended density distributions, resonant states, and the couplings with the continuum in exotic nuclei are shown. Finally, taking the neutron-rich odd nucleus $^{159}$Sn as an example, the halo structure is investigated by blocking the quasi-particle state $1p_{1/2}$. It is found that it is mainly the weakly bound states near the Fermi surface that contribute a lot for the extended density distributions at large coordinate space.
Background: The Density-constraint Time-dependent Hartree-Fock method is currently the tool of choice to predict fusion cross-sections. However, it does not include pairing correlations, which have been found recently to play an important role. Purpose: To describe the fusion cross-section with a method that includes the superfluidity and to understand the impact of pairing on both the fusion barrier and cross-section. Method: The density-constraint method is tested first on the following reactions without pairing, $^{16}$O+$^{16}$O and $^{40}$Ca+$^{40}$Ca. A new method is developed, the Density-constraint Time-dependent Hartree-Fock-Bogoliubov method. Using the Gogny-TDHFB code, it is applied to the reactions $^{20}$O+$^{20}$O and $^{44}$Ca+$^{44}$Ca. Results: The Gogny approach for systems without pairing reproduces the experimental data well. The DC-TDHFB method is coherent with the TDHFB fusion threshold. The effect of the phase-lock mechanism is shown for those reactions. Conclusions: The DC-TDHFB method is a useful new tool to determine the fusion potential between superfluid systems and to deduce their fusion cross-sections.
We solve the Hartree-Fock-Bogoliubov (HFB) equations for a spherical mean field and a pairing potential with the inverse Hamiltonian method, which we have developed for the solution of the Dirac equation. This method is based on the variational principle for the inverse Hamiltonian, and is applicable to Hamiltonians that are bound neither from above nor below. We demonstrate that the method works well not only for the Dirac but also for the HFB equations.
Weakly-bound deformed nuclei have been studied by the Skyrme Hartree-Fock-Bogoliubov (HFB) approach in large coordinate-space boxes. In particular, the box-size dependence of the HFB calculations of weakly-bound deformed nuclei are investigated, including the particle density and pairing density distributions at nuclear surfaces, the near-threshold resonant and continuum quasiparticle spectra, and energetic properties. The box size may have larger influences in pairing properties than in other bulk properties. We demonstrate that large-box calculations of weakly-bound nuclei are important to precisely describe exotic phenomena such as deformed halos and peninsulas of stability beyond drip lines.
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 coordinate space formulation of the Hartree-Fock-Bogoliubov (HFB) method enables self-consistent treatment of mean-field and pairing in weakly bound systems whose properties are affected by the particle continuum space. Of particular interest are neutron-rich, deformed drip-line nuclei which can exhibit novel properties associated with neutron skin. To describe such systems theoretically, we developed an accurate 2D lattice Skyrme-HFB solver {hfbax} based on B-splines. Compared to previous implementations, we made a number of improvements aimed at boosting the solvers performance. These include: explicit imposition of axiality and space inversion, use of the modified Broydens method to solve self-consistent equations, and a partial parallelization of the code. {hfbax} has been benchmarked against other HFB solvers, both spherical and deformed, and the accuracy of the B-spline expansion was tested by employing the multiresolution wavelet method. Illustrative calculations are carried out for stable and weakly bound nuclei at spherical and very deformed shapes, including constrained fission pathways. In addition to providing new physics insights, {hfbax} can serve as a useful tool to assess the reliability and applicability of coordinate-space and configuration-space HFB solvers, both existing and in development.