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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2023

One thousand soils for molecular understanding of belowground carbon cycling

While significant progress has been made in understanding global carbon (C) cycling, the mechanisms regulating belowground C fluxes and storage are still uncertain. New molecular technologies have the power to elucidate these processes, yet we have no widespread standardized implementation of molecular techniques. To address this gap, we introduce the Molecular Observation Network (MONet), a decadal vision from the Environmental Molecular Sciences Laboratory (EMSL), to develop a national network for understanding the molecular composition, physical structure, and hydraulic and biological properties of soil and water. These data are essential for advancing the next generation of multiscale Earth systems models. In this paper, we discuss the 1000 Soils Pilot for MONet, including a description of standardized sampling materials and protocols and a use case to highlight the utility of molecular-level and microstructural measurements for assessing the impacts of wildfire on soil. While the 1000 Soils Pilot generated a plethora of data, we focus on assessments of soil organic matter (SOM) chemistry via Fourier-transform ion cyclotron resonance-mass spectrometry and microstructural properties via X-ray computed tomography to highlight the effects of recent fire history in forested ecosystems on belowground C cycling. We observed decreases in soil respiration, microbial biomass, and potential enzyme activity in soils with high frequency burns. Additionally, the nominal oxidation state of carbon in SOM increased with burn frequency in surface soils. This results in a quantifiable shift in the molecular signature of SOM and shows that wildfire may result in oxidation of SOM and structural changes to soil pore networks that persist into deeper soils.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

A predictive model for failure properties of thermoset resins

A predictive model for the three-dimensional failure behavior of engineering polymers has been developed in a recent NASA-sponsored research program. This model acknowledges the underlying molecular deformation mechanisms and thus accounts for the effects of different chemical compositions, crosslink density, functionality of the curing agent, etc., on the complete nonlinear stress-strain response including yield. The material parameters required by the model can be determined from test-tube quantities of a new resin in only a few days. Thus, we can obtain a first-order prediction of the applicability of a new resin for an advanced aerospace application without synthesizing the large quantities of material needed for failure testing. This technology will effect order-of-magnitude reductions in the time and expense required to develop new engineering polymers.

Caruthers, James M.↗
DOE OSTI2025

Dislocation-Grain Boundary Interaction Dataset for FCC Cu

Interactions between dislocations and grain boundaries play a major role in controlling the strength and ductility of structural materials. Experimentally, assessing and probing geometric and stress-based criteria at the local level for dislocation transmission through grain boundaries remains challenging. Therefore, there have been many efforts to systematically generate datasets of dislocation-grain boundary interactions (DGI) via computational models such as molecular dynamics simulations. So far, most DGI datasets have focused only on the subset of nominal minimum-energy grain boundary structures, which limits their applicability, especially to materials processed far from equilibrium. We present a comprehensive database of dislocation-grain boundary interactions for edge, screw, and 60° mixed dislocation with 330 <110> and 257 <112> symmetric tilt grain boundaries (total of 587) in FCC Cu consisting of 73 minimum-energy grain boundary structures and 514 metastable structures. The dataset contains the outcomes for 5234 unique interactions for various dislocation types, grain boundary structures, and applied shear stresses.

36 MATERIALS SCIENCE↗
NASA NTRS1999

The Considere Condition and Rapid Stretching of Linear and Branched Polymer Melts

We analyze the onset of "necking" and subsequent filament failure during the transient uniaxial elongation of viscoelastic fluid samples in extensional rheometers. In the limit of rapid elongation (such that no molecular relaxation occurs), the external work applied is all stored elastically and the Considere criterion originally developed in solid mechanics can be used to quantitatively predict the critical Hencky strain to failure. By comparing the predictions of the Doi-Edwards model for linear homopolymer melts with those of the "Pom-Pom" model for prototypical branched melts we show that the critical strain to failure in rapid elongation of a rubbery material is intimately linked to the molecular topology of the chain, especially the degree of chain branching. The onset of necking instability is monotonically shifted to larger Hencky strains as the number of branches is increased. Numerical computations at finite Deborah numbers also show that there is an optimal range of deformation rates over which homogeneous extensions can be maintained to large strain. We also consider other rapid homogeneous stretching deformations, such as biaxial and planar stretching, and show that the degree of stabilization afforded by inclusion of material with long-chain branching is a sensitive function of the imposed mode of deformation.

McKinley, Gareth H.↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

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

Modeling dark- and light-induced crystal structures and single-crystal optical absorption spectra of ruthenium-based complexes that undergo SO 2 -linkage photoisomerization

