Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular dynamics computer 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Substrate binding-induced conformational transitions in the omega-3 fatty acid transporter MFSD2A

Major Facilitator Superfamily Domain containing 2 A (MFSD2A) is a transporter that is highly enriched at the blood-brain and blood-retinal barriers, where it mediates Na + -dependent uptake of ω-3 fatty acids in the form of lysolipids into the brain and eyes, respectively. Despite recent structural insights, it remains unclear how this process is initiated, and driven by Na + . Here, we perform Molecular Dynamics simulations which demonstrate that substrates enter outward facing MFSD2A from the outer leaflet of the membrane via lateral openings between transmembrane helices 5/8 and 2/11. The substrate headgroup enters first and engages in Na + -bridged interactions with a conserved glutamic acid, while the tail is surrounded by hydrophobic residues. This binding mode is consistent with a “trap-and-flip” mechanism and triggers transition to an occluded conformation. Furthermore, using machine learning analysis, we identify key elements that enable these transitions. These results advance our molecular understanding of the MFSD2A transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Completion of initial reduced-order model for flammable refrigerant dispersal in residential spaces

Environmental regulations aimed at reducing global warming impacts of HVAC&R refrigerants have resulted in the phase-out of chlorofluorocarbons (CFCs) in 2010 and hydrochlorofluorocarbons (HCFCs) by 2030 in developed countries with additional restrictions on the use of hydrofluorocarbons (HFCs) set to take effect in 2036. Many of the remaining alternative refrigerants that have lower global warming potential (GWP) and that are suitable for use in HVAC&R systems (e.g., propane, difluoromethane) are flammable to some degree. Flammable refrigerants introduce new challenges and hazards to property and personal safety related to potential deflagration during system maintenance or due to leakage of the refrigerant accumulating in the conditioned space. Standards have been developed to set maximum charge limits for flammable refrigerants in HVAC&R systems; however, existing safety standards and codes still restrict their use. The bodies that maintain and update these codes need publicly available, science-based information to enable credible guidelines for setting safe charge limits for different flammable refrigerants in different HVAC&R applications. In 2016, the Alliance for Responsible Atmospheric Policy, the Air-Conditioning Heating and Refrigeration Institute (AHRI), ASHRAE, the U.S. Department of Energy (DOE), and the State of California began efforts to develop such information. As part of this effort, Oak Ridge National Laboratory (ORNL) began the current, ongoing project to examine imposed charge limits for flammable refrigerants and identify reasonable adjustments to these limits when found appropriate. Past tasks under this project have included development of experimentally vetted, computational fluid dynamics (CFD) simulation approaches to study the results of leakage of flammable refrigerants from various HVAC&R systems into different types of commercial and residential spaces. In this report, we discuss recent efforts to develop a predictive model of the flammable volume fraction and accumulated refrigerant mass in a single-room residential space resulting from the leak of a flammable refrigerant from a small room air conditioning (RAC) unit. The eventual goal is development of a model suitable for public release which could simulate a range of scenarios. In discussions with the AHRTI Flammable Refrigerant Subcommittee (FRS) at the beginning of this effort, a total of 9 input parameters were chosen for consideration including room area, room/door opening area, ventilation fan flow rate, unit/leak height, leak area, leak rate, total refrigerant charge, refrigerant molecular weight, and state of the unit fan.

42 ENGINEERING↗
DOE OSTI2023

Constructing the Mechanism of Dinoflagellate Luciferin Bioluminescence Using Computation

Dinoflagellate luciferin bioluminescence is unique since it does not rely on decarboxylation but is poorly understood compared to that of firefly, bacteria, and coelenterata luciferins. Here, in this work, we computationally investigate possible protonation states, stereoisomers, a chemical mechanism, and the dynamics of the bioluminescence intermediate that is responsible for chemiexcitation. Using semiempirical dynamics, time-dependent density functional theory static calculations, and a correlation diagram, we find that the intermediate’s functional group that is likely responsible for chemiexcitation is a 4-member ring, a dioxetanol, that undergoes [2π + 2π] cycloreversion and the biolumiphore is the cleaved structure. The simulated emission spectra and luciferase-dependent absorbance spectra agree with the experimental data, giving support to our proposed mechanism and biolumiphore. We also compute circular dichroism spectra of the intermediate’s four stereoisomers to guide future experiments in differentiating them.

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

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗
NASA NTRS2012

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗
DOE OSTI2022

Smoothed particle hydrodynamics modeling and analysis of oxide reduction process for uranium oxides

