Neck Rupture and Scission Neutrons in Nuclear Fission.

Just before a nucleus undergoes fission, a neck is formed between the emerging fission fragments. It is widely accepted that this neck undergoes a rather violent rupture, despite the absence of unambiguous experimental evidence. The main difficulty in addressing the neck rupture and saddle-to-scission stages of fission is that both are highly nonequilibrium processes. Here, we present the first fully microscopic characterization of the scission mechanism, along with the spectrum and the spatial distribution of scission neutrons (SNs), and some upper limit estimates for the emission of charged particles. The spectrum of SNs has a distinct angular distribution, with neutrons emitted in roughly equal numbers in the equatorial plane and along the fission axis. They carry an average energy around 3±0.5 MeV for the fission of ^{236}U, ^{240}Pu, and ^{252}Cf, and a maximum of 16-18 MeV. We estimate a conservative lower bound of 9%-14% of the total emitted neutrons are produced at scission.

Fragment Intrinsic Spins and Fragments' Relative Orbital Angular Momentum in Nuclear Fission.

We present the first fully unrestricted microscopic calculations of the primary fission fragment intrinsic spins and of the fission fragments' relative orbital angular momentum for ^{236}U^{*}, ^{240}Pu^{*}, and ^{252}Cf using the time-dependent density functional theory framework. Within this microscopic approach, free of restrictions and unchecked assumptions and which incorporates the relevant physical observables for describing fission, we evaluate the triple distribution of the fission fragment intrinsic spins and of their fission fragments' relative orbital angular momentum and show that their dynamics is dominated by their bending collective modes in contradistinction to the predictions of the existing phenomenological models and some interpretations of experimental data.

Fission Fragment Intrinsic Spins and Their Correlations.

The intrinsic spins and their correlations are the least understood characteristics of fission fragments from both theoretical and experimental points of view. In many nuclear reactions, the emerging fragments are typically excited and acquire an intrinsic excitation energy and an intrinsic spin depending on the type of the reactions and interaction mechanism. Both the intrinsic excitation energies and the fragments' intrinsic spins and parities are controlled by the interaction mechanism and conservations laws, which lead to their correlations and determines the character of their deexcitation mechanism. We outline here a framework for the theoretical extraction of the intrinsic spin distributions of the fragments and their correlations within the fully microscopic real-time density-functional theory formalism and illustrate it on the example of induced fission of ^{236}U and ^{240}Pu, using two nuclear energy density functionals. These fission fragment intrinsic spin distributions display new qualitative features previously not discussed in literature. Within this fully microscopic framework, we extract for the first time the intrinsic spin distributions of fission fragments of ^{236}U and ^{240}Pu as well as the correlations of their intrinsic spins, which have been debated in literature for more than six decades with no definite conclusions so far.

Emergence of a Pseudogap in the BCS-BEC Crossover.

Strongly correlated Fermi systems with pairing interactions become superfluid below a critical temperature T_{c}. The extent to which such pairing correlations alter the behavior of the liquid at temperatures T>T_{c} is a subtle issue that remains an area of debate, in particular regarding the appearance of the so-called pseudogap in the BCS-BEC crossover of unpolarized spin-1/2 nonrelativistic matter. To shed light on this, we extract several quantities of crucial importance at and around the unitary limit, namely, the odd-even staggering of the total energy, the spin susceptibility, the pairing correlation function, the condensate fraction, and the critical temperature T_{c}, using a nonperturbative, constrained-ensemble quantum Monte Carlo algorithm.

Dynamics of Fragmented Condensates and Macroscopic Entanglement.

The relative phase of the order parameters in the collision of two condensates can influence the outcome of their collision in the case of weak coupling. With increasing interaction strength, however, the initially independent phases of the two order parameters in the colliding partners quickly become phase locked, as the strong coupling favors an overall phase rigidity of the entire condensate, and upon their separation the emerging superfluid fragments become entangled.

Induced Fission of (240)Pu within a Real-Time Microscopic Framework.

We describe the fissioning dynamics of ^{240}Pu from a configuration in the proximity of the outer fission barrier to full scission and the formation of the fragments within an implementation of density functional theory extended to superfluid systems and real-time dynamics. The fission fragments emerge with properties similar to those determined experimentally, while the fission dynamics appears to be quite complex, with many excited shape and pairing modes. The evolution is found to be much slower than previously expected, and the ultimate role of the collective inertia is found to be negligible in this fully nonadiabatic treatment of nuclear dynamics, where all collective degrees of freedom (CDOF) are included (unlike adiabatic treatments with a small number of CDOF).

Vortex Pinning and Dynamics in the Neutron Star Crust.

The nature of the interaction between superfluid vortices and the neutron star crust, conjectured by Anderson and Itoh in 1975 to be at the heart vortex creep and the cause of glitches, has been a long-standing question in astrophysics. Using a qualitatively new approach, we follow the dynamics as superfluid vortices move in response to the presence of "nuclei" (nuclear defects in the crust). The resulting motion is perpendicular to the force, similar to the motion of a spinning top when pushed. We show that nuclei repel vortices in the neutron star crust, and characterize the force per unit length of the vortex line as a function of the vortex element to the nucleus separation.

Quantized superfluid vortex rings in the unitary Fermi gas.

