No Arabic abstract
Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density approximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vortices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.
We discuss properties of the method based on time dependent superfluid local density approximation (TDSLDA) within an application to induced fission of 240Pu and surrounding nuclei. Various issues related to accuracy of time evolution and the determination of the properties of fission fragments are discussed.
The dynamic response of asymmetric nuclear matter is studied by using a Time-Dependent Local Isospin Density (TDLIDA) approximation approach. Calculations are based on a local density energy functional derived by an Auxiliary Field Diffusion Monte Carlo (AFDMC) calculation of bulk nuclear matter. Three types of excited states emerge: collective states, a continuum of quasi-particle-quasi-hole excitations and unstable solutions. These states are analyzed and discussed for different values of the nuclear density $rho$ and isospin asymmetry $xi=(N-Z)/A$. An analytical expression of the compressibility as a function of $rho$ and $xi$ is derived which show explicitly an instability of the neutron matter around $rhosimeq 0.09 fm^{-3}$ when a small fraction of protons is added to the system.
Here we describe the form of the Asymmetric Superfluid Local Density Approximation (ASLDA), a Density Functional Theory (DFT) used to model the two-component unitary Fermi gas. We give the rational behind the functional, and describe explicitly how we determine the form of the DFT from the to the available numerical and experimental data.
Firstly, a systematic procedure is derived for obtaining three-dimensional bound-state equations from four-dimensional ones. Unlike ``quasi-potential approaches this procedure does not involve the use of delta-function constraints on the relative four-momentum. In the absence of negative-energy states, the kernels of the three-dimensional equations derived by this technique may be represented as sums of time-ordered perturbation theory diagrams. Consequently, such equations have two major advantages over quasi-potential equations: they may easily be written down in any Lorentz frame, and they include the meson-retardation effects present in the original four-dimensional equation. Secondly, a simple four-dimensional equation with the correct one-body limit is obtained by a reorganization of the generalized ladder Bethe-Salpeter kernel. Thirdly, our approach to deriving three-dimensional equations is applied to this four-dimensional equation, thus yielding a retarded interaction for use in the three-dimensional bound-state equation of Wallace and Mandelzweig. The resulting three-dimensional equation has the correct one-body limit and may be systematically improved upon. The quality of the three-dimensional equation, and our general technique for deriving such equations, is then tested by calculating bound-state properties in a scalar field theory using six different bound-state equations. It is found that equations obtained using the method espoused here approximate the wave functions obtained from their parent four-dimensional equations significantly better than the corresponding quasi-potential equations do.
We propose a computationally efficient approach to the nonadiabatic time-dependent density functional theory (TDDFT) which is based on a representation of the frequency-dependent exchange correlation kernel as a response of a set of damped oscillators. The requirements to computational resources needed to implement our approach do not differ from those of the standard real-time TDDFT in the adiabatic local density approximation (ALDA). Thus, our result offers an exciting opportunity to take into account temporal nonlocality and memory effects in calculations with TDDFT in quantum chemistry and solid state physics for unprecedentedly low costs.