Source author record

Aurel Bulgac

Aurel Bulgac appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

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.

preprint2022arXiv

Sensitivity to the initial conditions of the Time-Dependent Density Functional Theory

Time-Dependent Density Functional Theory is mathematically formulated through non-linear coupled time-dependent 3-dimensional partial differential equations and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is however linear and mathematically the exact equations of the Density Functional Theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a Density Functional Theory framework. Whether for realistic implementations of the Time-Dependent Density Functional Theory the question of absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density functional theory is a unique tool allowing us the study of non-equilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom ${\cal O}(10^{10})$ and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

preprint2020arXiv

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 non-perturbative, constrained-ensemble quantum Monte Carlo algorithm.

preprint2020arXiv

Framework for Polarized Superfluid Fermion Systems

I discuss the advantages and disadvantages of several procedures, some known and some new, for constructing stationary states within the mean field approximation for a system with pairing correlations and unequal numbers spin-up and spin-down fermions, using the two chemical potentials framework. One procedure in particular appears to have significant physics advantages over previously suggested in the literature computational frameworks. Moreover, this framework is applicable to study strongly polarized superfluid fermion systems with arbitrarily large polarizations and with arbitrary total particle numbers. These methods are equally applicable to normal systems.

preprint2020arXiv

Nuclear Fission Dynamics: Past, Present, Needs, and Future

Recent developments in theoretical modeling and in computational power have allowed us to make significant progress on a goal not achieved yet in nuclear theory: a fully microscopic theory of nuclear fission. The complete microscopic description remains a computationally demanding task, but the information that can be provided by current calculations can be extremely useful to guide and constrain phenomenological approaches. First, a truly microscopic framework that can describe the real-time dynamics of the fissioning system can justify or rule out assumptions and approximations incompatible with an accurate quantum treatment or with our understanding of the inter nucleon interactions. Second, the microscopic approach can be used to obtain trends such as: the excitation energy sharing mechanism between fission fragments (FFs) with increasing excitation energy of the fissioning system, the angular momentum content of the FFs, or even to compute observables that cannot be otherwise calculated in phenomenological approaches or even measured, as in the case of astronomical environments. Merely the characterization of the trends would be of great importance for various application. We present here arguments that a truly microscopic approach to fission does not support the assumption of adiabaticity of the large amplitude collective motion in fission, particularly starting from the outer saddle down to the scission configuration.

preprint2016arXiv

Induced Fission of $^{240}$Pu within a Real-Time Microscopic Framework

We describe the fissioning dynamics of Pu240 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).

preprint2016arXiv

Towards Quantum Turbulence in Cold Atomic Fermionic Superfluids

Fermionic superfluids provide a new realization of quantum turbulence, accessible to both experiment and theory, yet relevant to phenomena from both cold atoms to nuclear astrophysics. In particular, the strongly interacting Fermi gas realized in cold-atom experiments is closely related to dilute neutron matter in neutron star crusts. Unlike the liquid superfluids 4He (bosons) and 3He (fermions) where quantum turbulence has been studied in the laboratory, superfluid Fermi gases stand apart for a number of reasons. They admit a reliable theoretical description based on a DFT called the TDSLDA that describes both static and dynamic phenomena. Cold atom experiments demonstrate exquisite control over particle number, spin polarization, density, temperature, and interaction strength. Topological defects such as domain walls and quantized vortices, which lie at the heart of quantum turbulence, can be created and manipulated with time-dependent external potentials, and agree with the time-dependent theoretical techniques. While similar experimental and theoretical control exists for weakly interacting Bose gases, the unitary Fermi gas is strongly interacting. The resulting vortex line density is extremely high, and quantum turbulence may thus be realized in small systems where classical turbulence is suppressed. Fermi gases also permit the study of exotic superfluid phenomena such as a 3D LOFF supersolid, and a finite temperature pseudo-gap in the regime of classical turbulence. The dynamics associated with these phenomena has only started to be explored. Finally, superfluid mixtures have recently been realized, providing experimental access to phenomena like Andreev-Bashkin entrainment. Superfluid Fermi gases thus provide a rich forum for addressing phenomena related to quantum turbulence with applications ranging from terrestrial superfluidity to astrophysical dynamics in neutron stars.

preprint2016arXiv

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 as a function of the vortex-nucleus separation.

preprint2015arXiv

Nuclear Energy Density Functionals: What do we really know?

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-Weizsäcker 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 Weizsäcker's 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.

preprint2015arXiv

Perfect-fluid behavior of a dilute Fermi gas near unitary

