Showing 1–50 of 176 results for author: Boronat, J

Equation of state of Bose gases beyond the universal regime

Authors: Marti Planasdemunt, Jordi Pera, Jordi Boronat

Abstract: The equation of state of dilute Bose gases, in which the energy only depends on the $s$-wave scattering length, is rather unknown beyond the universal limit. We have carried out a bunch of diffusion Monte Carlo calculations up to gas parameters of $10^{-2}$ to explore how the departure from the universality emerges. Using different model potentials, we calculate the energies of the gas in an exact… ▽ More The equation of state of dilute Bose gases, in which the energy only depends on the $s$-wave scattering length, is rather unknown beyond the universal limit. We have carried out a bunch of diffusion Monte Carlo calculations up to gas parameters of $10^{-2}$ to explore how the departure from the universality emerges. Using different model potentials, we calculate the energies of the gas in an exact way, within some statistical noise, and report the results as a function of the three relevant scattering parameters: the $s$-wave scattering length $a_0$, the $s$-wave effective range $r_0$, and the $p$-wave scattering length $a_1$. If the effective range is not large we observe universality in terms of $a_0$ and $r_0$ up to gas parameters of $10^{-2}$. If $r_0$ grows the regime of universality in these two parameters is reduced and effects of $a_1$ start to be observed. In the $(a_0,r_0)$ universal regime we propose an analytical law that reproduces fairly well the exact energies. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 16 pages

Interaction effects on the itinerant ferromagnetism phase transition

Authors: Jordi Pera, Joaquim Casulleras, Jordi Boronat

Abstract: Itinerant ferromagnetism is one of the most studied quantum phase transitions, the transition point and the nature of this phase transition being widely discussed. In dilute Fermi liquids, this analysis has been carried out up to second-order in the gas parameter, where the results for any spin degeneracy are universal in terms of only the s-wave scattering length $a_0$. We extend this analysis to… ▽ More Itinerant ferromagnetism is one of the most studied quantum phase transitions, the transition point and the nature of this phase transition being widely discussed. In dilute Fermi liquids, this analysis has been carried out up to second-order in the gas parameter, where the results for any spin degeneracy are universal in terms of only the s-wave scattering length $a_0$. We extend this analysis to third-order where energies depend, not only on $a_0$, but also on the s-wave effective range $r_0$ and the p-wave scattering length $a_1$. The introduction in the theory of these new parameters changes the transition point, with respect to the second-order estimation, and also can modify the nature of the phase transition itself. We analyze these interaction effects on the phase transition for different spin values. The emerging phase diagram shows that the type of ferromagnetic transition changes dramatically as a function of $r_0$ and $a_1$ and, importantly, that this classification is not solely determined by the spin value, as happens at second order. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: 9 pages

3He adsorbed on molecular hydrogen surfaces

Authors: M. C. Gordillo, J. Boronat

Abstract: Using a diffusion Monte Carlo (DMC) technique, we calculated the phase diagram of 3He adsorbed on a first solid layer of a molecular hydrogen isotope (H2,HD and D2) on top of graphite. The results are qualitatively similar in all cases: a two-dimensional gas spanning from the infinite dilution limit to a second-layer helium density of 0.048 +/- 0.004 Ang^{-2}. That gas is in equilibrium with a 7/1… ▽ More Using a diffusion Monte Carlo (DMC) technique, we calculated the phase diagram of 3He adsorbed on a first solid layer of a molecular hydrogen isotope (H2,HD and D2) on top of graphite. The results are qualitatively similar in all cases: a two-dimensional gas spanning from the infinite dilution limit to a second-layer helium density of 0.048 +/- 0.004 Ang^{-2}. That gas is in equilibrium with a 7/12 commensurate structure, more stable than any incommensurate triangular solid of similar density. These findings are in reasonably good agreement with available experimental data. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: To appear in J. Chem. Phys

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.

Comments: 6 + 1 pages text + Supplementary Material (appendix). 3 figures + 1 table

Self-Assembled Chains and Solids of Dipolar Atoms in a Multilayer

Authors: G. Guijarro, G. E. Astrakharchik, G. Morigi, J. Boronat

Abstract: We predict that ultracold bosonic dipolar gases, confined within a multilayer geometry, may undergo self-assembling processes, leading to the formation of chain gases and solids. These dipolar chains, with dipoles aligned across different layers, emerge at low densities and resemble phases observed in liquid crystals, such as nematic and smectic phases. We calculate the phase diagram using quantum… ▽ More We predict that ultracold bosonic dipolar gases, confined within a multilayer geometry, may undergo self-assembling processes, leading to the formation of chain gases and solids. These dipolar chains, with dipoles aligned across different layers, emerge at low densities and resemble phases observed in liquid crystals, such as nematic and smectic phases. We calculate the phase diagram using quantum Monte Carlo methods, introducing a newly devised trial wave function designed for describing the chain gas, where dipoles from different layers form chains without in-plane long-range order. We find gas, solid, and chain phases, along with quantum phase transitions between these states. Specifically, we predict a quantum phase transition from a gaseous to a self-ordered phase, which occurs at a critical interlayer distance. Remarkably, in the self-organized phases, the mean interparticle distance can significantly exceed the characteristic length of the interaction potential, yielding solids and chain gases with densities several orders of magnitude lower than those of conventional quantum solids. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: Main text and Supplemental Material: 9 pages, 7 figures. arXiv admin note: text overlap with arXiv:2012.14545

Temperature-induced miscibility of impurities in trapped Bose gases

Authors: G. Pascual, T. Wasak, A. Negretti, G. E. Astrakharchik, J. Boronat

Abstract: We study the thermal properties of impurities embedded in a repulsive Bose gas under a harmonic trapping potential. In order to obtain exact structural properties in this inhomogeneous many-body system, we resort to the path-integral Monte Carlo method. We find that, at low temperatures, a single impurity is expelled to the edges of the bath cloud if the impurity-boson coupling constant is larger… ▽ More We study the thermal properties of impurities embedded in a repulsive Bose gas under a harmonic trapping potential. In order to obtain exact structural properties in this inhomogeneous many-body system, we resort to the path-integral Monte Carlo method. We find that, at low temperatures, a single impurity is expelled to the edges of the bath cloud if the impurity-boson coupling constant is larger than the boson-boson one. However, when the temperature is increased, but still in the Bose-condensed phase, the impurity occupies the core of the trap and, thus, the system becomes miscible. This thermal-induced miscibility transition is also observed for a finite concentration of impurities in this inhomogeneous system. We find that the transition temperature for miscibility depends on the impurity-boson interaction and we indicate a novel nondestructive method to measure the temperature of a system based on the studied phenomenon. △ Less

Submitted 16 April, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 7 pages, 7 figures

Journal ref: Phys. Rev. Research 6, L022014, 2024

Ring solids and supersolids in spherical shell-shaped dipolar Bose-Einstein condensates

Authors: Juan Sánchez-Baena, Raúl Bombín, Jordi Boronat

