Simulation of the structure of oxygen hole centers in Mg2SiO4 and Mg2SiO4: Cr forsterite crystals by the interatomic potential method

The structure of oxygen hole centers in forsterite crystals has been simulated using the interatomic potential method. The energies of isolated oxygen hole centers, as well as the energies of their clusters with intrinsic and extrinsic defects of the crystal, have been estimated for different arrangements of point defects in the structure. It has been shown that the most energetically favorable position for isolated oxygen hole centers is the O3 position, in which the gain in the formation energy is equal to 0.17 eV as compared to the O2 position and 1.66 eV as compared to the O1 position. The maximum energy gain due to the association energy can be achieved when the oxygen hole centers are located at the vertices of the tetrahedron with a silicon vacancy. The presence of chromium in the forsterite crystal can increase the probability of the formation of silicon vacancies. The obtained results have been discussed in terms of the experimental investigations of the color centers generated in the Mg2SiO4 and Mg2SiO4: Cr crystals under ionizing radiation.