We present an ab initio calculation of the shear viscosity as a function of interaction strength in a two-component unpolarized Fermi gas near the unitary limit, within a finite temperature quantum Monte Carlo (QMC) framework and using the Kubo linear-response formalism. The shear viscosity decreases as we tune the interaction strength 1/ak_F from the Bardeen-Cooper-Schrieffer side of the Feshbach resonance towards Bose-Einstein condensation limit and it acquires the smallest value for 1/ak_F approx 0.4, with a minimum value of (eta/s)_min approx 0.2 hbar/k_B, which is about twice as small as the value reported for experiments in quark-gluon plasma (eta/s)_QGP lesssim 0.4 hbar/k_B. The Fermi gas near unitarity thus emerges as the most "perfect fluid" observed so far in nature. The clouds of dilute Fermi gas near unitarity exhibit the unusual attribute that, for the sizes realized so far in the laboratory or larger (less than 10^9 atoms), can sustain quantum turbulence below the critical temperature, but at the same time the classical turbulence is suppressed in the normal phase.

preprint2014arXiv

Life Cycle of Superfluid Vortices and Quantum Turbulence in the Unitary Fermi Gas

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.

preprint2014arXiv

Quantized Superfluid Vortex Rings in the Unitary Fermi Gas

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

preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

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

preprint2013arXiv

Cooper pairing above the critical temperature in a unitary Fermi gas

We present an ab initio determination of spin responses 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 presence of the Cooper pairs. The spin diffusion transport coefficient does not display the existence of a minimum in the vicinity of the critical temperature and it drops to very low values D_s approx 0.8hbar/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*, where the pseudogap regime disappears.

preprint2013arXiv

Quantum Friction: Cooling Quantum Systems with Unitary Time Evolution

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 friction 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 Schrödinger equation significantly simpler than the solution of its stationary counterpart.

preprint2013arXiv

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.

preprint2013arXiv

The temperature evolution of the shear viscosity in a unitary Fermi gas

We present an ab initio determination of the shear viscosity for the unitary Fermi gas based on finite temperature quantum Monte Carlo (QMC) calculations and the Kubo linear-response formalism. The results are confronted with the bound for the shear viscosity originating from hydrodynamic fluctuations. Assuming smoothness of the frequency dependent shear viscosity eta(omega), we show that the bound is violated in the low temperature regime and the violation occurs simultaneously with the onset of the Cooper paring in the system. In order to preserve the hydrodynamic bound in QMC eta(omega) has to possess a sharp structure located in the vicinity of zero frequency which is not resolved by an analytic continuation procedure.

preprint2013arXiv

Time-Dependent Density Functional Theory and the Real-Time Dynamics of Fermi Superfluids

I describe the Time-Dependent Superfluid Local Density Approximation, which is an adiabatic extension of the Density Functional Theory to superfluid Fermi systems and their real-time dynamics. This new theoretical framework has been applied to describe a number of phenomena in cold atomic gases and nuclear collective motion: excitation of the Higgs modes in strongly interacting Fermi superfluids, generation of quantized vortices, crossing and reconnection of vortices, excitation of the superflow at velocities above the critical velocity, excitation of quantum shock waves and domain walls in the collisions of superfluid atomic clouds, excitation of collective states in nuclei.

preprint2013arXiv

Time-Dependent Superfluid Local Density Approximation

The time-dependent superfluid local density approximation (TDSLDA) is an extension of the Hohenberg-Kohn density functional theory (DFT) to time-dependent phenomena in superfluid fermionic systems. Unlike linear response theory, which is only valid for weak external fields, the (TDSLDA) approach allows one to study non-linear excitations in fermionic superfluids, including large amplitude collective modes, and the response to strong external probes. Even in the case of weak external fields, the (TDSLDA) approach is technically easier to implement. We will illustrate the implementation of the (TDSLDA) for the unitary Fermi gas, where dimensional arguments and Galilean invariance simplify the form of the functional, and ab initio input from (QMC) simulations fix the coefficients to quite high precision.

preprint2013arXiv

Use of the Discrete Variable Representation Basis in Nuclear Physics

The discrete variable representation (DVR) basis is nearly optimal for numerically representing wave functions in nuclear physics: Suitable problems enjoy exponential convergence, yet the Hamiltonian remains sparse. We show that one can often use smaller basis sets than with the traditional harmonic oscillator basis, and still benefit from the simple analytic properties of the DVR bases which requires no overlap integrals, simply permit using various Jacobi coordinates, and admit straightforward analyses of the ultraviolet and infrared convergence properties.

preprint2012arXiv

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.

preprint2012arXiv

The Unitary Fermi Gas: From Monte Carlo to Density Functionals