Abstract: We study the interplay between the anisotropy of the dipole-dipole interaction and confinement in a curved geometry by means of the extended Gross-Pitaevskii equation, which allows us to characterize the ground state of a dipolar Bose gas under the confinement of a bubble trapping potential. We do so in terms of the scattering length a and the number of particles. We observe the emergence of a wid… ▽ More We study the interplay between the anisotropy of the dipole-dipole interaction and confinement in a curved geometry by means of the extended Gross-Pitaevskii equation, which allows us to characterize the ground state of a dipolar Bose gas under the confinement of a bubble trapping potential. We do so in terms of the scattering length a and the number of particles. We observe the emergence of a wide variety of dipolar solids, consisting on arrangements of different number of droplets along a ring over the equator of the spherical shell confinement. We also show that the transition between the different phases of the system can be engineered by varying a, the number of particles or the radius of the trap, parameters which can be experimentally tuned. Finally, we show the importance of working in microgravity conditions as gravity unstabilizes the observed dipolar solids. △ Less

Submitted 11 April, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: New version, 10 pages, 10 figures

Reply to "Comment on Absence of off-diagonal long-range order in hcp $^{\bf 4}$He dislocation cores"

Authors: Maurice de Koning, Wei Cai, Claudio Cazorla, Jordi Boronat

Abstract: We offer a Reply to Document arXiv:2303.12174 by M. Boninsegni et al., commenting on arXiv:2301.02854 . We offer a Reply to Document arXiv:2303.12174 by M. Boninsegni et al., commenting on arXiv:2301.02854 . △ Less

Submitted 1 December, 2023; originally announced December 2023.

Squeezing Bose-Bose liquid drops

Authors: A. Sanuy, R. Barač, P. Stipanović, L. Vranješ Markić, J. Boronat

Abstract: We explore ultradilute Bose-Bose liquid droplets squeezed by an external harmonic potential in one spatial direction. Our theoretical study is based on a functional that is built using quantum Monte Carlo results of the bulk phase and incorporates finite-range effects. A characteristic feature of these drops is the existence of a critical atom number, that is the minimum number of particles to hav… ▽ More We explore ultradilute Bose-Bose liquid droplets squeezed by an external harmonic potential in one spatial direction. Our theoretical study is based on a functional that is built using quantum Monte Carlo results of the bulk phase and incorporates finite-range effects. A characteristic feature of these drops is the existence of a critical atom number, that is the minimum number of particles to have a many-body bound state. We report results on the critical atom numbers for different magnetic fields and applying confinement strengths towards a quasi-two-dimensional setup. In the regime where the local density approximation is expected to be valid, we find that the critical atom number decreases linearly with the harmonic oscillator length of the confining potential. With the largest squeezing explored in our work, we predict stable drops at the level of one thousand atoms. Our functional reduces the critical numbers for any confinement and applied magnetic field with respect to the estimations based on the Lee-Huang-Yang model. We observe saturated drops when the number of atoms in the drop is much larger than the critical value, their central density being higher for the quantum Monte Carlo functional than for the Lee-Huang-Yang one. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 7 pages, 6 figures

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)

Phases of 4He and H2 adsorbed on doped graphene

Authors: M. C. Gordillo, J. Boronat

Abstract: The influence of attractive boron impurities, embedded on a graphene sheet, on the phase diagrams of $^4$He and H$_2$ adsorbed on top was studied using the diffusion Monte Carlo method. The doping of graphene was made by distributing the boron atoms following the same pattern found in an experimentally synthesized substrate. Our results show that while the different incommensurate solid phases of… ▽ More The influence of attractive boron impurities, embedded on a graphene sheet, on the phase diagrams of $^4$He and H$_2$ adsorbed on top was studied using the diffusion Monte Carlo method. The doping of graphene was made by distributing the boron atoms following the same pattern found in an experimentally synthesized substrate. Our results show that while the different incommensurate solid phases of both adsorbates remain largely unchanged after doping, the liquid/gas equations of state are significanty different from the ones on pristine graphene. Doping graphene produces new translationally invariant stable phases for $^4$He, depending on the concentration of boron impurities, but makes the H$_2$ ground state to remain solid. In addition, several new registered phases appear for both adsorbates. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 8 pages

