Dear Mark,<br><br>Thanks for your help. Actually I did separate runs for flexible and rigid models, and didn't use the equilibrated flexible structure to start with rigid model calculations. But the starting configurations are the same. So do you mean I need to reduce the T-coupling time? If the initial water geometry is not the same as the rigid model, does the program will "enforce" that geometry for water molecules, or to "relax" water molecure towards it? <br><br>Many thanks,<br><br>Cherry<br><br><br><br><br><br><br><br><b><i>Mark Abraham <Mark.Abraham@anu.edu.au></i></b> wrote:<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> Cherry Y. Yates wrote:<br>> Dear Mark,<br>> <br>> Thanks for your help, I followed you instruction, and I can get stable <br>> temperature after I remove the temperature coupling. In this case I used <br>> the flexible spc itp file. However if I use a rigid
spc.itp, after <br>> remveing the temperature coupling, the system temperature is shooting up <br>> from 300K to thousands degree, the total energy is going up also instead <br>> of fluctuating. I attached both the itp files and grompp.mdp files. Let <br>> me know what could be wrong. It is really wired to me.<br><br>The change in the model physics from flexible to rigid is severe one - <br>almost all of the water molecules won't fit the rigid geometry, and will <br>thus be in a high-potential-energy state. If the constraint-enforcing <br>algorithm conserves energy, then that potential energy has to be <br>converted to kinetic energy. Then if the simulation lacks temperature <br>coupling, it is at constant energy and the temperature will be <br>necessarily high. Actually you say you observe rising energy - this is <br>probably because your timestep is too large to conserve energy at this <br>very high temperature (and thus atomic velocity).<br><br>So the
take-home lesson is not to do two perturbations of your model <br>physics at one time. If you really want to start with a flexible water <br>model with temperature regulation and end up with rigid water model at <br>without temperature regulation, do the water transition with temperature <br>coupling still on to give the energy somewhere to flow - to the external <br>bath. I expect you'd need a short temperature coupling constant to make <br>this work in practice. Once that's stable, *now* try taking off <br>temperature coupling.<br><br>This suggests the question, why you'd want to do this flexible-rigid <br>transition in an equilibration? It would be simpler to start in the <br>rigid domain, and to do energy minimization to correct any geometries <br>before starting MD simulations.<br><br>Mark<br>_______________________________________________<br>gmx-users mailing list gmx-users@gromacs.org<br>http://www.gromacs.org/mailman/listinfo/gmx-users<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></blockquote><br><p> __________________________________________________<br>Do You Yahoo!?<br>Tired of spam? Yahoo! Mail has the best spam protection around <br>http://mail.yahoo.com