Effects of Na/K-Cl Salts on Hydrolysis of Aluminosilicate Glass Using Ab Initio Molecular Dynamics

The structural and chemical modifications on the surface of pure and alkali-doped aluminosilicate (AS) glasses due to hydrolysis are investigated using ab initio molecular dynamics. The effects of water on the glass network are fully elucidated by analyzing the short- and intermediate-range structural orders embedded in the pair distribution function, bond length and angle distribution, coordination number, and interatomic bonding. A novel concept of total bond order is used to quantify and compare the strength of bonds in hydrated and unhydrated glasses. We show that AS glass is hydrolyzed by water diffusion near the surface and by proton (H+) transfers into the bulk, which increases with time. Hence, a dissolved glass-water interface becomes rich in Si-OH and Al-OH. The alkali ions associated with the nonbridging oxygen accelerate the hydrolysis by facilitating water and H+ diffusion. Al is more impacted by hydrolysis than Si, resulting in greater variation in the Al-O bond order than Si-O. Doping of NaCl and KCl enhances the ionization of water and the hydrolysis of ASs with increased salt concentration. The KCl doping ionizes more water molecules and causes more degradation of the glass network than NaCl. Co-doping of Na and K results in a mixed alkali effect due to complex interatomic bonding from different-sized ions. These exceptionally detailed findings in highly complex glasses with varying salt compositions provide new and unprecedented atomistic insights that can help to understand the hydrolysis and dissolution mechanisms of ASs and other silicate glasses.

A Thorough Investigation of Electronic, Optical, Mechanical, and Thermodynamic Properties of Stable Glasslike Sodium Germanate under Compressive Hydrostatic Pressure: Ab Initio Study

In this paper, we have tried to elucidate the variation of structural, electronic, and thermodynamic properties of glasslike Na₂GeO₃ under compressive isotropic pressure within a framework of density functional theory (DFT). The result shows stable structural (orthorhombic → tetragonal) and electronic (indirect → direct) phase transitions at P ∼ 20 GPa. The electronic band gap transition plays a key role in the enhancement of optical properties. The results of the thermodynamic properties have shown that Na₂GeO₃ follows Debye's low-temperature specific heat law and the classical thermodynamic of the Dulong-Petit law at high temperature. The pressure sensitivity of the electronic properties led us to compute the piezoelectric tensor (both in relaxed and clamped ions). We have observed significant electric responses in the form of a piezoelectric coefficient under applied pressure. This property suggested that Na₂GeO₃ could be a potential material for energy harvest in future energy-efficient devices. As expected, Na₂GeO₃ becomes harder and harder under compressive pressure up to the phase transition pressure (∼20 GPa) which can be read from Pugh's ratio (k_H) > 1.75, however, at pressures above 20 GPa k_H < 1.75, which may be due to the formation of fractures at high pressure.

The Atomic Structure and Mechanical Properties of ZIF-4 under High Pressure: Ab Initio Calculations

The effects of pressure on the structural and electronic properties and the ionic configuration of ZIF-4 were investigated through the first-principles method based on the density functional theory. The elastic properties, including the isotropic bulk modulus K, shear modulus G, Young's modulus E, and Poisson's ratio ν of the orthorhombic-type structure ZIF-4 were determined using the Voigt-Reuss-Hill averaging scheme. The results show that the ZIF-4 phase is ductile according to the analysis of K/G and Cauchy pressure. The Debye temperatures obtained from the elastic stiffness constants increase with increasing pressure. Finally, the pressure-dependent behaviors of the density of states and ionic configuration are successfully calculated and discussed.

Structural Evolution of SiO_{2} Glass with Si Coordination Number Greater than 6

Pair distribution function measurement of SiO_{2} glass up to 120 GPa reveals changes in the first-, second-, and third-neighbor distances associated with an increase in Si coordination number C_{Si} to >6 above 95 GPa. Packing fractions of Si and O determined from the first- and second-neighbor distances show marked changes accompanied with the structural evolution from C_{Si}=6 to >6. Structural constraints in terms of ionic radius ratio of Si and O, and ratio of nonbonded radius to bonded Si─O distance support the structural evolution of SiO_{2} glass with C_{Si}>6 at high pressures.

Interfacial Interaction between Suolunite Crystal and Silica Binding Peptide for Novel Bioinspired Cement

