Using the Monte Carlo event generator tools Pythia and Herwig, we simulate the production of bottom/charm meson and antimeson pairs at hadron colliders in proton-proton/antiproton collisions. With these results, we derive an order-of-magnitude estimate for the production rates of the bottom analogs and the spin partner of the X (3872) as hadronic molecules at the LHC and Tevatron experiments. We find that the cross sections for these processes are at the nb level, so that the current and future data sets from the Tevatron and LHC experiments offer a significant discovery potential. We further point out that the X b / X b2 should be reconstructed in the γ ϒ(nS)(n = 1, 2, 3), ϒ(1S)π + π − π 0 , or χ bJ π + π − instead of the ϒ(nS)π + π − final states.