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.
Computational Modeling of Battery Materials
This chapter provides a broad overview of the various computational modeling techniques used to gain fundamental insights into coupled electrochemical processes that occur in battery materials at electronic-to-mesoscopic scales. Furthermore, representative successes of these techniques in modeling electrodes, electrolytes, and electrode-electrolyte interfaces are highlighted to establish the current state-of-the-art in the field.
Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions
The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.
New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions
Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.
Predicting the Dissolution Rate of Borosilicate Glasses using QSPR analysis based on molecular dynamics simulations
Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition-property relationships of glass and other materials with complex structures. A series of 20 borosilicate and boroaluminosilicate glasses have been modeled by using MD simulations with recently developed effective potentials. Short- and medium-range structures of these glasses were analyzed and, based on these structural information, QSPR analysis of the initial dissolution rate (r0) was made and compared with measured r0 at 90°C and pH 9 using various structural descriptors such as percentage of bridging oxygen species, network connectivity and average ring size. The structural descriptors, Fnet, containing energetic information such as single bond strength and other structural information were also used. It was found that overall network connectivity, average ring size and Fnet give reasonable predictions of the r0 of studied glasses, given the conditions that the glasses are homogeneous and dissolve congruently. Modifying glass compositions to account preferential release of modifiers gives a better prediction for incongruently dissolving glasses. The results were compared with our recent work of predicting glass dissolution behavior from compositions using the topological-constraints-based models.
Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading
Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.
Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators
Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.
Dramatic effect of electrode type on tunnel junction based molecular spintronic devices
A new class of molecular spintronic devices can be fabricated by chemically bonding magnetic molecular channels to the electrodes of a prefabricated tunnel junction with exposed side edges. Prior experimental studies showed that the cyanide-bridged octametallic molecular cluster, [(pzTp)FeIII(CN)3]4[NiII(L)]4¬[O3SCF3]4 [(pzTp) = tetra(pyrazol-1-yl)borate; L = 1-S(acetyl)tris(pyrazolyl)decane] molecule impact depended on the type of metallic electrodes used in the tunnel junction testbed. Experimental magnetization and transport studies showed a dramatic difference in molecule response on tunnel junctions with different combinations of metallic electrodes. Transport via paramagnetic molecular channels on a tunnel junction involving paramagnetic and ferromagnetic metal electrodes was dramatically different than the suppressed current state observed on tunnel junctions involving two ferromagnetic electrodes. We conducted theoretical studies to understand the experimental data and explore a wide range of electrode materials on tunnel junction-based molecular spintronics devices (TJMSD). Here, we report a Monte Carlo simulation study that focuses on understanding the effect of electrodes on the magnetic and physical properties of TJMSD. A 3D Heisenberg model of cross-junction-shaped TJMSD was used for the simulation study. We studied the effects of ferromagnetic, paramagnetic, and antiferromagnetic electrode materials. This study provides insights for designing and understanding futuristic molecular spintronics devices.
Overview of Ablative TPS Modeling at NASA Ames
Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.
Computationally Guided Design of Polymer-Coated Microparticles as Reusable Materials
Long-duration space exploration missions and sustained lunar or Martian surface operations present greater demands for multifunctional and reusable materials. By scaling down the amount of material to be launched from Earth, both mission cost and risk can be reduced. In this regard, leveraging in-space manufacturing capabilities with reusable feedstock materials is an attractive option, as it will allow for articles to be generated on demand, utilized, and then recycled for additional use. NASA’s Enabling Sustained Presence Using Recyclables (ESPUR) project aims to develop reusable materials using polymer-coated microparticles that are bonded via reversible Diels-Alder reactions, where only modest heat is needed to trigger the reverse reaction and enable reuse. For proof-of-concept demonstration, research is currently focused on the fabrication of epoxy microparticles that contain a copoly(carbonate urethane) coating with maleimide and furan functionalities. Here, we discuss the integration of computational materials modeling approaches to help navigate the large design space in this development effort. We perform molecular dynamics (MD) simulations with atomistic and coarse-grained models of the copolymer, which allow us to evaluate the effects of design parameters like the molecular weight and composition on the molecular interactions and chain dynamics. We show how the properties change with the reversible bonds. We also leverage discrete element method (DEM) simulations to assess how microparticle design parameters like the size ratio and volume fraction can be tuned to increase the packing density and number of microparticle contacts to improve the mechanical properties. Our results demonstrate how computational tools can be used in close collaboration with experimental efforts to accelerate material design.
Quantitative Effects of Disorder on Chemically Modified Amorphous Carbon Electrodes
Real materials are disordered. This disorder influences the properties of these materials and the chemical processes that occur at their interfaces. Gaining a molecular-level understanding of the underlying physical manifestations caused by disordered materials is crucial to unraveling and ultimately controlling the efficiency and performance of these materials in a range of energy-related devices. This understanding necessitates measurement techniques through which disorder can be detected, quantified, and monitored. However, such quantitative measurements are notoriously difficult, as effects often average out in ensemble measurements. In this work, we describe how a combination of electrochemical and spatially resolved surface spectroscopy measurements illuminate a molecular-level picture of disorder in materials. Using amorphous carbon as an intrinsically disordered material, we covalently attached a monolayer of ferrocene. Interfacial electron transfer across the amorphous carbon–ferrocene interface is highly sensitive to disruptions of order. Here, by systematically varying linker properties and surface loadings, the influence of lateral interactions between nonuniformly distributed ferrocene headgroups on ensemble electrochemical measurements is demonstrated. Electrochemical and imaging data collectively indicate that conformational flexibility of the ferrocene moieties provides a mechanism to elude repulsive and unbalanced lateral interactions, while rigid linkages provide direct information about the underlying disorder of the material. This study is the first of its kind to quantify and visualize molecular disorder and heterogeneity with an experimental model accessed through ensemble measurements.
A multiscale cohesive law for carbon fiber networks
Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.
Combining Molecular Simulations and Machine Learning to Comprehensively Explore Adsorption Space (Final Report)
This is the final report for the project "Combining Molecular Simulations and Machine Learning to Comprehensively Explore Adsorption Space". The central aim of this project was to develop modeling tools that vastly increase the diversity of molecules for which reliable adsorption isotherms in nanoporous materials are available. Within the project, molecular simulation data for a diverse set of examples was combined with machine-learning (ML) methods to rapidly predict adsorption equilibria in crystalline nanoporous materials (primarily MOFs).
Role of dislocations on martensitic transformation temperatures and microstructure: A molecular dynamics study
Microstructure and defects strongly affect martensitic transformations in metallic alloys. Significant progress has been made in understanding the atomic-level processes that control the role of grain boundaries and precipitates in these solid-to-solid phase transformations. Yet, the role of dislocations and their structures on martensitic transformation temperature and the resulting microstructure remains unclear. Therefore, we used large-scale molecular dynamics simulations to study the forward and reverse transformation of a martensitic material modeled after Ni63Al37 under cyclic thermal loading. The simulations reveal that dislocations in the austenite phase act as one-dimensional seeds for the martensite phase, which is present at temperatures significantly above the martensite start value. We find a reduction in the dislocation density during cyclic thermal loading, which results in the increase in martensite and austenite transition temperatures, in agreement with experiments. Importantly, we extracted a critical martensitic nuclei size for developing stable domains and found that relatively low dislocation densities are needed to grow independent martensitic variants resulting in a multi-domain structure.
Delocalized Excitation Transfer in Open Quantum Systems with Long-Range Interactions
The interplay between coherence and system-environment interactions is at the basis of a wide range of phenomena, from quantum information processing to charge and energy transfer in molecular systems, biomolecules, and photochemical materials. In this work, we use a Frenkel exciton model with long-range interacting qubits coupled to a damped collective bosonic mode to investigate vibrationally assisted transfer processes in donor-acceptor systems featuring internal substructures analogous to light-harvesting complexes. We find that certain delocalized excitonic states maximize the transfer rate and that the entanglement is preserved during the dissipative transfer over a wide range of parameters. We investigate the reduction in transfer caused by static disorder, white noise, and finite temperature and study how transfer efficiency scales as a function of the number of dimerized monomers and the component number of each monomer, finding which excitonic states lead to optimal transfer. Finally, we provide a realistic experimental setting to realize this model in analog trapped-ion quantum simulators. Analog quantum simulation of systems comprising many and increasingly complex monomers could offer valuable insights into the design of light-harvesting materials, particularly in the nonperturbative intermediate parameter regime examined in this study, where classical simulation methods are resource intensive.
Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model
Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.
Exploring a New Regime of Molecular Orbital Physics in 4d Cluster Magnets with Resonant Inelastic X-Ray Scattering
Molecular orbital systems with clusters of heavy transition metal (TM) ions are one of the most important classes of model materials for studying the interplay between local physics and effects of itinerancy. Despite a large number of candidates identified in the family of 4d TM materials, an understanding of their physics from competing microscopic energy scales is still missing. Here, we bridge this gap by reporting the first resonant inelastic x-ray scattering (RIXS) measurement on a well-known series of Ru cluster magnets with a 6H-perovskite structure Ba 3 MRu 2 O 9 (M 3+ = In 3+ , Y 3+ , La 3+ ) comprised of Ru dimers. In addition to providing a microscopic explanation for their anomalous magnetic properties, our RIXS measurements combined with theoretical modeling uncover a new regime of molecular orbital physics where the combined effect of large hopping and small spin-orbit coupling results in highly fragile electronic states in the Ru-dimer compounds directly manifested as an abrupt change in the RIXS spectrum accompanying a tiny change in the local structure tuned by the M-site ion. This unique combination of energy scales found only in the 4d but not the 5d cluster magnets highlights the value of these materials as a new platform for studying quantum phase transition involving molecular orbitals.