<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>You pull rate is ridiculously high (1 nm/ps).<br><br>Berk<br> <br><br><hr id="stopSpelling">Date: Tue, 17 Mar 2009 19:54:50 +0530<br>From: manik.mayur@gmail.com<br>To: gmx-users@gromacs.org<br>Subject: [gmx-users] mdrun hangs with the pull option<br><br>Hi,<br>I am trying to simulate a case where I am pulling a set of atoms with constant velocity with respect to another set of atoms, like:<br><br>XXXXXXXXXXXXXXXXXXXXXXXXXXX --> U (velocity, moving wall, pull group)<br>
ooooooooooooooooooooooooooo<br>ooooooooooooooooooooooooooo -> fluid<br>ooooooooooooooooooooooooooo<br>XXXXXXXXXXXXXXXXXXXXXXXXXXX --> fixed wall (reference group)<br><br>but while running mdrun for the case, it simply becomes unresponsive after 300 steps (no output message even with the -v option). Can somebody suggest what I might be missing. I am appending the md.log and .mdp file for reference however, for quick reference the excerpts are-><br>
<br>*) eql.mdp:<br><br>; COM PULLING <br>; Pull type: no, umbrella, constraint or constant_force<br>pull = constraint<br>pull_geometry = direction<br>pull_start = yes<br>pull_init = <br>
pull_nstxout = 10<br>pull_ngroups = 1<br>pull_group0 = WD<br>pull_group1 = WU<br>pull_vec1 = 1 0 0<br>pull_rate1 = 1<br><br>*) md.log:<br><br>pull = constraint<br>
pull_geometry = direction<br> pull_dim (3):<br> pull_dim[0]=1<br> pull_dim[1]=1<br> pull_dim[2]=1<br> pull_r1 = 1<br> pull_r0 = 1.5<br> pull_constr_tol = 1e-06<br>
pull_nstxout = 10<br> pull_nstfout = 1<br> pull_ngrp = 1<br> pull_group 0:<br> atom (672):<br> atom[0,...,671] = {672,...,1343}<br> weight: not available<br> pbcatom = 1007<br>
vec (3):<br> vec[0]= 0.00000e+00<br> vec[1]= 0.00000e+00<br> vec[2]= 0.00000e+00<br> init (3):<br> init[0]= 0.00000e+00<br> init[1]= 0.00000e+00<br> init[2]= 0.00000e+00<br>
rate = 0<br> k = 0<br> kB = 0<br> pull_group 1:<br> atom (672):<br> atom[0,...,671] = {0,...,671}<br> weight: not available<br> pbcatom = 335<br>
vec (3):<br> vec[0]= 1.00000e+00<br> vec[1]= 0.00000e+00<br> vec[2]= 0.00000e+00<br> init (3):<br> init[0]= 2.08088e+00<br> init[1]= 0.00000e+00<br> init[2]= 0.00000e+00<br>
rate = 1<br> k = 0<br> kB = 0<br><br>The other useful data might be:<br>OS- Ubuntu 8.04<br>Gromacs- 4.0.2<br clear="all"><br>Thanks,<br>Manik Mayur<br>-- <br>
Graduate student<br>Microfluidics Lab<br>Dept. of Mechanical Engg.<br>IIT Kharagpur<br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>