Showing 1–50 of 84 results for author: Magierski, P

Quantum vortices in fermionic superfluids: from ultracold atoms to neutron stars

Authors: Piotr Magierski, Andrea Barresi, Andrzej Makowski, Daniel Pęcak, Gabriel Wlazłowski

Abstract: Superfluid dilute neutron matter and ultracold gas, close to the unitary regime, exhibit several similarities. Therefore, to a certain extent, fermionic ultracold gases may serve as emulators of dilute neutron matter, which forms the inner crust of neutron stars and is not directly accessed experimentally. Quantum vortices are one of the most significant properties of neutron superfluid, essential… ▽ More Superfluid dilute neutron matter and ultracold gas, close to the unitary regime, exhibit several similarities. Therefore, to a certain extent, fermionic ultracold gases may serve as emulators of dilute neutron matter, which forms the inner crust of neutron stars and is not directly accessed experimentally. Quantum vortices are one of the most significant properties of neutron superfluid, essential for comprehending neutron stars' dynamics. The structure and dynamics of quantum vortices as a function of pairing correlations' strength are being investigated experimentally and theoretically in ultracold gases. Certain aspects of these studies are relevant to neutron stars. We provide an overview of the characteristics of quantum vortices in s-wave-type fermionic and electrically neutral superfluids. The main focus is on the dynamics of fermionic vortices and their intrinsic structure. △ Less

Submitted 16 September, 2024; v1 submitted 20 June, 2024; originally announced June 2024.

Journal ref: Eur. Phys. J. A 60, 186 (2024)

Stability of persistent currents in superfluid fermionic rings

Authors: Klejdja Xhani, Andrea Barresi, Marek Tylutki, Gabriel Wlazłowski, Piotr Magierski

Abstract: We investigate the stability of persistent currents in superfluid fermionic gases confined to a ring geometry. Our studies, conducted at zero temperature using time-dependent density functional theory, cover interaction regimes from strong (unitary Fermi gas) to weak (Bardeen-Cooper-Schrieffer regime) couplings. Stability is tested with respect to the presence of an external defect within the ring… ▽ More We investigate the stability of persistent currents in superfluid fermionic gases confined to a ring geometry. Our studies, conducted at zero temperature using time-dependent density functional theory, cover interaction regimes from strong (unitary Fermi gas) to weak (Bardeen-Cooper-Schrieffer regime) couplings. Stability is tested with respect to the presence of an external defect within the ring. The dissipation mechanism related to vortex generation is present in all interaction regimes. Interestingly, while the corresponding critical winding number is found to be independent of the regime, the flow energy dissipation and its origin strongly depend on it. Vortex emission is accompanied by Cooper pair breaking, which occurs even beyond the vortex core in the weakly interacting regime. The pair-breaking mechanism prevents the imprinting of a persistent current with a winding number above a threshold, which decreases as we approach the BCS regime. Our study reveals the existence of two types of critical winding numbers above which the currents cease to be persistent in Fermi superfluids: one related to the proliferation of quantum vortices and the other with the onset of the pair-breaking mechanism. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 8 pages, 4 figures

Time-dependent nuclear energy-density functional theory toolkit for neutron star crust: Dynamics of a nucleus in a neutron superfluid

Authors: Daniel Pȩcak, Agata Zdanowicz, Nicolas Chamel, Piotr Magierski, Gabriel Wlazłowski

Abstract: We present a new numerical tool designed to probe the dense layers of neutron star crusts. It is based on the Time-Dependent Hartree-Fock-Bogoliubov theory with generalized Skyrme nuclear energy density functionals, such as the Brussels-Montreal ones. We use it to study the time evolution of a nucleus accelerating through superfluid neutron medium in the inner crust of a neutron star. We extract a… ▽ More We present a new numerical tool designed to probe the dense layers of neutron star crusts. It is based on the Time-Dependent Hartree-Fock-Bogoliubov theory with generalized Skyrme nuclear energy density functionals, such as the Brussels-Montreal ones. We use it to study the time evolution of a nucleus accelerating through superfluid neutron medium in the inner crust of a neutron star. We extract an effective mass in the low velocity limit. We observe a threshold velocity and specify mechanisms of dissipation: phonon emission, Cooper pairs breaking, and vortex rings creation. The microscopic effects we study have impact on neutron star. Moreover, the mechanisms, we described, are general and apply also to other fermionic superfluid mixtures like liquid helium, or ultracold gases. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Generation and decay of Higgs mode in a strongly interacting Fermi gas

Authors: Andrea Barresi, Antoine Boulet, Gabriel Wlazłowski, Piotr Magierski

Abstract: We investigate the life cycle of the large amplitude Higgs mode in strongly interacting superfluid Fermi gas. Through numerical simulations with time-dependent density-functional theory and the technique of the interaction quench, we verify the previous theoretical predictions on the mode's frequency. Next, we demonstrate that the mode is dynamically unstable against external perturbation and qual… ▽ More We investigate the life cycle of the large amplitude Higgs mode in strongly interacting superfluid Fermi gas. Through numerical simulations with time-dependent density-functional theory and the technique of the interaction quench, we verify the previous theoretical predictions on the mode's frequency. Next, we demonstrate that the mode is dynamically unstable against external perturbation and qualitatively examine the emerging state after the mode decays. The post-decay state is characterized by spatial fluctuations of the order parameter and density at scales comparable to the superfluid coherence length scale. We identify similarities with FFLO states, which become more prominent at higher dimensionalities and nonzero spin imbalances. △ Less

Submitted 14 July, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: Ancillary files contain reproducibility packs of numerical calculations

Journal ref: Sci. Rep. 13, 11285 (2023)

Manifestation of pairing modes in nuclear collisions

Authors: A. Makowski, M. C. Barton, P. Magierski, K. Sekizawa, G. Wlazłowski

Abstract: We discuss the possible manifestation of pairing dynamics in nuclear collisions beyond the standard quasi-static treatment of pairing correlations. These involve solitonic excitations induced by pairing phase difference of colliding nuclei and pairing dynamic enhancement in the di-nuclear system formed by merging nuclei. We discuss the possible manifestation of pairing dynamics in nuclear collisions beyond the standard quasi-static treatment of pairing correlations. These involve solitonic excitations induced by pairing phase difference of colliding nuclei and pairing dynamic enhancement in the di-nuclear system formed by merging nuclei. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 2 figures, 56th Zakopane Conference On Nuclear Physics

Journal ref: Acta Phys. Pol. B Proc. Suppl. 16, 4-A38 (2023)

Disordered structures in ultracold spin-imbalanced Fermi gas

Authors: Buğra Tüzemen, Tomasz Zawiślak, Gabriel Wlazłowski, Piotr Magierski

Abstract: We investigate properties of spin-imbalanced ultracold Fermi gas in a large range of spin polarizations at low temperatures. We present results of microscopic calculations based on mean-field and density functional theory approaches, with no symmetry constraints. At low polarization values we predict the structure of the system as consisting of several spin-polarized droplets. As the polarization… ▽ More We investigate properties of spin-imbalanced ultracold Fermi gas in a large range of spin polarizations at low temperatures. We present results of microscopic calculations based on mean-field and density functional theory approaches, with no symmetry constraints. At low polarization values we predict the structure of the system as consisting of several spin-polarized droplets. As the polarization increases, the system self-organizes into a disordered structures similar to liquid crystals, and energetically they can compete with ordered structures such as grid-like domain walls. At higher polarizations the system starts to develop regularities that, in principle, can be called supersolid, where periodic density modulation and pairing correlations coexist. The robustness of the results has been checked with respect to temperature effects, dimensionality, and the presence of a trapping potential. Dynamical stability has also been investigated. △ Less

