Hello,
Thank you for all your helps. I solved my solvent box issue. It was actual bubbles caused by a sudden heat up, not wrapping problem. After I changed my heating procedure to a smoother one, I am able to eliminate those bubbles in NVT now.
Regards,
Sixue
-----Original Message-----
From: Daniel Roe [mailto:daniel.r.roe.gmail.com]
Sent: 2015Äê3ÔÂ9ÈÕ 9:59
To: AMBER Mailing List
Subject: Re: [AMBER] Solvent Box Issue When Running AMBER Dynamics
Hi,
You still have the same issue as before - your system's density is way, way too low in both restart files (~0.85 g/mL). You need to allow your system's density (i.e. box dimensions) to properly relax.
According to your output you ran 100 ps in NTV ensemble - this will
*not* get you an equilibrated system and in fact is what is letting those vacuum bubbles form in your simulation. You can run a short (<
15 ps) initial relaxation of solvent in NTV (holding solute fixed), but afterwards you're going to have to run in NTP to let your density equilibrate. If you're interested, I've previously posted the protocol I tend to use on the list, which you can view here:
http://archive.ambermd.org/201307/0550.html
Note of course this is by no means the only way to do it. However, whatever procedure you use to equilibrate, you need to make certain that your system's density and energy values look reasonable before considering running production dynamics.
Hope this helps,
-Dan
On Wed, Mar 4, 2015 at 3:46 PM, Zhang, Sixue <szhang74.illinois.edu> wrote:
> Hello,
>
> I brought out an solvent box issue to this mailing list a few weeks ago. I found that the solvent box became very weird in shape after dynamics (NVT or NPT). It looks like the water was boiling. I was not sure whether it was a wrapping problem or improper equilibration. I then did more tests according to people's suggestions. One Post-doc in my group who is very experienced in AMBER built and equilibrated my system from scratch. She also ended up with this bad box shape. She is pretty sure she equilibrated it properly. In one of my new test, I put high restraints on all atoms, 200 kcal/mol, which means everything almost can't move. If there is no wrapping problem, the solvent box should stay at the same shape after NVT. However, I found some part of the solvents still gradually evolved to a bad shape while the rest part stayed at their original positions as expected. Does anybody has any idea about this issue or has some suggestion on wrapping the system?
>
> Below is the dropbox links for one set of my test files in case you would like to take a look.
> https://www.dropbox.com/sh/wa7z4sapfr1o6od/AAAX1IdVaFsq_xx8GtMuui4Ga?d
> l=0
>
> Thanks a lot,
>
> Sixue Zhang
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
--
-------------------------
Daniel R. Roe, PhD
Department of Medicinal Chemistry
University of Utah
30 South 2000 East, Room 307
Salt Lake City, UT 84112-5820
http://home.chpc.utah.edu/~cheatham/
(801) 587-9652
(801) 585-6208 (Fax)
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Mar 09 2015 - 14:30:02 PDT
