Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics 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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A New Sodium Thioborate Fast Ion Conductor: Na 3 B 5 S 9

Abstract We report a new sodium fast‐ion conductor, Na 3 B 5 S 9 , that exhibits a high Na ion total conductivity of 0.80 mS cm −1 (sintered pellet; cold‐pressed pellet=0.21 mS cm −1 ). The structure consists of corner‐sharing B 10 S 20 supertetrahedral clusters, which create a framework that supports 3D Na ion diffusion channels. The Na ions are well‐distributed in the channels and form a disordered sublattice spanning five Na crystallographic sites. The combination of structural elucidation via single crystal X‐ray diffraction and powder synchrotron X‐ray diffraction at variable temperatures, solid‐state nuclear magnetic resonance spectra and ab initio molecular dynamics simulations reveal high Na‐ion mobility (predicted conductivity: 0.96 mS cm −1 ) and the nature of the 3D diffusion pathways. Notably, the Na ion sublattice orders at low temperatures, resulting in isolated Na polyhedra and thus much lower ionic conductivity. This highlights the importance of a disordered Na ion sublattice—and existence of well‐connected Na ion migration pathways formed via face‐sharing polyhedra—in dictating Na ion diffusion.

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

Atomically Local Electric Field Induced Interface Water Reorientation for Alkaline Hydrogen Evolution Reaction

The slow water dissociation process in alkaline electrolyte severely limits the kinetics of HER. The orientation of H 2 O is well known to affect the dissociation process, but H 2 O orientation is hard to control because of its random distribution. Herein, an atomically asymmetric local electric field was designed by IrRu dizygotic single-atom sites (IrRu DSACs) to tune the H 2 O adsorption configuration and orientation, thus optimizing its dissociation process. The electric field intensity of IrRu DSACs is over 4.00×10 10 N/C. The ab initio molecular dynamics simulations combined with in situ Raman spectroscopy analysis on the adsorption behavior of H 2 O show that the M–H bond length (M=active site) is shortened at the interface due to the strong local electric field gradient and the optimized water orientation promotes the dissociation process of interfacial water. This work provides a new way to explore the role of single atomic sites in alkaline hydrogen evolution reaction.

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

Aqueous Self‐Assembly of Cylindrical and Tapered Bottlebrush Block Copolymers

The self‐assembly of amphiphilic bottlebrush block copolymers (BCPs), featuring backbones densely grafted with two types of side chains, is less well understood compared to linear BCPs. In particular, the solution self‐assembly of tapered bottlebrush BCPs—cone‐shaped BCPs with hydrophilic or hydrophobic tips—remains unexplored. This study investigates eight tapered and four cylindrical bottlebrush BCPs with varied ratios of hydrophobic polystyrene (PS) and hydrophilic poly(acrylic acid) (PAA) side chains, synthesized via sequential addition of macromonomers using ring‐opening metathesis polymerization (SAM‐ROMP). Self‐assembled nanostructures formed in water were analyzed using cryogenic transmission electron microscopy, small‐angle neutron scattering, and dynamic light scattering. Most BCPs generated multiple nanostructures with surface protrusions, including spherical micelles, cylindrical micelles, and vesicles, alongside transitional forms like ellipsoids and semi‐vesicles. Coarse‐grained molecular dynamics simulations supported the experimental findings, which revealed two distinct self‐assembly pathways. The first involved micelle fusion, producing elliptical and cylindrical aggregates, sometimes forming Y‐junctions. The second pathway featured micelle maturation into semivesicles, which developed into vesicles or large compound vesicles. This work provides the first experimental evidence of vesicle formation via semivesicles in bottlebrush BCPs and demonstrates the significant influence of cone directionality on self‐assembly behavior in these cone‐shaped polymeric amphiphiles.

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

Is a Malleable Active Site Loop the Key to High Substrate Promiscuity? Hybrid, Biocatalytic Route to Structurally Diverse Taxoid Side Chains with Remarkable Dual Stereocontrol

