Search NASA⌕ Search

SEARCH · Search NASA

Results for “materials and molecular 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Impact of Confinement and Zwitterionic Ligand Chemistry on Ion–Ion Selectivity of Functionalized Nanopores

Membranes incorporating zwitterionic chemistries have recently emerged as promising candidates for facilitating challenging ion–ion separations. Transport of ions in such membranes predominantly occurs in hydrated nanopores lined with zwitterionic monomers. Here, to shed light on the physics of ion–ion selectivity underlying such materials, we conducted molecular dynamics simulations of sodium halide transport in model nanopores grafted with sulfobetaine methacrylate molecules. Our results reveal that in both functionalized and unfunctionalized nanopores smaller ions prefer to reside near the pore center, while the larger ions tend to reside near the pore walls. An enhancement in the selective transport of larger anions is observed within the unfunctionalized nanopores relative to that in salt-in-water solutions. Upon functionalization of the nanopores with zwitterions (ZIs), the disparities in the anionic distribution profiles within the pores coupled with differences in the anion-ZI interactions result in a slowdown of larger anions relative to smaller anions. Increasing the ZI grafting density exacerbates these effects, further promoting the selective transport of smaller anions. Our results suggest that selectivity toward large anions can be realized by using nanoporous membranes with ZI content that is high enough to facilitate ion/water partitioning into the pores while preserving the characteristic tendency of the unfunctionalized pores to facilitate faster transport of the larger anions. On the other hand, selectivity toward smaller anions can be achieved by targeting ZI content within the pores that is high enough to significantly slow down the transport of large anions but not high enough to hinder the partitioning of ions/water molecules into the pore due to steric effects.

anions↗
DOE OSTI2025

A foundation model for atomistic materials chemistry

Atomistic simulations of matter, especially those that leverage first-principles (ab initio) electronic structure theory, provide a microscopic view of the world, underpinning much of our understanding of chemistry and materials science. Over the last decade or so, machine-learned force fields have transformed atomistic modeling by enabling simulations of ab initio quality over unprecedented time and length scales. However, early machine-learning (ML) force fields have largely been limited by (i) the substantial computational and human effort required to develop and validate potentials for each particular system of interest and (ii) a general lack of transferability from one chemical system to the next. Here, we show that it is possible to create a general-purpose atomistic ML model, trained on a public dataset of moderate size, that is capable of running stable molecular dynamics for a wide range of molecules and materials. We demonstrate the power of the MACE-MP-0 model-and its qualitative and at times quantitative accuracy-on a diverse set of problems in the physical sciences, including properties of solids, liquids, gases, chemical reactions, interfaces, and even the dynamics of a small protein. The model can be applied out of the box as a starting or "foundation" model for any atomistic system of interest and, when desired, can be fine-tuned on just a handful of application-specific data points to reach ab initio accuracy. Establishing that a stable force-field model can cover almost all materials changes atomistic modeling in a fundamental way: experienced users obtain reliable results much faster, and beginners face a lower barrier to entry. Foundation models thus represent a step toward democratizing the revolution in atomic-scale modeling that has been brought about by ML force fields.

Batatia, Ilyes↗
DOE OSTI2024

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inverse design of nanoporous crystalline reticular materials with deep generative models

Reticular frameworks are crystalline porous materials that form via the self-assembly of molecular building blocks in different topologies, with many having desirable properties for gas storage, separation, catalysis, biomedical applications and so on. The notable variety of building blocks makes reticular chemistry both promising and challenging for prospective materials design. Here we propose an automated nanoporous materials discovery platform powered by a supramolecular variational autoencoder for the generative design of reticular materials. We demonstrate the automated design process with a class of metal–organic framework (MOF) structures and the goal of separating carbon dioxide from natural gas or flue gas. Our model shows high fidelity in capturing MOF structural features. Here, we show that the autoencoder has a promising optimization capability when jointly trained with multiple top adsorbent candidates identified for superior gas separation. MOFs discovered here are strongly competitive against some of the best-performing MOFs/zeolites ever reported.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

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

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Bound state dynamics in the long-range spin- 1 2 XXZ model

Experimental platforms based on trapped ions, cold molecules, and Rydberg atoms have made possible the investigation of highly nonlocal spin-1 / 2 Hamiltonians with long-range couplings. Here, we study the effects of such nonlocal couplings in the long-range spin-1 / 2 XXZ Heisenberg Hamiltonian. Here, we calculate explicitly the two-spin energy spectrum, which describes all possible energetic configurations of two spins pointing in a specific direction embedded in a background of spins with opposite orientation. For fast decay of the spin-spin couplings, we find that the two-spin energy spectrum is characterized by well-defined discrete values, corresponding to bound states, separated by a set of continuum states describing the scattering region. In the deep long-range regime instead, the bound states disappear as they get incorporated by the scattering region. The presence of two-spin bound states results to be crucial to determine both two- and many-spin dynamics. On one hand, radically different two-spin spreadings can be observed by tuning the decay of the spin couplings. On the other hand, two-spin bound states enable the dynamical stabilization of effective antiferromagnetic states in the presence of ferromagnetic couplings. Finally, we propose a novel scheme based on a trapped-ion quantum simulator to experimentally realize the long-range XXZ model and to study its out-of-equilibrium properties.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
NASA NTRS1975

