Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Molecular Dynamic Simulations of Interaction of an AFM Probe with the Surface of an SCN Sample

Molecular dynamic (MD) simulations is conducted in order to estimate forces of probe-substrate interaction in the Atomic Force Microscope (AFM). First a review of available molecular dynamic techniques is given. Implementation of MD simulation is based on an object-oriented code developed at the University of Delft. Modeling of the sample material - succinonitrile (SCN) - is based on the Lennard-Jones potentials. For the polystyrene probe an atomic interaction potential is used. Due to object-oriented structure of the code modification of an atomic interaction potential is straight forward. Calculation of melting temperature is used for validation of the code and of the interaction potentials. Various fitting parameters of the probe-substrate interaction potentials are considered, as potentials fitted to certain properties and temperature ranges may not be reliable for the others. This research provides theoretical foundation for an interpretation of actual measurements of an interaction forces using AFM.

Bune, Adris↗
DOE OSTI2025

Atomistic characterization of hydration-dependent fuel cell ionomer nanostructure: validation by vibrational spectroscopy

The development of Nafion alternatives for fuel cells and electrolyzers requires a fundamental understanding of hydration-dependent ion-exchange site acid/base chemistry. We present here reactive force field (ReaxFF) molecular dynamics (MD) simulations of Nafion at varying molar water/ion-exchange-site ratios (λ), which we correlate to our experimental and density functional theory-based vibrational spectra. ReaxFF describes the formation and breaking of covalent bonds, enabling simulations of proton exchange between sulfonic acid/sulfonate groups and water/hydronium. Our MD simulations determine the λ-dependent equilibrium proportions of protonated and deprotonated sites. We find that protonated sites persist across all λ from 0 to 20, challenging the widely accepted notion that all sites are ionized above a threshold λ value (e.g., 3 or 4). Our simulations generate hundreds of realistic exchange site environments, the characterization of which are based on 6 Å radii sulfur-centered ‘inner-spheres’. These inner- and outer-sphere regions elucidate an interplay of stereoelectronic factors that influence protonation states, including the number of inner-sphere waters (Λ). Our simulations produce broad λ-dependent distributions of Λ values, representing non-uniform exchange site hydration. Moreover, we demonstrate that these distributions are specific to both protonated and deprotonated sites, with significant overlap between the two distributions for all non-zero λ. These distributions underpin the IR spectra of hydrated membranes, with each exchange site contributing an IR spectrum characteristic of its protonation state and Λ. We expect that these nanostructural characterizations of Nafion exchange sites will contribute to the development of new ionomers.

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

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2023

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

Here, we use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100–200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Primary defect production from molecular dynamics simulations of high-energy displacement cascades in NbMoTaW alloys

In this work, we report on large-scale molecular-dynamics (MD) simulations of displacement cascades in equiatomic NbMoTaW alloys at PKA energies ranging from 0.15 to 150 keV. We find defect production to be strongly dependent on recoil energy, scaling sublinearly up to 10 keV, and linearly thereafter. We find the sublinear regime to be defined by low values of surviving Frenkel pairs, typically found as isolated point defects or small defect clusters, while at higher recoil energies dense cascades become more frequent, leading to splitting into subcascades and the production of relatively large prismatic-dislocation loops with ⟨111⟩ and ⟨001⟩ Burgers vectors. These loops immobilize large fractions of defects, leading to a rapid growth of the number of surviving defects in the linear regime. We also anneal post-cascade defect configurations using object-kinetic Monte Carlo (OKMC) simulations to account for intracascade recombination on time scales not accessible to MD simulations. Cascade annealing is strongly temperature dependent, with the OKMC simulations only showing significant recovery at 1000 K but not below. Our results are in general agreement with existing published data for refractory concentrated alloys.

Zhou, Xinran [Lawrence Berkeley National Laborator↗
DOE OSTI2024

Amantadine interactions with phase separated lipid membranes

Amantadine, a small amphilphic organic compound that consists of an adamantane backbone and an amino group, was first recognized as an antiviral in 1963 and received approval for prophylaxis against the type A influenza virus in 1976. Since then, it has also been used to treat Parkinson’s disease-related dyskinesia and is being considered as a treatment for corona viruses. Since amantadine usually targets membrane-bound proteins, its interactions with the membrane are also thought to be important. Biological membranes are now widely understood to be laterally heterogeneous and certain proteins are known to preferentially co-localize within specific lipid domains. Does amantadine, therefore, preferentially localize in certain lipid composition domains? To address this question, here, we studied amantadine’s interactions with phase separating membranes composed of cholesterol, DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine), POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine), and DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine), as well as single-phase DPhPC (1,2-diphytanoyl-sn-glycero-3-phos-phocholine) membranes. From Langmuir trough and differential scanning calorimetry (DSC) measurements, we determined, respectively, that amantadine preferentially binds to disordered lipids, such as POPC, and lowers the phase transition temperature of POPC/DSPC/cholesterol mixtures, implying that amantadine increases membrane disorder. Further, using droplet interface bilayers (DIBs), we observed that amantadine disrupts DPhPC membranes, consistent with its disordering properties. Finally, we carried out molecular dynamics (MD) simulations on POPC/DSPC/cholesterol membranes with varying amounts of amantadine. Consistent with experiment, MD simulations showed that amantadine prefers to associate with disordered POPC-rich domains, domain boundaries, and lipid glycerol backbones. Since different proteins co-localize with different lipid domains, our results have possible implications as to which classes of proteins may be better targets for amantadine.

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

