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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Kinetics of short-range order formation in GeSn alloy: MBE vs CVD

Recently, short-range order (SRO) has attracted significant attention, challenging the conventional view of the atomic positions in alloys being random. Furthermore, the presence of SRO has been predicted to have profound effects on the electronic and topological properties of group-IV alloys, offering a different direction in designing group-IV materials for photoelectronic and quantum devices. However, due to the limited understanding of the formation mechanisms, developing effective methods to manipulate SRO in epitaxy is still challenging. To address this, we propose a mechanism for the GeSn alloy, revealing that surface diffusion plays a key role in SRO formation. Building on this mechanism, we show that the distinct surface conditions in MBE and CVD lead to the formation of SRO with enhanced Sn–Sn pairing in MBE-grown samples, while CVD-grown samples remain random alloys. Furthermore, our findings provide an initial understanding of the kinetic process of SRO formation, providing guidance for the design of experiments to manipulate SRO.

Alloys↗

DOE OSTI2023

World beyond the nearest neighbors

The structure beyond the nearest neighbor atoms in liquid and glass is characterized by the medium-range order (MRO). In the conventional approach, the MRO is considered to result directly from the short-range order (SRO) in the nearest neighbors. To this bottom–up approach starting with the SRO, we propose to add a top–down approach in which global collective forces drive liquid to form density waves. The two approaches are in conflict with each other, and the compromise produces the structure with the MRO. The driving force to produce density waves provides the stability and stiffness to the MRO, and controls various mechanical properties. This dual framework provides a novel perspective for description of the structure and dynamics of liquid and glass.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Precision measurement of lifetime and branching ratios of the 4⁢𝑓 13⁡ 5⁢𝑑⁢6⁢𝑠 1 ⁢[5/2] 5/2 state in Yb + ions

We report spectroscopic and time-resolved experimental observations to characterize the [Xe]⁢4⁢𝑓 13 ⁡( 2 𝐹$^𝑜_{5/2}$)⁢5⁢𝑑⁢6⁢𝑠⁡( 1 𝐷)⁢ 1 ⁢[5/2]$^𝑜_{5/2}$ state in 172 Yb + ions. We access this state from the metastable 4⁢𝑓 14 ⁡5⁢𝑑⁡( 2 𝐷 3/2,5/2 ) manifold and observe an unexpectedly long lifetime of 𝜏 = 37.9⁢(9) ⁢µ⁢s that allows visible Rabi oscillations and resolved-sideband spectroscopy. Using a combination of coherent population dynamics, high-fidelity detection and heralded state preparation, and optical pumping methods, we measure the branching ratios to the 2 𝐷 3/2 , 2 𝐷 5/2 , and 2 𝑆 1/2 states to be 0.359⁢(2), 0.639(2), and 0.0023⁢(16), respectively. The branching ratio to the 4⁢𝑓 13⁡ 6⁢𝑠 2 ⁡( 2 𝐹 7/2 ) is compatible with zero within our experimental resolution. We also report measurements of Landé 𝑔-factor of the 1 ⁢[5/2]$^𝑜_{5/2}$ state. Further, the branching ratio of the 2 𝐷 5/2 to 2 𝑆 1/2 decay in 172 Yb + was measured to be 0.188(3), improving its relative uncertainty by an order of magnitude. Our measurements provide experimental benchmarks for better understanding the atomic structure of Yb + ions, which still lacks accurate numerical descriptions, and the use of high-lying excited states for partial detection and qubit manipulation in the omg architecture.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

DOE OSTI2020

Modeling Catalyzed Growth of Single-Wall Carbon Nanotubes (Final Report)