AESC multinodal free molecular contamination transport model

The film formation effects are examined of self-induced clouds and condensed films of volatile condensible materials (VCMs) which outgas from such sources as paints, coatings, adhesives, and other nonmetallic surfaces of satellites. A semi-empirical theory for the multiple scattering of free neutral molecules from dielectric surfaces which includes the effects of condensation and evaporation is defined and a mathematical model of the process is developed. The formation rate of thin surface films is considered to be the net accumulation of two simultaneous molecular processes; emission (evaporation) and adsorption (condensation). The net formation rate is considered to be dependent upon the temperature of the surface nodes, the mass and temperatures of the gas molecules, and the total pressure of the incident molecular flux. Apparatus and procedures for measuring these parameters are presented.

Zeiner, E. A.↗
DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2022

A Method for Incorporating Quartz Crystal Microbalance Thermo-Gravimetric Assessment Data in Contamination Modeling

The science objectives of increasingly ambitious space missions require meeting stringent contamination requirements. Contamination transport models are based on established principles such as Fick’s law for contamination diffusion within the source material, Arrhenius factors to account for temperature variations and free molecular transport for spacecraft operating in vacuum. In addition to the science principles, empirical data obtained from ASTM E1559 testing is used for engineering contamination modeling. The ASTM E1559 test provides outgassing rate data for materials at expected spacecraft temperatures during a mission. Multiple quartz crystal microbalances (QCMs) can be operated at different temperatures to measure the amount of contamination condensation that can occur on various spacecraft surfaces. At the conclusion of the ASTM E1559 test a thermo-gravimetric assessment (TGA) can be performed by heating the QCMs at a linear rate to obtain contamination desorption data. This TGA data provides a wealth of information that can be incorporated into contamination models. A method for applying this data is presented.

Anderson, John R↗
NASA NTRS2026

Gateway Element and Payload Materials Outgassing Analyses: HALO, HERMES, and ERSA

Gateway was intended to be humanity’s first space station around the Moon, but its development has been paused as the National Aeronautics and Space Administration (NASA) shifts focus to achieving the United States’ National Space Policy goals. Instead of an orbiting lunar outpost, NASA will now pursue the development of a lunar surface base to support a sustained human presence on the Moon. Before the program’s pause, Gateway’s Induced Environments team worked to ensure payloads and elements (i.e., modules) complied with induced environment requirements. Methods developed and insights gained from this work will have applicability to NASA’s Moon Base and the potential repurposing of Gateway elements and payloads, as well as to induced environments modeling for future space stations. The Gateway program’s induced environment included molecular contamination, electric thruster plume sputter and redeposition, and lunar dust transfer from the Human Landing System (HLS). Primary sources of external molecular contamination included materials outgassing, chemical thruster plume contamination, and vacuum venting. The focus of this paper will be on element- and payload-level materials outgassing analyses performed for Gateway Configuration 1, extending the previously-developed framework for Gateway system-level external molecular contamination modeling. Gateway Configuration 1 consisted of the Power and Propulsion Element (PPE) and the Habitation and Logistics Outpost (HALO). It also included payloads like the European Radiation Sensor Array (ERSA) attached to PPE and the Heliophysics Environmental and Radiation Measurement Experiment Suite (HERMES) attached to HALO. The element- and payload-level analyses to be introduced in this paper for HALO, HERMES, and ERSA enabled high-fidelity descriptions of Gateway’s external molecular contamination environment. Approaches to geometric modeling, meshing, outgassing rate assignment, molecular transport modeling, and analysis methodology will be presented. Element and payload contaminant deposition onto sensitive Gateway receiver surfaces will be summarized and results compared to induced environment requirements. While these results incorporate refinements made over the course of the program, they were not intended to be final. Therefore, modeling assumptions and inputs, potential improvements, and lessons-learned will be documented to inform future work on Moon Base, repurposed elements and payloads, and other space stations.

Gateway↗
DOE OSTI2024

Comparing plant litter molecular diversity assessed from proximate analysis and 13 C NMR spectroscopy

Accurate representation of the chemical diversity of litter in ecosystem-scale models is critical for improving predictions of decomposition rates and stabilization of plant material into soil organic matter. In this contribution, we conducted a systematic review to evaluate how conventional characterization of plant litter quality using proximate analysis compares with molecular-scale characterization using 13 C NMR spectroscopy. Using a molecular mixing model, we converted chemical shift regions from NMR into fractions of carbon (C) in five organic compound classes that are major constituents of plant material: carbohydrates, proteins, lignins, lipids, and carbonylic compounds. We found positive correlations between the acid soluble fraction and carbohydrates, and between the acid insoluble fraction and lignins. However, the acid-soluble fraction underestimated carbohydrates, and the acid insoluble fraction overestimated lignins by 243%. We identified two sources of uncertainties: i) disparities between litter chemical composition based on hydrolysability and actual chemical composition obtained from NMR and ii) conversion factors to translate proximate fractions into organic constituents. Both uncertainties are critical, potentially leading to misinterpretations of decay rates in litter decomposition models. Consequently, we recommend including explicit substrate chemistry data in the next generation of litter decomposition models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces

An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.

36 MATERIALS SCIENCE↗