We calculate the renormalized effective 2-, 3-, and 4-body interactions for N neutral ultracold bosons in the ground state of an isotropic harmonic trap, assuming 2-body interactions modeled with the combination of a zero-range and energy-dependent p
seudopotential. We work to third-order in the scattering length a defined at zero collision energy, which is necessary to obtain both the leading-order effective 4-body interaction and consistently include finite-range corrections for realistic 2-body interactions. The leading-order, effective 3- and 4-body interaction energies are U3 = -(0.85576...)(a/l)^2 + 2.7921(1)(a/l)^3 + O[(a/l)^4] and U4 = +(2.43317...)(a/l)^3 + O[(al)^4], where w and l are the harmonic oscillator frequency and length, respectively, and energies are in units of hbar*w. The one-standard deviation error 0.0001 for the third-order coefficient in U3 is due to numerical uncertainty in estimating a slowly converging sum; the other two coefficients are either analytically or numerically exact. The effective 3- and 4-body interactions can play an important role in the dynamics of tightly confined and strongly correlated systems. We also performed numerical simulations for a finite-range boson-boson potential, and it was comparison to the zero-range predictions which revealed that finite-range effects must be taken into account for a realistic third-order treatment. In particular, we show that the energy-dependent pseudopotential accurately captures, through third order, the finite-range physics, and in combination with the multi-body effective interactions gives excellent agreement with the numerical simulations, validating our theoretical analysis and predictions.
We show that there are effective three- and higher-body interactions generated by the two-body collisions of atoms confined in the lowest vibrational states of a 3D optical lattice. The collapse and revival dynamics of approximate coherent states loa
ded into a lattice are a particularly sensitive probe of these higher-body interactions; the visibility of interference fringes depend on both two-, three-, and higher-body energy scales, and these produce an initial dephasing that can help explain the surprisingly rapid decay of revivals seen in experiments. If inhomogeneities in the lattice system are sufficiently reduced, longer timescale partial and nearly full revivals will be visible. Using Feshbach resonances or control of the lattice potential it is possible to tune the effective higher-body interactions and simulate effective field theories in optical lattices.
