Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulations”

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 127 records · Page 7

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗

Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction

Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.

BISON↗

ChemGraph as an agentic framework for computational chemistry workflows

Atomistic simulations are essential in chemistry and materials science but remain challenging to run due to the expert knowledge required for the setup, execution, and validation stages of these calculations. We present ChemGraph, an agentic framework powered by artificial intelligence and state-of-the-art simulation tools to streamline and automate computational chemistry and materials science workflows. ChemGraph leverages graph neural network-based foundation models for accurate yet computationally efficient calculations and large language models (LLMs) for natural language understanding, task planning, and scientific reasoning to provide an intuitive and interactive interface. We evaluate ChemGraph across 13 benchmark tasks and demonstrate that smaller LLMs (GPT-4o-mini, Claude-3.5-haiku, Qwen-2.5-14B) perform well on simple workflows, while more complex tasks benefit from using larger models. Importantly, we show that decomposing complex tasks into smaller subtasks through a multi-agent framework enables GPT-4o to reach perfect accuracy and smaller LLMs to match or exceed single-agent GPT-4o's performance in these benchmarks.

Computational chemistry↗

Stress-dependent χ phase transformation in a Ni-based superalloy

The ongoing push to elevate operating temperatures in aerospace gas turbine engines – driven by goals of enhanced fuel efficiency and reduced CO 2 emissions – mandates advancements in the creep resistance of Ni- and Co-based superalloys, which are integral for critical engine components. This study elucidates the role of stress assisted localized phase transformations in the creep properties of these alloys. By leveraging chemo-mechanical coupling, self-healing γ′ precipitates are designed to immobilize planar defects, thereby increasing creep resistance. Employing advanced characterization techniques such as high-resolution Scanning Transmission Electron Microscopy (HR-STEM), in conjunction with atomistic simulations and thermodynamic calculations, novel deformation pathways facilitated by χ local phase transformation (LPT) strengthening have been uncovered; notably, the formation of χ nano-laths through microtwinning and superlattice intrinsic stacking fault (SISF) shearing. This study highlights critical insights into the compositional boundaries necessary for optimizing LPT strengthening while avoiding deleterious bulk formation of η/χ phases. These advancements will guide the design of new alloys maximizing high-temperature creep strength for advanced aerospace applications.