Cement and concrete have been important construction materials throughout human history. There is an urgent need to explore novel and untraditional cementitious materials to enhance the durability of building materials and structures in response to increased infrastructure demand worldwide. We report an exploratory study on a biocomposite cement based on a large-scale computational study using density functional theory. An explicitly solvated mixture of a mineral calcium silicate hydrate (C-S-H) crystal suolunite (Ca₂Si₂O₅(OH)₂·H₂O) and a silicon binding peptide with amino acid sequence PRO-PRO-PRO-TRP-LEU-PRO-TYR-MET-PRO-PRO-TRP-SER is constructed using ab initio molecular dynamics (AIMD). Detailed analysis on the interface structure, interatomic bonding, mechanical properties, and solvent effect of this model reveals a complex interplay of different types of covalent and ionic bonding, including ubiquitous hydrogen bonding which plays a crucial role in their properties. The use of the total bond order density (TBOD), a single quantum mechanical metric, for assessing the interfacial cohesion for this composite biocement is proposed. We find that the solvated model has a slightly larger TBOD than the dried one. These results could lead to a systematic search and rational design for different types of bioinspired and hybrid functional materials with other inorganic minerals and organic peptides.

Novel pressure-induced topological phase transitions of supercooled liquid and amorphous silicene

This molecular dynamics (MD) simulation carries a detailed analysis of a pressure-induced structural transition supercooled liquid and amorphous silicene (a-silicene). Low-density models of supercooled liquid and a-silicene containing 10 000 atoms are obtained by rapid cooling processes from the melts. Then, an a-silicene model at T  =  1000 K, a supercooled liquid model at T  =  1500 K and a liquid silicon model at T  =  2000 K have been isothermally compressed step by step up to a high density in order to observe the pressure-induced structural changes. Specifically 'Cairo tiling' pentagonal and square lattices of silicene are discovered in our calculations. Structural properties of those penta-silicene and tetra-silicene models have been carefully analyzed through the radial distribution functions, interatomic distances, bond-angle distributions under high-pressure condition. The dependence of pressure on formation behaviors is calculated via pressure-volume and energy-density relationships. The first order transition from low-density supercooled liquid/amorphous silicene to high-density penta-silicene and continuous transition from low-density liquid to high-density tetra-silicene are discussed. Atomic mechanism and sp³/sp² hybridization evolution are inspected whereas the role of low-membered ring defects/boundary promises remarkable application and advanced research in future.

Deformation behavior of an amorphous zeolitic imidazolate framework - from a supersoft material to a complex organometallic alloy

Zeolitic imidazolate frameworks (ZIFs)-a subset of metal-organic frameworks (MOFs)-have recently attracted immense attention. Many crystalline ZIFs (c-ZIFs) have highly porous zeolite structures that are ideal for molecular encapsulation. Recently emerging non-crystalline or amorphous ZIFs (a-ZIFs) with a similar short-range order are of interest because they can be converted from c-ZIFs for large-scale production. Here, we present a computational study of the deformation behavior of a unique a-ZIF model by simulating step-wise compression and expansion with strains between -0.389 and +0.376. An insulator-to-metal transition is observed at 51 GPa leading to a multicomponent light amorphous alloy of only 3.68 g (cm)-3. A high-density amorphous-to-amorphous phase transition is observed due to the sudden formation of N-N bond pairs. The systematic expansion of the a-ZIF retains the framework softness until it fractures at high strain. Based on the expansion data, we propose an empirical formula for super-soft materials, which is in line with available experimental data.

Complex interplay of interatomic bonding in a multi-component pyrophosphate crystal: K2Mg (H2P2O7)2·2H2O

The electronic structure and interatomic bonding of pyrophosphate crystal K₂Mg (H₂P₂O₇)₂·2H₂O are investigated for the first time showing complex interplay of different types of bindings. The existing structure from single-crystal X-ray diffraction is not sufficiently refined, resulting in unrealistic short O─H bonds which is rectified by high-precision density functional theory (DFT) calculation. K₂Mg (H₂P₂O₇)₂·2H₂O has a direct gap of 5.22 eV and a small electron effective mass of 0.14 m_e. Detailed bond analysis between every pair of atoms reveals the complexity of various covalent, ionic, hydrogen bonding and bridging bonding and their sensitive dependence on structural differences. The K--O bonds are much weaker than Mg--O bonds and contributions from the hydrogen bonds are non-negligible. Quantitative analysis of internal cohesion in terms of total bond order density and partial bond order density divulges the relative importance of different types of bonding. The calculated optical absorptions show multiple peaks and a sharp Plasmon peak at 23 eV and a refractive index of 1.44. The elastic and mechanical properties show features unique to this low-symmetry crystal. Phonon calculation gives vibrational frequencies in agreement with reported Raman spectrum. These results provide new insights indicating that acidic pyrophosphates could have a variety of unrealized applications in advanced technology.

Ab Initio Modeling of Structure and Properties of Single and Mixed Alkali Silicate Glasses

