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 determine the energetically lowest lying states in the BEC-BCS crossover regime of s-wave interacting two-component Fermi gases under harmonic confinement by solving the many-body Schrodinger equation using two distinct approaches. Essentially exa
ct basis set expansion techniques are applied to determine the energy spectrum of systems with N=4 fermions. Fixed-node diffusion Monte Carlo methods are applied to systems with up to N=20 fermions, and a discussion of different guiding functions used in the Monte Carlo approach to impose the proper symmetry of the fermionic system is presented. The energies are calculated as a function of the s-wave scattering length a_s for N=2-20 fermions and different mass ratios kappa of the two species. On the BEC and BCS sides, our energies agree with analytically-determined first-order correction terms. We extract the scattering length and the effective range of the dimer-dimer system up to kappa = 20. Our energies for the strongly-interacting trapped system in the unitarity regime show no shell structure, and are well described by a simple expression, whose functional form can be derived using the local density approximation, with one or two parameters. The universal parameter xi for the trapped system for various kappa is determined, and comparisons with results for the homogeneous system are presented.
