BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. Boronat

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

ResearcherUnclaimed source record
cond-mat.quant-gascond-mat.othercond-mat.mtrl-scicond-mat.mes-hallcond-mat.stat-mechnucl-thphysics.atm-clusphysics.atom-ph

Catalog footprint

What is connected

54works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

54 published item(s)

preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ultradilute quantum liquid of dipolar atoms in a bilayer

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Few-body bound states of two-dimensional bosons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

3He on preplated graphite

By using the diffusion Monte Carlo method, we obtained the full phase diagram of $^3$He on top of graphite preplated with a solid layer of $^4$He. All the $^4$He atoms of the substrate were explicitly considered and allowed to move during the simulation. We found that the ground state is a liquid of density 0.007 $\pm$ 0.001 Å$^{-2}$, in good agreement with available experimental data. This is significantly different from the case of $^3$He on clean graphite, in which both theory and experiment agree on the existence of a gas-liquid transition at low densities. Upon an increase in $^3$He density, we predict a first-order phase transition between a dense liquid and a registered 7/12 phase, the 4/7 phase being found metastable in our calculations. At larger second-layer densities, a final transition is produced to an incommensurate triangular phase.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Droplets of trapped quantum dipolar bosons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamic structure factor of liquid 4He across the normal-superfluid transition

We have carried out a microscopic study of the dynamic structure factor of liquid $^4$He across the normal-superfluid transition temperature using the path integral Monte Carlo method. The ill-posed problem of the inverse Laplace transform, from the imaginary-time intermediate scattering function to the dynamic response, is tackled by stochastic optimization. Our results show a quasi-particle peak and a small and broad multiphonon contribution. In spite of the lack of strength in the collective peaks, we clearly identify the rapid dropping of the roton peak amplitude when crossing the transition temperature $T_λ$. Other properties such as the static structure factor, static response, and one-phonon contribution to the response are also calculated at different temperatures. The changes of the phonon-roton spectrum with the temperature are also studied. An overall agreement with available experimental data is achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gapped spectrum in pair-superfluid bosons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Liquid and solid phases of 3He on graphite

Recent heat-capacity experiments show quite unambiguously the existence of a liquid $^3$He phase adsorbed on graphite. This liquid is stable at an extremely low density, possibly one of the lowest found in Nature. Previous theoretical calculations of the same system, and in strictly two dimensions, agree with the result that this liquid phase is not stable and the system is in the gas phase. We calculated the phase diagram of normal $^3$He adsorbed on graphite at $T=0$ using quantum Monte Carlo methods. Considering a fully corrugated substrate we observe that at densities lower that 0.006 Å$^{-2}$ the system is a very dilute gas, that at that density is in equilibrium with a liquid of density 0.014 Å$^{-2}$. Our prediction matches very well the recent experimental findings on the same system. On the contrary, when a flat substrate is considered, no gas-liquid coexistence is found, in agreement with previous calculations. We also report results on the different solid structures, and the corresponding phase transitions that appear at higher densities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Momentum distribution of liquid 4He across the normal-superfluid phase transition

We have carried out a study of the momentum distribution and of the spectrum of elementary excitations of liquid $^4$He across the normal-superfluid transition temperature, using the path integral Monte Carlo method. Our results for the momentum distribution in the superfluid regime show that a kink is present in the range of momenta corresponding to the roton excitation. This effect disappears when crossing the transition temperature to the normal fluid, in a behavior currently unexplained by theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

One Dimensional 1H, 2H and 3H

The ground-state properties of one-dimensional electron-spin-polarized hydrogen $^1$H, deuterium $^2$H, and tritium $^3$H are obtained by means of quantum Monte Carlo methods. The equations of state of the three isotopes are calculated for a wide range of linear densities. The pair correlation function and the static structure factor are obtained and interpreted within the framework of the Luttinger liquid theory. We report the density dependence of the Luttinger parameter and use it to identify different physical regimes: Bogoliubov Bose gas, super-Tonks-Girardeau gas, and quasi-crystal regimes for bosons; repulsive, attractive Fermi gas, and quasi-crystal regimes for fermions. We find that the tritium isotope is the one with the richest behaviour. Our results show unambiguously the relevant role of the isotope mass in the properties of this quantum system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Luttinger liquid behavior of one-dimensional 3He

