Ab Initio Study of the Vibrational Spectrum and Related Properties of Crystalline Compounds; the Case of CaCO3 Calcite

Abstract The static and high frequency dielectric tensors, Born effective charges, vibrational spectrum at the Γ point, TO-LO splitting and IR intensities of calcite CaCO3 have been calculated with the periodic ab initio CRYSTAL program, with five different basis sets of increasing size and four different Hamiltonians (HF, LDA, PBE, B3LYP). B3LYP is shown to perform better than the other options, in particular of LDA and PBE that are often used for the calculation of the vibrational spectrum of crystalline solids. When comparing B3LYP and experimental frequencies, the mean absolute difference is as small as 8.5 cm-1; this number reduces to 4.8 cm-1 if the two lowest experimental frequencies, that we suspect to be affected by a relatively large error, are excluded from statistics. Static and high frequency dielectric tensors, as well as IR intensities computed with the same hybrid scheme (B3LYP) compare quite favourably with experiment. The full set of modes is characterized by various tools including isotopic substitution, “freezing” one of the two subunits (Ca2+ or CO32-) and graphical representations. A general tool has been implemented, that permits the automatic generation of the animation of the full set of modes starting from the CRYSTAL output (available at www.crystal.unito.it/vibs/calcite).

[1]  P. Dirac.  Note on Exchange Phenomena in the Thomas Atom , 1930, Mathematical Proceedings of the Cambridge Philosophical Society.

[2]  M. Born,et al.  Dynamical Theory of Crystal Lattices , 1954 .

[3]  C. C. J. Roothaan,et al.  Self-Consistent Field Theory for Open Shells of Electronic Systems , 1960 .

[4]  R. Howie,et al.  An Introduction to the Rock-Forming Minerals , 1966 .

[5]  S. Naidu,et al.  Precision lattice parameters and thermal expansion of calcite , 1967 .

[6]  K. Hellwege,et al.  Zwei-Phononen-Absorptionsspektren und Dispersion der Schwingungszweige in Kristallen der Kalkspatstruktur , 1970 .

[7]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[8]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[9]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[10]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[11]  V. Streltsov,et al.  X-ray study of the electron density in calcite, CaCO3 , 1993 .

[12]  P. McMillan,et al.  Thermodynamic properties and isotopic fractionation of calcite from vibrational spectroscopy of 18O-substituted calcite , 1996 .

[13]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[14]  M. Böttcher,et al.  Characterization of inorganic and biogenic magnesian calcites by Fourier Transform infrared spectroscopy , 1997 .

[15]  C. Geiger,et al.  Low-temperature single-crystal Raman spectrum of pyrope , 2000 .

[16]  R. Dovesi,et al.  A general method to obtain well localized Wannier functions for composite energy bands in linear combination of atomic orbital periodic calculations , 2001 .

[17]  Alfredo Pasquarello,et al.  Raman scattering intensities in α-quartz: A first-principles investigation , 2001 .

[18]  Roberto Dovesi,et al.  Well localized crystalline orbitals obtained from Bloch functions: The case of KNbO 3 , 2001 .

[19]  Bartolomeo Civalleri,et al.  Hartree–Fock geometry optimisation of periodic systems with the Crystal code , 2001 .

[20]  Roberto Dovesi,et al.  Characterization of the electronic structure of crystalline compounds through their localized Wannier functions , 2002 .

[21]  Roberto Dovesi,et al.  Implementation of the finite field perturbation method in the CRYSTAL program for calculating the dielectric constant of periodic systems , 2003, J. Comput. Chem..

[22]  R. Orlando,et al.  The vibrational spectrum of calcite (CaCO3): an ab initio quantum-mechanical calculation , 2004 .

[23]  R. Orlando,et al.  Calculation of the vibration frequencies of α‐quartz: The effect of Hamiltonian and basis set , 2004, J. Comput. Chem..

[24]  Bartolomeo Civalleri,et al.  The calculation of the vibrational frequencies of crystalline compounds and its implementation in the CRYSTAL code , 2004, J. Comput. Chem..

[25]  Roberto Dovesi,et al.  Vibration Frequencies of Mg3Al2Si3O12 Pyrope. An ab initio study with the CRYSTAL code. , 2005, The journal of physical chemistry. B.

[26]  R. Dovesi,et al.  Influence of the exchange-correlation functional in all-electron calculations of the vibrational frequencies of corundum (α-Al2O3) , 2006 .