Search

Quantum mechanical few-body systems in reduced dimensionalities can exhibit many interesting properties such as scale-invariance and universality. Analytical descriptions are often available for integer dimensionality, however, numerical approaches are necessary for addressing dimensional transitions. The Fully-Correlated Gaussian method provides a variational description of the few-body real-space wavefunction. By placing the particles in a harmonic trap, the system can be described at various degrees of anisotropy by squeezing the confinement. Through this approach, configurations of two and three identical bosons as well as heteronuclear (Cs-Cs-Li and K-K-Rb) systems are described during a continuous deformation from three to one dimension. We find that the changes in binding energies between integer dimensional cases exhibit a universal behavior akin to that seen in avoided crossings or Zeldovich rearrangement., Comment: 11 pages, 5 figures, comments are most welcome

The quantum mechanical three-body problem is a source of continuing interest due to its complexity and not least due to the presence of fascinating solvable cases. The prime example is the Efimov effect where infinitely many bound states of identical bosons can arise at the threshold where the two-body problem has zero binding energy. An important aspect of the Efimov effect is the effect of spatial dimensionality; it has been observed in three dimensional systems, yet it is believed to be impossible in two dimensions. Using modern experimental techniques, it is possible to engineer trap geometry and thus address the intricate nature of quantum few-body physics as function of dimensionality. Here we present a framework for studying the three-body problem as one (continuously) changes the dimensionality of the system all the way from three, through two, and down to a single dimension. This is done by considering the Efimov favorable case of a mass-imbalanced system and with an external confinement provided by a typical experimental case with a (deformed) harmonic trap., Comment: 11 pages, 4 figures, comments are most welcome

We investigate the two lowest-lying weakly bound states of $N \leq 8$ bosons as functions of the strength of two-body Gaussian interactions. We observe the limit for validity of Efimov physics. We calculate energies and second radial moments as functions of scattering length. For identical bosons we find that two $(N-1)$-body states appear before the $N$-body ground states become bound. This pattern ceases to exist for $N \geq 7$ where the size of the ground state becomes smaller than the range of the two-body potential. All mean-square-radii for $N \geq 4$ remain finite at the threshold of zero binding, where they vary as $(N-1)^{p}$ with $p=-3/2,-3$ for ground and excited states, respectively. Decreasing the mass of one particle we find stronger binding and smaller radii. The identical particles form a symmetric system, while the lighter particle is further away in the ground states. In the excited states we find the identical bosons either surrounded or surrounding the light particle for few or many bosons, respectively. We demonstrate that the first excited states for all strengths resemble two-body halos of one particle weakly bound to a dense $N$-body system for $N=3,4$. This structure ceases to exist for $N \geq 5$., Comment: 13 pages, 10 figures, 2 tables

Matrix elements between shifted correlated Gaussians of various potentials with several form-factors are calculated analytically. Analytic matrix elements are of importance for the correlated Gaussian method in quantum few-body physics.

One-dimensional multi-component Fermi or Bose systems with strong zero-range interactions can be described in terms of local exchange coefficients and mapping the problem into a spin model is thus possible. For arbitrary external confining potentials the local exchanges are given by highly non-trivial geometric factors that depend solely on the geometry of the confinement through the single-particle eigenstates of the external potential. To obtain accurate effective Hamiltonians to describe such systems one needs to be able to compute these geometric factors with high precision which is difficult due to the computational complexity of the high-dimensional integrals involved. An approach using the local density approximation would therefore be a most welcome approximation due to its simplicity. Here we assess the accuracy of the local density approximation by going beyond the simple harmonic oscillator that has been the focus of previous studies and consider some double-wells of current experimental interest. We find that the local density approximation works quite well as long as the potentials resemble harmonic wells but break down for larger barriers. In order to explore the consequences of applying the local density approximation in a concrete setup we consider quantum state transfer in the effective spin models that one obtains. Here we find that even minute deviations in the local exchange coefficients between the exact and the local density approximation can induce large deviations in the fidelity of state transfer for four, five, and six particles., Comment: 12 pages, 7 figures, 1 table, final version

We study a quantum mechanical system consisting of up to three identical dipoles confined to move along a helical shaped trap. The long-range interactions between particles confined to move in this one dimension leads to an interesting effective two-particle potential with an oscillating behaviour. For this system we calculate the spectrum and the wave functions of the bound states. The full quantum solutions show clear imprints of the tendency for the system to form chains of dipoles along the helix, i.e. a configuration in which the dipoles are sitting approximately one winding of the helix apart so that they can take maximal advantage of the strong head-to-tail attraction that is a generic feature of the dipole-dipole interaction., Comment: 8 pages, 8 figures and 1 table

