Re: [AMBER] problem Finding entropy

From: Ramin Salimi <ramin.salimi01.utrgv.edu>
Date: Thu, 20 Jul 2017 15:48:01 +0000

Dr. Case

Awesome, sounds like a plan.

will do it.

As the final question, you said 4ns with 2000 snapshots shouldn't give convergence for entropy. why is that so? If the output file shows the energy converges by the end of the run, wouldn't that be enough? I mean, is convergence criterion any different for the entropy?


Regards

Ramin

________________________________
From: David A Case <david.case.rutgers.edu>
Sent: Thursday, July 20, 2017 10:16:47 AM
To: AMBER Mailing List
Subject: Re: [AMBER] problem Finding entropy

On Thu, Jul 20, 2017, Ramin Salimi wrote:
>
> Well, I was planning to calculate the entropy based on the method by Schlitter:

It's been a long time since this came up for me, but the entropy
calculation already done by Amber is more accurate than the Schlitter
approximation. See:

%A I. Andricioaei
%A M. Karplus
%T On the calculation of entropy from covariance matrices of the atomic
fluctuations
%J J. Chem. Phys.
%V 115
%P 6289-6292
%D 2001

It could be a good learning exercise to implement the Schlitter idea and make
comparisons for yourself, but you should first study the existing
quasi-harmonic code so that you understand what it is doing.

Note also that 4ns and only 2000 snapshots are almost certainly well below
what is needed to converge quasiharmonic entropies--but try some experiments
yourself. On the other hand, if the trajectories are not doing what you think
they should, that needs to be investigated. (It seems unlikely that the
strength of the restraint is very important, though.)

...dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=02%7C01%7Cramin.salimi01%40utrgv.edu%7C6cb6e832bdbd40b6dfd508d4cf82747c%7C990436a687df491c91249afa91f88827%7C0%7C0%7C636361606605228121&sdata=b6y3V%2FkhkNgwnJBUnqbkwwyT4X5AiUUj7TYmheynRIg%3D&reserved=0
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jul 20 2017 - 09:00:03 PDT
Custom Search