The problem of numerical accuracy in the calculation of vibrational frequencies of crystalline compounds from the hessian matrix is discussed with reference to alpha-quartz (SiO(2)) as a case study and to the specific implementation in the CRYSTAL code. The Hessian matrix is obtained by numerical differentiation of the analytical gradient of the energy with respect to the atomic positions. The process of calculating vibrational frequencies involves two steps: the determination of the equilibrium geometry, and the calculation of the frequencies themselves. The parameters controlling the truncation of the Coulomb and exchange series in Hartree-Fock, the quality of the grid used for the numerical integration of the Exchange-correlation potential in Density Functional Theory, the SCF convergence criteria, the parameters controlling the convergence of the optimization process as well as those controlling the accuracy of the numerical calculation of the Hessian matrix can influence the obtained vibrational frequencies to some extent. The effect of all these parameters is discussed and documented. It is concluded that with relatively economical computational conditions the uncertainty related to these parameters is smaller than 2-4 cm(-1). In the case of the Local Density Approximation scheme, comparison is possible with recent calculations performed with a Density Functional Perturbation Theory method and a plane-wave basis set.

The central-zone vibrational spectrum of alpha-quartz (SiO2) is calculated by building the Hessian matrix numerically from the analytical gradients of the energy with respect to the atomic coordinates. The nonanalytical part is obtained with a finite field supercell approach for the high-frequency dielectric constant and a Wannier function scheme for the evaluation of Born charges. The results obtained with four different Hamiltonians, namely Hartree-Fock, DFT in its local (LDA) and nonlocal gradient corrected (PBE) approximation, and hybrid B3LYP, are discussed, showing that B3LYP performs far better than LDA and PBE, which in turn provide better results than HF, as the mean absolute difference from experimental frequencies is 6, 18, 21, and 44 cm(-1), respectively, when a split valence basis set containing two sets of polarization functions is used. For the LDA results, comparison is possible with previous calculations based on the Density Functional Perturbation Theory and usage of a plane-wave basis set. The effects associated with the use of basis sets of increasing size are also investigated. It turns out that a split valence plus a single set of d polarization functions provides frequencies that differ from the ones obtained with a double set of d functions and a set of f functions on all atoms by on average less than 5 cm(-1).

No abstract

Vibrations / Calcite / IR Intensities / Ab Initio SimulationThe static and high frequency dielectric tensors, Born effective charges, vibrational spectrum at the Γ point, TO-LO splitting and IR intensities of calcite CaCO 3 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 (Ca 2+ or CO 3 2− ) 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).

The vibrational spectrum of Mg(3)Al(2)Si(3)O(12) pyrope is calculated at the Gamma point by using the periodic ab initio CRYSTAL program that adopts an all-electron Gaussian-type basis set and the B3LYP Hamiltonian. The full set of frequencies (17 IR active, 25 RAMAN active, 55 silent modes) is calculated. The effect of the basis set and of the computational parameters on the calculated frequencies is discussed. It is shown that the mean absolute difference with respect to the experimental IR and RAMAN data is as small as 6 and 8 cm(-1), respectively. The IR and RAMAN modes are fully characterized by various tools such as isotopic substitution, direct inspection of the eigenvectors, and graphical representation. The present calculation permits to clarify some of the assignment and interpretation problems raised by experiment and previous simulations with force fields.