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.
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.
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 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.
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.
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.
We study the low-temperature states of a system of spinless bosons trapped in a broad spherical shell. In our virtual (path-integral Monte Carlo) experiments, we specifically examine how the equilibrium structure and overall superfluidity of the system are affected by a change of the interparticle-potential core from hard to increasingly soft. In this manner we bridge the gap between the behaviors of two species of bosons: Helium-like atoms, which self-assemble into ordinary crystalline solids, and Rydberg-dressed atoms, which in simulation form cluster (super)solids. We also vary the particle mass, i.e., the quantumness of the system, so as to modify the relative stability of ``solid'' and fluid states. We thus find a whole spectrum of (finite-size) phases with a curved structure, including multilayer and cluster phases, as well as superfluids and single layer supersolids. In particular, we confirm that supersolids only occur in the presence of clusters. By considering three system sizes (from 60 to 180 particles), but rescaling the trap radius so that the density stays fixed, we verify that supersolid states only exist in the same range of particle masses. Our model serves as a platform to test how emergent properties are modified, in spherical geometry, in the crossover from atoms to dressed atoms.
By means of a mean-field method, we have studied the zero-temperature structure and excitation spectrum of a three-dimensional soft-core bosonic system for a value of the interaction strength that favors a crystal structure made of atomic nanoclusters arranged with fcc ordering. In addition to the longitudinal and transverse phonon branches expected for a normal crystal, the excitation spectrum shows a soft mode related to the breaking of gauge symmetry, which signals a partial superfluid character of the solid. Additional evidence of supersolidity is provided by the calculation of the superfluid fraction, which shows a first-order drop, from 1 to 0.4, at the liquid-supersolid transition and a monotonic decrease as the interaction strength parameter is increased. The conditions for the coexistence of the supersolid with the homogeneous superfluid are discussed, and the surface tension of a representative solid-liquid interface is calculated.