Recently, quantum embedding methods have emerged as efficient approaches for calculating optical transitions involving strongly correlated electrons. These methods involve mapping the infinite bulk system onto an impurity model embedded within a bath. By partitioning the full system into smaller fragments, these techniques address the computational complexity associated with extended systems featuring localized correlation. The fragment with the region of interest is treated with an expensive yet accurate method and is embedded in the rest of the system termed as the environment, which is solved at a low cost, mean field level of theory. Various quantum embedding theories have been reported in the literature leveraging the different central embedding variable. These variables include the density, the density matrix and the Green's function. In this study, we employ the periodic version of density matrix embedding theory (pDMET). Due to the localized nature of strongly correlated defects, such as F centers in MgO, they are well-suited for embedding methods.