<div dir="ltr">Dear CP2K experts:<div><br></div><div>I am running a system with QMMM method in CP2K. By checking the gaussian electrostatic coupling reference (http://pubs.acs.org/doi/abs/10.1021/ct050123f), I am setting the following:</div><div><br></div><div><div>    ECOUPL GAUSS</div><div>    NOCOMPATIBILITY .TRUE.</div><div>    USE_GEEP_LIB 10</div><div>    CENTER NEVER</div></div><div><br></div><div>and I am varying USE_GEEP_LIB with the values of 2, 5, 10 and 12. The problems I have are as the following:</div><div><br></div><div>1. When USE_GEEP_LIB is set to 10, the simulation is fastest. However, I expect the fewer number of levels of grid I use, the faster the simulation would be.</div><div><br></div><div>2. Even if the simulation got accelerated, it is still not enough. The fastest simulation gave me only 4 times speed up compared with pure ab initio MD with DFT for all the atoms.</div><div><br></div><div>3. In the reference, it also mentioned that the cutoff distance would affect the acceleration, while I did not find the corresponding set up in the manual.</div><div><br></div><div>I am attaching all my input files. I really want to know how could I set up the input files so that I could get faster simulations? (4 times faster is not good enough for me obviously).</div><div><br></div><div>Thanks very much for your help!</div><div><br></div><div>Tao</div></div>