Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conformation dynamics”

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 739 records · Page 41

ToPolyAgent: AI agents for coarse-grained bead-spring topological polymer simulations

We introduce ToPolyAgent, a multi-agent AI framework for performing coarse-grained molecular dynamics (MD) simulations of topological polymers through natural language instructions. By integrating large language models (LLMs) with domain-specific computational tools, ToPolyAgent supports both interactive and autonomous simulation workflows across diverse polymer architectures, including linear, ring, brush, and star polymers, as well as dendrimers. The system consists of four LLM-powered agents: a Config Agent for generating initial polymer–solvent configurations, a Simulation Agent for executing LAMMPS-based MD simulations and conformational analyses, a Report Agent for compiling markdown reports, and a Workflow Agent for streamlined autonomous operations. Interactive mode incorporates user feedback loops for iterative refinements, while autonomous mode enables end-to-end task execution from detailed prompts. We demonstrate ToPolyAgent's versatility through case studies involving diverse polymer architectures under varying solvent conditions, thermostats, and simulation lengths. Furthermore, we highlight its potential as a research assistant by directing it to investigate the effect of interaction parameters on the linear polymer conformation, and the influence of grafting density on the persistence length of the brush polymer. By coupling natural language interfaces with rigorous simulation tools, ToPolyAgent lowers barriers to complex computational workflows and advances AI-driven materials discovery in polymer science. It lays the foundation for autonomous and extensible multi-agent scientific research ecosystems.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗

Structure-Derived Proton-Transfer Mechanism of Action Human Pyruvate Dehydrogenase

The derivative of vitamin B1 thiamin pyrophosphate (TPP) is a cofactor of pyruvate dehydrogenase (E1p) that is involved in decarboxylation of pyruvate followed by reductive acetylation of lipoic acid covalently bound to a lysine residue of dihydrolipoamide acetyltransferase. The structure of E1p recently determined in our laboratory revealed patterns of association of foul subunits and specifics of two TPP binding sites. The mechanism of action in part includes a conserved hydrogen bond between the N1' atom of the aminopyrimidine ring of the cofactor and the carboxylate group of Glu59 from the beta subunits, and a V-conformation of the cofactor that brings the N4' atom of the aminopyrimidine ring to the distance of the intramolecular hydrogen bond formed with the C2-atom of the thiazolium moiety. The carboxylate group of Glu59 is the local proton acceptor that enables proton translocation within the aminopyrimidine ring and stabilization of the rare N4' - iminopyrimidine tautomer. Based on the analysis of E1p structure, we postulate that the protein environment drives N4' - amino/N4' - imino dynamics resulting in a concerted shuttle-like movement of the subunits. We also propose that this movement of the subunits is strictly coordinated with the two enzymatic reactions carried out in E1p by each of the two cofactor sites. It is proposed that these reactions are in alternating phases such that when one active site is involved in decarboxylation, the other is involved in acetylation of lipoyl noiety.

Ciszak, Ewa↗

The Acetylene Bridge in Intramolecular Singlet Fission: A Boon or A Nuisance?

Various analogues of the alkylsilylacetylene group are frequently used as auxiliary groups to enhance the solubility and stability of the acene dimer core, widely used as platforms to investigate intramolecular singlet fission (iSF) mechanisms. However, while in the 2,2’-linked dimers they are primarily auxiliary groups, these are essential fragments of the bridging units in 6,6’/5,5’-linked dimers, the two preferred choices for dimerization. The starkly different iSF dynamics observed in the two variants raise the question of what role the acetylene bridges play. Here, we systematically designed a set of (oligo-)para-phenylene bridged 2,2’-linked pentacene dimers with an additional acetylene fragment in the bridging unit to mimic the structure of 6,6’-linked dimers. Contrasting the results with previously reported analogous 2,2’-linked and 6,6’-linked pentacene dimers reveals that the acetylene bridges contribute to significant conformational freedom. This effect provides a mechanism to promote spin evolution within the triplet pair to achieve free triplets but also offers new parasitic pathways for triplet-pair recombination, revealing that this structural motif can be both a boon and a nuisance. Additionally, our analysis reveals that these bridges directly modify the electronic states, highlighting significant pitfalls of the standard chromophore-bridge-chromophore framework used to design and interpret photophysics of iSF materials.

