Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulation”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2020

Observation of Magnon Polarization

We measure the mode-resolved direction of the precessional motion of the magnetic order, i.e.,magnon polarization, via the chiral term of inelastic polarized neutron scattering spectra. The magnon polarization is a unique and unambiguous signature of magnets and is important in spintronics, affecting thermodynamic properties such as the magnitude and sign of the spin Seebeck effect. However, it has never been directly measured in any material until this work. The observation of both signs of magnon polarization in Y 3 Fe 5 O 12 also gives direct proof of its ferrimagnetic nature. The experiments agree very well with atomistic simulations of the scattering cross section.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ultrafast Neuromorphic Dynamics Using Hidden Phases in the Prototype of Relaxor Ferroelectrics

Materials possessing multiple states are promising to emulate synaptic and neuronic behaviors. Their operation frequency, typically in or below the GHz range, however, limits the speed of neuromorphic computing. Ultrafast THz electric field excitation has been employed to induce nonequilibrium states of matter, called hidden phases in oxides. One may wonder if there are systems for which THz pulses can generate neuronic and synaptic behavior, via the creation of hidden phases. Using atomistic simulations, we discover that relaxor ferroelectrics can emulate all the key neuronic and memristive synaptic features. Their occurrence originates from the activation of many hidden phases of polarization order, resulting from the response of nanoregions to THz pulses. We find such phases further possess different dielectric constants, which is also promising for memcapacitor devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

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↗
DOE OSTI2022

Predicting the location of shear band initiation in a metallic glass

We report atomistic simulation results which indicate that the location of shear banding in a metallic glass (MG) can be ascertained with reasonably high accuracy solely from the undeformed static structure. Correlation is observed between the location of the initiation of shear bands in a simulated MG and the initial distribution of the density of fertile sites (DFS) for stress-driven shear transformations identified a priori based on a deep learning model devised in our recent work. In addition, we demonstrated that one can judge whether a glass is brittle or ductile solely based upon its initial DFS distribution. These validate that shear bands in MG arise from nonlinear instabilities and that the as-quenched glass structure contains inhomogeneities that influence these instabilities. In this study, we also reveal an important subtlety regarding the nondeterministic nature of athermal quasistatic shear simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geometrically necessary dislocation fingerprints of dislocation loop absorption at grain boundaries

Here we present a numerical methodology to compute the Nye-tensor fingerprints of dislocation loop absorption at grain boundaries (GBs) for comparison with TEM observations of irradiated polycrystals. Our approach links atomistic simulations of self-interstitial atom (SIA) prismatic loops gliding toward and interacting with GBs in body-centered cubic iron with experimentally extracted geometrically necessary dislocation (GND) maps to facilitate the interpretation of damage processes. The Nye-tensor analysis is strongly mesh-size dependent—corresponding to resolution-dependent TEM observations. The method computes GND fingerprints from discretized dislocation line segments extracted from molecular dynamics simulations of dislocation loops being absorbed at a GB. Specifically, we perform MD simulation of prismatic loops of two diameters and monitor the three stages of the absorption process: loop glide, the partial, and full absorption of the loops at a [1 0 0] symmetric tilt GB. These methods provide a framework for future investigations of the nature of defect absorption by grain boundaries under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2022

Identifying the microstructural features associated with void nucleation during elevated–temperature deformation of copper

The microstructural-scale mechanisms that produce cracks in metals during deformation at elevated temperatures are relevant to applications that involve thermal exposure. Prior studies of cavitation during high-temperature deformation, for example, creep, suffered from an inability to directly observe the microstructural evolution that occurs during deformation and leads to void nucleation. Here the current study takes advantage of modern high-speed electron backscatter diffraction (EBSD) detectors to observe cavitation in oxygen-free, high-conductivity copper in situ during deformation at 300°C. Most voids formed at the triple junction between a twin boundary and a high-angle grain boundary (HAGB). This finding does not contradict previous studies that suggested that twins are resistant to cracking—it reveals that cracks in HAGBs originate at twin/HAGB triple junctions and that cracks preferentially grow along HAGBs rather than the accompanying twins. Atomistic simulations explored the origins of this observation and suggest that twin/HAGB triple junctions are microstructural weak points.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inversion, chemical complexity, and interstitial transport in spinels