The ground-state properties of one-dimensional 3He are studied using quantum Monte Carlo methods. The equation of state is calculated in a wide range of physically relevant densities and is well reproduced by a power-series fit. The Luttinger liquid theory is found to describe the long-range properties of the correlation function. The density dependence of the Luttinger parameter is explicitly found and interestingly it shows a non-monotonic behavior. Depending on the density, the static structure factor can be a smooth function of the momentum or might contain a peak of a finite or infinite height. Although no phase transitions are present in the system, we identify a number of physically different regimes, including an ideal Fermi gas, a "Bose-gas", a "super-Tonks-Girardeau" regime, and a "quasi-crystal".

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of a quantum Coulomb wire

We report the quantum phase diagram of a one-dimensional Coulomb wire obtained using the path integral Monte Carlo (PIMC) method. The exact knowledge of the nodal points of this system permits us to find the energy in an exact way, solving the sign problem which spoils fermionic calculations in higher dimensions. The results obtained allow for the determination of the stability domain, in terms of density and temperature, of the one-dimensional Wigner crystal. At low temperatures, the quantum wire reaches the quantum-degenerate regime, which is also described by the diffusion Monte Carlo method. Increasing the temperature the system transforms to a classical Boltzmann gas which we simulate using classical Monte Carlo. At large enough density, we identify a one-dimensional ideal Fermi gas which remains quantum up to higher temperatures than in two- and three-dimensional electron gases. The obtained phase diagram as well as the energetic and structural properties of this system are relevant to experiments with electrons in quantum wires and to Coulomb ions in one-dimensional confinement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distinguishability, degeneracy and correlations in three harmonically trapped bosons in one-dimension

We study a system of two bosons of one species and a third boson of a second species in a one-dimensional parabolic trap at zero temperature. We assume contact repulsive inter- and intra-species interactions. By means of an exact diagonalization method we calculate the ground and excited states for the whole range of interactions. We use discrete group theory to classify the eigenstates according to the symmetry of the interaction potential. We also propose and validate analytical ansatzs gaining physical insight over the numerically obtained wavefunctions. We show that, for both approaches, it is crucial to take into account that the distinguishability of the third atom implies the absence of any restriction over the wavefunction when interchanging this boson with any of the other two. We find that there are degeneracies in the spectra in some limiting regimes, that is, when the inter-species and/or the intra-species interactions tend to infinity. This is in contrast with the three-identical boson system, where no degeneracy occurs in these limits. We show that, when tuning both types of interactions through a protocol that keeps them equal while they are increased towards infinity, the systems's ground state resembles that of three indistinguishable bosons. Contrarily, the systems's ground state is different from that of three-identical bosons when both types of interactions are increased towards infinity through protocols that do not restrict them to be equal. We study the coherence and correlations of the system as the interactions are tuned through different protocols, which permit to built up different correlations in the system and lead to different spatial distributions of the three atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase diagram of dipolar bosons in 2D with tilted polarization

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum correlations and spatial localization in one-dimensional ultracold bosonic mixtures

We present the complete phase diagram for one-dimensional binary mixtures of bosonic ultracold atomic gases in a harmonic trap. We obtain exact results with direct numerical diagonalization for small number of atoms, which permits us to quantify quantum many-body correlations. The quantum Monte Carlo method is used to calculate energies and density profiles for larger system sizes. We study the system properties for a wide range of interaction parameters. For the extreme values of these parameters, different correlation limits can be identified, where the correlations are either weak or strong. We investigate in detail how the correlation evolve between the limits. For balanced mixtures in the number of atoms in each species, the transition between the different limits involves sophisticated changes in the one- and two-body correlations. Particularly, we quantify the entanglement between the two components by means of the von Neumann entropy. We show that the limits equally exist when the number of atoms is increased, for balanced mixtures. Also, the changes in the correlations along the transitions among these limits are qualitatively similar. We also show that, for imbalanced mixtures, the same limits with similar transitions exist. Finally, for strongly imbalanced systems, only two limits survive, i.e., a miscible limit and a phase-separated one, resembling those expected with a mean-field approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Single-particle vs. pair superfluidity in a bilayer system of dipolar bosons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality in molecular halo clusters