Abstract These studies reveal the first structure ofClostridium acetobutylicumalcohol dehydrogenase (CaADH), a protein exhibiting remarkable substrate promiscuity and stereochemical fidelity. The CaADH enzyme is utilized here for synthesizing 20 potential aryl isoserine side chains for the Taxotere family of tubulin‐binding chemotherapeutics. The approach involves dynamic reductive kinetic resolution (DYRKR) upon the corresponding α‐chloro‐β‐keto esters, showing high D‐synstereoselectivity, including those leading to the clinically relevant milataxel (Ar = 2‐furyl) and simotaxel (Ar = 2‐thienyl) side chains. Furthermore, various cross‐coupling chemistries performed on thep‐bromophenyl isoserine side chain significantly enhance the structural diversity of the taxoid side chain library obtained (16 additional taxoid side chains). The CaADH structure is notable: (i) the nicotinamide cofactor is bound in ananti‐conformation, with the amide carbonyl occupying the ketone binding pocket, and (ii) a flexible loop near the active site likely contributes to the remarkable substrate promiscuity observed in CaADH. We present our perspective on the dynamic nature of the CaADH active site through molecular dynamics simulation, proposing a halogen bonding model as a potential mechanism for the remarkable selectivity for an (S)‐configured C─Cl bond, in addition to the D‐facial selectivity, demonstrated across 20 diverse substrates by this remarkable short‐chain dehydrogenase enzyme.

Chemistry↗
DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗
DOE OSTI2022

Superoxide Dismutase 1 Folding Stability as a Target for Molecular Tweezers in SOD1-Related Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation are hallmarks of many severe neurodegenerative diseases including Alzheimer's, Parkinson's and Huntington's disease. As a supramolecular ligand that binds to lysine and arginine residues, the molecular tweezer CLR01 was found to modify the aggregation pathway of disease-relevant proteins in vitro and in vivo with beneficial effects on toxicity. However, the molecular mechanisms of how tweezers exert these effects remain mainly unknown, hampering further drug development. Here, we investigate the modulation mechanism of unfolding and aggregation pathways of SOD1, which are involved in amyotrophic lateral sclerosis (ALS), by CLR01. Using a truncated version of the wildtype SOD1 protein, SOD1 bar , we show that CLR01 acts on the first step of the aggregation pathway, the unfolding of the SOD1 monomer. CLR01 increases, by ~10 °C, the melting temperatures of the A4V and G41D SOD1 mutants, which are commonly observed mutations in familial ALS. Molecular dynamics simulations and binding free energy calculations as well as native mass spectrometry and mutational studies allowed us to identify K61 and K92 as binding sites for the tweezers to mediate the stability increase. The data suggest that the modulation of SOD1 conformational stability is a promising target for future developments of supramolecular ligands against neurodegenerative diseases.

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

Different Enzyme Conformations Induce Different Mechanistic Traits in HIV‐1 Protease

Abstract The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV‐1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol −1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active‐site water molecule is proposed to influence the mechanistic pathway.

Coimbra, João T. S.↗
DOE OSTI2022

Construction of Photo‐Responsive Pd 2 L 4 ‐Type Nanocages based on Feringa's Second‐Generation Motor and Its Guest Binding Ability for C 60

Abstract We present the construction of a M 2 L 4 ‐type metal‐organic nanocage featuring four endohedral Feringa's motor motifs and its adaptive encapsulation towards a C 60 guest molecule. The structure of the cage, though complicated on the 1 H NMR spectrum due to the adoption of mixed ligands, was unambiguously characterized with a combination of ESI‐MS, 2D DOSY, 13 C NMR and particularly the SAXS technique. The molecular motor within the cage demonstrated similar photophysical properties to the uncoordinated one, indicating the motor's function was not compromised when it was anchored in such a confined nanospace. Furthermore, the nanocage showed good guest encapsulation ability towards C 60 , and a guest induced‐fit behavior of the cage was revealed based on the extensive SAXS analysis and molecular dynamics simulation. The adaptive motorized nanocage reported here represents one of the very few examples of integrating individual motors into a discrete nanoconfined system and offers prospects to achieve its non‐equilibrium functions.

Chemistry↗
DOE OSTI2022

Discovery of Complex Binding and Reaction Mechanisms from Ternary Gases in Rare Earth Metal–Organic Frameworks