We consider two identical impurities immersed in a Fermi sea for a broad range of masses and for both interacting and non-interacting impurities. The interaction between the particles is described through attractive zero-range potentials and the problem is solved in momentum space. The two impurities can attach to a fermion from the sea and form three-body bound states. The energy of these states increase as function of the Fermi momentum $k_F$, leading to three-body bound states below the Fermi energy. The fate of the states depends highly on two- and three-body thresholds and we find evidence of medium-induced Borromean-like states in 2D. The corrections due to particle-hole fluctuations in the Fermi sea are considered in the three-body calculations and we show that in spite of the fact that they strongly affect both the two- and three-body systems, the correction to the point at which the three-body states cease to exist is small., Comment: 27 pages, 10 figures, including technical appendices, published version

We consider one particle confined to a deformed one-dimensional wire. The quantum mechanical equivalent of the classical problem is not uniquely defined. We describe several possible hamiltonians and corresponding solutions for a finite wire with fixed endpoints and non-vanishing curvature. We compute and compare the disparate eigenvalues and eigenfunctions obtained from different quantization prescriptions. The JWKB approximation without potential leads precisely to the square well spectrum and the coordinate dependent stretched or compressed box related eigenfunctions. The geometric potential arising from an adiabatic expansion in terms of curvature is at best only valid for very small curvature., Comment: 13 pages, 13 figures, 4 tables

We investigate three-body recombination rates into deep dimers in cold atomic gases with large scattering length within hyper-spherical adiabatic zero-range approach. We derive closed analytic expressions for the rates for one- and two-species gases. Although the deep dimers are beyond the zero-range theory the latter can still describe the recombination into deep dimers by use of one additional short-range absorption parameter. The recombination rate, as function of the scattering length, retains the known universal behavior --- the fourth power trend with characteristic log-periodic peaks --- however increasing the short-range absorption broadens the peaks until they are eventually completely smeared out. Increasing the heavy-to-light mass ratio in a two-species system decreases the distance between the peaks and increases the overal scale of the recombination rate., Comment: 8 pages, 5 figures, Conference "Systems on the verge of the stability" (Critical Stability 2014)

We consider three-body recombination into deep dimers in a mass-imbalanced two-component atomic gas. We use an optical model where a phenomenological imaginary potential is added to the lowest adiabatic hyper-spherical potential. The consequent imaginary part of the energy eigenvalue corresponds to the decay rate or recombination probability of the three-body system. The method is formulated in details and the relevant qualitative features are discussed as functions of scattering lengths and masses. We use zero-range model in analyses of recent recombination data. The dominating scattering length is usually related to the non-equal two-body systems. We account for temperature smearing which tends to wipe out the higher-lying Efimov peaks. The range and the strength of the imaginary potential determine positions and shapes of the Efimov peaks as well as the absolute value of the recombination rate. The Efimov scaling between recombination peaks is calculated and shown to depend on both scattering lengths. Recombination is predicted to be largest for heavy-heavy-light systems. Universal properties of the optical parameters are indicated. We compare to available experiments and find in general very satisfactory agreement., Comment: 16 pages, 12 figures, 3 tables, revised version

We describe a procedure to solve an up to $2N$ problem where the particles are separated topologically in $N$ groups with at most two particles in each. Arbitrary interactions are allowed between the (two) particles within one group. All other interactions are approximated by harmonic oscillator potentials. The problem is first reduced to an analytically solvable $N$-body problem and $N$ independent two-body problems. We calculate analytically spectra, wave functions, and normal modes for both the inverse square and delta-function two-body interactions. In particular, we calculate separation energies between two strings of particles. We find that the string separation energy increases with $N$ and interaction strength., Comment: 13 pages, 7 figures, revised and published version

Strongly interacting one-dimensional quantum systems often behave in a manner that is distinctly different from their higher-dimensional counterparts. When a particle attempts to move in a one-dimensional environment it will unavoidably have to interact and 'push' other particles in order to execute a pattern of motion, irrespective of whether the particles are fermions or bosons. A present frontier in both theory and experiment are mixed systems of different species and/or particles with multiple internal degrees of freedom. Here we consider trapped two-component bosons with short-range inter-species interactions much larger than their intra-species interactions and show that they have novel energetic and magnetic properties. In the strongly interacting regime, these systems have energies that are fractions of the basic harmonic oscillator trap quantum and have spatially separated ground states with manifestly ferromagnetic wave functions. Furthermore, we predict excited states that have perfect antiferromagnetic ordering. This holds for both balanced and imbalanced systems, and we show that it is a generic feature as one crosses from few- to many-body systems., Comment: 9 pages, 6 figures

