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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structural insights into opposing actions of neurosteroids on GABA A receptors

γ-Aminobutyric acid type A (GABA A ) receptors mediate fast inhibitory signaling in the brain and are targets of numerous drugs and endogenous neurosteroids. A subset of neurosteroids are GABA A receptor positive allosteric modulators; one of these, allopregnanolone, is the only drug approved specifically for treating postpartum depression. There is a consensus emerging from structural, physiological and photolabeling studies as to where positive modulators bind, but how they potentiate GABA activation remains unclear. Other neurosteroids are negative modulators of GABA A receptors, but their binding sites remain debated. Here we present structures of a synaptic GABA A receptor bound to allopregnanolone and two inhibitory sulfated neurosteroids. Allopregnanolone binds at the receptor-bilayer interface, in the consensus potentiator site. In contrast, inhibitory neurosteroids bind in the pore. MD simulations and electrophysiology support a mechanism by which allopregnanolone potentiates channel activity and suggest the dominant mechanism for sulfated neurosteroid inhibition is through pore block.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Rejuvenation as the origin of planar defects in the CrCoNi medium entropy alloy

Abstract High or medium- entropy alloys (HEAs/MEAs) are multi-principal element alloys with equal atomic elemental composition, some of which have shown record-breaking mechanical performance. However, the link between short-range order (SRO) and the exceptional mechanical properties of these alloys has remained elusive. The local destruction of SRO by dislocation glide has been predicted to lead to a rejuvenated state with increased entropy and free energy, creating softer zones within the matrix and planar fault boundaries that enhance the ductility, but this has not been verified. Here, we integrate in situ nanomechanical testing with energy-filtered four-dimensional scanning transmission electron microscopy (4D-STEM) and directly observe the rejuvenation during cyclic mechanical loading in single crystal CrCoNi at room temperature. Surprisingly, stacking faults (SFs) and twin boundaries (TBs) are reversible in initial cycles but become irreversible after a thousand cycles, indicating SF energy reduction and rejuvenation. Molecular dynamics (MD) simulation further reveals that the local breakdown of SRO in the MEA triggers these SF reversibility changes. As a result, the deformation features in HEAs/MEAs remain planar and highly localized to the rejuvenated planes, leading to the superior damage tolerance characteristic in this class of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Trapping of spermine, Kukoamine A, and polyamine toxin blockers in GluK2 kainate receptor channels

Abstract Kainate receptors (KARs) are a subtype of ionotropic glutamate receptor (iGluR) channels, a superfamily of ligand-gated ion channels which mediate the majority of excitatory neurotransmission in the central nervous system. KARs modulate neuronal circuits and plasticity during development and are implicated in neurological disorders, including epilepsy, depression, schizophrenia, anxiety, and autism. Calcium-permeable KARs undergo ion channel block, but the therapeutic potential of channel blockers remains underdeveloped, mainly due to limited structural knowledge. Here, we present closed-state structures of GluK2 KAR homotetramers in complex with ion channel blockers NpTx-8, PhTx-74, Kukoamine A, and spermine. We find that blockers reside inside the GluK2 ion channel pore, intracellular to the closed M3 helix bundle-crossing gate, with their hydrophobic heads filling the central cavity and positively charged polyamine tails spanning the selectivity filter. Molecular dynamics (MD) simulations of our structures illuminate interactions responsible for different affinity and binding poses of the blockers. Our structures elucidate the trapping mechanism of KAR channel block and provide a template for designing new blockers that can selectively target calcium-permeable KARs in neuropathologies.

Science & Technology - Other Topics↗
DOE OSTI2025

Non-monotonic size-dependent exciton radiative lifetime in CsPbBr3 nanocrystals

Lead halide perovskite nanocrystals have attracted intense interest due to their desirable optical properties, diverse structural features, and size-tunable excitonic structure. Here we show that, under ambient conditions, a non-monotonic trend in radiative lifetime emerges from the interplay of size, lattice symmetry and excitonic structure. Small nanocrystals exhibit long radiative lifetimes due to weakly emissive excitons, but the oscillator strength increases and shortens the lifetime for nanocrystals approaching intermediate confinement. For larger nanocrystals with higher exciton density of states (DOS), the radiative lifetime is lengthened due to depopulation of the bright exciton manifold into thermally accessible dim states. A size-dependent structural symmetry lowering transition from cubic to orthorhombic is observed by XRD and MD simulations, and the non-monotonic radiative lifetime trend emerges only in lower symmetry structures with an increased dim exciton DOS. These findings shed light on the impact of nanocrystal size and structure on radiative lifetime and pave the way for tailored optical materials in various optical applications.

