<div>Dear Gromacs Users,</div>
<div><br>I would like to compare the diffusion coefficient of water in different confined spheres with the bulk solvent, I have taken a solvent box of size 7nm, and I have created a sphere with radius 25 Å as lowest and 35 Å as highest spheres. I have not used the pressure coupling for MD, and temp coupling is applied for freezing and non freezing waters separately. In case of the bigger sphere (i. e, 35 Å) the diffusion coefficient is 2.4 10-5 cm2/s (comparable with bulk), but the system with 25 Å radius gives me almost 0 diffusion coefficient. When I open the log file of 25 Å system it shows the initial temperature as zero.</div>

<div><br>Is there any relation between size of freezing group and temperature coupling? How the sphere with 35 Å radius is working fine but not with the radius 25 Å?</div>
<p>Can anybody help me?</p>
<p>Here is my m.mdp file</p>
<p>+++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>title           = AB2130<br>cpp             = /usr/bin/cpp<br>constraints     = all-bonds<br>integrator      = md<br>dt              = 0.002 ; ps !<br>nsteps          = 1500000 ; total 3.0 ns.<br>
nstcomm         = 1<br>nstxout         = 1000 ; collect data     every 2.0 ps<br>nstvout         = 1000 ; collect velocity every 2.0 ps<br>nstfout         = 0<br>nstlog          = 0<br>nstenergy       = 1000 ; collect energy   every 2.0 ps<br>
nstlist         = 10<br>ns_type         = grid<br>rlist           = 1.0<br>coulombtype     = PME<br>rcoulomb        = 1.0<br>rvdw            = 1.0<br>rvdw_switch     = 0.9<br>fourierspacing  = 0.12<br>fourier_nx      = 0<br>
fourier_ny      = 0<br>fourier_nz      = 0<br>pme_order       = 4<br>ewald_rtol      = 1e-5<br>optimize_fft    = yes<br>; Berendsen temperature coupling is on in two groups<br>Tcoupl          = V-rescale<br>tau_t           = 0.1  0.1<br>
tc-grps         = Tmp1 Tmp2<br>ref_t           = 300  0<br>gen_vel         = yes<br>gen_temp        = 300.0<br>acc_grps        = Tmp2<br>accelerate      = 0.1 0.1 0.1<br>gen_seed        = 173529<br>freezegrps      = Freez<br>
freezedim       = Y Y Y<br>+++++++++++++++++++++++++++++++++++++++++++++++++++++++<br clear="all"><br>-- <br>*********************************************<br>J. Srinivasa Rao<br>Post-doctoral Research Associate<br>Computational Biophysics Group<br>
Department of Physics<br>Drexel University<br>3141 Chestnut St<br>Philadelphia, PA 19104<br>Ph:  Off:     215-895-1989<br>       Mob:  704-706-4191<br>Web:<a href="http://jsrao.web.officelive.com/default.aspx">http://jsrao.web.officelive.com/default.aspx</a><br>
**********************************************<br><br><br></p>