A study of the crystal chemistry, electron density distributions, and hydrogen incorporation in the Al₂SiO₅ polymorphs

The Al₂SiO₅ polymorphs have been examined to provide new insights into their chemical bonding, their crystal chemistry, their equations of state, and the incorporation of water in the form of hydroxyl in their structures. The Al₂SiO₅ polymorphs provide a unique structural assemblage for a crystal chemical examination due to the variation in Al coordination in the structures where Al is in 4-fold, 5-fold, and 6-fold in sillimanite, andalusite, and kyanite, respectively. Consequently, the Al₂SiO₅ polymorphs have been examined with a combination of experimental (high pressure X-ray diffraction and Polarized FTIR spectroscopy) and theoretical (VASP and Crystal 98) methods.

An experimental high pressure X-ray diffraction study on andalusite and sillimanite has constrained their equation of state and the pressure derivatives of their bulk modulus with pressure. Additionally, the effect of pressure on the crystal structures has been examined, where the main structural response is compression of the AlO₆ octahedra. Comparatively, compression of the AlO₆ octahedra in andalusite is more anisotropic, while the major direction of axial compressibility in both structures is dependent on the orientation of the AlO6 octahedra.

In order to better understand the crystal chemistry of the Al-O and Si-O bonds in the polymorphs, ELF isosurfaces were examined. ELF isosurfaces represent a graphical representation of the localized electron probability density. Six distinct types of ELF isosurfaces were observed in the Al₂SiO₅ polymorphs resulting from differences in the geometry, coordination, and coordinated cation atomic number surrounding the oxygens within the crystal structures. The ELF was also shown to be isostructurally related to electron density difference maps.

In a combined experimental and theoretical investigation of the Al₂SiO₅ polymorphs, potential protonation sites within the crystal structures were determined at an atomic level with polarized FTIR spectroscopy and analysis of (3,-3) critical points of the negative Laplacian. The polarized FTIR spectra indicate the orientation of the OH dipole in the three polymorphs and the (3,-3) critical points indicate regions of locally concentrated electron density. Potential protonation sites were determined based on the value of the negative Laplacian, the underbonded nature of the oxygens, and the number of surrounding cations.