Dear All,
I am working on a system with my own designed force field parameters for a linear pt-c7-n7 system, I determined that force field based on gaff force field and modify some of them to fit my own project. I redefined atom pt, c7, n7 in frcmod file as follows:
MASS
pt 195.09 1.00 atomic mass of pt
c7 12.01 0.360 gaff c1 (sp C)
n7 14.01 0.530 gaff n2 (aliphatic sp2 with two connected atoms)
BOND
pt-c7 366.0 2.042 based on x-ray
n7-c7 366.0 1.13 based on x-ray
ANGLE
pt-c7-n7 42.000 180.0 based on x-ray
DIHEDRAL
IMPROPER
NONBON
pt 4.88 0.4 literature
c7 1.9080 0.0860 same as gaff c1 (cp C)
n7 1.8240 0.1700 same as gaff n1 (OPLS)
=====================================================================
However, i found the following error message everything I run the min and md.
For the min.in, I used.
&cntrl
imin=1, maxcyc=500000, ncyc=40000, ntr=0, nmropt=0,
cut=16, ntb=0, igb=1,
The error message before final results were listed in the out file is as follows:
....Restarted due to Linmin failure...
***Repeated Linmin Failure***
For md.in, I used
&cntrl
imin=0, saltconc=0.3
nstlim=20000, dt=0.002, ntc=1, ntr=0, nmropt=0
ntpr=100, ntwx=100,
cut=16, ntb=0, igb=1,
ntt=3, gamma_ln=1.0,
tempi=300.0, temp0=300.0,
However, the temperature rise significantly during the md simulation
NSTEP=0, TEMP(K)=308.74
NSTEP=100, TEMP(K)=1317.96
NSTEP=200, TEMP(K)=2242.74
.....etc
finally, error message likes
vlimit exceeded for step .....
structure distorted,
TEMP(K)=******
Etot = *****
EPtot=*****
1-4NB=*****
occurs.
Would you please kindly help if I have done something obviously wrong here that lead to the error message above? Please kindly help. Many thanks
Best regards,
Catein
_________________________________________________________________
Explore the seven wonders of the world
http://search.msn.com/results.aspx?q=7+wonders+world&mkt=en-US&form=QBRE
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Fri Dec 05 2008 - 10:48:21 PST