Single-walled carbon nanotubes (SWCNTs) are among the most promising nanostructures because of their special and unusual physico-chemical properties that make them strong candidates for a variety of applications. Since their discovery, chemical vapor deposition (CVD) methods have become the preferred synthesis method. CVD uses a carbon-containing precursor species, which is decomposed over a catalyst yielding C atoms that eventually evolve in carbon nanostructures. CVD temperatures are in the order of 1000K or lower, and pressures in the range of 1-5 atm, although some processes use higher pressures. Depending on the nature of the catalyst, and on the conditions of pressure, temperature, and presence of other inert compounds, the resultant nanostructures may include nanofibers, nano-onions, multiple-, or single-wall carbon nanotubes, among others. When the CVD synthesis process is appropriately tuned, however, it is possible to obtain only SWCNTs, which are geometrically characterized by their diameter and by their chiral angle. The chiral indexes (n,m) determine the degree of helical twist of the graphite lattice along the nanotube axis, where the term chirality is derived from the left and right-handed helicity of the tubes defined by their chiral indexes (n,m) and (m,n) respectively. The main goal of this research is to build fundamental understanding that will help elucidating the nature and structure of the actual catalyst in a combined catalyst/support system, which we claim is essential especially when the catalyst is of nanodimensions. Given that we are considerably advanced in the understanding of selective single-walled carbon nanotube catalyst synthesis, we take this process as a case study. Thus, we focus on the understanding of the selective growth of SWCNTs on supported nanocatalysts in systems where the catalyst synthesis takes place in the same reactor preceding the nanotube synthesis. Once we understand what type of catalyst structure would be necessary for this specific growth, we will try to elucidate the catalytic synthesis process that leads to such catalytic structure. Thus, our aim of controlling nanostructure synthesis may start to be delineated.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗

NASA NTRS1973

Molecular collisions 21: Semiclassical approximation to atom-symmetric top rotational excitation

A distorted wave approximation to the T matrix for atom-symmetric top scattering was developed. The approximation is correct to first order in the part of the interaction potential responsible for transitions in the component of rotational angular momentum along the symmetry axis of the top. A semiclassical expression for this T matrix is derived by assuming large values of orbital and rotational angular momentum quantum numbers.

Russell, D.↗

NASA NTRS1975

Stability of localized plasma model in two and three dimensions

Stability against catastrophic radial collapse is demonstrated for initially well-behaved, localized solutions in a class of two- and three-dimensional wave equations characterized by saturating nonlinearities, which model various electrostatic and electromagnetic 'caviton' structures in an asymptotically uniform plasma. This result contrasts sharply with mathematical predictions of equations with low-order nonlinearities. Discrete classes of possible quasistationary modes exhibit features comparable to atomic wave functions and differ qualitatively from one-dimensional structures.

Wilcox, J. Z.↗

NASA NTRS1976

Theoretical L-shell Coster-Kronig energies 11 or equal to z or equal to 103

Relativistic relaxed-orbital calculations of L-shell Coster-Kronig transition energies have been performed for all possible transitions in atoms with atomic numbers. Hartree-Fock-Slater wave functions served as zeroth-order eigenfunctions to compute the expectation of the total Hamiltonian. A first-order approximation to the local approximation was thus included. Quantum-electrodynamic corrections were made. Each transition energy was computed as the difference between results of separate self-consistent-field calculations for the initial, singly ionized state and the final two-hole state. The following quantities are listed: total transition energy, 'electric' (Dirac-Hartree-Fock-Slater) contribution, magnetic and retardation contributions, and contributions due to vacuum polarization and self energy.

Chen, M. H.↗

NASA NTRS1978

A model of carbon production in a cometary coma

A model of the cometary ionosphere is developed in order to account for the large population of metastable C(D-1) atoms detected via the ultraviolet spectrum of Comet West (1976 VI). Dissociative recombination of CO(plus) ions and electrons is shown to be the dominant source of carbon atoms rather than photodissociation of CO so that the derived carbon production rate is only a lower limit to the evaporation rate of the carbon bearing mother molecule.

Feldman, P. D.↗

NASA NTRS1984

The calculation of molecular Eigen-frequencies

A method of determining molecular eigen-frequencies based on the function of Einstein expressing the variation of the atomic heat of various elements is proposed. It is shown that the same equation can be utilized to calculate both atomic heat and optically identifiably eigen-frequencies - at least to an order of magnitude - suggesting that in both cases the same oscillating structure is responsible.

Lindemann, F. A.↗

NASA NTRS1987

Stability relationship for water droplet crystallization with the NASA Lewis icing spray