14 SOLAR ENERGY↗

Structure and functional determinants of Rad6–Bre1 subunits in the histone H2B ubiquitin-conjugating complex

The conserved complex of the Rad6 E2 ubiquitin-conjugating enzyme and the Bre1 E3 ubiquitin ligase catalyzes histone H2B monoubiquitination (H2Bub1), which regulates chromatin dynamics during transcription and other nuclear processes. Here, we report a crystal structure of Rad6 and the non-RING domain N-terminal region of Bre1, which shows an asymmetric homodimer of Bre1 contacting a conserved loop on the Rad6 ‘backside’. This contact is distant from the Rad6 catalytic site and is the location of mutations that impair telomeric silencing in yeast. Mutational analyses validated the importance of this contact for the Rad6–Bre1 interaction, chromatin-binding dynamics, H2Bub1 formation and gene expression. Moreover, the non-RING N-terminal region of Bre1 is sufficient to confer nucleosome binding ability to Rad6 in vitro. Interestingly, Rad6 P43L protein, an interaction interface mutant and equivalent to a cancer mutation in the human homolog, bound Bre1 5-fold more tightly than native Rad6 in vitro, but showed reduced chromatin association of Bre1 and reduced levels of H2Bub1 in vivo. These surprising observations imply conformational transitions of the Rad6–Bre1 complex during its chromatin-associated functional cycle, and reveal the differential effects of specific disease-relevant mutations on the chromatin-bound and unbound states. Overall, our study provides structural insights into Rad6–Bre1 interaction through a novel interface that is important for their biochemical and biological responses.

59 BASIC BIOLOGICAL SCIENCES↗

Lower length scale informed improvements to Bison U-Pu-Zr fuel swelling model

Due to renewed interest in metallic fuels from the U-(Pu)-Zr material system, significant improvements have recently been made to BISON's capability to simulate metallic fuel, particularly in the area of swelling. In this report, lower length scale simulations that have been conducted to inform the engineering-scale models in BISON during FY20 are described. For the high-temperature regions of the fuel that are dominated by the $\gamma$ phase of U-(Pu)-Zr, two new models for swelling have implemented in BISON. Phase-field simulations were conducted to calculate parameters for these and other BISON models that control when gaseous swelling ceases and fission gas release begins. For the low-temperature regions of the fuel dominated by the phase of uranium, microstructurally resolved simulations of pore growth and crystal deformation were performed to understand the effect of irradiation-generated crystal shape changes. Conformally meshed microstructures with pores were simulated with a crystal plasticity model and with a burnup eigenstrain model. Plastic deformation will not cause pore size to increase; burnup eigenstrain (irradiation-induced crystal shape changes) will increase the volume of existing porosity and cause it to become anisotropic in shape. Additional model development work for early-stage fission gas cluster/bubble formation in the $\gamma$ phase were performed, by coupling with defect clustering within cluster dynamics framework. Sensitivity analysis were conducted to identify the critical thermo-kinetic parameters needed for rigorous assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Damage to Polystyrene Polymer Film by Shock Wave Induced Bubble Collapse

Metallic surfaces that are in contact with solutions are commonly used in numerous applications where these surfaces can be damaged by shock wave induced bubble collapse. Use of polymer films that coat such surfaces to prevent them from damage requires a better understanding of how much harm collapsing bubbles produce in the films. In this study, we report the results from coarse-grained molecular dynamics simulations to study the damage to polystyrene (PS) films coating a hard surface. The damage was caused by a collapsing nanobubble located in the proximity of the film and interacting with an impinging shock wave. This collapse produces a high-speed water jet that impacts the PS film with a greater force than the shock front and creates cavities/pits in the PS film. We observed that polymer molecules located in the jet vicinity undergo conformational extension in the direction perpendicular to the jet motion, while chain molecules in the rest of the film undergo compression. Furthermore, we also observed that damage to the film is sensitive to the strength of the shock wave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computation of separated flows by a vortex-tracing algorithm