Submitted 18 March, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: Ancillary files contain reproducibility packs of numerical calculations

Journal ref: New J. Phys. 25, 033013 (2023)

Dissipation Mechanisms in Fermionic Josephson Junction

Authors: Gabriel Wlazłowski, Klejdja Xhani, Marek Tylutki, Nikolaos P. Proukakis, Piotr Magierski

Abstract: We characterize numerically the dominant dynamical regimes in a superfluid ultracold fermionic Josephson junction. Beyond the coherent Josephson plasma regime, we discuss the onset and physical mechanism of dissipation due to the superflow exceeding a characteristic speed, and provide clear evidence distinguishing its physical mechanism across the weakly- and strongly-interacting limits, despite q… ▽ More We characterize numerically the dominant dynamical regimes in a superfluid ultracold fermionic Josephson junction. Beyond the coherent Josephson plasma regime, we discuss the onset and physical mechanism of dissipation due to the superflow exceeding a characteristic speed, and provide clear evidence distinguishing its physical mechanism across the weakly- and strongly-interacting limits, despite qualitative dynamics of global characteristics being only weakly sensitive to the operating dissipative mechanism. Specifically, dissipation in the strongly interacting regime occurs through the phase-slippage process, caused by the emission and propagation of quantum vortices, and sound waves -- similar to the Bose-Einstein condensation limit. Instead, in the weak interaction limit, the main dissipative channel arises through the pair-breaking mechanism. △ Less

Submitted 13 January, 2023; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: manuscript, supplementary material, ancillary files containing reproducibility packs

Journal ref: Phys. Rev. Lett. 130, 023003 (2023)

Dissipative Dynamics of Quantum Vortices in Fermionic Superfluid

Authors: Andrea Barresi, Antoine Boulet, Piotr Magierski, Gabriel Wlazłowski

Abstract: In a recent article, Kwon et al. [Nature (London) {\bf 600}, 64 (2021)] revealed nonuniversal dissipative dynamics of quantum vortices in a fermionic superfluid. The enhancement of the dissipative process is pronounced for the Bardeen-Cooper-Schrieffer interaction regime, and it was suggested that the effect is due to the presence of quasiparticles localized inside the vortex core. We test this hy… ▽ More In a recent article, Kwon et al. [Nature (London) {\bf 600}, 64 (2021)] revealed nonuniversal dissipative dynamics of quantum vortices in a fermionic superfluid. The enhancement of the dissipative process is pronounced for the Bardeen-Cooper-Schrieffer interaction regime, and it was suggested that the effect is due to the presence of quasiparticles localized inside the vortex core. We test this hypothesis through numerical simulations with time-dependent density-functional theory: a fully microscopic framework with fermionic degrees of freedom. The results of fully microscopic calculations expose the impact of the vortex-bound states on dissipative dynamics in a fermionic superfluid. Their contribution is too weak to explain the experimental measurements, and we identify that thermal effects, giving rise to mutual friction between superfluid and the normal component, dominate the observed dynamics. △ Less

Submitted 26 January, 2023; v1 submitted 2 July, 2022; originally announced July 2022.

Comments: main paper with supplemental material; ancillary files containing reproducibility packs and example movies

Journal ref: Phys. Rev. Lett. 130, 043001 (2023)

Local energy density functional for superfluid Fermi gases from effective field theory

Authors: Antoine Boulet, Gabriel Wlazłowski, Piotr Magierski