Egan, Ashton J. [Friedrich-Alexander University Er↗

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗

Structure and mechanical properties of grain boundaries in molybdenum disulfide (MoS 2 )

Molybdenum disulfide (MoS 2 ) is a two-dimensional material widely used as a lubricant in many applications involving mechanical loading under a wide range of operating temperatures. Since many synthesis and processing techniques yield MoS 2 in its polycrystalline form, establishing grain boundary (GB) structure–property relations is key to designing MoS 2 microstructures with tailored properties. Here, we employ classical, reactive atomistic simulations to investigate the structure and mechanical behavior of a wide range of GBs in MoS 2 as a function of temperature. Using a bicrystal MoS 2 geometry, we characterize the atomic structure and calculate the energy of several low-angle GBs. Then, we simulate the tensile deformation behavior of MoS 2 bicrystals at several temperatures. Our results reveal that at temperatures above 100 K, the deformation of MoS 2 bicrystals is characterized by the nucleation of shear bands from GBs that grow, with subsequent loading, into the MoS 2 crystals. At low temperatures, the tensile deformation is characterized by the nucleation and propagation of deformation fronts, resulting in altered bond angles and bond lengths. Quantitative analysis reveals a decrease in the ultimate tensile stress and ultimate failure strain of MoS 2 bicrystals with the increase in temperature. Furthermore, our simulations of the mechanical behavior of metastable GBs reveal that the strength and ductility decrease with the increase in energy of these boundary structures. In broad terms, our work provides future avenues to employ GB engineering as a strategy to tailor the properties of MoS 2 microstructures.

Moore, Robert D. [Lehigh Univ., Bethlehem, PA (Uni↗

An idealized model for the crystal structure of intermetallic compounds isostructural with Mg 3 Cr 2 Al 18

Mg 3 Cr 2 Al 18 (abbreviated in this report as MCA) is the parent phase for a large class of intermetallic compounds that belong to the cubic crystal space group, $Fd\overline{3}m$. The purpose of this paper is to introduce an ideal, unrelaxed crystal structure for compounds isostructural with MCA. There are five distinct atomic sublattices in MCA compounds, which can be denoted, $A, B, C, D,$ and $E$. With this, a general description for MCA structures can be written as $A^{8a}_{1}B^{16c}_{2}C^{16d}_{2}D^{48f}_{6}E^{96g}_{12}$, where the superscripts represent the Wyckoff special equipoints associated with the various sublattices in MCA, and the subscripts indicate the contributions of each sublattice to the stoichiometry of one formula unit in an any given MCA structured compound. Sublattices D and E are where deviations from ideality occur in real, MCA-like compounds. This paper examines MCA bond lengths, nearest-neighbour polyhedral arrangements, 3-D sublattice crystal structures, 2-D atom tessellation patterns, and crystal chemical effects associated with atomic relaxations on the $D$ and $E$ sublattices. The ideal MCA crystal structure developed in this report provides an appropriate initial structure for use as input to crystal structure refinements of diffraction data for MCA-like phases being examined experimentally, or as input for computational, atomistic simulations of the structures of such compounds.

36 MATERIALS SCIENCE↗

Reversible to irreversible transitions for ac driven skyrmions on periodic substrates

Abstract Using atomistic simulations, we investigate the dynamical behavior of magnetic skyrmions in dimer and trimer molecular crystal arrangements, as well as bipartite lattices at 3/2 and 5/2 fillings, under ac driving over a square array of anisotropy defects. For low ac amplitudes, at all fillings reversible motion appears in which the skyrmions return to their original positions at the end of each ac drive cycle and the diffusion is zero. We also identify two distinct irreversible regimes. The first is a translating regime in which the skyrmions form channels of flow in opposing directions and translate by one substrate lattice constant per ac drive cycle. The translating state appears in the dimer and trimer arrangements, and produces pronounced peaks in the diffusivity in the direction perpendicular to the external drive. For larger ac amplitudes, we find chaotic irreversible motion in which the skyrmions can randomly exchange places with each other over time, producing long-time diffusive behavior both parallel and perpendicular to the ac driving direction.

36 MATERIALS SCIENCE↗

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Temperature-dependent spin dynamics in Cr 2 O 3

Cr 2 ⁢O 3 emerges as a prominent candidate material for spintronics and magnetoelectronic applications. However, a comprehensive understanding of its temperature-dependent spin dynamics remains elusive, impeding the engineering of novel spintronics based on this material. Here, we delve into this through a combination of inelastic neutron scattering experiments and atomistic simulations. Our results unveiled the emergence of paramagnons above and below 𝑇 𝑁 . We demonstrated a significant softening of linear magnons upon heating in the antiferromagnetic state. Further analysis revealed that this softening primarily originated from four-magnon interactions, while thermal expansion played a minor role.

antiferromagnetism↗

Surface acoustic waves and lattice vibrations in two-dimensional Ti 3 ⁢C 2 ⁢𝑇 𝑥 (𝑇 =O, F) MXene films

We investigated surface acoustic wave (SAW) propagation and lattice vibrations in two-dimensional (2D) titanium carbide (Ti 3 ⁢C 2 ⁢𝑇 x ) MXene films as a function of surface termination and layer stacking, using atomistic simulations. We found that SAW propagation velocity is highly sensitive to both single-layer properties and interlayer bonding. Surface terminations significantly modulate wave behavior, with oxygen and fluorine terminations producing distinct effects on wave propagation, with oxygen-terminated monolayers exhibiting 20% higher wave speeds than fluorine counterparts due to strengthened intralayer bonds. Key observations include the transition from one to two layers causing wave speed variations, and the development of interlayer modes that generate more dispersed lattice vibrations. As the film layer thickness increases, SAW propagation becomes predominantly confined to the upper surface, with coherence of vibrational modes diminishing in multilayer structures. These findings suggest MXene terminations and layer stacking are crucial parameters for controlling SAW behavior, offering promising avenues for novel acoustic wave device applications.

2-dimensional systems↗

Stability and Mobility of Disconnections in Solute Atmospheres: Insights from Interfacial Defect Diagrams

This Letter explores the stability of disconnections (step-dislocation defects) at grain boundaries in binary alloys. We introduce interfacial defect diagrams, derived from atomistic simulations and segregation theory, to predict the stability of disconnections in the temperature-solute concentration phase space and relate it to governing segregation mechanisms. These diagrams reveal multiple stability regimes influenced by solute-induced clustering and pinning effects impacting the thermal migration of disconnections and offering insights into their thermodynamics and kinetic properties.

Crystal structure↗

Impurity induced confinement effects in size-separated Mn-doped CsPbCl 3 nanocrystals

Beyond quantum confinement, impurity doping of semiconductor nanocrystals (NCs) offers a degree of control over their optical and electronic properties, which are of interest for potential applications. Manganese (Mn)-doped CsPbCl 3 NCs, for example, show synergistic effects in the visible, with large enhancements in both blue (exciton) and red (dopant) emission over a narrow range of dopant concentration. Although Mn concentration is the primary parameter, NC size offers an additional degree of control over the number of dopants per NC. Size-resolved CsPbCl 3 NCs prepared in the limit of strong Mn doping highlight the emergence of defective NCs with deviations from anticipated confinement trends. In conclusion, the transition correlates with a redshift and increase in the amplitude of dopant emission, and we discuss these observations in the context of atomistic simulations of Ruddlesden-Popper defects in CsPbCl 3 .

36 MATERIALS SCIENCE↗

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗

Skyrmion-skyrmionium phase separation and laning transitions via spin-orbit torque currents

Many driven binary systems can exhibit laning transitions when the two species have different mobilities, such as colloidal particles with opposite charges in electric fields. Another example is pedestrian or active-matter systems, where particles moving in opposite directions form a phase-separated state that enhances the overall mobility. In this work, we use atomistic simulations to demonstrate that mixtures of skyrmions and skyrmioniums also exhibit pattern formation and laning transitions. Skyrmions move more slowly and at a finite skyrmion Hall angle compared to skyrmioniums, which move faster and without a skyrmion Hall effect. At low drives, the system forms a partially jammed phase where the skyrmioniums are dragged by the surrounding skyrmions, resulting in a finite angle of motion for the skyrmioniums. At higher drives, the system transitions into a laned state, but unlike colloidal systems, the lanes in the skyrmion-skyrmionium mixture are tilted relative to the driving direction due to the intrinsic skyrmion Hall angle. In the laned state, the skyrmionium angle of motion is reversed when it aligns with the tilted lane structure. At even higher drives, the skyrmioniums collapse into skyrmions. Below a critical skyrmion density, both textures can move independently with few collisions, but above this density, the laning state disappears entirely, and the system transitions to a skyrmion-only state. We map out the velocity and Hall responses of the different textures and identify three distinct phases: partially jammed, laned, and skyrmion-only moving crystal states. In conclusion, we compare our results to recent observations of tilted laning phases in pedestrian flows, where chiral symmetry breaking in the particle interactions leads to similar behavior.

36 MATERIALS SCIENCE↗

Grain boundary energy control in zinc aluminate nanoceramics

This study investigates the grain boundary energy dependence on segregated dopants in nanocrystalline zinc aluminate ceramics. Atomistic simulations of Σ3 and Σ9 grain boundaries showed that trivalent ions of varying ionic radii [Sc 3+ (74.5 pm), In 3+ (80.0 pm), Y 3+ (90.0 pm), and Nd 3+ (98.3 pm)] have a tendency to segregate to both interfaces, with Y 3+ presenting the highest segregation potentials. The connection between segregation and the reduction of interfacial energies was explored by measuring the grain boundary energy on nanoceramics fabricated via high‐pressure spark plasma sintering (HP‐SPS) using differential scanning calorimetry (DSC). The results revealed that Y 3+ doping at 0.5 mol% reduces the grain boundary energy in zinc aluminate nanoceramics from 1.1–1.3 J/m 2 to 0.6–0.8 J/m 2 ; the range correlates with the observed size dependence of the excess energy, with higher values observed for the smaller grain sizes (∼17 nm). The noted decrease in interfacial energies for doped samples suggests it is indeed possible to alter the stability of zinc aluminate grain boundaries via dopant segregation.

36 MATERIALS SCIENCE↗

Ordered and disordered skyrmion states on a square substrate

We examine the ordering of skyrmions interacting with a square substrate created from a modulation of anisotropy using atomistic simulations. We consider fillings of ƒ = 1.5, 2.0, 2.5, 3.0, and 3.5 skyrmions per potential minimum as a function of magnetic field and sample size. For a filling of ƒ = 2.0, we find various dimer orderings, such as tilted dimer states, as well as antiferromagnetic ordering that is similar to the colloidal dimer ordering seen on square substrates. Furthermore, the ability of the skyrmions to change shape or annihilate produces additional states that do not occur in the colloidal systems. For certain parameters at ƒ = 2.0, half of the skyrmions can annihilate to form a square lattice, or a superlattice of trimers and monomers containing skyrmions of different sizes can form. At lower fields, ordered stretched skyrmion states can appear, and for zero field, there can be ordered stripe states. For ƒ = 3.0, we find ferromagnetic ordered trimers, tilted lattices, columnar lattices, and stretched phases. For fillings of ƒ = 1.5 and 2.5, we find bipartite lattices, different ordered and disordered states, and several extended disordered regions produced by frustration effects.

36 MATERIALS SCIENCE↗