TY - JOUR

T1 - Ab initio simulation of the IR spectra of pyrope, grossular, and andradite

AU - Zicovich-Wilson, C. M.

AU - Torres, F. J.

AU - Pascale, F.

AU - Valenzano, L.

AU - Orlando, R.

AU - Dovesi, R.

PY - 2008/10

Y1 - 2008/10

N2 - IR spectra of pyrope Mg3Al2Si3O 12, grossular Ca3Al2Si3O 12 and andradite Ca3Fe2Si3O 12 garnets were simulated with the periodic ab initio CRYSTAL code by adopting an all-electron Gaussian-type basis set and the B3LYP Hamiltonian. Two sets of 17 F1u Transverse Optical (TO) and Longitudinal Optical (LO) frequencies were generated, together with their intensities. Because the generation of LO modes requires knowledge of the high frequency dielectric constant ε∞ and Born effective charges, they were preliminary evaluated by using a finite field saw-tooth model and well localized Wannier functions, respectively. As a by-product, the static dielectric constant ε0 was also obtained. The agreement of the present calculated wavenumbers (i.e. peak positions) with the available experimental data is excellent, in that the mean absolute difference for the full set of data smaller than 8 cm-1. Missing peaks in experimental spectra were found to correspond to modes with low calculated intensities. Correspondence between TO and LO modes was established on the basis of the overlap between the eigenvectors of the two sets and similarity of isotopic shifts; as result, the so called LO-TO splitting could be determined. Animation of the normal modes was employed to support the proposed pairing.

AB - IR spectra of pyrope Mg3Al2Si3O 12, grossular Ca3Al2Si3O 12 and andradite Ca3Fe2Si3O 12 garnets were simulated with the periodic ab initio CRYSTAL code by adopting an all-electron Gaussian-type basis set and the B3LYP Hamiltonian. Two sets of 17 F1u Transverse Optical (TO) and Longitudinal Optical (LO) frequencies were generated, together with their intensities. Because the generation of LO modes requires knowledge of the high frequency dielectric constant ε∞ and Born effective charges, they were preliminary evaluated by using a finite field saw-tooth model and well localized Wannier functions, respectively. As a by-product, the static dielectric constant ε0 was also obtained. The agreement of the present calculated wavenumbers (i.e. peak positions) with the available experimental data is excellent, in that the mean absolute difference for the full set of data smaller than 8 cm-1. Missing peaks in experimental spectra were found to correspond to modes with low calculated intensities. Correspondence between TO and LO modes was established on the basis of the overlap between the eigenvectors of the two sets and similarity of isotopic shifts; as result, the so called LO-TO splitting could be determined. Animation of the normal modes was employed to support the proposed pairing.

KW - Ab initio

KW - Dynamical matrix

KW - Gaussian basis sets

KW - Periodic calculations

KW - Vibrational spectrum

UR - http://www.scopus.com/inward/record.url?scp=50149088807&partnerID=8YFLogxK

U2 - 10.1002/jcc.20993

DO - 10.1002/jcc.20993

M3 - Artículo

AN - SCOPUS:50149088807

SN - 0192-8651

VL - 29

SP - 2268

EP - 2278

JO - Journal of Computational Chemistry

JF - Journal of Computational Chemistry

IS - 13

ER -