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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Water-induced surface ordering facilitating the microcutting of ductile metals

Metal cutting is a crucial process in modern manufacturing. Enhancing the machinability of metals can significantly improve their production efficiency and surface integrity. Coating surface-active media (SAM) on the free surface of the metals before cutting is an easy method to improve machinability, which usually pertains to the category of the renowned Rehbinder effect. However, the existing SAM is usually hazardous and complex materials. Besides, the effect of SAM on the local structure of the metal surface remains unclear. In this study, water is employed as a simple yet often overlooked SAM in the microcutting of copper. Using water as SAM also allows the employment of X-ray absorption fine structure spectroscopy (XAFS) to study the local structure of copper with and without water coating. Results show that water coating on the free surface of copper can significantly reduce the cutting force and chip thickness, and improve the surface finish. Interestingly, removing the water coating enables the recovery of the cutting force, demonstrating a reversible effect. Based on the XAFS results and molecular dynamics simulation, a water-induced surface ordering mechanism is proposed to explain the findings from the microcutting experiments. This mechanism suggests that water molecules can induce surface ordering in copper, resulting in reduced surface energy and fracture toughness of copper, thus enhancing machinability. In conclusion, this work provides valuable insights into the comprehension of the Rehbinder effect and shows that picometer-scale modifications of the surface atom arrangement can considerably alter the deformation mode of metals, paving the way for the development of new manufacturing processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-efficiency purification of CH 4 and H 2 energy sources enabled by a phosphotungstic acid-supported Os single-atom catalyst

Methane (CH 4 ) and hydrogen (H 2 ) show promise as low-carbon energy sources, but their impurities, including H 2 and CO, pose challenges for storage and use. To address these challenges, a robust purification protocol for CH 4 and/or H 2 , combined with the catalytic conversion of impurities into CO 2 and H 2 O, is a compelling solution. Here, in this work, we investigated 11 phosphotungstic acid (PTA)-supported single-atom catalysts (SACs) by density functional theory (DFT) computations. Os 1 /PTA SACs exhibited superior catalytic activity, and the ease of oxidation follows the CO > H 2 > CH 4 order. It facilitated efficient purification of CH 4 in solvents such as water, MeOH, and various others. For H 2 purification, Os 1 /PTA SACs demonstrated excellent performance in gas, water, and MeOH. Notably, in water and MeOH, it selectively removed CO without consuming H 2 with low free energy barriers. The strong Os-PTA interactions and charge transfer mechanism contributed to its exceptional catalytic activity. Our findings shed light on SAC behavior and their potential for efficient CH 4 and H 2 purification. By addressing impurity challenges and improving clean energy utilization, these findings contribute to the development of sustainable energy technologies.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Low-energy magneto-optics of Tb 2 Ti 2 O 7 in a [111] magnetic field

The pyrochlore magnet Tb 2 Ti 2 O 7 shows a lack of magnetic order to low temperatures and is considered to be a quantum spin liquid candidate. Here we perform time-domain THz spectroscopy on high-quality Tb 2 Ti 2 O 7 crystals and study the low-energy excitations as a function of [111] magnetic field with high energy resolution. The low-energy crystal-field excitations change their energies anomalously under magnetic field. Despite several sharp field-dependent changes, we show that the material's spectrum can be described not by phase transitions but by field-dependent hybridization between the low-energy crystal-field levels. We highlight the strong coupling between spin and lattice degrees of freedom in Tb 2 Ti 2 O 7 as evidenced by the magnetic-field tunable crystal-field environment. Calculations based on single ion physics with field-induced symmetry reduction of the crystal-field environment can reproduce our data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nature of Excitons and Their Ligand-Mediated Delocalization in Nickel Dihalide Charge-Transfer Insulators