optical materials↗
DOE OSTI2023

Multidimensional quantitative phenotypic and molecular analysis reveals neomorphic behaviors of p53 missense mutants

Abstract Mutations in the TP53 tumor suppressor gene occur in >80% of the triple-negative or basal-like breast cancer. To test whether neomorphic functions of specific TP53 missense mutations contribute to phenotypic heterogeneity, we characterized phenotypes of non-transformed MCF10A-derived cell lines expressing the ten most common missense mutant p53 proteins and observed a wide spectrum of phenotypic changes in cell survival, resistance to apoptosis and anoikis, cell migration, invasion and 3D mammosphere architecture. The p53 mutants R248W, R273C, R248Q, and Y220C are the most aggressive while G245S and Y234C are the least, which correlates with survival rates of basal-like breast cancer patients. Interestingly, a crucial amino acid difference at one position—R273C vs. R273H—has drastic changes on cellular phenotype. RNA-Seq and ChIP-Seq analyses show distinct DNA binding properties of different p53 mutants, yielding heterogeneous transcriptomics profiles, and MD simulation provided structural basis of differential DNA binding of different p53 mutants. Integrative statistical and machine-learning-based pathway analysis on gene expression profiles with phenotype vectors across the mutant cell lines identifies quantitative association of multiple pathways including the Hippo/YAP/TAZ pathway with phenotypic aggressiveness. Further, comparative analyses of large transcriptomics datasets on breast cancer cell lines and tumors suggest that dysregulation of the Hippo/YAP/TAZ pathway plays a key role in driving the cellular phenotypes towards basal-like in the presence of more aggressive p53 mutants. Overall, our study describes distinct gain-of-function impacts on protein functions, transcriptional profiles, and cellular behaviors of different p53 missense mutants, which contribute to clinical phenotypic heterogeneity of triple-negative breast tumors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Machine-learned interatomic potentials by active learning: amorphous and liquid hafnium dioxide

We propose an active learning scheme for automatically sampling a minimum number of uncorrelated configurations for fitting the Gaussian Approximation Potential (GAP). Our active learning scheme consists of an unsupervised machine learning (ML) scheme coupled with a Bayesian optimization technique that evaluates the GAP model. We apply this scheme to a Hafnium dioxide (HfO 2 ) dataset generated from a "melt-quench" ab initio molecular dynamics (AIMD) protocol. Our results show that the active learning scheme, with no prior knowledge of the dataset, is able to extract a configuration that reaches the required energy fit tolerance. Further, molecular dynamics (MD) simulations performed using this active learned GAP model on 6144 atom systems of amorphous and liquid state elucidate the structural properties of HfO 2 with near ab initio precision and quench rates (i.e., 1.0 K/ps) not accessible via AIMD. The melt and amorphous X-ray structural factors generated from our simulation are in good agreement with experiment. In addition, the calculated diffusion constants are in good agreement with previous ab initio studies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Probing multi-dimensional composition spaces in search of strong metallic alloys

Refractory complex concentrated alloys (RCCA) offer exceptionally high-temperature strength compared to pure metals and dilute alloys, but predictive theory for RCCA design is lacking. We present large-scale molecular Dynamics (MD) simulations of crystal plasticity to explore alloy compositions for maximum mechanical strength, focusing on Fe-Ta-W and Nb-Ta-Mo-W alloy families modeled with Embedded Atom Model (EAM) and Spectral Neighbor Analysis Potentials (SNAP). To efficiently guide the search for strong alloy compositions, we employ iterative optimization using Gaussian process regression. Many simulated RCCA compositions exhibit pronounced cocktail strengthening, with strengths surpassing their strongest constituent metal, tungsten. Contrary to expectations, the highest strength is found on binary edges of the RCCA composition space. Detailed analyses of atomistic simulations reveal that, similar to pure BCC metals, plastic response in RCCA is primarily governed by screw dislocations. However, at large strains, dislocation multiplication and interactions (Taylor hardening) become the dominant mechanisms contributing to RCCA strength.

