<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Hi,<br><br>Shake is not relevant for water and also a time step of 2 fs should be fine.<br>The cut-off's are the problem. You have a buffer size of 0.1 nm,<br>which is already smaller than 2 times the distance from the center of geometry<br>of a water molecule to a hydrogen. You need some additional distance<br>for water diffusion. I would use a buffer of 0.25 to 0.3 nm.<br>You don't list you coulombtype setting. Use pme (or if you want perfect energy<br>conservation: pme-switch), you can also use reaction-field-zero if you really don't<br>want to use PME.<br>Then use nstlist=-1, run a short simulation and check in at the end of your log file<br>that the neighbor list lifetime is somewhere between 5 and 20 steps.<br><br>We should have a wiki entry for such details. Maybe there is one, but I was too lazy<br>to check or make one.<br><br>Berk<br><br>> Date: Mon, 1 Mar 2010 09:16:36 +0100<br>> From: r.friedman@bioc.uzh.ch<br>> To: gmx-users@gromacs.org<br>> Subject: Re: [gmx-users] NVE of water<br>> <br>> Hi Andrea,<br>> Did you use double precision? Also, I'd try a lower dt (say 1fs) and<br>> shake tolerance (maybe 1E-8).<br>> Good luck,<br>> Ran<br>> <br>> Andrea wrote:<br>> > Dear users,<br>> > <br>> > for test purposes in order to set up a bigger system, I try to run NVE<br>> > simulations of SPC water, but the energy increases very rapidely. My<br>> > guess is that the cutoffs I use are not good for water. I that the<br>> > case ( I would be grateful for a good reference for suitable SPC water<br>> > parameters) or do I miss something else?<br>> > My parameter file for the NVE is:<br>> > <br>> > title = NVE<br>> > cpp = /lib/cpp<br>> > integrator = md<br>> > dt = 0.002 ; ps ! = 2 fs<br>> > nsteps = 50000 ; total 100 ps<br>> > nstxout = 5000<br>> > nstvout = 5000<br>> > nstxtcout = 0<br>> > nstlog = 5000<br>> > nstenergy = 5000<br>> > nstlist = 10<br>> > ns_type = grid<br>> > rlist = 1.1<br>> ><br>> > unconstrained-start = yes<br>> > constraints = all-bonds<br>> > constraint_algorithm = shake<br>> > shake_tol = 0.0001<br>> > ;VdW<br>> > vdwtype = Switch<br>> > rvdw = 1.0 ; rvdw+ (0.1:0.3)= rlist<br>> > rvdw_switch = 0.9<br>> > gen_vel = no ; yes<br>> > gen_temp = 300<br>> > gen_seed = -1<br>> > ;Temperature coupling<br>> > tc_grps = system<br>> > tcoupl = no ;nose-hoover<br>> > tau_t = 0.1<br>> > ref_t = 300<br>> > ;Pressure coupling<br>> > pcoupl = no<br>> > optimize_fft = yes<br>> > <br>> > Any suggesions are really welcome.<br>> > <br>> > Thank you.<br>> > <br>> > Regards,<br>> > Andrea Muntean<br>> <br>> <br>> -- <br>> ------------------------------------------------------<br>> Ran Friedman<br>> Postdoctoral Fellow<br>> Computational Structural Biology Group (A. Caflisch)<br>> Department of Biochemistry<br>> University of Zurich<br>> Winterthurerstrasse 190<br>> CH-8057 Zurich, Switzerland<br>> Tel. +41-44-6355559<br>> Email: r.friedman@bioc.uzh.ch<br>> Skype: ran.friedman<br>> ------------------------------------------------------<br>> <br>> -- <br>> gmx-users mailing list gmx-users@gromacs.org<br>> http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> Please search the archive at http://www.gromacs.org/search before posting!<br>> Please don't post (un)subscribe requests to the list. Use the <br>> www interface or send it to gmx-users-request@gromacs.org.<br>> Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>                                            <br /><hr />New Windows 7: Find the right PC for you. <a href='http://windows.microsoft.com/shop' target='_new'>Learn more.</a></body>
</html>