We study the spectrum of the family of Hamiltonian
$H_N = \sum_{j=1}^N \left[ -\Delta_{x_j} + V_\text{ext}(x_j) \right] + \sum_{1\leq i<j \leq N} N^2 V(N(x_i-x_j))$,
where $H_N$ acts on $L^2(\mathbb{R}^{3N})$. This family of Hamiltonians models a gas of bosons confined by an external potential $ V_\text{ext}$. These quantum-mechanical systems are interesting
as they exhibit a special state of matter, the so-called Bose-Einstein condensate. In
such Bose-Einstein condensates almost all particles can be described by the same one-
particle wave function. We will see how to make this notion mathematically precise and what the latest rigorous results are.
We show that in this regime low-energy states exhibit complete Bose-Einstein condensation with optimal rate of convergence. Building upon this condensation result we determine the next order correction of the asymptotics of the ground state energy for large number of particles and also the low-energy excitation spectrum.