Ground state of weakly bound dimers and trimers with a radius extending well into the classically forbidden region is explored, with the goal to test the predicted universality of quantum halo states. The focus of the study are molecules consisting of T$\downarrow$, D$\downarrow$, $^3$He, $^4$He and alkali atoms, where interaction between particles is much better known than in the case of nuclei, which are traditional examples of quantum halos. The study of realistic systems is supplemented by model calculations in order to analyze how low-energy properties depend on the interaction potential. The use of variational and diffusion Monte Carlo methods enabled very precise calculation of both size and binding energy of the trimers. In the quantum halo regime, and for large values of scaled binding energies, all clusters follow almost the same universal line. As the scaled binding energy decreases, Borromean states separate from tango trimers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamic structure function of a cold Fermi gas at unitarity

We present a theoretical study of the dynamic structure function of a resonantly interacting two-component Fermi gas at zero temperature. Our approach is based on dynamic many-body theory able to describe excitations in strongly correlated Fermi systems. The fixed-node diffusion Monte Carlo method is used to produce the ground-state correlation functions which are used as an input for the excitation theory. Our approach reproduces recent Bragg scattering data in both the density and the spin channel. In the BCS regime, the response is close to that of the ideal Fermi gas. On the BEC side, the Bose peak associated with the formation of dimers dominates the density channel of the dynamic response. When the fraction of dimers is large our theory departs from the experimental data, mainly in the spin channel.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Monte Carlo study of spin-polarized deuterium

The ground state properties of spin-polarized deuterium (D$\downarrow$) at zero temperature are obtained by means of the diffusion Monte Carlo calculations within the fixed-node approximation. Three D$\downarrow$ species have been investigated (D$\downarrow_1$, D$\downarrow_2$, D$\downarrow_3$), corresponding respectively to one, two and three equally occupied nuclear spin states. Influence of the backflow correlations on the ground state energy of the systems is explored. The equilibrium densities for D$\downarrow_2$ and D$\downarrow_3$ liquids are obtained and compared with ones obtained in previous approximate prediction. The density and the pressure at which the gas-liquid phase transition occurs at $T$=0 is obtained for D$\downarrow_1$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Second layer of H2 and D2 adsorbed on graphene

We report diffusion Monte Carlo calculations on the phase diagrams of para-H2 and ortho-D2 adsorbed on top of a first layer of the same substances on graphene. We found that the ground state of the second layer is a triangular incommensurate solid for both isotopes. The densities for promotion to a second layer and for the onset of a two-dimensional solid on that second layer compare favorably to available experimental data in both cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spin-polarized hydrogen adsorbed on the surface of superfluid 4He

The experimental realization of a thin layer of spin-polarized hydrogen adsorbed on top of the surface of superfluid 4He provides one of the best examples of a stable nearly two-dimensional quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid 4He slab and the adsorbed hydrogen layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly two-dimensional (2D) setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as $60\%$ at the largest coverage studied.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

4He adsorbed outside a single carbon nanotube

The phase diagrams of $^4$He adsorbed on the external surfaces of single armchair carbon nanotubes with radii in the range 3.42 -- 10.85 Å\ are calculated using the diffusion Monte Carlo method. For nanotubes narrower than a (10,10) one, the ground state is an incommensurate solid similar to the one found for H$_2$ on the same substrates. For wider nanotubes, the phase with the minimum energy per particle is a liquid layer. Curved $\sqrt 3 \times \sqrt 3$ registered solids similar to the ones found on graphene and graphite were unstable for all the tubes considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A microscopic description of vacancies in solid 4He

The changes that vacancies produce in the properties of hcp solid 4He are studied by means of quantum Monte Carlo methods. Our results show that the introduction of vacancies produces significant changes in the behavior of solid 4He, even when the vacancy concentration is very small. We show that there is an onset temperature where the properties of incommensurate 4He change significantly. Below this temperature, we observe the emergence of off-diagonal long range order and a complete spatial delocalization of the vacancies. This temperature is quite close to the temperature where non-classical rotational inertia has been experimentally observed. Finally, we report results on the influence of vacancies in the elastic properties of hcp 4He at zero temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ferromagnetic transition of a two-component Fermi gas of Hard Spheres

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ground state properties and excitation spectrum of a two dimensional gas of bosonic dipoles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stability of resonantly interacting heavy-light Fermi mixtures