Numerical solutions for two-dimensional, time-dependent, separated flows around bodies are obtained, using a new version of the vortex method. This method provides an efficient representation of flows involving large regions of separation. The modifications incorporated in the new version, which improve its accuracy, versatility, and computing speed, are described. The computer cost is only of the order of the 3/2 power on N, instead of N-squared, for each step with N vortices. Arbitrary shapes can be treated; a conformal mapping is not required. Special attention is paid to the viscous character of the solution and to the accurate computation of the pressure distribution at the body surface. The vortex solution for the outer flow is coupled to an inner solution for the attached part of the boundary layer. Numerical results are presented for several bluff bodies exhibiting dependence on Reynolds number, for stationary airfoils under steady or transient conditions and for oscillating airfoils, including dynamic stalls. These results are compared with other available results, analytical or experimental, and demonstrate the enhanced reliability and accuracy of the improved method.

Spalart, P. R.↗

Building bulk geometry from the tensor Radon transform

Using the tensor Radon transform and related numerical methods, we study how bulk geometries can be explicitly reconstructed from boundary entanglement entropies in the specific case of AdS 3 /CFT 2 . We find that, given the boundary entanglement entropies of a 2d CFT, this framework provides a quantitative measure that detects whether the bulk dual is geometric in the perturbative (near AdS) limit. In the case where a well-defined bulk geometry exists, we explicitly reconstruct the unique bulk metric tensor once a gauge choice is made. We then examine the emergent bulk geometries for static and dynamical scenarios in holography and in many-body systems. Apart from the physics results, our work demonstrates that numerical methods are feasible and effective in the study of bulk reconstruction in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Structural Constraints on Excited-State Properties in Dimeric Cu(I) Diimine Complexes

Copper(I) bis-diimine complexes have played important roles in light-activated processes that can lead to their potential applications in photocatalysis and chemical sensing. Their metal-to-ligand charge-transfer (MLCT) excited-state properties are tunable by various structural factors. Dimeric Cu(I) complexes with connecting diimine derivative ligands offer another structural tuning platform for the excited-state properties. Here, we investigate excited-state properties in two covalently connected dimeric Cu(I)'s with varying structural constraints exerted by the number of carbons in the polyethylene bridge (C0 and C4) connecting the two copper(I) diimine moieties. An interesting feature of Cu(I) diimine complexes is their ability to flatten following a photoinduced structural change. Herein, we observe larger structural constraints and more structural rearrangement required upon excitation of the longer bridged complex C4 to achieve a conformation toward a more flattened tetrahedral coordination geometry compared to the shorter bridged C0. Vibrational wavepacket analysis of these complexes further supports the effect of these structural constraints where we observe a more rapid dephasing of the C0 complex, as opposed to the C4 complex, despite similar normal mode vibrations. The experimental results were supplemented by TDDFT calculations. In conclusion, the studies provide insight into using metal-metal interactions through constraints to tune excited-state dynamics and pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗

Interactions of anesthetics with the membrane-water interface

Although the potency of conventional anesthetics correlates with lipophilicity, an affinity to water also is essential. It was recently found that compounds with very low affinities to water do not produce anesthesia regardless of their lipophilicity. This finding implies that clinical anesthesia might arise because of interactions at molecular sites near the interface of neuronal membranes with the aqueous environment and, therefore, might require increased concentrations of anesthetic molecules at membrane interfaces. As an initial test of this hypothesis, we calculated in molecular dynamics simulations the free energy profiles for the transfer of anesthetic 1,1,2-trifluoroethane and nonanesthetic perfluoroethane across water-membrane and water-hexane interfaces. Consistent with the hypothesis, it was found that trifluoroethane, but not perfluoroethane, exhibits a free energy minimum and, therefore, increased concentrations at both interfaces. The transfer of trifluoroethane from water to the nonpolar hexane or interior of the membrane is accompanied by a considerable, solvent-induced shift in the conformational equilibrium around the C-C bond.

