  We adopt the heavy baryon chiral perturbation theory (HBChPT) to calculate
the scattering lengths of $\phi B_{cc}^{(*)}$ up to $\mathcal{O}(p^3)$, where
$\phi$ is the pseudoscalar mesons. The recoil effect and the mass splitting
between the spin-$1\over 2$ and spin-$3\over 2$ doubly charmed baryons are
included. In order to give the numerical results, we construct the chiral
Lagrangians with heavy diquark-antiquark (HDA) symmetry in a formally covariant
approach. Then, we relate the low energy constants (LECs) of the doubly charmed
baryons to those of $D^{(*)}$ mesons. The LECs for the $\phi D^{(*)}$
scattering are estimated in two scenarios, fitting lattice QCD results and
using the resonance saturation model. The chiral convergence of the first
scenario is not good enough due to the the large strange quark mass and the
presence of the possible bound states, virtual states and resonance. The final
results for two scenarios are consistent with each other. The interaction for
the $[\pi\Xi^{(*)}_{cc}]^{(1/2)}$, $[K\Xi^{(*)}_{cc}]^{(0)}$,
$[K\Omega^{(*)}_{cc}]^{(1/2)}$, $[\eta\Xi^{(*)}_{cc}]^{(1/2)}$,
$[\eta\Omega^{(*)}_{cc}]^{(0)}$ and $[\bar{K}\Xi^{(*)}_{cc}]^{(0)}$ channels
are attractive. The most attractive channel $[\bar{K}\Xi^{(*)}_{cc}]^{(0)}$ may
help to form the partner states of the $D_{s0}^*(2317)$ ($D_{s1}(2460)$) in the
doubly heavy sector.