We investigate a two-component mixture of resonantly interacting Fermi gases as a function of the ratio κof the heavy to the light mass of the two species. The diffusion Monte Carlo method is used to calculate the ground-state energy and the pair correlation function starting from two different guiding wave functions, which describe respectively the superfluid and the normal state of the gas. Results show that the mixture is stable and superfluid for mass ratios smaller than the critical value κ_c = 13 +/- 1. For larger values of κsimulations utilizing the wave function of the normal state are unstable towards cluster formation. The relevant cluster states driving the instability appear to be formed by one light particle and two or more heavy particles within distances on the order of the range of the interatomic potential. The small overlap between the wave function of the trimer bound state and the guiding wave function used to describe the superfluid state produces the unphysical stability of the superfluid gas above κ_c.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Superfluidity of metastable bulk glass para-hydrogen at low temperature

Molecular para-hydrogen has been proposed theoretically as a possible candidate for superfluidity, but the eventual superfluid transition is hindered by its crystallization. In this work, we study a metastable non crystalline phase of bulk p-H2 by means of the Path Integral Monte Carlo method in order to investigate at which temperature this system can support superfluidity. By choosing accurately the initial configuration and using a non commensurate simulation box, we have been able to frustrate the formation of the crystal in the simulated system and to calculate the temperature dependence of the one-body density matrix and of the superfluid fraction. We observe a transition to a superfluid phase at temperatures around 1 K. The limit of zero temperature is also studied using the diffusion Monte Carlo method. Results for the energy, condensate fraction, and structure of the metastable liquid phase at T=0 are reported and compared with the ones obtained for the stable solid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Zero-temperature phase diagram of the second layer of $^{\bf 4}$He adsorbed on graphene

The phase diagram at zero temperature of $^4$He adsorbed on an helium incommensurate triangular solid on top of a single graphene sheet has been obtained using the diffusion Monte Carlo method. We have found that, in accordance with previous experimental and simulation results for graphite, the ground state of $^4$He on this setup is a liquid that, upon compression, transforms into a triangular solid. To define the stability limits of both liquid and solid phases, we considered not only the adsorption energies of the atoms located on the second layer but the average energy of the atoms in both layers. Our results show that the lower density limit for a stable liquid in the second layer is 0.163 $\pm$ 0.005 Å$^{-2}$ and that the lower limit for the existence of an incommensurate solid on the second layer is 0.186 $\pm$ 0.003 Å$^{-2}$. Both values are in overall agreement with the results of torsional oscillator experiments and heat capacity measurements on graphite. The 4/7 and 7/12 registered solids are found to be metastable with respect to triangular incommensurate arrangements of the same density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Zero-temperature phase diagram of Yukawa bosons

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bose-Einstein Condensation in liquid $^4$He near the liquid-solid transition line

We present precision neutron scattering measurements of the Bose-Einstein condensate fraction, n0(T), and the atomic momentum distribution, n\star(k), of liquid 4He at pressure p =24 bar. Both the temperature dependence of n0(T) and of the width of n\star(k) are determined. The n0(T) can be represented by n0(T) = n0(0)[1-(T/Tλ)γ] with a small n0(0) = 2.80\pm0.20% and large γ = 13\pm2 for T < Tλ indicating strong interaction. The onset of BEC is accompanied by a significant narrowing of the n\star(k). The narrowing accounts for 65% of the drop in kinetic energy below Tλ and reveals an important coupling between BEC and k > 0 states. The experimental results are well reproduced by Path Integral Monte Carlo calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ewald method for polytropic potentials in arbitrary dimensionality

The Ewald summation technique is generalised to power-law 1/|r|^k potentials in three-, two- and one-dimensional geometries with explicit formulae for all the components of the sums. The cases of short-range, long-range and "marginal" interactions are treated separately. The jellium model, as a particular case of a charge-neutral system, is discussed and the explicit forms of the Ewald sums for such system are presented. A generalised form of the Ewald sums for a noncubic (nonsquare) simulation cell for three- (two-) dimensional geometry is obtained and its possible field of application is discussed. A procedure for the optimisation of the involved parameters in actual simulations is developed and an example of its application is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ground-state properties and superfluidity of two- and quasi two-dimensional solid 4He

