Hi Nitol,
I'm not sure how you wanted to use VarD to integrate by parts and discard total derivtives but for that I use a simple rule like the following:
IBPrule1 = XX_*h[LI[1], aa_, bb_] CD[cc_]@CD[dd_]@h[LI[1], ee_, ff_] -> -CD[cc]@(XX*h[LI[1], aa, bb])*CD[dd]@h[LI[1], ee, ff]
or, more generally,
IBPrule2 = XX_*CD[cc_]@CD[dd_]@YY_ -> -CD[cc]@XX*CD[dd]@YY
Note that the indices are different (you use greek indices in your notebook so you'll have to change that)
Best,
Sergi