-
Dipolar droplets of strongly interacting molecules
Authors:
Tim Langen,
Jordi Boronat,
Juan Sánchez-Baena,
Raúl Bombín,
Tijs Karman,
Ferran Mazzanti
Abstract:
We simulate a molecular Bose-Einstein condensate in the strongly dipolar regime, observing the existence of self-bound droplets, as well as their splitting into multiple droplets by confinement-induced frustration. Our quantum Monte Carlo approach goes beyond the limits of the established effective mean-field theories for dipolar quantum gases, revealing small droplets produced by strong dipolar i…
▽ More
We simulate a molecular Bose-Einstein condensate in the strongly dipolar regime, observing the existence of self-bound droplets, as well as their splitting into multiple droplets by confinement-induced frustration. Our quantum Monte Carlo approach goes beyond the limits of the established effective mean-field theories for dipolar quantum gases, revealing small droplets produced by strong dipolar interactions outside known stable regimes. The simulations include realistic molecular interactions and therefore have direct relevance for current and future experiments.
△ Less
Submitted 15 July, 2024; v1 submitted 12 July, 2024;
originally announced July 2024.
-
Mean field initialization of the Annealed Importance Sampling algorithm for an efficient evaluation of the Partition Function of Restricted Boltzmann Machines
Authors:
A. Prat Pou,
E. Romero,
J. Martí,
F. Mazzanti
Abstract:
Probabilistic models in physics often require from the evaluation of normalized Boltzmann factors, which in turn implies the computation of the partition function Z. Getting the exact value of Z, though, becomes a forbiddingly expensive task as the system size increases. This problem is also present in probabilistic learning models such as the Restricted Boltzmann Machine (RBM), where the situatio…
▽ More
Probabilistic models in physics often require from the evaluation of normalized Boltzmann factors, which in turn implies the computation of the partition function Z. Getting the exact value of Z, though, becomes a forbiddingly expensive task as the system size increases. This problem is also present in probabilistic learning models such as the Restricted Boltzmann Machine (RBM), where the situation is even worse as the exact learning rules implies the computation of Z at each iteration. A possible way to tackle this problem is to use the Annealed Importance Sampling (AIS) algorithm, which provides a tool to stochastically estimate the partition function of the system. So far, the standard application of the AIS algorithm starts from the uniform probability distribution and uses a large number of Monte Carlo steps to obtain reliable estimations of Z following an annealing process. In this work we show that both the quality of the estimation and the cost of the computation can be significantly improved by using a properly selected mean-field starting probability distribution. We perform a systematic analysis of AIS in both small- and large-sized problems, and compare the results to exact values in problems where these are known. As a result of our systematic analysis, we propose two successful strategies that work well in all the problems analyzed. We conclude that these are good starting points to estimate the partition function with AIS with a relatively low computational cost.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
A quantum Monte Carlo based density functional for Dysprosium dipolar system
Authors:
Raúl Bombín,
Viktor Cikojević,
Ferran Mazzanti,
Jordi Boronat
Abstract:
We present a quantum Monte Carlo based density functional to describe droplet formation and supersolidity in dipolar systems. The usual Lee-Huang-Yang term, accounting for quantum correlations in the conventional extended Gross-Pitaievskii equation (eGPE), has been substituted by the correlation energy evaluated with Quantum Monte Carlo. We demonstrate the ability of the new functional to reproduc…
▽ More
We present a quantum Monte Carlo based density functional to describe droplet formation and supersolidity in dipolar systems. The usual Lee-Huang-Yang term, accounting for quantum correlations in the conventional extended Gross-Pitaievskii equation (eGPE), has been substituted by the correlation energy evaluated with Quantum Monte Carlo. We demonstrate the ability of the new functional to reproduce existing experimental data for the minimum critical number of atoms $N_\mathrm{c}$ required for droplet formation. $N_\mathrm{c}$ is a challenging quantity for theoretical predictions, and the eGPE provides only a qualitative description of it, mainly when it is applied to Dysprosium. We also use the new approach to characterize the BEC-supersolid transition. The quantum Monte Carlo based functional can be easily implemented in any existing eGPE code, improving the description of dipolar systems without increasing the computational cost.
△ Less
Submitted 8 March, 2024; v1 submitted 8 November, 2023;
originally announced November 2023.
-
Off-diagonal long-range order in arrays of dipolar droplets
Authors:
R. Bombin,
F. Mazzanti,
J. Boronat
Abstract:
We report quantum Monte Carlo results of harmonically confined quantum Bose dipoles within a range of interactions covering the evolution from a gas phase to the formation of an array of droplets. Scaling the experimental setup to a computationally accessible domain we characterize that evolution in qualitative agreement with experiments. Our microscopic approach generates ground-state results fre…
▽ More
We report quantum Monte Carlo results of harmonically confined quantum Bose dipoles within a range of interactions covering the evolution from a gas phase to the formation of an array of droplets. Scaling the experimental setup to a computationally accessible domain we characterize that evolution in qualitative agreement with experiments. Our microscopic approach generates ground-state results free from approximations, albeit with some controlled statistical noise. The simultaneous estimation of the static structure factor and the one-body density matrix allows for a better knowledge of the quantum coherence between droplets. Our results show a narrow window of interaction strengths where diagonal and off-diagonal long-range order can coexist. This domain, which is the key signal of a supersolid state, is reduced with respect to the one predicted by the extended Gross-Pitaevskii equation. Differences are probably due to an increase of attraction in our model, observed previously in the calculation of critical atom numbers for single dipolar drops.
△ Less
Submitted 6 February, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Universal properties of dipolar Bose polarons in two dimensions
Authors:
Juan Sánchez-Baena,
Luis A. Peña Ardila,
Grigory Astrakharchik,
Ferran Mazzanti
Abstract:
We study the quasiparticle properties of a dipolar impurity immersed in a two-dimensional dipolar bath. We use the ab-initio Diffusion Monte Carlo technique to determine the polaron energy, effective mass and quasiparticle residue. We find that these quantities follow a universal behaviour when properly scaled in terms of the polarization angle and scattering length. This trend is maintained over…
▽ More
We study the quasiparticle properties of a dipolar impurity immersed in a two-dimensional dipolar bath. We use the ab-initio Diffusion Monte Carlo technique to determine the polaron energy, effective mass and quasiparticle residue. We find that these quantities follow a universal behaviour when properly scaled in terms of the polarization angle and scattering length. This trend is maintained over a wide range of values of the gas parameter, even in the highly correlated regime. Additionally, we show that the anisotropy of the impurity-bath interaction leads to an anisotropic effective mass that is unexpectedly larger in the direction of minimum repulsion of the impurity-bath interaction. Finally, we use our Monte Carlo results to check the accuracy of perturbative approaches and determine their range of validity in terms of the gas parameter.
△ Less
Submitted 7 June, 2023; v1 submitted 31 May, 2023;
originally announced May 2023.
-
Helium bubbles in liquid lithium: a potential issue for ITER
Authors:
Edgar Alvarez-Galera,
Jordi Marti,
Ferran Mazzanti,
Lluis Batet
Abstract:
Future fusion nuclear reactors will produce sustainable energy form the fusion of deuterium and tritium. In order to do so, the reactors will need to produce their own tritium through the neutron bombardment of lithium. Such reaction will produce tritium and helium inside the breeding blanket of the reactor. Helium can trigger nucleation mechanisms due to its very low solubility inside liquid meta…
▽ More
Future fusion nuclear reactors will produce sustainable energy form the fusion of deuterium and tritium. In order to do so, the reactors will need to produce their own tritium through the neutron bombardment of lithium. Such reaction will produce tritium and helium inside the breeding blanket of the reactor. Helium can trigger nucleation mechanisms due to its very low solubility inside liquid metals. Consequently, the knowledge and understanding of the microscopic processes of helium nucleation is crucial to improve the efficiency, sustainability and safety of the fusion energy production. The formation of helium bubbles inside the liquid metal used as breeding material may be a serious issue that has yet to be fully understood. We provide further insight on the behavior of lithium and helium mixtures at experimentally corresponding operating conditions (800~K and pressures between 1 and 100 bar) using a suitable microscopic model able to describe the helium and lithium atomic interactions, in excellent agreement with available experimental data. The simulations predict the formation of helium bubbles with radii around 10 Angstroem at ambient pressure and with surface tension values between 0.6-1.0 N/m, with a dependency of the concentration of helium. We also report cohesive energies of helium as well as a quantitative estimation of the Hildebrand and Kumar cohesion parameters.
△ Less
Submitted 1 April, 2023;
originally announced April 2023.
-
Striped Ultradilute Liquid of Dipolar Bosons in Two Dimensions
Authors:
Clemens Staudinger,
Diana Hufnagl,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
We investigate the phases of a Bose-Einstein condensate of dipolar atoms restricted to move in a two-dimensional plane. The dipole moments are all aligned in a direction tilted with respect to the plane normal. As a result of the attractive and repulsive components of the dipole-dipole interaction, the dipolar gas has a self-bound phase, which is stabilized by quantum fluctuations. Furthermore, ti…
▽ More
We investigate the phases of a Bose-Einstein condensate of dipolar atoms restricted to move in a two-dimensional plane. The dipole moments are all aligned in a direction tilted with respect to the plane normal. As a result of the attractive and repulsive components of the dipole-dipole interaction, the dipolar gas has a self-bound phase, which is stabilized by quantum fluctuations. Furthermore, tilting the dipoles tunes the anisotropy of the dipole-dipole interaction, which can trigger a spatial density modulation. In this work we study these two aspects and investigate the conditions for the formation of a self-bound and striped phase, which has been realized in experiments with dipolar droplets. We use a variational method based on the hypernetted-chain Euler-Lagrange optimization of a Jastrow-Feenberg ansatz for the many-body wave function to study the ground state properties. This method takes into account quantum fluctuations in a non-perturbative way and thus can be used also for strongly correlated systems.
△ Less
Submitted 20 March, 2023; v1 submitted 9 August, 2022;
originally announced August 2022.
-
Formal Modeling and Initial Analysis of the 4SECURail Case Study
Authors:
Franco Mazzanti,
Dimitri Belli
Abstract:
We present the case study developed in the context of the 4SECURail project and the approach used for its formal modeling and analysis. Starting from a simple SysML/UML behavioral model of the system requirements, three formal models have been developed using three different frameworks, namely UMC, ProB, and CADP/LNT. The paper shows how the different ways to represent and analyze the system from…
▽ More
We present the case study developed in the context of the 4SECURail project and the approach used for its formal modeling and analysis. Starting from a simple SysML/UML behavioral model of the system requirements, three formal models have been developed using three different frameworks, namely UMC, ProB, and CADP/LNT. The paper shows how the different ways to represent and analyze the system from the three different points of view allow us to take advantage of the resulting diversity.
△ Less
Submitted 18 March, 2022;
originally announced March 2022.
-
Time-dependent variational Monte Carlo study of the dynamic response of bosons in an optical lattice
Authors:
Mathias Gartner,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
We study the dynamics of a one-dimensional Bose gas at unit filling in both shallow and deep optical lattices and obtain the dynamic structure factor ${S(k,ω)}$ by monitoring the linear response to a weak probe pulse. We introduce a new procedure, based on the time-dependent variational Monte Carlo method (tVMC), which allows to evolve the system in real time, using as a variational model a Jastro…
▽ More
We study the dynamics of a one-dimensional Bose gas at unit filling in both shallow and deep optical lattices and obtain the dynamic structure factor ${S(k,ω)}$ by monitoring the linear response to a weak probe pulse. We introduce a new procedure, based on the time-dependent variational Monte Carlo method (tVMC), which allows to evolve the system in real time, using as a variational model a Jastrow-Feenberg wave function that includes pair correlations. Comparison with exact diagonalization results of ${S(k,ω)}$ obtained on a lattice in the Bose-Hubbard limit shows good agreement of the dispersion relation for sufficiently deep optical lattices, while for shallow lattices we observe the influence of higher Bloch bands. We also investigate non-linear response to strong pulses. From the power spectrum of the density fluctuations we obtain the excitation spectrum, albeit broadened, by higher harmonic generation after a strong pulse with a single low wave number. As a remarkable feature of our simulations we furthermore demonstrate that the full excitation spectrum can be retrieved from the power spectrum of the density fluctuations due to the stochastic noise inherent in any Monte Carlo method, without applying an actual perturbation.
△ Less
Submitted 27 May, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
Universal Properties of Anisotropic Dipolar Bosons in Two Dimensions
Authors:
J. Sánchez-Baena,
L. A. Peña Ardila,
G. Astrakharchik,
F. Mazzanti
Abstract:
The energy of ultra-dilute quantum many-body systems is known to exhibit a universal dependence on the gas parameter $x=n a_0^d$, with $n$ the density, $d$ the dimensionality of the space ($d=1,2,3$) and $a_0$ the $s$-wave scattering length. The universal regime typically extends up to $x\approx 0.001$, while at larger values specific details of the interaction start to be relevant and different m…
▽ More
The energy of ultra-dilute quantum many-body systems is known to exhibit a universal dependence on the gas parameter $x=n a_0^d$, with $n$ the density, $d$ the dimensionality of the space ($d=1,2,3$) and $a_0$ the $s$-wave scattering length. The universal regime typically extends up to $x\approx 0.001$, while at larger values specific details of the interaction start to be relevant and different model potentials lead to different results. Dipolar systems are peculiar in this regard since the anisotropy of the interaction makes $a_0$ depend on the polarization angle $α$, so different combinations of $n$ and $α$ can lead to the same value of the gas parameter $x$. In this work we analyze the scaling properties of dipolar bosons in two dimensions as a function of the density and polarization dependent scattering length up to very large values of the gas parameter $x$. Using Quantum Monte Carlo (QMC) methods we study the energy and the main structural and coherence properties of the ground state of a gas of dipolar bosons by varying the density and scattering length for fixed gas parameter. We find that the dipolar interaction shows relevant scaling laws up to unusually large values of $x$ that hold almost to the boundaries in the phase diagram where a transition to a stripe phase takes place.
△ Less
Submitted 13 May, 2022; v1 submitted 21 December, 2021;
originally announced December 2021.
-
Nucleation of helium in liquid lithium
Authors:
J. Martí,
F. Mazzanti,
G. E. Astrakharchik,
L. Batet,
L. Portos-Amill,
B. Pedreño
Abstract:
Fusion energy stands out as a promising alternative for a future decarbonised energy system. To be sustainable, future fusion nuclear reactors will have to produce their own tritium. In the so-called breeding blanket of a reactor, the neutron bombardment of lithium will produce the desired tritium, but also helium, which can trigger nucleation mechanisms owing to the very low solubility of helium…
▽ More
Fusion energy stands out as a promising alternative for a future decarbonised energy system. To be sustainable, future fusion nuclear reactors will have to produce their own tritium. In the so-called breeding blanket of a reactor, the neutron bombardment of lithium will produce the desired tritium, but also helium, which can trigger nucleation mechanisms owing to the very low solubility of helium in liquid metals. An understanding of the underlying microscopic processes is important for improving the efficiency, sustainability and reliability of the fusion energy conversion process. A spontaneous creation of helium drops or bubbles in the liquid metal used as breeding material in some designs may be a serious issue for the performance of the breeding blankets. This phenomenon has yet to be fully studied and understood. This work aims to provide some insight on the behavior of lithium and helium mixtures at experimentally corresponding operating conditions (843 K and pressures between 0.1 and 7 GPa). We report a microscopic study of the thermodynamic, structural and dynamical properties of lithium-helium mixtures, as a first step to the simulation of the environment in a nuclear fusion power plant. We introduce a microscopic model devised to describe the formation of helium drops in the thermodynamic range considered. A transition from a miscible homogeneous mixture to a phase-separated one, in which helium drops are nucleated, is observed as the pressure is increased above 0.175 GPa. The diffusion coefficient of lithium (2 Å 2 /ps) is in excellent agreement with reference experimental data, whereas the diffusion coefficient of helium is in the range of 1 Å 2 /ps and tends to decrease as pressure increases. The radii of helium drops have been found to be between 1 and 2 Å.
△ Less
Submitted 19 July, 2021;
originally announced July 2021.
-
Systematic Evaluation and Usability Analysis of Formal Tools for Railway System Design
Authors:
Alessio Ferrari,
Franco Mazzanti,
Davide Basile,
Maurice H. ter Beek
Abstract:
Formal methods and supporting tools have a long record of success in the development of safety-critical systems. However, no single tool has emerged as the dominant solution for system design. Each tool differs from the others in terms of the modeling language used, its verification capabilities and other complementary features, and each development context has peculiar needs that require differen…
▽ More
Formal methods and supporting tools have a long record of success in the development of safety-critical systems. However, no single tool has emerged as the dominant solution for system design. Each tool differs from the others in terms of the modeling language used, its verification capabilities and other complementary features, and each development context has peculiar needs that require different tools. This is particularly problematic for the railway industry, in which formal methods are highly recommended by the norms, but no actual guidance is provided for the selection of tools. To guide companies in the selection of the most appropriate formal tools to adopt in their contexts, a clear assessment of the features of the currently available tools is required.
To address this goal, this paper considers a set of 13 formal tools that have been used for railway system design, and it presents a systematic evaluation of such tools and a preliminary usability analysis of a subset of 7 tools, involving railway practitioners. The results are discussed considering the most desired aspects by industry and earlier related studies. While the focus is on the railway domain, the overall methodology can be applied to similar contexts. Our study thus contributes with a systematic evaluation of formal tools and it shows that despite the poor graphical interfaces, usability and maturity of the tools are not major problems, as claimed by contributions from the literature. Instead, support for process integration is the most relevant obstacle for the adoption of most of the tools.
△ Less
Submitted 13 May, 2021; v1 submitted 27 January, 2021;
originally announced January 2021.
-
Efficient Evaluation of the Partition Function of RBMs with Annealed Importance Sampling
Authors:
Ferran Mazzanti,
Enrique Romero
Abstract:
Probabilistic models based on Restricted Boltzmann Machines (RBMs) imply the evaluation of normalized Boltzmann factors, which in turn require from the evaluation of the partition function Z. The exact evaluation of Z, though, becomes a forbiddingly expensive task as the system size increases. This even worsens when one considers most usual learning algorithms for RBMs, where the exact evaluation…
▽ More
Probabilistic models based on Restricted Boltzmann Machines (RBMs) imply the evaluation of normalized Boltzmann factors, which in turn require from the evaluation of the partition function Z. The exact evaluation of Z, though, becomes a forbiddingly expensive task as the system size increases. This even worsens when one considers most usual learning algorithms for RBMs, where the exact evaluation of the gradient of the log-likelihood of the empirical distribution of the data includes the computation of Z at each iteration. The Annealed Importance Sampling (AIS) method provides a tool to stochastically estimate the partition function of the system. So far, the standard use of the AIS algorithm in the Machine Learning context has been done using a large number of Monte Carlo steps. In this work we show that this may not be required if a proper starting probability distribution is employed as the initialization of the AIS algorithm. We analyze the performance of AIS in both small- and large-sized problems, and show that in both cases a good estimation of Z can be obtained with little computational cost.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
Supersolid Striped Droplets in a Raman Spin-Orbit-Coupled System
Authors:
J. Sánchez-Baena,
J. Boronat,
F. Mazzanti
Abstract:
We analyze the role played by quantum fluctuations on a Raman Spin-Orbit Coupled system in the stripe phase. We show that beyond mean-field effects stabilize the collapse predicted by mean-field theory and induce the emergence of two phases: a gas and a liquid, which also show spatial periodicity along a privileged direction. We show that the energetically favored phase is determined by the Raman…
▽ More
We analyze the role played by quantum fluctuations on a Raman Spin-Orbit Coupled system in the stripe phase. We show that beyond mean-field effects stabilize the collapse predicted by mean-field theory and induce the emergence of two phases: a gas and a liquid, which also show spatial periodicity along a privileged direction. We show that the energetically favored phase is determined by the Raman coupling and the spin-dependent scattering lengths. We obtain the ground-state solution of the finite system by solving the extended Gross-Pitaevskii equation and find self-bound, droplet-like solutions that feature internal structure through a striped pattern. We estimate the critical number for binding associated to these droplets and show that their value is experimentally accessible. We report an approximate energy functional in order to ease the evaluation of the Lee-Huang-Yang correction in practical terms.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Reply to the Comment on "Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipolar Stripes"
Authors:
Raúl Bombín,
Ferran Mazzanti,
Jordi Boronat
Abstract:
This is a Reply to the Comment from F. Cinti and M. Boninsegni on our recent work on the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in a two-dimensional dipolar system [R.Bombín, F. Mazzanti and J. Boronat, Physical Review A 100, 063614 (2019)]. The main criticism about our work, expressed in that Comment, is that we did not explicitly report the two spatial contributions to the total…
▽ More
This is a Reply to the Comment from F. Cinti and M. Boninsegni on our recent work on the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in a two-dimensional dipolar system [R.Bombín, F. Mazzanti and J. Boronat, Physical Review A 100, 063614 (2019)]. The main criticism about our work, expressed in that Comment, is that we did not explicitly report the two spatial contributions to the total superfluid fraction. Here, we analyze our results for a point of the phase diagram corresponding to the stripe phase, close to the gas to stripe transition line, and for a temperature below the BKT critical temperature. The scaling with the system size of the contribution to the superfluid fraction, coming from the direction in which spatial order appears, shows that it remains finite in the thermodynamic limit, as we already stated in our original work. This allow us to state that the stripe phase is superfluid at low temperatures. Furthermore, we offer some comments that help to understand where the differences between the results of Cinti and Boninsegni and ours comes from.
△ Less
Submitted 20 October, 2020; v1 submitted 3 July, 2020;
originally announced July 2020.
-
Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipolar Stripes
Authors:
Raul Bombin,
Ferran Mazzanti,
Jordi Boronat
Abstract:
A two-dimensional quantum system of dipoles, with a polarization angle not perpendicular to the plane, shows a transition from a gas to a stripe phase. We have studied the thermal properties of these two phases using the path integral Monte Carlo (PIMC) method. By simulating the thermal density matrix, PIMC provides exact results for magnitudes of interest such as the superfluid fraction and the o…
▽ More
A two-dimensional quantum system of dipoles, with a polarization angle not perpendicular to the plane, shows a transition from a gas to a stripe phase. We have studied the thermal properties of these two phases using the path integral Monte Carlo (PIMC) method. By simulating the thermal density matrix, PIMC provides exact results for magnitudes of interest such as the superfluid fraction and the one-body density matrix. As it is well known, in two dimensions the superfluid-to-normal phase transition follows the Berezinskii-Kosterlitz-Thouless (BKT) scenario. Our results show that both the anisotropic gas and the stripe phases follow the BKT scaling laws. At fixed density and increasing the tilting angle, the transition temperature decreases in going from the gas to the stripe phase. Superfluidity in the perpendicular direction to the stripes is rather small close to the critical temperature but it becomes larger at lower temperatures, mainly close to the transition to the gas. Our results are in qualitative agreement with the supersolidity observed recently in a quasi-one-dimensional array of dipolar droplets.
△ Less
Submitted 9 December, 2019;
originally announced December 2019.
-
Two-dimensional repulsive Fermi polarons with short and long-range interactions
Authors:
Raúl Bombín,
Tommaso Comparin,
Gianluca Bertaina,
Ferran Mazzanti,
Stefano Giorgini,
Jordi Boronat
Abstract:
We study the repulsive polaron problem in a two-component two-dimensional system of fermionic atoms. We use two different interaction models: a short-range (hard-disk) potential and a dipolar potential. In our approach, all the atoms have the same mass and we consider the system to be composed of a uniform bath of a single species and a single atomic impurity. We use the diffusion Monte Carlo meth…
▽ More
We study the repulsive polaron problem in a two-component two-dimensional system of fermionic atoms. We use two different interaction models: a short-range (hard-disk) potential and a dipolar potential. In our approach, all the atoms have the same mass and we consider the system to be composed of a uniform bath of a single species and a single atomic impurity. We use the diffusion Monte Carlo method to evaluate polaron properties such as its chemical potential and pair distribution functions, together with a discussion on the deficit of volume induced by the impurity. We also evaluate observables that allow us to determine the validity of the quasi-particle picture: the quasi-particle residue and the effective mass of the polaron. Employing two different potentials allows us to identify the universality regime, where the properties depend only on the gas parameter $n a_s^2$ fixed by the bath density and the two-dimensional scattering length.
△ Less
Submitted 18 July, 2019; v1 submitted 24 May, 2019;
originally announced May 2019.
-
Dilute dipolar quantum droplets beyond the extended Gross-Pitaevskii equation
Authors:
Fabian Böttcher,
Matthias Wenzel,
Jan-Niklas Schmidt,
Mingyang Guo,
Tim Langen,
Igor Ferrier-Barbut,
Tilman Pfau,
Raúl Bombín,
Joan Sánchez-Baena,
Jordi Boronat,
Ferran Mazzanti
Abstract:
Dipolar quantum droplets are exotic quantum objects that are self-bound due to the subtle balance of attraction, repulsion and quantum correlations. Here we present a systematic study of the critical atom number of these self-bound droplets, comparing the experimental results with extended mean-field Gross-Pitaevskii equation (eGPE) and quantum Monte-Carlo simulations of the dilute system. The res…
▽ More
Dipolar quantum droplets are exotic quantum objects that are self-bound due to the subtle balance of attraction, repulsion and quantum correlations. Here we present a systematic study of the critical atom number of these self-bound droplets, comparing the experimental results with extended mean-field Gross-Pitaevskii equation (eGPE) and quantum Monte-Carlo simulations of the dilute system. The respective theoretical predictions differ, questioning the validity of the current theoretical state-of-the-art description of quantum droplets within the eGPE framework and indicating that correlations in the system are significant. Furthermore, we show that our system can serve as a sensitive testing ground for many-body theories in the near future.
△ Less
Submitted 4 November, 2019; v1 submitted 23 April, 2019;
originally announced April 2019.
-
Two-dimensional Mixture of Dipolar Fermions: Equation of State and Magnetic Phases
Authors:
Tommaso Comparin,
Raul Bombin,
Markus Holzmann,
Ferran Mazzanti,
Jordi Boronat,
Stefano Giorgini
Abstract:
We study a two-component mixture of fermionic dipoles in two dimensions at zero temperature, interacting via a purely repulsive $1/r^3$ potential. This model can be realized with ultracold atoms or molecules, when their dipole moments are aligned in the confinement direction orthogonal to the plane. We characterize the unpolarized mixture by means of the Diffusion Monte Carlo technique. Computing…
▽ More
We study a two-component mixture of fermionic dipoles in two dimensions at zero temperature, interacting via a purely repulsive $1/r^3$ potential. This model can be realized with ultracold atoms or molecules, when their dipole moments are aligned in the confinement direction orthogonal to the plane. We characterize the unpolarized mixture by means of the Diffusion Monte Carlo technique. Computing the equation of state, we identify the regime of validity for a mean-field theory based on a low-density expansion and compare our results with the hard-disk model of repulsive fermions. At high density, we address the possibility of itinerant ferromagnetism, namely whether the ground state can be fully polarized in the fluid phase. Within the fixed-node approximation, we show that the accuracy of Jastrow-Slater trial wave functions, even with the typical two-body backflow correction, is not sufficient to resolve the relevant energy differences. By making use of the iterative-backflow improved trial wave functions, we observe no signature of a fully-polarized ground state up to the freezing density.
△ Less
Submitted 9 April, 2019; v1 submitted 19 December, 2018;
originally announced December 2018.
-
Diffusion Monte Carlo methods for Spin-Orbit-Coupled ultracold Bose gases
Authors:
J. Sanchez-Baena,
J. Boronat,
F. Mazzanti
Abstract:
We present two Diffusion Monte Carlo (DMC) algorithms for systems of ultracold quantum gases featuring synthetic spin-orbit interactions. The first one is a discrete spin generalization of the T- moves spin-orbit DMC, which provides an upper bound to the fixed-phase energy. The second is a spin-integrated DMC method which recovers the fixed-phase property by avoiding the definition of the effectiv…
▽ More
We present two Diffusion Monte Carlo (DMC) algorithms for systems of ultracold quantum gases featuring synthetic spin-orbit interactions. The first one is a discrete spin generalization of the T- moves spin-orbit DMC, which provides an upper bound to the fixed-phase energy. The second is a spin-integrated DMC method which recovers the fixed-phase property by avoiding the definition of the effective Hamiltonian involved in the T-moves approach. The latter is a more accurate method but it is restricted to spin-independent two-body interactions. We report a comparison between both algorithms for different systems. As a check of the efficiency of both methods, we compare the DMC energies with results obtained with other numerical methods, finding agreement between both estimation
△ Less
Submitted 29 July, 2018; v1 submitted 21 July, 2018;
originally announced July 2018.
-
Self-bound Bose mixtures
Authors:
Clemens Staudinger,
Ferran Mazzanti,
Robert E. Zillich
Abstract:
Recent experiments confirmed that fluctuations beyond the mean-field approximation can lead to self-bound liquid droplets of ultra-dilute binary Bose mixtures. We proceed beyond the beyond-mean-field approximation, and study liquid Bose mixtures using the variational hypernetted-chain Euler Lagrange method, which accounts for correlations non-perturbatively. Focusing on the case of a mixture of un…
▽ More
Recent experiments confirmed that fluctuations beyond the mean-field approximation can lead to self-bound liquid droplets of ultra-dilute binary Bose mixtures. We proceed beyond the beyond-mean-field approximation, and study liquid Bose mixtures using the variational hypernetted-chain Euler Lagrange method, which accounts for correlations non-perturbatively. Focusing on the case of a mixture of uniform density, as realized inside large saturated droplets, we study the conditions for stability against evaporation of one of the components (both chemical potentials need to be negative) and against liquid-gas phase separation (spinodal instability), the latter being accompanied by a vanishing speed of sound. Dilute Bose mixtures are stable only in a narrow range near an optimal ratio $ρ_1/ρ_2$ and near the total energy minimum. Deviations from a universal dependence on the s-wave scattering lengths are significant despite the low density.
△ Less
Submitted 16 May, 2018;
originally announced May 2018.
-
Ten Diverse Formal Models for a CBTC Automatic Train Supervision System
Authors:
Franco Mazzanti,
Alessio Ferrari
Abstract:
Communications-based Train Control (CBTC) systems are metro signalling platforms, which coordinate and protect the movements of trains within the tracks of a station, and between different stations. In CBTC platforms, a prominent role is played by the Automatic Train Supervision (ATS) system, which automatically dispatches and routes trains within the metro network. Among the various functions, an…
▽ More
Communications-based Train Control (CBTC) systems are metro signalling platforms, which coordinate and protect the movements of trains within the tracks of a station, and between different stations. In CBTC platforms, a prominent role is played by the Automatic Train Supervision (ATS) system, which automatically dispatches and routes trains within the metro network. Among the various functions, an ATS needs to avoid deadlock situations, i.e., cases in which a group of trains block each other. In the context of a technology transfer study, we designed an algorithm for deadlock avoidance in train scheduling. In this paper, we present a case study in which the algorithm has been applied. The case study has been encoded using ten different formal verification environments, namely UMC, SPIN, NuSMV/nuXmv, mCRL2, CPN Tools, FDR4, CADP, TLA+, UPPAAL and ProB. Based on our experience, we observe commonalities and differences among the modelling languages considered, and we highlight the impact of the specific characteristics of each language on the presented models.
△ Less
Submitted 27 March, 2018;
originally announced March 2018.
-
Dipolar Bose Superstripes
Authors:
R. Bombin,
J. Boronat,
F. Mazzanti
Abstract:
We study the superfluid properties of a system of fully polarized dipolar bosons moving in the $xy$ plane. We focus on the general case where the polarization field forms an arbitrary angle $α$ with respect to the $z$ axis, while the system is still stable. We use the diffusion Monte Carlo and the path integral ground state methods to evaluate the one-body density matrix and the superfluid fractio…
▽ More
We study the superfluid properties of a system of fully polarized dipolar bosons moving in the $xy$ plane. We focus on the general case where the polarization field forms an arbitrary angle $α$ with respect to the $z$ axis, while the system is still stable. We use the diffusion Monte Carlo and the path integral ground state methods to evaluate the one-body density matrix and the superfluid fractions in the region of the phase diagram where the system forms stripes. Despite its oscillatory behavior, the presence of a finite large-distance asymptotic value in the $s$-wave component of the one-body density matrix indicates the existence of a Bose condensate. The superfluid fraction along the stripes direction is always close to 1, while in the $y$ direction decreases to a small value that is nevertheless different from zero. These two facts confirms that the stripe phase of the dipolar Bose gas in 2D is superfluid.
△ Less
Submitted 11 December, 2017; v1 submitted 25 August, 2017;
originally announced August 2017.
-
Optical lattices as a tool to study defect-induced superfluidity
Authors:
Grigory E. Astrakharchik,
Konstantin V. Krutitsky,
Maciej Lewenstein,
Ferran Mazzanti,
Jordi Boronat
Abstract:
We study the superfluid response, the energetic and structural properties of a one-dimensional ultracold Bose gas in an optical lattice of arbitrary strength. We use the Bose-Fermi mapping in the limit of infinitely large repulsive interaction and the diffusion Monte Carlo method in the case of finite interaction. For slightly incommensurate fillings we find a superfluid behavior which is discusse…
▽ More
We study the superfluid response, the energetic and structural properties of a one-dimensional ultracold Bose gas in an optical lattice of arbitrary strength. We use the Bose-Fermi mapping in the limit of infinitely large repulsive interaction and the diffusion Monte Carlo method in the case of finite interaction. For slightly incommensurate fillings we find a superfluid behavior which is discussed in terms of vacancies and interstitials. It is shown that both the excitation spectrum and static structure factor are different for the cases of microscopic and macroscopic fractions of defects. This system provides a extremely well-controlled model for studying defect-induced superfluidity.
△ Less
Submitted 6 September, 2017; v1 submitted 22 December, 2016;
originally announced December 2016.
-
Gapped spectrum in pair-superfluid bosons
Authors:
G. E. Astrakharchik,
R. E. Zillich,
F. Mazzanti,
J. Boronat
Abstract:
We study the ground state of a bilayer system of dipolar bosons with dipoles oriented by an external field perpendicularly to the two parallel planes. By decreasing the interlayer distance, for a fixed value of the strength of the dipolar interaction, the system undergoes a quantum phase transition from an atomic to a pair superfluid. We investigate the excitation spectrum across this transition b…
▽ More
We study the ground state of a bilayer system of dipolar bosons with dipoles oriented by an external field perpendicularly to the two parallel planes. By decreasing the interlayer distance, for a fixed value of the strength of the dipolar interaction, the system undergoes a quantum phase transition from an atomic to a pair superfluid. We investigate the excitation spectrum across this transition by using microscopic approaches. Quantum Monte Carlo methods are employed to obtain the static structure factors and intermediate scattering functions in imaginary time. The dynamic response is calculated using both the correlated basis functions (CBF) method and the approximate inversion of the Laplace transform of the quantum Monte Carlo imaginary time data. In the atomic phase, both density and spin excitations are gapless. However, in the pair-superfluid phase a gap opens in the excitation energy of the spin mode. For small separation between layers, the minimal spin excitation energy equals the binding energy of a dimer and is twice the gap value.
△ Less
Submitted 25 August, 2016;
originally announced August 2016.
-
Droplets of trapped quantum dipolar bosons
Authors:
A. Macia,
J. Sánchez-Baena,
J. Boronat,
F. Mazzanti
Abstract:
Strongly interacting systems of dipolar bosons in three dimensions confined by harmonic traps are analyzed using the exact Path Integral Ground State Monte Carlo method. By adding a repulsive two-body potential, we find a narrow window of interaction parameters leading to stable ground- state configurations of droplets in a crystalline arrangement. We find that this effect is entirely due to the i…
▽ More
Strongly interacting systems of dipolar bosons in three dimensions confined by harmonic traps are analyzed using the exact Path Integral Ground State Monte Carlo method. By adding a repulsive two-body potential, we find a narrow window of interaction parameters leading to stable ground- state configurations of droplets in a crystalline arrangement. We find that this effect is entirely due to the interaction present in the Hamiltonian without resorting to additional stabilizing mechanisms or specific three-body forces. We analyze the number of droplets formed in terms of the Hamiltonian parameters, relate them to the corresponding s-wave scattering length, and discuss a simple scaling model for the density profiles. Our results are in qualitative agreement with recent experiments showing a quantum Rosensweig instability in trapped Dy atoms.
△ Less
Submitted 25 July, 2016;
originally announced July 2016.
-
One-dimensional Bose gas in optical lattices of arbitrary strength
Authors:
Grigory E. Astrakharchik,
Konstantin V. Krutitsky,
Maciej Lewenstein,
Ferran Mazzanti
Abstract:
One-dimensional Bose gas with contact interaction in optical lattices at zero temperature is investigated by means of the exact diffusion Monte Carlo algorithm. The results obtained from the fundamental continuous model are compared with those obtained from the lattice (discrete) Bose-Hubbard model, using exact diagonalization, and from the quantum sine-Gordon model. We map out the complete phase…
▽ More
One-dimensional Bose gas with contact interaction in optical lattices at zero temperature is investigated by means of the exact diffusion Monte Carlo algorithm. The results obtained from the fundamental continuous model are compared with those obtained from the lattice (discrete) Bose-Hubbard model, using exact diagonalization, and from the quantum sine-Gordon model. We map out the complete phase diagram of the continuous model and determine the regions of applicability of the Bose-Hubbard model. Various physical quantities characterizing the systems are calculated and it is demonstrated that the sine-Gordon model used for shallow lattices is inaccurate.
△ Less
Submitted 8 February, 2016; v1 submitted 4 September, 2015;
originally announced September 2015.
-
A Neighbourhood-Based Stopping Criterion for Contrastive Divergence Learning
Authors:
E. Romero,
F. Mazzanti,
J. Delgado
Abstract:
Restricted Boltzmann Machines (RBMs) are general unsupervised learning devices to ascertain generative models of data distributions. RBMs are often trained using the Contrastive Divergence learning algorithm (CD), an approximation to the gradient of the data log-likelihood. A simple reconstruction error is often used as a stopping criterion for CD, although several authors \cite{schulz-et-al-Conve…
▽ More
Restricted Boltzmann Machines (RBMs) are general unsupervised learning devices to ascertain generative models of data distributions. RBMs are often trained using the Contrastive Divergence learning algorithm (CD), an approximation to the gradient of the data log-likelihood. A simple reconstruction error is often used as a stopping criterion for CD, although several authors \cite{schulz-et-al-Convergence-Contrastive-Divergence-2010-NIPSw, fischer-igel-Divergence-Contrastive-Divergence-2010-ICANN} have raised doubts concerning the feasibility of this procedure. In many cases the evolution curve of the reconstruction error is monotonic while the log-likelihood is not, thus indicating that the former is not a good estimator of the optimal stopping point for learning. However, not many alternatives to the reconstruction error have been discussed in the literature. In this manuscript we investigate simple alternatives to the reconstruction error, based on the inclusion of information contained in neighboring states to the training set, as a stopping criterion for CD learning.
△ Less
Submitted 24 July, 2015;
originally announced July 2015.
-
Quantum Monte Carlo estimation of complex-time correlations for the study of the ground-state dynamic structure function
Authors:
Riccardo Rota,
Joaquim Casulleras,
Ferran Mazzanti,
Jordi Boronat
Abstract:
We present a method based on the Path Integral Monte Carlo formalism for the calculation of ground-state time correlation functions in quantum systems. The key point of the method is the consideration of time as a complex variable whose phase $δ$ acts as an adjustable parameter. By using high-order approximations for the quantum propagator, it is possible to obtain Monte Carlo data all the way fro…
▽ More
We present a method based on the Path Integral Monte Carlo formalism for the calculation of ground-state time correlation functions in quantum systems. The key point of the method is the consideration of time as a complex variable whose phase $δ$ acts as an adjustable parameter. By using high-order approximations for the quantum propagator, it is possible to obtain Monte Carlo data all the way from purely imaginary time to $δ$ values near the limit of real time. As a consequence, it is possible to infer accurately the spectral functions using simple inversion algorithms. We test this approach in the calculation of the dynamic structure function $S(q,ω)$ of two one-dimensional model systems, harmonic and quartic oscillators, for which $S(q,ω)$ can be exactly calculated. We notice a clear improvement in the calculation of the dynamic response with respect to the common approach based on the inverse Laplace transform of the imaginary-time correlation function.
△ Less
Submitted 4 March, 2015;
originally announced March 2015.
-
Phase diagram of dipolar bosons in 2D with tilted polarization
Authors:
A. Macia,
J. Boronat,
F. Mazzanti
Abstract:
We analyze the ground state of a system of dipolar bosons moving in the $XY$ plane and such that their dipolar moments are all aligned in a fixed direction in space. We focus on the general case where the polarization field forms a generic angle $α$ with respect to the $Z$ axis. We use the Path Integral Ground State method to analyze the static properties of the system as both $α$ and the density…
▽ More
We analyze the ground state of a system of dipolar bosons moving in the $XY$ plane and such that their dipolar moments are all aligned in a fixed direction in space. We focus on the general case where the polarization field forms a generic angle $α$ with respect to the $Z$ axis. We use the Path Integral Ground State method to analyze the static properties of the system as both $α$ and the density $n$ vary over a wide range were the system is stable. We use the maximum of the static structure function as an order parameter to characterize the different phases and the transition lines among them. We find that aside of a superfluid gas and a solid phase, the system reaches a stripe phase at large tilting angles that is entirely induced by the anisotropic character of the interaction. We also show that the quantum phase transition from the gas to the stripe phase is of second order, and report approximate values for the critical exponents.
△ Less
Submitted 25 July, 2014;
originally announced July 2014.
-
Single-particle vs. pair superfluidity in a bilayer system of dipolar bosons
Authors:
A. Macia,
G. E. Astrakharchik,
F. Mazzanti,
S. Giorgini,
J. Boronat
Abstract:
We consider the ground state of a bilayer system of dipolar bosons, where dipoles are oriented by an external field in the direction perpendicular to the parallel planes. Quantum Monte Carlo methods are used to calculate the ground-state energy, the one-body and two-body density matrix, and the superfluid response as a function of the separation between layers. We find that by decreasing the inter…
▽ More
We consider the ground state of a bilayer system of dipolar bosons, where dipoles are oriented by an external field in the direction perpendicular to the parallel planes. Quantum Monte Carlo methods are used to calculate the ground-state energy, the one-body and two-body density matrix, and the superfluid response as a function of the separation between layers. We find that by decreasing the interlayer distance for fixed value of the strength of the dipolar interaction, the system undergoes a quantum phase transition from a single-particle to a pair superfluid. The single-particle superfluid is characterized by a finite value of both the atomic condensate and the super-counterfluid density. The pair superfluid phase is found to be stable against formation of many-body cluster states and features a gap in the spectrum of elementary excitations.
△ Less
Submitted 25 July, 2014;
originally announced July 2014.
-
Stopping Criteria in Contrastive Divergence: Alternatives to the Reconstruction Error
Authors:
David Buchaca,
Enrique Romero,
Ferran Mazzanti,
Jordi Delgado
Abstract:
Restricted Boltzmann Machines (RBMs) are general unsupervised learning devices to ascertain generative models of data distributions. RBMs are often trained using the Contrastive Divergence learning algorithm (CD), an approximation to the gradient of the data log-likelihood. A simple reconstruction error is often used to decide whether the approximation provided by the CD algorithm is good enough,…
▽ More
Restricted Boltzmann Machines (RBMs) are general unsupervised learning devices to ascertain generative models of data distributions. RBMs are often trained using the Contrastive Divergence learning algorithm (CD), an approximation to the gradient of the data log-likelihood. A simple reconstruction error is often used to decide whether the approximation provided by the CD algorithm is good enough, though several authors (Schulz et al., 2010; Fischer & Igel, 2010) have raised doubts concerning the feasibility of this procedure. However, not many alternatives to the reconstruction error have been used in the literature. In this manuscript we investigate simple alternatives to the reconstruction error in order to detect as soon as possible the decrease in the log-likelihood during learning.
△ Less
Submitted 9 April, 2014; v1 submitted 20 December, 2013;
originally announced December 2013.
-
Ground state properties and excitation spectrum of a two dimensional gas of bosonic dipoles
Authors:
A. Macia,
F. Mazzanti,
J. Boronat
Abstract:
We present a quantum Monte Carlo study of two-dimensional dipolar Bose gases in the limit of zero temperature. The analysis is mainly focused on the anisotropy effects induced in the homogeneous gas when the polarization angle with respect to the plane is changed. We restrict our study to the regime where the dipolar interaction is strictly repulsive, although the strength of the pair repulsion de…
▽ More
We present a quantum Monte Carlo study of two-dimensional dipolar Bose gases in the limit of zero temperature. The analysis is mainly focused on the anisotropy effects induced in the homogeneous gas when the polarization angle with respect to the plane is changed. We restrict our study to the regime where the dipolar interaction is strictly repulsive, although the strength of the pair repulsion depends on the vector interparticle distance. Our results show that the effect of the anisotropy in the energy per particle scales with the gas parameter at low densities as expected, and that this scaling is preserved for all polarization angles even at the largest densities considered here. We also evaluate the excitation spectrum of the dipolar Bose gas in the context of the Feynman approximation and compare the results obtained with the Bogoliubov ones. As expected, we find that these two approximations agree at very low densities, while they start to deviate from each other as the density increases. For the largest densities studied, we observe a significant influence of the anisotropy of the dipole-dipole interaction in the excitation spectrum.
△ Less
Submitted 15 October, 2012;
originally announced October 2012.
-
Excitations and Stripe Phase Formation in a 2D Dipolar Bose Gas with Tilted Polarization
Authors:
Adrian Macia,
Diana Hufnagl,
Ferran Mazzanti,
Jordi Boronat,
Robert E. Zillich
Abstract:
We present calculations of the ground state and excitations of an anisotropic dipolar Bose gas in two dimensions, realized by a non-perpendicular polarization with respect to the system plane. For sufficiently high density an increase of the polarization angle leads to a density instability of the gas phase in the direction where the anisotropic interaction is strongest. Using a dynamic many-body…
▽ More
We present calculations of the ground state and excitations of an anisotropic dipolar Bose gas in two dimensions, realized by a non-perpendicular polarization with respect to the system plane. For sufficiently high density an increase of the polarization angle leads to a density instability of the gas phase in the direction where the anisotropic interaction is strongest. Using a dynamic many-body theory, we calculate the dynamic structure function in the gas phase which shows the anisotropic dispersion of the excitations. We find that the energy of roton excitations in the strongly interacting direction decreases with increasing polarization angle and almost vanishes close to the instability. Exact path integral ground state Monte Carlo simulations show that this instability is indeed a quantum phase transition to a stripe phase, characterized by long-range order in the strongly interacting direction.
△ Less
Submitted 7 September, 2012;
originally announced September 2012.
-
Ferromagnetic transition of a two-component Fermi gas of Hard Spheres
Authors:
F. Arias de Saavedra,
F. Mazzanti,
J. Boronat,
A. Polls
Abstract:
We use microscopic many-body theory to analyze the problem of itinerant ferromagnetism in a repulsive atomic Fermi gas of Hard Spheres. Using simple arguments, we show that the available theoretical predictions for the onset of the ferromagnetic transition predict a transition point at a density ($k_F a \sim 1$) that is too large to be compatible with the universal low-density expansion of the ene…
▽ More
We use microscopic many-body theory to analyze the problem of itinerant ferromagnetism in a repulsive atomic Fermi gas of Hard Spheres. Using simple arguments, we show that the available theoretical predictions for the onset of the ferromagnetic transition predict a transition point at a density ($k_F a \sim 1$) that is too large to be compatible with the universal low-density expansion of the energy. We present new variational calculations for the hard-sphere Fermi gas, in the framework of Fermi hypperneted chain theory, that shift the transition to higher densities ($k_F a \sim 1.8$). Backflow correlations, which are mainly active in the unpolarized system, are essential for this shift.
△ Less
Submitted 12 March, 2012;
originally announced March 2012.
-
Zero-temperature phase diagram of Yukawa bosons
Authors:
O. N. Osychenko,
G. E. Astrakharchik,
F. Mazzanti,
J. Boronat
Abstract:
We study the zero-temperature phase diagram of bosons interacting via screened Coulomb (Yukawa) potential by means of the diffusion Monte Carlo method. The Yukawa potential is used as a model interaction in the neutron matter, dusty plasmas and charged colloids. As shown by D. S. Petrov et al. [Phys. Rev. Lett. 99, 130407 (2007)], interactions between weakly bound molecules of heavy and light ferm…
▽ More
We study the zero-temperature phase diagram of bosons interacting via screened Coulomb (Yukawa) potential by means of the diffusion Monte Carlo method. The Yukawa potential is used as a model interaction in the neutron matter, dusty plasmas and charged colloids. As shown by D. S. Petrov et al. [Phys. Rev. Lett. 99, 130407 (2007)], interactions between weakly bound molecules of heavy and light fermionic atoms are described by an effective Yukawa potential with a strength related to the heavy-light mass ratio M/m which might lead to crystallization in a two-dimensional geometry if the mass ratio of heavy-light fermions exceeds a certain critical value. In the present work we do a thorough study of the quantum three-dimensional Yukawa system. For strong interactions (equivalently, large mass ratios) the system experiences several phase transitions as the density is increased, passing from gas to solid and to gas phase again. Weakly interacting Yukawa particles do not crystallize at any density. We find the minimal interaction strength at which the crystallization happens. In terms of the two-component fermionic system, this strength corresponds to a heavy-light mass ratio of M/m ~ 180, so that it is impossible to realize the gas-crystal transition in a conventional bulk system. For the Yukawa model of fermionic mixtures we also analyze the possibility of building molecular systems with very large effective mass ratios by confining the heavy component to a sufficiently deep optical lattice. We show how the effective mass of the heavy component can be made arbitrarily large by increasing the lattice depth, thus leading to a tunable effective mass ratio that can be used to realize a molecular superlattice.
△ Less
Submitted 4 June, 2012; v1 submitted 29 December, 2011;
originally announced December 2011.
-
Microscopic description of anisotropic low-density dipolar Bose gases in two dimensions
Authors:
A. Macia,
F. Mazzanti,
J. Boronat,
R. E. Zillich
Abstract:
A microscopic description of the zero energy two-body ground state and many-body static properties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is presented and discussed. By changing the polarization angle with respect to the plane, we study the impact of the anisotropy, present in the dipole--dipole interaction, on the energy per particle, comparing the…
▽ More
A microscopic description of the zero energy two-body ground state and many-body static properties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is presented and discussed. By changing the polarization angle with respect to the plane, we study the impact of the anisotropy, present in the dipole--dipole interaction, on the energy per particle, comparing the results with mean field predictions. We restrict the analysis to the regime where the interaction is always repulsive, although the strength of the repulsion depends on the orientation with respect to the polarization field. We present a series expansion of the solution of the zero energy two-body problem which allows us to find the scattering length of the interaction and to build a suitable Jastrow factor that we use as a trial wave function for both a variational and diffusion Monte Carlo simulation of the infinite system. We find that the anisotropy has an almost negligible impact on the ground state properties of the many-body system in the universal regime where the scattering length governs the physics of the system. We also show that scaling in the gas parameter persists in the dipolar case up to values where other isotropic interactions with the same scattering length yield different predictions.
△ Less
Submitted 30 May, 2011;
originally announced May 2011.
-
Elementary Excitations in Bose-Einstein Condensates at Large Scattering Lengths
Authors:
R. Sarjonen,
M. Saarela,
F. Mazzanti
Abstract:
We present a theoretical analysis of excitation modes in Bose-Einstein condensates in ultracold alkali-metal gases for large scattering lengths and momenta where corrections to the mean field approximation become important. We assume that the effective interaction in the metastable, single channel, gaseous phase has a well defined Fourier transform that scales with the scattering length. Based o…
▽ More
We present a theoretical analysis of excitation modes in Bose-Einstein condensates in ultracold alkali-metal gases for large scattering lengths and momenta where corrections to the mean field approximation become important. We assume that the effective interaction in the metastable, single channel, gaseous phase has a well defined Fourier transform that scales with the scattering length. Based on this we show that for increasing scattering lengths or equivalently increasing densities the system becomes less correlated and that at large values of the scattering length Bragg scattering measures directly the Fourier transform of the effective two-body potential. We construct model potentials which fit the recently measured line shifts in $^{85}$Rb by Papp et al. (Phys. Rev. Lett. {\bf 101}, 135301 (2008)), and show that they fix the low momentum expansion of the effective range function. We find excellent agreement with the experimental data when the effective range is $\ll 1$ and the coefficient of the $k^4$-term is $-7.5 \pm 0.5$ in scattering length units. The resolution in Bragg scattering experiments so far does not reveal details of the frequency dependence in the dynamic structure function $S(k,ω)$ and we show that the Feynman spectrum determines the measured line shifts. We propose the possibility of a transition to a novel density wave state.
△ Less
Submitted 17 December, 2009;
originally announced December 2009.
-
Low-dimensional weakly interacting Bose gases: non-universal equations of state
Authors:
G. E. Astrakharchik,
J. Boronat,
I. L. Kurbakov,
Yu. E. Lozovik,
F. Mazzanti
Abstract:
The zero-temperature equation of state is analyzed in low-dimensional bosonic systems. In the dilute regime the equation of state is universal in terms of the gas parameter, i.e. it is the same for different potentials with the same value of the s-wave scattering length. Series expansions of the universal equation of state are reported for one- and two- dimensional systems. We propose to use the…
▽ More
The zero-temperature equation of state is analyzed in low-dimensional bosonic systems. In the dilute regime the equation of state is universal in terms of the gas parameter, i.e. it is the same for different potentials with the same value of the s-wave scattering length. Series expansions of the universal equation of state are reported for one- and two- dimensional systems. We propose to use the concept of energy-dependent s-wave scattering length for obtaining estimations of non-universal terms in the energy expansion. We test this approach by making a comparison to exactly solvable one-dimensional problems and find that the generated terms have the correct structure. The applicability to two-dimensional systems is analyzed by comparing with results of Monte Carlo simulations. The prediction for the non-universal behavior is qualitatively correct and the densities, at which the deviations from the universal equation of state become visible, are estimated properly. Finally, the possibility of observing the non-universal terms in experiments with trapped gases is also discussed.
△ Less
Submitted 30 July, 2009;
originally announced July 2009.
-
High-order Time Expansion Path Integral Ground State
Authors:
R. Rota,
J. Casulleras,
F. Mazzanti,
J. Boronat
Abstract:
The feasibility of path integral Monte Carlo ground state calculations with very few beads using a high-order short-time Green's function expansion is discussed. An explicit expression of the evolution operator which provides dramatic enhancements in the quality of ground-state wave-functions is examined. The efficiency of the method makes possible to remove the trial wave function and thus obta…
▽ More
The feasibility of path integral Monte Carlo ground state calculations with very few beads using a high-order short-time Green's function expansion is discussed. An explicit expression of the evolution operator which provides dramatic enhancements in the quality of ground-state wave-functions is examined. The efficiency of the method makes possible to remove the trial wave function and thus obtain completely model-independent results still with a very small number of beads. If a single iteration of the method is used to improve a given model wave function, the result is invariably a shadow-type wave function, whose precise content is provided by the high-order algorithm employed.
△ Less
Submitted 23 July, 2009;
originally announced July 2009.
-
Off-diagonal Ground State Properties of a 1D Gas of Fermi Hard Rods
Authors:
F. Mazzanti,
G. E. Astrakharchik,
J. Boronat,
J. Casulleras
Abstract:
A variational Monte Carlo calculation of the one-body density matrix and momentum distribution of a system of Fermi hard rods (HR) is presented and compared with the same quantities for its bosonic counterpart. The calculation is exact within statistical errors since we sample the exact ground state wave function, whose analytical expression is known. The numerical results are in good agreement…
▽ More
A variational Monte Carlo calculation of the one-body density matrix and momentum distribution of a system of Fermi hard rods (HR) is presented and compared with the same quantities for its bosonic counterpart. The calculation is exact within statistical errors since we sample the exact ground state wave function, whose analytical expression is known. The numerical results are in good agreement with known asymptotic expansions valid for Luttinger liquids. We find that the difference between the absolute value of the bosonic and fermionic density matrices becomes marginally small as the density increases. In this same regime, the corresponding momentum distributions merge into a common profile that is independent of the statistics. Non-analytical contributions to the one--body density matrix are also discussed and found to be less relevant with increasing density.
△ Less
Submitted 21 December, 2007;
originally announced December 2007.
-
Ground-State Properties of a One-Dimensional System of Hard Rods
Authors:
F. Mazzanti,
G. E. Astrakharchik,
J. Boronat,
J. Casulleras
Abstract:
A quantum Monte Carlo simulation of a system of hard rods in one dimension is presented and discussed. The calculation is exact since the analytical form of the wavefunction is known, and is in excellent agreement with predictions obtained from asymptotic expansions valid at large distances. The analysis of the static structure factor and the pair distribution function indicates that a solid-lik…
▽ More
A quantum Monte Carlo simulation of a system of hard rods in one dimension is presented and discussed. The calculation is exact since the analytical form of the wavefunction is known, and is in excellent agreement with predictions obtained from asymptotic expansions valid at large distances. The analysis of the static structure factor and the pair distribution function indicates that a solid-like and a gas-like phases exist at high and low densities, respectively. The one-body density matrix decays following a power-law at large distances and produces a divergence in the low density momentum distribution at k=0 which can be identified as a quasi-condensate.
△ Less
Submitted 30 May, 2007;
originally announced May 2007.
-
Ground state properties of a dilute homogeneous Bose gas of hard disks in two dimensions
Authors:
F. Mazzanti,
A. Polls,
A. Fabrocini
Abstract:
The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a variational many-body approach based on a Jastrow correlated ground state wave function. The asymptotic behaviors of the radial distribution function and the one-body density matrix are analyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain energy functional. Our…
▽ More
The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a variational many-body approach based on a Jastrow correlated ground state wave function. The asymptotic behaviors of the radial distribution function and the one-body density matrix are analyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain energy functional. Our results show important deviations from those of the available low density expansions, already at gas parameter values $x\sim 0.001$. The condensate fraction in 2D is also computed and found generally lower than the 3D one at the same $x$.
△ Less
Submitted 27 October, 2004;
originally announced October 2004.
-
Energy and Structure of Hard-Sphere Bose Gases in three and two dimensions
Authors:
F. Mazzanti,
A. Polls,
A. Fabrocini
Abstract:
The energy and structure of dilute gases of hard spheres in three dimensions is discussed, together with some aspects of the corresponding 2D systems. A variational approach in the framework of the Hypernetted Chain Equations (HNC) is used starting from a Jastrow wavefunction that is optimized to produce the best two--body correlation factor with the appropriate long range. Relevant quantities d…
▽ More
The energy and structure of dilute gases of hard spheres in three dimensions is discussed, together with some aspects of the corresponding 2D systems. A variational approach in the framework of the Hypernetted Chain Equations (HNC) is used starting from a Jastrow wavefunction that is optimized to produce the best two--body correlation factor with the appropriate long range. Relevant quantities describing static properties of the system are studied as a function of the gas parameter $x=ρa^d$ where $ρ$, $a$ and $d$ are the density, $s$--wave scattering length of the potential and dimensionality of the space, respectively. The occurrence of a maximum in the radial distribution function and in the momentum distribution is a natural effect of the correlations when $x$ increases. Some aspects of the asymptotic behavior of the functions characterizing the structure of the systems are also investigated.
△ Less
Submitted 26 October, 2004;
originally announced October 2004.
-
High momentum response of liquid 3He
Authors:
F. Mazzanti,
A. Polls,
J. Boronat,
J. Casulleras
Abstract:
A final-state-effects formalism suitable to analyze the high-momentum response of Fermi liquids is presented and used to study the dynamic structure function of liquid $^3$He. The theory, developed as a natural extension of the Gersch-Rodriguez formalism, incorporates the Fermi statistics explicitely through a new additive term which depends on the semi-diagonal two-body density matrix. The use…
▽ More
A final-state-effects formalism suitable to analyze the high-momentum response of Fermi liquids is presented and used to study the dynamic structure function of liquid $^3$He. The theory, developed as a natural extension of the Gersch-Rodriguez formalism, incorporates the Fermi statistics explicitely through a new additive term which depends on the semi-diagonal two-body density matrix. The use of a realistic momentum distribution, calculated using the diffusion Monte Carlo method, and the inclusion of this additive correction allows for a good agreement with available deep-inelastic neutron scattering data.
△ Less
Submitted 9 January, 2004;
originally announced January 2004.
-
Energy and structure of dilute hard- and soft-sphere gases
Authors:
F. Mazzanti,
A. Polls,
A. Fabrocini
Abstract:
The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, as the variational correlated theory, the Bogoliubov model and the uniform limit approximation, valid in the weak interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. A Jastrow type correlation provides…
▽ More
The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, as the variational correlated theory, the Bogoliubov model and the uniform limit approximation, valid in the weak interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. A Jastrow type correlation provides a good description of the systems, both hard- and soft-spheres, if the hypernetted chain energy functional is freely minimized and the resulting Euler equation is solved. The study of the soft-spheres potentials confirms the appearance of a dependence of the energy on the shape of the potential at gas paremeter values of $x \sim 0.001$. For quantities other than the energy, such as the radial distribution functions and the momentum distributions, the dependence appears at any value of $x$. The occurrence of a maximum in the radial distribution function, in the momentum distribution and in the excitation spectrum is a natural effect of the correlations when $x$ increases. The asymptotic behaviors of the functions characterizing the structure of the systems are also investigated. The uniform limit approach results very easy to implement and provides a good description of the soft-sphere gas. Its reliability improves when the interaction weakens.
△ Less
Submitted 21 May, 2003;
originally announced May 2003.
-
Finite Temperatur Dynamic Susceptibility of the Free Bose Gas
Authors:
F. Mazzanti,
A. Polls
Abstract:
A detailed calculation of the real part of the finite temperature dynamic susceptibility of the free Bose gas is presented. After a short discussion on the different ways in which it can be calculated for temperatures above and below the Bose--Einstein transition temperature, its main properties and its evolution with $q$ and $T$ are analyzed. Finally, expressions for the lowest order energy wei…
▽ More
A detailed calculation of the real part of the finite temperature dynamic susceptibility of the free Bose gas is presented. After a short discussion on the different ways in which it can be calculated for temperatures above and below the Bose--Einstein transition temperature, its main properties and its evolution with $q$ and $T$ are analyzed. Finally, expressions for the lowest order energy weighted sum rules are also derived and studied.
△ Less
Submitted 20 March, 2001;
originally announced March 2001.
-
High-momentum dynamic structure function of liquid 3He-4He mixtures: a microscopic approach
Authors:
F. Mazzanti,
A. Polls,
J. Boronat
Abstract:
The high-momentum dynamic structure function of liquid 3He-4He mixtures has been studied introducing final state effects. Corrections to the impulse approximation have been included using a generalized Gersch-Rodriguez theory that properly takes into account the Fermi statistics of 3He atoms. The microscopic inputs, as the momentum distributions and the two-body density matrices, correspond to a…
▽ More
The high-momentum dynamic structure function of liquid 3He-4He mixtures has been studied introducing final state effects. Corrections to the impulse approximation have been included using a generalized Gersch-Rodriguez theory that properly takes into account the Fermi statistics of 3He atoms. The microscopic inputs, as the momentum distributions and the two-body density matrices, correspond to a variational (fermi)-hypernetted chain calculation. The agreement with experimental data obtained at $q=23.1$ Å$^{-1}$ is not completely satisfactory, the comparison being difficult due to inconsistencies present in the scattering measurements. The significant differences between the experimental determinations of the 4He condensate fraction and the 3He kinetic energy, and the theoretical results, still remain unsolved.
△ Less
Submitted 24 October, 2000;
originally announced October 2000.
-
Finite temperature dynamic structure function of the free Bose gas
Authors:
F. Mazzanti,
A. Polls
Abstract:
A detailed calculation of the finite temperature dynamic structure function of the free Bose gas is presented and discussed. After a short derivation of the expressions describing the exact response $S(q,ω;T)$ and the Impulse Approximation (IA) $S_{IA}(q,ω;t)$, their main properties and their evolution with q and T are analyzed. The lowest order energy weighted sum rules of both $S(q,ω;T)$ and…
▽ More
A detailed calculation of the finite temperature dynamic structure function of the free Bose gas is presented and discussed. After a short derivation of the expressions describing the exact response $S(q,ω;T)$ and the Impulse Approximation (IA) $S_{IA}(q,ω;t)$, their main properties and their evolution with q and T are analyzed. The lowest order energy weighted sum rules of both $S(q,ω;T)$ and $S_{IA}(q,ω;T)$ are also derived and commented. Finally, the $q\to\infty$ asymptotic behavior is analyzed and discussed in terms of scaling laws.
△ Less
Submitted 12 September, 2000;
originally announced September 2000.
-
Finite Temperature Dynamic Structure Function of the Free Fermi Gas
Authors:
F. Mazzanti,
A. Polls
Abstract:
A detailed discussion of the coherent and incoherent dynamic structure function of the free Fermi gas at finite temperature is presented. Their behavior and evolution with the momentum transfer and the temperature is analyzed, while particular attention is paid to the way in which their relative contribution changes with respect to the $\tilde{T}=0$ case. The influence of thermal effects on the…
▽ More
A detailed discussion of the coherent and incoherent dynamic structure function of the free Fermi gas at finite temperature is presented. Their behavior and evolution with the momentum transfer and the temperature is analyzed, while particular attention is paid to the way in which their relative contribution changes with respect to the $\tilde{T}=0$ case. The influence of thermal effects on the lowest order sum rules of the coherent and incoherent responses is also discussed. Finally, the scaling properties of the responses at high momentum transfer as a function of the temperature are also analyzed.
△ Less
Submitted 9 March, 2000;
originally announced March 2000.