Quantum-Mechanical ab Initio Simulation of the Raman and IR Spectra of Fe3Al2Si3O12 Almandine

The CRYSTAL code, 1976-2020 and beyond, a long story

CRYSTAL is a periodic ab initio code that uses a Gaussian-type basis set to express crystalline orbitals (i.e., Bloch functions). The use of atom-centered basis functions allows treating 3D (crystals), 2D (slabs), 1D (polymers), and 0D (molecules) systems on the same grounds. In turn, all-electron calculations are inherently permitted along with pseudopotential strategies. A variety of density functionals are implemented, including global and range-separated hybrids of various natures and, as an extreme case, Hartree-Fock (HF). The cost for HF or hybrids is only about 3-5 times higher than when using the local density approximation or the generalized gradient approximation. Symmetry is fully exploited at all steps of the calculation. Many tools are available to modify the structure as given in input and simplify the construction of complicated objects, such as slabs, nanotubes, molecules, and clusters. Many tensorial properties can be evaluated by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, first and second hyperpolarizabilities, etc. The calculation of infrared and Raman spectra is available, and the intensities are computed analytically. Automated tools are available for the generation of the relevant configurations of solid solutions and/or disordered systems. Three versions of the code exist: serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated on each core, whereas in the third one, the Fock matrix is distributed for diagonalization. All the relevant vectors are dynamically allocated and deallocated after use, making the code very agile. CRYSTAL can be used efficiently on high performance computing machines up to thousands of cores.

Multiscale Computational Simulation of Amorphous Silicates’ Structural, Dielectric, and Vibrational Spectroscopic Properties

Silicates are among the most abundant and important inorganic materials, not only in the Earth’s crust, but also in the interstellar medium in the form of micro/nanoparticles or embedded in the matrices of comets, meteorites, and other asteroidal bodies. Although the crystalline phases of silicates are indeed present in nature, amorphous forms are also highly abundant. Here, we report a theoretical investigation of the structural, dielectric, and vibrational properties of the amorphous bulk for forsterite (Mg2SiO4) as a silicate test case by a combined approach of classical molecular dynamics (MD) simulations for structure evolution and periodic quantum mechanical Density Functional Theory (DFT) calculations for electronic structure analysis. Using classical MD based on an empirical partial charge rigid ionic model within a melt-quenching scheme at different temperatures performed with the GULP 4.0 code, amorphous bulk structures for Mg2SiO4 were generated using the crystalline phase as the initial guess. This has been done for bulk structures with three different unit cell sizes, adopting a super-cell approach; that is, 1 × 1 × 2, 2 × 1 × 2, and 2 × 2 × 2. The radial distribution functions indicated a good degree of amorphization of the structures. Periodic B3LYP-geometry optimizations performed with the CRYSTAL14 code on the generated amorphous systems were used to analyze their structure; to calculate their high-frequency dielectric constants (ε∞); and to simulate their IR, Raman, and reflectance spectra, which were compared with the experimental and theoretical crystalline Mg2SiO4. The most significant changes of the physicochemical properties of the amorphous systems compared to the crystalline ones are presented and discussed (e.g., larger deviations in the bond distances and angles, broadening of the IR bands, etc.), which are consistent with their disordered nature. It is also shown that by increasing the unit cell size, the bulk structures present a larger degree of amorphization.

B3LYP periodic study of the physicochemical properties of the nonpolar (010) Mg-pure and fe-containing olivine surfaces

B3LYP periodic simulations have been carried out to study some physicochemical properties of the bulk structures and the corresponding nonpolar (010) surfaces of Mg-pure and Fe-containing olivine systems; i.e., Mg2SiO4 (Fo) and Mg1.5Fe0.5SiO4 (Fo75). A detailed structural analysis of the (010) Fo and Fo75 surface models shows the presence of coordinatively unsaturated metal cations (Mg(2+) and Fe(2+), respectively) with shorter metal-O distances compared to the bulk ones. Energetic analysis devoted to the Fe(2+) electronic spin configuration and to the ion position in the surfaces reveals that Fe(2+) in its quintet state and placed at the outermost positions of the slab constitutes the most stable Fe-containing surface, which is related to the higher stability of high spin states when Fe(2+) is coordinatively unsaturated. Comparison of the simulated IR and the corresponding reflectance spectra indicates that Fe(2+) substitution induces an overall bathochromic shift of the spectra due to the larger mass of Fe compared to Mg cation. In contrast, the IR spectra of the surfaces are shifted to upper values and exhibit more bands compared to the corresponding bulk systems due to the shorter metal-O distances given in the coordinatively unsaturated metals and to symmetry reduction which brings nonequivalent motions between the outermost and the internal modes, respectively.

Raman spectrum of pyrope garnet. A quantum mechanical simulation of frequencies, intensities, and isotope shifts

The Raman spectrum of pyrope garnet is simulated in ab initio quantum mechanical calculations, using an all-electron Gaussian-type basis set and the hybrid B3LYP functional. Frequencies calculated for the 25 Raman-active modes are in excellent agreement with the several sets of experimental data, with the mean absolute difference ranging from 4 to 8 cm(-1). Comparison of the computed and experimental spectrum shows excellent agreement for most of the intensities as well. Modes missing from experiment are shown to be characterized by low (computed) intensity. Spurious peaks in the experimental spectra are also identified. The isotopic effect has been simulated for (24)Mg → (26)Mg substitution and shows excellent agreement with shifts reported in one of the experiments. Agreement is excellent for all but one mode, which turns out to be attributed to the wrong symmetry in the experiment.