Abstract Spinels with the generic chemical formula AB 2 O 4 have potential applications in nuclear energy and batteries. In both cases, their functionality is related to mass transport through the crystal. Here, using long‐time atomistic simulations, we examine the impact of the cation structure on interstitial transport in two spinel chemistries, inverse MgGa 2 O 4 and double MgAlGaO 4 . We emphasize two aspects of the transport properties: the unit mechanisms that are described by individual barriers, for which we introduce pole‐figure‐like plots, and the aggregate behavior of those unit mechanisms. Compared to previous work on normal spinels, we find that inversion significantly reduces the rate of interstitial transport in these structures and has an impact on the stability of defects as they move through the lattice. In particular, B cation interstitials are found to be kinetically stable only in the inverse MgGa 2 O 4 . These results provide new insight into relationship between structure, chemistry, and transport in spinels.

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

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↗
DOE OSTI2022

Irradiation-induced grain boundary facet motion: In situ observations and atomic-scale mechanisms

Metals subjected to irradiation environments undergo microstructural evolution and concomitant degradation, yet the nanoscale mechanisms for such evolution remain elusive. Here, we combine in situ heavy ion irradiation, atomic resolution microscopy, and atomistic simulation to elucidate how radiation damage and interfacial defects interplay to control grain boundary (GB) motion. While classical notions of boundary evolution under irradiation rest on simple ideas of curvature-driven motion, the reality is far more complex. Focusing on an ion-irradiated Pt Σ3 GB, we show how this boundary evolves by the motion of 120° facet junctions separating nanoscale {112} facets. Our analysis considers the short- and mid-range ion interactions, which roughen the facets and induce local motion, and longer-range interactions associated with interfacial disconnections, which accommodate the intergranular misorientation. We suggest how climb of these disconnections could drive coordinated facet junction motion. These findings emphasize that both local and longer-range, collective interactions are important to understanding irradiation-induced interfacial evolution.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Maximum strength and dislocation patterning in multi–principal element alloys

Multi–principal element alloys (MPEAs) containing three or more components in high concentrations render a tunable chemical short-range order (SRO). Leveraging large-scale atomistic simulations, we probe the limit of Hall-Petch strengthening and deformation mechanisms in a model CrCoNi alloy and unravel chemical ordering effects. The presence of SRO appreciably increases the maximum strength and lowers the propensity for faulting and structure transformation, accompanied by intensification of planar slip and strain localization. Deformation grains exhibit notably different microstructures and dislocation patterns that prominently depend on their crystallographic orientation and the number of active slip planes. Grain of single-planar slip attains the highest volume fraction of deformation-induced structure transformation, and grain with double-slip planes develops the densest dislocation network. These results advancing the fundamental understanding of deformation mechanisms and dislocation patterning in MPEAs suggest a mechanistic strategy for tuning mechanical behavior through simultaneously tailoring grain texture and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transonic dislocation propagation in diamond

Here, the motion of line defects (dislocations) has been studied for more than 60 years, but the maximum speed at which they can move is unresolved. Recent models and atomistic simulations predict the existence of a limiting velocity of dislocation motion between the transonic and subsonic ranges at which the self-energy of dislocation diverges, though they do not deny the possibility of the transonic dislocations. We used femtosecond x-ray radiography to track ultrafast dislocation motion in shock-compressed single-crystal diamond. By visualizing stacking faults extending faster than the slowest sound wave speed of diamond, we show the evidence of partial dislocations at their leading edge moving transonically. Understanding the upper limit of dislocation mobility in crystals is essential to accurately model, predict, and control the mechanical properties of materials under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dissolution enables dolomite crystal growth near ambient conditions

Crystals grow in supersaturated solutions. A mysterious counterexample is dolomite CaMg(CO 3 ) 2 , a geologically abundant sedimentary mineral that does not readily grow at ambient conditions, not even under highly supersaturated solutions. Using atomistic simulations, we show that dolomite initially precipitates a cation-disordered surface, where high surface strains inhibit further crystal growth. However, mild undersaturation will preferentially dissolve these disordered regions, enabling increased order upon reprecipitation. Furthermore, our simulations predict that frequent cycling of a solution between supersaturation and undersaturation can accelerate dolomite growth by up to seven orders of magnitude. We validated our theory with in situ liquid cell transmission electron microscopy, directly observing bulk dolomite growth after pulses of dissolution. This mechanism explains why modern dolomite is primarily found in natural environments with pH or salinity fluctuations. More generally, it reveals that the growth and ripening of defect-free crystals can be facilitated by deliberate periods of mild dissolution.

Kim, Joonsoo↗
DOE OSTI2025

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↗