On 6/14/17 5:18 PM, Mohsen Ramezanpour wrote:
Hi Justin,

Thanks for your reply,

Here are the outputs for EM step:

Steepest Descents converged to Fmax < 100 in 5198 steps
Potential Energy  = -2.6418572e+05
Maximum force     =  9.5666359e+01 on atom 54770
Norm of force     =  5.7127428e+00

It seems okay to me.

I chose the setting you mentioned.

Interestingly, it also works for MD with dt=1 fs. Not for dt=2 fs, though.
Based on visualization, some lipid chains are crossed which I think this
cause the problem. Not sure how can I fix them.


Surprising. With such a low force, I wouldn't expect there to be any overlap or detrimental close contact. An initial equilibration with dt = 1 fs is a possibility, or a more robust method for constructing the system, whatever that may be.

-Justin

--
==================================================

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==================================================
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.

Reply via email to