Abstract: Over the past two decades, many studies in the Density Functional Theory context revealed new aspects and properties of strongly correlated superfluid quantum systems in numerous configurations that can be simulated in experiments. This was made possible by the generalization of the Local Density Approximation to superfluid systems by Bulgac in [Phys. Rev. C 65, 051305, (2002), Phys. Rev. A 76, 04… ▽ More Over the past two decades, many studies in the Density Functional Theory context revealed new aspects and properties of strongly correlated superfluid quantum systems in numerous configurations that can be simulated in experiments. This was made possible by the generalization of the Local Density Approximation to superfluid systems by Bulgac in [Phys. Rev. C 65, 051305, (2002), Phys. Rev. A 76, 040502, (2007)]. In the presented work, we propose an extension of the Superfluid Local Density Approximation systematically improvable and applicable to a large range of many-body quantum problems getting rid of the fitting procedures of the functional parameters. It turns out that only the knowledge of the density dependence of the quasi-particle properties, namely, the chemical potential, the effective mass, and the pairing gap function, are enough to obtain an explicit and accurate local functional of the densities without any adjustment a posteriori. This opens the way toward an Effective Field Theory formulation of the Density Functional Theory in the sense that we obtain a universal expansion of the functional parameters entering in the theory as a series in pairing gap function. Finally, we discuss possible applications of the developed approach allowing precise analysis of experimental observations. In that context, we focus our applications on the static structure properties of superfluid vortices. △ Less

Submitted 7 July, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Journal ref: Phys. Rev. A 106, 013306 (2022)

Pairing dynamics and solitonic excitations in collisions of medium-mass, identical nuclei

Authors: Piotr Magierski, Andrzej Makowski, Matthew C. Barton, Kazuyuki Sekizawa, Gabriel Wlazłowski

Abstract: We present results of collisions of $^{90}$Zr+$^{90}$Zr and $^{96}$Zr+$^{96}$Zr obtained within time-dependent density functional theory (TDDFT) extended to superfluid systems, known as time-dependent superfluid local density approximation (TDSLDA). We discuss qualitatively new features occurring in collisions of two superfluid nuclei at energies in the vicinity of the Coulomb barrier. We show tha… ▽ More We present results of collisions of $^{90}$Zr+$^{90}$Zr and $^{96}$Zr+$^{96}$Zr obtained within time-dependent density functional theory (TDDFT) extended to superfluid systems, known as time-dependent superfluid local density approximation (TDSLDA). We discuss qualitatively new features occurring in collisions of two superfluid nuclei at energies in the vicinity of the Coulomb barrier. We show that a \textit{solitonic excitation} -- an abrupt pairing phase distortion -- reported previously [P.~Magierski et al., Phys. Rev. Lett. \textbf{119}, 042501 (2017)], increases the barrier for capture generating effective repulsion between colliding nuclei. Moreover we demonstrate that pairing field leads to qualitatively different dynamics at the Coulomb barrier which manifests itself in a slower evolution of deformation towards a compact shape. Last but not least, we show that magnitude of pairing correlations can be dynamically enhanced after collision. We interpret it as a dynamically induced $U(1)$ symmetry breaking, which leads to large-amplitude oscillations of pairing field and bear similarity to the pairing Higgs mechanism. △ Less

Submitted 8 June, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 13 pages (including supplemental material), 6 figures

Journal ref: Phys. Rev. C 105, 064602, 2022

Properties of a quantum vortex in neutron matter

Authors: Daniel Pęcak, Nicolas Chamel, Piotr Magierski, Gabriel Wlazłowski

Abstract: We have studied systematically microscopic properties of a quantum vortex in neutron matter at finite temperatures and densities corresponding to different layers of the inner crust of a neutron star. To this end and in preparation of future simulations of the vortex dynamics, we have carried out fully self-consistent 3D Hartree-Fock-Bogoliubov calculations, using one of the latest nuclear energy-… ▽ More We have studied systematically microscopic properties of a quantum vortex in neutron matter at finite temperatures and densities corresponding to different layers of the inner crust of a neutron star. To this end and in preparation of future simulations of the vortex dynamics, we have carried out fully self-consistent 3D Hartree-Fock-Bogoliubov calculations, using one of the latest nuclear energy-density functionals from the Brussels-Montreal family, which has been developed specifically for applications to neutron superfluidity in neutron-star crusts. By analyzing the flow around the vortex, we have determined the effective radius relevant for the vortex filament model. We have also calculated the specific heat in the presence of the quantum vortex and have shown that it is substantially larger than for a uniform system at low temperatures. The low temperature limit of the specific heat has been identified as being determined by Andreev states inside the vortex core. We have shown that the specific heat in this limit does not scale linearly with temperature. The typical energy scale associated with Andreev states is defined by the minigap, which we have extracted for various neutron-matter densities. Our results suggest that vortices may be spin-polarized in the crust of magnetars. Finally, we have obtained a lower bound for the specific heat of a collection of vortices with given surface density, taking into account both the contributions from the vortex core states and from the hydrodynamic flow. △ Less

Submitted 5 November, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. C 104, 055801 (2021)

Dynamics of spin-polarized impurity in ultracold Fermi gas

Authors: Piotr Magierski, Buğra Tüzemen, Gabriel Wlazłowski

Abstract: We show that the motion of spin-polarized impurity (ferron) in ultracold atomic gas is characterized by a certain critical velocity which can be traced back to the amount of spin imbalance inside the impurity. We have calculated the effective mass of ferron in two dimensions. We show that the effective mass scales with the surface of the ferron. We discuss the impact of these findings; in particul… ▽ More We show that the motion of spin-polarized impurity (ferron) in ultracold atomic gas is characterized by a certain critical velocity which can be traced back to the amount of spin imbalance inside the impurity. We have calculated the effective mass of ferron in two dimensions. We show that the effective mass scales with the surface of the ferron. We discuss the impact of these findings; in particular, we demonstrate that ferrons become unstable in the vicinity of a vortex. △ Less

Submitted 14 September, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Journal ref: Phys. Rev. A 104, 033304 (2021)

Spin-polarized vortices with reversed circulation

Authors: Piotr Magierski, Gabriel Wlazłowski, Andrzej Makowski, Konrad Kobuszewski

Abstract: We present the analysis of the structure of fermionic vortices with the spin-polarized core from a weak coupling limit to the unitary regime. We show the mechanism for the generation of the {\it reversed circulation} in the vortex core induced by an excess of majority spin particles. We introduce the classification of the polarized vortices based on the number of Fermi circles where the minigap va… ▽ More We present the analysis of the structure of fermionic vortices with the spin-polarized core from a weak coupling limit to the unitary regime. We show the mechanism for the generation of the {\it reversed circulation} in the vortex core induced by an excess of majority spin particles. We introduce the classification of the polarized vortices based on the number of Fermi circles where the minigap vanishes. This provides a unique description of the vortex as one cannot smoothly map wave functions into one another corresponding to vortices differing by the number of Fermi circles. The effective mass of quasiparticles along the vortex core is analyzed and its role in the propagation of spin-polarization along the vortex line is discussed. △ Less

Submitted 30 September, 2022; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 11 pages, 8 figures

Journal ref: Phys. Rev. A 106, 033322 (2022)

Rotating quantum turbulence in the unitary Fermi gas

Authors: Khalid Hossain, Konrad Kobuszewski, Michael McNeil Forbes, Piotr Magierski, Kazuyuki Sekizawa, Gabriel Wlazłowski

Abstract: Quantized vortices carry the angular momentum in rotating superfluids, and are key to the phenomenon of quantum turbulence. Advances in ultra-cold atom technology enable quantum turbulence to be studied in regimes with both experimental and theoretical control, unlike the original contexts of superfluid helium experiments. While much work has been performed with bosonic systems, detailed studies o… ▽ More Quantized vortices carry the angular momentum in rotating superfluids, and are key to the phenomenon of quantum turbulence. Advances in ultra-cold atom technology enable quantum turbulence to be studied in regimes with both experimental and theoretical control, unlike the original contexts of superfluid helium experiments. While much work has been performed with bosonic systems, detailed studies of fermionic quantum turbulence are nascent, despite wide applicability to other contexts such as rotating neutron stars. In this paper, we present the first large-scale study of quantum turbulence in rotating fermionic superfluids using an accurate orbital based time-dependent density functional theory (DFT) called the superfluid local density approximation (SLDA). We identify two different modes of turbulent decay in the dynamical equilibration of a rotating fermionic superfluid, and contrast these results with a computationally simpler orbital-free DFT, which we find can qualitatively reproduce these decay mechanisms if dissipation is explicitly included. These results demonstrate that one-body dissipation mechanisms intrinsic to fermionic superfluids play a key role differentiating fermionic from bosonic turbulence, but also suggest that simpler orbital-free theories may be corrected so that these more efficient techniques can be used to model extended physical systems such as neutron superfluids in neutron stars. △ Less

Submitted 14 January, 2022; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 13 pages, 8 figures. Accompanying movies can be found at http://wlazlowski.fizyka.pw.edu.pl/supplement/sm-rotating-turbulence/

Report number: INT-PUB-20-040

Journal ref: Phys. Rev. A 105, 013304 (2022)

Future of Nuclear Fission Theory

Authors: Michael Bender, Remi Bernard, George Bertsch, Satoshi Chiba, Jacek Dobaczewski, Noel Dubray, Samuel Giuliani, Kouichi Hagino, Denis Lacroix, Zhipan Li, Piotr Magierski, Joachim Maruhn, Witold Nazarewicz, Junchen Pei, Sophie Peru, Nathalie Pillet, Jorgen Randrup, David Regnier, Paul-Gerhard Reinhard, Luis Robledo, Wouter Ryssens, Jhilam Sadhukhan, Guillaume Scamps, Nicolas Schunck, Cedric Simenel , et al. (8 additional authors not shown)

Abstract: There has been much recent interest in nuclear fission, due in part to a new appreciation of its relevance to astrophysics, stability of superheavy elements, and fundamental theory of neutrino interactions. At the same time, there have been important developments on a conceptual and computational level for the theory. The promising new theoretical avenues were the subject of a workshop held at the… ▽ More There has been much recent interest in nuclear fission, due in part to a new appreciation of its relevance to astrophysics, stability of superheavy elements, and fundamental theory of neutrino interactions. At the same time, there have been important developments on a conceptual and computational level for the theory. The promising new theoretical avenues were the subject of a workshop held at the University of York in October 2019; this report summarises its findings and recommendations. △ Less

Submitted 2 November, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 86 LaTeX pages, 3 figures, 331 references, version published as Topical Review in Journal of Physics G

Journal ref: J. Phys. G: Nucl. Part. Phys. 47 (2020) 113002

Properties of spin-polarized impurities -- ferrons, in the unitary Fermi gas

Authors: Buğra Tüzemen, Paweł Kukliński, Piotr Magierski, Gabriel Wlazłowski

Abstract: A new excitation mode has been predicted to exist in the unitary Fermi gas. It has a form of a spin-polarized impurity, which was dubbed as ferron. It is characterized by a closed nodal surface of the pairing field surrounding a partially spin-polarized superfluid region, where the phase differs by $π$. In this paper, we discuss the effect of temperature on the generation of the ferron and the adi… ▽ More A new excitation mode has been predicted to exist in the unitary Fermi gas. It has a form of a spin-polarized impurity, which was dubbed as ferron. It is characterized by a closed nodal surface of the pairing field surrounding a partially spin-polarized superfluid region, where the phase differs by $π$. In this paper, we discuss the effect of temperature on the generation of the ferron and the adiabaticity of the spin-polarizing potential together with ferron's ground state properties. △ Less

Submitted 11 March, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Journal ref: Acta Physics Polonica B 51, 595 (2020)

Pairing dynamics in low energy nuclear collisions

Authors: M. C. Barton, S. Jin, P. Magierski, K. Sekizawa, G. Wlazłowski, A. Bulgac

Abstract: Superfluidity is a generic feature of various quantum systems at low temperatures and it is in particular important for the description of dynamics of low energy nuclear reactions. The time-dependent density functional theory (TDDFT) is, to date, the only microscopic method which takes into account in a consistent way far from equilibrium dynamics of pairing field and single-particle degrees of fr… ▽ More Superfluidity is a generic feature of various quantum systems at low temperatures and it is in particular important for the description of dynamics of low energy nuclear reactions. The time-dependent density functional theory (TDDFT) is, to date, the only microscopic method which takes into account in a consistent way far from equilibrium dynamics of pairing field and single-particle degrees of freedom. The local version of TDDFT, so called TDSLDA, is particularly useful for the description of nuclear reactions and is well suited for leadership class computers of hybrid (CPU+GPU) architecture. The preliminary results obtained for collisions involving both medium-mass and heavy nuclei at the energies around the Coulomb barrier are presented. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 6 Pages, 1 Figure, XXXVI Mazurian Lakes Conference

Journal ref: Acta Physics Polonica B 51, 605 (2020)

Spin-polarized droplets in the unitary Fermi gas

Authors: Piotr Magierski, Buğra Tüzemen, Gabriel Wlazłowski

Abstract: We demonstrate the existence of a new type of spatially localized excitations in the unitary Fermi gas: spin polarized droplets with a peculiar internal structure involving the abrupt change of the pairing phase at the surface of the droplet. It resembles the structure of the Josephson-$π$ junction occurring when a slice of a ferromagnet is sandwiched between two superconductors. The stability of… ▽ More We demonstrate the existence of a new type of spatially localized excitations in the unitary Fermi gas: spin polarized droplets with a peculiar internal structure involving the abrupt change of the pairing phase at the surface of the droplet. It resembles the structure of the Josephson-$π$ junction occurring when a slice of a ferromagnet is sandwiched between two superconductors. The stability of the impurity is enhanced by the mutual interplay between the polarization effects and the pairing field, resulting in an exceptionally long-lived state. The prospects for its realization in experiment are discussed. △ Less

Submitted 17 September, 2019; v1 submitted 1 November, 2018; originally announced November 2018.

Comments: Manuscript with included Supplementary Material. Published version. Accompanying movies are available at http://wlazlowski.fizyka.pw.edu.pl/supplement/2/

Journal ref: Phys. Rev. A 100, 033613 (2019)

Pairing dynamics and time dependent density functional theory

Authors: P. Magierski, J. Grineviciute, K. Sekizawa

Abstract: We discuss issues related to pairing dynamics in nuclear large amplitude collective motion. The examples of effects which are not properly described within BCS theory are presented. In the second part we review properties of TDDFT and in particular we discuss the time-dependent local density approximation (TDSLDA) starting from the stationary action principle. We discuss issues related to pairing dynamics in nuclear large amplitude collective motion. The examples of effects which are not properly described within BCS theory are presented. In the second part we review properties of TDDFT and in particular we discuss the time-dependent local density approximation (TDSLDA) starting from the stationary action principle. △ Less

Submitted 2 July, 2018; v1 submitted 26 December, 2017; originally announced December 2017.

Comments: 11 pages

Journal ref: Acta Phys. Pol. B49, 281 (2018)

Suppressed solitonic cascade in spin-imbalanced superfluid Fermi gas

Authors: Gabriel Wlazłowski, Kazuyuki Sekizawa, Maciej Marchwiany, Piotr Magierski

Abstract: Cold atoms experiments offer invaluable information on superfluid dynamics, including decay cascades of topological defects. While the cascade properties are well established for Bose systems, our understanding of their behavior in Fermi counterparts is very limited, in particular in spin-imbalanced systems, where superfluid (paired) and normal (unpaired) particles naturally coexist giving rise to… ▽ More Cold atoms experiments offer invaluable information on superfluid dynamics, including decay cascades of topological defects. While the cascade properties are well established for Bose systems, our understanding of their behavior in Fermi counterparts is very limited, in particular in spin-imbalanced systems, where superfluid (paired) and normal (unpaired) particles naturally coexist giving rise to complex spatial structure of the atomic cloud. Here we show, based on a newly developed microscopic approach, that the decay cascades of topological defects are dramatically modified by the spin-polarization. We demonstrate that decay cascades end up at different stages: "dark soliton", "vortex ring" or "vortex line", depending on the polarization. We reveal that it is caused by sucking of unpaired particles into the soliton's internal structure. As a consequence vortex reconnections are hindered and we anticipate that quantum turbulence phenomenon can be significantly affected, indicating new physics induced by polarization effects. △ Less

Submitted 26 June, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: Published version. Manuscript with included Supplementary Material. Movies available at http://www.if.pw.edu.pl/~gabrielw/supplement/1/supplement.html

Journal ref: Phys. Rev. Lett. 120, 253002 (2018)

Accuracy of fission dynamics within the time dependent superfluid local density approximation

Authors: J. Grineviciute, P. Magierski, A. Bulgac, S. Jin, I. Stetcu

Abstract: 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. 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. △ Less

Submitted 20 November, 2017; v1 submitted 6 November, 2017; originally announced November 2017.

Comments: Contribution to Proceedings of XXXV Mazurian Lakes Conference on Physics, Piaski, Poland, September 3-9, 2017

Journal ref: Acta Physica Polonica B 49 (2018) 591

Solitonic excitations in collisions of superfluid nuclei: a qualitatively new phenomenon distinct from the Josephson effect

Authors: Kazuyuki Sekizawa, Gabriel Wlazłowski, Piotr Magierski

Abstract: Recently, we have reported a novel role of pairing in low-energy heavy ion reactions at energies above the Coulomb barrier, which may have a detectable impact on reaction outcomes, such as the kinetic energy of fragments and the fusion cross section [arXiv:1611.10261, arXiv:1702.00069]. The phenomenon mimics the one studied experimentally with ultracold atomic gases, where two clouds of fermionic… ▽ More Recently, we have reported a novel role of pairing in low-energy heavy ion reactions at energies above the Coulomb barrier, which may have a detectable impact on reaction outcomes, such as the kinetic energy of fragments and the fusion cross section [arXiv:1611.10261, arXiv:1702.00069]. The phenomenon mimics the one studied experimentally with ultracold atomic gases, where two clouds of fermionic superfluids with different phases of the pairing fields are forced to merge, inducing various excitation modes of the pairing field. Although it originates from the phase difference of the pairing fields, the physics behind it is markedly different from the so-called Josephson effect. In this short contribution, we will briefly outline the results discussed in our recent papers and explain relations with the field of ultracold atomic gases. △ Less

Submitted 13 May, 2017; originally announced May 2017.

Comments: 5 pages, 3 figures, Account of an invited talk given at the International Conference on heavy-ion collisions at near-barrier energies (FUSION17), Hobart, Tasmania, Australia, Feb. 20-24, 2017

Journal ref: EPJ Web of Conf. 163, 00051 (2017)

Nuclear Fission: from more phenomenology and adjusted parameters to more fundamental theory and increased predictive power

Authors: A. Bulgac, S. Jin, P. Magierski, K. Roche, N. Schunck, I. Stetcu

Abstract: Two major recent developments in theory and computational resources created the favorable conditions for achieving a microscopic description of nuclear fission almost eighty years after its discovery in 1939 by Hahn and Strassmann (1930). The first major development was in theory, the extension of the Time-Dependent Density Functional Theory (TDDFT) to superfluid fermion systems. The second develo… ▽ More Two major recent developments in theory and computational resources created the favorable conditions for achieving a microscopic description of nuclear fission almost eighty years after its discovery in 1939 by Hahn and Strassmann (1930). The first major development was in theory, the extension of the Time-Dependent Density Functional Theory (TDDFT) to superfluid fermion systems. The second development was in computing, the emergence of powerful enough supercomputers capable of solving the complex systems of equations describing the time evolution in three dimensions without any restrictions of hundreds of strongly interacting nucleons. Even though the available nuclear energy density functionals (NEDFs) are phenomenological still, their accuracy is improving steadily and the prospects of being able to perform calculations of the nuclear fission dynamics and to predict many properties of the fission fragments, otherwise not possible to extract from experiments, are within reach, all without making recourse anymore to uncontrollable assumptions and simplified phenomenological models. △ Less

Submitted 28 April, 2017; originally announced May 2017.

Comments: 6 pages, account of invited talk given at FUSION17, Hobart, Tasmania, February 20-24, 2017

Journal ref: EPJ Web of Conferences, 163, 00007 (2017)

Microscopic Theory of Nuclear Fission

Authors: Aurel Bulgac, Shi Jin, Piotr Magierski, Kenneth J. Roche, Ionel Stetcu

Abstract: We describe the fission dynamics of $^{240}$Pu within an implementation of the Density Functional Theory (DFT) extended to superfluid systems and real-time dynamics. We demonstrate the critical role played by the pairing correlations, which even though are not the driving force in this complex dynamics, are providing the essential lubricant, without which the nuclear shape evolution would come to… ▽ More We describe the fission dynamics of $^{240}$Pu within an implementation of the Density Functional Theory (DFT) extended to superfluid systems and real-time dynamics. We demonstrate the critical role played by the pairing correlations, which even though are not the driving force in this complex dynamics, are providing the essential lubricant, without which the nuclear shape evolution would come to a screeching halt. The evolution is found to be much slower than previously expected in this fully non-adiabatic treatment of nuclear dynamics, where there are no symmetry restrictions and all collective degrees of freedom (CDOF) are allowed to participate in the dynamics. △ Less

Submitted 5 April, 2017; v1 submitted 3 April, 2017; originally announced April 2017.

Comments: 8 pages, 3 figures, Account of invited talk given at The 26th International Nuclear Physics Conference, 11-16 September, 2016, Adelaide, Australia (corrected a couple of silly typos)

Journal ref: PoS(INPC2016)225

Induced fission of 240Pu

Authors: A. Bulgac, S. Jin, P. Magierski, K. J. Roche, I. Stetcu

Abstract: We study the fission dynamics of 240Pu within an implementation of the Density Functional Theory (DFT) extended to superfluid systems and real-time dynamics. We demonstrate the critical role played by the pairing correlations. The evolution is found to be much slower than previously expected in this fully non-adiabatic treatment of nuclear dynamics, where there are no symmetry restrictions and all… ▽ More We study the fission dynamics of 240Pu within an implementation of the Density Functional Theory (DFT) extended to superfluid systems and real-time dynamics. We demonstrate the critical role played by the pairing correlations. The evolution is found to be much slower than previously expected in this fully non-adiabatic treatment of nuclear dynamics, where there are no symmetry restrictions and all collective degrees of freedom (CDOF) are allowed to participate in the dynamics. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 8 pages, 4 figures, talk given at The 6th International Conference on Fission and Properties of Neutron-Rich Nuclei, Sanibel Island, Florida, November 6-2 (2016)

Report number: NT@UW-17-03

Journal ref: pages 564-572, in Fission and Properties of Neutron-Rich-Nuclei, eds. J.H. Hamilton, A.V. Ramayya, and P. Talou, World Scientific, Singapore (2017)

Solitonic Excitations in Collisions Of Superfluid Nuclei

Authors: Kazuyuki Sekizawa, Piotr Magierski, Gabriel Wlazłowski

Abstract: We investigate the role of the pairing field dynamics in low-energy heavy ion reactions within the nuclear time-dependent density functional theory extended to superfluid systems. Recently, we have reported on unexpectedly large effects associated with the relative phase of the pairing field of colliding nuclei on the reaction outcomes, such as the total kinetic energy and the fusion cross section… ▽ More We investigate the role of the pairing field dynamics in low-energy heavy ion reactions within the nuclear time-dependent density functional theory extended to superfluid systems. Recently, we have reported on unexpectedly large effects associated with the relative phase of the pairing field of colliding nuclei on the reaction outcomes, such as the total kinetic energy and the fusion cross section [P. Magierski, K. Sekizawa, and G. Wlazłowski, arXiv:1611.10261 [nucl-th]]. We have elucidated that the effects are due to creation of a "domain wall" or a "solitonic structure" of the pairing field in the neck region, which hinders energy dissipation as well as the neck formation, leading to significant changes of the reaction dynamics. The situation nicely mimics the one extensively studied experimentally with ultracold atomic gases, where two clouds of superfluid atoms possessing different phases of the pairing field are forced to merge, creating various topological excitations, quantum vortices and solitons, as well as Josephson currents. In this paper, we present unpublished results for a lighter system, namely, $^{44}$Ca+$^{44}$Ca. It is shown that the pairing effects on the fusion hindrance are rather small in lighter systems, due to a strong tendency towards fusion, which is consistent with an earlier study. △ Less

Submitted 31 January, 2017; originally announced February 2017.

Comments: 8 pages, 2 figures, Proceedings of the 26th International Nuclear Physics Conference (INPC2016), Adelaide, Australia, Sep. 11-16, 2016

Journal ref: PoS(INPC2016)214

Novel Role of Superfluidity in Low-Energy Nuclear Reactions

Authors: Piotr Magierski, Kazuyuki Sekizawa, Gabriel Wlazłowski

Abstract: We demonstrate, within symmetry unrestricted time-dependent density functional theory, the existence of new effects in low-energy nuclear reactions which originate from superfluidity. The dynamics of the pairing field induces solitonic excitations in the colliding nuclear systems, leading to qualitative changes in the reaction dynamics. The solitonic excitation prevents collective energy dissipati… ▽ More We demonstrate, within symmetry unrestricted time-dependent density functional theory, the existence of new effects in low-energy nuclear reactions which originate from superfluidity. The dynamics of the pairing field induces solitonic excitations in the colliding nuclear systems, leading to qualitative changes in the reaction dynamics. The solitonic excitation prevents collective energy dissipation and effectively suppresses fusion cross section. We demonstrate how the variations of the total kinetic energy of the fragments can be traced back to the energy stored in the superfluid junction of colliding nuclei. Both contact time and scattering angle in non-central collisions are significantly affected. The modification of the fusion cross section and possibilities for its experimental detection are discussed. △ Less

Submitted 25 July, 2017; v1 submitted 30 November, 2016; originally announced November 2016.

Comments: Manuscript with included Supplemental online material, published version

Journal ref: Phys. Rev. Lett. 119, 042501 (2017)

Microscopic Calculations of Vortex-Nucleus Interaction in the Neutron Star Crust

Authors: Kazuyuki Sekizawa, Gabriel Wlazłowski, Piotr Magierski, Aurel Bulgac, Michael McNeil Forbes

Abstract: We investigate the dynamics of a quantized vortex and a nuclear impurity immersed in a neutron superfluid within a fully microscopic time-dependent three-dimensional approach. The magnitude and even the sign of the force between the quantized vortex and the nuclear impurity have been a matter of debate for over four decades. We determine that the vortex and the impurity repel at neutron densities,… ▽ More We investigate the dynamics of a quantized vortex and a nuclear impurity immersed in a neutron superfluid within a fully microscopic time-dependent three-dimensional approach. The magnitude and even the sign of the force between the quantized vortex and the nuclear impurity have been a matter of debate for over four decades. We determine that the vortex and the impurity repel at neutron densities, 0.014 fm$^{-3}$ and 0.031 fm$^{-3}$, which are relevant to the neutron star crust and the origin of glitches, while previous calculations have concluded that the force changes its sign between these two densities and predicted contradictory signs. The magnitude of the force increases with the density of neutron superfluid, while the magnitude of the pairing gap decreases in this density range. △ Less

Submitted 11 September, 2016; originally announced September 2016.

Comments: 4 pages, 2 figures, Talk given at the 14th International Symposium on "Nuclei in the Cosmos" (NIC-XIV), June 19-24, 2016, Toki Messe, Niigata, Japan

Report number: NT@UW-16-11

Journal ref: JPS Conf. Proc. 14, 010807 (2017)

Vortex pinning and dynamics in the neutron star crust

Authors: Gabriel Wlazłowski, Kazuyuki Sekizawa, Piotr Magierski, Aurel Bulgac, Michael McNeil Forbes

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 2016; v1 submitted 8 June, 2016; originally announced June 2016.

Comments: supplemental materials included

Report number: NT@UW-16-06

Journal ref: Phys. Rev. Lett. 117, 232701 (2016)

Nuclear Reactions and Superfluid Time Dependent Density Functional Theory

Authors: Piotr Magierski

Abstract: The extension of Time Dependent Density Functional Theory (TDDFT) to superfluid systems is discussed in the context of nuclear reactions and large amplitude collective motion. The extension of Time Dependent Density Functional Theory (TDDFT) to superfluid systems is discussed in the context of nuclear reactions and large amplitude collective motion. △ Less

Submitted 23 June, 2019; v1 submitted 7 June, 2016; originally announced June 2016.

Comments: 15 pages, to be published in the ebook "Progress of time-dependent nuclear reaction theory" (Betham Science Publishers 2019) honoring Prof. Joachim Maruhn's retirement

Journal ref: Frontiers in Nuclear and Particle Physics, 2019, Vol. 2, 57-71

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

Authors: Aurel Bulgac, Piotr Magierski, Kenneth J. Roche, Ionel Stetcu

Abstract: 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… ▽ More 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). △ Less

