Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atomistic modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Atomistic Modeling of the Interphase Structure between Alumina and Al12Fe2Cr Quasicrystal Approximant

Molecular dynamics simulations are performed to study the stability of an fcc metallic interlayer phase between an iron aluminide quasicrystal approximant (QCA) phase (84.2 Al, 5.7 Fe, 5.4 Ni, and 4.7 at.% Cr) and aluminum oxides (α-Al 2 O 3 and γ-Al 2 O 3 ). The presence of the interlayer phase was experimentally observed in some regions of the inner diameter surface of a TPBAR cladding tube. The simulations employed existing many-body (MEAM and EAM) interatomic potentials to describe atomic interactions in the interlayer and QCA phases, and Buckingham+Coulomb pair interactions for the alumina and between the alumina and the QCA (or fcc). Interface systems of QCA/alumina, fec/alumina, and QCA/fcc are modeled and the fracture energy of the interfaces is calculated. Based on the fracture energies, for the γ-Al 2 O 3 , both MEAM and EAM suggest that the fcc interlayer phase would not form. For the α-Al 2 O 3 , the MEAM predicts that the fcc phase would form, however, the EAM predicts that it forms only if the α-Al 2 O 3 is Al-terminated at the interface. For the O-terminated α-Al 2 O 3 , the EAM predicts the fcc phase would not form, consistent with the results for γ-Al 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Solvated Ionic Liquid-based Electrolytes: Experimental Characterization, Atomistic Modeling and Applications in Energy Storage

This study aimed on understanding the molecular structure and transport mechanisms of high-concentration lithium electrolytes, with a particular focus on solvate ionic liquids (SILs) based on glyme solvents, for use in next-generation, high-voltage energy storage systems. The hypothesis is that the composition of the solvents in these electrolytes can be designed and adjusted to control cation solvation, enabling charge transport via an ion-hopping mechanism. This mechanism resembles Grotthuss conduction of protons in water and is independent of electrolyte viscosity. The experimental studies employed linear and nonlinear infrared spectroscopy (FTIR and 2D IR), electrochemical measurements, and vibrational probes.

25 ENERGY STORAGE↗
NASA NTRS2002

Atomistic Modeling of Surface and Bulk Properties of Cu, Pd and the Cu-Pd System

The BFS (Bozzolo-Ferrante-Smith) method for alloys is applied to the study of the Cu-Pd system. A variety of issues are analyzed and discussed, including the properties of pure Cu or Pd crystals (surface energies, surface relaxations), Pd/Cu and Cu/Pd surface alloys, segregation of Pd (or Cu) in Cu (or Pd), concentration dependence of the lattice parameter of the high temperature fcc CuPd solid solution, the formation and properties of low temperature ordered phases, and order-disorder transition temperatures. Emphasis is made on the ability of the method to describe these properties on the basis of a minimum set of BFS universal parameters that uniquely characterize the Cu-Pd system.

Bozzolo, Guillermo↗
NASA NTRS2002

Growth of Pt/Cu(100): An Atomistic Modeling Comparison with the Pd/Cu(100) Surface Alloy

The Bozzolo, Ferrante, and Smith (BFS) method for alloys is applied to the study of Pt deposition on Cu(100). The formation of a Cu-Pt surface alloy is discussed within the framework of previous results for Pd/Cu(100). In spite of the fact that both Pd and Pt share the same basic behavior when deposited on Cu, it is seen that subtle differences become responsible for the differences in growth observed at higher cover-ages. In agreement with experiment, all the main features of Pt/Cu(100) and Pd/Cu(100) are obtained by means of a simple modeling scheme, and explained in terms of a few basic ingredients that emerge from the BFS analysis.

Demarco, Gustavo↗
NASA NTRS2003

Atomistic Modeling of Nanostructures via the BFS Quantum Approximate Method

Ideally, computational modeling techniques for nanoscopic physics would be able to perform free of limitations on the type and number of elements, while providing comparable accuracy when dealing with bulk or surface problems. Computational efficiency is also desirable, if not mandatory, for properly dealing with the complexity of typical nano-strucured systems. A quantum approximate technique, the BFS method for alloys, which attempts to meet these demands, is introduced for the calculation of the energetics of nanostructures. The versatility of the technique is demonstrated through analysis of diverse systems, including multi-phase precipitation in a five element Ni-Al-Ti-Cr-Cu alloy and the formation of mixed composition Co-Cu islands on a metallic Cu(III) substrate.

Bozzolo, Guillermo↗
NASA NTRS2004

Atomistic Modeling of Pd Site Preference in NiTi

