Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2025

Cross-scale covariance for material property prediction

A simulation can stand its ground against an experiment only if its prediction uncertainty is known. The unknown accuracy of interatomic potentials (IPs) is a major source of prediction uncertainty, severely limiting the use of large-scale classical atomistic simulations in a wide range of scientific and engineering applications. Here we explore covariance between predictions of metal plasticity, from 178 large-scale (~10 8 atoms) molecular dynamics (MD) simulations, and a variety of indicator properties computed at small-scales (≤10 2 atoms). All simulations use the same 178 IPs. In a manner similar to statistical studies in public health, we analyze correlations of strength with indicators, identify the best predictor properties, and build a cross-scale “strength-on-predictors” regression model. This model is then used to estimate regression error over the statistical pool of IPs. Small-scale predictors found to be highly covariant with strength are computed using expensive quantum-accurate calculations and used to predict flow strength, within the statistical error bounds established in our study.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Atomistic Simulation of Glasses and Amorphous Materials: Challenges and Opportunities for the Next Decade

Atomistic simulations have become indispensable tools for understanding glass structure, dynamics, and properties, yet persistent challenges limit their predictive power. This perspective examines three interconnected issues, namely glass formation procedures, interatomic potential development, and machine learning applications, which emerged from the 5th International Workshop on Challenges of Atomistic Simulations of Glasses and Amorphous Materials. We identify convergent community priorities for (i) standardized validation protocols, (ii) curated benchmark datasets with complete metadata, and (iii) open repositories for glasses. A systematic was forward is provided by a hierarchical validation framework for assessing the structural fidelity, property prediction, and behavioral realism of simulation techniques. Looking ahead, transformative advances are promised by the fusion of classical techniques with machine learning based approaches, for instance, by integrating swap Monte Carlo with machine-learning (ML) potentials, leveraging foundation models through transfer learning, and finetuning ML potentials with experimental data. Progress depends on the community committing to validated models, reproducible protocols, and sustained data sharing.

Krishnan, N. M. Anoop↗
DOE OSTI2025

First-principles effective Hamiltonian for finite-temperature modeling of nonperovskite ferroelectrics

First-principles-based effective Hamiltonian techniques have been widely employed for over three decades to investigate ferroelectricity and related phenomena in perovskite materials. These techniques offer high accuracy, transferability, compatibility with various finite-temperature algorithms, computational efficiency, and ease in incorporating interactions with external fields. They have been adapted to study diverse phenomena, ranging from topological dipole patterns in ferroelectric nanostructures to multicaloric effects. In this work, we develop an effective Hamiltonian for the nonperovskite ferroelectric HfO 2 (hafnia). Applying this methodology to explore the finite-temperature and finite-electric-field properties of ferroelectric hafnia revealed (1) exceptionally large intrinsic coercive fields, an order of magnitude higher than those observed in perovskite ferroelectrics; (2) their atomistic origin; and (3) the existence of a regime where the relationship between the coercive field and the energy barrier for polarization reversal is counterintuitive. Here, these developments could accelerate progress both in methodological advancements for simulating ferroics and in the atomistic understanding of a broad range of ferroelectrics.

Electric polarization↗
DOE OSTI2022

Revealing hidden defects through stored energy measurements of radiation damage

With full knowledge of a material’s atomistic structure, it is possible to predict any macroscopic property of interest. In practice, this is hindered by limitations of the chosen characterization techniques. For example, electron microscopy is unable to detect the smallest and most numerous defects in irradiated materials. Instead of spatial characterization, we propose to detect and quantify defects through their excess energy. Differential scanning calorimetry of irradiated Ti measures defect densities five times greater than those determined using transmission electron microscopy. Our experiments also reveal two energetically distinct processes where the established annealing model predicts one. Molecular dynamics simulations discover the defects responsible and inform a new mechanism for the recovery of irradiation-induced defects. The combination of annealing experiments and simulations can reveal defects hidden to other characterization techniques and has the potential to uncover new mechanisms behind the evolution of defects in materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2021

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of↗
NASA NTRSDate not supplied

