AMBER Archive (2006)Subject: Re: AMBER: how to include EXTRA POINTS in AMBER8
From: Rachel (comeonsos_at_googlemail.com)
Date: Mon Dec 04 2006 - 11:03:11 CST
Dear David,
Thanks for your reply. I have tried to put a finite mass on that extra point
(I used mass = 1.008 as the normal hydrogen atom), and the result is even I
used very very small timestep (e.g. 0.1fs), I don't have any problem with
vlimit, and the error message is:
####################################################################
NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 0.00 PRESS =
0.0
Etot = -217128.3786 EKtot = 0.0000 EPtot = -
217128.3786
BOND = 400.3220 ANGLE = 2165.2216 DIHED =
6798.3874
1-4 NB = 2513.3831 1-4 EEL = 34285.2873 VDWAALS =
25111.1035
EELEC = -288402.0835 EHBOND = 0.0000 RESTRAINT =
0.0000
Ewald error estimate: 0.8763E-04
------------------------------------------------------------------------------
Coordinate resetting (SHAKE) was not accomplished
within 3000 iterations
Note: This is usually a symptom of some deeper
problem with the energetics of the system.
#####################################################################
And about the name, I used 'EPH' and in xleap, I defined its type as 'EP',
and element 'Lp', which is the same as TIP4P, TIP5P water.
Thanks a lot for your kind concern, and I really appreciate your help and
suggestions.
Best regards,
Rachel
On 12/4/06, David A. Case <case_at_scripps.edu> wrote:
>
> On Mon, Dec 04, 2006, Rachel wrote:
> >
> > I have been trying to include some massless extra points in my system,
> > however not successful, and the program kept complaining that 'vlimit
> > exceeded for step 0, vmax=infinity',
>
> First, did you try making the masses non-zero to see if the error changes
> or
> goes away?
>
> Second, do the names of your extra point atoms begin with the letters
> "EP"?
>
> ...dac
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|