Materials science↗
DOE OSTI2020

Real-time mass spectrometric characterization of the solid–electrolyte interphase of a lithium-ion battery

Solid-electrolyte interphase (SEI) is a key component that dictates the performance of most electrochemical devices, but the understanding of its chemistry and structure has been limited by the lack of in-situ tools. In this work, we report SEI live-formation in Li-ion battery under operando condition with liquid secondary ion mass spectrometry (SIMS), which, in combination with molecular dynamics (MD) simulation, presents the first dynamic picture with molecular-accuracy about this key component. Before any interphasial chemistry occurs during the initial charging, an electric double layer forms at the electrode-electrolyte interface via the self-assembly of solvent molecules directed by both Li+ and potential of the electrified surface. Such double layer structure predicts the eventual interphasial chemistries formed from the electrolytes, where the negatively-charged electrode surface repels salt anions from the inner-Helmholtz layer and results in an inner SEI that is thin, dense and inorganic in nature but LiF-depleted. It was this dense layer that exercises the major functions of an interphase, i.e., conducting Li+ but insulating electrons. An outer layer that is electrolyte-permeable and organic-rich appears after formation of the inner SEI. In highly-concentrated electrolyte, the presence of LiF in the inner-SEI is raised due to the high population of salt anions in the double layer. These real-time observations at nano-scale reveal in general sense how an interphase inherits chemistry from inner-Helmholtz layer structure, and establish guidelines in engineering desirable interphases for the future battery chemistries.

Zhou, Yufan↗
DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interface microstructure effects on dynamic failure behavior of layered Cu/Ta microstructures

Abstract Structural metallic materials with interfaces of immiscible materials provide opportunities to design and tailor the microstructures for desired mechanical behavior. Metallic microstructures with plasticity contributors of the FCC and BCC phases show significant promise for damage-tolerant applications due to their enhanced strengths and thermal stability. A fundamental understanding of the dynamic failure behavior is needed to design and tailor these microstructures with desired mechanical responses under extreme environments. This study uses molecular dynamics (MD) simulations to characterize plasticity contributors for various interface microstructures and the damage evolution behavior of FCC/BCC laminate microstructures. This study uses six model Cu/Ta interface systems with different orientation relationships that are as- created, and pre-deformed to understand the modifications in the plasticity contributions and the void nucleation/evolution behavior. The results suggest that pre-existing misfit dislocations and loading orientations (perpendicular to and parallel to the interface) affect the activation of primary and secondary slip systems. The dynamic strengths are observed to correlate with the energy of the interfaces, with the strengths being highest for low-energy interfaces and lowest for high-energy interfaces. However, the presence of pre-deformation of these interface microstructures affects not only the dynamic strength of the microstructures but also the correlation with interface energy.

42 ENGINEERING↗
DOE OSTI2022

The allosteric mechanism leading to an open-groove lipid conductive state of the TMEM16F scramblase

TMEM16F is a Ca 2+ -activated phospholipid scramblase in the TMEM16 family of membrane proteins. Unlike other TMEM16s exhibiting a membrane-exposed hydrophilic groove that serves as a translocation pathway for lipids, the experimentally determined structures of TMEM16F shows the groove in a closed conformation even under conditions of maximal scramblase activity. It is currently unknown if/how TMEM16F groove can open for lipid scrambling. Here we describe the analysis of ~400 µs all-atom molecular dynamics (MD) simulations of the TMEM16F revealing an allosteric mechanism leading to an open-groove, lipid scrambling competent state of the protein. The groove opens into a continuous hydrophilic conduit that is highly similar in structure to that seen in other activated scramblases. The allosteric pathway connects this opening to an observed destabilization of the Ca 2+ ion bound at the distal site near the dimer interface, to the dynamics of specific protein regions that produces the open-groove state to scramble phospholipids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Mechanistic insights into a heterobifunctional degrader-induced PTPN2/N1 complex

