<div dir="ltr"><div><br></div><div>Hello all,</div><div>I noticed that torsions evaluated by CP2K when using AMBER parameters are exactly half of the correct answer. Here is the description below.</div><div><br></div><div><span style="font-size:13px">I do single point energies of a small dipeptide in vacuum using amber forcefield. (attached input files)</span><br></div><div><span style="font-size:13px">Gromacs and OPENMM give the same intramolecular energy components.</span></div><div>_____________________________________________</div><div><div>However CP2K does not:</div><div><span style="font-size:13px">more 1.out | grep -A 2 "FIST energy"</span></div><div><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span style="font-size:13px"> FIST energy contributions in kcal/mol:</span><br></p><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span> BOND    =       17.8009  ANGLE   =       34.4220  UBRAD   =        0.0000</span></p><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span> TORSION =       10.0356  IMPTORS =        0.0000  OPBEND  =        0.0000</span></p></div></div><div><div>BOND = 17.8009* 4.184 = 74.478 kJ/mol</div><div>ANGLE = 34.4220 <span style="font-size:13px">* 4.184 = 144.021 kJ/mol</span></div><div>TORSION = 10.0356 <span style="font-size:13px">* 4.184 = <b><u><font color="#ff0000">41.988 kJ/mol</font></u></b></span></div></div><div><span style="font-size:13px">_____________________________________________________</span></div><div><br></div><div><font size="2">GROMACS (after converting amber.prmtop -> </font><font size="2">Gromacs topology)</font></div><div><span style="font-size:13px">BOND = </span><span style="font-size:13px">74.478 kJ/mol</span></div><div><span style="font-size:13px">ANGLE = </span><span style="font-size:13px">144.021 kJ/mol</span></div><div><span style="font-size:13px">TOTAL TORSIONS =  80.691(proper) + 3.2858 (improper) = <b><u><font color="#ff0000">83.3549 kJ/mol</font> </u></b></span></div><div><span style="font-size:13px"><b> </b></span><span style="font-size:13px"><b>(Exactly double as that from CP2K!) </b></span></div><div><span style="font-size: 13px;">TOTAL PE = ~251 kJ/mol. (This matches exactly with OPENMM when using amber.prmtop )</span></div><div>So Gromacs and openmm results agree with each other. </div><div>______________________________________________________</div><div><span style="font-size:13px"><br></span></div><div>Then saw this post: </div><div><u><a href="https://groups.google.com/forum/#!searchin/cp2k/amber$20torsion%7Csort:date/cp2k/PeMpumD2oXU/1uWBed2IfhkJ" target="_blank">https://groups.google.com/<wbr>forum/#!searchin/cp2k/amber$<wbr>20torsion%7Csort:date/cp2k/<wbr>PeMpumD2oXU/1uWBed2IfhkJ</a></u><br></div><div><br></div><div><div><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span>mpirun -n 1 cp2k.popt --version</span></p><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span> CP2K version 6.1</span></p><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span> SVN source code revision svn:18464</span></p><p style="font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(242,242,242);background-color:rgb(0,0,0)"><span> cp2kflags: libint fftw3 libxc parallel mpi3 scalapack xsmm libderiv_max_am1=5 libint_max_am=6 plumed2 mkl</span></p></div></div><div><span><br></span></div><div><span style="font-size:13px"><b>(Seems like TORSIONS is still broken in svn:18464)</b></span></div><div><br></div></div>