Atomic and Electronic Structure in MgO-SiO2.

Understanding disordered structure is difficult due to insufficient information in experimental data. Here, we overcome this issue by using a combination of diffraction and simulation to investigate oxygen packing and network topology in glassy (g-) and liquid (l-) MgO-SiO2 based on a comparison with the crystalline topology. We find that packing of oxygen atoms in Mg2SiO4 is larger than that in MgSiO3, and that of the glasses is larger than that of the liquids. Moreover, topological analysis suggests that topological similarity between crystalline (c)- and g-(l-) Mg2SiO4 is the signature of low glass-forming ability (GFA), and high GFA g-(l-) MgSiO3 shows a unique glass topology, which is different from c-MgSiO3. We also find that the lowest unoccupied molecular orbital (LUMO) is a free electron-like state at a void site of magnesium atom arising from decreased oxygen coordination, which is far away from crystalline oxides in which LUMO is occupied by oxygen's 3s orbital state in g- and l-MgO-SiO2, suggesting that electronic structure does not play an important role to determine GFA. We finally concluded the GFA of MgO-SiO2 binary is dominated by the atomic structure in terms of network topology.

Intermolecular correlations of liquid and glassy CS2 studied by synchrotron radiation x-ray diffraction.

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. We have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS2, which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16-25.7 Å-1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS2 were also recorded in a wide temperature range of 160-300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS2 crystal, indicating that the structure of glassy CS2 is close to that of crystalline CS2.

Illustrated formalisms for total scattering data: a guide for new practitioners.

The total scattering method is the simultaneous study of both the real- and reciprocal-space representations of diffraction data. While conventional Bragg-scattering analysis (employing methods such as Rietveld refinement) provides insight into the average structure of the material, pair distribution function (PDF) analysis allows for a more focused study of the local atomic arrangement of a material. Generically speaking, a PDF is generated by Fourier transforming the total measured reciprocal-space diffraction data (Bragg and diffuse) into a real-space representation. However, the details of the transformation employed and, by consequence, the resultant appearance and weighting of the real-space representation of the system can vary between different research communities. As the worldwide total scattering community continues to grow, these subtle differences in nomenclature and data representation have led to conflicting and confusing descriptions of how the PDF is defined and calculated. This paper provides a consistent derivation of many of these different forms of the PDF and the transformations required to bridge between them. Some general considerations and advice for total scattering practitioners in selecting and defining the appropriate choice of PDF in their own research are presented. This contribution aims to benefit people starting in the field or trying to compare their results with those of other researchers.

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions.

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson's theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs-Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born-Oppenheimer molecular dynamics performed at a single temperature. This advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

Local structure and distortions of mixed methane-carbon dioxide hydrates.

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

A high temperature gas flow environment for neutron total scattering studies of complex materials.

We present the design and capabilities of a high temperature gas flow environment for neutron diffraction and pair distribution function studies available at the Nanoscale Ordered Materials Diffractometer instrument at the Spallation Neutron Source. Design considerations for successful total scattering studies are discussed, and guidance for planning experiments, preparing samples, and correcting and reducing data is defined. The new capabilities are demonstrated with an in situ decomposition study of a battery electrode material under inert gas flow and an in operando carbonation/decarbonation experiment under reactive gas flow. This capability will aid in identifying and quantifying the atomistic configurations of chemically reactive species and their influence on underlying crystal structures. Furthermore, studies of reaction kinetics and growth pathways in a wide variety of functional materials can be performed across a range of length scales spanning the atomic to the nanoscale.

Molecular Dynamics Simulations of Hydration Effects on Solvation, Diffusivity, and Permeability in Chitosan/Chitin Films.

The effects of hydration on the solvation, diffusivity, solubility, and permeability of oxygen molecules in sustainable, biodegradable chitosan/chitin food packaging films were studied via molecular dynamics and confined random walk simulations. With increasing hydration, the membrane has a more homogeneous water distribution with the polymer chains being fully solvated. The diffusivity increased by a factor of 4 for oxygen molecules and by an order of magnitude for water with increasing the humidity. To calculate the Henry's constant and solubility of oxygen in the membranes with changing hydration, the excess chemical potential was calculated via free energy perturbation, thermodynamic integration and direct particle deletion methods. The simulations predicted a higher solubility and permeability for the lower humidity, in contradiction to experimental results. All three methods for calculating the solubility were in good agreement. It was found that the Coulombic interactions in the potential caused the oxygen to bind too strongly to the protonated amine group. Insight from this work will help guide molecular modeling of chitosan/chitin membranes, specifically permeability measurements for small solute molecules. Efforts to chemically tailor chitosan/chitin membranes to favor discrete as opposed to continuous aqueous domains could reduce oxygen permeability.

Ab Initio Molecular Dynamics Simulations of an Excess Proton in a Triethylene Glycol-Water Solution: Solvation Structure, Mechanism, and Kinetics.

We investigate the solvation shell structures, the distribution of protonic defects, mechanistic details, kinetics, and dynamics of proton transfer for an excess proton in bulk water and for an excess proton in an aqueous solution of triethylene glycol (TEG) via Car-Parrinello molecular dynamics simulations. The PW91, PBE, and PBE with the Tkatchenko-Scheffler (TS) density-dependent dispersion functionals were used and compared for bulk water and the TEG-water mixtures. The excess proton is found to reside predominantly on water molecules but also resides on hydroxyl groups of TEG. The lifetimes associated with structural diffusion time scales of the protonated water were found to be on the order of ∼1 ps. All three functionals studied support the presolvation requirement for structural diffusion. The highest level of theory shows a reduction in the free energy barrier for water-water proton transfer in TEG-water mixtures compared to bulk water. The effect of TEG shows no strong change in the kinetics for TEG-water mixtures compared to bulk water for this same level of theory. The excess proton displays burst-rest behavior in the presence of TEG, similar to that found in bulk water. We find that the TEG chain disrupts the hydrogen-bond network, causing the solvation shell around water to be populated by TEG chain groups instead of other waters, reducing the rigidity of the hydrogen-bond network. Methylene is a dominant hydrogen bond donor for the protonated water in hydrogen-bond networks associated with proton transfer and structural diffusion. This is consistent with previous studies that have found the hydronium ion to be amphiphilic in nature and to have higher proton mobility at oil-water interfaces.