No Arabic abstract
Based on density-functional perturbation theory we have computed the photoelastic tensor of a model of sodium silicate glass of composition (Na$_2$O)$_{0.25}$(SiO$_2$)$_{0.75}$ (NS3). The model (containig 84 atoms) is obtained by quenching from the melt in combined classical and Car-Parrinello molecular dynamics simulations. The calculated photoelastic coefficients are in good agreement with experimental data. In particular, the calculation reproduces quantitatively the decrease of the photoelastic response induced by the insertion of Na, as measured experimentally. The extension to NS3 of a phenomenological model developed in a previous work for pure a-SiO$_2$ indicates that the modulation upon strain of other structural parameters besides the SiOSi angles must be invoked to explain the change in the photoelstic response induced by Na.
We use textit{ab initio} molecular dynamics simulations to investigate the properties of the dry surface of pure silica and sodium silicate glasses. The surface layers are defined based on the atomic distributions along the direction ($z-$direction) perpendicular to the surfaces. We show that these surfaces have a higher concentration of dangling bonds as well as two-membered (2M) rings than the bulk samples. Increasing concentration of Na$_2$O reduces the proportion of structural defects. From the vibrational density of states, one concludes that 2M rings have a unique vibrational signature at a frequency $approx850$~cm$^{-1}$, compatible with experimental findings. We also find that, due to the presence of surfaces, the atomic vibration in the $z-$direction is softer than for the two other directions. The electronic density of states shows clear the differences between the surface and interior and we can attribute these to specific structural units. Finally, the analysis of the electron localization function allows to get insight on the influence of local structure and the presence of Na on the nature of chemical bonding in the glasses.
Based on density-functional perturbation theory we have computed from first principles the photoelastic tensor of few crystalline phases of silica at normal conditions and high pressure (quartz, $alpha$-cristobalite, $beta$-cristobalite) and of models of amorphous silica (containig up to 162 atoms), obtained by quenching from the melt in combined classical and Car-Parrinello molecular dynamics simulations. The computational framework has also been checked on the photoelastic tensor of crystalline silicon and MgO as prototypes of covalent and ionic systems. The agreement with available experimental data is good. A phenomenological model suitable to describe the photoelastic properties of different silica polymorphs is devised by fitting on the ab-initio data.
We have performed first-principles calculations of thick slabs of Ti-doped sodium alanate (NaAlH_4), which allows to study the system energetics as the dopant progresses from the surface to the bulk. Our calculations predict that Ti stays on the surface, substitutes for Na, and attracts a large number of H atoms to its vicinity. Molecular dynamics simulations suggest that the most likely product of the Ti-doping is the formation of H-rich TiAl_n (n>1) compounds on the surface, and hint at the mechanism by which Ti enhances the reaction kinetics of NaAlH_4.
Using molecular dynamics simulations we investigate the dependence of the structural and vibrational properties of the surfaces of sodo-silicate glasses on the sodium content as well as the nature of the surface. Two types of glass surfaces are considered: A melt-formed surface (MS) in which a liquid with a free surface has been cooled down into the glass phase and a fracture surface (FS) obtained by tensile loading of a glass sample. We find that the MS is more abundant in Na and non-bridging oxygen atoms than the FS and the bulk glass, whereas the FS has higher concentration of structural defects such as two-membered rings and under-coordinated Si than the MS. We associate these structural differences to the production histories of the glasses and the mobility of the Na ions. It is also found that for Na-poor systems the fluctuations in composition and local atomic charge density decay with a power-law as a function of distance from the surface while Na-rich systems show an exponential decay with a typical decay length of $approx2.3$~AA. The vibrational density of states shows that the presence of the surfaces leads to a decrease of the characteristic frequencies in the system. The two-membered rings give rise to a pronounce band at $approx880$~cm$^{-1}$ which is in good agreement experimental observations.
Calcium silicate perovskite (CaSiO$_3$) is one of the major mineral components of the lower mantle, but has been the subject of relatively little work compared to the more abundant Mg-based materials. One of the major problems related to CaSiO$_3$ that is still the subject of research is its crystal structure under lower mantle conditions - a cubic Pm$bar{3}$m structure is accepted in general, but some have suggested that lower-symmetry structures may be relevant. In this paper, we use a fully first-principles vibrational self-consistent field method to perform high accuracy anharmonic vibrational calculations on several candidate structures at a variety of points along the geotherm near the base of the lower mantle to investigate the stability of the cubic structure and related distorted structures. Our results show that the cubic structure is the most stable throughout the lower mantle, and that this result is robust against the effects of thermal expansion.