Structural Elucidation of a Polypeptoid Chain in a Crystalline Lattice Reveals Key Morphology-Directing Role of the N-Terminus
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.
Molecular basis for lipid recognition by the prostaglandin D 2 receptor CRTH2
Significance Many amphipathic lipids, including cannabinoids, lysophospholipids, prostanoids, and leukotrienes, bind to G protein–coupled receptors (GPCRs) through the lipid bilayer. How these lipid ligands get anchored to their receptors and enter ligand-binding pockets is largely elusive. Our results revealed a distinctive lipid recognition mechanism for the prostaglandin D 2 (PGD 2 ) receptor CRTH2 that is not shared by other prostanoid receptors. Based on the structural analysis, we propose two representative binding modes of lipids in GPCRs. We also identified critical structural motifs and conformational changes of CRTH2 for capturing lipids from the lipidic environment and lipid transition to the binding pocket, thereby offering unprecedented molecular insight into the dynamic lipid-binding processes for GPCRs.
59 BASIC BIOLOGICAL SCIENCES↗
Radiation Damage in Materials—Helium Effects
Despite its scarcity in terrestrial life, helium effects on microstructure evolution and thermo-mechanical properties can have a significant impact on the operation and lifetime of applications, including: advanced structural steels in fast fission reactors, plasma facing and structural materials in fusion devices, spallation neutron target designs, energetic alpha emissions in actinides, helium precipitation in tritium-containing materials, and nuclear waste materials. The small size of a helium atom combined with its near insolubility in almost every solid makes the helium–solid interaction extremely complex over multiple length and time scales. This Special Issue, “Radiation Damage in Materials—Helium Effects”, contains review articles and full-length papers on new irradiation material research activities and novel material ideas using experimental and/or modeling approaches. These studies elucidate the interactions of helium with various extreme environments and tailored nanostructures, as well as their impact on microstructural evolution and material properties.
36 MATERIALS SCIENCE↗
Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel
Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.
36 MATERIALS SCIENCE↗
Fluctuating hydrodynamics and the Rayleigh–Plateau instability
The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.
97 MATHEMATICS AND COMPUTING↗
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↗
Comparative Energetics of Various Membrane Distillation Configurations and Guidelines for Design and Operation
This paper presents a comparative performance study of single-stage desalination processes with major configurations of membrane distillation (MD) modules. MD modules covered in this study are (a) direct contact MD (DCMD), (b) vacuum MD (VMD), (c) sweeping gas MD (SGMD), and (d) air gap MD (AGMD). MD-based desalination processes are simulated with rigorous theoretical MD models supported by molecular thermodynamic property models for the accurate calculation of performance metrics. The performance metrics considered in MD systems are permeate flux and energy efficiency, i.e., gained output ratio (GOR). A general criterion is established to determine the critical length of these four MDs (at fixed width) for the feasible operation of desalination in a wide range of feed salinities. The length of DCMD and VMD is restricted by the feed salinity and permeate flux, respectively, while relatively large AGMD and SGMD are allowed. The sensitivity of GOR flux with respect to permeate conditions is investigated for different MD configurations. AGMD outperforms other configurations in terms of energy efficiency, while VMD reveals the highest permeate production. With larger MD modules, utilization of thermal energy supplied by the hot feed for evaporation is in the order of VMD > AGMD > SGMD > DCMD. Simulation results highlight that energy efficiency of the overall desalination process relies on the efficient recovery of spent for evaporation, suggesting potential improvement in energy efficiency for VMD-based desalination.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation
The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).
25 ENERGY STORAGE↗
Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction
Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Enhanced mobility of dislocation network nodes and its effect on dislocation multiplication and strain hardening
Understanding plastic deformation of crystals in terms of the fundamental physics of dislocations has remained a grand challenge in materials science for decades. To overcome this, the Discrete Dislocation Dynamics (DDD) method has been developed, but its lack of atomistic resolution leaves open the possibility that certain key mechanisms may be overlooked. Here, by comparing large-scale Molecular Dynamics (MD) with DDD simulations performed under identical conditions we uncover significant discrepancies in the predicted strength and microstructure evolution in BCC crystals under high-strain rate conditions. These are traced to unexpected behaviors of dislocation network nodes forming at dislocation intersections, that can move in ways not previously anticipated as revealed by MD. Once these newfound freedoms of nodal motion are incorporated, DDD simulations begin to closely match plastic evolution observed in MD. This additional mechanism of motion whereby non-screw dislocations can change their glide plane profoundly affects fundamental processes of dislocation multiplication, recovery and storage that define strength of metals.
36 MATERIALS SCIENCE↗
Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction
A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.
Kittell, David E. [Sandia National Laboratories (S↗
Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores
Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Machine learning the electric field response of condensed phase systems using perturbed neural network potentials
Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.
Science & Technology - Other Topics↗
Molecular Dynamics Simulations of Helium Bubble Formation in Palladium Tritides [Slides]
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
74 ATOMIC AND MOLECULAR PHYSICS↗
Molecular Dynamics Simulations of Helium Bubble Formation in Palladium Tritides [Slides]
Pd readily absorbs hydrogen and its isotopes, and can be used to purify gas mixtures involving tritium. Tritium decays to He, forming He bubbles. Bubbles causes possible PCT effects swelling, He release, all leading to failures. Radioactive decay experiments take many years. Molecular dynamics (MD) studies can be quickly done. No previous MD methods can simulate He bubble nucleation and growth.
74 ATOMIC AND MOLECULAR PHYSICS↗
MD-HM: Memoization-based Molecular Dynamics Simulations on Big Memory System
Explore the source record for details and available documents.
Computer science↗
Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes
Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.
Duignan, Timothy T.↗
Design of High Efficient Mid-Wavelength Infrared Polarizer on ORMOCHALC Polymer
While an organically modified chalcogenide (ORMOCHALC) can be used to fabricate a polymeric mid-wavelength infrared (MWIR) polarizer with competitive extinction ratio compared to the commercial wire-grid polarizers, which are composed of fragile inorganic materials, there is still a knowledge gap regarding the systematic design process to obtain high transmission efficiency and extinction ratio. As such, a computational parameter study for design optimization is conducted with the geometric parameters of the bilayer grating ORMOCHALC polarizer. The computational study shows that the Fabry–Pérot cavity is the primary mechanism that determines the transmission behaviors and the extinction ratio. A bilayer grating design, guided by the parameter study, is realized through the thermal nanoimprint and metal deposition processes. The extinction ratios measured with the Fourier-transform infrared are 245, 304, and 351 at the wavelength of 3, 4, and 5 μm, respectively. Compared to the state-of-the-art of the polymeric MWIR linear polarizers, the extinction ratio is improved by 1.4 times, and the transmission efficiency is increased by 2.5 times. Theoretical analysis with the multiple-layer model based on the transfer matrix method predicts a matched transmission behavior with the experiment and a full-wave electromagnetic simulation.
42 ENGINEERING↗