<div>Hi gmx-users and developers,</div>
<div>&nbsp;</div>
<div>I met a strange error when using wall function combined with user defined potential. Following is a part of the mdp:</div>
<div>&nbsp;</div>
<div>------------------------------------------------------------------------</div>
<div>vdw-type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = user</div>
<div>energygrps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=&nbsp;Surf&nbsp; &nbsp;SOL</div>
<div>&nbsp;</div>
<div>; 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;&nbsp;&nbsp; = 2<br>wall_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = table<br>wall_r_linpot&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1<br>wall_atomtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;= W1&nbsp; W2<br>
wall_density&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;= 1.5&nbsp; 1.5<br>wall_ewald_zfac&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2</div>
<div>energygrp_table&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=&nbsp;Surf &nbsp;wall0&nbsp; Surf&nbsp; wall1&nbsp; SOL&nbsp; wall0&nbsp; SOL&nbsp; wall1 </div>
<div>table-extension&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0</div>
<div>------------------------------------------------------------------------</div>
<div>1) With the command bellow,</div>
<div>&nbsp;</div>
<div>mdrun -table table_Surf_wall0.xvg table_Surf_wall1.xvg table_SOL_wall0.xvg table_SOL_wall1.xvg -deffnm prod</div>
<div>&nbsp;</div>
<div>A fatal error was generated: Library file table_Surf_wall0_Surf_wall0.xvg not found in current dir nor in default directories.</div>
<div>&nbsp;</div>
<div>Then I copy the table6-12.xvg from the top dir, rename it&nbsp;by&nbsp;table.xvg and add it in the command line,</div>
<div>&nbsp;</div>
<div>mdrun -table table6-12.xvg table_Surf_wall0.xvg table_Surf_wall1.xvg table_SOL_wall0.xvg table_SOL_wall1.xvg -deffnm prod</div>
<div>&nbsp;</div>
<div>Am I right? It seems to be ok now. BTW, only the parameters in the replusive colunm in table_SOL_wall1.xvg are none zero (when z&gt;1, these values are aslo set to zero.)</div>
<div>&nbsp;</div>
<div>2) But when I make a parallel run, a crash was found: </div>
<div>&nbsp;</div>
<div>One of the processes started by mpirun has exited with a nonzero exit<br>code.&nbsp; This typically indicates that the process finished in error.<br>If your process did not finish in error, be sure to include a &quot;return<br>
0&quot; or &quot;exit(0)&quot; in your C code before exiting the application.</div>
<div>PID 5798 failed on node n0 (<a href="http://127.0.0.1">127.0.0.1</a>) due to signal 11.<br>-----------------------------------------------------------------------------<br>mpirun failed with exit status 11</div>
<div>&nbsp;</div>
<div>3) Moreover, if table-extension set to 0.0, a fatal error was caught: <br>&#39;An atom is beyond the wall: coordinates -5371173273600.000000 -20382570512384.000000 2778490208256.000000, distance -1.000000<br>You might want to use the mdp option wall_r_linpot&#39;</div>

<div>&nbsp;</div>
<div>I have tested several values, crash is found only when this value is less than 0.2.</div>
<div>&nbsp;</div>
<div>If you have any suggestions, that would be much appreciated. Thanks in<br>advance.<br></div>