Search NASA⌕ Search

SEARCH · Search NASA

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 217 records · Page 12

Atomic-Scale Dynamics of Five-Fold Twin Mediated Coalescence: Pathway-Dependent and Defect-Governed Nonclassical Growth Mechanisms

Defective crystals with distinct properties have been discovered in many systems. However, the growth mechanism of defective crystals is still poorly understood. Here, in this work, using a 5-fold twinned gold (Au) nanocrystal (NC) as a model system, three new coalescence pathways involving detwinning or twinning have been uncovered through atomic-scale dynamic observations in an aberration-corrected transmission electron microscope coupled with atomistic simulations. This demonstrates that beyond crystal size, coalescence growth dynamics involving 5-fold twins (5-FTs) are highly dependent on crystal defect density and the approach pathways of the crystals. When a 5-FT encounters a smaller 5-FT or a smaller NC in a face-to-face way, a new, larger 5-FT is produced at a relatively fast coalescence growth rate; while in a corner-to-corner way, the coalescence dynamics are more retarded and sluggish, which is conducive to the formation of complex multitwined structures rather than 5-FTs. This highlights that the planar defect density and crystal approach pathway influence the coalescence dynamics of crystals containing 5-FT. Moreover, a column-by-column grain boundary (GB) migration mechanism, which results in bent GBs, was also unveiled during the crystal coalescence process. These results enrich the general understanding of the crystal growth theory and provide new insights into the controllable fabrication of 5-FTs by crystal coalescence mechanisms.

crystallization↗

Ab Initio Structures and Energetics of Hydrated Flat and Terrace-Step Surfaces of Forsterite (Mg 2 SiO 4 )

Forsterite (Mg 2 SiO 4 ), a model divalent metal silicate mineral, has been extensively studied in the context of mineral carbonation. Although dissolution is a key step in this process, the mechanisms by which forsterite dissolves under high CO 2 conditions remain poorly understood. Atomistic simulations could aid in exploring these mechanisms, but it is essential first to understand the structures and energetics of the relevant forsterite surfaces. We present an ab initio study of the structure and surface energy at 0 K of the flat $(010), (110), (001), (111), (021), (101)$ and $(120)$ faces of forsterite using the density functional PBE Hamiltonian and a plane-wave basis set. Dry surfaces became stabilized upon hydration through the formation of bonds between surface Mg and O from water, as well as by the formation of hydrogen bonds. According to surface energy values, the stability order of the hydrated forsterite faces was found to be $(120) < (101) < (021) < (111) < (001) < (110) < (010)$. We also investigated the energetics of the terrace-step $(0\bar{41})$ surface as a model site for forsterite dissolution. Among all the facets, the $(0\bar{41})$ surface is the least stable termination in water. Hydration of Mg atoms on the $(0\bar{41})$ surface increases their susceptibility to dissolution. The presence of a step and its hydration destabilizes the terraces, making step retreat more likely than a dissolution front advancing along the [010] direction. This research will support future simulations to investigate forsterite dissolution in water under CO 2 -rich conditions.

PBE Hamiltonian↗

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION↗

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↗

Polariton-assisted incoherent to coherent excitation energy transfer between colloidal nanocrystal quantum dots

Here, we explore the dynamics of energy transfer between two nanocrystal quantum dots placed within an optical microcavity. By adjusting the coupling strength between the cavity photon mode and the quantum dots, we have the capacity to fine-tune the effective coupling between the donor and acceptor. Introducing a non-adiabatic parameter, γ, governed by the coupling to the cavity mode, we demonstrate the system’s capability to shift from the overdamped Förster regime (γ $\ll$ 1) to an underdamped coherent regime (γ $\gg$ 1). In the latter regime, characterized by swift energy transfer rates, the dynamics are influenced by decoherence time. To illustrate this, we study the exciton energy transfer dynamics between two closely positioned CdSe/CdS core/shell quantum dots with sizes and separations relevant to experimental conditions. Employing an atomistic approach, we calculate the excitonic level arrangement, exciton–phonon interactions, and transition dipole moments of the quantum dots within the microcavity. These parameters are then utilized to define a model Hamiltonian. Subsequently, we apply a generalized non-Markovian quantum Redfield equation to delineate the dynamics within the polaritonic framework.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗

Computational materials assessment of the D/Li-stripping neutron source as a prototypical facility for fusion materials testing

