Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic level 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Interfacial Defect of Lithium Metal in Solid-State Batteries

All-solid-state battery with Li metal anode is a promising rechargeable battery technology with high energy density and improved safety. Currently, the application of Li metal anode is plagued by the failure at the interfaces between lithium metal and solid electrolyte (SE). However, little is known about the defects at Li–SE interfaces and their effects on Li cycling, impeding further improvement of Li metal anodes. Herein, by performing large-scale atomistic modeling of Li metal interfaces with common SEs, we discover that lithium metal forms an interfacial defect layer of nanometer-thin disordered lithium at the Li–SE interfaces. Further, this interfacial defect Li layer is highly detrimental, leading to interfacial failure such as pore formation and contact loss during Li stripping. By systematically studying and comparing incoherent, coherent, and semi-coherent Li–SE interfaces, we find that the interface with good lattice coherence has reduced Li defects at the interface and has suppressed interfacial failure during Li cycling. Our finding discovered the critical roles of atomistic lithium defects at interfaces for the interfacial failure of Li metal anode, and motivates future atomistic-level interfacial engineering for Li metal anode in solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

CHARMM-GUI Bicelle Builder : An Extension of Membrane Builder for Modeling and Simulation of Bicelle Systems

Membrane mimetics, such as detergent micelles, nanodiscs, and amphipol complexes, which can provide membrane-like environments while retaining small and soluble features, have been utilized to study membrane proteins. A bicelle, composed of varying lipids and detergents, is a useful membrane mimetic because the lipid-to-detergent ratio, the q-value, can be adjusted to alter the properties of the aggregate, including the thickness and size of the bicelle. However, building a bicelle model for modeling and simulation studies requires nontrivial efforts, even for experts. We introduce CHARMM-GUI Bicelle Builder, a web-based platform that can generate various all-atom bicelle systems via a graphical user interface with all available lipids and detergents in Membrane Builder. To illustrate and validate Bicelle Builder with practical systems, we have modeled and simulated pure bicelles consisting of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids with 1,2-dihexanoyl-sn-glycero-3-phosphocholine (C6DHPC) detergents and protein–bicelle complexes, composed of DMPC with C6DHPC, foscholine-10 (FOS10), and lysophosphatidylcholine-12 (LPC12) detergents. Our simulation results indicate that Bicelle Builder can generate reliable and robust bicelle models with and without proteins that retain DMPC bilayer characteristics. Bicelle Builder is expected to help researchers better understand not only bicelles themselves but also atomistic-level structures of protein–bicelle complexes that are often difficult to access through experimental approaches.

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

Ion and Water Dynamics in the Transition from Dry to Wet Conditions in Salt-Doped PEG

The influence of the water content on ion and water transport mechanisms in polymer membranes under low to moderate hydration conditions remains poorly understood. In this study, we combine ion and water diffusivity (PFG-NMR) measurements with atomistic molecular dynamics simulations to better understand transport processes in hydrated salt-doped poly(ethylene glycol). Above the water percolation threshold, the experimental and simulated diffusivities are in good agreement with the free volume transport models. At low hydration levels, unlike dry systems, ion diffusion cannot be described by polymer segmental dynamics alone. We rationalize such observations by the interplay between ion–water and ion–polymer solvation of cations and between ion–water and cation–anion interactions for anions. Further, we demonstrate that a two-state model combining ion–water solvation and free volume transport can describe water dynamics across the entire hydration range of interest. Here, our findings provide a more encompassing analysis of ion and water transport in hydrated polyelectrolytes, specifically in the low hydration regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dislocation loop evolution in Kr-irradiated ThO 2

Here, the early stage of microstructural evolution of ThO 2 , under krypton irradiation at 600, 800, and 1000°C, was investigated using in situ transmission electron microscopy (TEM). Dislocation loops grew faster, whereas their number density decreased with increasing irradiation temperature. Loop density was found to decrease with ion dose. Interstitial dislocation loops, including Frank loops with Burgers vector of $\textit{a}/3\langle111\rangle$ and perfect loops with Burgers vector of $\textit{a}/2\langle110\rangle$, were determined by traditional TEM and atomic resolution–scanning TEM techniques. Atomistic and mesoscale level modeling are performed to interpret experimental observations. The migration energy barriers of defects in ThO 2 were calculated using density-functional theory. The energetics of different dislocation loop types were studied using molecular dynamics simulations. Loop density and diameter were analyzed using a kinetic rate theory model that considers stoichiometric loop evolution. This analysis reveals that loop growth is governed by the mobility of cation interstitials, whereas loop nucleation is determined by the mobility of anion defects. Lastly, a rate theory model was used to extract the diffusion coefficients of thorium interstitials, oxygen interstitials, and vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2022

The heterogeneous nature of mechanically accelerated grain growth

Abstract While grain growth is traditionally viewed as a purely thermally driven process, nanocrystalline metals can undergo grain growth under mechanical loads, even at room temperature. We performed a detailed atomistic study of the heterogeneous nature of mechanically accelerated grain growth in a polycrystalline Pt nanowire. Using molecular dynamics simulations, we compared the grain-growth behavior of individual grains during tensile and shear cyclic loading, for three different equivalent strain levels, and at two temperatures. Pure thermal grain growth with no mechanical loading provided a baseline reference case. On average, grains that were already susceptible to thermal grain growth were stimulated to grow faster with mechanical loading, as expected. However, when analyzed on a grain-by-grain basis, the results were far more complex: grains that grew fastest under one stimuli were less accelerated under other stimuli. Even when the magnitude of loading changed, the relative growth of individual grains was distorted. We interpret this complexity from the perspective of superimposed growth mechanisms.