We investigate a two-dimensional system of with two values of the internal (spin) degree of freedom. It is confined by a deformed harmonic trap and subject to a Zeeman field, Rashba or Dresselhaus one-body spin-orbit couplings and two-body short range repulsion. We obtain self-consistent mean-field $N$-body solutions as functions of the interaction parameters. Single-particle Spectra and total energies are computed and compared to the results without interaction. We perform a statistical analysis for the distributions of nearest neighbor energy level spacings and show that quantum signatures of chaos are seen in certain parameters regimes. Furthermore, the effects of two-body repulsion on the nearest neighbor distributions are investigated. This repulsion can either promote or destroy the signatures of potential chaotic behavior depending on relative strengths of parameters. Our findings support the suggestion that cold atoms may be used to study quantum chaos both in the presence and absence of interactions., Comment: 12 pages, 9 figures, revised version

We examine a one-dimensional two-component fermionic system in a trap, assuming that all particles have the same mass and interact through a strong repulsive zero-range force. First we show how a simple system of three strongly interacting particles in a harmonic trap can be treated using the hyperspherical formalism. Next we discuss the behavior of the energy for the N-body system., Comment: 5 pages. Original paper for EPJ ST in connection with the workshop BEC2014 28-31 May 2014 in Levico Terme, Italy

We consider a one-dimensional gas of cold atoms with strong contact interactions and construct an effective spin-chain Hamiltonian for a two-component system. The resulting Heisenberg spin model can be engineered by manipulating the shape of the external confining potential of the atomic gas. We find that bosonic atoms offer more flexibility for tuning independently the parameters of the spin Hamiltonian through interatomic (intra-species) interaction which is absent for fermions due to the Pauli exclusion principle. Our formalism can have important implications for control and manipulation of the dynamics of few- and many-body quantum systems; as an illustrative example relevant to quantum computation and communication, we consider state transfer in the simplest non-trivial system of four particles representing exchange-coupled qubits., Comment: 10 pages including appendix, 3 figures, revised version