The fundamental optical excitations of correlated transition-metal compounds are typically identified with multielectronic transitions localized at the transition-metal site, such as d d transitions. In this vein, intense interest has surrounded the appearance of sharp, below-band-gap optical transitions, i.e., excitons, within the magnetic phase of correlated Ni 2 + van der Waals magnets. The interplay of magnetic and charge-transfer insulating ground states in Ni 2 + systems raises intriguing questions on the roles of long-range magnetic order and of metal-ligand charge transfer in the exciton nature, which inspired microscopic descriptions beyond typical d d excitations. Here we study the impact of charge transfer and magnetic order on the excitation spectrum of the nickel dihalides ( Ni X 2 , X = Cl , Br, and I) using Ni- L 3 edge resonant inelastic x-ray scattering (RIXS). In all compounds, we detect sharp excitations, analogous to the recently reported excitons, and assign them to spin-singlet multiplets of octahedrally coordinated Ni 2 + stabilized by intra-atomic Hund’s exchange. Additionally, we demonstrate that these excitons are dispersive using momentum-resolved RIXS. Our data evidence a ligand-mediated multiplet dispersion, which is tuned by the charge-transfer gap and independent of the presence of long-range magnetic order. This reveals the mechanisms governing nonlocal interactions of on-site d d excitations with the surrounding crystal or magnetic structure, in analogy to ground-state superexchange. These measurements thus establish the roles of magnetic order, self-doped ligand holes, and intersite-coupling mechanisms for the properties of d d excitations in charge-transfer insulators. Published by the American Physical Society 2024

2-dimensional systems↗
DOE OSTI2023

Atomic masses with machine learning for the astrophysical r process

The astrophysical r process plays a vital role in the production of heavy elements. Modeling of the r process is sensitive to masses and further requires knowledge of masses beyond current experimental reach. Therefore, simulations of the r process offer a unique test bed for predicting mass extrapolations. We take a Machine-Learning (ML) approach to model the masses across the entire chart of nuclides. For the first time, we simulate r-process nucleosynthesis with a physics-based ML mass model. We compare simulated abundances to solar data in order to evaluate the model's performance far from stability. The resulting r-process abundances up to thorium and uranium qualitatively match those of the observed solar system abundance pattern, with the characteristic peaks well positioned. We propagate the mass uncertainties obtained from the ML model to r-process abundance yields to estimate an uncertainty band associated with our approach. The size of the uncertainty band is approximately one order of magnitude which aligns with the uncertainty reported using alternative techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Theory of nuclear fission

Atomic nuclei are quantum many-body systems of protons and neutrons held together by strong nuclear forces. Under the proper conditions, nuclei can break into two (sometimes three) fragments which will subsequently decay by emitting particles. This phenomenon is called nuclear fission. Since different fission events may produce different fragmentations, the end-products of all fissions that occurred in a small chemical sample of matter comprise hundreds of different isotopes, including α particles, together with a large number of emitted neutrons, photons, electrons and antineutrinos. The extraordinary complexity of this process, which happens at length scales of the order of a femtometer, mostly takes less than a femtosecond but is not entirely over until all the lingering β decays have completed – which can take years – is a fascinating window into the physics of atomic nuclei. While fission may be more naturally known in the context of its technological applications, it also plays a crucial role in the synthesis of heavy elements in astrophysical environments. In both cases, simulations are needed for the many systems or energies inaccessible to experiments in the laboratory. In this context, the level of accuracy and precision required poses formidable challenges to nuclear theory. Overall, the goal of this article is to provide a comprehensive overview of the theoretical methods employed in the description of nuclear fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Helium impurities and interactions in lithium