Submitted 25 March, 2016; v1 submitted 2 November, 2015; originally announced November 2015.

Comments: updated version in agreement with published paper, updated and new references, 7 pages, additional material at http://faculty.washington.edu/bulgac/Pu240/

Report number: NT@UW-15-14, PNNL-SA-114173, LA-UR-15-28092

Journal ref: Phys. Rev. Lett. 116, 122504 (2016)

Relativistic Coulomb excitation within Time Dependent Superfluid Local Density Approximation

Authors: I. Stetcu, C. Bertulani, A. Bulgac, P. Magierski, K. J. Roche

Abstract: Within the framework of the unrestricted time-dependent density functional theory, we present for the first time an analysis of the relativistic Coulomb excitation of the heavy deformed open shell nucleus $^{238}$U. The approach is based on Superfluid Local Density Approximation (SLDA) formulated on a spatial lattice that can take into account coupling to the continuum, enabling self-consistent st… ▽ More Within the framework of the unrestricted time-dependent density functional theory, we present for the first time an analysis of the relativistic Coulomb excitation of the heavy deformed open shell nucleus $^{238}$U. The approach is based on Superfluid Local Density Approximation (SLDA) formulated on a spatial lattice that can take into account coupling to the continuum, enabling self-consistent studies of superfluid dynamics of any nuclear shape. We have computed the energy deposited in the target nucleus as a function of the impact parameter, finding it to be significantly larger than the estimate using the Goldhaber-Teller model. The isovector giant dipole resonance, the dipole pygmy resonance and giant quadrupole modes were excited during the process. The one body dissipation of collective dipole modes is shown to lead a damping width $Γ_\downarrow \approx 0.4$ MeV and the number of pre-equilibrium neutrons emitted has been quantified. △ Less