The spectrum and properties of quantum bound states is strongly dependent on the dimensionality of space. How this comes about and how one may theoretically and experimentally study the interpolation between different dimensions is a topic of great interest in different fields of physics. In this paper we study weakly bound states of non-relativistic two and three boson systems when passing continuously from a three (3D) to a two-dimensional (2D) regime within a 'squeezed dimension' model. We use periodic boundary conditions to derive a surprisingly simple form of the three-boson Schr{\"o}dinger equation in momentum space that we solve numerically. Our results show a distinct dimensional crossover as three-boson states will either disappear into the continuum or merge with a 2D counterpart, and also a series of sharp transitions in the ratios of three-body and two-body energies from being purely 2D to purely 3D., Comment: 5 pages, 2 figures, revised version

We compute single-particle energy spectra for a one-body hamiltonian consisting of a two-dimensional deformed harmonic oscillator potential, the Rashba spin-orbit coupling and the Zeeman term. To investigate the statistical properties of the obtained spectra as functions of deformation, spin-orbit and Zeeman strengths we examine the distributions of the nearest neighbor spacings. We find that the shapes of these distributions depend strongly on the three potential parameters. We show that the obtained shapes in some cases can be well approximated with the standard Poisson, Brody and Wigner distributions. The Brody and Wigner distributions characterize irregular motion and help identify quantum chaotic systems. We present a special choices of deformation and spin-orbit strengths without the Zeeman term which provide a fair reproduction of the fourth-power repelling Wigner distribution. By adding the Zeeman field we can reproduce a Brody distribution, which is known to describe a transition between the Poisson and linear Wigner distributions., Comment: 12 pages, 10 figures, published version

We consider crystal formation of particles with dipole-dipole interactions that are confined to move in a one-dimensional helical geometry with their dipole moments oriented along the symmetry axis of the confining helix. The stable classical lowest energy configurations are found to be chain structures for a large range of pitch-to-radius ratios for relatively low density of dipoles and a moderate total number of particles. The classical normal mode spectra support the chain interpretation both through structure and the distinct degeneracies depending discretely on the number of dipoles per revolution. A larger total number of dipoles leads to a clusterization where the dipolar chains move closer to each other. This implies a change in the local density and the emergence of two length scales, one for the cluster size and one for the inter-cluster distance along the helix. Starting from three dipoles per revolution, this implies a breaking of the initial periodicity to form a cluster of two chains close together and a third chain removed from the cluster. This is driven by the competition between in-chain and out-of-chain interactions, or alternatively the side-by-side repulsion and the head-to-tail attraction in the system. The speed of sound propagates along the chains. It is independent of the number of chains although depending on geometry., Comment: 13 pages, 13 figures, revised version

The study of quantum mechanical bound states is as old as quantum theory itself. Yet, it took many years to realize that three-body borromean systems that are bound when any two-body subsystem is unbound are abundant in nature. Here we demonstrate the existence of borromean systems of spin-polarized (spinless) identical fermions in two spatial dimensions. The ground state with zero orbital (planar) angular momentum exists in a borromean window between critical two- and three-body strengths. The doubly degenerate first excited states of angular momentum one appears only very close to the two-body threshold. They are the lowest in a possible sequence of so-called super-Efimov states. While the observation of the super-Efimov scaling could be very difficult, the borromean ground state should be observable in cold atomic gases and could be the basis for producing a quantum gas of three-body states in two dimensions., Comment: 9 pages, 3 figures, published version

We study two-component bosonic systems with strong inter-species and vanishing intra-species interactions. A new class of exact eigenstates is found with energies that are {\it not} sums of the single-particle energies with wave functions that have the characteristic feature that they vanish over extended regions of coordinate space. This is demonstrated in an analytically solvable model for three equal mass particles, two of which are identical bosons, which is exact in the strongly-interacting limit. We numerically verify our results by presenting the first application of the stochastic variational method to this kind of system. We also demonstrate that the limit where both inter- and intra-component interactions become strong must be treated with extreme care as these limits do not commute. Moreover, we argue that such states are generic also for general multi-component systems with more than three particles. The states can be probed using the same techniques that have recently been used for fermionic few-body systems in quasi-1D., Comment: 6 pages, 4 figures, published version

We study the two dimensional three-body problem in the general case of three distinguishable particles interacting through zero-range potentials. The Faddeev decomposition is used to write the momentum-space wave function. We show that the large-momentum asymptotic spectator function has the same functional form as derived previously for three identical particles. We derive analytic relations between the three different Faddeev components for three distinguishable particles. We investigate the one-body momentum distributions both analytically and numerically and analyze the tail of the distributions to obtain two- and three-body contact parameters. We specialize from the general cases to examples of two identical, interacting or non-interacting, particles. We find that the two-body contact parameter is not a universal constant in the general case and show that the universality is recovered when a subsystem is composed of two identical non-interacting particles. We also show that the three-body contact parameter is negligible in the case of one non-interacting subsystem compared to the situation where all subsystem are bound. As example, we present results for mixtures of Lithium with two Cesium or two Potassium atoms, which are systems of current experimental interest., Comment: 24 pages, 7 figures, published version

We investigate three-boson recombination of equal mass systems as function of (negative) scattering length, mass, finite energy, and finite temperature. An optical model with an imaginary potential at short distance reproduces experimental recombination data and allows us to provide a simple parametrization of the recombination rate as function of scattering length and energy. Using the two-body van der Waals length as unit we find that the imaginary potential range and also the potential depth agree to within thirty percent for Lithium and Cesium atoms. As opposed to recent studies suggesting universality of the threshold for bound state formation, our results suggest that the recombination process itself could have universal features., Comment: 5 pages, 5 figures

In one dimension, the study of magnetism dates back to the dawn of quantum mechanics when Bethe solved the famous Heisenberg model that describes quantum behaviour in magnetic systems. In the last decade, one-dimensional systems have become a forefront area of research driven by the realization of the Tonks-Girardeau gas using cold atomic gases. Here we prove that one-dimensional fermionic and bosonic systems with strong short-range interactions are solvable in arbitrary confining geometries by introducing a new energy-functional technique and obtaining the full spectrum of energies and eigenstates. As a first application, we calculate spatial correlations and show how both ferro- and anti-ferromagnetic states are present already for small system sizes that are prepared and studied in current experiments. Our work demonstrates the enormous potential for quantum manipulation of magnetic correlations at the microscopic scale., Comment: 11 pages, 2 figures, including methods, final version

We solve the three-body bound state problem in three dimensions for mass imbalanced systems of two identical bosons and a third particle in the universal limit where the interactions are assumed to be of zero-range. The system displays the Efimov effect and we use the momentum-space wave equation to derive formulas for the scaling factor of the Efimov spectrum for any mass ratio assuming either that two or three of the two-body subsystems have a bound state at zero energy. We consider the single-particle momentum distribution analytically and numerically and analyse the tail of the momentum distribution to obtain the three-body contact parameter. Our finding demonstrate that the functional form of the three-body contact term depends on the mass ratio and we obtain an analytic expression for this behavior. To exemplify our results, we consider mixtures of Lithium with either two Caesium or Rubium atoms which are systems of current experimental interest., Comment: 16 pages, 9 figures, 1 appendix, revised version

We consider one-dimensional tubes containing bosonic polar molecules. The long-range dipole-dipole interactions act both within a single tube and between different tubes. We consider arbitrary values of the externally aligned dipole moments with respect to the symmetry axis of the tubes. The few-body structures in this geometry are determined as function of polarization angles and dipole strength by using both essentially exact stochastic variational methods and the harmonic approximation. The main focus is on the three, four, and five-body problems in two or more tubes. Our results indicate that in the weakly-coupled limit the inter-tube interaction is similar to a zero-range term with a suitable rescaled strength. This allows us to address the corresponding many-body physics of the system by constructing a model where bound chains with one molecule in each tube are the effective degrees of freedom. This model can be mapped onto one-dimensional Hamiltonians for which exact solutions are known., Comment: 22 pages, 7 figures, revised version

We consider a spin-orbit coupled system of particles in an external trap that is represented by a deformed harmonic oscillator potential. The spin-orbit interaction is a Rashba interaction that does not commute with the trapping potential and requires a full numerical treatment in order to obtain the spectrum. The effect of a Zeeman term is also considered. Our results demonstrate that variable spectral gaps occur as a function of strength of the Rashba interaction and deformation of the harmonic trapping potential. The single-particle density of states and the critical strength for superfluidity vary tremendously with the interaction parameter. The strong variations with Rashba coupling and deformation implies that the few- and many-body physics of spin-orbit coupled systems can be manipulated by variation of these parameters., Comment: 18 pages, 6 figures, published version

We investigate finite-range effects in systems with three identical bosons. We calculate recombination rates and bound state spectra using two different finite-range models that have been used recently to describe the physics of cold atomic gases near Feshbach resonances where the scattering length is large. The models are built on contact potentials which take into account finite range effects; one is a two-channel model and the other is an effective range expansion model implemented through the boundary condition on the three-body wave function when two of the particles are at the same point in space. We compare the results with the results of the ubiquitous single-parameter zero-range model where only the scattering length is taken into account. Both finite range models predict variations of the well-known geometric scaling factor 22.7 that arises in Efimov physics. The threshold value at negative scattering length for creation of a bound trimer moves to higher or lower values depending on the sign of the effective range compared to the location of the threshold for the single-parameter zero-range model. Large effective ranges, corresponding to narrow resonances, are needed for the reduction to have any noticeable effect., Comment: 16 pages, 7 figures, published version

We consider three-body systems in two dimensions with zero-range interactions for general masses and interaction strengths. The momentum-space Schr\"odinger equation is solved numerically and in the Born-Oppenheimer (BO) approximation. The BO expression is derived using separable potentials and yields a concise adiabatic potential between the two heavy particles. The BO potential is Coulomb-like and exponentially decreasing at small and large distances, respectively. While we find similar qualitative features to previous studies, we find important quantitative differences. Our results demonstrate that mass-imbalanced systems that are accessible in the field of ultracold atomic gases can have a rich three-body bound state spectrum in two dimensional geometries. Small light-heavy mass ratios increase the number of bound states. For 87Rb-87Rb-6Li and 133Cs-133Cs-6Li we find respectively 3 and 4 bound states., Comment: 17 pages, 8 figures, revised version

We search for Borromean three-body systems of identical bosons in two dimensional geometry, i.e. we search for bound three-boson system without bound two-body subsystems. Unlike three spatial dimensions, in two-dimensional geometry the two- and three-body thresholds often coincide ruling out Borromean systems. We show that Borromean states can only appear for potentials with substantial attractive and repulsive parts. Borromean states are most easily found when a barrier is present outside an attractive pocket. Extensive numerical search did not reveal Borromean states for potentials without an outside barrier. We outline possible experimental setups to observe Borromean systems in two spatial dimensions., Comment: 9 pages, 5 figures, published version

The momentum distribution is a powerful probe of strongly-interacting systems that are expected to display universal behavior. This is contained in the contact parameters which relate few- and many-body properties. Here we consider a Bose gas in two dimensions and explicitly show that the two-body contact parameter is universal and then demonstrate that the momentum distribution at next-to-leading order has a logarithmic dependence on momentum which is vastly different from the three-dimensional case. Based on this, we propose a scheme for measuring the effective dimensionality of a quantum many-body system by exploiting the functional form of the momentum distribution., Comment: 8 pages, 5 figures, final version

We calculate shallow three-body bound states in the universal regime, defined by Efimov, with inclusion of both scattering length and effective range parameters. The universal spectrum is recovered for the least bound states, whereas for larger binding energies we find corrections to the universal scaling laws. We recover known results for broad Feshbach resonances with small effective range, whereas in the case of narrow resonances we find a distinct non-monotonic behavior of the threshold at which the lowest Efimov trimer merges with the three-body continuum. To address the issue of the physical origin of the three-body parameter we provide a physically clear model for the relation between three-body physics and typical two-body atom-atom interactions. Our results demonstrate that experimental information from narrow Feshbach resonances and/or mixed systems are of vital importance to pin down the relation of two- and three-body physics in atomic systems., Comment: 6 pages, 6 figures, 1 appendix, accepted version

The virial expansion method is applied within a harmonic approximation to an interacting N-body system of identical fermions. We compute the canonical partition functions for two and three particles to get the two lowest orders in the expansion. The energy spectrum is carefully interpolated to reproduce ground state properties at low temperature and the non-interacting large temperature limit of constant virial coefficients. This resembles the smearing of shell effects in finite systems with increasing temperature. Numerical results are discussed for the second and third virial coefficients as function of dimension, temperature, interaction, and the transition temperature between low and high energy limits., Comment: 11 pages, 7 figures, published version

We consider energetics and structural properties of a many particle system in one dimension with pairwise contact interactions confined in a parabolic external potential. To render the problem analytically solvable, we use the harmonic approximation scheme at the level of the Hamiltonian. We investigate the scaling with particle number of the ground state energies for systems consisting of identical bosons or fermions. We then proceed to focus on bosonic systems and make a detailed comparison to known exact results in the absence of the parabolic external trap for three-body systems. We also consider the thermodynamics of the harmonic model which turns out to be similar for bosons and fermions due to the lack of degeneracy in one dimension., Comment: 7+epsilon pages, 7 figures, published version

We consider few-body bound state systems and provide precise definitions of Borromean and Brunnian systems. The initial concepts are more than a hundred years old and originated in mathematical knot-theory as purely geometric considerations. About thirty years ago they were generalized and applied to the binding of systems in nature. It now appears that recent generalization to higher order Brunnian structures may potentially be realized as laboratory made or naturally occurring systems. With the binding energy as measure, we discuss possibilities of physical realization in nuclei, cold atoms, and condensed matter systems. Appearance is not excluded. However, both the form and the strengths of the interactions must be rather special. The most promising subfields for present searches would be in cold atoms because of external control of effective interactions, or perhaps in condensed-matter systems with non-local interactions. In nuclei, it would only be by sheer luck due to a lack of tunability., Comment: 8 pages, 5 figures, revised version

We describe a method to compute thermodynamic quantities in the harmonic approximation for identical bosons and fermions in an external confining field. We use the canonical partition function where only energies and their degeneracies enter. The number of states of given energy and symmetry is found by separating the center of mass motion, and counting the remaining states of given symmetry and excitation energy of the relative motion. The oscillator frequencies that enter the harmonic Hamiltonian can be derived from realistic model parameters and the method corresponds to an effective interaction approach based on harmonic interactions. To demonstrate the method, we apply it to systems in two dimensions. Numerical calculations are compared to a brute force method that is considerably more computationally intensive., Comment: 12 pages, 8 figures

The thermodynamics of a quantum system of layers containing perpendicularly oriented dipolar molecules is studied within an oscillator approximation for both bosonic and fermionic species. The system is assumed to be built from chains with one molecule in each layer. We consider the effects of the intralayer repulsion and quantum statistical requirements in systems with more than one chain. Specifically, we consider the case of two chains and solve the problem analytically within the harmonic Hamiltonian approach which is accurate for large dipole moments. The case of three chains is calculated numerically. Our findings indicate that thermodynamic observables, such as the heat capacity, can be used to probe the signatures of the intralayer interaction between chains. This should be relevant for near future experiments on polar molecules with strong dipole moments., Comment: 15 pages, 5 figures, final version

Ultracold polar molecules in multilayered systems have been experimentally realized very recently. While experiments study these systems almost exclusively through their chemical reactivity, the outlook for creating and manipulating exotic few- and many-body physics in dipolar systems is fascinating. Here we concentrate on few-body states in a multilayered setup. We exploit the geometry of the interlayer potential to calculate the two- and three-body chains with one molecule in each layer. The focus is on dipoles that are aligned at some angle with respect to the layer planes by means of an external eletric field. The binding energy and the spatial structure of the bound states are studied in several different ways using analytical approaches. The results are compared to stochastic variational calculations and very good agreement is found. We conclude that approximations based on harmonic oscillator potentials are accurate even for tilted dipoles when the geometry of the potential landscape is taken into account., Comment: 10 pages, 6 figures. Submitted to Few-body Systems special issue on Critical Stability, revised version

We consider bound states of asymmetric three-body systems confined to two dimensions. In the universal regime, two energy ratios and two mass ratios provide complete knowledge of the three-body energy measured in units of one two-body energy. The lowest number of stable bound states is produced when one mass is larger than two similar masses. We focus on selected asymmetric systems of interest in cold atom physics. The scaled three-body energy and the two scaled two-body energies are related through an equation for a supercircle whose radius increases almost linearly with three-body energy. The exponents exhibit an increasing behavior with three-body energy. The mass dependence is highly non-trivial. We give a simple relation that predicts the universal three-body energy., Comment: 5 pages, 4 figures, final version

Bound structures among dipolar molecules in multilayers are a topic of great interest in the light of recent experiments that have demonstrated the feasibility of the setup. While it is known that two molecules in two adjacent layers will always bind, larger complexes have only been scarcely addressed thus far. Here we prove rigorously that three- and four-body states will never be bound when the dipoles are oriented perpendicular to the layers. The technique employed is general and can be used for more molecules/layers and other geometries. Our analytical findings are supported by numerical calculations for both fermionic and bosonic molecules. Furthermore, we calculate the reduction in intralayer repulsion necessary to bind large complexes and estimate the influence of bound complexes in systems with many layers., Comment: 5 pages, 4 figures, final version

The momentum space zero-range model is used to investigate universal properties of three interacting particles confined to two dimensions. The pertinent equations are first formulated for a system of two identical and one distinct particle and the two different two-body subsystems are characterized by two-body energies and masses. The three-body energy in units of one of the two-body energies is a universal function of the other two-body energy and the mass ratio. We derive convenient analytical formulae for calculations of the three-body energy as function of these two independent parameters and exhibit the results as universal curves. In particular, we show that the three-body system can have any number of stable bound states. When the mass ratio of the distinct to identical particles is greater than 0.22 we find that at most two stable bound states exist, while for two heavy and one light mass an increasing number of bound states is possible. The specific number of stable bound states depends on the ratio of two-body bound state energies and on the mass ratio and we map out an energy-mass phase-diagram of the number of stable bound states. Realizable systems of both fermions and bosons are discussed in this framework., Comment: 21 pages, 8 figures, 4 tables, final version

We consider the N-body problem in a layered geometry containing cold polar molecules with dipole moments that are polarized perpendicular to the layers. A harmonic approximation is used to simplify the hamiltonian and bound state properties of the two-body inter-layer dipolar potential are used to adjust this effective interaction. To model the intra-layer repulsion of the polar molecules, we introduce a repulsive inter-molecule potential that can be parametrically varied. Single chains containing one molecule in each layer, as well as multi-chain structures in many layers are discussed and their energies and radii determined. We extract the normal modes of the various systems as measures of their volatility and eventually of instability, and compare our findings to the excitations in crystals. We find modes that can be classified as either chains vibrating in phase or as layers vibrating against each other. The former correspond to acoustic and the latter to optical phonons. Instabilities can occur for large intra-layer repulsion and produce diverging amplitudes of molecules in the outer layers. Lastly, we consider experimentally relevant regimes to observe the structures., Comment: 17 pages, 20 figures, accepted version

We consider a bilayer setup with two parallel planes of cold fermionic polar molecules when the dipole moments are oriented perpendicular to the planes. The binding energy of two-body states with one polar molecule in each layer is determined and compared to various analytic approximation schemes in both coordinate- and momentum-space. The effective interaction of two bound dimers is obtained by integrating out the internal dimer bound state wave function and its robustness under analytical approximations is studied. Furthermore, we consider the effect of the background of other fermions on the dimer state through Pauli blocking, and discuss implications for the zero-temperature many-body phase diagram of this experimentally realizable system., Comment: 18 pages, 10 figures, accepted version

We calculate the energy and wave functions of two particles confined to two spatial dimensions interacting via arbitrary anisotropic potentials with negative or zero net volume. The general rigorous analytic expressions are given in the weak coupling limit where universality or model independence are approached. The monopole part of anisotropic potentials is crucial in the universal limit. We illustrate the universality with a system of two arbitrarily polarized cold dipolar molecules in a bilayer. We discuss the transition to universality as function of polarization and binding energy, and compare analytic and numerical results obtained by the stochastic variational method. The universal limit is essentially reached for experimentally accessible strengths., Comment: 4.1 pages, 3 figures, published version

The exploration of cold polar molecules in different geometries is a rapidly developing experimental and theoretical pursuit. Recently, the implementation of optical lattices has enabled confinement in stacks of planes, the number of which is also controllable. Here we consider the bound state structure of two polar molecules confined in two adjacent planes as function of the polarization angle of the dipole moment of the molecules. We prove analytically and present numerical evidence for the existence of bound states for arbitrary dipole moments and polarization directions in this two-dimensional geometry. The spatial structure of the bound states is dominated by two-dimensional s- and p-waves, where the latter exceeds 40 percent over a large range of polarization angles for intermediate or strong dipole strength. Finally, we consider the influence of the dimer bound states on the potential many-body ground-state of the system., Comment: 25 pages, 9 figures, final version

The recent experimental realization of cold polar molecules in the rotational and vibrational ground state opens the door to the study of a wealth of phenomena involving long-range interactions. By applying an optical lattice to a gas of cold polar molecules one can create a layered system of planar traps. Due to the long-range dipole-dipole interaction one expects a rich structure of bound complexes in this geometry. We study the bilayer case and determine the two-body bound state properties as a function of the interaction strength. The results clearly show that a least one bound state will always be present in the system. In addition, bound states at zero energy show universal behavior and extend to very large radii. These results suggest that non-trivial bound complexes of more than two particles are likely in the bilayer and in more complicated chain structures in multi-layer systems., Comment: 6 pages, 5 figures. Revised version to be published

We consider a system of three identical bosons near a Feshbach resonance in the universal regime with large scattering length usually described by model independent zero-range potentials. We employ the adiabatic hyperspherical approximation and derive the rigorous large-distance equation for the adiabatic potential for finite-range interactions. The effective range correction to the zero-range approximation must be supplemented by a new term of the same order. The non-adiabatic term can be decisive. Efimov physics is always confined to the range between effective range and scattering length. The analytical results agree with numerical calculations for realistic potentials., Comment: 4 pages, 2 figures

We demonstrate the possibility of existence of meta-stable N-body Efimov states in trapped Bose systems with large scattering length. We calculate spectra of trapped systems of N=3,4,5,6, and 7 bosons using a stochastic variational method with a restricted correlated Gaussian basis. For each system the calculations reveal a series of Efimov states where the energy and the r.m.s. radius exhibit the characteristic exponential dependence upon the state number. We also estimate the contribution of these states to the recombination rate of Bose-Einstein condensates.

We calculate the energy and the condensate fraction of a system of trapped bosons interacting via a short-range two-body potential with positive scattering length. The potential is attractive and has a two-body bound state. When the scattering length is small compared to the trap length the system is model independent: all potential models -- attractive, repulsive and zero-range -- provide similar results. When the scattering length is large the attractive model differs qualitatively from the repulsive and zero-range models. In this regime the system with attractive potential becomes independent of the scattering length, with both the energy and the condensate fraction converging towards finite constants.

We propose a coordinate-space regularization of the three-body problem with zero-range potentials. We include the effective range and the shape parameter in the boundary condition of the zero-range potential. The proposed extended zero-range model is tested against atomic helium trimers and is shown to provide an adequate quantitative description of these systems.

The zero-range potential is customarily employed in various mean-field calculations of many-body systems in atomic and nuclear physics within, correspondingly, Gross-Pitaevskii and Skyrme-Hartree-Fock approach. We argue, however, that a many-body system with zero-range potentials is unstable against clusterization into collapsed three-body subsystems. We show that neither the density dependence of the potential nor an additional repulsive three-body potential can prevent this unexpected correlational collapse if the potentials are of zero range. Therefore the zero-range potential can only be used in many-body calculations where all three-body correlations are explicitly excluded., Comment: Corrected reference