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 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↗
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↗
DOE OSTI2021

27 Al NMR diffusometry of Al 13 Keggin nanoclusters

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, in this work, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes–Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

27Al↗
DOE OSTI2024

Synthesizing realistic sand assemblies with denoising diffusion in latent space

Abstract The shapes and morphological features of grains in sand assemblies have far‐reaching implications in many engineering applications, such as geotechnical engineering, computer animations, petroleum engineering, and concentrated solar power. Yet, our understanding of the influence of grain geometries on macroscopic response is often only qualitative, due to the limited availability of high‐quality 3D grain geometry data. In this paper, we introduce a denoising diffusion algorithm that uses a set of point clouds collected from the surface of individual sand grains to generate grains in the latent space. By employing a point cloud autoencoder, the three‐dimensional point cloud structures of sand grains are first encoded into a lower‐dimensional latent space. A generative denoising diffusion probabilistic model is trained to produce synthetic sand that maximizes the log‐likelihood of the generated samples belonging to the original data distribution measured by a Kullback‐Leibler divergence. Numerical experiments suggest that the proposed method is capable of generating realistic grains with morphology, shapes and sizes consistent with the training data inferred from an F50 sand database. We then use a rigid contact dynamic simulator to pour the synthetic sand in a confined volume to form granular assemblies in a static equilibrium state with targeted distribution properties. To ensure third‐party validation, 50,000 synthetic sand grains and the 1542 real synchrotron microcomputed tomography (SMT) scans of the F50 sand, as well as the granular assemblies composed of synthetic sand grains are made available in an open‐source repository.

Vlassis, Nikolaos N.↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