Submitted 13 January, 2015; v1 submitted 11 March, 2014; originally announced March 2014.

Comments: 5 pages, 3 figures and 21 pages, 17 figures of supplementary material; published version

Report number: LA-UR-14-21210, NT@UW-14-06

Journal ref: Phys. Rev. Lett. 114, 012701 (2015)

Properties of N=90 Isotones within the Mean Field Perspective

Authors: E. Ganioğlu, R. Wyss, P. Magierski

Abstract: In recent years, the N=90 isotones have been investigated to a large extent in relation to studies of quantum phase transitions. In this paper, we use the mean field approach with pairing-deformation self consistent Total Routhian Surface (TRS) calculations to study the N=90 isotones and neigh- bouring nuclei. The important probes, such as moments of inertia, quadrupole moments, the energy ratio o… ▽ More In recent years, the N=90 isotones have been investigated to a large extent in relation to studies of quantum phase transitions. In this paper, we use the mean field approach with pairing-deformation self consistent Total Routhian Surface (TRS) calculations to study the N=90 isotones and neigh- bouring nuclei. The important probes, such as moments of inertia, quadrupole moments, the energy ratio of E(4+1)/E(2+1), octupole and hexadecapole degrees of freedom are considered and the cal- culated results are compared with the available experimental data. From a microscopic point of view, the N=90 isotones characterize the onset of the deformed region and are very well described by mean field calculations. The results are compared with those from other studies in beyond mean-field approximations. Shape coexistence phenomena in the region of interest are discussed. △ Less

