The properties of a macroscopic assembly of weakly repulsive bosons at zero temperature are well described by Gross-Pitaevskii mean-field theory. According to this formalism the system exhibits a quantum transition from superfluid to cluster supersolid as a function of pressure. We develop a thermodynamically rigorous treatment of the different phases of the system by adopting a variational formulation of the condensate wave function---represented as a sum of Gaussians---that is amenable to exact manipulations. Not only is this description quantitatively accurate, but it is also capable to predict the order (and sometimes even the location) of the transition. We consider a number of crystal structures in two and three dimensions and determine the phase diagram. Depending on the lattice, the transition from fluid to solid can be first order or continuous, a lower coordination entailing a milder transition. In two dimensions, crystallization would occur at the same pressure on three distinct lattices (square, honeycomb, and stripes), all providing metastable phases with respect to the triangular crystal. A similar scenario holds in three dimensions, where the simple-cubic and diamond crystals also share a common melting point; however, the stable crystal at low pressure is typically fcc. Upon compression and depending on the shape of the potential, the fcc crystal may transform into hcp. We conclude by sketching a theory of the solid-fluid interface and of quantum nucleation of the solid from the fluid.
The properties of a macroscopic assembly of weakly repulsive bosons at zero temperature are well described by Gross-Pitaevskii mean-field theory. According to this formalism the system exhibits a quantum transition from superfluid to cluster supersolid as a function of pressure. We develop a thermodynamically rigorous treatment of the different phases of the system by adopting a variational formulation of the condensate wave function---represented as a sum of Gaussians---that is amenable to exact manipulations. Not only is this description quantitatively accurate, but it is also capable to predict the order (and sometimes even the location) of the transition. We consider a number of crystal structures in two and three dimensions and determine the phase diagram. Depending on the lattice, the transition from fluid to solid can be first order or continuous, a lower coordination entailing a milder transition. In two dimensions, crystallization would occur at the same pressure on three distinct lattices (square, honeycomb, and stripes), all providing metastable phases with respect to the triangular crystal. A similar scenario holds in three dimensions, where the simple-cubic and diamond crystals also share a common melting point; however, the stable crystal at low pressure is typically fcc. Upon compression and depending on the shape of the potential, the fcc crystal may transform into hcp. We conclude by sketching a theory of the solid-fluid interface and of quantum nucleation of the solid from the fluid.
We work out the ground-state diagram of weakly-repulsive penetrable bosons, using mean-field theory with a Gaussian ansatz on the single-particle wave function. Upon compression, the fluid transforms into a cluster supersolid, whose structure is characterized for various choices of the embedding space. In Euclidean space, the stable crystals are those with the most compact structure, i.e., triangular and fcc in two and three dimensions, respectively. For particles confined in a spherical surface, as the sphere radius increases we observe a sequence of transitions between different cluster phases, all having a regular or semiregular polyhedron as supporting frame for the clusters. The present results are relevant for the behavior of ultracold bosons weakly coupled to a Rydberg state.
Abstract We study a system of penetrable bosons on a line, focusing on the high-density/weak-interaction regime, where the ground state is, to a good approximation, a condensate. Under compression, the system clusterizes at zero temperature, i.e. particles gather together in separate, equally populated bunches. We compare predictions from the Gross–Pitaevskii (GP) equation with those of two distinct variational approximations of the single-particle state, written as either a sum of Gaussians or the square root of it. Not only the wave functions in the three theories are similar, but also the phase-transition density is the same for all. In particular, clusterization occurs together with the softening of roton excitations in GP theory. Compared to the latter theory, Gaussian variational theory has the advantage that the mean-field energy functional is written in (almost) closed form, which enables us to extract the phase-transition and high-density behaviors in fully analytic terms. We also compute the superfluid fraction of the clustered system, uncovering its exact behavior close, as well as very far away from, the transition.
We study a system of penetrable bosons embedded in a spherical surface. Under the assumption of weak interaction between the particles, the ground state of the system is, to a good approximation, a pure condensate. We employ thermodynamic arguments to investigate, within a variational ansatz for the single-particle state, the crossover between distinct finite-size ``phases'' in the parameter space spanned by the sphere radius and the chemical potential. In particular, for radii up to a few interaction ranges we examine the stability of the fluid phase with respect to a number of crystal-like arrangements having the symmetry of a regular or semiregular polyhedron. We find that, while quantum fluctuations keep the system fluid at low density, upon compression it eventually becomes inhomogeneous, i.e., particles gather together in clusters. As the radius increases, the nature of the high-density aggregate varies and we observe a sequence of transitions between different cluster phases (``solids''), whose underlying rationale is to maximize the coordination number of clusters, while ensuring at the same time the proper distance between each neighboring pair. We argue that, at least within our mean-field description, every cluster phase is supersolid.
The supersolid phase of soft-core bosons in two dimensions is investigated using the self-consistent Hartree-Fock and quantum Monte Carlo methods. An approximate phase diagram at finite temperatures is initially constructed using the mean-field approach, which is subsequently validated through precise path-integral simulations, enabling a microscopic characterization of the various phases. Superfluid and melting/freezing transitions are analyzed through the superfluid density and the long-range behavior of correlation functions associated with positional and orientational order, in accordance with the general picture of Berezinskii-Kosterlitz-Thouless transitions. A broad region at low temperatures is identified where the supersolid phase exists, separating the uniform superfluid phase from the normal quasicrystal phase. Additionally, a potential intermediate hexatic phase with quasi-long-range orientational order is identified in a narrow region between the normal solid and fluid phases. These findings establish self-consistent Hartree-Fock theory beyond the local density approximation as an effective tool, complementary to computationally intensive quantum Monte Carlo simulations, for investigating the melting of the supersolid phase and the possible emergence of the hexatic superfluid phase in bosonic systems with various interaction potentials.
The ground state of a two-dimensional (2D) system of Bose particles of spin zero, interacting via a repulsive Gaussian-core potential, has been investigated by means of quantum Monte Carlo simulations. The quantum phase diagram is qualitatively identical to that of 2D Yukawa bosons. While the system is a fluid at all densities for weak coupling, in the strong-coupling regime it transitions upon compression from a low-density superfluid to a crystal, and then into a reentrant superfluid phase. No evidence of a (supersolid) cluster crystal phase is seen.
The low temperature phase diagram of an assembly of Bose particles, interacting via a simple soft disc potential, is studied numerically by means of Monte Carlo simulations. At sufficiently high densities, the system features a cluster crystal phase, as observed in its classical counterpart. In the low temperature limit, and in a relatively wide region of parameter space, the quantum cluster crystal undergoes a transition to a supersolid phase, displaying both a finite superfluid fraction and crystalline order. We argue that soft-core Bosons ought be regarded as a ‘minimal model’ of continuous-space supersolidity. With the aim of elucidating the interplay between superfluid response and solid order, we compute the excitation spectrum of the system. A comparison of that of the insulating cluster crystal with that of the supersolid, shows that the latter features an additional acoustic mode.