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

Abstract

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.