Engineering Papers⌕ Search

Engineering topics

Tucker, Matthew G.

Publications and source records attributed to Tucker, Matthew G..

Integrating machine learning interatomic potentials with hybrid reverse Monte Carlo structure refinements in RMCProfile

Structure refinement with reverse Monte Carlo (RMC) is a powerful tool for interpreting experimental diffraction data. To ensure that the under-constrained RMC algorithm yields reasonable results, the hybrid RMC approach applies interatomic potentials to obtain solutions that are both physically sensible and in agreement with experiment. To expand the range of materials that can be studied with hybrid RMC, we have implemented a new interatomic potential constraint in RMCProfile that grants flexibility to apply potentials supported by the Large-scale Atomic/Molecular Massively Parallel Simulator ( LAMMPS ) molecular dynamics code. This includes machine learning interatomic potentials, which provide a pathway to applying hybrid RMC to materials without currently available interatomic potentials. To this end, we present a methodology to use RMC to train machine learning interatomic potentials for hybrid RMC applications.

Cuillier, Paul↗

Atomic and Electronic Structure in MgO–SiO 2

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–SiO 2 based on a comparison with the crystalline topology. We find that packing of oxygen atoms in Mg 2 SiO 4 is larger than that in MgSiO 3 , 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-) Mg 2 SiO 4 is the signature of low glass-forming ability (GFA), and high GFA g-(l-) MgSiO 3 shows a unique glass topology, which is different from c-MgSiO 3 . 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–SiO 2 , suggesting that electronic structure does not play an important role to determine GFA. We finally concluded the GFA of MgO–SiO 2 binary is dominated by the atomic structure in terms of network topology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Achieving giant electrostrain of above 1% in (Bi,Na)TiO 3 -based lead-free piezoelectrics via introducing oxygen-defect composition

Piezoelectric ceramics have been extensively used in actuators, where the magnitude of electrostrain is key indicator for large-stroke actuation applications. Here, we propose an innovative strategy based on defect chemistry to form a defect-engineered morphotropic phase boundary and achieve a giant strain of 1.12% in lead-free Bi 0.5 Na 0.5 TiO 3 (BNT)–based ceramics. The incorporation of the hypothetical perovskite BaAlO 2.5 with nominal oxygen defect into BNT will form strongly polarized directional defect dipoles, leading to a strong pinning effect after aging. The large asymmetrical strain is mainly attributed to two factors: The defect dipoles along crystallographic [001] direction destroy the long-range ordering of the ferroelectric and activate a reversible phase transition while promoting polarization rotation when the dipoles are aligned along the applied electric field. Our results not only demonstrate the potential application of BNT-based materials in low-frequency, large-stroke actuators but also provide a general methodology to achieve large strain.

36 MATERIALS SCIENCE↗

Lorentz factor for time-of-flight neutron Bragg and total scattering

We report the three fundamental origins of the Lorentz factor for neutron time-of-flight powder diffraction are revisited. A detailed derivation of the Lorentz factor is presented in the context of diffuse scattering modelling in reciprocal space when perfect periodicity is assumed, and the total scattering pattern is constructed in its discrete form – the factor in this case becomes 1/Q 2 (or d 2 ). Discussion is also presented with respect to practical data reduction where a vanadium measurement is usually taken as the normalization factor (to account for various factors such as detector efficiency), and it is shown that the existence of the Lorentz factor is independent of such a normalization process.

36 MATERIALS SCIENCE↗

In situ inelastic neutron scattering of mixed CH 4 –CO 2 hydrates

An abundant source of CH 4 can be found in natural hydrate deposits. Recent demonstration of CH 4 recovery from hydrates via CO 2 exchange has revealed the potential as a fuel source that also provides a medium for carbon sequestration. It is vital to understand the structural and dynamic impacts of guest variation in CH 4 , CO 2 , and mixed hydrates and link the results to the stability of various deposits in nature, harvesting methane, and sequestering CO 2 . Molecular vibrations are examined in CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates at 5 and 190 K and Xe hydrates for comparison. Inelastic neutron scattering (INS) is an ideal spectroscopy technique to observe the dynamic modes in the hydrate structure and enclathrated CH 4 , as it is extremely sensitive to 1 H. The presence of CO 2 in hydrates tightens the lattice. It introduces more active librational modes to the host lattice, while hindering the motion of CH 4 in mixed CH 4 -CO 2 hydrate at 5 K. At 190 K, a large broadening of the CH 4 librational modes indicates disorder in the structure leading to dissociation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural analyses of amorphous calcium carbonate before and after removing strontium ions from an aqueous solution