A family of coordination complexes of the type [Ru(SO 2 )(NH 3 ) 4 X] m+ Y n – (m, n = 1 or 2) exhibit optical switching capabilities in their single-crystal states. This striking effect is caused by the light-induced formation of SO 2 -linkage photoisomers, which are metastable if kept at suitably cool temperatures. We modeled the dark- and light-induced states of these large crystalline complexes via plane-wave (PW)- and molecular-orbital (MO)-based density functional theory (DFT) and time-dependent DFT in order to calculate their structural and optical properties; the calculated results are compared with experimental data. We show that the PW-DFT-based periodic models replicate the structural properties of these complexes more effectively than the MO-DFT-based molecular-fragment models, observing only small deviations in key bond lengths relative to the experimentally derived crystal structures. The periodic models were also found to more effectively simulate trends seen in experimental optical absorption spectra, with optical absorbance and coverage of the visible region increasing with the formation of the photoinduced geometries. The contribution of the metastable photoisomeric species more heavily focuses on the lower-energy end of the spectra. Spectra generated from the molecular-fragment models are limited by the geometry of the fragment used and the number of excited-state roots considered in those calculations. In general, periodic models outperform the molecular-fragment models owing to their ability to better appreciate the periodic phenomena that are present in these crystalline materials as opposed to MO approaches, which are finite methods. We thus demonstrate that PW-DFT-based periodic models should be considered as a more than viable method for simulating the optical and electronic properties of these single-crystal optical switches.

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

Programmable exploration of magnetic states in Lieb-kagome interpolated lattices

We investigate a hybrid modeling framework in which a quantum annealer is used to simulate magnetic interactions in molecular qubit lattices inspired by experimentally realizable systems. Using phthalocyanine assemblies as a structurally constrained prototype, we model a continuous deformation from a Lieb to a kagome lattice, revealing frustration-driven disorder and magnetic field-induced reordering in the spin structure. Here, the goal is to show how a quantum annealer can operate as a physically instantiated, programmable platform to emulate experimentally relevant lattice deformations and produce observables in a manner analogous to an experimental measurement, enabling the characterization of magnetic arrangements beyond the reach of current molecular architectures. This surrogate modeling approach offers a pathway to explore and iteratively design tunable magnetic states in synthetic materials. The synthetic design, structural characterization, and quantum simulation framework established here defines a modular and scalable paradigm for probing the limits of engineered matter across chemistry, condensed matter, and quantum information science.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Multiple-objective Reinforcement Learning for Inverse Design and Identification

The aim of inverse chemical design is to develop new molecules with given optimized molecular properties or objectives. Recently, generative deep learning (DL) networks are considered as the state-of-the-art in inverse chemical design and have achieved early success in generating molecular structures with desired properties in the pharmaceutical and material chemistry fields. However, satisfying a large number (> 10 objectives) of molecular objectives is a limitation of current generative models. To improve the model’s ability to handle a large number of molecule design objectives, we developed a Reinforcement Learning (RL) based generative framework to optimize chemical molecule generation. Our use of Curriculum Learning (CL) to fine-tune the pre-trained generative network allowed the model to satisfy up to 21 objectives and increase the generative network’s robustness. The experiments show that the proposed multiple-objective RL-based generative model can correctly identify unknown molecules with an 83% to 100% success rate, compared to the baseline approach of 0%. Additionally, this proposed generative model is not limited to just chemistry research challenges; we anticipate that problems that utilize RL with multiple objectives will benefit from this framework.

reinforcement learning, generative model, inverse ↗
DOE OSTI2022

Exploration of the Application of Machine Learning to the improvement of Interatomic Potentials

Current methods for atomistic simulations of material systems suffer from limitations which restrict the ability of the simulation to correctly characterize certain material behavior and physical phenomena. Small scale ab initio (AIMD) modeling is highly accurate but is computationally expensive. Classical molecular dynamics (CMD) simulations use interatomic potentials (IPs) to describe larger systems at a reduced computational cost, but with a reduced accuracy. Creating more robust IPs enable more precise CMD simulations. In this thesis, the application of a specific machine learning process,artifical neural networks (ANNs), to the improvement of IPs and MD simulation is discussed.

O'Connor, Sean↗
DOE OSTI2024

Modeling and Optimization of Zeolites for Contaminant Removal from Coal Combustion Impoundment Leachates

