<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Hello</p>
<p>there is a separate residue for charged GLN (QLN) in the OPLS-forcefield file aminoacids.rtp.</p>
<p>This residue uses the preprocessor define &quot;<span>dih_GLN_chi4_C_C_O_H</span>&quot; for a dihedral parameter set.</p>
<p>However, this macro is not defined anywhere (esp. not in ffbonded.itp with all the others).</p>
<p>The reason why grompp does not complain, is because the function push_bond() uses sscanf</p>
<p><span>&nbsp;&nbsp; nread = sscanf(line,format,&amp;cc[0],</span>.....</p>
<p>expecting double values. When strings are read instead, the variable nread is zero, and grompp assumes</p>
<p>default values for this dihedral.</p>
<p>If I'm correct, we should either delete the non-existing macro, or have to add values for this dihedral.</p>
<p>Thank you for your time,</p>
<p>regards Martin<br>
</p>
</div>
</body>
</html>