ddcMD converter

This Python package is developed for converting GROMACS files to ddcMD inputs. The code makes it easier for ddcMD users to convert GROMACS MD simulations into ddcMD. It is used in the MuMMI framework for multiscale MD simulations.

Glosli, JamesN↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pretreatment with Sodium Methyl Mercaptide Increases Carbohydrate Yield during Kraft Pulping

Kraft pulping is the predominant technology in the pulp and paper industry for removing lignin from wood carbohydrates to produce paper, board, packaging, tissue, and specialty cellulose. However, the kraft process is energy intensive and expensive, and its yield is limited by the degradation of carbohydrates. Pretreatment can increase carbohydrate yield by limiting degradation via primary peeling of reducing end groups. However, protection of galactoglucomannan (GGM), the primary hemicellulose component of softwood, is minimal when conventional pretreatments are used. Here, we investigated the effectiveness of sodium methyl mercaptide pretreatment on southern pine wood chips under a range of experimental conditions. We found that pretreatment of biomass with 4.38% sodium methyl mercaptide at pH 12 and 105 °C for 60 min provided small but significant increases in xylan and cellulose yields relative to control conditions, but preservation of GGM was minimal. Furthermore, to provide insight into molecular-scale details of primary peeling, pretreatment, and alkaline hydrolysis, we performed classical molecular dynamics (MD) simulations under selected process conditions and quantum mechanical (QM) calculations of selected reactions. MD simulations showed that C1 of the GGM reducing end is more readily accessible by HO – and CH 3 S – ions than in cellulose. The free energy barrier for peeling calculated with QM is lower for GGM than for cellulose, indicating increased susceptibility to peeling. In addition, we found that GGM may be more susceptible to internal chain cleavage than cellulose. Thus, even though reducing end groups may be protected initially through pretreatment, new unprotected reducing end groups may be generated through alkaline hydrolysis. Taken together, these findings show the promise of methyl mercaptide as a pretreatment technology for cellulose retention and also provide molecular insight for improving its effectiveness toward GGM.

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

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗
DOE OSTI2021

Prediction of solubility parameters of lignin and ionic liquids using multi-resolution simulation approaches

The solubility parameter (SP) of a molecular species is a vital feature that indicates polarity and quantifies the ‘like-seeks-like’ principle, which is used in chemistry to screen solvents for dissolution. Recent studies have demonstrated that ionic liquids (ILs) and deep eutectic solvents (DESs) efficiently solubilize lignocellulosic biomass and promote enzymatic saccharification into sugars used for the production of biofuels and value-added chemicals. Understanding the solubility of plant biopolymers, particularly lignin, in ILs and DESs is critical for selecting candidate ILs and DESs for biomass pretreatment; however, experimentally measuring SPs is challenging. Thus, the present study investigates lignin dissolution mechanisms in IL/DES and prediction of the solubility parameters (Hildebrand and Hansen) of lignin, ILs, and DESs using multi-resolution simulation approaches. Here, solubility parameters of the studied compounds were predicted using molecular dynamics (MD) simulations, and the SP of lignin was determined to be 23–27 MPa 1/2 , which was close to the polymeric lignin solubility parameter (24.3–25.5 MPa 1/2 ). The SPs of ILs namely [Ch][Lys], [Ch][Oct], and [Emim][Lys] were predicted to be ~26 MPa 1/2 , which is close to lignin's SPs and resulted in increased biomass delignification. The MD simulated SPs were validated by both the COSMO-RS model and experimental investigations, with the results showing a close agreement between the predicted and experimentally obtained SPs. In addition, the enthalpy of vaporization (ΔH vap ) of ILs/DESs was predicted based on the potential energy of the system, and the ΔH vap of ILs/DESs was around 40–65 kcal mol –1 , which is 5–8 times higher than that of traditional organic solvents.

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

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating

The voltage-dependent anion channel (VDAC) is a β-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by “gating,” i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of β-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC’s voltage sensitivity. Using the analysis of over 40 μs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the β-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the β-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 μs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that β-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