<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns:p="urn:schemas-microsoft-com:office:powerpoint" xmlns:a="urn:schemas-microsoft-com:office:access" xmlns:dt="uuid:C2F41010-65B3-11d1-A29F-00AA00C14882" xmlns:s="uuid:BDC6E3F0-6DA3-11d1-A2A3-00AA00C14882" xmlns:rs="urn:schemas-microsoft-com:rowset" xmlns:z="#RowsetSchema" xmlns:b="urn:schemas-microsoft-com:office:publisher" xmlns:ss="urn:schemas-microsoft-com:office:spreadsheet" xmlns:c="urn:schemas-microsoft-com:office:component:spreadsheet" xmlns:odc="urn:schemas-microsoft-com:office:odc" xmlns:oa="urn:schemas-microsoft-com:office:activation" xmlns:html="http://www.w3.org/TR/REC-html40" xmlns:q="http://schemas.xmlsoap.org/soap/envelope/" xmlns:rtc="http://microsoft.com/officenet/conferencing" xmlns:D="DAV:" xmlns:Repl="http://schemas.microsoft.com/repl/" xmlns:mt="http://schemas.microsoft.com/sharepoint/soap/meetings/" xmlns:x2="http://schemas.microsoft.com/office/excel/2003/xml" xmlns:ppda="http://www.passport.com/NameSpace.xsd" xmlns:ois="http://schemas.microsoft.com/sharepoint/soap/ois/" xmlns:dir="http://schemas.microsoft.com/sharepoint/soap/directory/" xmlns:ds="http://www.w3.org/2000/09/xmldsig#" xmlns:dsp="http://schemas.microsoft.com/sharepoint/dsp" xmlns:udc="http://schemas.microsoft.com/data/udc" xmlns:xsd="http://www.w3.org/2001/XMLSchema" xmlns:sub="http://schemas.microsoft.com/sharepoint/soap/2002/1/alerts/" xmlns:ec="http://www.w3.org/2001/04/xmlenc#" xmlns:sp="http://schemas.microsoft.com/sharepoint/" xmlns:sps="http://schemas.microsoft.com/sharepoint/soap/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:udcs="http://schemas.microsoft.com/data/udc/soap" xmlns:udcxf="http://schemas.microsoft.com/data/udc/xmlfile" xmlns:udcp2p="http://schemas.microsoft.com/data/udc/parttopart" xmlns:wf="http://schemas.microsoft.com/sharepoint/soap/workflow/" xmlns:dsss="http://schemas.microsoft.com/office/2006/digsig-setup" xmlns:dssi="http://schemas.microsoft.com/office/2006/digsig" xmlns:mdssi="http://schemas.openxmlformats.org/package/2006/digital-signature" xmlns:mver="http://schemas.openxmlformats.org/markup-compatibility/2006" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns:mrels="http://schemas.openxmlformats.org/package/2006/relationships" xmlns:spwp="http://microsoft.com/sharepoint/webpartpages" xmlns:ex12t="http://schemas.microsoft.com/exchange/services/2006/types" xmlns:ex12m="http://schemas.microsoft.com/exchange/services/2006/messages" xmlns:pptsl="http://schemas.microsoft.com/sharepoint/soap/SlideLibrary/" xmlns:spsl="http://microsoft.com/webservices/SharePointPortalServer/PublishedLinksService" xmlns:Z="urn:schemas-microsoft-com:" xmlns:st="" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.25in 1.0in 1.25in;}
div.WordSection1
        {page:WordSection1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=WordSection1>
<p class=MsoNormal>Dear GROMACS users,<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>I am encountering a couple of issues when trying to perform
normal mode analysis in an implicit solvent (GBSA) setting.<o:p></o:p></p>
<p class=MsoNormal>I am using version 4.5.1 with double precision;
unfortunately I do not have the single precision version installed for comparison.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The starting point is a small protein which was energy minimized
in two stages, also in double precision and with GBSA, producing “em2”
files. The mdp file used for the normal mode analysis is the following:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>----------------- nm.mdp -------------------------<o:p></o:p></p>
<p class=MsoNormal>integrator = nm<o:p></o:p></p>
<p class=MsoNormal>nsteps
= 1<o:p></o:p></p>
<p class=MsoNormal>implicit_solvent = GBSA<o:p></o:p></p>
<p class=MsoNormal>gb_algorithm = Still ; the default<o:p></o:p></p>
<p class=MsoNormal>rgbradii =
1.0 ; must be equal to rlist <o:p></o:p></p>
<p class=MsoNormal>rlist
= 1.0<o:p></o:p></p>
<p class=MsoNormal>coulombtype = cut-off <o:p></o:p></p>
<p class=MsoNormal>rcoulomb =
1.0<o:p></o:p></p>
<p class=MsoNormal>vdwtype
= cut-off<o:p></o:p></p>
<p class=MsoNormal>rvdw
= 1.0<o:p></o:p></p>
<p class=MsoNormal>------------------------------------------------------- <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The command used is <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>grompp_d -f nm.mdp -p topol.top -c em2.gro -t em2.trr -o
nm.tpr<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>which ran fine except for the following note:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>NOTE 1 [file nm.mdp]:<o:p></o:p></p>
<p class=MsoNormal> You are using a plain Coulomb cut-off, which might
produce artifacts.<o:p></o:p></p>
<p class=MsoNormal> You might want to consider using PME electrostatics.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Then, I tried to run<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>mdrun_d -v -nice 0 -deffnm nm<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>However, the command seems to be stuck, and the following
serious warning is produced:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Warning: 1-4 interaction between 64 and 71 at distance 3.114
which is larger than the 1-4 table size 2.000 nm<o:p></o:p></p>
<p class=MsoNormal>These are ignored for the rest of the simulation<o:p></o:p></p>
<p class=MsoNormal>This usually means your system is exploding,<o:p></o:p></p>
<p class=MsoNormal>if not, you should increase table-extension in your mdp file<o:p></o:p></p>
<p class=MsoNormal>or with user tables increase the table size<o:p></o:p></p>
<p class=MsoNormal>Maximum force: 3.48709e+06<o:p></o:p></p>
<p class=MsoNormal>Maximum force probably not small enough to ensure that you
are in an <o:p></o:p></p>
<p class=MsoNormal>energy well. Be aware that negative eigenvalues may occur
when the<o:p></o:p></p>
<p class=MsoNormal>resulting matrix is diagonalized.<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The em2.gro is definitely not having such a large distance
between atoms 64 and 71:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal> 6GLN
CB 64 1.879 1.895 2.423<o:p></o:p></p>
<p class=MsoNormal>.<o:p></o:p></p>
<p class=MsoNormal>.<o:p></o:p></p>
<p class=MsoNormal> 6GLN OE1
71 2.047 1.715 2.642<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>which are only about 0.346 nm apart. Also, the energy
minimization, which used the same options, mutatis mutandis,<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>---------------------- em2.mdp
----------------------------------<o:p></o:p></p>
<p class=MsoNormal>integrator = cg<o:p></o:p></p>
<p class=MsoNormal>nsteps
= 1000<o:p></o:p></p>
<p class=MsoNormal>nstcgsteep = 40 <o:p></o:p></p>
<p class=MsoNormal>implicit_solvent = GBSA<o:p></o:p></p>
<p class=MsoNormal>gb_algorithm = Still ; the default<o:p></o:p></p>
<p class=MsoNormal>rgbradii =
1.0 ; must be equal to rlist <o:p></o:p></p>
<p class=MsoNormal>nstlist
= 10<o:p></o:p></p>
<p class=MsoNormal>rlist
= 1.0<o:p></o:p></p>
<p class=MsoNormal>coulombtype = cut-off<o:p></o:p></p>
<p class=MsoNormal>rcoulomb =
1.0<o:p></o:p></p>
<p class=MsoNormal>vdwtype
= cut-off<o:p></o:p></p>
<p class=MsoNormal>rvdw
= 1.0<o:p></o:p></p>
<p class=MsoNormal>nstenergy = 10<o:p></o:p></p>
<p class=MsoNormal>---------------------------------------------------------------------<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>converged to machine precision having a much smaller maximum
force:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Stepsize too small, or no change in energy.<o:p></o:p></p>
<p class=MsoNormal>Converged to machine precision,<o:p></o:p></p>
<p class=MsoNormal>but not to the requested precision Fmax < 10<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Polak-Ribiere Conjugate Gradients converged to machine
precision in 0 steps,<o:p></o:p></p>
<p class=MsoNormal>but did not reach the requested Fmax < 10.<o:p></o:p></p>
<p class=MsoNormal>Potential Energy = -2.71868395521758e+04<o:p></o:p></p>
<p class=MsoNormal>Maximum force =
4.63436548427712e+02 on atom 584<o:p></o:p></p>
<p class=MsoNormal>Norm of force =
1.25296847735530e+02<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The other problem arises when I try to follow the grompp_d note
above and change in nm.mdp to<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>coulombtype = pme<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>I then have a fatal error:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>ERROR 1 [file nm.mdp]:<o:p></o:p></p>
<p class=MsoNormal> With GBSA, coulombtype must be equal to Cut-off<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Am I doing something wrong or are there still some problems
in the new GBSA option of version 4.5.1?<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Thanks,<o:p></o:p></p>
<p class=MsoNormal>Ehud Schreiber.<o:p></o:p></p>
<p class=MsoNormal> <o:p></o:p></p>
</div>
</body>
</html>