In a recent study we have reported a new type of trial wave function symmetric under the exchange of particles and which is able to describe a supersolid phase. In this work, we use the diffusion Monte Carlo method and this model wave function to study the properties of solid 4He in two- and quasi two-dimensional geometries. In the purely two-dimensional case, we obtain results for the total ground-state energy and freezing and melting densities which are in good agreement with previous exact Monte Carlo calculations performed with a slightly different interatomic potential model. We calculate the value of the zero-temperature superfluid fraction ρ_{s} / ρof 2D solid 4He and find that it is negligible in all the considered cases, similarly to what is obtained in the perfect (free of defects) three-dimensional crystal using the same computational approach. Interestingly, by allowing the atoms to move locally in the perpendicular direction to the plane where they are confined to zero-point oscillations (quasi two-dimensional crystal) we observe the emergence of a finite superfluid density that coexists with the periodicity of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microscopic description of anisotropic low-density dipolar Bose gases in two dimensions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram of Rydberg atoms with repulsive van der Waals interaction

We report a quantum Monte Carlo calculation of the phase diagram of bosons interacting with a repulsive inverse sixth power pair potential, a model for assemblies of Rydberg atoms in the local van der Waals blockade regime. The model can be parametrized in terms of just two parameters, the reduced density and temperature. Solidification happens to the fcc phase. At zero temperature the transition density is found with the diffusion Monte Carlo method at density $ρ= 3.9 (\hbar^2/m C_6)^{3/4} $, where $C_6$ is the strength of the interaction. The solidification curve at non-zero temperature is studied with the path integral Monte Carlo approach and is compared with transitions in corresponding harmonic and classical crystals. Relaxation mechanisms are considered in relation to present experiments, especially pertaining to hopping of the Rydberg excitation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase transitions of H2 adsorbed on the surface of single carbon nanotubes

By means of Diffusion Monte Carlo calculations, we obtained the complete phase diagrams of H$_2$ adsorbed on the outer surface of isolated armchair carbon nanotubes of radii ranging from 3.42 to 10.85 Å. We only considered density ranges corresponding to the filling of the first adsorption layer in these curved structures. In all cases, the zero-temperature ground state was found to be an incommensurate solid, except in the widest tube, in which the structure with lowest energy is an analogous of the $\sqrt{3} \times \sqrt{3}$ phase found in planar substrates. Those incommensurate solids result form the arrangement of the hydrogen molecules in circumferences whose plane is perpendicular to the main axis of the carbon nanotube. For each tube, there is only one of such phases stable in the density range considered, except in the case of the (5,5) and (6,6) tubes, in which two of these incommensurate solids are separated by novel first order phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Solidification of small para-H2 clusters at zero temperature

We have determined the ground-state energies of para-H$_2$ clusters at zero temperature using the diffusion Monte Carlo method. The liquid or solid character of each cluster is investigated by restricting the phase through the use of proper importance sampling. Our results show inhomogeneous crystallization of clusters, with alternating behavior between liquid and solid phases up to N=55. From there on, all clusters are solid. The ground-state energies in the range N=13--75 are established and the stable phase of each cluster is determined. In spite of the small differences observed between the energy of liquid and solid clusters, the corresponding density profiles are significantly different, feature that can help to solve ambiguities in the determination of the specific phase of H$_2$ clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Supersolidity in quantum films adsorbed on graphene and graphite

Using quantum Monte Carlo we have studied the superfluid density of the first layer of $^4$He and H$_2$ adsorbed on graphene and graphite. Our main focus has been on the equilibrium ground state of the system, which corresponds to a registered $\sqrt3 \times \sqrt3$ phase. The perfect solid phase of H$_2$ shows no superfluid signal whereas $^4$He has a finite but small superfluid fraction (0.67%). The introduction of vacancies in the crystal makes the superfluidity increase, showing values as large as 14% in $^4$He without destroying the spatial solid order.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The elastic constants of solid 4He under pressure: a diffusion Monte Carlo study

We study the elasticity of perfect 4He at zero-temperature using the diffusion Monte Carlo method and a realistic semi-empirical pairwise potential to describe the He-He interactions. Specifically, we calculate the value of the elastic constants of hcp helium C_{ij} as a function of pressure up to 110 bar. It is found that the pressure dependence of all five non-zero C_{ij} is linear and we provide accurate parametrization of each of them. Our elastic constants results are compared to previous variational calculations and low-temperature measurements and in general notably good agreement is found among them. Furthermore, we report T = 0 results for the Gruneisen parameters, sound velocities and Debye temperature of hcp 4He. This work represents the first of a series of computational studies aimed at thoroughly characterizing the response of solid helium to external stress-strain.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground state of small mixed helium and spin-polarized tritium clusters: a quantum Monte Carlo stud

