We present the simplest nuclear energy density functional (NEDF) to date, determined by only 4 significant phenomenological parameters, yet capable of fitting measured nuclear masses with better accuracy than the Bethe-Weizsacker mass formula, while
also describing density structures (charge radii, neutron skins etc.) and time-dependent phenomena (induced fission, giant resonances, low energy nuclear collisions, etc.). The 4 significant parameters are necessary to describe bulk nuclear properties (binding energies and charge radii); an additional 2 to 3 parameters have little influence on the bulk nuclear properties, but allow independent control of the density dependence of the symmetry energy and isovector excitations, in particular the Thomas-Reiche-Kuhn sum rule. This Hohenberg-Kohn-style of density functional theory successfully realizes Weizsackers ideas and provides a computationally tractable model for a variety of static nuclear properties and dynamics, from finite nuclei to neutron stars, where it will also provide a new insight into the physics of the r-process, nucleosynthesis, and neutron star crust structure. This new NEDF clearly separates the bulk geometric properties - volume, surface, symmetry, and Coulomb energies which amount to 8MeV per nucleon or up to 2000MeV per nucleus for heavy nuclei - from finer details related to shell effects, pairing, isospin breaking, etc. which contribute at most a few MeV for the entire nucleus. Thus it provides a systematic framework for organizing various contributions to the NEDF. Measured and calculated physical observables - symmetry and saturation properties, the neutron matter equation of state, and the frequency of giant dipole resonances - lead directly to new terms not considered in current NEDF parameterizations.
The unitary Fermi gas (UFG) offers an unique opportunity to study quantum turbulence both experimentally and theoretically in a strongly interacting fermionic superfluid. It yields to accurate and controlled experiments, and admits the only dynamical
microscopic description via time-dependent density functional theory (DFT) - apart from dilute bosonic gases - of the crossing and reconnection of superfluid vortex lines conjectured by Feynman in 1955 to be at the origin of quantum turbulence in superfluids at zero temperature. We demonstrate how various vortex configurations can be generated by using well established experimental techniques: laser stirring and phase imprinting. New imagining techniques demonstrated by the MIT group [Ku et al. arXiv:1402.7052] should be able to directly visualize these crossings and reconnections in greater detail than performed so far in liquid helium. We demonstrate the critical role played by the geometry of the trap in the formation and dynamics of a vortex in the UFG and how laser stirring and phase imprint can be used to create vortex tangles with clear signatures of the onset of quantum turbulence.
We present a simple Hamiltonian description of the dynamics of a quantized vortex ring in a trapped superfluid, compare this description with dynamical simulations, and characterize the dependence of the dynamics of the shape of the trap.
In a recent article, Yefsah et al. [Nature 499, 426 (2013)] report the observation of an unusual excitation in an elongated harmonically trapped unitary Fermi gas. After phase imprinting a domain wall, they observe oscillations almost an order of mag
nitude slower than predicted by any theory of domain walls which they interpret as a heavy soliton of inertial mass some 200 times larger than the free fermion mass or 50 times larger than expected for a domain wall. We present compelling evidence that this soliton is instead a quantized vortex ring by showing that the main aspects of the experiment can be naturally explained within the framework of time-dependent superfluid DFT.
We introduce a type of quantum dissipation -- local quantum friction -- by adding to the Hamiltonian a local potential that breaks time-reversal invariance so as to cool the system. Unlike the Kossakowski-Lindblad master equation, local quantum frict
ion directly effects unitary evolution of the wavefunctions rather than the density matrix: it may thus be used to cool fermionic many-body systems with thousands of wavefunctions that must remain orthogonal. In addition to providing an efficient way to simulate quantum dissipation and non-equilibrium dynamics, local quantum friction coupled with adiabatic state preparation significantly speeds up many-body simulations, making the solution of the time-dependent Schrodinger equation significantly simpler than the solution of its stationary counterpart.
The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists ca
n lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay.
We present an efficient and general method to compute vortex-pinning interactions - which arise in neutron stars, superconductors, and trapped cold atoms - at arbitrary separations using real-time dynamics. This method overcomes uncertainties associa
ted with matter redistribution by the vortex position and the related choice of ensemble that plague the typical approach of comparing energy differences between stationary pinned and unpinned configurations: uncertainties that prevent agreement in the literature on the sign and magnitude of the vortex-nucleus interaction in the crust of neutron stars. We demonstrate and validate the method with Gross-Pitaevskii-like equations for the unitary Fermi gas, and demonstrate how the technique of adiabatic state preparation with time-dependent simulation can be used to calculate vortex-pinning interactions in fermionic systems such as the vortex-nucleus interaction in the crust of neutron stars.
A unitary Fermi gas has a surprisingly rich spectrum of large amplitude modes of the pairing field alone, which defies a description within a formalism involving only a reduced set of degrees of freedom, such as quantum hydrodynamics or a Landau-Ginz
burg-like description. These modes are very slow, with oscillation frequencies well below the pairing gap, which makes their damping through quasiparticle excitations quite ineffective. In atomic traps these modes couple naturally with the density oscillations, and the corresponding oscillations of the atomic cloud are an example of a new type of collective mode in superfluid Fermi systems. They have lower frequencies than the compressional collective hydrodynamic oscillations, have a non-spherical momentum distribution, and could be excited by a quick time variation of the scattering length.
We present strong theoretical evidence that a Larkin-Ovchinnikov (LOFF/FFLO) pairing phase is favoured over the homogeneous superfluid and normal phases in three-dimensional unitary Fermi systems. Using a Density Functional Theory (DFT) based on the
latest quantum Monte-Carlo calculations and experimental results, we show that this phase is competitive over a large region of the phase diagram. The oscillations in the number densities and pairing field have a substantial amplitude, and a period some 3 to 10 times the average interparticle separation. Within the DFT, the transition to a normal polarized Fermi liquid at large polarizations is smooth, while the transition to a fully-paired superfluid is abrupt.
The Quantum Monte Carlo method for spin 1/2 fermions at finite temperature is formulated for dilute systems with an s-wave interaction. The motivation and the formalism are discussed along with descriptions of the algorithm and various numerical issu
es. We report on results for the energy, entropy and chemical potential as a function of temperature. We give upper bounds on the critical temperature T_c for the onset of superfluidity, obtained by studying the finite size scaling of the condensate fraction. All of these quantities were computed for couplings around the unitary regime in the range -0.5 le (k_F a)^{-1} le 0.2, where a is the s-wave scattering length and k_F is the Fermi momentum of a non-interacting gas at the same density. In all cases our data is consistent with normal Fermi gas behavior above a characteristic temperature T_0 > T_c, which depends on the coupling and is obtained by studying the deviation of the caloric curve from that of a free Fermi gas. For T_c < T < T_0 we find deviations from normal Fermi gas behavior that can be attributed to pairing effects. Low temperature results for the energy and the pairing gap are shown and compared with Green Function Monte Carlo results by other groups.
