Hydrogenlike molecules composed of $D_1D_1$, $D_1D^*_2$ and $D^*_2D^*_2$

We systematically explore the S-wave $D_1D_1$, $D_1D^*_2$ and $D^*_2D^*_2$ states with various isospin-spin-orbit ($ISL$) configurations in the quark model. We propose nine stable dimeson states with the $ISL$ configurations, $ISL=001$, $010$, $012$, $100$, $102$, $110$, $112$, $120$, and $122$, against dissociation into their constituent mesons. Those bound states are hydrogenlike molecular states, where the two subclusters are moderately overlapped and the QCD covalent bond is formed due to the delocalization of light quarks. The QCD covalent bond serves as the primary binding mechanism in the bound states with $I=1$. However, the exchange of $\pi$ and $\sigma$-meson plays a pivotal role in the bound states with $I=0$. The coupled-channel effect is essential in the formation of the bound states with $ISL=001$, $010$, $012$, $100$, and $102$.