NASA Center ARC↗

Controlling Structural, Electronic, and Energy Flow Dynamics of Catalytic Processes through Tailored Nanostructures

MoS 2 (molybdenum disulfide) is a highly-versatile catalyst material for support of numerous reactions from hydrodesulfurization and denitrogenation to the focus of this renewal proposal: hydrogenation of CO/CO 2 towards (higher) alcohols. At the same time, MoS 2 is a non-toxic, environmentally-benign and rather inert material- which under ambient conditions for some time has served as a lubricant and, more recently, as next-generation electronic material. The apparent contrast between inertness and stability in ambient, and catalytic activity under reactive conditions is puzzling and calls for a synergistic theoretical and experimental investigation with the long-term objective of enabling the rational design of MoS 2 -based catalysts for alcohol-formation reactions by providing a microscopic understanding of the environmental factors that determines site activity and selectivity. Our research project seeks answers to the questions (a) what conformation does MoS 2 adopt under reaction conditions (as opposed to that under ultrahigh vacuum and low temperatures)?; (b) what reaction pathways exist on such a material?; (c) how can the local environment of the active sites be manipulated so as to make MoS 2 an efficient catalyst for production of higher alcohol from syngas? In particular, research strategies will explore how the basal plane composed of sulfur atoms can be activated so as to exhibit a reactivity of its own, by addressing three research targets and building on extensive preliminary and enabling work: (1) vacancies and vacancy aggregates on the basal plane; (2) non-local catalyst transformation through alkali doping, hydrogenation and phase transition; (3) fabrication of a metal-nanoparticle-activated MoS 2 system, in which particle anchoring, reactive sites and pathways as well as selectivity are controlled by design, as an example of predictive development of a catalyst material,. All strategies are directed to improve the efficacy of the key reactive sites and selectivity of chemical pathways by design, to replace the inefficient methodology of trial and error in catalyst development. This research project represents a synergistic combination of computational guidance, foundational surface-science-based experiments and validation under reactive conditions that aims at transformative new insights into the working of MoS 2 -based hydrogenation-catalysts, a topic squarely at the center of the interest of DOE BES. Guided and led by Talat Rahman, a computational physicist, this project will apply density functional theory to understand structure, reaction pathways and chemical potentials associated with MoS 2 -based CO/CO 2 hydrogenation, augmented by kinetic Monte Carlo methods for reaction rates and prefactors as well as ab-initio molecular dynamics for evaluation of thermal stability. Complementary experimental input and validation will originate from co-PI Ludwig Bartels, a physical chemist and materials scientist, whose group focuses on local imaging and preparation of MoS 2 materials, from co-PI Peter Dowben, an experimental physicist, whose group is expert in the spectroscopy of occupied and unoccupied electronic states, and from senior collaborator Michael White of Brookhaven National Laboratory, whose group generates high-resolution electronic and activity information on size-selected well-defined metal chalcogenide clusters. This collaborative effort will enable a comprehensive understanding of the correlation of structural integrity and catalytic activity of MoS 2 in forms ranging from extended films to individual particles with known geometries and binding sites. Alcohol formation from syngas is a rapidly emerging application that has great potential through facile, economic, and decentralized biomass gasification. CO 2 activation is one of the most pressing concerns of our time: increasing CO 2 levels in the atmosphere change the climate and expose the globe to environmental transformations with the potential for enormous economic and societal impact. We will investigate CO/CO 2 hydrogenation via an interdisciplinary research collaboration with established synergy – one that, in accordance with the mission of the DOE, involves accredited Hispanic-Serving Institutions and that, through student exchange with international collaborators and National Labs, will directly benefit a broad spectrum of communities and generate human resources in sciences essential to their future.

2D materials↗

Generative $β$-hairpin design using a residue-based physicochemical property landscape

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the -hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial -hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The confluence of machine learning and multiscale simulations

