We present a systematic investigation of the possible molecular states composed of a pair of doubly charmed baryons ($\Xi_{cc}\Xi_{cc}$) or one doubly charmed baryon and one doubly charmed antibaryon $(\Xi_{cc}\bar{\Xi}_{cc})$ within the framework of the one-boson-exchange-potential model. For the spin-triplet systems, we take into account the mixing between the ${}^3S_1$ and ${}^3D_1$ channels. For the baryon-baryon system $\Xi_{cc}\Xi_{cc}$ with $(R,I) = (\bar{3}, 1/2)$ and $(\bar{3}, 0)$, where $R$ and $I$ represent the group representation and the isospin of the system, respectively, there exist loosely bound molecular states. For the baryon-antibaryon system $\Xi_{cc}\bar{\Xi}_{cc}$ with $(R,I) = (8, 1)$, $(8, 1/2)$ and $(8,0)$, there also exist deuteron-like molecules. The $B_{cc}\bar{B}_{cc}$ molecular states may be produced at LHC. The proximity of their masses to the threshold of two doubly charmed baryons provides a clean clue to identify them.