Coal Combustion Residual (CCR) landfills in the U.S. can lead to elevated contaminant concentrations in groundwater and runoff, particularly of arsenic, boron, and selenium. Zeolites can be readily synthesized from materials contained within the coal combustion ash impoundments and can be cation-exchanged to enhance adsorption capacity, selectivity, and reactivity. However, the optimization of zeolites becomes a daunting task when accounting for the variety of Si:Al ratios, the species of extra-framework cations present, and the zeolite pore structure. Molecular simulations provide methods to study and guide the design of zeolites for the sorption of contaminants from aqueous solution. In this work, models that predict the sorption of arsenic, boron, and selenium from water using cation-exchanged zeolites. Because of the lack of experimental adsorption data for these species of contaminants, models were parametrized to reproduce the results of first-principles calculations and then used to predict sorption for zeolites for a dataset containing 6000 combinations of zeolites and sorption conditions. Machine learning was used to train a model to predict sorption for materials in this database based on the results of the molecular simulations. Next, a genetic algorithm was used to optimize zeolites for the removal of each contaminant from aqueous solution for individual impoundment sites based on contaminant concentrations reported by the Electric Power Research Institute.

Findley, John↗
DOE OSTI2021

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

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

What dictates soft clay-like lithium superionic conductor formation from rigid salts mixture

Abstract Soft clay-like Li-superionic conductors, integral to realizing all-solid-state batteries, have been recently synthesized by mixing rigid-salts. Here, through computational and experimental analysis, we clarify how a soft clay-like material can be created from a mixture of rigid-salts. Using molecular dynamics simulations with a deep learning-based interatomic potential energy model, we uncover the microscopic features responsible for soft clay-formation from ionic solid mixtures. We find that salt mixtures capable of forming molecular solid units on anion exchange, along with the slow kinetics of such reactions, are key to soft-clay formation. Molecular solid units serve as sites for shear transformation zones, and their inherent softness enables plasticity at low stress. Extended X-ray absorption fine structure spectroscopy confirms the formation of molecular solid units. A general strategy for creating soft clay-like materials from ionic solid mixtures is formulated.

25 ENERGY STORAGE↗
DOE OSTI2021

Phase evolution in two-phase alloys during severe plastic deformation

Herein, phase evolution in FCC metals with strongly interacting alloy components during severe plastic deformation is investigated using molecular dynamics simulations. Specifically, we study the alloy microstructure in steady state, nucleation and growth of precipitates in supersaturated alloys, and the decomposition of precipitates in undersaturated alloys. The results are compared to a modified effective temperature model, providing a physical understanding for the atomic processes underlying the model and a perspective on its strengths and weaknesses. Key observations in this work are nucleation and growth of precipitates during SPD at a temperature of 100 K; Gibbs-Thomson-like behavior relating solubility to precipitate size under steady-state shearing; a direct relationship between the effective temperature and the shear modulus; and the importance of cluster agglomeration during precipitate growth. The study also reveals that the mechanisms of forced chemical mixing depends on precipitate size, adding complications for effective temperature models describing inhomogeneous systems. The simulations are shown to provide good semiquantitative agreement with experimental findings reported in the literature.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Brown dwarfs and Jovian planets: A comparison

The recent detection of a subluminous companion to the M dwarf star VB8 has renewed interest in the characteristics of objects spanning the mass range from Jupiter to hydrogen burning stars. Atmospheric and interior models were constructed for objects in this mass regime, up to 30 Jupiter masses, with emphasis on understanding the relationship of brown dwarfs such as the VB8 companion to the better-studied Jovian planets. The atmospheric model solves the equation of radiative transfer assuming frequency dependent molecular opacity sources H2, He, H2O, CO, and CH4 which are important by virtue of the high cosmic abundance of their constituent atoms. Condensation of cosmochemically important materials, iron and silicates, in the atmosphere is possible, and the effect of such grains as opacity sources is assessed. The luminosity of the object is presumed due to degenerate cooling following a collapse phase and possibly deuterium burning and an interior model is constructed using as an outer boundary condition the temperature and pressure level at which the atmosphere becomes convective. The interior model is analogous to Jupiter, with a large liquid metallic-hydrogen core and a thinner molecular-hydrogen envelope. The oxidation state of carbon in the outer envelope of a brown dwarf of similar age to Jupiter is a function of the object's mass. This makes the wavelength dependence of the atmospheric opacity sensitive to the carbon to oxygen ratio, since the abundance of the primary source of molecular opacity, H2O, decreases as more oxygen is tied up as CO.

Lunine, J. I.↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