Polish Institute of Nuclear Physics Modeling SiC Crystal Defects

Kristallmodell von Siliziumkarbid mit Stufenversetzungen an rot markierten Stellen. Eine einzige kristallografische Ebene ist unten dargestellt. Die Stellen, an denen elektrische Ladungen in benachbarte Schichten gelangen können, sind gelb markiert.
Silicon carbide crystal model with edge dislocations introduced in places marked in red. A single crystallographic plane is presented at the bottom. The places where electric charges can ‘leak’ to neighboring layers are marked in yellow.

Imperfections of crystal structure deeply modify basic properties of the entire material and, in consequence, drastically limit its applications. Using silicon carbide as an example, Polish physicists have shown that even such computationally demanding defects can be examined with atomic accuracy.

For a long time, theoretical physicists have tried to model crystals incorporating realistic defects and predict their influence on the physical properties of a material. The models, based on the results of various experiments, have described the changes in basic properties of a material without explaining the real causes and effects of the occurring phenomena.

A new model of silicon carbide (SiC), built by physicists from the Institute of Nuclear Physics of the Polish Academy of Sciences (IFJ PAN) in Cracow, has allowed them to demonstrate that now it is possible to study crystals ab initio with such complex defects as edge dislocations and to explain their characteristics by processes occurring on an atomic scale. For this, the physicists of the IFJ PAN cooperated with the Institute of Fundamental Research and the Institute of High Pressure Physics of the Polish Academy of Sciences, both located in Warsaw.

»We tried to find the mechanisms responsible at the atomic level for lowering the breakdown voltage in silicon carbide crystals. Our ab initio calculations lead to a qualitative understanding of the problem and contribute to explaining the details of this phenomenon«, says Dr. Jan Łażewski, professor at the IFJ PAN.

Ab initio calculations have been known for some time – Walter Kohn and John Pople were awarded the Nobel Prize for Chemistry in 1998 –, but they were only recently introduced for simulating linear crystal defects. This term is used to describe calculations carried out using quantum mechanics equations, supported only by knowledge about the structure of the atom and the symmetry of crystals. There is no direct information from experiments in such models, which means that they can also be used to analyze materials that have never been studied or even synthesized before. Because of relatively substantial complication of the issue, so far ab initio calculations worked, at most, in the case of point defects, related to vacancies (missing atoms or holes in the crystal structure) as well as admixtures introduced into the crystal.

For their study, the researchers used silicon carbide. At the atomic level, these crystals are composed of many flat layers arranged one on top of each other. Each layer resembles a honeycomb: it consists of hexagonal cells in which the silicon carbide molecules are located vertically in the corners. Each two adjacent layers can be combined in three ways. The multilayer ‘sandwiches’ with different layouts create so-called polytypes, of which there exist more than 250 in the case of silicon carbide. The group from IFJ PAN used the 4H-SiC polymorph.

»When modelling such structures, one of the main problems is computational complexity. A model of pure crystal, devoid of admixtures or dislocations, is characterized by high symmetry and can be calculated even in a few minutes. In order to carry out a calculation for a material with dislocation, we need months working on a high power computer«, emphasizes Dr. Paweł Jochym, professor at the IFJ PAN.

The problems with edge dislocations result from the scale of their influence on the crystal structure of the material. As an illustration, they can be compared to the problem of disguising a gap in a row of tiles on a floor. The gap can be ‘camouflaged’ by moving the tiles of adjacent rows, but the defect will always remain visible. Edge dislocations resulting from the lack of whole lengths or regions of atoms/molecules in individual crystal layers act similarly, affecting the positions of atoms and molecules in many adjacent layers. And since the dislocations can extend over long distances, in practice the disturbances caused by them include the entire crystal.

The most interesting phenomena take place in the dislocation core, i.e. in the vicinity of the edge of the damaged layer of the crystal network. In order to eliminate long-range effects caused by a single dislocation, and thus significantly reduce the number of atoms under consideration, a trick was employed: a second dislocation of the opposite effect was introduced. In this way, the impact of the first dislocation over longer distances was compensated for.

The SiC crystal model consisted of about 400 atoms. The simulations showed that in the layers of crystals, along the edge of the core of the defect, ‘tunnels’ appear in the form of channels with reduced charge density. They lower the potential barrier locally and cause electric charges to ‘leak’ from the valence band. In addition, in the forbidden gap, which in the insulator guarantees a lack of electrical conductivity, conditions appear which reduce its width and effectiveness in limiting the flow of charge. It was shown that these states originate from atoms located in the dislocation core.

»The situation can be compared to a deep, steep ravine that a squirrel is trying to cross. If the bottom of the ravine is empty, the squirrel will not get to the other side. However, if there are a number of trees at the bottom that are high enough, the squirrel can jump over their tops to the other side of the ravine. In the crystal we modelled, the squirrels are the electrical charges, the valence band is one edge of the ravine, the conduction band is the other, and the trees are the aforementioned states associated with the atoms of the dislocation core«, says Prof. Łażewski.

Now that the mechanisms responsible for lowering the threshold of the energy barrier have become known at the atomic level, there is a huge scope for experimentation. The proposed mechanism will have to be verified in order to be able to use it to limit the negative influence of the tested defects. Fortunately, there are already technical possibilities for this.

»The future will verify whether our ideas will be confirmed in their entirety. However, we are confident about the fate of our model and the presented approach to simulating edge dislocations. We already know that the ab initio model has proved its worth in confrontation with certain experimental data«, concludes Prof. Jochym.

Scientific papers

J. Łażewski, et al., Destructive influence of the edge dislocation on 4H-SiC properties – DFT study

J. Łażewski, et al., DFT modelling of the edge dislocation in 4H-SiC, Journal of Materials Science (2019) 54:10737-10745, DOI: 10.1007/s10853-019-03630-5