Energy scan of the $e^+e^- \to h_b(nP)\pi^+\pi^-$ $(n=1,2)$ cross sections and evidence for $\Upsilon(11020)$ decays into charged bottomonium-like states

Using data collected with the Belle detector at the KEKB asymmetric-energy $e^+e^-$ collider, we measure the energy dependence of the $e^+e^- \to h_b(nP)\pi^+\pi^-$ $(n=1,2)$ cross sections from thresholds up to $11.02\,$GeV. We find clear $\Upsilon(10860)$ and $\Upsilon(11020)$ peaks with little or no continuum contribution. We study the resonant substructure of the $\Upsilon(11020) \to h_b(nP)\pi^+\pi^-$ transitions and find evidence that they proceed entirely via the intermediate isovector states $Z_b(10610)$ and $Z_b(10650)$. The relative fraction of these states is loosely constrained by the current data: the hypothesis that only $Z_b(10610)$ is produced is excluded at the level of 3.3 standard deviations, while the hypothesis that only $Z_b(10650)$ is produced is not excluded at a significant level.