Journal ref: Phys. Rev. B 108, 134505 (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.

Comments: 7 pages

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

Phases of $^4$He and H$_2$ adsorbed on a single carbon nanotube

Authors: M. C. Gordillo, J. Boronat

Abstract: Using a diffusion Monte Carlo (DMC) technique, we calculated the phase diagrams of $^4$He and H$_2$ adsorbed on a single (5,5) carbon nanotube, one of the narrowest that can be obtained experimentally. For a single monolayer, when the adsorbate density increases, both species undergo a series of first order solid-solid phase transitions between incommensurate arrangements. Remarkably, the $^4$He l… ▽ More Using a diffusion Monte Carlo (DMC) technique, we calculated the phase diagrams of $^4$He and H$_2$ adsorbed on a single (5,5) carbon nanotube, one of the narrowest that can be obtained experimentally. For a single monolayer, when the adsorbate density increases, both species undergo a series of first order solid-solid phase transitions between incommensurate arrangements. Remarkably, the $^4$He lowest-density solid phase shows supersolid behavior in contrast with the normal solid that we found for H$_2$. The nature of the second-layer is also different for both adsorbates. Contrarily to what happens on graphite, the second-layer of $^4$He on that tube is a liquid, at least up to the density corresponding to a third-layer promotion on a flat substrate. However, the second-layer of H$_2$ is a solid that, at its lowest stable density, has a small but observable superfluid fraction. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Journal ref: Phys. Rev. B 107, 174518 (2023)

H2 superglass on an amorphous carbon substrate

Authors: M. C. Gordillo, J. Boronat

Abstract: The phase diagram of a para-H2 monolayer absorbed on a experimentally syntetized amorphous carbon sheet was calculated using a diffusion Monte Carlo technique. We found that the ground state of that system changed drastically from a perfectly flat substrate to a situation in which the carbon atoms were allowed a certain degree of disorder in the $z$ direction. In the first case, at zero pressure w… ▽ More The phase diagram of a para-H2 monolayer absorbed on a experimentally syntetized amorphous carbon sheet was calculated using a diffusion Monte Carlo technique. We found that the ground state of that system changed drastically from a perfectly flat substrate to a situation in which the carbon atoms were allowed a certain degree of disorder in the $z$ direction. In the first case, at zero pressure we have a glass of density 0.056 $\pm$ 0.003 Å$^{-2}$ in equilibrium with an incommensurate solid of 0.068 $\pm$ 0.002 Å$^{-2}$. At the equilibrium density, the glass was found to have a tiny, but non-negligible superfluid faction of less than 1 \% (0.44 $\pm$ 0.05 \%). In the $z$-disordered substrate, we observe a significant enhancement of the superfluid fraction in the glass phase as well as a smaller but not zero value in the incommensurate crystal. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Journal ref: Phys. Rev. B 107, L060505 (2023)

Thermal-induced Local Imbalance in Repulsive Binary Bose Mixtures

Authors: Gerard Pascual, Gabriele Spada, Sebastiano Pilati, Stefano Giorgini, Jordi Boronat

Abstract: We study repulsive two-component Bose mixtures with equal populations and confined in a finite-size box through path-integral Monte Carlo simulations. For different values of the $s$-wave scattering length of the interspecies potential, we calculate the local population imbalance in a region of fixed volume inside the box at different temperatures. We find two different behaviors: for phase-separa… ▽ More We study repulsive two-component Bose mixtures with equal populations and confined in a finite-size box through path-integral Monte Carlo simulations. For different values of the $s$-wave scattering length of the interspecies potential, we calculate the local population imbalance in a region of fixed volume inside the box at different temperatures. We find two different behaviors: for phase-separated states at $T=0$, thermal effects induce a diffusion process which reduces the local imbalance whereas, for miscible states at $T=0$, a maximum in the local population imbalance appears at a certain temperature, below the critical one. We show that this intriguing behavior is strongly related to the bunching effect associated with the Bose-Einstein statistics of the particles in the mixture and to an unexpected behavior of the cross pair distribution function not reported before. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 7 pages, 7 figures

Journal ref: Phys. Rev. Research 5, L032041 (2023)

Low-energy scattering parameters: A theoretical derivation of the effective range and scattering length for arbitrary angular momentum

Authors: Jordi Pera, Jordi Boronat

Abstract: The most important parameters in the study of low-energy scattering are the s-wave and p-wave scattering lengths and the s-wave effective range. We solve the scattering problem and find two useful formulas for the scattering length and the effective range for any angular momentum, as long as the Wigner threshold law holds. Using that formalism, we obtain a set of useful formulas for the angular-mo… ▽ More The most important parameters in the study of low-energy scattering are the s-wave and p-wave scattering lengths and the s-wave effective range. We solve the scattering problem and find two useful formulas for the scattering length and the effective range for any angular momentum, as long as the Wigner threshold law holds. Using that formalism, we obtain a set of useful formulas for the angular-momentum scattering parameters of four different model potentials: hard-sphere, soft-sphere, spherical well, and well-barrier potentials. The behavior of the scattering parameters close to Feshbach resonances is also analyzed. Our derivations can be useful as hands-on activities for learning scattering theory. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Main document with two supplementary materials

Journal ref: American Journal of Physics 91, 90 (2023)

Thermal fading of the $1/k^4$-tail of the momentum distribution induced by the hole anomaly

Authors: Giulia De Rosi, Grigori E. Astrakharchik, Maxim Olshanii, Jordi Boronat

Abstract: We study the thermal behavior of correlations in a one-dimensional Bose gas with tunable interaction strength, crossing from weakly-repulsive to Tonks-Girardeau regime. A reference temperature in this system is that of the hole anomaly, observed as a peak in the specific heat and a maximum in the chemical potential. We find that at large momenta $k$ and temperature above the anomaly threshold, the… ▽ More We study the thermal behavior of correlations in a one-dimensional Bose gas with tunable interaction strength, crossing from weakly-repulsive to Tonks-Girardeau regime. A reference temperature in this system is that of the hole anomaly, observed as a peak in the specific heat and a maximum in the chemical potential. We find that at large momenta $k$ and temperature above the anomaly threshold, the tail $\mathcal{C}/k^4$ of the momentum distribution (proportional to the Tan contact $\mathcal{C}$) is screened by the $1/|k|^3$-term due to a dramatic thermal increase of the internal energy emerging from the thermal occupation of spectral excitation states. The same fading is consistently revealed in the behavior at short distances $x$ of the one-body density matrix (OBDM) where the $|x|^3$-dependence disappears for temperatures above the anomaly. We obtain a new general analytic tail for the momentum distribution and a minimum $k$ fixing its validity range, both calculated with exact Bethe-Ansatz method and valid in all interaction and thermal regimes, crossing from the quantum to the classical gas limit. Our predictions are confirmed by comparison with ab-initio Path Integral Monte Carlo calculations for the momentum distribution and the OBDM exploring a wide range of interaction strength and temperature. Our results unveil a novel connection between excitations and correlations. We expect them to be of interest to any cold atomic, nuclear, solid-state, electronic and spin system exhibiting an anomaly or a thermal second-order phase transition. △ Less

Submitted 20 February, 2024; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: Main Text: 7 pages, 4 figures; Supplemental Material: 10 pages, 8 figures

Journal ref: Phys. Rev. A 109, L031302 (2024)

Correlation properties of a one-dimensional repulsive Bose gas at finite temperature

Authors: Giulia De Rosi, Riccardo Rota, Grigori E. Astrakharchik, Jordi Boronat

Abstract: We present a comprehensive study shedding light on how thermal fluctuations affect correlations in a Bose gas with contact repulsive interactions in one spatial dimension. The pair correlation function, the static structure factor, and the one-body density matrix are calculated as a function of the interaction strength and temperature with the exact ab-initio Path Integral Monte Carlo method. We e… ▽ More We present a comprehensive study shedding light on how thermal fluctuations affect correlations in a Bose gas with contact repulsive interactions in one spatial dimension. The pair correlation function, the static structure factor, and the one-body density matrix are calculated as a function of the interaction strength and temperature with the exact ab-initio Path Integral Monte Carlo method. We explore all possible gas regimes from weak to strong interactions and from low to high temperatures. We provide a detailed comparison with a number of theories, such as perturbative (Bogoliubov and decoherent classical), effective (Luttinger liquid) and exact (ground-state and thermal Bethe Ansatz) ones. Our Monte Carlo results exhibit an excellent agreement with the tractable limits and provide a fundamental benchmark for future observations which can be achieved in atomic gases, cavity quantum-electrodynamic and superconducting-circuit platforms. △ Less

Submitted 6 April, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: Main Text: 13 pages, 9 figures. Appendix: 1 page

Journal ref: New J. Phys. 25 043002 (2023)

Excitonic Insulator to Superconductor Phase Transition in Ultra-Compressed Helium

Authors: Cong Liu, Ion Errea, Chris Pickard, Lewis J. Conway, Bartomeu Monserrat, Yue-Wen Fang, Chi Ding, Qing Lu, Jian Sun, Jordi Boronat, Claudio Cazorla

Abstract: Helium, the second most abundant element in the universe, exhibits an extremely large electronic band gap of about $20$ eV at low pressures ($\le 0.1$ GPa). While the metallization pressure of hcp helium has been accurately predicted, thus far little attention has been paid to the specific mechanisms driving the band-gap closure and electronic properties of this quantum crystal in the terapascal r… ▽ More Helium, the second most abundant element in the universe, exhibits an extremely large electronic band gap of about $20$ eV at low pressures ($\le 0.1$ GPa). While the metallization pressure of hcp helium has been accurately predicted, thus far little attention has been paid to the specific mechanisms driving the band-gap closure and electronic properties of this quantum crystal in the terapascal regime (1 TPa $= 1,000$ GPa). Here, we employ state-of-the-art density functional theory and many-body perturbation theory calculations to fill up this knowledge gap. It is found that prior to reaching metallicity bulk solid helium becomes an excitonic insulator (EI), an exotic state of matter typically observed in low-dimensional systems in which electrostatically bound electron-hole pairs form spontaneously. Furthermore, it is shown that electron-phonon coupling (EPC) is significantly enhanced across the EI to metal phase transition as signaled by prominent phonon softening and giant EPC strength values ($λ\sim 10-100$) estimated at specific reciprocal space points. Accordingly, we predict metallic helium to be a superconductor with a critical temperature of $\approx 30$ K at $20$ TPa and of $\approx 100$ K at $100$ TPa. These unforeseen phenomena have important consequences on the elastic, thermodynamic and transport properties of metallic helium hence may be critical for improving our fundamental understanding and modelling of celestial bodies. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 25 pages, 18 figures

Absence of off-diagonal long-range order in hcp $^{\bf 4}$He dislocation cores

Authors: Maurice de Koning, Wei Cai, Claudio Cazorla, Jordi Boronat

Abstract: The mass transport properties along dislocation cores in hcp $^4$He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation c… ▽ More The mass transport properties along dislocation cores in hcp $^4$He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation core structures and their off-diagonal long-range order properties. It is found that the Bose-Einstein condensate fraction of such defective $^4$He systems is practically null ($\le 10^{-6}$), just as in the bulk defect-free crystal. These results provide compelling evidence for the absence of intrinsic superfluidity in dislocation cores in hcp $^4$He and challenge the superfluid dislocation-network interpretation of the mass-flux-experiment observations, calling for further experimental investigation. △ Less

Submitted 7 January, 2023; originally announced January 2023.

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

Theoretical approaches to liquid Helium

Authors: Jordi Boronat

Abstract: In this article, we review the main theoretical methods applied to the study of liquid Helium adopting a microscopic approach, that is, starting from the many-particle Hamiltonian of the system. Following an introduction on the first early approaches, we discuss two main issues. In the first one, we report the main ingredients of a theoretical approach based on the variational method, with the dis… ▽ More In this article, we review the main theoretical methods applied to the study of liquid Helium adopting a microscopic approach, that is, starting from the many-particle Hamiltonian of the system. Following an introduction on the first early approaches, we discuss two main issues. In the first one, we report the main ingredients of a theoretical approach based on the variational method, with the discussion of the high accuracy obtained with them and the related progress in the design of highly accurate trial wave functions. In the second part, the main stochastic methods used in this study are briefly discussed. Altogether reflects the strong links between the study of liquid Helium and the progress in the development of new and extremely powerful approaches to deal with one of the most strongly quantum many-body systems in Nature. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: To appear in Encyclopedia of Condensed Matter Physics, 2nd edition

Phase separation in binary Bose mixtures at finite temperature

Authors: G. Spada, L. Parisi, G. Pascual, N. G. Parker, T. P. Billam, S. Pilati, J. Boronat, S. Giorgini

Abstract: We investigate the magnetic behavior of finite-temperature repulsive two-component Bose mixtures by means of exact path-integral Monte-Carlo simulations. Novel algorithms are implemented for the free energy and the chemical potential of the two components. Results on the magnetic susceptibility suggest that the conditions for phase separation are not modified from the zero temperature case. This c… ▽ More We investigate the magnetic behavior of finite-temperature repulsive two-component Bose mixtures by means of exact path-integral Monte-Carlo simulations. Novel algorithms are implemented for the free energy and the chemical potential of the two components. Results on the magnetic susceptibility suggest that the conditions for phase separation are not modified from the zero temperature case. This contradicts previous predictions based on approximate theories. We also determine the temperature dependence of the chemical potential and the contact parameters for experimentally relevant balanced mixtures. △ Less

Submitted 20 October, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: v2: 19 pages, 10 figures, matches published version; v1: 6 pages, 5 figures plus supplemental material of 8 pages, 5 figures)

