When running AIMD, I usually add DFT-D3(BJ) dispersion correction to account for the lack of dispersion when running at the level of BLYP/TZVP-MOLOPT-GTH.<div><br></div><div><div>&VDW_POTENTIAL</div><div>        POTENTIAL_TYPE PAIR_POTENTIAL</div><div>        &PAIR_POTENTIAL</div><div>          PARAMETER_FILE_NAME dftd3.dat</div><div>          TYPE DFTD3(BJ)</div><div>          REFERENCE_FUNCTIONAL BLYP</div><div>        &END PAIR_POTENTIAL</div><div>      &END VDW_POTENTIAL</div></div><div><br></div><div>However, I wonder if I would still need DFT-D3(BJ) during QM/MM calculation? The dispersion between the QM region and the MM region or between the atoms in the MM region is taken into account by Lennard-Jones potential.</div><div>For cp2k, is the dispersion between the atoms in the QM region taken into account by the Lennard-Jones potential, or one would need to add dispersion correction via DFT-D3(BJ)?</div><div>Thank you.</div>