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

G. E. Astrakharchik

G. E. Astrakharchik 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.otherphysics.atom-phcond-mat.mes-hallcond-mat.stat-mechcond-mat.str-elmath.DSquant-ph

Catalog footprint

What is connected

41works
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

41 published item(s)

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

The origin of the period-$2T/7$ quasi-breathing in disk-shaped Gross-Pitaevskii breathers

We address the origins of the quasi-periodic breathing observed in [Phys. Rev.\ X vol. 9, 021035 (2019)] in disk-shaped harmonically trapped two-dimensional Bose condensates, where the quasi-period $T_{\text{quasi-breathing}}\sim$~$2T/7$ and $T$ is the period of the harmonic trap. We show that, due to an unexplained coincidence, the first instance of the collapse of the hydrodynamic description, at $t^{*} = \arctan(\sqrt{2})/(2π) T \approx T/7$, emerges as a `skillful impostor' of the quasi-breathing half-period $T_{\text{quasi-breathing}}/2$. At the time $t^{*}$, the velocity field almost vanishes, supporting the requisite time-reversal invariance. We find that this phenomenon persists for scale-invariant gases in all spatial dimensions, being exact in one dimension and, likely, approximate in all others. In $d$ dimensions, the quasi-breathing half-period assumes the form $T_{\text{quasi-breathing}}/2 \equiv t^{*} = \arctan(\sqrt{d})/(2π) T$. Remaining unresolved is the origin of the period-$2T$ breathing, reported in the same experiment.

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

Strong coupling Bose polarons in a two-dimensional gas

We study the properties of Bose polarons in two dimensions using quantum Monte Carlo techniques. Results for the binding energy, the effective mass, and the quasiparticle residue are reported for a typical strength of interactions in the gas and for a wide range of impurity-gas coupling strengths. A lower and an upper branch of the quasiparticle exist. The lower branch corresponds to an attractive polaron and spans from the regime of weak coupling where the impurity acts as a small density perturbation of the surrounding medium to deep bound states which involve many particles from the bath and extend as far as the healing length. The upper branch corresponds to an excited state where due to repulsion a low-density bubble forms around the impurity but might be unstable against decay into many-body bound states. Interaction effects strongly affect the quasiparticle properties of the polaron. In particular, in the strongly correlated regime, the impurity features a vanishing quasiparticle residue, signaling the transition from an almost free quasiparticle to a bound state involving many atoms from the bath.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The dynamics of digits: Calculating pi with Galperin's billiards

In Galperin billiards, two balls colliding with a hard wall form an analog calculator for the digits of the number $π$. This classical, one-dimensional three-body system (counting the hard wall) calculates the digits of $π$ in a base determined by the ratio of the masses of the two particles. This base can be any integer, but it can also be an irrational number, or even the base can be $π$ itself. This article reviews previous results for Galperin billiards and then pushes these results farther. We provide a complete explicit solution for the balls' positions and velocities as a function of the collision number and time. We demonstrate that Galperin billiard can be mapped onto a two-particle Calogero-type model. We identify a second dynamical invariant for any mass ratio that provides integrability for the system, and for a sequence of specific mass ratios we identify a third dynamical invariant that establishes superintegrability. Integrability allows us to derive some new exact results for trajectories, and we apply these solutions to analyze the systematic errors that occur in calculating the digits of $π$ with Galperin billiards, including curious cases with irrational number bases.

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

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
preprint2016arXiv

One-dimensional multicomponent Fermi gas in a trap: quantum Monte Carlo study

One-dimensional world is very unusual as there is an interplay between quantum statistics and geometry, and a strong short-range repulsion between atoms mimics Fermi exclusion principle, fermionizing the system. Instead, a system with a large number of components with a single atom in each, on the opposite acquires many bosonic properties. We study the ground-state properties a multi-component Fermi gas trapped in a harmonic trap by fixed-node diffusion Monte Carlo method. We investigate how the energetic properties (energy, contact) and correlation functions (density profile and momentum distribution) evolve as the number of components is changed. It is shown that the system fermionizes in the limit of strong interactions. Analytical expression are derived in the limit of weak interactions within the local density approximation for arbitrary number of components and for one plus one particle using an exact solution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ultradilute low-dimensional liquids

We calculate the energy of one- and two-dimensional weakly interacting Bose-Bose mixtures analytically in the Bogoliubov approximation and by using the diffusion Monte Carlo technique. We show that in the case of attractive inter- and repulsive intraspecies interactions the energy per particle has a minimum at a finite density corresponding to a liquid state. We derive the Gross-Pitaevskii equation to describe droplets of such liquids and solve it analytically in the one-dimensional case.

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
preprint2015arXiv

Quantum versus mean-field collapse in a many-body system

The recent analysis, based on the mean-field approximation (MFA), has predicted that the critical quantum collapse of the bosonic wave function, pulled to the center by the inverse-square potential in the three-dimensional space, is suppressed by the repulsive cubic nonlinearity in the bosonic gas, the collapsing ground state being replaced by a regular one. We demonstrate that a similar stabilization acts in a quantum many-body system, beyond the MFA. While the collapse remains possible, repulsive two-particle interactions give rise to a metastable gaseous state, which is separated by a potential barrier from the collapsing regime. The stability of this state improves with the increase of the number of particles. The results are produced by calculations of the variational energy, with the help of the Monte Carlo method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reentrant behavior of the breathing-mode-oscillation frequency in a one-dimensional Bose gas

We calculate the breathing mode frequency $ω$ in a one-dimensional Bose gas confined to a harmonic trap of frequency $ω_z$. We predict Exciting temporal oscillations of the density distribution is a high-precision method for probing ultracold trapped atomic gases. Interaction effects in their many-body dynamics are particularly puzzling and counter-intuitive in one spatial dimension (1D) due to enhanced quantum correlations. We consider 1D quantum Bose gas in a parabolic trap at zero temperature and explain, analytically and numerically, how oscillation frequency depends on the number of particles, their repulsion and the trap strength. We identify the frequency with the energy difference between the ground state and a particular excited state. This way we avoided resolving the dynamical evolution of the system, simplifying the problem immensely. We find an excellent quantitative agreement of our results with the data from the Innsbruck experiment [Science 325, 1224 (2009)].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Breakdown of the scale invariance in the vicinity of Tonks-Girardeau gas

In this article, we consider the monopole excitations of the harmonically trapped Bose gas in the vicinity of the Tonks-Girardeau limit. Using Girardeau's Fermi-Bose duality and subsequently an effective fermion-fermion odd-wave interaction, we obtain the dominant correction to the scale-invariance-protected value of the excitation frequency, for microscopically small excitation amplitudes. We produce a series of diffusion Monte Carlo results that confirm our analytic prediction for three particles. And less expectedly, our result stands in excellent agreement with the result of a hydrodynamic simulation of the microscopically large but macroscopically small excitations.

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 quasi-two-dimensional bosons in laser speckle potential

We have studied the phase diagram of a quasi-two-dimensional interacting Bose gas at zero temperature in the presence of random potential created by laser speckles. The superfluid fraction and the fraction of particles with zero momentum are obtained within the mean-field Gross-Pitaevskii theory and in diffusion Monte Carlo simulations. We find a transition from the superfluid to the insulating state, when the strength of the disorder grows. Estimations of the critical parameters are compared with the predictions of the percolation theory in the Thomas-Fermi approximation. Analytical expressions for the zero-momentum fraction and the superfluid fraction are derived in the limit of weak disorder and weak interactions within the framework of the Bogoliubov theory. Limits of validity of various approximations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Properties of a Bose Gas in the Presence of Disorder (Laurea thesis)

The phenomenon of Bose-Einstein condensation and superfluidity in a Bose gas with disorder is investigated. Diffusion Monte Carlo (DMC) method is used to calculate superfluid and condensate fraction of the system as a function of density and strength of disorder at zero temperature. The algorithm and implementation of the Diffusion Monte Carlo method is explained in details. Bogoliubov theory is developed for the analytical description of the problem. Ground state energy, superfluid fraction and condensate fraction are calculated. It is shown that same results for the superfluid fraction can be obtained in a perturbative manner from Gross-Pitaevskii equation. Ground state energy, obtained from DMC calculations, is compared to predictions of Bogoliubov theory, which are found to be valid in the regime, when the strength of disorder is small. It is shown that "unusual" situation, when the superfluid fraction is smaller than the condensate fraction, can be realized in this system.

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

Quantum Monte Carlo study of ultracold gases (PhD thesis)

This Dissertation presents results of a thorough study of ultracold bosonic and fermionic gases in three-dimensional and quasi-one-dimensional systems. Although the analyses are carried out within various theoretical frameworks (Gross-Pitaevskii, Bethe ansatz, local density approximation, etc.) the main tool of the study is the Quantum Monte Carlo method in different modifications (variational Monte Carlo, diffusion Monte Carlo, fixed-node Monte Carlo methods). We benchmark our Monte Carlo calculations by recovering known analytical results (perturbative theories in dilute limits, exactly solvable models, etc.) and extend calculations to regimes, where the results are so far unknown. In particular we calculate the equation of state and correlation functions for gases in various geometries and with various interatomic interactions.

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
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

Lieb's soliton-like excitations in harmonic trap

We study the solitonic Lieb II branch of excitations in one-dimensional Bose-gas in homogeneous and trapped geometry. Using Bethe-ansatz Lieb's equations we calculate the "effective number of atoms" and the "effective mass" of the excitation. The equations of motion of the excitation are defined by the ratio of these quantities. The frequency of oscillations of the excitation in a harmonic trap is calculated. It changes continuously from its "soliton-like" value ω_h/\sqrt{2} in the high density mean field regime to ω_h in the low density Tonks-Girardeau regime with ω_h the frequency of the harmonic trapping. Particular attention is paid to the effective mass of a soliton with velocity near the speed of sound.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Trapped one-dimensional ideal Fermi gas with a single impurity

Properties of a single impurity in a one-dimensional Fermi gas are investigated in homogeneous and trapped geometries. In a homogeneous system we use McGuire's expression [J. B. McGuire, J. Math. Phys. 6, 432 (1965)] to obtain interaction and kinetic energies, as well as the local pair correlation function. The energy of a trapped system is obtained (i) by generalizing McGuire expression (ii) within local density approximation (iii) using perturbative approach in the case of a weakly interacting impurity and (iv) diffusion Monte Carlo method. We demonstrate that a closed formula based on the exact solution of the homogeneous case provides a precise estimation for the energy of a trapped system for arbitrary coupling constant of the impurity even for a small number of fermions. We analyze energy contributions from kinetic, interaction and potential components, as well as spatial properties such as the system size. Finally, we calculate the frequency of the breathing mode. Our analysis is directly connected and applicable to the recent experiments in microtraps.

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

Super Tonks-Girardeau state in an attractive one-dimensional dipolar gas

The ground state of a one-dimensional (1D) quantum gas of dipoles oriented perpendicular to the longitudinal axis, with a strong 1/x^3 repulsive potential, is studied at low 1D densities $n$. Near contact the dependence of the many-body wave function on the separation x_{jl} of two particles reduces to a two-body wave function Ψ_{rel}(x_{jl}). Immediately after a sudden rotation of the dipoles so that they are parallel to the longitudinal axis, this wave function will still be that of the repulsive potential, but since the potential is now that of the attractive potential, it will not be stationary. It is shown that as nd^2 -> 0 the rate of change of this wave function approaches zero. It follows that for small values of nd^2, this state is metastable and is an analog of the super Tonks-Girardeau state of bosons with a strong zero-range attraction. The dipolar system is equivalent to a spinor Fermi gas with spin $z$ components σ_{\uparrow}=\perp (perpendicular to the longitudinal axis) and σ_{\downarrow}=|| (parallel to the longitudinal axis). A Fermi-Fermi mapping from spinor to spinless Fermi gas followed by the standard 1960 Fermi-Bose mapping reduces the Fermi system to a Bose gas. Potential experiments realizing the sudden spin rotation with ultracold dipolar gases are discussed, and a few salient properties of these states are accurately evaluated by a Monte Carlo method.

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

Condensate fraction of cold gases in non-uniform external potential

Exact calculation of the condensate fraction in multi-dimensional inhomogeneous interacting Bose systems which do not possess continuous symmetries is a difficult computational problem. We have developed an iterative procedure which allows to calculate the condensate fraction as well as the corresponding eigenfunction of the one-body density matrix. We successfully validate this procedure in diffusion Monte Carlo simulations of a Bose gas in an optical lattice at zero temperature. We also discuss relation between different criteria used for testing coherence in cold Bose systems, such as fraction of particles that are superfluid, condensed or are in the zero-momentum state.

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

Mesoscopic supersolid of dipoles in a trap

A mesoscopic system of indirect dipolar bosons trapped by a harmonic potential is considered. The system has a number of physical realizations including dipole excitons, atoms with large dipolar moment, polar molecules, Rydberg atoms in inhomogenious electric field. We carry out a diffusion Monte Carlo simulation to define the quantum properties of a two-dimensional system of trapped dipoles at zero temperature. In dimensionless units the system is described by two control parameters, namely the number of particles and the strength of the interparticle interaction. We have shown that when the interparticle interaction is strong enough a mesoscopic crystal is formed. As the strength of interactions is decreased a multi-stage melting takes place. Off-diagonal order in the system is tested using natural orbitals analysis. We have found that the system might be Bose-condensed even in the case of strong interparticle interactions. There is a set of parameters for which a spatially ordered structure is formed while simultaneously the fraction of Bose condensed particles is non zero. This might be considered as a realization of a mesoscopic supersolid.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical study of Kapitza pendulum

A driven pendulum with vertical oscillations of pendulum support (Kapitza pendulum) possesses a number of unusual properties and is a popular object of both analytical and numerical studies. Although some spectacular results can be obtained, such as the vertical position of the pendulum under certain conditions might become stable, no explicit analytical solution for the pendulum trajectory is known. We carry out a numerical study of Kapitza pendulum for a number of different physical regimes. Comparison is made with the limiting cases where the exact solution is known.

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
preprint2010arXiv

Exact ground state properties of the one-dimensional Coulomb gas

The ground state properties of a single-component one-dimensional Coulomb gas are investigated. We use Bose-Fermi mapping for the ground state wave function which permits to solve the Fermi sign problem in the following respects (i) the nodal surface is known, permitting exact calculations (ii) evaluation of determinants is avoided, reducing the numerical complexity to that of a bosonic system, thus allowing simulation of a large number of fermions. Due to the mapping the energy and local properties in one-dimensional Coulomb systems are exactly the same for Bose-Einstein and Fermi-Dirac statistics. The exact ground state energy has been calculated in homogeneous and trapped geometries by using the diffusion Monte Carlo method. We show that in the low-density Wigner crystal limit an elementary low-lying excitation is a plasmon, which is to be contrasted with the large-density ideal Fermi gas/Tonks-Girardeau limit, where low lying excitations are phonons. Exact density profiles are confronted to the ones calculated within the local density approximation which predicts a change from a semicircular to inverted parabolic shape of the density profile as the value of the charge is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground state of a mixture of two bosonic Calogero-Sutherland gases with strong odd-wave interspecies attraction

A model of two Calogero-Sutherland Bose gases A and B with strong odd-wave AB attractions induced by a p-wave AB Feshbach resonance is studied. The ground state wave function is found analytically by a Bose-Bose duality mapping, which permits one to accurately determine static physical properties by a Monte Carlo method. The condensation of particles or particle pairs (molecules) is tested by analyzing the presence of the off-diagonal long-range order in one- or two-body density matrices. The p-wave symmetry of AB interaction makes possible quasi-condensation of type A particles at the Fermi momentum of the B component. The zero-temperature phase diagram is drawn in terms of densities and interaction strengths.

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

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
preprint2010arXiv

Wave functions of the super Tonks-Girardeau gas and the trapped 1D hard sphere Bose gas

Recent theoretical and experimental results demonstrate a close connection between the super Tonks-Girardeau (sTG) gas and a 1D hard sphere Bose (HSB) gas with hard sphere diameter nearly equal to the 1D scattering length $a_{1D}$ of the sTG gas, a highly excited gas-like state with nodes only at interparticle separations $|x_{j\ell}|=x_{node}\approx a_{1D}$. It is shown herein that when the coupling constant $g_B$ in the Lieb-Liniger interaction $g_Bδ(x_{j\ell})$ is negative and $|x_{12}|\ge x_{node}$, the sTG and HSB wave functions for $N=2$ particles are not merely similar, but identical; the only difference between the sTG and HSB wave functions is that the sTG wave function allows a small penetration into the region $|x_{12}|<x_{node}$, whereas for a HSB gas with hard sphere diameter $a_{h.s.}=x_{node}$, the HSB wave function vanishes when all $|x_{12}|<a_{h.s.}$. Arguments are given suggesting that the same theorem holds also for $N>2$. The sTG and HSB wave functions for N=2 are given exactly in terms of a parabolic cylinder function, and for $N\ge 2$, $x_{node}$ is given accurately by a simple parabola. The metastability of the sTG phase generated by a sudden change of the coupling constant from large positive to large negative values is explained in terms of the very small overlap between the ground state of the Tonks-Girardeau gas and collapsed cluster states.

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

Thermal and quantum fluctuations in chains of ultracold polar molecules

Ultracold polar molecules, in highly anisotropic traps and interacting via a repulsive dipolar potential, may form one-dimensional chains at high densities. According to classical theory, at low temperatures there exists a critical value of the density at which a second order phase transition from a linear to a zigzag chain occurs. We study the effect of thermal and quantum fluctuations on these self-organized structures using classical and quantum Monte Carlo methods, by means of which we evaluate the pair correlation function and the static structure factor. Depending on the parameters, these functions exhibit properties typical of a crystalline or of a liquid system. We compare the thermal and the quantum results, identifying analogies and differences. Finally, we discuss experimental parameter regimes where the effects of quantum fluctuations on the linear - zigzag transition can be observed.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Number fluctuations in cold quantum gases

In ultracold gases many experiments use atom imaging as a basic observable. The resulting image is averaged over a number of realizations and mostly only this average is used. Only recently the noise has been measured to extract physical information. In the present paper we investigate the quantum noise arising in these gases at zero temperature. We restrict ourselves to the homogeneous situation and study the fluctuations in particle number found within a given volume in the gas, and more specifically inside a sphere of radius $R$. We show that zero-temperature fluctuations are not extensive and the leading term scales with sphere radius $R$ as $R^2\ln R$ (or $\ln R$) in three- (or one-) dimensional systems. We calculate systematically the next term beyond this leading order. We consider first the generic case of a compressible superfluid. Then we investigate the whole Bose-Einstein-condensation (BEC)-BCS crossover crossover, and in particular the limiting cases of the weakly interacting Bose gas and of the free Fermi gas.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Equation of state and collective frequencies of a trapped Fermi gas along the BEC-unitarity crossover

We show that the study of the collective oscillations in a harmonic trap provides a very sensitive test of the equation of state of a Fermi gas near a Feshbach resonance. Using a scaling approach, whose high accuracy is proven by comparison with exact hydrodynamic solutions, the frequencies of the lowest compressional modes are calculated at T=0 in terms of a dimensionless parameter characterizing the equation of state. The predictions for the collective frequencies, obtained from the equations of state of mean field BCS theory and of recent Monte-Carlo calculations, are discussed in detail.

0 citations0 reviewsNo code linkedOpen access