Journal ref: SciPost Phys. 15, 171 (2023)

Van der Waals five-body size-energy universality

Authors: Petar Stipanović, Leandra Vranješ Markić, Jordi Boronat

Abstract: A universal relationship between scaled size and scaled energy is explored in five-body self-bound quantum systems. The ground-state binding energy and structure properties are obtained by means of the diffusion Monte Carlo method. We use pure estimators to eliminate any residual bias in the estimation of the cluster size. Strengthening the inter-particle interaction, we extend the exploration fro… ▽ More A universal relationship between scaled size and scaled energy is explored in five-body self-bound quantum systems. The ground-state binding energy and structure properties are obtained by means of the diffusion Monte Carlo method. We use pure estimators to eliminate any residual bias in the estimation of the cluster size. Strengthening the inter-particle interaction, we extend the exploration from the halo region to classical systems. Universal scaled size-scaled energy line, which does not depend on the short-range potential details and binding strength, is found for homogeneous pentamers with interaction potentials decaying at long range predominantly as $r^{-6}$. For mixed pentamers, we discuss under which conditions the universal line can approximately describe the size-energy ratio. Our data is compatible with generalized Tjon lines, which assume a linear dependence between the binding energy of the pentamers and the one of tetramers, when both are divided by the trimer energies. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: This version of the article has been accepted for publication, after peer review but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections. The Version of Record is available online at: https://dx.doi.org/10.1038/s41598-022-13630-2

ACM Class: J.2; I.6; G.1.8

Journal ref: Scientific Reports 12, 10368 (2022)

Beyond universality in repulsive SU(N) Fermi gases

Authors: Jordi Pera, Joaquim Casulleras, Jordi Boronat

Abstract: Itinerant ferromagnetism in dilute Fermi gases is predicted to emerge at values of the gas parameter where second-order perturbation theory is not accurate enough to properly describe the system. We have revisited perturbation theory for SU(N) fermions and derived its generalization up to third order both in terms of the gas parameter and the polarization. Our results agree satisfactorily with qua… ▽ More Itinerant ferromagnetism in dilute Fermi gases is predicted to emerge at values of the gas parameter where second-order perturbation theory is not accurate enough to properly describe the system. We have revisited perturbation theory for SU(N) fermions and derived its generalization up to third order both in terms of the gas parameter and the polarization. Our results agree satisfactorily with quantum Monte Carlo results for hard-sphere and soft-sphere potentials for $S = 1/2$. Although the nature of the phase transition depends on the interaction potential, we find that for a hard-sphere potential a phase transition is guaranteed to occur. While for $S= 1/2$ we observe a quasi-continuous transition, for spins $3/2$ and $5/2$, a first-order phase transition is found. For larger spins, a double transition (combination of continuous and discontinuous) occurs. The critical density reduces drastically when the spin increases, making the phase transition more accessible to experiments with ultracold dilute Fermi gases. Estimations for Fermi gases of Yb and Sr with spin $5/2$ and $9/2$, respectively, are reported. △ Less

Submitted 30 May, 2024; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 24 pages. Some minor corrections introduced

Itinerant ferromagnetism in dilute SU(N) Fermi gases

Authors: Jordi Pera, Joaquim Casulleras, Jordi Boronat

