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

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗
DOE OSTI2024

Development of Transformational Solvents for CO2 Separations

The primary objective of this CRADA activity is to use a combined molecular modeling and experimental validation approach to refine and develop transformational solvents for carbon capture. PNNL’s role on this project is currently funded by the Department of Energy’s (DOE) Office of Fossil Energy (OFE). PNNL is developing advanced molecular modeling based on their CO2BOLs solvent platform as a demonstration solvent for the activity; the model was developed and compared against measured data for CO2BOL derivatives. Here, a CRADA with PNNL and GE will leverage their current molecular models and apply them to solvent classes that operate on carbamate chemistry, specifically GE’s aminosilicone solvent class. The molecular models will be used to predict physical and thermodynamic properties, such as viscosity, as a means to predict advanced formulations with reduced viscosity compared to current aminosilicone derivatives, enabling optimized thermodynamic and kinetic metrics for economical carbon capture for this class of materials. Together, PNNL and GE will develop a comprehensive means of linking molecular modeling parameters to intermediate physical properties as a means to improve solvent performance.

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

Structure and ultrafast dynamics of tri-nuclear Ag-/Tl–Pt 2 POP 4 complexes in solution

The energetics and dynamics of ion assembly in solution has broad influence in nanomaterials and inorganic synthesis. To investigate the fundamental processes involved, we present a time-resolved x-ray solution scattering (TR-XSS) study of the trinuclear silver and thallium complexes of the diplatinum ion PtPOP [Pt 2 (H 2 P 2 O 5 )$_4^{4−}$] in aqueous solution. These complexes, their structural properties, and their electronic structure are not well understood and afford a unique opportunity to study the metal–metal bond formation that influences molecular and material assembly in solution. We present model-independent analysis of the observed dynamics as well as an analysis incorporating time-resolved structural refinements of key bond lengths with $<$100 fs time resolution. We find that upon photoexcitation, the Pt atoms contract ∼0.25 Å toward the center of both the Ag- and the Tl-PtPOP complexes, as previously observed for the PtPOP anion. For the AgPtPOP system, an ultrafast Ag-Pt bond expansion of ∼0.2 Å is observed, whereas in contrast, the TlPtPOP system exhibits a Tl-Pt bond contraction of ∼0.3 Å upon photoexcitation. For both complexes, the change in electronic state leads to coherent (“wave-packet”) oscillations along the metal–Pt coordinates. Based on these structural dynamics, we propose an electronic structure model that describes the metal–metal bonding behavior in both the ground and excited state for both complexes.

