[gmx-users] Energy minimization problem in Choline Chloride+Urea mixture in Gromacs

Suman Das sumandas212 at gmail.com
Fri Jan 2 17:21:24 CET 2015


Hi,
 I am trying to simulate a mixture of Choline Chloride+Urea mixture with
33:67 ratio at 298 K. I have taken the force field of choline from this
paper :http://pubs.acs.org/doi/full/10.1021/jp404619x
<https://www.researchgate.net/go.Deref.html?url=http%3A%2F%2Fpubs.acs.org%2Fdoi%2Ffull%2F10.1021%2Fjp404619x>,
chloride from http://pubs.acs.org/doi/full/10.1021/jp063901o
<https://www.researchgate.net/go.Deref.html?url=http%3A%2F%2Fpubs.acs.org%2Fdoi%2Ffull%2F10.1021%2Fjp063901o>
and urea from http://pubs.acs.org/doi/full/10.1021/jp030534x
<https://www.researchgate.net/go.Deref.html?url=http%3A%2F%2Fpubs.acs.org%2Fdoi%2Ffull%2F10.1021%2Fjp030534x>.
During the minimization process Fmax < 200 but potential energy was quite
high~5.4465270e+04. And finally the simulation crushes.  I have used this
force field of urea previously which matches well with experimental data.
So I have tried with another force field of choline (
http://pubs.acs.org/doi/full/10.1021/jp301270y
<https://www.researchgate.net/go.Deref.html?url=http%3A%2F%2Fpubs.acs.org%2Fdoi%2Ffull%2F10.1021%2Fjp301270y>)
but exactly the same thing happens. I have also tried the freeze group
option. But the result is same. I built the box using packmol. I took the
initial box length 5nm.

regards,
Suman


More information about the gromacs.org_gmx-users mailing list