Abstract: We present exact analytic results for the energy of a SU(N) repulsive Fermi gas as a function of the spin-channel occupation at second order in the gas parameter. This is an extension of an old result that now incorporates the degree of polarization of the system. Therefore, the magnetic properties of the gas can be obtained, free from numerical uncertainties. For spin 1/2 we find that second-orde… ▽ More We present exact analytic results for the energy of a SU(N) repulsive Fermi gas as a function of the spin-channel occupation at second order in the gas parameter. This is an extension of an old result that now incorporates the degree of polarization of the system. Therefore, the magnetic properties of the gas can be obtained, free from numerical uncertainties. For spin 1/2 we find that second-order corrections change the itinerant ferromagnetic transition from continuous to first-order. Instead, for spin larger than 1/2 the phase transition is always of first-order type. The transition critical density reduces when the spin increases, making the phase transition more accessible to experiments with ultracold dilute Fermi gases. Estimations for Fermi gases of Yb and Sr with spin 5/2 and 9/2, respectively, are reported. △ Less

Submitted 2 November, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: 12 pages, 10 figures. Minor corrections

Journal ref: SciPost Phys. 14, 038 (2023)

Supersolidity in the second layer of $para$-H$_2$ adsorbed on graphite

Authors: M. C. Gordillo, J. Boronat

Abstract: We calculated the phase diagram of the second layer of $para$-H$_2$ adsorbed on graphite using quantum Monte Carlo methods. The second layer shows an incommensurate triangular crystal structure. By using a symmetric wave function, that makes possible molecule exchanges, we observed that this nearly two-dimensional crystal shows a finite superfluid density around a total density of 0.1650 Å$^{-2}$.… ▽ More We calculated the phase diagram of the second layer of $para$-H$_2$ adsorbed on graphite using quantum Monte Carlo methods. The second layer shows an incommensurate triangular crystal structure. By using a symmetric wave function, that makes possible molecule exchanges, we observed that this nearly two-dimensional crystal shows a finite superfluid density around a total density of 0.1650 Å$^{-2}$. The superfluid fraction of this supersolid phase was found to be small, 0.41$\pm$0.05 \%, but still experimentally accessible. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 6 pages

Journal ref: Physical Review B 105, 094501 (2022)

Mixtures of dipolar gases in two dimensions: a quantum Monte Carlo study

Authors: Sergi Pradas, Jordi Boronat

Abstract: We studied the miscibility of two dipolar quantum gases in the limit of zero temperature. The system under study is composed by a mixture of two Bose gases with dominant dipolar interaction in a two-dimensional harmonic confinement. The dipolar moments are considered all to be perpendicular to the plane, turning the dipolar potential in a purely repulsive and isotropic model. Our analysis is carri… ▽ More We studied the miscibility of two dipolar quantum gases in the limit of zero temperature. The system under study is composed by a mixture of two Bose gases with dominant dipolar interaction in a two-dimensional harmonic confinement. The dipolar moments are considered all to be perpendicular to the plane, turning the dipolar potential in a purely repulsive and isotropic model. Our analysis is carried out by using the diffusion Monte Carlo method which allows for an exact solution to the many-body problem within some statistical noise. Our results show that the miscibility between the two species is rather constrained as a function of the relative dipolar moments and masses of the two components. A narrow regime is predicted where both species mix and we introduce an adimensional parameter whose value predicts quite accurately the miscibility of the two dipolar gases. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 14 pages, 7 figures

Dilute quantum liquid in a K-Rb Bose mixture

Authors: Viktor Cikojević, Elena Poli, Francesco Ancilotto, Leandra Vranješ Markić, Jordi Boronat

Abstract: A quantum liquid in a heterogeneous mixture of $^{41}$K and $^{87}$Rb atoms is studied using the diffusion Monte Carlo method and Density Functional Theory. The perturbative Lee-Huang-Yang term for a heterogeneous mixture is verified and it is proved to be valid only near the gas-liquid transition. Based on the equations of state of the bulk mixture, calculated with diffusion Monte Carlo, extensio… ▽ More A quantum liquid in a heterogeneous mixture of $^{41}$K and $^{87}$Rb atoms is studied using the diffusion Monte Carlo method and Density Functional Theory. The perturbative Lee-Huang-Yang term for a heterogeneous mixture is verified and it is proved to be valid only near the gas-liquid transition. Based on the equations of state of the bulk mixture, calculated with diffusion Monte Carlo, extensions to Lee-Huang-Yang corrected mean-field energy functionals (MF+LHY) are presented. Using Density Functional Theory, a systematic comparison between different functionals is performed, focusing on the critical atom number, surface tension, surface width, Tolman length, and compressibility. These results are given as a function of the inter-species interaction strength, within the stability domain of the liquid mixture. △ Less

Submitted 13 July, 2021; originally announced July 2021.

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

Ultradilute quantum liquid of dipolar atoms in a bilayer

Authors: G. Guijarro, G. E. Astrakharchik, J. Boronat

Abstract: We show that ultradilute quantum liquids can be formed with ultracold bosonic dipolar atoms in a bilayer geometry. Contrary to previous realizations of ultradilute liquids, there is no need for stabilizing the system with an additional repulsive short-range potential. The advantage of the proposed system is that dipolar interactions on their own are sufficient for creation of a self-bound state an… ▽ More We show that ultradilute quantum liquids can be formed with ultracold bosonic dipolar atoms in a bilayer geometry. Contrary to previous realizations of ultradilute liquids, there is no need for stabilizing the system with an additional repulsive short-range potential. The advantage of the proposed system is that dipolar interactions on their own are sufficient for creation of a self-bound state and no additional short-range potential is needed for the stabilization. We perform quantum Monte Carlo simulations and find a rich ground-state phase diagram that contains quantum phase transitions between liquid, solid, atomic gas, and molecular gas phases. The stabilization mechanism of the liquid phase is consistent with the microscopic scenario in which the effective dimer-dimer attraction is balanced by an effective three-dimer repulsion. The equilibrium density of the liquid, which is extremely small, can be controlled by the interlayer distance. From the equation of state, we extract the spinodal density, below which the homogeneous system breaks into droplets. Our results offer a new example of a two-dimensional interacting dipolar liquid in a clean and highly controllable setup. △ Less

Submitted 11 February, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: 7 pages, 5 figures; Supplemental Materials included

Journal ref: Phys. Rev. Lett. 128, 063401 (2022)

Quasiparticle nature of the Bose polaron at finite temperature

Authors: Gerard Pascual, Jordi Boronat

Abstract: The Bose polaron has attracted theoretical and experimental interest because the mobile impurity is surrounded by a bath that undergoes a superfluid-to-normal phase transition. Although many theoretical works have studied this system in its ground state, only few analyze its behavior at finite temperature. We have studied the effect of temperature on a Bose polaron system performing ab-initio Path… ▽ More The Bose polaron has attracted theoretical and experimental interest because the mobile impurity is surrounded by a bath that undergoes a superfluid-to-normal phase transition. Although many theoretical works have studied this system in its ground state, only few analyze its behavior at finite temperature. We have studied the effect of temperature on a Bose polaron system performing ab-initio Path Integral Monte Carlo simulations. This method is able to approach the critical temperature without losing accuracy, in contrast with perturbative approximations. We have calculated the polaron energy for the repulsive and attractive branches and we have observed an asymmetric behavior between the two branches. When the potential is repulsive, the polaron energy decreases when the temperature increases, and contrariwise for the attractive branch. Our results for the effective mass and the dynamical structure factor of the polaron show unambiguously that its quasiparticle nature disappears close to the critical temperature, in agreement with recent experimental findings. Finally, we have also estimated the fraction of bosons in the condensate as well as the superfluid fraction, and we have concluded that the impurity hinders the condensation of the rest of bosons. △ Less

Submitted 11 November, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: 9 pages, 6 figures

Hole-induced anomaly in the thermodynamic behavior of a one-dimensional Bose gas

Authors: Giulia De Rosi, Riccardo Rota, Grigori E. Astrakharchik, Jordi Boronat

Abstract: We reveal an intriguing anomaly in the temperature dependence of the specific heat of a one-dimensional Bose gas. The observed peak holds for arbitrary interaction and remembers a superfluid-to-normal phase transition in higher dimensions, but phase transitions are not allowed in one dimension. The presence of the anomaly signals a region of unpopulated states which behaves as an energy gap and is… ▽ More We reveal an intriguing anomaly in the temperature dependence of the specific heat of a one-dimensional Bose gas. The observed peak holds for arbitrary interaction and remembers a superfluid-to-normal phase transition in higher dimensions, but phase transitions are not allowed in one dimension. The presence of the anomaly signals a region of unpopulated states which behaves as an energy gap and is located below the hole branch in the excitation spectrum. The anomaly temperature is found to be of the same order of the energy of the maximum of the hole branch. We rely on the Bethe Ansatz to obtain the specific heat exactly and provide interpretations of the analytically tractable limits. The dynamic structure factor is computed with the Path Integral Monte Carlo method for the first time. We notice that at temperatures similar to the anomaly threshold, the energy of the thermal fluctuations become comparable with the maximal hole energy, leading to a qualitative change in the structure of excitations. This excitation pattern experiences the breakdown of the quasi-particle description for any value of the interaction strength at the anomaly, similarly to any superfluid phase transition at the critical temperature. We provide indications for future observations and how the hole anomaly can be employed for in-situ thermometry, identifying different collisional regimes and understanding other anomalies in atomic, solid-state, electronic, spin-chain and ladder systems. △ Less

Submitted 23 May, 2022; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: Main Text: 16 pages, 6 figures. Appendices: 11 pages, 4 figures. Open-access datasets can be found at 10.5821/data-2117-353128-1

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

Dynamics of equilibration and collisions in ultradilute quantum droplets

Authors: Viktor Cikojević, Leandra Vranješ Markić, Martí Pi, Manuel Barranco, Francesco Ancilotto, Jordi Boronat

Abstract: Employing time-dependent density-functional theory, we have studied dynamical equilibration and binary head-on collisions of quantum droplets made of a $^{39}$K-$^{39}$K Bose mixture. The phase space of collision outcomes is extensively explored by performing fully three-dimensional calculations with effective single-component QMC based and two-components LHY-corrected mean-field functionals. We e… ▽ More Employing time-dependent density-functional theory, we have studied dynamical equilibration and binary head-on collisions of quantum droplets made of a $^{39}$K-$^{39}$K Bose mixture. The phase space of collision outcomes is extensively explored by performing fully three-dimensional calculations with effective single-component QMC based and two-components LHY-corrected mean-field functionals. We exhaustively explored the important effect -- not considered in previous studies -- of the initial population ratio deviating from the optimal mean-field value $N_2/N_1 = \sqrt{a_{11} / a_{22}}$. Both stationary and dynamical calculations with an initial non-optimal concentration ratio display good agreement with experiments. Calculations including three-body losses acting only on the $\left|F, m_{F}\right\rangle=|1,0\rangle$ state show dramatic differences with those obtained with the three-body term acting on the total density. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Quantum halo states in two-dimensional dipolar clusters

Authors: G. Guijarro, G. E. Astrakharchik, J. Boronat

Abstract: A halo is an intrinsically quantum object defined as a bound state of a spatial size which extends deeply into the classically forbidden region. Previously, halos have been observed in bound states of two and less frequently of three atoms. Here, we propose a realization of halo states containing as many as six atoms. We report the binding energies, pair correlation functions, spatial distribution… ▽ More A halo is an intrinsically quantum object defined as a bound state of a spatial size which extends deeply into the classically forbidden region. Previously, halos have been observed in bound states of two and less frequently of three atoms. Here, we propose a realization of halo states containing as many as six atoms. We report the binding energies, pair correlation functions, spatial distributions, and sizes of few-body clusters composed by bosonic dipolar atoms in a bilayer geometry. We find two very distinct halo structures, for large interlayer separation the halo structure is roughly symmetric and we discover an unusual highly anisotropic shape of halo states close to the unbinding threshold. Our results open avenues of using ultracold gases for the experimental realization of halos with the largest number of atoms ever predicted before. △ Less

Submitted 15 January, 2021; originally announced January 2021.

Comments: 8 pages, 4 figures

Journal ref: Scientific Reports 11, 19437 (2021)

Structural Superfluid-Mott Insulator Transition for a Bose Gas in Multi-Rods

Authors: Omar Abel Rodríguez-López, M. A. Solís, J. Boronat

Abstract: We report on a novel structural Superfluid-Mott Insulator (SF-MI) quantum phase transition for an interacting one-dimensional Bose gas within permeable multi-rod lattices, where the rod lengths are varied from zero to the lattice period length. We use the ab-initio diffusion Monte Carlo method to calculate the static structure factor, the insulation gap, and the Luttinger parameter, which we use t… ▽ More We report on a novel structural Superfluid-Mott Insulator (SF-MI) quantum phase transition for an interacting one-dimensional Bose gas within permeable multi-rod lattices, where the rod lengths are varied from zero to the lattice period length. We use the ab-initio diffusion Monte Carlo method to calculate the static structure factor, the insulation gap, and the Luttinger parameter, which we use to determine if the gas is a superfluid or a Mott insulator. For the Bose gas within a square Kronig-Penney (KP) potential, where barrier and well widths are equal, the SF-MI coexistence curve shows the same qualitative and quantitative behavior as that of a typical optical lattice with equal periodicity but slightly larger height. When we vary the width of the barriers from zero to the length of the potential period, keeping the height of the KP barriers, we observe a new way to induce the SF-MI phase transition. Our results are of significant interest, given the recent progress on the realization of optical lattices with a subwavelength structure that would facilitate their experimental observation. △ Less

Submitted 10 January, 2021; originally announced January 2021.

Comments: 6 pages, 4 figures, accepted for publication in Phys. Rev. A

Trapped Bose-Bose mixtures at finite temperature: a quantum Monte Carlo approach

Authors: K. Dželalija, V. Cikojević, J. Boronat, L. Vranješ Markić