Here we investigate helium interactions with lithium using density functional theory. Like other body-centered cubic (bcc) metals, the lowest-energy site for interstitial helium is a tetrahedral site. However, helium in lithium shows a higher substitutional formation energy than either the tetrahedral or octahedral interstitial formation energies, which is unique. The calculated migration energies of helium in lithium are also very low ( ~ 1 meV), an order of magnitude lower than those in other bcc metals. We find an increase in the binding energy to a vacancy as the number of helium atoms bound to it increases. The extremely low migration energies suggest that helium transport in lithium will be very fast.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a 1 mm radius vitreous carbon pellet (~3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cold gas mass measurements for the era of large optical spectroscopic surveys

ABSTRACT Gas plays an important role in many processes in galaxy formation and evolution, but quantifying the importance of gas has been hindered by the challenge to measure gas masses for large samples of galaxies. Data sets of direct atomic and molecular gas measurements are sufficient to establish simple scaling relations, but often not large enough to quantify three-parameter relations, or second-order dependences. As an alternative approach, we derive here indirect cold gas measurements from optical emission lines using photoionization models for galaxies in the Sloan Digital Sky Survey (SDSS) main galaxy sample and the PHANGS-MUSE survey. We calibrate the gas surface density measurements using xCOLD GASS and PHANGS-ALMA molecular gas measurements to ensure that our measurements are reliable. We demonstrate the importance of taking into account the scale dependence of the relation between optical depth (τ V ) and gas surface density (Σ gas ) and provide a general prescription to estimate Σ gas from τ V , metallicity, and the dust-to-metal ratio, at any arbitrary physical resolution. To demonstrate that the indirect cold gas masses are accurate enough to quantify the role of gas in galaxy evolution, we study the mass–metallicity relation of SDSS galaxies and show that as a third parameter gas mass is better than star formation rate at reducing the scatter of the relation, as predicted by models and simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

First-principles study of bilayer polymeric manganese phthalocyanine

Here we study bilayer manganese phthalocyanine (MnPc) molecules and MnPc polymeric sheets using first-principles simulations with a focus on the magnetic interactions between Mn atoms. We find that the most stable position of the upper layer with respect to the lower layer is shifted about 1/8 of a lattice vector from the center of the bottom layer along the direction toward a nearest-neighbor N atom. The magnetic ground state is the Néel antiferromagnetic (AFM) configuration within a layer and ferromagnetic (FM) between Mn atoms in adjacent layers. In this state, the system becomes a semiconductor with an indirect band gap of 11meV. The strongest interaction is the interlayer coupling between the closest Mn atoms. A maximally localized Wannier analysis suggests that the dominant coupling pathway is Mn-(N,C)-Mn rather than a direct Mn-Mn coupling. The maximum calculated magnetic anisotropy energy is found to be 1.0meV per Mn atom. We also find that the bilayer molecule shows a significant stacking angle change from FM to AFM configurations accompanied by a change of orbital filling ordering.

2-dimensional systems↗
DOE OSTI2023

Strain-tuned quantum criticality in electronic Potts-nematic systems

Motivated by recent observations of threefold rotational symmetry breaking in twisted moiré systems, cold-atom optical lattices, quantum Hall systems, and triangular antiferromagnets, we phenomenologically investigate the strain-temperature phase diagram of the electronic 3-state Potts-nematic order. While in the absence of strain the quantum Potts-nematic transition is first-order, quantum critical points (QCP) emerge when uniaxial strain is applied, whose nature depends on whether the strain is compressive or tensile. In one case, the nematic amplitude jumps between two non-zero values while the nematic director remains pinned, leading to a symmetry-preserving metanematic transition that terminates at a quantum critical end-point. For the other type of strain, the nematic director unlocks from the strain direction and spontaneously breaks an in-plane twofold rotational symmetry, which in twisted moiré superlattices triggers an electric polarization. Such a piezoelectric transition changes from first to second-order upon increasing strain, resulting in a quantum tricritical point. Using a Hertz-Millis approach, we show that these QCPs share interesting similarities with the widely studied Ising-nematic QCP. The existence of three minima in the nematic action also leaves fingerprints in the strain-nematic hysteresis curves, which display multiple loops. At non-zero temperatures, because the upper critical dimension of the 3-state Potts model is smaller than three, the Potts-nematic transition is expected to remain first-order in 3D, but to change to second-order in 2D. As a result, the 2D strain-temperature phase diagram displays two first-order transition wings bounded by lines of critical end-points or tricritical points, reminiscent of the phase diagram of metallic ferromagnets. Furthermore, we discuss how our results can be used to unambiguously identify spontaneous Potts-nematic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. Here, these methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantifying Stern layer water alignment before and during the oxygen evolution reaction

While water’s oxygen is the electron source in the industrially important oxygen evolution reaction, the strong absorber problem clouds our view of how the Stern layer water molecules orient themselves in response to applied potentials. Here, we report nonlinear optical measurements on nickel electrodes held at pH 13 indicating a disorder-to-order transition in the Stern layer water molecules before the onset of Faradaic current. A full water monolayer (1.1 × 10 15 centimeter −2 ) aligns with oxygen atoms pointing toward the electrode at +0.8 volt and the associated work is 80 kilojoule per mole. Our experiments identify water flipping energetics as a target for understanding overpotentials, advance molecular electrochemistry, provide benchmarks for electrical double layer models, and serve as a diagnostic tool for understanding electrocatalysis.

Science & Technology - Other Topics↗