In a recent article, Yefsah et al. [Nature (London) 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 magnitude 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 density functional theories.

Strength of the vortex-pinning interaction from real-time dynamics.

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 associated 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.

Cooper pairing above the critical temperature in a unitary Fermi gas.

We present an ab initio determination of the spin response of the unitary Fermi gas. Based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism, we determine the temperature dependence of the spin susceptibility and the spin conductivity. We show that both quantities exhibit suppression above the critical temperature of the superfluid-to-normal phase transition due to Cooper pairing. The spin diffusion transport coefficient does not display a minimum in the vicinity of the critical temperature and drops to very low values D(s)≈0.8h/m in the superfluid phase. All these spin observables show a smooth and monotonic behavior with temperature when crossing the critical temperature T(c), until the Fermi liquid regime is attained at the temperature T(*), above which the pseudogap regime disappears.

Quantum shock waves and domain walls in the real-time dynamics of a superfluid unitary Fermi gas.

We show that in the collision of two superfluid fermionic atomic clouds one observes the formation of quantum shock waves as discontinuities in the number density and collective flow velocity. Domain walls, which are topological excitations of the superfluid order parameter, are also generated and exhibit abrupt phase changes by π and slower motion than the shock waves. The domain walls are distinct from the gray soliton train or number density ripples formed in the wake of the shock waves and observed in the collisions of superfluid bosonic atomic clouds. Domain walls with opposite phase jumps appear to collide elastically.

Onset of a pseudogap regime in ultracold Fermi gases.

We show, using an ab initio approach based on Quantum Monte Carlo technique, that the pseudogap regime emerges in ultracold Fermi gases close to the unitary point. We locate the onset of this regime at a value of the interaction strength corresponding to (k(F)a)(-1)≈-0.05 (a-scattering length). We determine the evolution of the gap as a function of temperature and interaction strength in the Fermi gas around the unitary limit and show that our results exhibit a remarkable agreement with the recent wave-vector resolved radio frequency spectroscopy data. Our results indicate that a finite temperature structure of the Fermi gas around unitarity is more complicated and involves the presence of the phase with preformed Cooper pairs, which, however, do not contribute to the long range order.

Finite-temperature pairing gap of a unitary fermi gas by quantum monte carlo calculations.

We calculate the one-body temperature Green's (Matsubara) function of the unitary Fermi gas via quantum Monte Carlo, and extract the spectral weight function A(p,omega) using the methods of maximum entropy and singular value decomposition. From A(p,omega) we determine the quasiparticle spectrum, which can be accurately parametrized by three functions of temperature: an effective mass m{*}, a mean-field potential U, and a gap Delta. Below the critical temperature T{c}=0.15 epsilon{F} the results for m{*}, U, and Delta can be accurately reproduced using an independent quasiparticle model. We find evidence of a pseudogap in the fermionic excitation spectrum for temperatures up to T{*} approximately 0.20 epsilon{F}> T{c}.

Real-time dynamics of quantized vortices in a unitary Fermi superfluid.

We introduce a comprehensive theoretical framework for the fermionic superfluid dynamics, grounded on a local extension of the time-dependent density functional theory. With this approach, we describe the generation and the real-time evolution and interaction of quantized vortices, the large-amplitude collective modes, as well as the loss of superfluidity at high flow velocities. We demonstrate the formation of vortex rings and provide a microscopic description of the crossing and reconnection of quantized vortex lines in a fermion superfluid, which provide the mechanism for the emergence of quantum turbulence at very low temperatures. We observe that superfluidity often survives when these systems are stirred with velocities far exceeding the speed of sound.

Large amplitude dynamics of the pairing correlations in a unitary fermi gas.

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-Ginzburg-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 nonspherical momentum distribution, and could be excited by a quick time variation of the scattering length.

Unitary Fermi supersolid: the Larkin-Ovchinnikov phase.

We present strong theoretical evidence that a Larkin-Ovchinnikov (LOFF/FFLO) pairing phase is favored 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.

Thermodynamics of a trapped unitary Fermi gas.

We present the first model-independent comparison of recent measurements of the entropy and of the critical temperature of a unitary Fermi gas, performed by Luo et al., with the most complete results currently available from finite temperature Monte Carlo calculations. The measurement of the critical temperature in a cold fermionic atomic cloud is consistent with a value T(c) = 0.23(2)epsilon(F) in the bulk, as predicted by the present authors in their Monte Carlo calculations.

Spin 1/2 fermions in the unitary regime: a superfluid of a new type.

We study, in a fully nonperturbative calculation, a dilute system of spin 1/2 interacting fermions, characterized by an infinite scattering length at finite temperatures. Various thermodynamic properties and the condensate fraction are calculated and we also determine the critical temperature for the superfluid-normal phase transition in this regime. The thermodynamic behavior appears as a rather surprising and unexpected mélange of fermionic and bosonic features. The thermal response of a spin 1/2 fermion at the BCS-BEC crossover should be classified as that of a new type of superfluid.

Induced P-wave superfluidity in asymmetric fermi gases.

We show that two new intraspecies P-wave superfluid phases appear in two-component asymmetric Fermi systems with short-range S-wave interactions. In the BEC limit, phonons of the molecular BEC induce P-wave superfluidity in the excess fermions. In the BCS limit, density fluctuations induce P-wave superfluidity in both the majority and the minority species. These phases may be realized in experiments with spin-polarized Fermi gases.

Specific heat of a fermionic atomic cloud in the unitary regime.

In the unitary regime, when the scattering amplitude greatly exceeds in magnitude the average interparticle separation, and below the critical temperature thermal properties of an atomic fermionic cloud are governed by the collective modes, specifically the Bogoliubov-Anderson sound modes. The specific heat of an atomic cloud in an elongated trap, in particular, has a rather complex temperature dependence, which changes from an exponential behavior at very low temperatures (T << h omega(parallel)), to proportional T for h omega(parallel) << T << h omega(perpendicular) and then continuously to proportional T4 at temperatures just below the critical temperature, when the surface modes play a dominant role. Only the low (h omega(parallel) << T << h omega(perpendicular)) and high (h omega(perpendicular) << T << T(c)) temperature power laws are well defined. For the intermediate temperatures one can introduce at most a gradually increasing with temperature exponent.