Abstract: We study thermal properties of a trapped Bose-Bose mixture in a dilute regime using quantum Monte Carlo methods. Our main aim is to investigate the dependence of the superfluid density and the condensate fraction on temperature, for the mixed and separated phases. To this end, we use the diffusion Monte Carlo method, in the zero-temperature limit, and the path-integral Monte Carlo method for finit… ▽ More We study thermal properties of a trapped Bose-Bose mixture in a dilute regime using quantum Monte Carlo methods. Our main aim is to investigate the dependence of the superfluid density and the condensate fraction on temperature, for the mixed and separated phases. To this end, we use the diffusion Monte Carlo method, in the zero-temperature limit, and the path-integral Monte Carlo method for finite temperatures. The results obtained are compared with solutions of the coupled Gross-Pitaevskii equations for the mixture at zero temperature. We notice the existence of an anisotropic superfluid density in some phase-separated mixtures. Our results also show that the temperature evolution of the superfluid density and condensate fraction is slightly different, showing noteworthy situations where the superfluid fraction is smaller than the condensate fraction. △ Less

Submitted 8 November, 2020; originally announced November 2020.

Comments: accepted for publication in Phys. Rev. A

Journal ref: Phys. Rev. A 102, 063304 (2020)

Finite-range effects in the two-dimensional repulsive Fermi polaron

Authors: Raúl Bombín, Viktor Cikojević, Juan Sánchez-Baena, Jordi Boronat

Abstract: We study the repulsive Fermi polaron in a two-component, two-dimensional system of fermionic atoms inspired by the results of a recent experiment with $^{173}$Yb atoms [N. Darkwah Oppong \textit{et al.}, Phys. Rev. Lett. \textbf{122}, 193604 (2019)]. We use the diffusion Monte Carlo method to report properties such as the polaron energy and the quasi-particle residue that have been measured in tha… ▽ More We study the repulsive Fermi polaron in a two-component, two-dimensional system of fermionic atoms inspired by the results of a recent experiment with $^{173}$Yb atoms [N. Darkwah Oppong \textit{et al.}, Phys. Rev. Lett. \textbf{122}, 193604 (2019)]. We use the diffusion Monte Carlo method to report properties such as the polaron energy and the quasi-particle residue that have been measured in that experiment. To provide insight on the quasi-particle character of the problem, we also report results for the effective mass. We show that the effective range, together with the scattering length, is needed in order to reproduce the experimental results. Using different model potentials for the interaction between the Fermi sea and the impurity, we show that it is possible to establish a regime of universality, in terms of these two parameters, that includes the whole experimental regime. This illustrates the relevance of quantum fluctuations and beyond mean-field effects to correctly describe the Fermi polaron problem. △ Less

Submitted 8 April, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Journal ref: Phys. Rev. A 103, 041302 (2021)

Estimation of the condensate fraction from the static structure factor

Authors: Yu. E. Lozovik, I. L. Kurbakov, G. E. Astrakharchik, J. Boronat

Abstract: We present an analytical method to estimate the condensate fraction $n_0/n$ in strongly correlated systems for which the zero-temperature static structure factor $S({\bf p})$ is known. The advantage of the proposed method is that it allows one to predict the long-range behavior of the one-body density matrix (i) in macroscopic and mesoscopic systems, (ii) in three- and two-dimensional geometry, (i… ▽ More We present an analytical method to estimate the condensate fraction $n_0/n$ in strongly correlated systems for which the zero-temperature static structure factor $S({\bf p})$ is known. The advantage of the proposed method is that it allows one to predict the long-range behavior of the one-body density matrix (i) in macroscopic and mesoscopic systems, (ii) in three- and two-dimensional geometry, (iii) at zero and low finite temperature, and (iv) in weakly and strongly correlated regimes. Our method is tested against exact values obtained with various quantum Monte Carlo methods in a number of strongly correlated systems showing an excellent agreement. The proposed technique is also useful in numerical simulations as it allows one to extrapolate the condensate fraction to the thermodynamic limit for particle numbers as small as tens to hundreds. Our method is especially valuable for extracting the condensate fraction from the experimentally measured static structure factor $S({\bf p})$, thus providing a new simple alternative technique for the estimation of $n_0/n$. We analyze available experimental data for $S({\bf p})$ of superfluid helium and find an excellent agreement with the experimental value of $n_0/n$. △ Less

Submitted 7 April, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 11 pages, 10 figures, published version

Journal ref: Phys. Rev. B 103, 094511 (2021)

Towards a QMC-based density functional including finite-range effects: excitation modes of a $^{39}$K quantum droplet

Authors: V. Cikojević, L. Vranješ Markić, M. Pi, M. Barranco, J. Boronat

Abstract: Some discrepancies between experimental results on quantum droplets made of a mixture of $^{39}$K atoms in different hyperfine states and their analysis within extended Gross-Pitaevskii theory (which incorporates beyond mean-field corrections) have been recently solved by introducing finite-range effects into the theory. Here, we study the influence of these effects on the monopole and quadrupole… ▽ More Some discrepancies between experimental results on quantum droplets made of a mixture of $^{39}$K atoms in different hyperfine states and their analysis within extended Gross-Pitaevskii theory (which incorporates beyond mean-field corrections) have been recently solved by introducing finite-range effects into the theory. Here, we study the influence of these effects on the monopole and quadrupole excitation spectrum of extremely dilute quantum droplets using a density functional built from first-principles quantum Monte Carlo calculations, which can be easily introduced in the existing Gross-Pitaevskii numerical solvers. Our results show differences of up to $20\%$ with those obtained within the extended Gross-Pitaevskii theory, likely providing another way to observe finite-range effects in mixed quantum droplets by measuring their lowest excitation frequencies. △ Less

Submitted 14 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.

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)

Superfluid and supersolid phases of 4He on the second layer of graphite

Authors: M. C. Gordillo, J. Boronat

Abstract: We revisited the phase diagram of the second layer of 4He on top of graphite using quantum Monte Carlo methods. Our aim was to explore the existence of the novel phases suggested recently in experimental works, and determine their properties and stability limits. We found evidence of a superfluid quantum phase with hexatic correlations, induced by the corrugation of the first Helium layer, and a q… ▽ More We revisited the phase diagram of the second layer of 4He on top of graphite using quantum Monte Carlo methods. Our aim was to explore the existence of the novel phases suggested recently in experimental works, and determine their properties and stability limits. We found evidence of a superfluid quantum phase with hexatic correlations, induced by the corrugation of the first Helium layer, and a quasi-two-dimensional supersolid corresponding to a 7/12 registered phase. The 4/7 commensurate solid was found to be unstable, while the triangular incommensurate crystals, stable at large densities, were normal. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. Let. 124, 205301 (2020)

Finite-range effects in ultradilute quantum drops

Authors: Viktor Cikojević, Leandra Vranješ Markić, Jordi Boronat

