  In this paper we study $I = 0$ bottomonium in $S$, $P$, $D$ and $F$ waves
considering five coupled channels, one confined quarkonium and four open
$B^{(*)} \bar B^{(*)}$ and $B_s^{(*)} \bar B_s^{(*)}$ meson-meson channels. To
this end we use and extend a recently developed novel approach utilizing
lattice QCD string breaking potentials for the study of quarkonium bound states
and resonances. This approach is based on the Born Oppenheimer approximation
and the unitary emergent wave method and allows to compute the poles of the
$\mbox{T}$ matrix. We compare our results to existing experimental results for
$I = 0$ bottomonium and discuss masses, decay widths and the assignment of
angular momentum quantum numbers. Moreover, we determine the quarkonium and
meson-meson composition of these states to clarify, which of them are ordinary
quarkonium, and which of them should rather be interpreted as tetraquarks.
