I already have the free energy TOTEN of the structure after optimization with GGA-PBE and now I want to calculate vdW corrections. I modified the INCAR file with
IVDW = 11 [zero damping DFT-D3 method of Grimme (available as of VASP.5.3.4)]
or
IVDW = 12 [DFT-D3 method with Becke-Jonson damping (available as of VASP.5.3.4)]
Made use of CONTCAR, KPOINTS, POTCAR and executable from the previous calculation.
However, I realized that its not working at all. After running for several hours, still there is not much of change and hence no convergence.
Can someone please suggest the steps for calculating vdW using DFT-D3.
Thank you.