Dear Hilal Kucuk, van der Waals correction can be switched on in VASP with the IVDW tag in INCAR, for example IVDW=2 (Tkatchenko-Scheffler method) or IVDW=22 (Many-body dispersion energy method), see the Vaspwiki for more information. Steric repulsion is already part of DFT-calculated forces in any case, so is any interaction due to static charge densities. So the basic components of hydrogen bonding should be included in DFT + vdW correction. But I am not sure whether the accuracy you cat get with that would be sufficient of if you need to use some beyond DFT method.

Thanks a lot Dear Martin. Can I find any tool to calculate them?

You can get the vdW correction directly from VASP. If you need it as a separate component of the energy, you can find it in OUTCAR as an entry titled "Estimated vdW energy (eV)". The final forces and TOTEN (final total energy) following after the vdW entry should include the vdW component then. Of course, you need to do the calculation with the IVDW tag set (to something nonzero) in INCAR. The steric repulsion and electrostatic interaction are also accounted for in the total energy, but if you need to know their contributions as separate components, that may be tricky. There is a kind of decomposition of TOTEN into several components in OUTCAR, namely components called "alpha Z" or PSCENC, "Ewald energy" or TEWEN, "-1/2 Hartree" or DENC, "-exchange" or EXHF, "-V(xc)+E(xc)" or XCENC, "PAW double counting", "eigenvalues" or EBANDS, and "atomic energy" or EATOM, but honestly, I do not know how to get anything useful from that.

Dear Martin,

Sorry for the late answer and thanks a lot for the answer. It is quite helpful.

Best wishes