Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enabling Efficient Water Splitting with Advanced Materials Designed for High pH Membrane Interface

This project was focused on developing the durable, high-performance materials and interfaces for advanced water splitting, enabling a clear pathway for achieving <$2/Kg H2 (on scale) with efficiency of 43 kWh/kg H 2 via anion exchange membrane (AEM)-based electrolysis. We aimed to advance these final goals via an improved fundamental understanding of both hydrogen and oxygen evolution reactions (HER/OER) leading to novel platinum group metal (PGM)-free catalyst materials in conjunction with critical improvements in membrane and ionomers and gas evolution electrodes with corresponding characterization and testing. Northeastern University (NU) lead this effort focusing on catalyst development and characterization (both in situ and ex situ) while project partners lead improvements in ionomer and membrane materials and will aid in the development of specialized electrode and membrane electrode assemblies. In addition, close collaboration occured with the HydroGEN Energy Materials Network (EMN) National Laboratory consortium including efforts related to use of advanced ionomers, durability protocols and validation of electrolyzer materials (e.g. NREL), multiscale modeling and computation (e.g. LBNL), and molecular dynamics (MD) simulations of the membrane catalyst interface (e.g. SNL). The interactions with HydroGEN included exchange of data and materials as needed to facilitate project success.

08 HYDROGEN↗
DOE OSTI2022

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

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

Development of bubble evolution model for new mechanistic transient fission gas release capability in BISON

This report summarizes efforts within NEAMS to investigate the mechanisms that govern fission gas behavior in UO 2 . In particular, the focus is on understanding how fission gas behavior causes transient fission gas release and fragmentation/pulverization of high burnup structure (HBS) in UO 2 . HBS forms in the periphery of the pellet where temperatures are relatively low. Previously, MD simulations were performed to determine the reaction energies for various Xe and U defects with bubbles, as a function of Xe to vacancy ratio or, equivalently, pressure. As had been shown in FY22, it was found that the unmodified version of the Simple Integrated Fission Gas Release and Swelling (SIFGRS) model within BISON greatly over-predicted the number gas atoms per vacancy in the bubbles in the outer rim of the pellet (a ratio of > 1 million). This was due to slow grain boundary vacancy diffusivity and not accounting for the pressure-dependent reaction energy for Xe interstitials with bubbles. The application of the pressure dependent reaction energies was able to restrict Xe to vacancy ratios to 2:1, which is far more realistic than those originally obtained from SIFGRS.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterize and control the multi-range structure of solution-phase systems with resonant x-ray scattering

This work aims at demonstrating the feasibility and impact of Anomalous X-ray Scattering (AXS) in characterizing the multi-range structures of solution-state systems. We will show preliminary investigation of long-range correlations in concentrated aqueous electrolytes with a combination of AXS and molecular dynamics (MD) simulations. We will also start exploring the capability of AXS to capture intra- and inter-molecular structure of dilute molecular systems, and specifically its sensitivity to the chemical environment surrounding active metal sites. We anticipate that the development of this method will help us to understand ion solvation and transport, which affect the performances of, for instance, electrocatalytic cells, as well as to identify/control the intrinsic structural factors that lead to selectivity, efficiency, and stability control during catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2025

Establishing defect-property relationships for 2D-nanomaterials (Final technical report)

Studies of new families of two-dimensional nanomaterials (2DNMs) have established that they possess unique properties that diverge from those of their bulk counterparts. This project aimed to determine how tolerant 2DNMs are to extreme photon and particle fluxes and to identify the mechanisms governing their radiation response. Early work indicated that graphene is not as representative of other 2DNMs as previously assumed; however, the origin of this difference remained unclear. To address these knowledge gaps, this project investigated the physical processes occurring at multiple length scales in transition metal dichalcogenides (TMDs) and quantified their structural stability and property evolution under far-from-equilibrium conditions. In-situ and ex-situ ion and electron irradiations were performed to directly control and monitor defect formation in TMDs. Complementary density functional theory (DFT) and molecular dynamics (MD) simulations were employed to elucidate the mechanisms of defect generation and evolution. The outcomes of this work established mechanistic understanding of irradiation-induced defect formation in 2DNMs, quantified the radiation tolerance of TMDs, and built a fundamental knowledge base for correlating defect structures with material properties. Collectively, these results provide new insights into the stability of low-dimensional materials in extreme environments and enable the predictive design of radiation-tolerant 2DNMs.

2D materials↗
DOE OSTI2025

Bridging the length scales in ionic separations via data-driving machine learning

We pursued a data science driven machine learning (ML) approach that blended molecular scale attributes informed from molecular dynamics (MD) simulation and materials properties to the selectivity and energy efficiency in targeted ionic separations using electric fields. The model mixtures investigated for ionic separations are pH sensitive and include organic acids, silica and boron, transition metals, such as copper and chromium. There were two major research thrusts of this project. Firstly, we investigated surrogate models and deep learning that relate material chemistries and structures to selective transport of ionic species under applied electric fields. Secondly we investigated how the bipolar junction interfacial design and water dissociation catalyst in bipolar membranes affect reverse bias polarization behavior and pH modulation in deionization platforms as a function of the platform operating parameters (e.g., cell voltage, residence time, and salt feed concentration). As a result of this work, we also were able to start a new direction, namely ML models for molecular design of surfactants.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Role of Surface Hydrophobicity on the Structure and Dynamics of CO2 and CH4 Confined in Silica Nanopores

Advancing a portfolio of technologies that range from the storage of excess renewable natural gas for distributed use to the capture and storage of CO 2 in geological formation are essential for meeting our energy needs while responding to challenges associated with climate change. Delineating the surface interactions and the organization of these gases in nanoporous environments is one of the less explored approaches to ground advances in novel materials for gas storage or predict the fate of stored gases in subsurface environments. To this end, the molecular scale interactions underlying the organization and transport behavior of CO 2 and CH 4 molecules in silica nanopores need to be investigated. To probe the influence of hydrophobic surfaces, a series of classical molecular dynamics (MD) simulations are performed to investigate the structure and dynamics of CO 2 and CH 4 confined in OH-terminated and CH 3 -terminated silica pores with diameters of 2, 4, 6, 8, and 10 nm at 298 K and 10 MPa. Higher adsorption extents of CO 2 compared to CH 4 are noted on OH-terminated and CH 3 -terminated pores. The adsorbed extents increase with the pore diameter. Further, the interfacial CO 2 and CH 4 molecules reside closer to the surface of OH-terminated pores compared to CH 3 -terminated pores. The lower adsorption extents of CH 4 on OH-terminated and CH 3 -terminated pores result in higher diffusion coefficients compared to CO 2 molecules. The diffusivities of both gases in OH-terminated and CH 3 -terminated pores increase systematically with the pore diameter. The higher adsorption extents of CO 2 on OH-terminated and CH 3 -terminated pores are driven by higher van der Waals and electrostatic interactions with the pore surfaces, while CH 4 adsorption is mainly due to van der Waals interactions with the pore walls. These findings provide the interfacial chemical basis underlying the organization and transport behavior of pressurized CO 2 and CH 4 gases in confinement.

Mohammed, Sohaib↗
DOE OSTI2023

Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants

The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.

Guo, Hao↗
DOE OSTI2021

The Role of Hydrophobic Nodes in the Dynamics of Class A β-Lactamases

Class A β-lactamases are known for being able to rapidly gain broad spectrum catalytic efficiency against most β-lactamase inhibitor combinations as a result of elusively minor point mutations. The evolution in class A β-lactamases occurs through optimisation of their dynamic phenotypes at different timescales. At long-timescales, certain conformations are more catalytically permissive than others while at the short timescales, fine-grained optimisation of free energy barriers can improve efficiency in ligand processing by the active site. Free energy barriers, which define all coordinated movements, depend on the flexibility of the secondary structural elements. The most highly conserved residues in class A β-lactamases are hydrophobic nodes that stabilize the core. To assess how the stable hydrophobic core is linked to the structural dynamics of the active site, we carried out adaptively sampled molecular dynamics (MD) simulations in four representative class A β-lactamases (KPC-2, SME-1, TEM-1, and SHV-1). Using Markov State Models (MSM) and unsupervised deep learning, we show that the dynamics of the hydrophobic nodes is used as a metastable relay of kinetic information within the core and is coupled with the catalytically permissive conformation of the active site environment. Finally, our results collectively demonstrate that the class A enzymes described here, share several important dynamic similarities and the hydrophobic nodes comprise of an informative set of dynamic variables in representative class A β-lactamases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling bacterial microcompartment architectures for enhanced cyanobacterial carbon fixation

The carboxysome is a bacterial microcompartment (BMC) which plays a central role in the cyanobacterial CO 2 -concentrating mechanism. These proteinaceous structures consist of an outer protein shell that partitions Rubisco and carbonic anhydrase from the rest of the cytosol, thereby providing a favorable microenvironment that enhances carbon fixation. The modular nature of carboxysomal architectures makes them attractive for a variety of biotechnological applications such as carbon capture and utilization. In silico approaches, such as molecular dynamics (MD) simulations, can support future carboxysome redesign efforts by providing new spatio-temporal insights on their structure and function beyond in vivo experimental limitations. However, specific computational studies on carboxysomes are limited. Fortunately, all BMC (including the carboxysome) are highly structurally conserved which allows for practical inferences to be made between classes. Here, we review simulations on BMC architectures which shed light on (1) permeation events through the shell and (2) assembly pathways. These models predict the biophysical properties surrounding the central pore in BMC-H shell subunits, which in turn dictate the efficiency of substrate diffusion. Meanwhile, simulations on BMC assembly demonstrate that assembly pathway is largely dictated kinetically by cargo interactions while final morphology is dependent on shell factors. Overall, these findings are contextualized within the wider experimental BMC literature and framed within the opportunities for carboxysome redesign for biomanufacturing and enhanced carbon fixation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

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

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young’s modulus, Poisson’s ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young’s modulus and Poisson’s ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a “safe zone” regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

bimodular materials↗