Understanding the selectivity of metal–organic frameworks (MOFs) to complex acid gas streams will enable their use in industrial applications. In this study, ab initio molecular dynamic simulations (AIMD) were used to simulate ternary gas mixtures (H 2 O-NO 2 -SO 2 ) in rare earth 2,5-dihydroxyterephthalic acid (RE-DOBDC) MOFs. Stronger H 2 O gas-metal binding arose from thermal vibrations in the MOF sterically hindering access of SO 2 and NO 2 molecules to the metal sites. Gas-gas and gas-linker interactions within the MOF framework resulted in the formation of multiple secondary gas species including HONO, HNO 2 , NOSO, and HNO 3 ⁻. Four gas adsorption sites were identified along with a new de-protonation reaction mechanism not observable through experiment. This study not only provides valuable information on competitive gas binding energies in the MOF, but it also provides important chemical insights into transient chemical reactions and mechanisms.

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

Allosteric Tuning of Caspase–7: Establishing the Nexus of Structure and Catalytic Power

Caspase-7 (C7), a cysteine protease involved in apoptosis, is a valuable drug target for its role in human diseases (e. g., Parkinson's, Alzheimer's, sepsis). The C7 allosteric site has great potential for small-molecule targeting, but numerous drug discovery efforts have identified precious few allosteric inhibitors. Here we present the first selective, drug-like inhibitor of C7 along with several other improved inhibitors based on our previous fragment hit. We also provide a rational basis for the impact of allosteric binding on the C7 catalytic cycle by using an integrated approach including X-ray crystallography, stopped-flow kinetics, and molecular dynamics simulations. Our findings suggest allosteric binding disrupts C7 pre-acylation by neutralization of the catalytic dyad, displacement of substrate from the oxyanion hole, and altered dynamics of substrate binding loops. This work advances drug targeting efforts and bolsters our understanding of allosteric structure–activity relationships (ASARs).

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

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗
DOE OSTI2025

Unfolding of the Villin Headpiece Domain: Revealing Structural Heterogeneity with Time‐Resolved X‐Ray Solution Scattering and Markov State Modeling

Understanding protein folding pathways is crucial to deciphering the principles of protein structure and function. Here, the unfolding dynamics of the 35‐residue villin headpiece (HP35) and a norleucine‐substituted variant (2F4K) using a combination of experimental and computational techniques is investigated. Time‐resolved X‐ray solution scattering coupled with equilibrium molecular dynamics simulations and Markov state modeling reveals distinct unfolding mechanisms between the two variants: HP35 and 2F4K. Specifically, HP35 exhibits a two‐state unfolding process, whereas an intermediate state is identified for the 2F4K mutant. A Markov state model constructed from simulations is used to map atomic‐level transitions to experimental observations, providing insights into the role of sequence variations in modulating folding pathways. The findings underscore the importance of integrating experimental and computational approaches to unravel protein unfolding mechanisms between heterogenous structural ensembles.

