Effective renormalized multi-body interactions of harmonically confined ultracold neutral bosons

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 pseudopotential. 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[(a\l)^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.