<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div dir="ltr">Dear gmx users<br>I am simulation a system including carbon nanotube+water.I have done these things:<br>1. I copied oplsaa.ff folder in my working directory<br>2. I added following lines to atomname2type.n2t<br>C&nbsp;&nbsp;&nbsp; opls_995&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12.011&nbsp; 2&nbsp;&nbsp;&nbsp; C&nbsp; 0.142&nbsp; C 0.142<br>C&nbsp;&nbsp;&nbsp; opls_996&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12.011&nbsp; 3&nbsp;&nbsp;&nbsp; C&nbsp; 0.142&nbsp; C 0.142&nbsp; C 0.142<br>C&nbsp;&nbsp;&nbsp; opls_997&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12.011&nbsp; 1&nbsp;&nbsp;&nbsp; C&nbsp; 0.142&nbsp; <br>3.I added these to ffbonded.itp<br>[ bondtypes ]<br>&nbsp;C&nbsp;&nbsp; C&nbsp; 1&nbsp;&nbsp; 0.1418&nbsp;&nbsp; 478900<br><br>&nbsp;[ angletypes ]<br>&nbsp;C&nbsp;&nbsp;
 C&nbsp;&nbsp; C&nbsp;&nbsp; 1&nbsp; 120.000&nbsp; 562.2<br><br>[ dihedraltypes ]<br>&nbsp;C&nbsp;&nbsp; C&nbsp;&nbsp; 1&nbsp;&nbsp; 0.000 167.360&nbsp; 1<br>4.I added these to ffnonbonded.itp<br>&nbsp;opls_995&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp; 12.01100&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A&nbsp;&nbsp;&nbsp; 3.85100e-01&nbsp; 4.39600e-01<br>&nbsp;opls_996&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp; 12.01100&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A&nbsp;&nbsp;&nbsp; 3.85100e-01&nbsp; 4.39600e-01<br>&nbsp;opls_997&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp; 12.01100&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A&nbsp;&nbsp;&nbsp; 3.85100e-01&nbsp; 4.39600e-01<br>&nbsp;5.then I run g_x2top with this command:<br>g_x2top -f cnt.pdb -o cnt.top -name CNT -nexcl 3 -kb 478900 -kt 562.2<br>6.then I used editconf and
 genbox:<br>editconf -f cnt.pdb -o cnt.gro -d 1 -bt cubic<br>genbox -cp cnt.gro -cs spc216.gro -o sol.gro -p cnt.top<br>7.I run grompp for EM with this command:<br>grompp -f em.mdp -c sol.gro -p cnt.top -o em.tpr -maxwarn 1<br>mdrun -deffnm em<br>8.again I run grompp for equilibration phase:<br>grompp -f nvt.mdp -c em.gro -p cnt.top -o nvt.tpr -n index.ndx -maxwarn 1<br>mdrun -deffnm nvt<br>However as I looked into the nvt.gro I realized that the cylindrical shape of carbon nanotube<br>&nbsp;has been changed.I am not sure what I have done is correct or not?and how to keep nanotube cylindrical ? any help would be really appreciated.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Best regards<br></div></div></body></html>