Showing 1–48 of 48 results for author: Mazzanti, F

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.

Comments: 4 pages, 3 figures

Journal ref: Physical Review A 109 033312 (2024)

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.

Comments: 7 pages

Journal ref: New Journal of Physics 26, 013052 (2024)

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.

Comments: 9 pages, 7 figures

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.

Comments: 32 pages, 10 figures

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.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. A 108, 033303 (2023)

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.

Comments: 12 pages, 3 figures

Journal ref: SciPost Phys. 13, 025 (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.

Journal ref: SciPost Phys. 13, 031 (2022)

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.

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.

Comments: 10 pages, 6 Figures

Journal ref: Phys. Rev. A 102, 053308 (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.

Journal ref: Phys. Rev. A 102, 047302 (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.

Journal ref: Physical Review A 100, 063614 (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.

Journal ref: Phys. Rev. A 100, 023608 (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.

Journal ref: Phys. Rev. Research 1, 033088 (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.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. A 99, 043609 (2019)

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.

Journal ref: Phys. Rev. A 98, 053632 (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.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. A 98, 023633 (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.

Comments: 5 pages, 4 figures. Te appear in Physical Review Letters

Journal ref: Phys. Rev. Lett. 119, 250402 (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.

Comments: 14 pages, 13 figures, published version

Journal ref: Phys. Rev. A 96, 033606 (2017)

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.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. A 94, 063630 (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.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 117, 205301 (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.

Comments: 5 figures

Journal ref: Phys. Rev. A 93, 021605 (2016)

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.

Comments: Accepted for publication on "Jornal of Chemical Physics"

Journal ref: J. Chem. Phys. 142, 114114 (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.

Comments: 4 figures

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.

Comments: 4 figures

Journal ref: Phys. Rev. A 90, 043623 (2014)

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.

Comments: 6 pages, 6 figures

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.

Comments: 5 pages, 4 figures

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.

Comments: added figure with finite-size dependence of the energy; comments and references added; title changed

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

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.

Comments: 9 figures, 1 table

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.

Comments: 5 pages containing 5 figures

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.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. A 81, 013612 (2010)

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.

Comments: 4 pages

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.

Comments: 5 figures

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.

Comments: 4 pages, 4 figures

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.

Comments: Submitted to PRA. 7 pages and 8 figures

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.

Comments: Proceedings of the QFS2004 conference in Trento. To appear in JLTP

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.

Comments: 4 pages, accepted in Phys. Rev. Lett

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.

Comments: Accepted in Phys. Rev. A

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.

Comments: To be published in Phys.Lett.A

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.

Comments: 18 pages, 11 figures, to appear in Phys. Rev. B

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.

Comments: To appear in Phys.Lett.A, seven figures

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.

Comments: LaTeX, 6 figures

Journal ref: Phys. Lett. A, Vol.263, pag.416 (1999)

Coherent and Incoherent Dynamic Structure Function of the Free Fermi Gas

Authors: F. Mazzanti, A. Polls, J. Boronat

Abstract: A detailed calculation of the coherent and incoherent dynamic structure functions of the free Fermi gas, starting from their expressions in terms of the one- and semi-diagonal two-body density matrices, is derived and discussed. Their behavior and evolution with the momentum transfer is analyzed, and particular attention is devoted to the contributions that both functions present at negative ene… ▽ More A detailed calculation of the coherent and incoherent dynamic structure functions of the free Fermi gas, starting from their expressions in terms of the one- and semi-diagonal two-body density matrices, is derived and discussed. Their behavior and evolution with the momentum transfer is analyzed, and particular attention is devoted to the contributions that both functions present at negative energies. Finally, an analysis of the energy weighted sum rules satisfied by both responses is also performed. Despite of the simplicity of the model, some of the conclusions can be extended to realistic systems. △ Less

Submitted 9 March, 2000; originally announced March 2000.

Comments: LaTeX, 3 figures

Journal ref: Phys.Lett.A, Vol.220, pag.251 (1996)

Final State Effects in the high q response of $^3$He-$^4$He mixtures

Authors: F. Mazzanti

Abstract: A modified Gersch-Rodriguez formalism describing the leading Final State Effects in the high momentum transfer response of low concentration $^3$He-$^4$He mixtures is presented and discussed. The leading corrections to the Impulse Approximation are expressed in terms of the interatomic potentials and the semidiagonal two-body density matrices of both the mixture and its boson-boson approximation… ▽ More A modified Gersch-Rodriguez formalism describing the leading Final State Effects in the high momentum transfer response of low concentration $^3$He-$^4$He mixtures is presented and discussed. The leading corrections to the Impulse Approximation are expressed in terms of the interatomic potentials and the semidiagonal two-body density matrices of both the mixture and its boson-boson approximation, in which the $^3$He atoms are replaced by bosons of the same mass and at the same partial density. Numerical calculations of the Final State Effects functions of $^3$He and $^4$He are finally presented and discussed. △ Less

Submitted 13 April, 2000; v1 submitted 8 June, 1999; originally announced June 1999.

Comments: RevTeX, 1 figure. To appear in Phys.Lett.A

Description of recent large-$q$ neutron inclusive scattering data from liquid $^4$He

Authors: A. S. Rinat, M. F. Taragin, F. Mazzanti, A. Polls

Abstract: We report dynamical calculations for large-$q$ structure functions of liquid $^4$He at $T$=1.6 and 2.3 K and compare those with recent MARI data. We extend those calculations far beyond the experimental range $q\le 29\Ain$ in order to study the approach of the response to its asymptotic limit for a system with interactions having a strong short-range repulsion. We find only small deviations from… ▽ More We report dynamical calculations for large-$q$ structure functions of liquid $^4$He at $T$=1.6 and 2.3 K and compare those with recent MARI data. We extend those calculations far beyond the experimental range $q\le 29\Ain$ in order to study the approach of the response to its asymptotic limit for a system with interactions having a strong short-range repulsion. We find only small deviations from theoretical $1/q$ behavior, valid for smooth $V$. We repeat an extraction by Glyde et al of cumulant coefficients from data. We argue that fits determine the single atom momentum distribution, but express doubt as to the extraction of meaningful Final State Interaction parameters. △ Less

Submitted 22 June, 1997; originally announced June 1997.

Comments: 37 pages, 13 postscript figs

Report number: WIS-97/01/Jan-PH

A microscopic approach to the response of $^{\bf 3}$He -$^{\bf 4}$He mixtures

Authors: A. Fabrocini, L. Vichi, F. Mazzanti, A. Polls

Abstract: Correlated Basis Function perturbation theory is used to evaluate the zero temperature response $S(q,ω)$ of $^3$He-$^4$He mixtures for inelastic neutron scattering, at momentum transfers $q$ ranging from $1.1$ to $1.7 Å^{-1}$. We adopt a Jastrow correlated ground state and a basis of correlated particle-hole and phonon states. We insert correlated one particle-one hole and one-phonon states to c… ▽ More Correlated Basis Function perturbation theory is used to evaluate the zero temperature response $S(q,ω)$ of $^3$He-$^4$He mixtures for inelastic neutron scattering, at momentum transfers $q$ ranging from $1.1$ to $1.7 Å^{-1}$. We adopt a Jastrow correlated ground state and a basis of correlated particle-hole and phonon states. We insert correlated one particle-one hole and one-phonon states to compute the second order response. The decay of the one-phonon states into two-phonon states is accounted for in boson-boson approximation. The full response is splitted into three partial components $S_{αβ}(q,ω)$, each of them showing a particle-hole bump and a one phonon, delta shaped peak, which stays separated from the multiphonon background. The cross term $S_{34}(q,ω)$ results to be of comparable importance to $S_{33}(q,ω)$ in the particle-hole sector and to $S_{44}(q,ω)$ in the phonon one. Once the one-phonon peak has been convoluted with the experimental broadening, the computed scattering function is in semiquantitative agreement with recent experimental measurements. △ Less

Submitted 22 December, 1995; originally announced December 1995.

Comments: 26 pages, RevTex 3.0, 8 figures available upon request

Report number: ECM-UB-70/95

Final state effects on superfluid $^{\bf 4}$He in the deep inelastic regime

Authors: F. Mazzanti, J. Boronat, A. Polls

Abstract: A study of Final State Effects (FSE) on the dynamic structure function of superfluid $^4$He in the Gersch--Rodriguez formalism is presented. The main ingredients needed in the calculation are the momentum distribution and the semidiagonal two--body density matrix. The influence of these ground state quantities on the FSE is analyzed. A variational form of $ρ_2$ is used, even though simpler forms… ▽ More A study of Final State Effects (FSE) on the dynamic structure function of superfluid $^4$He in the Gersch--Rodriguez formalism is presented. The main ingredients needed in the calculation are the momentum distribution and the semidiagonal two--body density matrix. The influence of these ground state quantities on the FSE is analyzed. A variational form of $ρ_2$ is used, even though simpler forms turn out to give accurate results if properly chosen. Comparison to the experimental response at high momentum transfer is performed. The predicted response is quite sensitive to slight variations on the value of the condensate fraction, the best agreement with experiment being obtained with $n_0=0.082$. Sum rules of the FSE broadening function are also derived and commented. Finally, it is shown that Gersch--Rodriguez theory produces results as accurate as those coming from other more recent FSE theories. △ Less

Submitted 21 December, 1995; originally announced December 1995.

Comments: 20 pages, RevTex 3.0, 11 figures available upon request, to be appear in Phys. Rev. B

Report number: ECM-UB-69/95

Dynamic Structure Function in 3he-4he Mixtures

Authors: J. Boronat, F. Dalfovo, F. Mazzanti, A. Polls

Abstract: Relevant features of the dynamic structure function $S(q,ω)$ in $^3$He-$^4$He mixtures at zero temperature are investigated starting from known properties of the ground state. Sum rules are used to fix rigorous constraints to the different contributions to $S(q,ω)$, coming from $^3$He and $^4$He elementary excitations, as well as to explore the role of the cross term $S^{(3,4)}(q,ω)$. Both the l… ▽ More Relevant features of the dynamic structure function $S(q,ω)$ in $^3$He-$^4$He mixtures at zero temperature are investigated starting from known properties of the ground state. Sum rules are used to fix rigorous constraints to the different contributions to $S(q,ω)$, coming from $^3$He and $^4$He elementary excitations, as well as to explore the role of the cross term $S^{(3,4)}(q,ω)$. Both the low-$q$ (phonon-roton $^4$He excitations and 1p-1h $^3$He excitations) and high-$q$ (deep inelastic scattering) ranges are discussed. △ Less

Submitted 2 March, 1993; originally announced March 1993.

Comments: 29 pages, Plain TeX, 11 figures available by request from dalfovo@itnvax.cineca.it