<p>&nbsp;Dear all,</p>
<p>I have sent this e-mail to the gmx-users mailing list for a few days, but without any reply until now. I'm not sure if it has something to do with the format, because the link can&rsquo;t be opened when I search the mailing list. So this time I use the plain text and send it once more, would anyone please help me with this?</p>
<p>I&rsquo;m working on some simulations about the adsorption of protein on solid surfaces, which have slab geometry in the x-y plane. In order to reduce the amount of water molecules and at the same time to decrease the unphysical Coulomb interaction between periodic images in the z-direction, an empty layer should be added in the z-direction. To prevent the water molecules from evaporating into the vacuum layer, I&rsquo;m going to use the wall option. However, there are some details that I&rsquo;m not pretty sure about.</p>
<p>What are your suggestions about the choice of wall_atomtype and wall_type? Is it necessary to leave some room for the two walls (or it may lead to high interaction energies between the walls and the system) and how to determine its height (in my simulations, I leave about 1.5 Angstrom, respectively)? Does it have something to do with the wall_atomtype and wall_type?</p>
<p>When choosing nwall=2, pressure coupling and Ewald summation can be used (it is usually best to use semi-isotropic pressure coupling with x/y compressibility set to 0). It means the wall can move in the z-direction? Can the pressure coupling be used in system with fixed atoms and how can we control/calculate the pressure of the mobile phase (as it has been discussed in the paper [Biointerphases 5, 85 (2010)] using the CHARMM package)?</p>
<p>When combining walls with the PME method, it is suggested the eward_geometry be set to 3dc and the wall_eward_zfac be 3. Does this mean there will be an empty layer whose height is 3 times the slab height added to increase the z-dimension of the box? And I&rsquo;m not sure about the exact meaning of the &ldquo;slab height&rdquo;; it seems to be the length/width of the slab (as described in the paper [J. Chem. Phys. 111, 3155 (1999)]).</p>
<p>I&rsquo;m sorry for troubling you with so many questions. But I really need your help. Anyone could help me? Thank you!!</p>
<p>Below is an attachment of mdp file used in my simulation work. If there is anything wrong, please be kind to point it out. Thanks again!</p>
<p>md.mdp<br>
------------------------------------------------------------------------------<br>
title&nbsp;&nbsp;= cyt-c on Au MD <br>
; Run parameters<br>
integrator&nbsp;= md&nbsp;&nbsp;; leap-frog integrator<br>
nsteps&nbsp;&nbsp;= 10000000&nbsp;; 2 * 10000000 = 20000 ps, 20 ns<br>
dt&nbsp;&nbsp;= 0.002&nbsp;&nbsp;; 2 fs<br>
comm-mode&nbsp;&nbsp;&nbsp;&nbsp; = Linear&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; mode for center of mass motion removal</p>
<p>; Output control<br>
nstxout&nbsp;&nbsp;= 10000&nbsp;&nbsp;; save coordinates every 20 ps<br>
nstvout&nbsp;&nbsp;= 10000&nbsp;&nbsp;; save velocities every 20 ps<br>
nstenergy&nbsp;= 10000&nbsp;&nbsp;; save energies every 20 ps<br>
nstlog&nbsp;&nbsp;= 10000&nbsp;&nbsp;; update log file every 20 ps</p>
<p>; Selection of energy groups<br>
energygrps&nbsp;&nbsp;&nbsp; = Protein_HEM&nbsp; GLD ;Group(s) to write to energy file</p>
<p>; Bond parameters<br>
continuation&nbsp;= yes&nbsp;&nbsp;; Restarting after NPT <br>
constraint_algorithm = lincs&nbsp;; holonomic constraints, GolP has been tested with lincs only <br>
constraints&nbsp;= hbonds&nbsp;; bonds with H-atoms constrained<br>
lincs_iter&nbsp;= 1&nbsp;&nbsp;; accuracy of LINCS<br>
lincs_order&nbsp;= 4&nbsp;&nbsp;; also related to accuracy</p>
<p>; Neighborsearching<br>
ns_type&nbsp;&nbsp;= grid&nbsp;&nbsp;; search neighboring grid cells<br>
nstlist&nbsp;&nbsp;= 10&nbsp;&nbsp;; 20 fs<br>
rlist&nbsp;&nbsp;= 1.1&nbsp;&nbsp;; short-range neighborlist cutoff (in nm)<br>
; Periodic boundary conditions<br>
pbc&nbsp;&nbsp;= xy&nbsp;&nbsp;; 2-D PBC</p>
<p>; Method for doing Van der Waals<br>
vdw-type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = switch<br>
rvdw-switch&nbsp;&nbsp; = 0.9<br>
rvdw&nbsp;&nbsp;= 1.0&nbsp;&nbsp;; short-range van der Waals cutoff (in nm)</p>
<p>; Electrostatics<br>
coulombtype&nbsp;= PME&nbsp;&nbsp;; Particle Mesh Ewald for long-range electrostatics<br>
rcoulomb&nbsp;= 1.1&nbsp;&nbsp;; short-range electrostatic cutoff (in nm)<br>
pme_order&nbsp;= 4&nbsp;&nbsp;; cubic interpolation<br>
fourierspacing&nbsp;= 0.12&nbsp;&nbsp;; grid spacing for FFT<br>
ewald_rtol&nbsp;&nbsp;&nbsp;&nbsp; = 1e-5<br>
ewald_geometry&nbsp; = 3dc</p>
<p>; FFT grid size, when a value is 0 fourierspacing will be used<br>
fourier_nx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>
fourier_ny&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>
fourier_nz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>
optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes</p>
<p>; Temperature coupling is on<br>
tcoupl&nbsp;&nbsp;= Nose-Hoover&nbsp;&nbsp;; Nose-Hoover thermostat<br>
tc-grps&nbsp;= Protein_HEM GLD&nbsp;Water_and_ions&nbsp; ; three coupling groups - more accurate<br>
tau_t&nbsp;&nbsp;= 0.5&nbsp;0.5&nbsp;0.5&nbsp;; time constant, in ps<br>
ref_t&nbsp;&nbsp;= 300 &nbsp;300 &nbsp;300&nbsp;; reference temperature, one for each group, in K</p>
<p>; Pressure coupling is on<br>
pcoupl&nbsp;&nbsp;= Parrinello-Rahman&nbsp;; Pressure coupling on in NPT<br>
pcoupltype&nbsp;= semiisotropic&nbsp;; nonuniform scaling of box vectors<br>
tau_p&nbsp;&nbsp;= 1.0&nbsp;1.0&nbsp;&nbsp;; time constant, in ps<br>
ref_p&nbsp;&nbsp;= 1.0&nbsp;1.0&nbsp;&nbsp;; reference pressure, in bar<br>
compressibility = 0&nbsp;4.5e-5&nbsp;&nbsp;; isothermal compressibility of water, bar^-1</p>
<p>; Velocity generation<br>
gen_vel&nbsp;&nbsp;= no&nbsp;&nbsp;; Velocity generation is off</p>
<p>; Non-equilibrium MD stuff<br>
freezegrps&nbsp;&nbsp;&nbsp; = LOCK<br>
freezedim&nbsp;&nbsp;&nbsp;&nbsp; = Y&nbsp; Y&nbsp; Y</p>
<p>; WALLS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
; Number of walls, type, atom types, densities and box-z scale factor for Ewald<br>
nwall&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<br>
wall_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 9-3<br>
wall_r_linpot&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1&nbsp; -1<br>
wall_atomtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = OWT3&nbsp; OWT3&nbsp;&nbsp;&nbsp;&nbsp; ; oxygen atom of TIP3P water in charmm27.ff<br>
wall_density&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 33.4&nbsp; 33.4<br>
wall_ewald_zfac&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 3<br>
---------------------------------------------------------------------------</p>
<p>Yours, sincerely!<br>
Chunwang Peng</p>
<p>Chemistry &amp; Chemical Engineering£¬<br>
South China University of Technology£¬<br>
Tianhe District, Guangzhou, China.<br>
&nbsp;</p>