<div dir="ltr">Dear all,<br><br>I am calculating the binding energy of graphene bilayer of the size 3nmX3nm (2X416 atoms). It turns out that the result from calculations with classical force field is quite different from that of DFTs:<br><br>CHARMM FF:                                       -0.579 Hartree<br>DFT-D3:                                                -1.056 Hartree<br>DFT-D3 with BSSE correction:           -0.968 Hartree<br><br>It is quite strange that the value of CHARMM is nearly one half of DFT results. <br><br>The binding energy is defined here as: dE = E_bilayer - 2*E_single_layer.<br><br>All my input files are attached in the tar file. Please help me to find out where the problem lies in. Thank you!<br><br>Best wishes,<br>Robert<br></div>