We report results for the ground-state energy and structural properties of small $^4$He-T$\downarrow$ clusters consisting of up to 4 T$\downarrow$ and 8 $^4$He atoms. These results have been obtained using very well-known $^4$He-$^4$He and T$\downarrow$-T$\downarrow$ interaction potentials and several models for the $^4$He-T$\downarrow$ interatomic potential. All the calculations have been performed with variational and diffusion Monte Carlo methods. It takes at least three atoms to form a mixed bound state. In particular, for small clusters the binding energies are significantly affected by the precise form of the $^4$He-T$\downarrow$ interatomic potential but the stability limits remain unchanged. The only exception is the $^4$He$_2$T$\downarrow$ trimer whose stability in the case of the weakest $^4$He-T$\downarrow$ interaction potential is uncertain, while it seems stable for other potentials. The mixed trimer $^4$He(T$\downarrow$)$_2$, a candidate for Borromean state, is not bound. All other studied clusters are stable. Some of the weakest bound clusters can be classified as quantum halos, as a consequence of having high probability of being in a classically forbidden region.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagram of H2 adsorbed on graphene

The phase diagram of the first layer of H$_2$ adsorbed on top of a single graphene sheet has been calculated by means of a series of diffusion Monte Carlo (DMC) simulations. We have found that, as in the case of $^4$He, the ground state of molecular hydrogen is a $\sqrt3 \times \sqrt3$ commensurate structure, followed, upon a pressure increase, by an incommensurate triangular solid. A striped phase of intermediate density was also considered, and found lying on top of the equilibrium curve separating both commensurate and incommensurate solids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Properties of Vacancy Formation in He-4 hcp Crystals at Zero Temperature and Fixed Pressure

Equation of state of He-4 hcp crystals with vacancies is determined at zero temperature using the diffusion Monte Carlo technique, an exact ground state zero-temperature method. This allows us to extract the formation enthalpy and isobaric formation energy of a single vacancy in otherwise perfect helium solid. Results were obtained for pressures up to 160 bar. The isobaric formation energy is found to reach a minimum near 57 bar where it is equal to $10.5\pm 1.2$ K. At the same pressure, the vacancy formation volume exhibits a maximum and reaches the volume of the unit cell. This pressure coincides with the pressure interval over which a peak in the supersolid fraction of He-4 was observed in a recent experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo study of large spin-polarized tritium clusters

This work expands recent investigations in the field of spin-polarized tritium (T$\downarrow$) clusters . We report the results for the ground state energy and structural properties of large T$\downarrow$ cl usters consisting of up to 320 atoms. All calculations have been performed with variational and diffusi on Monte Carlo methods, using an accurate {\it ab initio} interatomic potential. Our results for $N \le q 40$ are in good agreement with results obtained by other groups. Using a liquid-drop expression for t he energy per particle, we estimate the liquid equilibrium density, which is in good agreement with our recently obtained results for bulk T$\downarrow$. In addition, the calculations of the energy for larg e clusters have allowed for an estimation of the surface tension. From the mean-square radius of the dr op, determined using unbiased estimators, we determine the dependence of the radii on the size of the c luster and extract the unit radius of the T$\downarrow$ liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quasiequilibrium supersolid phase of a two-dimensional dipolar crystal

We have studied the possible existence of a supersolid phase of a two-dimensional dipolar crystal using quantum Monte Carlo methods at zero temperature. Our results show that the commensurate solid is not a supersolid in the thermodynamic limit. The presence of vacancies or interstitials turns the solid into a supersolid phase even when a tiny fraction of them are present in a macroscopic system. The effective interaction between vacancies is repulsive making a quasiequilibrium dipolar supersolid possible.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High order Chin actions in path integral Monte Carlo

High order actions proposed by Chin have been used for the first time in path integral Monte Carlo simulations. Contrarily to the Takahashi-Imada action, which is accurate to fourth order only for the trace, the Chin action is fully fourth order, with the additional advantage that the leading fourth and sixth order error coefficients are finely tunable. By optimizing two free parameters entering in the new action we show that the time step error dependence achieved is best fitted with a sixth order law. The computational effort per bead is increased but the total number of beads is greatly reduced, and the efficiency improvement with respect to the primitive approximation is approximately a factor of ten. The Chin action is tested in a one-dimensional harmonic oscillator, a H$_2$ drop, and bulk liquid $^4$He. In all cases a sixth-order law is obtained with values of the number of beads that compare well with the pair action approximation in the stringent test of superfluid $^4$He.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High-order Time Expansion Path Integral Ground State

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Low-dimensional weakly interacting Bose gases: non-universal equations of state

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum Fluids in Nanotubes: a Quantum Monte Carlo Approach

