In the present work, we study a numerical approach to one dimensional finite volume three-body interaction, the method is demonstrated by considering a toy model of three spinless particles interacting with pair-wise $\ensuremath{\delta}$-function potentials. The numerical results are compared with the exact solutions of three spinless bosons interaction when the strength of short-range interactions are set equal for all pairs.
In the present work, we study a numerical approach to one dimensional finite volume three-body interaction, the method is demonstrated by considering a toy model of three spinless particles interacting with pair-wise \ensuremathδ-function potentials. The numerical results are compared with the exact solutions of three spinless bosons interaction when the strength of short-range interactions are set equal for all pairs.
We explore a variational approach to the finite-volume N-body problem. The general formalism for N nonrelativistic spinless particles interacting with periodic pairwise potentials yields N-body secular equations. The solutions depend on the infinite-volume N-body wave functions. Given that the infinite-volume N-body dynamics may be solved by the standard Faddeev approach, the variational N-body formalism can provide a convenient numerical framework for finding discrete energy spectra in periodic lattice structures.
We extend a method previously presented for the three-body problem with attractive interparticle potentials to the case of repulsive interparticle potentials, with periodic boundary conditions on each particle. As before, we decompose the wave function into three parts, and from the Schrödinger equation write equations for these parts. We observe that, in one-dimension and for δ-function potentials,these equations can be easily solved numerically, and we present these solutions for potential“strengths” ranging from zero to infinity. We then discuss a more general method for solving these equations, which method involves expansions in a certain set of two-body functions. This general method is not necessarily limited to either one dimension or δ-function potentials, but, as a check on it, we do apply it to that case and get good agreement with the previous numerical results. As a further exploration of the method, we apply it to square well potentials in one dimension. In an Appendix we discuss the set of two-body functions that we use.
The low-energy scattering of three bosons or distinguishable particles with short-range interactions is characterized by a fundamental parameter, the three-body scattering hypervolume, which is responsible for the nonuniversal effects in dilute Bose-Einstein condensates. We derive an analytical formula of the three-body scattering hypervolumes for weak interactions. When the interaction supports two-body bound states, the three-body scattering hypervolume gains a negative imaginary part, which is directly related to the rate constant for three-body recombinations. We develop a numerical method to calculate the three-body scattering hypervolumes for some model potentials with variable strengths. For real atoms with van der Waals potential, the three-body scattering hypervolumes are much harder to compute, because the three-body wave function is highly oscillatory at smaller inter-atomic distances. However, they may be extractable from precision data such as the collective-mode frequencies of trapped Bose-Einstein condensates. In many numerical simulations, the system being simulated is put into a large but finite volume, such as a large periodic box of size L. To extract the low-energy scattering properties of two particles in infinite space from such simulations, L\\"uscher's formula must be used. In the second part of this thesis, we generalize L\\"uscher's formula to d spatial dimensions. We obtain its s-wave and p-wave approximations and the systematic expansions of the energies of the low-lying states. At a s-wave resonance in dge4 dimensions, we identify two low-lying states close to the threshold with nearly opposite energies, Esimpm1/Ld/2 when dge5, or Esimpm1/L2sqrtlnL when d=4. These calculations provide important insights to the physics of three particles at a three-body resonance in finite volumes in two or three dimensions. Three-body resonances are important phenomena not yet completely explored in the experiments on cold atoms.
A simple model of a boson system is considered for which the exact Schr\"odinger equation reduces to a difference equation which can be solved numerically. Application of the Bogoliubov method and comparison of its results with the numerical results then yield some insight into the validity of the Bogoliubov approximations. The simple model considered has, in the absence of interactions, only a zero energy state and two states of unit energy available to each boson. Initially there is assumed to be no interaction between bosons in the degenerate excited states and all existing interactions are assumed to be repulsive. The ground-state energy is calculated in the Bogoliubov approximation, with and without the depletion effect considered, and in another numerical approximation. These results and the results of the numerical solution of the exact Schr\"odinger equation converge to the Bogoliubov result with depletion ignored when the number of particles in the system approaches infinity. It is surprising to note that, for the interaction strengths considered, the Bogoliubov result is within 3% of the numerical result for as few as 32 particles in the system, and within 10% for only 4 particles. A modification of the above system is considered in which there exists an additional two-body interaction between particles in the degenerate excited states which may be attractive or repulsive. It is shown that the ground-state energy, with this additional interaction present, deviates from the Bogoliubov value linearly with the strength of the added interaction, i.e., $W={W}_{\mathrm{Bog}}+\ensuremath{\alpha}F$ (with F the strength parameter of the interaction). The factor $\ensuremath{\alpha}$ is found to be small and constant over a wide range of F and in this range the deviation of W from the Bogoliubov value is small. There is, however, a transition point FT, beyond which $\ensuremath{\alpha}$ is constant and large and the usual Bogoliubov approximation is invalid for interactions more attractive than that characterized by FT. A Bogoliubov-like approximation is shown to be quite accurate in this region. It is then shown that $F_{T}^{\ensuremath{-}1}$ is a linear function of the number of particles in the system.
We consider the evolution of N bosons, where N is large, with two-body interactions of the form N3βv(Nβ⋅), 0≤β≤1. The parameter β measures the strength of interactions. We compare the exact evolution with an approximation which considers the evolution of a mean field coupled with an appropriate description of pair excitations, see [25, 26]. For 0≤β<1/2, we derive an error bound of the form p(t)/Nα, where α>0 and p(t) is a polynomial, which implies a specific rate of convergence as N→∞.
We consider either 3 spinless bosons or 3 equal mass spin-1/2 fermions, interacting via a short-range potential of infinite scattering length and trapped in an isotropic harmonic potential. For a zero-range model, we obtain analytically the exact spectrum and eigenfunctions: for fermions all the states are universal; for bosons there is a coexistence of decoupled universal and efimovian states. All the universal states, even the bosonic ones, have a tiny 3-body loss rate. For a finite range model, we numerically find for bosons a coupling between zero angular momentum universal and efimovian states; the coupling is so weak that, for realistic values of the interaction range, these bosonic universal states remain long-lived and observable.