AMBER Archive (2006)

Subject: RE: AMBER: the total energy is not converge as the cutoff increase

From: Ross Walker (ross_at_rosswalker.co.uk)
Date: Mon Mar 06 2006 - 09:39:51 CST


Hi Thanyarat,
 
What you are seeing is to be expected. You are not using periodic boundaries
and I assume are doing a gas phase simulation without a dielectric or any
implicit solvent treatment. Then what you are seeing is correct.
Electrostatic interactions in gas phase are extremely long range. So long
range in fact that I would argue that you cannot justify the use of a cutoff
with a gas phase simulation.
 
Why do you want to simulate your system in gas phase? If you do an explicit
solvent calculation with PME you should see the electrostatic energy
converge around 8 angstroms as these types of calculations include the full
infinite electrostatic interactions via a combination of real and reciprocal
space sums.
 
All the best
Ross

/\
\/
|\oss Walker

| Department of Molecular Biology TPC15 |
| The Scripps Research Institute |
| Tel: +1 858 784 8889 | EMail:- ross_at_rosswalker.co.uk |
| http://www.rosswalker.co.uk <http://www.rosswalker.co.uk/> | PGP Key
available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.

 

  _____

From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf Of
thanyarat Udom
Sent: Monday, March 06, 2006 00:52
To: amber_at_scripps.edu
Subject: AMBER: the total energy is not converge as the cutoff increase

Dear Amber users
        
          I want to run an energy minimization of protein/water/counter ions
system. First, I have tested cutoff the value of non-bonded. I summarized
the output from the first step (NSTEP = 1) of energy minimization as
follows.
cut off Total Energy(NSTEP=1) BOND ANGLE DIHED
VDWAALS EEL
  10 -145910 2629.283 5045.8942 10443.8164
30238.3127 -243907.946
 12 -262380 2629.283 5045.8942 10443.8164
29518.1267 -356388.753
 15 -302340 2629.283 5045.8942 10443.8164
29052.2034 -395395.948
 18 -260040 2629.283 5045.8942 10443.8164
28855.7353 -355521.277
 20 -571410 2629.283 5045.8942 10443.8164
28731.883 -667743.075

 Surprisingly, the total energy is not converge as the cutoff increase.
Particularly, the results at cutoff = 20 Å look very strange. The
minimization has not been employed PBC (ntb =0). Please help me.

                                                             Best Regard

                                                  Thanyarat
Udommaneethanakit

  _____

Yahoo! Mail
Bring photos to life! New
<http://pa.yahoo.com/*http://us.rd.yahoo.com/evt=39174/*http://photomail.mai
l.yahoo.com> PhotoMail makes sharing a breeze.

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu