RE: AMBER: les/nudged elastic band

From: Ross Walker <>
Date: Thu, 31 May 2007 16:59:25 -0700

Hi Tamara,
The thing to remember when running QM/MM simulations is that the actual
bonding information in the prmtop file is not used for the QM region. Hence
while the two end points have to have identical topology files this in no
way prevents you breaking a bond along the reaction coordinate in NEB
assuming that that bond is within the QM region. You will see no difference
in the results if you don't have the bond in both end states or do have the
bond specified in both end states - the choice is yours.
All the best

|\oss Walker

| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- |
| <> | 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: [] On Behalf Of
Tamara Rogers
Sent: Thursday, May 31, 2007 16:20
Subject: AMBER: les/nudged elastic band

Dear Amber users,
I have a bond breaking reaction that I which to follow using the combination
of QMMD and LES/nudged elastic band. I followed the online tutorial and have
a general sense of how to to perform such a simulation. However, in the
tutorial it states that LES requires the all snapshots have the same
topology. Does anyone have any experience using this approach when there is
a change of topology during the reaction (i.e. bond breaking/forming). I
intend to model the region of the changing topology using QM. In particular,
how would I perform a classic transition state search using QMMD in Amber 9?


Luggage? GPS? Comic books?
Check out fitting gifts
mail&p=graduation+gifts&cs=bz> for grads at Yahoo! Search.

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Sun Jun 03 2007 - 06:07:30 PDT
Custom Search