Molten Salt Electrolytes: Computational Analysis of Transport, Electrochemical Properties and Designing New Mixtures

Electric aircraft propulsion has gained a traction over the last decade due to possible high-energy density electrochemistries with reasonable cycle life and identification of non-flammable chemistries that can guarantee much better safety. Commercial Li-ion batteries cannot reach such high capacities because of weight limitations that arise from the widely used intercalation electrodes. Also, use of organic electrolytes make them highly susceptible to fire on exposure to air and humidity. Currently, use of Lithium metal anode along with conversion electrochemistries such as Li-O2 and Li-S are being pursued to achieve such high energy densities. Safer inorganic solid-state electrolytes, recently discovered Water-in-Salt electrolytes, molten salt electrolytes are some of the alternatives being pursued for a safer/non-flammable battery. Of these safer alternatives, molten-salt electrolytes offer some attractive properties: liquid at operating temperatures resulting in better electrode-wetting, stable interface with Li-metal anodes and good ionic conductivity; their only drawback being high operating temperatures. In this talk, we will discuss some of the computational studies, driven via atomistic simulations, of the properties of molten-salt electrolytes including transport mechanism and interface stability. We will also discuss predicting electrochemical properties of these high-temperature electrolytes. Further, we will discuss a thermodynamic approach to predicting new molten-salt eutectics with lower melting points.

Molten Salt Electrolytes↗
DOE OSTI2023

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale . Here, our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

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

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗
DOE OSTI2025

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Domain-wall-mediated polarization switching in ferroelectric AlScN: Strain relief and field-dependent dynamics

While scandium-doped aluminum nitride (AlScN) exhibits robust ferroelectricity and excellent thermal stability, its utility is limited by an exceptionally high coercive field (𝐸 𝑐 ) for polarization switching. Unraveling the atomistic switching dynamics is therefore critical for tailoring 𝐸 𝑐 . Here, in this study, we combine density functional theory and machine-learning molecular dynamics to elucidate the polarization switching mechanisms in AlScN over various Sc concentrations and applied electric fields. We find that excessive lattice strain strictly prohibits collective polarization switching, but the preexisting domain walls relieve strain and lead to a distinct switching dynamics—dictating a field-dependent switching mechanism. At low electric fields, switching occurs via gradual domain-wall propagation consistent with the Kolmogorov-Avrami-Ishibashi model. In contrast, high fields stimulate additional nucleation, driving a rapid, homogeneous reversal process described by the simultaneous nonlinear nucleation and growth model. These findings highlight the critical role of domain-wall dynamics and suggest domain engineering as a viable strategy to tailor coercive fields in AlScN and related ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

SCOMAP-XD : atomistic deuterium contrast matching for small-angle neutron scattering in biology

The contrast-variation method in small-angle neutron scattering (SANS) is a uniquely powerful technique for determining the structure of individual components in biomolecular systems containing regions of different neutron scattering length density ρ . By altering the ρ of the target solute and the solvent through judicious incorporation of deuterium, the scattering of desired solute features can be highlighted. Most contrast-variation methods focus on highlighting specific bulk solute elements, but not on how the scattering at specific scattering vectors q , which are associated with specific structural distances, changes with contrast. Indeed, many systems exhibit q -dependent contrast effects. Here, a method is presented for calculating both bulk contrast-match points and q -dependent contrast using 3D models with explicit solute and solvent atoms and SASSENA , an explicit-atom SANS calculator. The method calculates the bulk contrast-match points within 2.4% solvent D 2 O accuracy for test protein–nucleic acid and lipid nanodisc systems. The method incorporates a general model for the incorporation of deuterium at non-exchangeable sites that was derived by performing mass spectrometry on green fluorescent protein. The method also decomposes the scattering profile into its component parts and identifies structural features that change with contrast. The method is readily applicable to a variety of systems, will expand the understanding of q -dependent contrast matching and will aid in the optimization of next-generation neutron scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