An analysis of the site subsitution behavior of Pd in NiTi was performed using the BFS method for alloys. Through a combination of Monte Carlo simulations and detailed atom-by-atom energetic analyses of various computational cells, representing compositions of NiTi with up to 10 at% Pd, a detailed understanding of site occupancy of Pd in NiTi was revealed. Pd subsituted at the expense of Ni in a NiTi alloy will prefer the Ni-sites. Pd subsituted at the expense of Ti shows a very weak preference for Ti-sites that diminishes as the amount of Pd in the alloy increases and as the temperature increases.

Bozzolo, Guillermo↗
NASA NTRS2007

Quantum Approximate Methods for the Atomistic Modeling of Multicomponent Alloys

This chapter describes the role of quantum approximate methods in the understanding of complex multicomponent alloys at the atomic level. The need to accelerate materials design programs based on economical and efficient modeling techniques provides the framework for the introduction of approximations and simplifications in otherwise rigorous theoretical schemes. As a promising example of the role that such approximate methods might have in the development of complex systems, the BFS method for alloys is presented and applied to Ru-rich Ni-base superalloys and also to the NiAI(Ti,Cu) system, highlighting the benefits that can be obtained from introducing simple modeling techniques to the investigation of such complex systems.

Bozzolo, Guillermo↗
NASA NTRS2006

Electrostatic Properties of PE and PTFE Subjected to Atmospheric Pressure Plasma Treatment; Correlation of Experimental Results with Atomistic Modeling

The use of an atmospheric pressure glow discharge (APGD) plasma was used at KSC to increase the hydrophilicity of spaceport materials to enhance their surface charge dissipation and prevent possible ESD in spaceport operations. Significant decreases in charge decay times were observed after tribocharging the materials using the standard KSC tribocharging test. The polarity and amount of charge transferred was dependent upon the effective work function differences between the respective materials. In this study, polyethylene (PE) and polytetrafluoroethylene (PTFE) were exposed to a He+O2 APGD. The pre and post treatment surface chemistry was analyzed by X-ray photoelectron spectroscopy and contact angle measurements. Semi-empirical and ab initio calculations were performed to correlate the experimental results with some plausible molecular and electronic structure features of the oxidation process. For the PE, significant surface oxidation was observed, as indicated by XPS showing C-O, C=O, and O-C=O bonding, and a decrease in the surface contact angle from 98.9 deg to 61.2 deg. For the PTFE, no C-O bonding appeared and the surface contact angle increased indicating the APGD only succeeded in cleaning the PTFE surface without affecting the surface structure. The calculations using the PM3 and DFT methods were performed on single and multiple oligomers to simulate a wide variety of oxidation scenarios. Calculated work function results suggest that regardless of oxidation mechanism, e.g. -OH, =0 or a combination thereof, the experimentally observed levels of surface oxidation are unlikely to lead to a significant change in the electronic structure of PE and that its increased hydrophilic properties are the primary reason for the observed changes in its electrostatic behavior. The calculations for PTFE argue strongly against significant oxidation of that material, as confirmed by the XPS results.

Trigwell, Steve↗
NASA NTRS2007

Electrostatic Properties of Polymers Subjected to Atmospheric Pressure Plasma Treatment; Correlation of Experimental Results with Atomistic Modeling

this study, PE, PTFE, PS and PMMA were exposed to a He+O2, APGD and pre and post treatment surface chemistries were analyzed by X-ray photoelectron spectroscopy and contact angle measurements. Semi-empirical and ab-initio calculations were performed to correlate the experimental results with sonic plausible molecular and electronic structure features of the oxidation process. For the PE and PS, significant surface oxidation showing C-O, C=O, and O-C=O bonding, and a decrease in the surface contact angles was observed. For the PTFE and PM MA, little change in the surface composition was observed. The molecular modeling calculations were performed on single and multiple oligomers and showed regardless of oxidation mechanism, e.g. -OH, =O or a combination thereof, experimentally observed levels of surface oxidation were unlikely to lead to a significant change in the electronic structure of PE and PS, and that the increased hydrophilic properties are the primary reason for the observed changes in its electrostatic behavior. Calculations for PTFE and PMMA argue strongly against significant oxidation of those materials, as confirmed by the XPS results.

Trigwell, S.↗
NASA NTRS2019

Physically-Informed Artificial Neural Networks for Atomistic Modeling of Materials

A new approach is presented for the development of classical interatomic potentials using physically-informed neural networks (PINN) combined with an analytical bond-order atomic interaction model. Due to the strong physical underpinnings, the PINN potentials demonstrate much better transferability than the existing machine-learning potentials while drastically improving the accuracy in comparison with traditional potentials. PINN potentials can be constructed for both metallic and covalent materials in a unified manner. A number of applications of PINN potentials to large-scale molecular dynamics and Monte Carlo simulations and calculation of thermal and mechanical properties of diverse materials are demonstrated. Some of the specific materials systems include silicon and aluminum, as well as alloys and compounds. Computational aspects of PINN potentials are discussed and future developments in this field are outlined.

J Hickman↗