In this chapter, we describe three related studies of the universal physics of two-component unitary Fermi gases with resonant short-ranged interactions. First we discuss an ab initio auxiliary field quantum Monte Carlo technique for calculating thermodynamic properties of the unitary gas from first principles. We then describe in detail a Density Functional Theory (DFT) fit to these thermodynamic properties: the Superfluid Local Density Approximation (SLDA) and its Asymmetric (ASLDA) generalization. We present several applications, including vortex structure, trapped systems, and a supersolid Larkin-Ovchinnikov (FFLO/LOFF) state. Finally, we discuss the time-dependent extension to the density functional (TDDFT) which can describe quantum dynamics in these systems, including non-adiabatic evolution, superfluid to normal transitions and other modes not accessible in traditional frameworks such as a Landau-Ginzburg, Gross-Pitaevskii, or quantum hydrodynamics.

preprint2011arXiv

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}\approx -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.

preprint2010arXiv

Generation and Dynamics of Quantized Vortices in a Unitary Fermi Superfluid

Superfluidity and superconductivity are remarkable manifestations of quantum coherence at a macroscopic scale. The dynamics of superfluids has dominated the study of these systems for decades now, but a comprehensive theoretical framework is still lacking. We introduce a local extension of the time-dependent density functional theory to describe the dynamics of fermionic superfluids. Within this approach one can correctly represent vortex quantization, generation, and dynamics, the transition from a superfluid to a normal phase and a number of other large amplitude collective modes which are beyond the scope of two-fluid hydrodynamics, Ginzburg-Landau and/or Gross-Pitaevskii approaches. We illustrate the power of this approach by studying the generation of quantized vortices, vortex rings, vortex reconnection, and transition from a superfluid to a normal state in real time for a unitary Fermi gas. We predict the emergence of a new qualitative phenomenon in superfluid dynamics of gases, the existence of stable superfluidity when the systems are stirred with velocities significantly exceeding the nominal Landau critical velocity in these systems.

preprint2010arXiv

The Long Journey from Ab Initio Calculations to Density Functional Theory for Nuclear Large Amplitude Collective Motion

At present there are two vastly different ab initio approaches to the description of the the many-body dynamics: the Density Functional Theory (DFT) and the functional integral (path integral) approaches. On one hand, if implemented exactly, the DFT approach can allow in principle the exact evaluation of arbitrary one-body observable. However, when applied to Large Amplitude Collective Motion (LACM) this approach needs to be extended in order to accommodate the phenomenon of surface-hoping, when adiabaticity is strongly violated and the description of a system using a single (generalized) Slater determinant is not valid anymore. The functional integral approach on the other hand does not appear to have such restrictions, but its implementation does not appear to be straightforward endeavor. However, within a functional integral approach one seems to be able to evaluate in principle any kind of observables, such as the fragment mass and energy distributions in nuclear fission. These two radically approaches can likely be brought brought together by formulating a stochastic time-dependent DFT approach to many-body dynamics.

preprint2006arXiv

Spin 1/2 Fermions in the Unitary Regime: A Superfluid of a New Type

We have studied, in a fully non-perturbative 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 were calculated and we have also determined the critical temperature for the superfluid-normal phase transition in this regime. The thermodynamic behavior appears as a rather surprising and unexpected melange 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.

preprint2004arXiv

Collective oscillations of a trapped Fermi gas near the unitary limit

We calculate the oscillation frequencies of trapped Fermi condensate with particular emphasis on the equation of state of the interacting Fermi system. We confirm Stringari's finding that the frequencies are independent of the interaction in the unitary limit, and we extend the theory away from that limit, where the interaction does affect the frequencies of the compressional modes only.

preprint2002arXiv

Dilute Quantum Droplets

In the limit when the two-body scattering length $a$ is negative and much larger than the effective two-body interaction radius the contribution to the ground state energy due to the three-body correlations is given by the Efimov effect. For particular values of the diluteness parameter $ρ|a|^3$ the three-body contribution can become the dominant term of the energy density functional. Under these conditions both Bose-Einstein and Fermi-Dirac systems could become self-bound and either boson droplets or fermion ``designer nuclei'' of various sizes and densities could be manufactured.

preprint1999arXiv

Proton Single Particle Energy Shifts due to Coulomb Correlations

A theoretically consistent approach to the calculation of the Coulomb correlation corrections to the single--particle energies is presented. New contributions to the single--particle energies previously overlooked in the literature are now identified and taken into account. We show that the interplay between the Coulomb interaction and the strong interaction, which is enhanced in the nuclear surface, leads to an upward shift of the proton single--particle levels. This shift affects the position of the calculated proton drip line, a shift towards decreasing $Z$. We describe briefly a similar mechanism which is at work for neutron levels. The same mechanism is responsible for significant corrections to the mass difference of the mirror nuclei (Nolen--Schiffer anomaly) and to the effective proton mass.