Abstract: In the first experimental realization of dilute Bose-Bose liquid drops using two hyperfine states of $^{39}$K some discrepancies between theory and experiment were observed. The standard analysis of the data using the Lee-Huang-Yang beyond mean-field theory predicted critical numbers which were significantly off the experimental measurements. Also, the radial size of the drops in the experiment pr… ▽ More In the first experimental realization of dilute Bose-Bose liquid drops using two hyperfine states of $^{39}$K some discrepancies between theory and experiment were observed. The standard analysis of the data using the Lee-Huang-Yang beyond mean-field theory predicted critical numbers which were significantly off the experimental measurements. Also, the radial size of the drops in the experiment proved to be larger than expected from this theory. Using a new functional, which is based on quantum Monte Carlo results of the bulk phase incorporating finite-range effects, we can explain the origin of the discrepancies in the critical number. This result proves the necessity of including finite-range corrections to deal with the observed properties in this setup. The controversy on the radial size is reasoned in terms of the departure from the optimal concentration ratio between the two species of the mixture. △ Less

Submitted 24 January, 2020; originally announced January 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)

Few-body bound states of two-dimensional bosons

Authors: G. Guijarro, G. E. Astrakharchik, J. Boronat, B. Bazak, D. S. Petrov

Abstract: We study clusters of the type A$_N$B$_M$ with $N\leq M\leq 3$ in a two-dimensional mixture of A and B bosons, with attractive AB and equally repulsive AA and BB interactions. In order to check universal aspects of the problem, we choose two very different models: dipolar bosons in a bilayer geometry and particles interacting via separable Gaussian potentials. We find that all the considered cluste… ▽ More We study clusters of the type A$_N$B$_M$ with $N\leq M\leq 3$ in a two-dimensional mixture of A and B bosons, with attractive AB and equally repulsive AA and BB interactions. In order to check universal aspects of the problem, we choose two very different models: dipolar bosons in a bilayer geometry and particles interacting via separable Gaussian potentials. We find that all the considered clusters are bound and that their energies are universal functions of the scattering lengths $a_{AB}$ and $a_{AA}=a_{BB}$, for sufficiently large attraction-to-repulsion ratios $a_{AB}/a_{BB}$. When $a_{AB}/a_{BB}$ decreases below $\approx 10$, the dimer-dimer interaction changes from attractive to repulsive and the population-balanced AABB and AAABBB clusters break into AB dimers. Calculating the AAABBB hexamer energy just below this threshold, we find an effective three-dimer repulsion which may have important implications for the many-body problem, particularly for observing liquid and supersolid states of dipolar dimers in the bilayer geometry. The population-imbalanced ABB trimer, ABBB tetramer, and AABBB pentamer remain bound beyond the dimer-dimer threshold. In the dipolar model, they break up at $a_{AB}\approx 2 a_{BB}$ where the atom-dimer interaction switches to repulsion. △ Less

Submitted 7 April, 2020; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: 3 figures

Journal ref: Phys. Rev. A 101, 041602 (2020)

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)

Layering transitions in superfluid helium adsorbed on a carbon nanotube mechanical resonator

Authors: Adrien Noury, Jorge Vergara-Cruz, Pascal Morfin, Bernard Plaçais, Maria Carmen Gordillo, Jordi Boronat, Sébastien Balibar, Adrian Bachtold

Abstract: Helium is recognized as a model system for the study of phase transitions. Of particular interest is the superfluid phase in two dimensions. We report measurements on superfluid helium films adsorbed on the surface of a suspended carbon nanotube. We measure the mechanical vibrations of the nanotube to probe the adsorbed helium film. We demonstrate the formation of helium layers up to five atoms th… ▽ More Helium is recognized as a model system for the study of phase transitions. Of particular interest is the superfluid phase in two dimensions. We report measurements on superfluid helium films adsorbed on the surface of a suspended carbon nanotube. We measure the mechanical vibrations of the nanotube to probe the adsorbed helium film. We demonstrate the formation of helium layers up to five atoms thickness. Upon increasing the vapour pressure, we observe layer-by-layer growth with discontinuities in both the number of adsorbed atoms and the speed of sound in the adsorbed film. These hitherto unobserved discontinuities point to a series of first-order layering transitions. Our results show that helium multilayers adsorbed on a nanotube are of unprecedented quality compared to previous works. They pave the way to new studies of quantized superfluid vortex dynamics on cylindrical surfaces, of the Berezinskii-Kosterlitz-Thouless phase transition in this new geometry, perhaps also to supersolidity in crystalline single layers as predicted in quantum Monte Carlo calculations. △ Less

Submitted 29 April, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 10 pages, 3 figures

Journal ref: Phys. Rev. Lett. 122, 165301 (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)

Universality in ultradilute liquid Bose-Bose mixtures

Authors: V. Cikojević, L. Vranješ Markić, G. E. Astrakharchik, J. Boronat

Abstract: We have studied dilute Bose-Bose mixtures of atoms with attractive interspecies and repulsive intraspecies interactions using quantum Monte Carlo methods at $T=0$. Using a number of models for interactions, we determine the range of validity of the universal equation of state of the symmetric liquid mixture as a function of two parameters: the $s$-wave scattering length and the effective range of… ▽ More We have studied dilute Bose-Bose mixtures of atoms with attractive interspecies and repulsive intraspecies interactions using quantum Monte Carlo methods at $T=0$. Using a number of models for interactions, we determine the range of validity of the universal equation of state of the symmetric liquid mixture as a function of two parameters: the $s$-wave scattering length and the effective range of the interaction potential. It is shown that the Lee-Huang-Yang correction is sufficient only for extremely dilute liquids with the additional restriction that the range of the potential is small enough. Based on the quantum Monte Carlo equation of state we develop a new density functional which goes beyond the Lee-Huang-Yang term and use it together with local density approximation to determine density profiles of realistic self-bound drops. △ Less

Submitted 11 November, 2018; originally announced November 2018.

Comments: "A related code, authored by Viktor Cikojevic, can be found at https://github.com/viktorcikojevic/Square-well-range-R-a12-reff-/tree/master"

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

Harmonically trapped Bose-Bose mixtures: a quantum Monte Carlo study

Authors: V. Cikojevic, L. Vranjes Markic, J. Boronat

Abstract: We study a harmonically confined Bose-Bose mixture using quantum Monte Carlo methods. Our results for the density profiles are systematically compared with mean-field predictions derived through the Gross-Pitaevskii equation in the same conditions. The phase space as a function of the interaction strengths and the relation between masses is quite rich. The miscibility criterion for the homogeneous… ▽ More We study a harmonically confined Bose-Bose mixture using quantum Monte Carlo methods. Our results for the density profiles are systematically compared with mean-field predictions derived through the Gross-Pitaevskii equation in the same conditions. The phase space as a function of the interaction strengths and the relation between masses is quite rich. The miscibility criterion for the homogeneous system applies rather well to the system, with some discrepancies close to the critical line for separation. We observe significant differences between the mean-field results and the Monte Carlo ones, that magnify when the asymmetry between masses increases. In the analyzed interaction regime, we observe universality of our results which extend beyond the applicability regime for the Gross-Pitaevskii equation. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: 9 pages

Journal ref: New J. Phys. 20, 085002 (2018)