We review quantum Monte Carlo results on energetic and structure properties of quantum fluids adsorbed in a bundle of carbon nanotubes. Using realistic interatomic interactions the different adsorption sites that a bundle offer are accurately studied and compared in some cases with strictly one-dimensional geometries. The study is performed quite extensively for $^4$He and restricted to the inner part of a single nanotube for H$_2$ and D$_2$. From a theoretical point of view, nanotubes open the real possibility of a quasi-one-dimensional confinement where to study quantum fluids in extremely reduced dimensionality. The results obtained show that in the narrowest configurations the system is nearly one-dimensional reinforcing the interest on the physics of one-dimensional quantum fluids. Experimental confirmation of the theoretical results obtained is still not in a satisfactory situation due to the difficulties on extracting from the data the dominant adsorption sites.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum Monte Carlo simulation of spin-polarized tritium

The ground-state properties of spin-polarized tritium T$\downarrow$ at zero temperature are obtained by means of diffusion Monte Carlo calculations. Using an accurate {\em ab initio} T$\downarrow$-T$\downarrow$ interatomic potential we have studied its liquid phase, from the spinodal point until densities above its freezing point. The equilibrium density of the liquid is significantly higher and the equilibrium energy of $-3.664(6)$ K significantly lower than in previous approximate descriptions. The solid phase has also been studied for three lattices up to high pressures, and we find that hcp lattice is slightly preferred. The liquid-solid phase transition has been determined using the double-tangent Maxwell construction; at zero temperature, bulk tritium freezes at a pressure of $P=9(1)$ bar.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Zero-temperature equation of state of solid 4He at low and high pressures

We study the zero-temperature equation of state (EOS) of solid 4He in the hexagonal closed packet (hcp) phase over the 0-57 GPa pressure range by means of the Diffusion Monte Carlo (DMC) method and the semi-empirical Aziz pair potential HFD-B(HE). In the low pressure regime (P ~ 0-1 GPa) we assess excellent agreement with experiments and we give an accurate description of the atomic kinetic energy, Lindemann ratio and Debye temperature over a wide range of molar volumes (22-6 cm^{3}/mol). However, on moving to higher pressures our calculated P-V curve presents an increasingly steeper slope which ultimately provides differences within ~40 % with respect to measurements. In order to account for many-body interactions arising in the crystal with compression which are not reproduced by our model, we perform additional electronic density-functional theory (DFT) calculations for correcting the computed DMC energies in a perturbative way. We explore both generalized gradient and local density approximations (GGA and LDA, respectively) for the electronic exchange-correlation potential. By proceeding in this manner, we show that discrepancies with respect to high pressure data are reduced to 5-10 % with few computational extra cost. Further comparison between our calculated EOSs and ab initio curves deduced for the perfect crystal and corrected for the zero-point motion of the atoms enforces the reliability of our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Equation of state of an interacting Bose gas at finite temperature: a Path Integral Monte Carlo study

By using exact Path Integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Momentum distributions in ^3He-^4He liquid mixtures

We present variational calculations of the one-body density matrices and momentum distributions for ^3He-^4He mixtures in the zero temperature limit, in the framework of the correlated basis functions theory. The ground-state wave function contains two- and three-body correlations and the matrix elements are computed by (Fermi)Hypernetted Chain techniques. The dependence on the ^3He concentration (x_3) of the ^4He condensate fraction $(n_0^{(4)})$ and of the ^3He pole strength (Z_F) is studied along the P=0 isobar. At low ^3He concentration, the computed ^4He condensate fraction is not significantly affected by the ^3He statistics. Despite of the low x_3 values, Z_F is found to be quite smaller than that of the corresponding pure ^3He because of the strong ^3He-^4He correlations and of the overall, large total density ρ. A small increase of $n_0^{(4)}$ along x_3 is found, which is mainly due to the decrease of ρrespect to the pure ^4He phase.

0 citations0 reviewsNo code linkedOpen access