How to define complex defect

Hi all,

I am trying to define a defect that is composed of an interstitial and a vacancy (a color center).

How can I achieve this in pymatgen? Do I have to define multiple defects and then combine them somehow?

Thank you for your help!

Hi AxLamelas, sorry for the late reply.

Currently, pymatgen supports the creation of point defects, but not point defect clusters, unfortunately.

The best you can do is, as you say, to create two defects and then combine them. What pymatgen can do is help you identify the unique sites. Perhaps you can you use the InterstitialGenerator located in pymatgen.analysis.defects.generators in order to find the different possible sites, then you can simply remove an atom to make the vacancy?

Part of the reason why we do not, yet, have a generator for these clusters is the number of permutations can quickly become massive. An interstitial site, for example, can be identified for a primitive cell and has N sites. A di-interstitial, however, does not have 2N possibilities, rather it the number of possibilities depends on how large of a supercell you want, because the separation of the two interstitials would be a new variable. Hopefully we will have a generator to take all of this into account some time in the future, but it isn’t in our immediate plans.

Good luck on your work and let me know if you have further questions.

1 Like

Thank you very much for the answer.
I created a custom class for dealing with the structure generations as this was my immediate problem. I’m defining particular defects by hand so the generators are not of immediate concern.
I’ll implement additional functionality when the need arises based on pymatgen’s implementation.
I’ll also have a look at your suggestion.

Sorry for the delay in my response.