Multiscale modeling has a long history of use in structural biology, as computational biologists strive to overcome the time- and length-scale limits of atomistic molecular dynamics. Contemporary machine learning techniques, such as deep learning, have promoted advances in virtually every field of science and engineering and are revitalizing the traditional notions of multiscale modeling. Deep learning has found success in various approaches for distilling information from fine-scale models, such as building surrogate models and guiding the development of coarse-grained potentials. However, perhaps its most powerful use in multiscale modeling is in defining latent spaces that enable efficient exploration of conformational space. In conclusion, this confluence of machine learning and multiscale simulation with modern high-performance computing promises a new era of discovery and innovation in structural biology.

59 BASIC BIOLOGICAL SCIENCES↗

Crash Simulation of a Boeing 737 Fuselage Section Vertical Drop Test

A 30-ft/s vertical drop test of a fuselage section of a Boeing 737 aircraft was conducted in October of 1999 at the FAA Technical Center in Atlantic City, NJ. This test was performed to evaluate the structural integrity of a conformable auxiliary fuel tank mounted beneath the floor and to determine its effect on the impact response of the airframe structure and the occupants. The test data were used to compare with a finite element simulation of the fuselage structure and to gain a better understanding of the impact physics through analytical/experimental correlation. To perform this simulation, a full-scale 3-dimensional finite element model of the fuselage section was developed using the explicit, nonlinear transient-dynamic finite element code, MSC.Dytran. The emphasis of the simulation was to predict the structural deformation and floor-level acceleration responses obtained from the drop test of the B737 fuselage section with the auxiliary fuel tank.

Fasanella, Edwin L.↗

A Design Methodology for Optimizing and Integrating Composite Materials in Gear Structures

The application of composite materials to gear structures is complex because of the gear shape, the need for precise dimensional tolerances, and the complex dynamic load condition. Methods are presented in this work to design and optimize an integrated composite hub-web structure that can be used as part of a hybrid composite-steel gear. The composite hub-web structure is a planar structure with a large decrease in thickness from the hub to the rim. Methods for design and optimization of this variable-thickness structure are presented along with a method for forming braided prepreg material to conform to the shape of this structure. A layered approach is presented that integrates cut plies or filler materials for thickness buildup with continuous-fiber layers for the primary load path. An additional gear design concept is presented that has an axially extended continuous-fiber composite structure that can be combined with the planar structure. A proposed optimization methodology is introduced where an optimized design is output for each type of optimization simulation. Through this process, composite knowledge can be incorporated in the optimization, and more control is given over the design during the process. The methods in this study provide a tool for designing and optimizing composite structures for gears and other high-power-density applications.

optimization↗

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune “bond” lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol’s role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA’s conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Altogether, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational design of closely related proteins that adopt two well-defined but structurally divergent folds

The plasticity of naturally occurring protein structures, which can change shape considerably in response to changes in environmental conditions, is critical to biological function. While computational methods have been used for de novo design of proteins that fold to a single state with a deep free-energy minimum [P.-S. Huang, S. E. Boyken, D. Baker, Nature 537, 320–327 (2016)], and to reengineer natural proteins to alter their dynamics [J. A. Davey, A. M. Damry, N. K. Goto, R. A. Chica, Nat. Chem. Biol.13, 1280–1285 (2017)] or fold [P. A. Alexander, Y. He, Y. Chen, J. Orban, P. N. Bryan, Proc. Natl. Acad. Sci. U.S.A. 106, 21149–21154 (2009)], the de novo design of closely related sequences which adopt well-defined but structurally divergent structures remains an outstanding challenge. We designed closely related sequences (over 94% identity) that can adopt two very different homotrimeric helical bundle conformations—one short (~66 Å height) and the other long (~100 Å height)—reminiscent of the conformational transition of viral fusion proteins. Crystallographic and NMR spectroscopic characterization shows that both the short- and long-state sequences fold as designed. We sought to design bistable sequences for which both states are accessible, and obtained a single designed protein sequence that populates either the short state or the long state depending on the measurement conditions. The design of sequences which are poised to adopt two very different conformations sets the stage for creating large-scale conformational switches between structurally divergent forms.

59 BASIC BIOLOGICAL SCIENCES↗