As part of ADIS2020 we had two hands-on programming sessions with the results collected here.
We tried a few methods to find vacancies in bulk structures as well tried to get some intuitive understanding of steinhardt parameters and what subset is most useful for structure identification.
See e.g. here We looked a bit into elastic behaviour and crack formation under hydrogen influence using a very simple decohesion zone model.