Lenzen, Philipp [Technical Univ. of Denmark, Lyngb↗
DOE OSTI2020

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

On the Prospect of Chemically Transferable Coarse-Grained Electronic Models for Soft Materials

Electronic coarse-graining (ECG) methods predict quantum-mechanical electronic properties directly from coarse-grained (CG) molecular configurations, enabling electronic predictions at mesoscale length scales. Here, we present a diagnostic assessment of the feasibility of chemically transferable ECG models across a broad polymer-relevant chemical space using all-atom, united-atom, and Martini-scale representations. While high-resolution ECG models achieve near-quantitative accuracy, we show that chemically transferable ECG at the Martini resolution fails because the CG force field does not sample the same configurational distribution of local molecular structure as that underlying the DFT-parameterized ECG model. We demonstrate that our proposed Element-Count-Label (ECL) representation, which augments Martini beads with explicit stoichiometric data, significantly improves chemical generalization across diverse polymer chemistries. However, we find that even with improved chemical resolution, the model cannot recover electronic property distributions that are absent from the configurational space sampled by the CG force field. These results demonstrate that chemically transferable ECG requires future Martini-like force fields to explicitly preserve quantum chemistry–compatible local molecular structure in addition to thermodynamic and structural fidelity.

Kidder, Katherine M [Department of Chemistry; Univ↗
DOE OSTI2020

Artificial Intelligence for Dynamic Phase Transitions in Materials [PowerPoint]

The presentation covers the following areas: emergent phenomena and non-equilibrium transitions in granular materials; materials science applications, including defects, interfaces and emergent phenomena; modeling with molecular dynamics simulations; principal component analysis (PCA); previous works; PCA applied to nonequilibrium transitions; PCA for varied quenched disorder; PCA incorporating particle velocity; and ending with a summary and impacts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Using anion photoelectron spectroscopy of cluster models to gain insights into mechanisms of catalyst-mediated H 2 production from water

We report that metal oxide cluster models of catalyst materials offer a powerful platform for probing the molecular-scale features and interactions that govern catalysis. This perspective gives an overview of studies implementing the combination of anion photoelectron (PE) spectroscopy and density functional theory calculations toward exploring cluster models of metal oxides and metal-oxide supported Pt that catalytically drive the hydrogen evolution reaction (HER) or the water–gas shift reaction. The utility in the combination of these experimental and computational techniques lies in our ability to unambiguously determine electronic and molecular structures, which can then connect to results of reactivity studies. In particular, we focus on the activity of oxygen vacancies modeled by suboxide clusters, the critical mechanistic step of forming proximal metal hydride and hydroxide groups as a prerequisite for H2 production, and the structural features that lead to trapped dihydroxide groups. The pronounced asymmetric oxidation found in heterometallic group 6 oxides and near-neighbor group 5/group 6 results in higher activity toward water, while group 7/group 6 oxides form very specific stoichiometries that suggest facile regeneration. Studies on the trans-periodic combination of cerium oxide and platinum as a model for ceria supported Pt atoms and nanoparticles reveal striking negative charge accumulation by Pt, which, combined with the ionic conductivity of ceria, suggests a mechanism for the exceptionally high activity of this system towards the water–gas shift reaction.

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

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗
DOE OSTI2022

Molecular structure models of amorphous bismuth and cerium carboxylate catalyst precursors

As our societal need for materials and energy has grown, so has our need for catalyst processes in hydrogen production. A major function in these applications, for both homogenous and heterogeneous catalysis processes, is the synthesis of an active metal catalyst. It must first be soluble to control the physical properties of the metal being used. Recent work in metal precursors has begun to turn toward these metal carboxylate types of material. Here, structural models are proposed for bismuth 2-ethylhexanoate and 2,2-dimethyloctanoate and cerium 2-ethylhexanoate. The bismuth compounds have been characterized at different ratios of bismuth to carboxylate as solutions of the free acids. Their structures are most consistent with a Bi 4 (RCO 2 ) 12 motif where the Bi ions are arranged in a flattened tetrahedron with Bi – Bi distances of about 4.3 Å. There is evidence for Bi – O – Bi linkages at low free acid concentrations. The cerium compound is most consistent with a linear tetracerium molecule where the Ce – Ce distances repeat at about 4.3 Å out to 16.4 Å. The models were generated by analogy with known crystal structures and compared to high-energy x-ray scattering data. To further evaluate the models, DFT calculations were made, and the equilibrium geometries were compared. The vibrational spectra calculated from those geometries are presented and compared to the experimental results. Magnetization vs. temperature data was collected on the cerium compound, and its behavior was consistent with the proposed model. A geometrical approach to determining the dimensionality and relative positions of the metal ions in these structures is presented.

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

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

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

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree–Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Chemistry↗
DOE OSTI2023

Planar, curved and twisted molecular nanographenes: Reduction-induced alkali metal coordination

Planar and curved polycyclic aromatic hydrocarbons (PAHs) attract significant attention as molecular models of fullerenes, carbon nanotubes, and graphene, thus stimulating broad investigation of chemical reactivity and materials applications of designed nanocarbon π-systems. Non-planar molecular nanographenes (NGs) recently emerge as advanced anode materials in energy storage, showing high reduction limits and enhanced alkali metal intercalation levels. However, the lack of direct structure–property correlations in such complex hybrid systems impedes their further development and utilization. With a focus on alkali-metal-induced reduction of selected PAHs, we herein review original metal binding and intercalation trends, site specific coordination, and distinct carbon framework responses to stepwise electron uptake. Small planar graphene fragments, like triphenylene and coronene, are compared to π-expanded hexabenzocoronenes, followed by the discussion of bowl-shaped corannulene, sumanene and other carbon bowls, as well as bent, warped, and twisted molecular nanographenes. The effect of size, symmetry, and framework topology along with the structural deformation of carbon backbones upon reduction are analyzed, using recent crystallographic examples of alkali-metal intercalated products. In conclusion, the revealed insights into the structures, binding, and metal intercalation in molecular nanographenes should stimulate their future applications as new functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