In order to produce small droplets for icing cloud simulation, high pressure air atomizing nozzles are used. For certain icing testing applications, median drop sizes as small as 5 mm are needed, which require air atomizing pressures greater than 3000 kPa. Isentropic expansion of the ambient temperature atomizing air to atmospheric pressure can result in air stream temperatures of -160 C which results in ice crystals forming in the cloud. To avoid such low temperatures, it is necessary to heat the air and water to high initial temperatures. An icing spray research program was conducted to map the temperatures below which ice crystals form. A soot slide technique was used to determine the presence of crystals in the spray.

Marek, C. John↗

NASA NTRS1988

Stability relationship for water droplet crystallization with the NASA Lewis icing spray nozzle

Marek, C. John↗

NASA NTRS1988

Elemental abundance analyses with coadded DAO spectrograms. IV - Revision of previous analyses. V - The mercury-manganese stars Phi Herculis, 28 Herculis and HR 7664

Changes in chromium, manganese, and nickel abundances derived from singly ionized lines are incorporated into the elemental abundance of Adelman and Hill (1987) in order to provide more accurate gf values and damping constants for several atomic species. An improved agreement with the values from neutral lines of the same element is found. In the second part, the method is applied to an elemental abundance analysis of three mercury-manganese stars, and correlations are found between the derived abundances and the effective temperature.

Adelman, Saul J.↗

NASA NTRS1991

Detection of global state predicates

The problem addressed here arises in the context of Meta: how can a set of processes monitor the state of a distributed application in a consistent manner? For example, consider the simple distributed application as shown here. Each of the three processes in the application has a light, and the control processes would each like to take an action when some specified subset of the lights are on. The application processes are instrumented with stubs that determine when the process turns its lights on or off. This information is disseminated to the control processes, each of which then determines when its condition of interest is met. Meta is built on top of the ISIS toolkit, and so we first built the sensor dissemination mechanism using atomic broadcast. Atomic broadcast guarantees that all recipients receive the messages in the same order and that this order is consistent with causality. Unfortunately, the control processes are somewhat limited in what they can deduce when they find that their condition of interest holds.

Marzullo, Keith↗

NASA NTRS1992

The molecular interstellar medium in two giant H II regions in M33 - NGC 604 and NGC 595

This study reports single-dish and interferometric CO observations for the two largest H II regions in M33, NGC 604, and NGC 595. Six individual molecular clouds are identified in the H II regions. Molecular and virial masses for three resolved clouds in NGC 604 agree very well, from which it is inferred that the CO-to-H2 conversion factor in this region is the same as it is in the inner disk of M33 and in the Galaxy, despite a metallicity in NGC 604 that is three times below solar. The total gas mass of the molecular complex associated with NGC 604 estimated from single-dish data is dissimilar to those of molecular complexes in the inner disk of M33, while the molecular gas mass of NGC 505 is roughly an order or magnitude smaller. Both regions contain a large mass fraction of atomic hydrogen, which suggests that photodissociation of molecular hydrogen is an important process in these energetic H II regions.

Wilson, Christine D.↗

NASA NTRS1998

Analysis of Phase Separation in Czochralski Grown Single Crystal Ilmenite

Ilmenite (FeTiOs) is a wide bandgap semiconductor with an energy gap of 2.58 eV. Ilmenite has properties suited for radiation tolerant applications, as well as a variety of other electronic applications. Single crystal ilmenite has been grown from the melt using the Czochralski method. Growth conditions have a profound effect on the microstructure of the samples. Here we present data from a variety of analytical techniques which indicate that some grown crystals exhibit distinct phase separation during growth. This phase separation is apparent for both post-growth annealed and unannealed samples. Under optical microscopy, there appear two distinct areas forming a matrix with an array of dots on order of 5 pm diameter. While appearing bright in the optical micrograph, atomic force microscope (AFM) shows the dots to be shallow pits on the surface. Magnetic force microscope (MFM) shows the dots to be magnetic. Phase identification via electron microprobe analysis (EMPA) indicates two major phases in the unannealed samples and four in the annealed samples, where the dots appear to be almost pure iron. This is consistent with micrographs taken with a scanning probe microscope used in the magnetic force mode. Samples that do not exhibit the phase separation have little or no discernible magnetic structure detectable by the MFM.

Wilkins, R.↗