  Intensity mapping of 21cm emission from neutral hydrogen (HI) promises to be
a powerful probe of large-scale structure in the post-reionisation epoch.
However, HI intensity mapping (IM) experiments will suffer the loss of
long-wavelength line-of-sight HI modes in the galactic foreground subtraction
process. The loss of these modes is particularly problematic for detecting HI
IM cross-correlations with projected large-scale structure tracers, such as CMB
secondary anisotropies. Here, we propose a cross-bispectrum estimator to
recover the cross-correlation of the HI IM field, $\delta T_{21},$ with the CMB
lensing field, $\kappa,$ constructed by correlating the position-dependent HI
power spectrum with the mean overdensity traced by CMB lensing.
  We study the cross-bispectrum estimator, $B^{\bar \kappa \delta T_{21} \delta
T_{21}},$ in the squeezed limit and forecast its detectability based on HI IM
measurements from HIRAX and CMB lensing measurements from AdvACT. The
cross-bispectrum improves constraints on cosmological parameters; in
particular, the constraint on the dark energy equation-of-state parameter,
$w_0,$ improves on the HI IM auto-power spectra constraint by 44\% (to 0.014),
while the constraint on $w_a$ improves by 33\% (to 0.08), assuming Planck
priors in each case. These results are robust to HI IM foreground removal
because they largely derive from small-scale HI modes. The HI-HI-$\kappa$
cross-bispectrum thus provides a novel way to recover HI correlations with CMB
lensing and constrain cosmological parameters at a level that is competitive
with next-generation galaxy redshift surveys. As a striking example of this, we
find that the combined constraint on the sum of the neutrino masses, while
varying all redshift and standard cosmological parameters within a
$w_0w_a\Omega_K$CDM model, is 5.5 meV.
