Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic ordering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deconstructing the Local Intermolecular Ordering and Dynamics of Liquid Chloroform and Bromoform

Local intermolecular structure and dynamics of the polar molecular liquids chloroform and bromoform are studied by molecular dynamics simulation. Structural distribution functions, including 1- and 2-D pair correlations and dipole contour plots allow direct comparison and show agreement with recent analyses of diffraction experiments. Studies of the haloforms’ reorientational dynamics and longevity of structural features resulting from intermolecular interaction extend previous work toward deeper understanding of the factors controlling these features. Analyses of ensemble average structures and dynamical properties isolate mass, electrostatics, and steric packing as driving forces or contributing factors for the observed ordering and dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bottlebrush Block Copolymer Thin Films

The self-assembly of block copolymers (BCPs) in thin films is governed by interfacial interactions (enthalpy) between the blocks and interfaces (at both the substrate and the surface). Advantage can also be taken of entropy to control the polymer orientation. Here, we synthesized a series of bottlebrush block copolymers (BBCPs), with deuterated polystyrene (DPS) as the core block and poly­(solketal acrylate) (PSA) as the corona block, where the backbone length (N BB) and grafting densities (GDs) were varied. The hydrophobic PSA block was converted to a hydrophilic poly­(glyceryl acrylate) (PGA) block by solid-state hydrolysis, bringing the BBCP from a disordered state into an ordered state with a lamellar microdomain morphology. The orientation of the morphology as a function of distance from the interfaces was systematically investigated by atomic force microscopy (AFM), interference microscopy, grazing-incidence small-angle X-ray scattering (GISAXS), and grazing-incidence small angle neutron scattering (GISANS). For N BB = 1 (diblock BCP), a mixed lamellar orientation was found. For N BB = 5 (star-like BBCP), a vertical orientation of the lamellar microdomains was found dominantly at the interfaces that propagated into the bulk of the film due to the entropic penalty associated with constraints on the junction points. Increasing the N BB to 50 (rod-like BBCP) and 100 (worm-like BBCP) yielded a random orientation, as the side chains became overcrowded when aligned in the same direction. By reducing the grafting density (GD) of a N BB = 50 BBCP from 100% to 75% and 50% by the incorporation of spacers into the backbone, the steric hindrance of side chains is reduced markedly, facilitating a vertical orientation of the lamellar microdomains. However, the vertical orientation becomes more random within the bulk of the thin film for different BBCP architectures. By modifying the substrate with a layer of PS, the preferential interactions between the core block of the BBCP and the interfaces induce a reconfiguration of the BBCP, drawing the core block out to contact substrate directly driven by enthalpy. Under these conditions, N BB = 5 showed a predominantly parallel orientation of the lamellar microdomains with the corona block minimizing contact with the substrate. Higher N BB values of 50 and 100 showed a more random orientation of the lamellar microdomains at different depths on the PS-modified substrate. Overall, thin-film morphology was regulated by tuning the configuration of the BBCP, varying the entropic contribution to the microdomain orientation.

Chen, Zhan↗

Epitaxial Stabilization and Oxygen Evolution Reaction Activity of Metastable Columbite Iridium Oxide

Non-rutile polymorphs of binary iridium oxide such as columbite IrO 2 (α-IrO 2 ) are promising candidates for highly active acid-stable oxygen evolution reaction (OER) catalysts, yet their synthesis has been challenging due to the dominant thermodynamic stability of rutile IrO 2 (R-IrO 2 ). Here, we report the growth of α-IrO 2 via epitaxial thin films using pulsed laser deposition. We observe that, in competition with R-IrO 2 (100), the films can be optimized to be predominantly (100)-oriented α-IrO 2 . Surprisingly, the activity of α-IrO 2 shows a large discrepancy of ~0.2 V in the overpotential compared to its predicted activity, which is resolved via theoretical calculations to be a crystal orientation effect. This demonstrates that the electrocatalytic activity can be significantly varied upon crystal orientation, a parameter that is difficult to control in conventional polycrystalline systems but accessible in epitaxial thin films. In total, this study demonstrates epitaxial thin film growth as a powerful technique, which can overcome large energetic instabilities on the order of ~300 meV to stabilize metastable material structures inaccessible by bulk synthesis. As a result, this provides unique opportunities to effectively identify the atomic structure of active catalysts by combining investigations of metastable materials with theoretical predictions.

36 MATERIALS SCIENCE↗

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗

Magnetic metamaterials by ion-implantation

We present a method for the additive fabrication of planar magnetic nanoarrays with minimal surface roughness. Synthesis is accomplished by combining electron-beam lithography, used to generate nanometric patterned masks, with ion implantation in thin films. By implanting 56 Fe + ions, we are able to introduce magnetic functionality in a controlled manner into continuous Pd thin films, achieving 3D spatial resolution down to a few tens of nanometers. Our results demonstrate the application of this technique in fabricating square artificial spin ice lattices, which exhibit well-defined magnetization textures and interactions among the patterned magnetic elements.

36 MATERIALS SCIENCE↗

Evaluation of current density measurement methods for high-current, low-beta electron beams

Recent current density measurements with a low-beta electron beam (β=0.5-0.75) are strongly affected by electron scatter and Cherenkov limits. These invasive measurements are heavily influenced by the electron energy and the intercepting material properties: density, index of refraction, and atomic number. We evaluate various measurement methods, including x-ray scintillation and Cherenkov emission. To optimize the x-ray scintillation technique for this energy regime, metal foils of differing atomic number (Z) and thickness were placed upstream of a scintillation screen. We have selected Cherenkov emitters, based on a material refractive index, in order to differentiate the contributions from electron scatter and produced fluorescence, Cherenkov emission, and total internal reflection. The observed distributions for both approaches are verified through MCNP6®. Additionally, the simulations provide insight as to what physical processes dominate the measured distribution. We conclude by determining the limits and optimal measurement range for each measurement technique.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗

