I am getting the error “pair distance < table inner cutoff” during the first timestep of a tribology simulation using a potential table, obviously - with an inner cutoff. The simulation has under 80 atoms total, with periodic boundary conditions. The strange part is, using a debugger I have found the the atoms who are too close to one another are labeled i=48 and j=1695. The only reasoning I can find for an atom to be index 1695 is if it is a periodic image, but then they certainly shouldn’t be this close to one another. Overall, we are pretty baffled. Any tips? Thanks so much!
input file
3-d 2 metal surfaces coming in contact and then sliding
dimension 3
boundary p p p
units metal
atom_style atomic
Can you print out the distance the code thinks the
2 atoms are apart, and their coordinates? I.e.
add a print statement before the error message
is invoked. You are using a triclinic box,
so periodic images of atoms are also shifted,
so are you sure your input config accounts for that?
Here are the coordinates and such:
i 48 x -14.830252 y 6.634883 z 6.931973
j 1695 x -15.092447 y 7.468805 z 6.935740
delx 0.262194 dely -0.833922 delz -0.003767 rsq 0.764186 tb->innersq 1.234921
The strangest part is that j=1695, when there are fewer than 100 atoms in the simulation. Perhaps it is possible that the table inner cutoff is simply too high. I would really like to make sure that this j value makes sense before proceeding.
Actually, the part I said about the inner cutoff can’t be correct, as with an even shorter time step I’ve gotten values for rsq < 0.1, which is completely unrealistic given the potential (equilibrium ~ 3 )
This output doesn't make sense to me. You said this
is for the 1st timestep, so I assume the coords will
be the same or nearly what is in the input data file.
In the data file the z extent is 0 to 48.96. None of
the 79 atoms has a z coord anywhere close to 6.9.
In fact atom 49 (i = 48) is at 6.2, 4.5, 2.9 in the data
file, but 20 Angs away in the print out below.
The number of ghosts is determined by the outer
cutoff of the table potential. Is it really long,
thus creating a large number of perioidic images?
Thanks Steve. You are completely right. We’ve been able to isolate the problem to something wrong with the potential file (ironically, the only thing I didn’t supply on this listserv - or perhaps, expectedly).