Critical temperature of trapped interacting bosons from large-N -based theories

Ultracold atoms provide clues to an important many-body problem regarding the dependence of the Bose-Einstein condensation transition temperature T_c on interactions. However, cold atoms are trapped in harmonic potentials and theoretical evaluations of the T_c shift of trapped interacting Bose gases are challenging. While previous predictions of the leading-order shift have been confirmed, more recent experiments exhibit higher-order corrections beyond available mean-field theories. By implementing two large-N -based theories with the local density approximation (LDA), we extract next-order corrections of the T_c shift. The leading-order large-N theory produces results quantitatively different from the latest experimental data. The leading-order auxiliary-field (LOAF) theory, containing both normal and anomalous density fields, captures the T_c shift accurately in the weak-interaction regime. However, the LOAF theory shows incompatible behavior with the LDA, and forcing the LDA leads to density discontinuities in the trap profiles. We present a phenomenological model based on the LOAF theory, which repairs the incompatibility and provides a prediction of the T_c shift in the stronger-interaction regime.