In present work, we study an 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 $delta$-function potentials. The numerical results are compared with the exact solutions of three spinless bosons interaction when strength of short-range interactions are set equal for all pairs.
In this work, based on consideration of periodicity and asymptotic forms of wave function, we propose a novel approach to the solution of finite volume three-body problem by mapping a three-body problem into a higher dimensional two-body problem. The idea is demonstrated by an example of two light spinless particles and one heavy particle scattering in one spatial dimension. This 1D three-body problem resembles a two-body problem in two spatial dimensions mathematically, and quantization condition of 1D three-body problem is thus derived accordingly.
The volume-dependence of a shallow three-particle bound state in the cubic box with a size $L$ is studied. It is shown that, in the unitary limit, the energy-level shift from the infinite-volume position is given by $Delta E=c (kappa^2/m),(kappa L)^{-3/2}|A|^2 exp(-2kappa L/sqrt{3})$, where $kappa$ is the bound-state momentum and $|A|^2$ denotes the three-body analog of the asymptotic normalization constant, which encodes the information about the short-range interactions in the three-body system.
Using the general formalism presented in Refs. [1,2], we study the finite-volume effects for the $mathbf{2}+mathcal{J}tomathbf{2}$ matrix element of an external current coupled to a two-particle state of identical scalars with perturbative interactions. Working in a finite cubic volume with periodicity $L$, we derive a $1/L$ expansion of the matrix element through $mathcal O(1/L^5)$ and find that it is governed by two universal current-dependent parameters, the scalar charge and the threshold two-particle form factor. We confirm the result through a numerical study of the general formalism and additionally through an independent perturbative calculation. We further demonstrate a consistency with the Feynman-Hellmann theorem, which can be used to relate the $1/L$ expansions of the ground-state energy and matrix element. The latter gives a simple insight into why the leading volume corrections to the matrix element have the same scaling as those in the energy, $1/L^3$, in contradiction to earlier work, which found a $1/L^2$ contribution to the matrix element. We show here that such a term arises at intermediate stages in the perturbative calculation, but cancels in the final result.
We discuss signatures of bound-state formation in finite volume via the Luscher finite size method. Assuming that the phase-shift formula in this method inherits all aspects of the quantum scattering theory, we may expect that the bound-state formation induces the sign of the scattering length to be changed. If it were true, this fact provides us a distinctive identification of a shallow bound state even in finite volume through determination of whether the second lowest energy state appears just above the threshold. We also consider the bound-state pole condition in finite volume, based on Luschers phase-shift formula and then find that the condition is fulfilled only in the infinite volume limit, but its modification by finite size corrections is exponentially suppressed by the spatial lattice size L. These theoretical considerations are also numerically checked through lattice simulations to calculate the positronium spectrum in compact scalar QED, where the short-range interaction between an electron and a positron is realized in the Higgs phase.
Using the framework of non-relativistic effective field theory, the finite-volume ground-state energy shift is calculated up-to-and-including $O(L^{-6})$ for the system of three pions in the channel with the total isospin $I=1$. The relativistic corrections are included perturbatively, up to the same order in the inverse of the box size $L$. The obtained explicit expression, together with the known result for the system with maximal isospin $I=3$, can be used for the extraction of two independent effective three-body couplings from the measured ground-state spectrum of three pions.