Submitted 2 September, 2013; originally announced September 2013.

Comments: 7 pages, 10 figures

Journal ref: Phys. Rev. C 89, 014311,14 January 2014

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

Authors: Gabriel Wlazłowski, Piotr Magierski, Aurel Bulgac, Kenneth J. Roche

Abstract: 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 boun… ▽ More 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. △ Less

Submitted 31 July, 2013; v1 submitted 8 April, 2013; originally announced April 2013.

Report number: NT@UW-13-12

Journal ref: Phys. Rev. A 88, 013639 (2013)

Cooper pairing above the critical temperature in a unitary Fermi gas

Authors: Gabriel Wlazłowski, Piotr Magierski, Joaquín E. Drut, Aurel Bulgac, Kenneth J. Roche

Abstract: 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 transitio… ▽ More 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. △ Less

Submitted 26 February, 2013; v1 submitted 6 December, 2012; originally announced December 2012.

Comments: 9 pages, 8 figures; supplemental materials included; published version

Journal ref: Phys. Rev. Lett. 110, 090401 (2013)

Shear Viscosity of a Unitary Fermi Gas

Authors: Gabriel Wlazłowski, Piotr Magierski, Joaquín E. Drut

Abstract: We present the first ab initio determination of the shear viscosity eta 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 shear viscosity to entropy density ratio eta/s. The minimum of eta/s appears to be located above the critical temperature for the superfluid-to-normal pha… ▽ More We present the first ab initio determination of the shear viscosity eta 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 shear viscosity to entropy density ratio eta/s. The minimum of eta/s appears to be located above the critical temperature for the superfluid-to-normal phase transition with the most probable value being eta/s approx 0.2 hbar/kB, which almost saturates the Kovtun-Son-Starinets universal value hbar/(4 pi kB). △ Less