Amorphous calcium carbonate (ACC) is a precursor of crystalline calcium carbonate; hence, its structural information at the atomic level is very important for controlling the morphology of crystalline calcium carbonate. In this study, we attempted to elucidate the process of Sr extraction from aqueous solution by using ACC for the purpose of removal of radioactive Sr from the contaminated water leaked after the Fukushima Daiichi nuclear accident. Here, the pair distribution functions, g(r) obtained by X-ray and neutron diffraction (ND) measurements show that ACC has a structure partially similar to that of monohydrocalcite, suggesting that ACC is transferred to the crystalline calcium carbonate starting from its crystal nucleus. Rietveld analysis of the ND data showed that the ACC that removed Sr was crystallized to calcite. However, the Sr–O coordination analyzed using extended X-ray absorption fine structure (EXAFS) implies that the local environment of O around Sr is similar to that in crystalline calcium carbonate aragonite.

36 MATERIALS SCIENCE↗

Intermolecular correlations of liquid and glassy CS 2 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. Here, 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 CS 2 , 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 CS 2 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 CS 2 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 CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Broad Distribution of Local Polar States Generates Large Electrothermal Properties in Pb-Free Relaxor Ferroelectrics

Electrothermal energy conversion provides attractive solutions for global energy management, such as energy harvesting from waste heat using pyroelectric energy conversion (PEC) and efficient cooling of portable electronics or data servers using the electrocaloric effect. Relaxor ferroelectrics are attractive for electrothermal energy conversion because of their large pyroelectric coefficients over a wide temperature range. Although Pb-based relaxors are well-known, toxicity concerns have mandated the intense search for Pb-free alternatives. Here, we engineered (Ba,Ca)TiO 3 -based relaxors based on a multisite doping strategy, which show promising electrothermal performance, viz. a maximum PEC efficiency of 14% and electrocaloric refrigeration capacity of 115 J/kg. Using local-scale structural analysis, we provide an atomistic model for large electrothermal properties in the newly designed Pb-free ferroelectrics, whereby a temperature-independent continuous distribution of cation displacement directions creates easy pathways for microscopic polarization reorientation. This research provides key structural insight for future atomic-scale engineering of environmentally sustainable ferroelectrics in energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dual Orbital Degeneracy Lifting in a Strongly Correlated Electron System

The local structure of NaTiSi 2 O 6 is examined across its Ti-dimerization orbital-assisted Peierls transition at 210 K. An atomic pair distribution function approach evidences local symmetry breaking preexisting far above the transition. The analysis unravels that, on warming, the dimers evolve into a short range orbital degeneracy lifted (ODL) state of dual orbital character, persisting up to at least 490 K. The ODL state is correlated over the length scale spanning ~ 6 sites of the Ti zigzag chains. Results imply that the ODL phenomenology extends to strongly correlated electron systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A cloud platform for atomic pair distribution function analysis: PDFitc

A cloud web platform for analysis and interpretation of atomic pair distribution function (PDF) data ( PDFitc ) is described. The platform is able to host applications for PDF analysis to help researchers study the local and nanoscale structure of nanostructured materials. The applications are designed to be powerful and easy to use and can, and will, be extended over time through community adoption and development. The currently available PDF analysis applications, structureMining, spacegroupMining and similarityMapping , are described. In the first and second the user uploads a single PDF and the application returns a list of best-fit candidate structures, and the most likely space group of the underlying structure, respectively. In the third, the user can upload a set of measured or calculated PDFs and the application returns a matrix of Pearson correlations, allowing assessment of the similarity between different data sets. structureMining is presented here as an example to show the easy-to-use workflow on PDFitc . In the future, as well as using the PDFitc applications for data analysis, it is hoped that the community will contribute their own codes and software to the platform.

