jats:titleAjats:scbstract</jats:sc> </jats:title>jats:pWe explore the spectrum of excited and exotic bottomonia using lattice QCD. Highly excited states are identified with masses up to 11,000 MeV, many of which can be grouped into supermultiplets matching those of the quark model while exotic spin-parity-charge-conjugation quantum numbers jats:italicJ</jats:italic>jats:supPC</jats:sup> = 0jats:sup+jats:italic−</jats:italic></jats:sup>jats:italic,</jats:italic> 1jats:supjats:italic−</jats:italic>+</jats:sup>jats:italic,</jats:italic> 2jats:sup+jats:italic−</jats:italic></jats:sup> that cannot be formed from jats:inline-formulajats:alternativesjats:tex-math$$ \overline{q}q $$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:mover mml:miq</mml:mi> mml:mo¯</mml:mo> </mml:mover> mml:miq</mml:mi> </mml:math></jats:alternatives></jats:inline-formula> alone are also identified. Single-meson operator constructions are used that have good jats:italicJ</jats:italic>jats:supPC</jats:sup> in the continuum, these are found to overlap well onto heavy quark states up to jats:italicJ</jats:italic> ≤ 4. A continuum jats:italicJ</jats:italic>jats:supPC</jats:sup> is assigned to each level, based on the distribution amongst lattice irreps and dominant operator overlaps. States with a dominant gluonic component are identified and form a hybrid supermultiplet with jats:italicJ</jats:italic>jats:supPC</jats:sup> = (0jats:italic,</jats:italic> 1jats:italic,</jats:italic> 2)jats:supjats:italic−</jats:italic>+</jats:sup>jats:italic,</jats:italic> 1jats:supjats:italic−−</jats:italic></jats:sup>, approximately 1500 MeV above the ground-state jats:italicη</jats:italic>jats:subjats:italicb</jats:italic></jats:sub>, similar to previous computations with light, strange and charm quark systems.</jats:p>