PTPN2 (protein tyrosine phosphatase non-receptor type 2, or TC-PTP) and PTPN1 are attractive immuno-oncology targets, with the deletion of Ptpn1 and Ptpn2 improving response to immunotherapy in disease models. Targeted protein degradation has emerged as a promising approach to drug challenging targets including phosphatases. We developed potent PTPN2/N1 dual heterobifunctional degraders (Cmpd-1 and Cmpd-2) which facilitate efficient complex assembly with E3 ubiquitin ligase CRL4 CRBN , and mediate potent PTPN2/N1 degradation in cells and mice. To provide mechanistic insights into the cooperative complex formation introduced by degraders, we employed a combination of structural approaches. Our crystal structure reveals how PTPN2 is recognized by the tri-substituted thiophene moiety of the degrader. We further determined a high-resolution structure of DDB1-CRBN/Cmpd-1/PTPN2 using single-particle cryo-electron microscopy (cryo-EM). This structure reveals that the degrader induces proximity between CRBN and PTPN2, albeit the large conformational heterogeneity of this ternary complex. The molecular dynamic (MD)-simulations constructed based on the cryo-EM structure exhibited a large rigid body movement of PTPN2 and illustrated the dynamic interactions between PTPN2 and CRBN. Together, our study demonstrates the development of PTPN2/N1 heterobifunctional degraders with potential applications in cancer immunotherapy. Furthermore, the developed structural workflow could help to understand the dynamic nature of degrader-induced cooperative ternary complexes.

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

Decoupling copolymer, lipid and carbon nanotube interactions in hybrid, biomimetic vesicles

Bilayer vesicles that mimic a real biological cell can be tailored to carry out a specific function by manipulating the molecular composition of the amphiphiles. These bio-inspired and bio-mimetic structures are increasingly being employed for a number of applications from drug delivery to water purification and beyond. Complex hybrid bilayers are the key building blocks for fully synthetic vesicles that can mimic biological cell membranes, which often contain a wide variety of molecular species. While the assembly and morpholgy of pure phospholid bilayer vesicles is well understood, the functionality and structure dramaticlly changes when copolymer and/or carbon nanotube porins (CNTP) are added. The aim of this study is to understand how the collective molecular interactions within hybrid vesicles affect their nanoscale structure and properties. In situ small and wide angle X-ray scattering (SAXS/WAXS) and molecular dynamics simulations (MD) are used to investigate the morphological effect of molecular interactions between polybutadiene polyethylene oxide, lipids and carbon nanotubes (CNT) within the hybrid vesicle bilayer. Within the lipid/copolymer system, the hybrid bilayer morphology transitions from phase separated lipid and compressed copolymer at low copolymer loadings to a mixed bilayer where opposing lipids are mostly separated from the inner region. This transition begins between 60 wt% and 70 wt%, with full homogenization observed by 80 wt% copolymer. The incorporation of CNT into the hybrid vesicles increases the bilayer thickness and enhances the bilayer symmetry. Analysis of the WAXS and MD indicate that the CNT–dioleoyl interactions are much stronger than the CNT–polybutadiene.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Graphene-extracted membrane lipids facilitate the activation of integrin α v β 8

Despite the remarkable electrochemical properties of graphene, strong van der Waals attraction between graphene and biomolecules often causes cytotoxicity, which hinders its applications in the biomedical field. Unfortunately, surface passivation of graphene might stimulate undesired immune response as the nanomaterial triggers cytokine production through membrane receptor activation. In this work, we use all-atom Molecular Dynamics (MD) simulations to unravel the underlying mechanism of graphene-induced inside-out activation of integrin αvβ8, a prominent membrane receptor expressed in immune cells. We model the transmembrane (TM) domains of integrin αvβ8 in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and observe the structural changes in the integrin–membrane complex when interacting with a graphene nanosheet across the membrane. We find that the β8 TM domain interacts with the graphene nanosheet directly or indirectly through extracted lipids, facilitating the pulling of a β8 subunit away from an αv subunit and thus leading to the disruption of the TM domain association by breaking the hydrophobic cluster in the cytoplasmic domains of the αv and β8 subunits. Alanine substitution of two conserved phenylalanine residues on the αv subunit at this hydrophobic cluster further reveals the importance of a stable T-shaped structure in retaining integrin in its inactive state. Our results agree with previous studies on the interactions between other integrin subtypes and their endogenous activators, suggesting an intriguing role that the graphene nanosheet may play in the integrin-related signal transduction during its interaction with the membrane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