97 MATHEMATICS AND COMPUTING↗

Structure and properties of densified silica glass: characterizing the order within disorder

The broken symmetry in the atomic-scale ordering of glassy versus crystalline solids leads to a daunting challenge to provide suitable metrics for describing the order within disorder, especially on length scales beyond the nearest neighbor that are characterized by rich structural complexity. Here, we address this challenge for silica, a canonical network-forming glass, by using hot versus cold compression to (i) systematically increase the structural ordering after densification and (ii) prepare two glasses with the same high-density but contrasting structures. The structure was measured by high-energy X-ray and neutron diffraction, and atomistic models were generated that reproduce the experimental results. The vibrational and thermodynamic properties of the glasses were probed by using inelastic neutron scattering and calorimetry, respectively. Traditional measures of amorphous structures show relatively subtle changes upon compacting the glass. The method of persistent homology identifies, however, distinct features in the network topology that change as the initially open structure of the glass is collapsed. The results for the same high-density glasses show that the nature of structural disorder does impact the heat capacity and boson peak in the low-frequency dynamical spectra. Densification is discussed in terms of the loss of locally favored tetrahedral structures comprising oxygen-decorated SiSi 4<:sub> tetrahedra.

36 MATERIALS SCIENCE↗

Molecular Rotational Dynamics in Mixed CH 4 –CO 2 Hydrates: Insights from Molecular Dynamics Simulations

The impact of guest molecule composition on the rotational dynamics in CH 4 , CO 2 , and mixed CH 4 –CO 2 gas hydrates is investigated with classical molecular dynamics simulations. Rotational autocorrelation functions are calculated for the guest and host molecules in each hydrate composition from simulation trajectories at 10, 40, 190, and 270 K. Analysis of these functions for each molecule is further decomposed into cage type for each CH 4 and CO 2 guest and cage face for the H 2 O host. CH 4 becomes more constrained, and CO 2 gains freedom in the mixed guest systems. Mixing guest species in gas hydrates alters the intermolecular interaction environment, impacting the rotational motion of the guest molecules. This effect is also seen in the host lattice H 2 O molecules, as evidence of molecular rotations is seen in the CO 2 hydrate at 270 K during longer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reverse Monte Carlo modeling for low-dimensional systems

Reverse Monte Carlo (RMC) is one of the commonly used approaches for modeling total scattering data. However, to extend the capability of the RMC method for refining the structure of nanomaterials, the dimensionality and finite size need to be considered when calculating the pair distribution function (PDF). To achieve this, the simulation box must be set up to remove the periodic boundary condition in one, two or three of the dimensions. This then requires a correction to be applied for the difference in number density between the real system and the simulation box. In certain circumstances an analytical correction for the uncorrelated pairings of atoms is also applied. The validity and applicability of our methodology is demonstrated by applying the algorithms to simulate the PDF patterns of carbon systems with various dimensions, and also by using them to fit experimental data of CuO nanoparticles. This alternative approach for characterizing the local structure of nano-systems with the total scattering technique will be made available via the RMCProfile package. The theoretical formulation and detailed explanation of the analytical corrections for low-dimensional systems – 2D nanosheets, 1D nanowires and 0D nanoparticles – is also given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Local magnetic cluster size identified by neutron total scattering in the site-diluted spin glass S n x F e 4 - x N ( x = 0.88 )

A detailed structure analysis for the site-diluted Sn x Fe 4- x N ( x = 0.25, 0.41, and 0.88) has been carried out through complex modeling of the neutron total scattering data. We present quantitative evidence showing the local ferromagnetic cluster size extending to ~8 Å on average when Sn 0.88 Fe 3.12 N undergoes the spin-glass transition (the other two not showing such a transition). In conclusion, the modeling methodology used in this work involving the corefinement of the nuclear and magnetic structure in both real and reciprocal space can potentially be applied generally to explore a variety of spin-glass material problems.

36 MATERIALS SCIENCE↗