On 11/21/2012 11:41 AM, Ali Alizadeh wrote:
Dear Ivan
Thank you for your reply,
Dear Ali
please keep the discussion open to everybody...
In this way you will be sure that someone else will read and interact
with us, correcting me if I say something wrong to you..;-)
1- In your opinion, Can i simulate that system?
In my (humble) opinion:
1)Of course you can simulate that system...however I doubt that,
without starting from the exact initial configuration with the exactly
same set-up, you can get the same results (i.e. see the nucleation).
The onset of ice nucleation is a random process and requires very long
simulation (the paper that you posted was analysing micro-second
trajectories!!).
There is the risk that you could try several different initial
configurations at several temperature without getting anything.
However, read carefully that paper, I do not remember all the details.
If you are interested in ice crystal growth, I would suggest to start
with an initial water/ice system: at temperature below the melting one
you will see formation of new ice starting from the initial ice matrix.
There are several works in literature on this.
2- How can i use rigid TIP4P model of water?
2) I would rather use other water model that have been explicitly
tested for Ice
(e.g. TIP4P/2005, TIP4P/Ice, TIP5P-Ew, NE6)
Best
Ivan
Sincerely
Ali Alizadeh
On 11/21/2012 09:43 AM, Ivan Gladich wrote:
Dear Ali
the paper that you are citing is using a rigid TIP4P water model
As far as I know, emtol is relevant only for minimization or molecular
dynamics with shell particle or flexible constraints.
Therefore, as Justin told you, the emtol value should be irrelevant.
Concerning this paper, I would like to warn you that h_omogeneous ice
nucleation from bulk water_ with explicit water molecule is very rare
event...
It's depends from the initial condition and it requires very long
simulations.
Indeed, if I remember well this paper, they observed homogeneous ice
crystal formation in only one of their trajectories...read carefully
the paper!
Best
Ivan
On 11/21/2012 07:47 AM, Ali Alizadeh wrote:
On 11/20/12 3:38 PM, Ali Alizadeh wrote:
Dear Justin
Thank you ,
This link is the paper that i want to reproduce its results,
Title: Molecular dynamics simulation of
the ice nucleation and growth
process leading to water freezing
http://www.nature.com/nature/journal/v416/n6879/pdf/416409a.pdf
I am using MD instead of steep, and i have a problem, I do not know
value of emtol so that energy of system reaches to -51 kj/mol, i
determine it then run it but it reaches to -30 kj/mol(only for 1ns)
original results :
http://alichemical.persiangig.com/document/11.jpg
Justin wrote:
> The value of emtol is only relevant in EM and has no
relationship whatsoever to
> MD. In fact, this keyword is ignored if an EM algorithm
is not chosen as the
>integrator.
This means When i use the md-vv integrator, emtol is ignored.
>No one can help you if you don't fully describe what
you're doing. That appears
>to be a figure from some paper - what's going on? If it
is from a paper, what
>were the methods? How was the plot produced? Otherwise
no one will have any
idea what's going on. You can't just force a system to
adopt some series of
>configurations that yield a given potential energy.
>Please describe in full what you are doing, and provide
answers to all the
>questions posed above, or else I'm afraid I'm not going to
have anything else
>useful to say.
-Justin
--
------
Ivan Gladich, Ph.D.
Postdoctoral Fellow
Academy of Sciences of the Czech Republic
Institute of Organic Chemistry and Biochemistry AS CR, v.v.i.
Flemingovo nám. 2.
166 10 Praha 6
Czech Republic
Tel: +420775504164
e-mail: ivan.glad...@uochb.cas.cz
web page:http://www.molecular.cz/~gladich/
-----
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists