Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular 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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

Here, understanding the size- and shape-dependent properties of platinum nanoparticles is critical for enabling the design of nanoparticle-based applications with optimal and potentially tunable functionality. Toward this goal, we evaluated nine different empirical potentials with the purpose of accurately modeling faceted platinum nanoparticles using molecular dynamics simulation. First, the potentials were evaluated by computing bulk and surface properties - surface energy, lattice constant, stiffness constants, and the equation of state - and comparing these to prior experimental measurements and quantum mechanics calculations. Then, the potentials were assessed in terms of the stability of cubic and icosahedral nanoparticles with faces in the {100} and {111} planes, respectively. Although none of the force fields predicts all the evaluated properties with perfect accuracy, one potential - the embedded atom method formalism with a specific parameter set - was identified as best able to model platinum in both bulk and nanoparticle forms.

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

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Assessment of the impact of reactor residence time distribution on non-equilibrium product selectivity of polypropylene pyrolysis using reactive molecular dynamics simulations

Mass residence time distribution (RTD) is considered to be an important factor controlling the product selectivity in the pyrolysis of biomass and plastic wastes along with the pyrolysis chemistry. However, due to the complex pyrolysis chemistry of biomass and plastic waste, the coupling between the reaction chemistry, RTD, and product selectivity is challenging to understand. Here, we introduce a reaction molecular dynamics-based method to examine pyrolysis chemistry and species timescales to assess the impact of RTD on product selectivity and yield. To validate this method, reactive molecular dynamics simulations were conducted for polypropylene pyrolysis and its non-equilibrium product selectivity using 6 different RTDs. We find that the RTD and the reaction chemistry control the peak non-equilibrium product concentrations. The peak monomer (C 3 H 6 ) concentration during pyrolysis can be increased by up to 25 % by using a narrow RTD in the case of polypropylene pyrolysis. We also find that product selectivity is strongly affected by the average residence time and RTD. This coupling between the reaction chemistry, RTD, and product selectivity highlights the need to understand detailed reaction chemistry to control RTD and optimize non-equilibrium product selectivity during polymer and biomass pyrolysis. The present method provides a new way to design RTD for reactors to reach maximized product selectivity of plastic waste and biomass.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Deformation behavior of annealed Cu64Zr36 metallic glass via molecular dynamics simulations

In the present work, we investigate the isothermal annealing process of Cu 64 Zr 36 metallic glass (MG) by means of molecular dynamics (MD) simulations, and characterize the compression properties under different temperatures and compressive rates. The pair distribution functions, the bond-orientational order, and the cluster-type index method are modeled and analyzed to characterize changes in the structure. Results of the modeling and analysis reveal that the MgZn 2 -type Laves, MgCu 2 -type Laves and five-fold local symmetry structures can be formed during the isothermal holding at 950 K. Study on the compression properties shows that at a fixed strain rate, the yield strength and yield drop all increase with decreasing temperature. At a fixed temperature, however, the increase of the strain rate leads to a noticeable increase in the yield strength at 600 K, but has little effect on yield strength at 10 K. Moreover, the modeling and analysis of the structure at a temperature of 10 K and strain-rate of 5 × 10 7 s -1 demonstrate that the order-disorder transformation initiates the shear band.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2025

Unveiling the lithium-ion transport mechanism in Li{sub 2}ZrCl{sub 6} Solid-State Electrolyte {ital via} deep learning-accelerated molecular dynamics simulations.

Lithium zirconium chlorides (LZCs) present a promising class of cost-effective solid electrolytes for next-generation all-solid-state batteries. The unique crystal structure of LZCs plays a crucial role in facilitating lithium-ion mobility, which further affects the electrochemical performance. To understand the underlying mechanism governing ion transport, we employed deep learning-accelerated molecular dynamics simulation on Li2ZrCl6 (trigonal alpha- and monoclinic beta-LZC), focusing specifically on the zirconium coordination environment. Our results reveal that disordered alpha-LZC exhibits the highest ionic conductivity, while beta-LZC demonstrates significantly lower conductivity, closely aligning with experimental findings. The study confirms that across all phases, lithium migration proceeds via the site-to-site hopping mechanism, where variations in site residence times critically impact the overall ionic conductivity. In alpha-LZCs, lithium ions prefer to anisotropically diffuse across interlayers as the result of a lower energy barrier, driven primarily by collective diffusion. In contrast, lithium ions in beta-LZC primarily isotropically diffuse within the intralayer, hindered by higher energy barriers and determined by individual diffusion. The variation in ZrCl6 2- octahedral unit softening, induced by the specific layered arrangement of zirconium atoms, emerges as a critical determinant of the energy barriers across the LZC phases. These atomic-scale insights into the transport processes provide valuable guidance for the rational design and optimization of LZCs-based electrolytes, accelerating their practical application in advanced energy storage technologies.

Guo, Hanzeng↗
DOE OSTI2022

Non-equilibrium molecular dynamics simulation to evaluate the effect of confinement on fluid flow in silica nanopores

The flow of fluids in nano-confinement has applications in separations, water purification, medical systems and in the recovery of fluids in petroleum systems. It is believed that fluids do not obey continuum laws when flowing in nanopores (or in confinement). This paper attempts to shed light on various nano confinement effects such as fluid-wall interactions, pore size and molecular geometry, using molecular dynamic simulation. Here, water, hexane, and methanol flow behaviors were simulated for pore diameters ranging from 1 to 8 nm. In addition to the density analysis of the confined fluids, the fluid flows in saturated nanopores were simulated by the sectional flow method. Water and methanol molecules were completely stabilized in the 1 nm pore. Hexane molecules could not enter the 1 nm pore, due to geometric considerations. For the 2 – 8 nm pores, all fluids showed reduced flow rate compared to the Hagen–Poiseuille flow, due to an interfacial molecular layer stabilized at the pore surface. Flow reduction observed in these studies is contrary to significant flow enhancements observed for fluid flow in carbon nanotubes of similar dimensions. Water flow showed almost constant stick length for all the pore sizes. Methanol flow had the largest stick length. Hexane flow was reduced because of overcrowding of molecules at the pore surface. The effect of confinement diminished with an increase in pore diameter.

42 ENGINEERING↗
DOE OSTI2024

Determination of the Diffusion Coefficients of Binary CH4 and C2H6 in a Supercritical CO2 Environment (500–2000 K and 100–1000 atm) by Molecular Dynamics Simulations

The self-diffusion coefficients of carbonaceous fuels in a supercritical CO2 environment provide transport information that can help us understand the Allam Cycle mechanism at a high pressure of 300 atm. The diffusion coefficients of pure CO2 and binary CO2/CH4 and CO2/C2H6 at high temperatures (500 K~2000 K) and high pressures (100 atm~1000 atm) are determined by molecular dynamics simulations in this study. Increasing the temperature leads to an increase in the diffusion coefficient, and increasing the pressure leads to a decrease in the diffusion coefficients for both methane and ethane. The diffusion coefficient of methane at 300 atm is approximately 0.012 cm2/s at 1000 K and 0.032 cm2/s at 1500 K. The diffusion coefficient of ethane at 300 atm is approximately 0.016 cm2/s at 1000 K and 0.045 cm2/s at 1500 K. The understanding of diffusion coefficients potentially leads to the reduction in fuel consumption and minimization of greenhouse gas emissions in the Allam Cycle.

Energy & Fuels↗
DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗
DOE OSTI2021

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Predicting the Mechanical Response of Polyhydroxyalkanoate Biopolymers Using Molecular Dynamics Simulations

Polyhydroxyalkanoates (PHAs) have emerged as a promising class of biosynthesizable, biocompatible, and biodegradable polymers to replace petroleum-based plastics for addressing the global plastic pollution problem. Although PHAs offer a wide range of chemical diversity, the structure–property relationships in this class of polymers remain poorly established. In particular, the available experimental data on the mechanical properties is scarce. In this contribution, we have used molecular dynamics simulations employing a recently developed forcefield to predict chemical trends in mechanical properties of PHAs. Specifically, we make predictions for Young’s modulus, and yield stress for a wide range of PHAs that exhibit varying lengths of backbone and side chains as well as different side chain functional groups. Deformation simulations were performed at six different strain rates and six different temperatures to elucidate their influence on the mechanical properties. Our results indicate that Young’s modulus and yield stress decrease systematically with increase in the number of carbon atoms in the side chain as well as in the polymer backbone. In addition, we find that the mechanical properties were strongly correlated with the chemical nature of the functional group. The functional groups that enhance the interchain interactions lead to an enhancement in both the Young’s modulus and yield stress. Finally, we applied the developed methodology to study composition-dependence of the mechanical properties for a selected set of binary and ternary copolymers. Overall, our work not only provides insights into rational design rules for tailoring mechanical properties in PHAs, but also opens up avenues for future high throughput atomistic simulation studies geared towards identifying functional PHA polymer candidates for targeted applications.

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

Investigating the time dependence of neutron-proton equilibration using molecular dynamics simulations

Neutron-proton equilibration has previously been studied experimentally in dynamically deformed nuclei in heavy ion collisions [Phys. Rev. Lett. 118, 062501 (2017); Phys. Rev. C 95, 044604 (2017)] using the NIMROD detector array at the Cyclotron Institute at Texas A&M University. Results indicated the composition of the two heaviest fragments originating from the excited projectile-like fragment evolved exponential with respect to their orientation angle to be more similar. Constrained Molecular Dynamics and Antisymmetrized Molecular Dynamics simulations were performed for different formulations of the density dependence of the asymmetry energy term of the nuclear equation of state. The simulations are compared to experimental results, which indicate a better agreement with a softer interaction in the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effect of Tethering Anions in Block Copolymer Electrolytes via Molecular Dynamics Simulations

Ion-containing block copolymers (BCPs) are promising solid-state electrolytes as they can combine ion conduction and mechanical robustness via multiple microphases. However, salt-doped BCPs can have low ion conduction and transference number. One strategy to increase transference number is to tether anions to the polymer backbone (creating a single-ion conductor), though this may impact cation mobility. Leveraging coarse-grained molecular dynamics simulations, we compare the effects of ion concentration and polymer dielectric constant on BCPs with and without bonding anions to the chain (without changing any other chemical properties). We find that single-ion BCPs have higher cation conductivity than analogous salt-doped BCPs, especially at moderate ion concentrations. Furthermore, we also compare several structural features of these systems; notably, tethering anions causes the ions to spread more evenly throughout the conducting phase, which is related to improved cation conduction.

25 ENERGY STORAGE↗
DOE OSTI2021

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Molecular dynamics simulation of effects of solutes on dislocation propagation in Ni-based superalloys

Ni-based superalloys are used in the hot sections of jet turbine engines because of their high strength, stability and resistance to oxidation. The properties of these alloys can be further optimized by adding alloying elements. Therefore, a fundamental knowledge on the effect of different elements on properties of Ni-based superalloys is required. Molecular dynamics simulation could shed light here but its application is hindered by the absence of reliable and computationally cheap semi-empirical potential of the interatomic interaction for 4 and more element alloys. We will present a new Ni-Al-Cr-Nb Finnis-Sinclair (FS) potential specially designed to simulate the dislocation propagation from to  phase. In order to construct this potential, we designed a special algorithm to incorporate the data on element partitioning in the potential development procedure. For example, it is known from experiment, that Cr is mostly present in the gamma phase. Figure 1 shows a snapshot obtained after equilibration of the model of the Ni68Al17Cr15 alloy at T=1000 K using the hybrid Monte-Carlo (MC)/molecular dynamics (MD) simulation with the developed semi-empirical potential. One can clearly see that the Cr partitioning is in agreement with the experimental data. We will discuss the developed algorithm to incorporate the solute partition data in details. Using the developed semi-empirical potential, we first investigated the effect of anti-site defects in the  phase on the single dislocation propagation. It was found that the dislocation velocity increases with the increasing of the anti-site defect concentration. This effect was attributed to smaller number of Al-Al pairs forming during the dislocation migration in the presence of the anti-site defects. Next, we investigated the effect of Nb on the dislocation pair propagation in the Kolbe mechanism. It was found that the addition of Nb leads to considerable decrease in the dislocation propagation rate. This is in agreement with the experimental data on the effect of Nb on the creep resistance of the Ni-based superalloys. We will discuss the origin of this effect.

Mikhail I. Mendelev↗
DOE OSTI2020

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

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

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