Submitted 12 July, 2012; v1 submitted 1 April, 2012; originally announced April 2012.

Comments: 8 pages, 9 figures; supplemental materials included

Report number: LA-UR-12-20206/LA-UR-12-20207; NT@UW-12-04

Journal ref: Phys. Rev. Lett. 109:020406, 2012

The Equation of State of the Unitary Fermi Gas: An Update on Lattice Calculations

Authors: Joaquín E. Drut, Timo A. Lähde, Gabriel Wlazłowski, Piotr Magierski

Abstract: The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below T/eF ~ 0.20, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) resul… ▽ More The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below T/eF ~ 0.20, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data for the energy, chemical potential, and density. These AFQMC results have been obtained using methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the context of lattice QCD. △ Less

Submitted 7 May, 2012; v1 submitted 21 November, 2011; originally announced November 2011.

Comments: 4 pages, 3 figures

Report number: LA-UR-11-12070; INT-PUB-11-058

Journal ref: Phys. Rev. A 85, 051601(R) (2012)

LINPRO: linear inverse problem library for data contaminated by statistical noise

Authors: Piotr Magierski, Gabriel Wlazlowski

Abstract: The library LINPRO which provides solution to the linear inverse problem for data contaminated by a statistical noise is presented. The library makes use of two methods: Maximum Entropy Method and Singular Value Decomposition. As an example it has been applied to perform an analytic continuation of the imaginary time propagator obtained within the Quantum Monte Carlo method. The library LINPRO which provides solution to the linear inverse problem for data contaminated by a statistical noise is presented. The library makes use of two methods: Maximum Entropy Method and Singular Value Decomposition. As an example it has been applied to perform an analytic continuation of the imaginary time propagator obtained within the Quantum Monte Carlo method. △ Less

Submitted 27 June, 2012; v1 submitted 25 October, 2011; originally announced October 2011.

Comments: The associated computer program is available at: http://tja.if.pw.edu.pl/linpro/

Journal ref: Comput. Phys. Commun. 183 (2012) 2264-2271

Isovector Giant Dipole Resonance from the 3D Time-Dependent Density Functional Theory for Superfluid Nuclei

Authors: I. Stetcu, A. Bulgac, P. Magierski, K. J. Roche

Abstract: A fully symmetry unrestricted Time-Dependent Density Functional Theory extended to include pairing correlations is used to calculate properties of the isovector giant dipole resonances of the deformed open-shell nuclei 172Yb (axially deformed), 188Os (triaxially deformed), and 238U (axially deformed), and to demonstrate good agreement with experimental data on nuclear photo-absorption cross-sectio… ▽ More A fully symmetry unrestricted Time-Dependent Density Functional Theory extended to include pairing correlations is used to calculate properties of the isovector giant dipole resonances of the deformed open-shell nuclei 172Yb (axially deformed), 188Os (triaxially deformed), and 238U (axially deformed), and to demonstrate good agreement with experimental data on nuclear photo-absorption cross-sections for two different Skyrme force parametrizations of the energy density functional: SkP and SLy4. △ Less

Submitted 21 November, 2011; v1 submitted 15 August, 2011; originally announced August 2011.

Comments: 5 pages, 3 figures, published version

Journal ref: Phys. Rev. C 84, 051309(R) (2011)

Onset of a Pseudogap Regime in Ultracold Fermi Gases

Authors: Piotr Magierski, Gabriel Wlazłowski, Aurel Bulgac

Abstract: 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 i… ▽ More 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. △ Less

Submitted 3 October, 2011; v1 submitted 22 March, 2011; originally announced March 2011.

Comments: 5 pages, 4 figures + 5 pages of supplementary online material

Journal ref: Phys. Rev. Lett. 107: 145304,2011

Generation and Dynamics of Quantized Vortices in a Unitary Fermi Superfluid

Authors: Aurel Bulgac, Yuan-Lung, Luo, Piotr Magierski, Kenneth J. Roche, Yongle Yu

Abstract: 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 a… ▽ More 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. △ Less

Submitted 27 November, 2010; originally announced November 2010.

Comments: 5 pages main text + 20 pages supplemental material. 32 accompanying videos can be viewed and downloaded from http://www.phys.washington.edu/groups/qmbnt/UFG/ (Please follow instructions in supplemental material.)

Journal ref: Science 332:1288,2011

Superfluid properties of dilute neutron matter

Authors: Gabriel Wlazlowski, Piotr Magierski

Abstract: We report results of fully non-perturbative calculations, based on Auxiliary Field Quantum Monte Carlo (AFQMC) approach, for the dilute neutron matter at the density $ρ=0.003\fm^{-3}$. Fundamental quantities which characterize the superfluid state: the single particle energy gap $Δ(T)$, and the critical temperature $\Tc$ have been determined. The large value of $Δ(0)/\Tc\approx 3.2$ indicates that… ▽ More We report results of fully non-perturbative calculations, based on Auxiliary Field Quantum Monte Carlo (AFQMC) approach, for the dilute neutron matter at the density $ρ=0.003\fm^{-3}$. Fundamental quantities which characterize the superfluid state: the single particle energy gap $Δ(T)$, and the critical temperature $\Tc$ have been determined. The large value of $Δ(0)/\Tc\approx 3.2$ indicates that the system is not a BCS-type superfluid at low temperatures. △ Less

Submitted 16 November, 2010; originally announced November 2010.

Comments: Proceeding: XVII Nuclear Physics Workshop in Kazimierz Dolny, Poland 2010. Accepted to publish in the International Journal of Modern Physics E, vol. 20

Journal ref: Int. J. Mod. Phys. E20, 569 (2011)

The Unitary Fermi Gas: From Monte Carlo to Density Functionals

