Modeling the polychromism of oxide minerals: The case of alexandrite and cordierite

In this work, we investigate the spectroscopic properties of photochromic alexandrite and cordierite by TD-DFT. The objective is to assess the TD-DFT for the simulation of pleochroism (change of color depending on the crystallographic direction of the observation) and the change of color as a function of the light source. For these simulations, we compared an embedding where dangling bonds are saturated by hydrogen atoms and an electrostatic embedding. The electrostatic embedding provided numerically more stable results and allowed a good reproduction of the pleochroism of cordierite, based on a Fe²⁺-Fe³⁺ intervalence charge transfer transition. However, the pleochroism of alexandrite is not as well reproduced, suggesting that TD-DFT has some difficulties to reproduce the anisotropy of the transition dipole moment, an aspect that is not deeply documented in the literature.