We report on a technique to determine the van der Waals coefficients of lithium ͑Li͒ atoms based on relativistic coupled-cluster theory. These quantities are determined using the imaginary parts of the scalar dipole and quadrupole polarizabilities, which are evaluated using an approach that we have proposed earlier ͓B. K. Sahoo, Chem. Phys. Lett. 448, 144 ͑2007͔͒. Our procedure is fully ab initio, and avoids the sum-over-thestates approach. We present the dipole and quadrupole polarizabilities of many of the low-lying excited states of Li. Also, the off-diagonal dipole and quadrupole polarizabilites between some of the low-lying states of Li have been calculated.