Chen, Elton Y.↗
DOE OSTI2025

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
DOE OSTI2025

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

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

Theory of Defect-Induced Crystal Field Perturbations in Rare-Earth Magnets

We present a theory describing the single-ion anisotropy of rare-earth (RE) magnets in the presence of point defects. Taking the RE-lean 1 ∶ 12 magnet class as a prototype, we use first-principles calculations to show how the introduction of Ti substitutions into SmFe 12 perturbs the crystal field, generating new coefficients due to the lower symmetry of the RE environment. We then demonstrate that these perturbations can be described extremely efficiently using a screened point charge model. We provide analytical expressions for the anisotropy energy that can be straightforwardly implemented in atomistic spin dynamics simulations, meaning that such simulations can be carried out for an arbitrary arrangement of point defects. The significant crystal field perturbations calculated here demonstrate that a sample that is single phase from a structural point of view can nonetheless have a dramatically varying anisotropy profile at the atomistic level if there is compositional disorder, which may influence localized magnetic objects like domain walls or skyrmions. Published by the American Physical Society 2024

25 ENERGY STORAGE↗
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2021

Screw dislocation impingement and slip transfer at fcc-bcc semicoherent interfaces

Understanding of dislocation-interface interactions at the atomistic level could provide useful insights into the mechanisms behind slip transfer, a key information for higher length-scale formulations. Here, a systematic study of screw dislocation impingement on both Fe/Al and Fe/Cu interfaces with Kurdjumov-Sachs orientation relationship has been carried out using atomistic simulations. For the “strong” Fe/Al interface, there is appreciable dislocation delocalization, in which the change in dislocation content provided by Nye tensor analysis. Interestingly, the delocalization is significant enough that large simulation cells must be employed in order to avoid boundary effects. When the impinging dislocations are far away from the boundaries, the large accumulation of shear strains nearby the impinging dislocations at the interface causes dislocation nucleation/emission at these places. While for the “weak” Fe/Cu interface, the dislocation delocalization at the interface is notably long-range, extending more than 20 nm away from where the initial dislocations entered the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lithiation of Sulfur-Graphene Compounds Using Reactive Force-Field Molecular Dynamics Simulations

We performed molecular dynamics simulations of lithium-sulfur-graphene compounds using reactive force fields, providing a time scale to observe atomistic features relevant to the microscopic behavior of the of the bulk of sulfur-based cathodes to be used beyond our present Li-ion batteries. The samples we used were set to realistic geometries through sophisticated protocols to simulate ultrafast reactions that occur within the picosecond range, thus allowing us to get some insights into the characteristics of the bulk material in working cathodes of Li–S batteries, which are mixed with carbon to increase the poor electronic conductivity of S. We report chemical speciation and geometrical data at atomistic levels. We observed that slowly lithiated cathodes were more stable and with higher density than those that were suddenly fully-lithiated. We did not observe molecular Li 2 S formation; however, we observed an amorphous solid arrangement with the same stoichiometry of Li and S, with S–Li–S angles of ~111° and smaller ones due to the interaction between polysulfides that did not reacted totally. In addition, graphene keeps its planar shape; however, S8 changes its shape from rings to chains. Lithiated structures are more stable with lower energies, and more close-packed structures than structures with Li already inserted.

25 ENERGY STORAGE↗
DOE OSTI2022

Structure of Water Adsorbed on Nanocrystalline Calcium Silicate Hydrate Determined from Neutron Scattering and Molecular Dynamics Simulations

Calcium silicate hydrate (C-S-H) is a disordered, nanocrystalline material that acts as a primary binding phase in Portland cement. Thin films of water are present on the surfaces and in nanopores of C-S-H, impacting many of its chemical and mechanical properties, such as ion transport, creep, or thermal behavior. Despite decades of research, a full understanding of the structural details of adsorbed, confined, and bulk water in C-S-H remains elusive. In this work, we applied a multitechnique study involving molecular dynamics (MD) simulations validated by neutron diffraction with isotopic substitution (NDIS) and X-ray scattering methods to investigate the structure of water in C-S-H and C-A-S-H (an Al-bearing, low-CO 2 C-S-H substitute). Direct comparison of NDIS data with the MD results reveals that the structure of confined and interfacial water differs significantly from the bulk water and exhibits a larger degree of mesoscale ordering for more hydrated C-S-H structures. As a result, this observation suggests an important role of water as a stabilizer of the atomistic-level structure of C-S-H.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elastic Properties and Deformation Mechanisms in the van der Waals Single-Crystalline Indium Selenide

The bulk van der Waals (vdW) single-crystalline indium selenide (β-InSe) displays exceptional plasticity in a layered crystalline form at both micro- and macroscale. However, the nanoscale origin of plasticity remains unclear. Herein, an atomic-level study on the deformation mechanisms of InSe by using first-principles calculations is reported. Remarkable anisotropic elasticity is observed in the vdW InSe layered crystal, and the stiffness is dramatically softened because of the vdW gap in the layered structure. The simulations capture the distinct fracture modes in the uniaxial tensile deformation, depending on the loading directions—brittle fracture in the [100] and [110] directions while ductile failure in the [001] direction. The InSe layered crystal structure exhibits superplastic deformability under uniaxial compression. Different transition pathways, including interlayer tangling, amorphization, and cross-linking, are tracked along respective deformation directions. The unprecedented plasticity of InSe layered crystals can be attributed to the phase transition coupled with interlayer gliding and cross-layer dislocation slipping. Here, this study deepens our understanding of the deformation mechanisms of layered materials at the atomistic level and provides insights into tailoring material properties for low-dimensional material design based on its deformation mechanisms.

2D materials↗
DOE OSTI2020

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