<DIV>HI everyone,<BR> I want to make a test using tables potential, I choose full-atoms benzene molecules for example.firstly ,I transform the 12-6 LJ potential into table potential, but when mdrun step, I got a error message as follows:<BR>******************************************************************<BR>Getting Loaded...<BR>Reading file ben_md.tpr, VERSION 3.3.3 (single precision)<BR>-----------------------------------------------------------------------------<BR>One of the processes started by mpirun has exited with a nonzero exit<BR>code. This typically indicates that the process finished in error.<BR>If your process did not finish in error, be sure to include a "return<BR>0" or "exit(0)" in your C code before exiting the application.<BR>PID 16534 failed on node n0 (127.0.0.1) due to signal 11.<BR>-----------------------------------------------------------------------------<BR>**************************************************************************<BR>here is my itp and table potential file:<BR>[ defaults ] <BR>; nbfunc comb-rule gen-pairs fudgeLJ fudgeQQ<BR> 1 1 no 1.0 1.0<BR> [ atomtypes ]<BR>; type mass charge ptype c6 c12<BR> BEN 12 0 A 1.00 0.0000</DIV>
<DIV>[ moleculetype ]<BR>; Name nrexcl <BR>BEN 6 <BR> <BR>[ atoms ] <BR>; nr type resnr resid atom cgnr charge mass<BR> 1 BEN 1 BEN C1 1 0 12.0110<BR> 2 BEN 1 BEN C2 2 0 12.0110<BR> 3 BEN 1 BEN C3 3 0 12.0110<BR> 4 BEN 1 BEN C4 4 0 12.0110<BR> 5 BEN 1 BEN C5 5 0 12.0110<BR> 6 BEN 1 BEN C6 6 0 12.0110</DIV>
<DIV>[ bonds ]<BR>; ai aj fu c0, c1, ...<BR> 1 2 2 0.139 10800000.0 0.139 10800000.0 ; CAA CAB<BR> 2 3 2 0.139 10800000.0 0.139 10800000.0 ; CAA CAD<BR> 3 4 2 0.139 10800000.0 0.139 10800000.0 ; CAB CAC<BR> 4 5 2 0.139 10800000.0 0.139 10800000.0 ; CAC CAF<BR> 5 6 2 0.139 10800000.0 0.139 10800000.0 ; CAF CAE<BR> 6 1 2 0.139 10800000.0 0.139 10800000.0 ; CAE CAD</DIV>
<DIV>[ angles ]<BR>; ai aj ak fu c0, c1, ...<BR> 1 2 3 2 120.0 505.0 120.0 505.0 ; HAA CAA CAB<BR> 2 3 4 2 120.0 505.0 120.0 505.0 ; HAA CAA CAD<BR> 3 4 5 2 120.0 505.0 120.0 505.0 ; CAB CAA CAD<BR> 4 5 6 2 120.0 505.0 120.0 505.0 ; CAA CAB HAB<BR> 5 6 1 2 120.0 505.0 120.0 505.0 ; CAA CAB CAC<BR> 6 1 2 2 120.0 505.0 120.0 505.0 ; HAB CAB CAC</DIV>
<DIV>[ dihedrals ]<BR>; ai aj ak al fu c0, c1, m, ...<BR> 1 2 3 4 2 0.0 167.4 0.0 167.4 ; imp CAA CAD CAB HAA<BR> 2 3 4 5 2 0.0 167.4 0.0 167.4 ; imp CAB CAA CAC HAB<BR> 3 4 5 6 2 0.0 167.4 0.0 167.4 ; imp CAC CAB CAF HAC<BR> 4 5 6 1 2 0.0 167.4 0.0 167.4 ; imp CAF CAC CAE HAF<BR> 5 6 1 2 2 0.0 167.4 0.0 167.4 ; imp CAE CAF CAD HAE<BR> 6 1 2 3 2 0.0 167.4 0.0 167.4 ; imp CAD CAA CAE HAD<BR> <BR>part of table potential as follows:</DIV>
<DIV> <BR> 0.36 0 0 0.65757 0 0 0<BR> 0.362 0 0 0.53518 0 0 0<BR> 0.364 0 0 0.42386 0 0 0<BR> 0.366 0 0 0.32271 0 0 0<BR> 0.368 0 0 0.23088 0 0 0<BR> 0.37 0 0 0.14758 0 0 0<BR> 0.372 0 0 0.07213 0 0 0<BR> 0.374 0 0 0.00385 0 0 0<BR> 0.376 0 0 -0.05783 0 0 0<BR> 0.378 0 0 -0.11347 0 0 0<BR> 0.38 0 0 -0.16358 0 0 0<BR> 0.382 0 0 -0.2086 0 0 0<BR> 0.384 0 0 -0.24897 0 0 0<BR> 0.386 0 0 -0.28507 0 0 0<BR> 0.388 0 0 -0.31727 0 0 0<BR> 0.39 0 0 -0.34588 0 0 0<BR> 0.392 0 0 -0.37121 0 0 0<BR> 0.394 0 0 -0.39354 0 0 0<BR> 0.396 0 0 -0.41312 0 0 0<BR> 0.398 0 0 -0.43018 0 0 0<BR> 0.4 0 0 -0.44495 0 0 0<BR> <BR>*****************************<BR>then ,I change somethings in itp files, add a new type C3,but set its parameters the same as BEN, then em, md ,this time, it words well,but the density of system are 10% less than the full-atom simulation with LJ12-6 potential. I want to know why? table potential can't give the accurate results? another quetion is that, when using tables, if the system has only one type of atoms, Gromas can't work well?<BR> <BR> <BR> changes as follows:<BR>********************************<BR>[ defaults ]<BR>; nbfunc comb-rule gen-pairs fudgeLJ fudgeQQ<BR> 1 1 no 1.0 1.0</DIV>
<DIV>[ atomtypes ]<BR>; type mass charge ptype c6 c12<BR> BEN 12 0 A 1.00 0.0000<BR> C3 12 0 A 1.00 0.0000<BR> [ nonbond_params ]<BR>; i j func c6 c12<BR> BEN C3 1 1.00 0.0000</DIV>
<DIV>[ moleculetype ]<BR>; Name nrexcl<BR>BEN 6<BR>[ atoms ]<BR>; nr type resnr resid atom cgnr charge mass<BR> 1 BEN 1 BEN C1 1 0 12.0110<BR> 2 BEN 1 BEN C2 2 0 12.0110<BR> 3 BEN 1 BEN C3 3 0 12.0110<BR> 4 C3 1 C3 C4 4 0 12.0110<BR> 5 C3 1 C3 C5 5 0 12.0110<BR> 6 C3 1 C3 C6 6 0 12.0110<BR> *************************<BR>table_BEN_BEN.xvg,table_BEN_C3.xvg,table_BEN_C3.xvg are all the same.<BR> <BR>Best!<BR></DIV><br><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net/?from=footer">没有广告的终身免费邮箱,www.yeah.net</a>
</span>