A common kinetic feature for oxide reduction chemical/electrochemical processes is oxygen transport via a porous metallic layer, which has been considered as a rate-determining step for reducing uranium oxides to metallic uranium. Accounting this kinetic behavior must involve the resolution of the moving reactive interface between shrinking oxide and expanding metal phases. This study presents a numerical model using smoothed particles hydrodynamics (SPH) to effectively deal with the evolution of the shrinking core reaction interface and oxygen transport via mass transfer of lithium oxide (Li 2 O) species in multiple mass transfer domains. We successfully validated the proposed model against a theoretical derivation for the oxide reduction process handling a shrinking oxide core with molten salt and metal ash medium on a simple planar geometry. Armed with successful validation results, the model examined a realistic reactant geometry to extend the arguments beyond the one-dimensional analyses, allowing the proposed model to apply to general application scenarios with multi-dimensional geometry. Here, this study demonstrated that the proposed model could simulate and evaluate an arbitrary reaction basket design without iterative experimental trials, which is prohibitive for a scaled high-temperature molten salt study in an inert environment. The construct potentially provides not only deep insights on multiphysics behaviors governing the process dynamics but also a robust framework for evaluating and screening candidate basket designs in the most cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Utilizing CellBox to Describe Varied Dynamical Systems

CellBox is a hybrid modeling framework that integrates machine-learning approaches with explicit mathematical modeling to determine the behavior of dynamical systems, intended for determining the complex interactions of molecule and phenotype within a cellular environment using perturbation/response-based data. The original work utilized analysis of a melanoma cell line, SK-Mel-133, which we sought to expand to and evaluate effectiveness in describing other nonlinear systems both within and outside of biology. We evaluated with two models. One model is based upon the free-fall of an object under a velocity-dependent drag force; this produced inherent perturbation/response data. The other model is based upon a form of mass-action molecular kinetics; for this, we ran the simulation to steady state in each species for each perturbation and select the steady state as input for CellBox. We determined that CellBox, outside of its original use case, can still strongly predict even non-biological dynamical systems in the small-perturbation regime, with r = 0.99977 to 0.85123, p < 0.05 for freefall. However, CellBox cannot accurately predict behavior within large-perturbation regimes, p >> 0.05 for freefall. In addition, we performed a naïve comparison of CellBox network structure to known network structure in the mass-action model to assess accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Relating Geometric Nanoconfinement and Local Molecular Environment to Diffusion in Ionic Polymer Membranes

Understanding structure-diffusion relationships in ionic polymer membranes not only illuminates fundamental mechanisms for membrane transport, but also informs further membrane development. In this study, we employ NMR diffusometry and molecular dynamics (MD) simulations to investigate diffusion coefficients and activation energies of diffusion (Ea) for both a lithium-form perfluorosulfonic acid membrane and aqueous solutions of lithium triflate. NMR shows that, at high water/ion molar ratio (?), Ea for water diffusion is slightly lower in the membrane than in solution. Conversely, membrane transport exhibits much higher Ea at low hydration as compared to solutions. MD simulations of a model system consisting of carbon nanotubes with varying diameter reveal that Ea of diffusion clearly relates to both the geometric nanoconfinement of the hydrophilic pathways in the membrane as well as the local molecular environment. These results demonstrate that Ea of diffusion can serve as a revelatory tool for the study of molecular transport processes and the coupling of morphology with transport at the nanoscale. This study thus provides new insights and new experimental and computational models for understanding transport in a wide range of polymer membranes, such as those used in molecular separations applications.

Zhang, Rui↗
DOE OSTI2021

Mechanism of Rad26-assisted rescue of stalled RNA polymerase II in transcription-coupled repair

Transcription-coupled repair is essential for the removal of DNA lesions from the transcribed genome. The pathway is initiated by CSB protein binding to stalled RNA polymerase II. Mutations impairing CSB function cause severe genetic disease. Yet, the ATP-dependent mechanism by which CSB powers RNA polymerase to bypass certain lesions while triggering excision of others is incompletely understood. Here we build structural models of RNA polymerase II bound to the yeast CSB ortholog Rad26 in nucleotide-free and bound states. This enables simulations and graph-theoretical analyses to define partitioning of this complex into dynamic communities and delineate how its structural elements function together to remodel DNA. We identify an allosteric pathway coupling motions of the Rad26 ATPase modules to changes in RNA polymerase and DNA to unveil a structural mechanism for CSB-assisted progression past less bulky lesions. Our models allow functional interpretation of the effects of Cockayne syndrome disease mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Hub stability in the calcium calmodulin-dependent protein kinase II

