We discuss a comprehensive study of bottomonium bound states and resonances based on static-static-light-light potentials. We use the Born-Oppenheimer diabatic approximation [1]}, which is a two-step approach. First, static quark-antiquark potentials in presence of a light quark-antiquark pair are computed with lattice QCD (here we use existing results from Ref. [2]}). Then, in a second step, these potentials are used in a coupled channel Schrödinger equation. This approach was successfully applied to \(\bar{b}\bar{b}ud\)  systems (see e.g. Refs. [3]}, [4]}) and to \(I = 0\)  bottomonium in an \(S\)  wave [5]}, [6]}. There are also ongoing efforts to study \(I = 1\)  bottomonium in a similar way [7]}.