A density functional theory (DFT)-based ab initio molecular dynamics (AIMD) has been applied to simulate models of single and mixed alkali silicate glasses with two different molar concentrations of alkali oxides. The structural environments and spatial distributions of alkali ions in the 10 simulated models with 20% and 30% of Li, Na, K and equal proportions of Li-Na and Na-K are studied in detail for subtle variations among the models. Quantum mechanical calculations of electronic structures, interatomic bonding, and mechanical and optical properties are carried out for each of the models, and the results are compared with available experimental observation and other simulations. The calculated results are in good agreement with the experimental data. We have used the novel concept of using the total bond order density (TBOD), a quantum mechanical metric, to characterize internal cohesion in these glass models. The mixed alkali effect (MAE) is visible in the bulk mechanical properties but not obvious in other physical properties studied in this paper. We show that Li doping deviates from expected trend due to the much stronger Li-O bonding than those of Na and K doping. The approach used in this study is in contrast with current studies in alkali-doped silicate glasses based only on geometric characterizations.

Beyond sixfold coordinated Si in SiO2 glass at ultrahigh pressures

We investigated the structure of SiO2 glass up to 172 GPa using high-energy X-ray diffraction. The combination of a multichannel collimator with diamond anvil cells enabled the measurement of structural changes in silica glass with total X-ray diffraction to previously unachievable pressures. We show that SiO2 first undergoes a change in Si-O coordination number from fourfold to sixfold between 15 and 50 GPa, in agreement with previous investigations. Above 50 GPa, the estimated coordination number continuously increases from 6 to 6.8 at 172 GPa. Si-O bond length shows first an increase due to the fourfold to sixfold coordination change and then a smaller linear decrease up to 172 GPa. We reconcile the changes in relation to the oxygen-packing fraction, showing that oxygen packing decreases at ultrahigh pressures to accommodate the higher than sixfold Si-O coordination. These results give experimental insight into the structural changes of silicate glasses as analogue materials for silicate melts at ultrahigh pressures.

Ion irradiation induced structural modifications and increase in elastic modulus of silica based thin films

Ion irradiation is an alternative to heat treatment for transforming organic-inorganic thin films to a ceramic state. One major shortcoming in previous studies of ion-irradiated films is the assumption that constituent phases in ion-irradiated and heat-treated films are identical and that the ion irradiation effect is limited to changes in composition. In this study, we investigate the effects of ion irradiation on both the composition and structure of constituent phases and use the results to explain the measured elastic modulus of the films. The results indicated that the microstructure of the irradiated films consisted of carbon clusters within a silica matrix. It was found that carbon was present in a non-graphitic sp²-bonded configuration. It was also observed that ion irradiation caused a decrease in the Si-O-Si bond angle of silica, similar to the effects of applied pressure. A phase transformation from tetrahedrally bonded to octahedrally bonded silica was also observed. The results indicated the incorporation of carbon within the silica network. A combination of the decrease in Si-O-Si bond angle and an increase in the carbon incorporation within the silica network was found to be responsible for the increase in the elastic modulus of the films.

Computer simulations of glasses: the potential energy landscape

We review the current state of research on glasses, discussing the theoretical background and computational models employed to describe them. This article focuses on the use of the potential energy landscape (PEL) paradigm to account for the phenomenology of glassy systems, and the way in which it can be applied in simulations and the interpretation of their results. This article provides a broad overview of the rich phenomenology of glasses, followed by a summary of the theoretical frameworks developed to describe this phenomonology. We discuss the background of the PEL in detail, the onerous task of how to generate computer models of glasses, various methods of analysing numerical simulations, and the literature on the most commonly used model systems. Finally, we tackle the problem of how to distinguish a good glass former from a good crystal former from an analysis of the PEL. In summarising the state of the potential energy landscape picture, we develop the foundations for new theoretical methods that allow the ab initio prediction of the glass-forming ability of new materials by analysis of the PEL.

Temperature-dependent densification of sodium borosilicate glass

We provide a comprehensive understanding of the temperature-dependent changes in the network topology, structure, and properties of densified borosilicate glass.

Elastic moduli of permanently densified silica glasses

Modelling the mechanical response of silica glass is still challenging, due to the lack of knowledge concerning the elastic properties of intermediate states of densification. An extensive Brillouin Light Scattering study on permanently densified silica glasses after cold compression in diamond anvil cell has been carried out, in order to deduce the elastic properties of such glasses and to provide new insights concerning the densification process. From sound velocity measurements, we derive phenomenological laws linking the elastic moduli of silica glass as a function of its densification ratio. The found elastic moduli are in excellent agreement with the sparse data extracted from literature, and we show that they do not depend on the thermodynamic path taken during densification (room temperature or heating). We also demonstrate that the longitudinal sound velocity exhibits an anomalous behavior, displaying a minimum for a densification ratio of 5%, and highlight the fact that this anomaly has to be distinguished from the compressibility anomaly of a-SiO₂ in the elastic domain.

A molecular dynamics study of ambient and high pressure phases of silica: structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, α-PbO2-type, and pyrite-type - for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase.