Berkovich indenter simulation


I was trying to simulate nanoindentation process on ZnO surface, using a Berkovich indenter. However lammps doesn’t include options for indenters other than spherical and cylindrical. I tried exploring the plane option and also creating a smaller rigid indenter (of the Berkovich shape) and then obtaining the values as per the interaction between rigid indenter and ZnO surface. However, I would like to know if there is a better way to do it, so that I am able to simulate the interaction forces and activation volume using the berkovich indenter shape. (Depth of indentation 25nm and loading rate 0.1mN/s for the simulation).

I don't know what a Berkovich indenter is. You could

a) write your own region class to use whatever shape you want
b) write your own fix to do indent however you wish