The calcium calmodulin protein kinase II (CaMKII) is a multi-subunit ring assembly with a central hub formed by the association domains. There is evidence for hub polymorphism between and within CaMKII isoforms, but the link between polymorphism and subunit exchange has not been resolved. Here, we present near-atomic resolution cryogenic electron microscopy (cryo-EM) structures revealing that hubs from the α and β isoforms, either standalone or within an β holoenzyme, coexist as 12 and 14 subunit assemblies. Single-molecule fluorescence microscopy of Venus-tagged holoenzymes detects intermediate assemblies and progressive dimer loss due to intrinsic holoenzyme lability, and holoenzyme disassembly into dimers upon mutagenesis of a conserved inter-domain contact. Molecular dynamics (MD) simulations show the flexibility of 4-subunit precursors, extracted in-silico from the β hub polymorphs, encompassing the curvature of both polymorphs. The MD explains how an open hub structure also obtained from the β holoenzyme sample could be created by dimer loss and analysis of its cryo-EM dataset reveals how the gap could open further. An assembly model, considering dimer concentration dependence and strain differences between polymorphs, proposes a mechanism for intrinsic hub lability to fine-tune the stoichiometry of αβ heterooligomers for their dynamic localization within synapses in neurons.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Mechanics↗
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 OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2024

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗
NASA NTRS1997

The rate of the reaction between C2H and C2H2 at interstellar temperatures

The reaction between the radical C2H and the stable hydrocarbon C2H2 is one of the simplest neutral-neutral hydrocarbon reactions in chemical models of dense interstellar clouds and carbon-rich circumstellar shells. Although known to be rapid at temperatures > or = 300 K, the reaction has yet to be studied at lower temperatures. We present here ab initio calculations of the potential surface for this reaction and dynamical calculations to determine its rate at low temperature. Despite a small potential barrier in the exit channel, the calculated rate is large, showing that this reaction and, most probably, more complex analogs contribute to the formation of complex organic molecules in low-temperature sources.

Non-NASA Center↗
DOE OSTI2020

Direct observation of 3D atomic packing in monatomic amorphous materials

Liquids and solids are two fundamental states of matter. Although the structure of crystalline solids has long been solved by crystallography, our understanding of the 3D atomic structure of liquids and amorphous materials remained speculative. In 1952, Frank hypothesized that icosahedral order is the prevalent atomic motif in monatomic liquids. Over the past six decades, there have been a great deal of experimental, computational and theoretical studies to understand the structure of liquids and amorphous materials. A polytetrahedral packing model was proposed to explain the 3D atomic structure of monatomic liquids and amorphous materials, in which icosahedral order is a key feature. The icosahedral order has also been found to play an important role in the structure of metallic glasses and quasicrystals. Here we advance atomic electron tomography to determine for the first time the 3D atomic positions in monatomic amorphous materials, including a Ta thin film and two Pd nanoparticles. Despite the different chemical composition and synthesis methods, we observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Contrary to the traditional understanding, the pentagonal bipyramids do not assemble full icosahedra, but are closely connected with each other to form networks that extend to medium-range scale. Our molecular dynamic simulations, coupled with experimental results, further reveal that the 3D atomic structure of monatomic liquids is similar to the experimental amorphous material. Moreover, we find that the pentagonal bipyramid networks, prevalent in monatomic liquids, rapidly grow in size and form a fraction of full icosahedra during the quench from a liquid to a metallic glass state, providing a possible new picture of structural evolution through the glass transition.

Yuan, Yakun↗
DOE OSTI2026

Revealing the Hidden Third Dimension of Point Defects in Two-Dimensional MXenes

Point defects govern many important functional properties of two-dimensional (2D) materials. However, resolving the three-dimensional (3D) arrangement of these defects in multi-layer 2D materials remains a fundamental challenge, hindering rational defect engineering. Here, we overcome this limitation using an artificial intelligence-guided electron microscopy workflow to map the 3D topology and clustering of atomic vacancies in Ti3C2TX MXene. Our approach reconstructs the 3D coordinates of vacancies across hundreds of thousands of lattice sites, generating robust statistical insight into their distribution that can be correlated with specific synthesis pathways. This large-scale data enables us to classify a hierarchy of defect structures-from isolated vacancies to nanopores-revealing their preferred formation and interaction mechanisms, as corroborated by molecular dynamics simulations. This work provides a generalizable framework for understanding and ultimately controlling point defects across large volumes, paving the way for the rational design of defect-engineered functional 2D materials.

2D materials↗