  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.