Crystal and magnetic structures of magnetic topological insulators MnBi 2 Te 4 and MnBi 4 Te 7

Using single-crystal neutron diffraction, we present a systematic investigation of the crystal structure and magnetism of the van der Waals topological insulators MnBi 2 Te 4 and MnBi 4 Te 7 , where rich topological quantum states have been recently predicted and observed. Structural refinements reveal that considerable Bi atoms occupied on the Mn sites in both materials, distinct from the previously reported antisite disorder. We show unambiguously that MnBi 2 Te 4 orders antiferromagnetically below 24 K, featured by a magnetic symmetry R I - 3 c , while MnBi 4 Te 7 is antiferromagnetic below 13 K with a magnetic space group P c - 3 c 1 . They both present antiferromagnetically coupled ferromagnetic layers with spins along the c axis. We put forward a stacking rule for the crystal structure of an infinitely adaptive series MnBi 2 n Te 3 n + 1 ( n ≥ 1 ) with a building unit of [ Bi 2 Te 3 ] . Overall, by comparing the magnetic properties between MnBi 2 Te 4 and MnBi 4 Te 7 , together with recent density-functional theory calculations, we concluded that a two-dimensional magnetism limit might be realized in the derivatives. Our work may promote theoretical studies of topological magnetic states in the series of MnBi 2 n Te 3 n + 1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

4 d element induced improvement of structural disorder and development of weakly reentrant spin-glass behavior in NiRuMnSn

The pursuit of efficient spin polarization in quaternary Heusler alloys with the general formula XX'YZ (where X, X', and Y are transition metals and Z is a p-block element) has been a subject of significant scientific interest. While previous studies showed that isoelectronic substitution of a 4d element in place of a 3d element in quaternary Heusler alloys improves the half-metallic ferromagnetic characteristics, our research here, on the quaternary Heusler alloy NiRuMnSn, reveals a strikingly different scenario. In this study, we present a detailed structural analysis of the material using x-ray absorption fine structure and neutron diffraction techniques which confirms the formation of a single-phase compound with 50:50 site disorder between Ni and Ru atoms at 4c and 4d sites. Contrary to expectations, our density functional theory calculations suggest a considerable decrease in spin polarization even in the ordered structure. Additionally, we report on the compound's exceptional behavior, displaying a rare reentrant spin-glass property below ~ 60 K, a unique and intriguing feature for quaternary Heusler-type compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗

CEGANN: CRYSTAL EDGE GRAPH ATTENTION NEURAL NETWORK

SF-22-156 Machine learning (ML) models and applications in materials design and discovery typically involve the use of feature representations or descriptors followed by a learning algorithm that maps them to user desired properties of interest. Most popular mathematical formulation-based descriptors are not unique across atomic environments and suffer from transferability issues across different application domains and/or material classes. The CEGANN code provides a unified interface to facilitate material characterization across materials across multiple scales (from atomic to mesoscale) and diverse classes of materials ranging from metals oxides, non-metals, and even hierarchical materials such as zeolites and semi ordered materials such as mesophases. CEGANN implements a Graph Attention Network (GAT) type convolution architecture. The details of network architecture can be found in the paper https://doi.org/10.48550/arXiv.2207.10168. The software comes with pretrained examples and dataset for the classification of the following representative systems: (1) Structure-level representation such as space group (2) Structural dimensionality (e.g., bulk, 2D, clusters etc.) (3) Grain boundary identification (4) Nucleation and growth of a zeolite polymorph (5) Characterization of binary mesophases and their phase transitions (6) Growth of ice. The code is written in python programming language.

CHAN, HENRYT↗

Machine Learning Predictions of Transition Probabilities in Atomic Spectra

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

74 ATOMIC AND MOLECULAR PHYSICS↗

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Finite-volume pionless effective field theory for few-nucleon systems with differentiable programming

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wave functions optimized using differentiable programming and via solution of a generalized eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wave functions, yielding comparably accurate representations of the ground-state wave functions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. Further, the method is demonstrated through calculations of the binding energies of nuclei with atomic number A ϵ {2,3,4} in finite volume, matched to lattice QCD calculations at quark masses corresponding to m π = 806 MeV, and infinite-volume effective field theory calculations of A ϵ {2,3,4,5,6} systems based on this matching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamical phase transitions in the collisionless pre-thermal states of isolated quantum systems: theory and experiments

Abstract We overview the concept of dynamical phase transitions (DPTs) in isolated quantum systems quenched out of equilibrium. We focus on non-equilibrium transitions characterized by an order parameter, which features qualitatively distinct temporal behavior on the two sides of a certain dynamical critical point. DPTs are currently mostly understood as long-lived prethermal phenomena in a regime where inelastic collisions are incapable to thermalize the system. The latter enables the dynamics to substain phases that explicitly break detailed balance and therefore cannot be encompassed by traditional thermodynamics. Our presentation covers both cold atoms as well as condensed matter systems. We revisit a broad plethora of platforms exhibiting pre-thermal DPTs, which become theoretically tractable in a certain limit, such as for a large number of particles, large number of order parameter components, or large spatial dimension. The systems we explore include, among others, quantum magnets with collective interactions, ϕ 4 quantum field theories, and Fermi–Hubbard models. A section dedicated to experimental explorations of DPTs in condensed matter and AMO systems connects this large variety of theoretical models.

Physics↗

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