Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Auxetic Two‐Dimensional Nanostructures from DNA**

Abstract Architectured materials exhibit negative Poisson's ratios and enhanced mechanical properties compared with regular materials. Their auxetic behaviors emerge from periodic cellular structures regardless of the materials used. The majority of such metamaterials are constructed by top‐down approaches and macroscopic with unit cells of microns or larger. There are also molecular auxetics including natural crystals which are not designable. There is a gap from few nanometers to microns, which may be filled by biomolecular self‐assembly. Herein, we demonstrate two‐dimensional auxetic nanostructures using DNA origami. Structural reconfigurations are performed by two‐step DNA reactions and complemented by mechanical deformation studies using molecular dynamics simulations. We find that the auxetic behaviors are mostly defined by geometrical designs, yet the properties of the materials also play an important role. From elasticity theory, we introduce design principles for auxetic DNA metamaterials.

Li, Ruixin↗
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.

Zhou, Laidong↗
DOE OSTI2024

The Influence of Alkyl Spacers and Molecular Weight on the Charge Transport and Storage Properties of Oxy‐Bithiophene‐Based Conjugated Polymers

Abstract Conjugated polymers (CPs) with polar side chains can conduct electronic and ionic charges simultaneously, making them promising for bioelectronics, electrocatalysis and energy storage. Recent work showed that adding alkyl spacers between CP backbones and polar side chains improved electronic charge carrier mobility, reduced swelling and enhanced stability, without compromising ion transport. However, how alkyl spacers impact polymer backbone conformation and, subsequently, electronic properties remain unclear. In this work, we design two oxy‐bithiophene‐based CP series, each featuring progressively extended alkyl spacer lengths and two distinct molecular weight (MW) distributions. Using operando characterisations, we evaluate the (spectro)electrochemical and swelling properties of the polymer thin films, and their performance in organic field‐effect transistors and organic electrochemical transistors. Surprisingly, alkyl spacers negatively impact the hole mobility of our polymers, with higher MW amplifying this effect. Using molecular dynamics simulations, we show that it is thermodynamically favourable for adjacent non‐polar alkyl spacers to aggregate in polar electrolytes, leading to backbone twisting. Further spectroscopic measurements corroborate this prediction. Our findings demonstrate the active interactions between side chain structure, MW and electrolyte/solvent polarity in influencing polymer performance, underscoring the importance of considering solvation environment effects on polymer conformation when designing new mixed conducting CPs for electrochemical applications.

Yu, Hang [Department of Physics and Centre for Pro↗
DOE OSTI2021

Auxetic Two-Dimensional Nanostructures from DNA

Architectured materials exhibit negative Poisson's ratios and enhanced mechanical properties compared with regular materials. Their auxetic behaviors emerge from periodic cellular structures regardless of the materials used. The majority of such metamaterials are constructed by top-down approaches and macroscopic with unit cells of microns or larger. There are also molecular auxetics including natural crystals which are not designable. There is a gap from few nanometers to microns, which may be filled by biomolecular self-assembly. In this work, we demonstrate two-dimensional auxetic nanostructures using DNA origami. Structural reconfigurations are performed by two-step DNA reactions and complemented by mechanical deformation studies using molecular dynamics simulations. We find that the auxetic behaviors are mostly defined by geometrical designs, yet the properties of the materials also play an important role. From elasticity theory, we introduce design principles for auxetic DNA metamaterials.

59 BASIC BIOLOGICAL SCIENCES↗
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 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 OSTI2025

The Role of Asparagine as a Gatekeeper Residue in the Selective Binding of Rare Earth Elements by Lanthanide‐Binding Peptides

Abstract Lanthanide‐binding tag (LBT) peptides selectively complex lanthanide cations (Ln 3+ ) in their binding pockets and are promising for lanthanide separation. However, designing LBTs that selectively target specific Ln 3+ cations remains a challenge due to limited molecular‐level understanding and control of interactions within the lanthanide‐binding pocket. In this study, we reveal that the N5 asparagine residue acts as a gatekeeper in the binding pocket, resulting in a 100‐fold selectivity for smaller Lu 3+ over larger La 3+ cations. Nuclear magnetic resonance spectroscopy and molecular dynamics simulations show that the N5 residue weakly binds to the larger La 3+ cation, permitting H 2 O molecules inside the pocket. For the smaller Lu 3+ cations, the N5 residue forms an inter‐arm hydrogen bond with the E14 glutamic acid residue, locking the Lu 3+ cation in the pocket and preventing H 2 O infiltration. Mutating the N5 asparagine to a D5 aspartic acid prevents such a hydrogen bond, eliminating the gatekeeping mechanism and precipitously reducing selectivity. The resulting binding affinity to Ln 3+ cations is non‐monotonic but generally increases with cation size. These results suggest a molecular design paradigm: the reduced affinity for larger lanthanides is due to open pocket conformations, while the selectivity of smaller Ln 3+ cations over larger ones is due to the gatekeeping hydrogen bond.

Chemistry↗
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↗