As the US fusion materials community awaits the selection and design of a fusion prototypical neutron source (FPNS), a risk reduction exercise has been conducted to (i) provide an updated materials performance evaluation using state-of-the-art computational materials modeling, (ii) expand on legacy analysis based on pure Fe to other relevant fusion structural materials types, and (iii) ensure that materials response under FPNS operational conditions is consistent with referential fusion reactor conditions. The current paper describes the efforts undertaken to assemble a comprehensive computational methodology that includes neutronics, primary damage calculations, atomistic simulations of displacement cascades, chemical inventory evolution calculations, and a computational thermodynamic analysis of emerging phases during irradiation. Our work extends existing studies in pure Fe to reduced-activation ferritic/martensitic steels, tungsten, silicon carbide, and vanadium alloys. We focus on the single-beam deuteron/lithium-stripping neutron source behind the IFMIF-DONES concept, which we assess against ITER, two DEMO designs, and an ideal pure 14-MeV flux. Our analysis indicates that, within standard uncertainties inherent to the models employed, the DONES concept adequately captures fusion conditions in the four materials analyzed. Our work is intended as a comprehensive irradiation damage analysis of fusion-representative neutron sources, to be used for further neutron source evaluation and fusion facility operation.

Marian, Jaime [Univ. of California, Los Angeles, C↗

The evolution of model Rh/Fe 3 O 4 (001) catalysts in hydrogen environments

Single metal atoms dispersed on oxides are a new emerging class of catalysts owing to their unique electronic and chemical properties. Here, in this study, we have prepared a series of model single-atom catalysts possessing well-characterized Rh sites that include Rh adatoms (Rh ad ), mixed surface layers with octahedrally-coordinated Rh (Rh oct ), as well as metallic Rh clusters and nanoparticles (Rh met ) on Fe 3 O 4 (001). Using X-ray photoelectron spectroscopy (XPS) and scanning tunneling microscopy (STM), we investigated the activity of such model systems towards H2 and their stability in reducing environments. Our results show that the atomically dispersed Rhad and Rhoct species do not activate H 2 , which would result in the formation of surface hydroxyls on Fe 3 O 4 (001). In contrast, the presence of Rh met in H 2 results in the formation of hydroxyls and subsequent etching of the Fe 3 O 4 (001) at higher temperatures (≥ 500 K) due to water formation via the Mars-van Krevelen mechanism. Additionally, such surface etching leads to the release of the Rh oct from the surface lattice and their sintering to Rh met . To bridge the material gap between the surface science models and high surface area catalysts, we perform parallel studies on powder Rh/Fe 3 O 4 catalysts. The XPS characterization shows remarkable similarities between these systems. Further, our surface science studies provide an atomistic picture of the behavior of high surface area catalysts in the H 2 atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chromium segregation-induced oxide evolution in Ni-10Cr alloys during high-temperature oxidation

The oxidation behavior of a Ni-10(wt%)Cr alloy under high-temperature O 2 conditions is investigated using transmission electron microscopy and first-principles calculations. Results reveal that chromium segregation plays a central role in driving the evolution of complex oxide phase structures during oxidation. At low Cr concentrations, Cr preferentially segregates to NiO grain boundaries or internal pores, substituting for Ni atoms and forming Ni(Cr)O solid solutions. As Cr content increases, enhanced diffusion promotes Cr penetration into the NiO lattice, leading to the formation of multiphase oxide structures. First-principles modeling corroborates these findings: at low Cr concentrations, Cr atoms favor surface and grain-boundary segregation, while higher concentrations lead to Cr aggregation within the NiO bulk. Furthermore, the integrated experimental-theoretical approach provides atomistic insights into Cr-mediated mass transport mechanisms during alloy oxidation and offers valuable guidance for controlling oxide growth kinetics and phase stability in Ni-Cr alloys, with implications for improving oxidation resistance in high-temperature structural applications.

36 MATERIALS SCIENCE↗

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗

Atomistic Mechanisms of Stress-Dependent Molten Salt Corrosion in NiCr Alloys

Ni-based structural alloys in molten salt environments often experience simultaneous mechanical loading and corrosive attack, yet the mechanisms governing stress-corrosion interactions remain unclear. Prior studies largely emphasize tensile stress, while the role of compressive stress has received limited attention. Here, reactive molecular dynamics simulations are used to investigate the coupled effects of applied strain and corrosion in Ni 0.75 Cr 0.25 exposed to molten FLiNaK at 800 °C. A Σ5(210) grain boundary model is subjected to tensile (+4%) to compressive (−4%) uniaxial strains, and corrosion behavior is evaluated through fluorine adsorption, charge redistribution, and grain boundary evolution. Tensile strain accelerates intergranular corrosion susceptibility by reducing local atomic packing through elastic dilation and increasing excess free volume at the grain boundary, which enhances atomic mobility and salt infiltration. In contrast, compressive strain can suppress corrosion by promoting the formation of a ridge-like surface layer along the grain boundary, limiting salt access to the underlying alloy. These results provide atomistic insight into how stress states influence grain boundary corrosion in molten salts.

36 - MATERIALS SCIENCE↗

Energetics of the nucleation and glide of disconnection modes in symmetric tilt grain boundaries

Grain boundaries (GBs) evolve by the nucleation and glide of disconnections, which are dislocations with a step character. In this work, motivated by recent success in predicting GB properties such as the shear coupling factor and mobility from the intrinsic properties of disconnections, we develop a systematic method to calculate the energy barriers for the nucleation and glide of individual disconnection modes under arbitrary driving forces and a quasi-2D setting. This method combines tools from bicrystallography to enumerate disconnection modes and the Nudged elastic band (NEB) method to calculate their energetics, yielding minimum energy paths and atomistic mechanisms for the nucleation and glide of each disconnection mode. We apply the method to accurately predict shear coupling factors of $[001]$ symmetric tilt grain boundaries in Cu. Particular attention is paid to the boundaries where the dislocation-based disconnection nucleation model produces incorrect nucleation barriers. We demonstrate that the method can accurately compute energy barriers and predict shear-coupling factors in the low-temperature regime. For certain disconnection modes in which the assumptions underlying our method do not hold, we report upper bounds on the energy barriers for disconnection nucleation and glide. In addition, the NEB trajectories reveal interesting phenomena such as the dissociation of a higher energy mode into lower energy modes, and in some cases, shear coupling being mediated by partial disconnections, wherein the GB structure temporarily changes to a metastable state before reverting back to its original structure. Graphical abstract

36 MATERIALS SCIENCE↗

Learning from metastable symmetric-tilt grain boundaries using physics-based descriptors

Grain boundaries (GBs) govern critical properties of polycrystalline materials. Although significant advancements have been made in characterizing minimum energy and ordered GBs, real GBs are seldom found in such well-defined states. This diversity of atomic arrangements in metastable states makes it challenging to establish structure-property relationships with physical insights. Here, to address this challenge, we use data-driven methods to explore these relationships and examine the underlying physics. In this study, we utilize a large atomistic database (~5000) of minimum energy and metastable states of symmetric-tilt copper GBs, combined with physically motivated local atomic environment (LAE) descriptors [strain functional descriptors (SFDs)], to predict GB properties and gain physical insights. Our regression models exhibit robust predictive capabilities using only 19 descriptors, generalizing to atomic environments in nanocrystals. A significant highlight of our work is the integration of an unsupervised method with SFDs to elucidate LAEs at GBs and their role in determining properties. The model, trained on these minimum energy and metastable GBs using SFDs, predicts the properties of unseen nanocrystals with good accuracy. Our research underscores the role of a physics-based representation of LAEs and the efficacy of data-driven methods in establishing GB structure-property relationships.

36 MATERIALS SCIENCE↗

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology↗

OpenFerro v0.1.0

OpenFerro is a Python package for on-lattice atomistic dynamics simulation of ferroic materials. OpenFerro is based on JAX, a high-performance linear algebra package supporting auto-differentiation and GPU acceleration. OpenFerro is designed to minimize the effort required to build on-lattice Hamiltonian models, and to perform molecular dynamics (MD) and Landau-Lifshitz-Gilbert simulations. Unlike existing codes, OpenFerro provides a unified interface to model different types of local order parameters.

Xie, Pinchen [Lawrence Berkeley National Laborator↗

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗

Simulation of ultrafast transient absorption spectra of a perylene-based light harvesting antenna

Atomistic simulations of photo-induced responses in artificial light-harvesting molecular systems help to reveal the mechanisms of ultrafast intramolecular energy transfer between individual chromophores. These light-induced processes mimic the primary events occurring in natural photosynthesis. Modeling studies contribute to the design of more efficient molecular architectures enabling performance optimization for applications in light harvesting, energy conversion, and optoelectronics. Within this context, the direct comparison between simulated and experimental transient absorption pump–probe (TA-PP) spectra are especially valuable for validating theoretical approaches and deepening mechanistic understanding. Herein, we investigate the photoinduced dynamics of an antenna system composed of two naphthalene monoimides donor units covalently linked to a perylene derived acceptor. Following photoexcitation, the exciton rapidly self-traps on one of the donor units. Thereafter, efficient ultrafast energy transfer to the acceptor unit takes place via two possible pathways: either through transient exciton localization on the second donor unit or by direct transfer to the acceptor. The simulated TA-PP spectra clearly capture these distinct energy transfer pathways and enable a detailed comparison of their relative efficiencies. This highlights the system's potential for tunable exciton dynamics towards advancing light-harvesting and optoelectronic molecular materials.

36 MATERIALS SCIENCE↗

Selected magnetoelectrical phonon stiffening in Cr 2 ⁢O 3

A thorough understanding of lattice dynamics, particularly how they couple with spin through spin-phonon interactions, is crucial for pioneering new spin-caloritronic applications. Here, despite extensive studies on Cr 2 ⁢O 3 , the origin of nonlinear thermal expansions and the unusual stiffening of optical phonon modes is still elusive. Through inelastic neutron scattering and atomistic calculations, we found that these behaviors can be attributed to the renormalization of electron states owing to the magnetic transition, and we ruled out the effects from thermal expansion, phonon anharmonicity, magnetostriction, or electron-phonon interactions. Importantly, our quantitative modeling suggests that the common belief that dynamic spin-phonon interactions are the origin of anomalous phonon energy stiffening in magnetic compounds does not apply in Cr 2 ⁢O 3 .

DFT+U↗