I am trying to calculate a melting temperature based on two phases.
The first thing I want to do, is melt one part of the crystal.
Right now, I have atoms flying out of the box so I was thinking of forcing them to stay in the box by having a flat indenter.
According to the manual, the fix should have this form:
fix 6 all indent 10 plane y 0.0 18. 0.0 200 vel 0. 0. 0. units box
However, I received the error message:
ERROR: Illegal fix indent command
I am using the 9 January version so flat indenter should already be implemented in it.
Do you have any comments?