Abstract:Background: Complex many-body systems, such as triaxial and reflection-asymmetric nuclei, weakly bound halo states, cluster configurations, nuclear fragments produced in heavy-ion fusion reactions, cold Fermi gases, and pasta phases in neutron star crust, are all characterized by large sizes and complex topologies in which many geometrical symmetries characteristic of ground-state configurations are broken. A tool of choice to study such complex forms of matter is an adaptive multi-resolution wavelet analysis. This method has generated much excitement since it provides a common framework linking many diversified methodologies across different fields, including signal processing, data compression, harmonic analysis and operator theory, fractals, and quantum field theory.Purpose: To describe complex superfluid many-fermion systems, we introduce an adaptive pseudospectral method for solving self-consistent equations of nuclear density functional theory in three dimensions, without symmetry restrictions.Methods: The numerical method is based on the multi-resolution and computational harmonic analysis techniques with a multi-wavelet basis. The application of state-of-the-art parallel programming techniques include sophisticated object-oriented templates which parse the high-level code into distributed parallel tasks with a multi-thread task queue scheduler for each multi-core node. The internode communications are asynchronous. The algorithm is variational and is capable of solving coupled complex-geometric systems of equations adaptively, with functional and boundary constraints, in a finite spatial domain of very large size, limited by existing parallel computer memory. For smooth functions, user-defined finite precision is guaranteed.Results: The new adaptive multi-resolution Hartree-Fock-Bogoliubov (HFB) solver MADNESS-HFB is bench-marked against a two-dimensional coordinate-space solver HFB-AX that is based on the B-spline technique and a three-dimensional solver HFODD that is based on the harmonic-oscillator basis expansion. Several examples are considered, including the self-consistent HFB problem for spin-polarized trapped cold fermions and the Skyrme-Hartree-Fock (+BCS) problem for triaxial deformed nuclei.Conclusions: The new MADNESS-HFB framework has many attractive features when applied to nuclear and atomic problems involving many-particle superfluid systems. Of particular interest are weakly bound nuclear configurations close to particle drip lines, strongly elongated and dinuclear configurations such as those present in fission and heavy-ion fusion, and exotic pasta phases that appear in neutron star crust.

Iterative solutions of the ATDHFB equations to determine the nuclear collective inertia

Efficient Method for Computing the Thouless-Valatin Inertia Parameters

Adaptive Multi-Resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure

Ab Initio Validation on the Connection Between Atomistic and Hydrodynamic Description to Unravel the Ion Dynamics of Warm Dense Matter

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method

Unified First-Principles Equations of State of Deuterium-Tritium Mixtures in the Global Inertial Confinement Fusion Region

A iterative finite method approach to Kohn-Sham equations of light atoms: Density Functional Theory tutorial for undergraduate students

Two-center harmonic oscillator basis for Skyrme-DFT calculations (I): formalism and Proof of Principle

A program for 3D nuclear static and time-dependent density-functional theory with full Skyrme energy density functional: HIT3D

Iterative subspace algorithms for finite-temperature solution of Dyson equation

Solution of universal nonrelativistic nuclear DFT equations in the Cartesian deformed harmonic-oscillator basis. (IX) HFODD (v3.06h): a new version of the program

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

A Well-Tempered Hybrid Method for Solving Challenging Time-Dependent Density Functional Theory (TDDFT) Systems

Effective and Efficient Algorithm for the Wigner Rotation Matrix at High Angular Momenta

Five-dimensional collective Hamiltonian with improved inertial functions

An Efficient Implementation of Analytical Nuclear Gradients for Linear-Response Time-Dependent Density Functional Theory in the Plane Wave Basis

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian

Deformed Coordinate-Space Hartree-Fock-Bogoliubov Approach to Weakly Bound Nuclei and Large Deformations

Multiconfigurational time-dependent density functional theory for atomic nuclei: Technical and numerical aspects

Time-dependent coupled cluster with orthogonal adaptive basis functions: General formalism and application to the vibrational problem