Authors: Aurel Bulgac, Michael McNeil Forbes, Piotr Magierski

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2012; v1 submitted 23 August, 2010; originally announced August 2010.

Comments: 73 pages, 20 figures. Chapter 9 in "The BCS-BEC Crossover and the Unitary Fermi Gas" edited by W. Zwerger (Springer, 2012). Updated to match published version

Report number: LA-UR 10-05509; INT-PUB-10-037; NW@UW-10-18

Journal ref: Lecture Notes in Physics, Vol. 836, Chap. 9, p.305-373 (2012)

Quantum Monte Carlo study of dilute neutron matter at finite temperatures

Authors: Gabriel Wlazlowski, Piotr Magierski

Abstract: We report results of fully non-perturbative, Path Integral Monte Carlo (PIMC) calculations for dilute neutron matter. The neutron-neutron interaction in the s channel is parameterized by the scattering length and the effective range. We calculate the energy and the chemical potential as a function of temperature at the density $\dens=0.003\fm^{-3}$. The critical temperature $\Tc$ for the superflui… ▽ More We report results of fully non-perturbative, Path Integral Monte Carlo (PIMC) calculations for dilute neutron matter. The neutron-neutron interaction in the s channel is parameterized by the scattering length and the effective range. We calculate the energy and the chemical potential as a function of temperature at the density $\dens=0.003\fm^{-3}$. The critical temperature $\Tc$ for the superfluid-normal phase transition is estimated from the finite size scaling of the condensate fraction. At low temperatures we extract the spectral weight function $A(p,ω)$ from the imaginary time propagator using the methods of maximum entropy and singular value decomposition. We determine the quasiparticle spectrum, which can be accurately parameterized by three parameters: an effective mass $m^*$, a mean-field potential $U$, and a gap $Δ$. Large value of $Δ/\Tc$ indicates that the system is not a BCS-type superfluid at low temperatures. △ Less

Submitted 1 February, 2011; v1 submitted 2 December, 2009; originally announced December 2009.

Comments: 4 pages, 3 figures

Journal ref: Phys.Rev.C83:012801,2011

Effective Range Expansion for the Interaction Defined on the Lattice

Authors: Gabriel Wlazłowski, Piotr Magierski

Abstract: The relation between the interaction parameters for fermions on the spatial lattice and the two-body $T$ matrix is discussed. The presented method allows determination of the interaction parameters through the relatively simple computational scheme which include the effect of finite lattice spacing. In particular the relation between the interaction parameters and the effective range expansion p… ▽ More The relation between the interaction parameters for fermions on the spatial lattice and the two-body $T$ matrix is discussed. The presented method allows determination of the interaction parameters through the relatively simple computational scheme which include the effect of finite lattice spacing. In particular the relation between the interaction parameters and the effective range expansion parameters is derived in the limit of large lattices. △ Less

Submitted 7 November, 2009; originally announced November 2009.

Comments: Proceedings from XVI Nuclear Physics Workshop in Kazimierz Dolny, Poland. Accepted to publish in the International Journal of Modern Physics E, vol. 19

Journal ref: Int.J.Mod.Phys.E19:781-786,2010

Quantum Monte Carlo method applied to strongly correlated dilute fermi gases with finite effective range

Authors: G. Wlazlowski, P. Magierski

Abstract: We discuss the Auxiliary Field Quantum Monte Carlo (AFQMC) method applied to dilute neutron matter at finite temperatures. We formulate the discrete Hubbard-Stratonovich transformation for the interaction with finite effective range which is free from the sign problem. The AFQMC results are compared with those obtained from exact diagonalization for a toy model. Preliminary calculations of energ… ▽ More We discuss the Auxiliary Field Quantum Monte Carlo (AFQMC) method applied to dilute neutron matter at finite temperatures. We formulate the discrete Hubbard-Stratonovich transformation for the interaction with finite effective range which is free from the sign problem. The AFQMC results are compared with those obtained from exact diagonalization for a toy model. Preliminary calculations of energy and chemical potential as a function of temperature are presented. △ Less

Submitted 4 December, 2008; originally announced December 2008.

Comments: 7 pages, 1 figure, proceedings from XIV Nuclear Physics Workshop in Kazimierz Dolny, Poland

Journal ref: Int.J.Mod.Phys.E18:919-925,2009

Nuclear tetrahedral configurations at spin zero

Authors: Krzysztof Zberecki, Paul-Henri Heenen, Piotr Magierski

Abstract: The possibility of the existence of stable tetrahedral deformations at spin zero is investigated using the Skyrme-HFBCS approach and the generator coordinate method (GCM). The study is limited to nuclei in which the tetrahedral mode has been predicted to be favored on the basis of non self-consistent models. Our results indicate that a clear identification of tetrahedral deformations is unlikely… ▽ More The possibility of the existence of stable tetrahedral deformations at spin zero is investigated using the Skyrme-HFBCS approach and the generator coordinate method (GCM). The study is limited to nuclei in which the tetrahedral mode has been predicted to be favored on the basis of non self-consistent models. Our results indicate that a clear identification of tetrahedral deformations is unlikely as they are strongly mixed with the axial octupole mode. However, the excitation energies related to the tetrahedral mode are systematically lower than those of the axial octupole mode in all the nuclei included in this study. △ Less

Submitted 4 February, 2009; v1 submitted 5 June, 2008; originally announced June 2008.

Comments: 10 pages, 4 figures

Journal ref: Phys.Rev.C79:014319,2009

Quantum Monte Carlo Simulations of the BCS-BEC Crossover at Finite Temperature

Authors: Aurel Bulgac, Joaquin E. Drut, Piotr Magierski

Abstract: 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 issues. We report on results for the energy, entropy and chemical potential as a function of temperature. We give upper bounds on the critical temperatu… ▽ More 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 issues. 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. △ Less

Submitted 21 August, 2008; v1 submitted 21 March, 2008; originally announced March 2008.

Comments: 24 pages, 20 figures, published version, additional material and figures

Journal ref: Phys.Rev.A78:023625,2008

The Finite Temperature Pairing Gap of a Unitary Fermi Gas by Quantum Monte Carlo Calculations

Authors: P. Magierski, G. Wlazlowski, A. Bulgac, J. E. Drut

Abstract: 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,ω)$ using the methods of maximum entropy and singular value decomposition. From $A(p,ω)$ we determine the quasiparticle spectrum, which can be accurately parametrized by three functions of temperature: an effective mass $m^*$, a mean-f… ▽ More 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,ω)$ using the methods of maximum entropy and singular value decomposition. From $A(p,ω)$ 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 $Δ$. Below the critical temperature $T_c=0.15\varepsilon_F$ the results for $m^*$, $U$ and $Δ$ 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^*\approx 0.20\varepsilon_{F} > T_c$}. △ Less

Submitted 22 November, 2009; v1 submitted 9 January, 2008; originally announced January 2008.

Comments: 4 pages, 5 figures. Published version

Journal ref: Phys. Rev. Lett. 103, 210403 (2009)

Thermodynamics of a Trapped Unitary Fermi Gas

Authors: Aurel Bulgac, Joaquin E. Drut, Piotr Magierski

Abstract: 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$ i… ▽ More 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. △ Less

Submitted 18 September, 2007; v1 submitted 31 January, 2007; originally announced January 2007.

Comments: 5 pages, 4 figures, published version

Journal ref: Phys.Rev.Lett.99:120401,2007