Nijhawan, Adam K. [Department of Chemistry Northwe↗
DOE OSTI2021

Computational Approach for Rapidly Predicting Temperature–Dependent Polymer Solubilities Using Molecular–Scale Models

One promising approach to recycle multicomponent plastic waste (e. g., multilayer plastic films) is selective dissolution. Selective dissolution is a solvent-mediated process in which differences in polymer solubility in a carefully chosen solvent system are exploited to recover a target polymer. Here, a computational approach was developed that rapidly predicts temperature-dependent polymer solubilities to guide the design of solvent systems for solvent-mediated polymer recycling. Polymer conformations were obtained from molecular dynamics simulations by modeling the polymer as a short oligomer and then used as input to the conductor-like screening model for real solvents (COSMO-RS) for solubility predictions. Using polyethylene (PE) and ethylene vinyl alcohol (EVOH) as representative polymers, the effect of simulation parameters was systematically studied, and predicted solubilities were found to be in good agreement with experimental measurements. Furthermore, the applicability of the approach was demonstrated by identifying selective solvents for PE and EVOH dissolution from a library of 524 solvents.

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

Tailoring Chemical Absorption-Precipitation to Lower the Regeneration Energy of a CO 2 Capture Solvent

Solvent-based CO 2 capture consumes significant amounts of energy for solvent regeneration. To improve energy efficiency, this study investigates CO 2 fixation in a solid form through solvation, followed by ionic self-assembly-aided precipitation. Based on the hypothesis that CO 3 2- ions may bind with monovalent metal ions, we introduced Na + into an aqueous hexane-1,6-diamine solution where CO 2 forms carbamate and bicarbonate. Then, Na + ions in the solvent act as a seed for ionic self-assembly with diamine carbamate to form an intermediate ionic complex. The recurring chemical reactions lead to the formation of an ionic solid from a mixture of organic carbamate/carbonate and inorganic sodium bicarbonate (NaHCO 3 ), which can be easily removed from the aqueous solvent through sedimentation or centrifugation and heated to release the captured CO 2 . Mild-temperature heating of the solids at 80–150 °C causes decomposition of the solid CO 2 -diamine-Na molecular aggregates and discharge of CO 2 . This sorbent regeneration process requires 6.5–8.6 GJ/t CO 2 . It was also found that the organic carbamate/carbonate solid, without NaHCO 3 , contains a significant amount of CO 2 , up to 6.2 mmol CO 2 /g-sorbent, requiring as low as 2.9–5.8 GJ/t CO 2 . In conclusion, molecular dynamic simulations support the hypothesis of using Na + to form relatively less stable, yet sufficiently solid, complexes for the least energy-intensive recovery of diamine solvents compared to bivalent carbonate–forming ions.

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

The homogeneous mixing of MgO and H 2 O at extreme conditions

Investigating water worlds presents a unique opportunity to understand the fundamental processes of planetary formation and evolution. One key aspect is characterizing the interactions between water and rock under the pressures and temperatures present within these worlds. Investigating the conditions for the homogeneous mixing of these materials is imperative to characterizing bulk properties and evolution of water-rich exoplanets. Here, in this study, we use density functional molecular dynamics simulations to study MgO-H 2 O mixtures at high pressure–temperature conditions where H 2 O occurs in solid, superionic, or liquid form. MgO, the representative rocky material, can be either solid or liquid. We start from 500 K at 120 GPa, increasing the temperature step by step up to 8000 K. By inspection, we determine the temperature at which MgO-HO homogeneously mix in our simulations. At 6000 K and 174 GPa is when we find the system to homogeneously mix. This heat-until-it-mixes approach provides us with an upper bound on the temperature for the mixing of MgO and H 2 O. We find that homogeneous mixing occurs at sufficiently low temperatures to be relevant for the collisional growth of a water-rich planet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing the electrical double layer by operando X‐ray photoelectron spectroscopy through a graphene‐carbon nanotube composite window

Abstract The electrical double layer is known to spontaneously form at the electrode‐electrolyte interface, impacting many important chemical and physical processes as well as applications including electrocatalysis, electroorganic synthesis, nanomaterial preparation, energy storage, and even emulsion stabilization. However, it has been challenging to study this fundamental phenomenon at the molecular level because the electrical double layer is deeply “buried” by the bulk electrolyte solution. Here, we report a quantitative probing of the electrical double layer of ionic liquids from the solid side of a photoelectron‐transparent graphene‐carbon nanotube hybrid membrane electrode using X‐ray photoelectron spectroscopy. The membrane window is ultrathin (1‐1.5 nm), large (~1 cm 2 ), and robust, enabling a tight seal of the electrolyte and quantitative measurement with excellent photoelectron signals. By operando monitoring the population changes of cations and anions in response to the applied electrical potentials, we experimentally resolve the chemical structure and dynamics of the electrical double layer, which corroborate results from molecular dynamics simulations. image

Wang, Peng↗
DOE OSTI2022

Manipulating Conjugated Polymer Backbone Dynamics through Controlled Thermal Cleavage of Alkyl Side Chains

The morphological stability of an organic photovoltaic (OPV) device is greatly affected by the dynamics of donors and acceptors occurring near the device's operational temperature. These dynamics can be quantified by the glass transition temperature (T g ) of conjugated polymers (CPs). Because flexible side chains possess much faster dynamics, the cleavage of the alkyl side chains will reduce chain dynamics, leading to a higher T g . In this work, the T g s for CPs are systematically studied with controlled side chain cleavage. Isothermal annealing of polythiophenes featuring thermally cleavable side chains at 140 °C, is found to remove more than 95% of alkyl side chains in 24 h, and raise the backbone T g from 23 to 75 °C. Coarse grain molecular dynamics simulations are used to understand the T g dependence on side chain cleavage. X-ray scattering indicates that the relative degree of crystallization remains constantduring isothermal annealing process. The effective conjugation length is not influenced by thermal cleavage; however, the density of chromophore is doubled after the complete removal of alkyl side chains. Here, the combined effect of enhancing T g and conserving crystalline structures during the thermal cleavage process can provide a pathway to improving the stability of optoelectronic properties in future OPV devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