Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCoordinate-Space Hartree-Fock-Bogoliubov Solvers for Superfluid Fermi Systems in Large Boxes
119
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The self-consistent Hartree-Fock-Bogoliubov problem in large boxes can be solved accurately in the coordinate space with the recently developed solvers HFB-AX (2D) and MADNESS-HFB (3D). This is essential for the description of superfluid Fermi systems with complicated topologies and significant spatial extend, such as fissioning nuclei, weakly-bound nuclei, nuclear matter in the neutron star rust, and ultracold Fermi atoms in elongated traps. The HFB-AX solver based on B-spline techniques uses a hybrid MPI and OpenMP programming model for parallel computation for distributed parallel computation, within a node multi-threaded LAPACK and BLAS libraries are used to further enable parallel calculations of large eigensystems. The MADNESS-HFB solver uses a novel multi-resolution analysis based adaptive pseudo-spectral techniques to enable fully parallel 3D calculations of very large systems. In this work we present benchmark results for HFB-AX and MADNESS-HFB on ultracold trapped fermions.
rate research
Read More
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.
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.
Condensed Fermi systems with an odd number of particles can be described by means of polarizing external fields having a time-odd character. We illustrate how this works for Fermi gases and atomic nuclei treated by density functional theory or Hartree-Fock-Bogoliubov (HFB) theory. We discuss the method based on introducing two chemical potentials for different superfluid components, whereby one may change the particle-number parity of the underlying quasiparticle vacuum. Formally, this method is a variant of non-collective cranking, and the procedure is equivalent to the so-called blocking. We present and exemplify relations between the two-chemical-potential method and the cranking approximation for Fermi gases and nuclei.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
