Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural 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 649 records · Page 36

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

Structure of the Diphtheria Toxin at Acidic pH: Implications for the Conformational Switching of the Translocation Domain

Diphtheria toxin, an exotoxin secreted by Corynebacterium that causes disease in humans by inhibiting protein synthesis, enters the cell via receptor-mediated endocytosis. The subsequent endosomal acidification triggers a series of conformational changes, resulting in the refolding and membrane insertion of the translocation (T-)domain and ultimately leading to the translocation of the catalytic domain into the cytoplasm. Here, we use X-ray crystallography along with circular dichroism and fluorescence spectroscopy to gain insight into the mechanism of the early stages of pH-dependent conformational transition. For the first time, we present the high-resolution structure of the diphtheria toxin at a mildly acidic pH (5–6) and compare it to the structure at neutral pH (7). We demonstrate that neither catalytic nor receptor-binding domains change their structure upon this acidification, while the T-domain undergoes a conformational change that results in the unfolding of the TH2–3 helices. Surprisingly, the TH1 helix maintains its conformation in the crystal of the full-length toxin even at pH 5. This contrasts with the evidence from the new and previously published data, obtained by spectroscopic measurements and molecular dynamics computer simulations, which indicate the refolding of TH1 upon the acidification of the isolated T-domain. The overall results imply that the membrane interactions of the T-domain are critical in ensuring the proper conformational changes required for the preparation of the diphtheria toxin for the cellular entry.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular dynamics simulations for glass transition temperature predictions of polyhydroxyalkanoate biopolymers

Polyhydroxyalkanoates (PHAs) represent an emerging class of biosynthetic and biodegradable polyesters that exhibit considerable potential to replace petroleum-based plastics towards a sustainable future. Despite the promise, general structure–property mappings within this class of polymers remain largely unexplored. An efficient exploration of this vast chemical space calls for the development and validation of predictive methods for accurate estimation of a diverse range of properties for PHA-based polymers. Towards this aim, we present and validate in this work the results of our molecular dynamics (MD) simulation based approach aimed at predicting glass transition temperatures (T g ) of PHA-based polymers. Since generally available and widely used polymer forcefields exhibit a relatively poor performance for T g predictions, we have developed a new forcefield by modifying the polymer consistent force field (PCFF) via refining a selected set of torsion potentials of the polymer backbone using accurate density functional theory (DFT) computations. After carefully assessing the dependence of critical simulation parameters, such as, polymer chain length, number of polymer chains, supercell size, and thermal quenching rate used in the simulation, the applicability and transferability of the modified PCFF (mPCFF) is demonstrated by directly comparing the computed T g predictions of various polymers with different chemistries, polymer side chain lengths and functional groups forming the polymer side chains against the respective experimentally measured values. Furthermore, the transport properties such as self-diffusion coefficient and viscosity are computationally determined and their well-known correlation with the target properties is demonstrated. Lastly, we have employed the developed approach to predict T g values for a number of yet-to-be-synthesized PHA-based polymers with a diverse set of functional groups in the polymer side chains. The results are further rationalized by correlating the predicted T g values with the inter-chain H-bond formation tendencies of the different side chain functional groups. This work represents an important first step towards computationally guided design of PHA-based functional polymers and opens up new directions for a systematic investigation of composition- and configuration-dependent structure–property relationships in more complex binary and ternary copolymer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alkali Cation Inhibition of Imidazolium-Mediated Electrochemical CO 2 Reduction on Silver

Imidazolium-based ionic liquids have led to enhanced CO 2 electroreduction activity due to cation effects at the cathode surface, stabilizing the reaction intermediates and decreasing the activation energy. In aqueous media, alkali cations are also known to improve CO 2 reduction activity on metals such as Ag, with the enhancement attributed to electrical double layer effects and trending with the size of the alkali cation. However, the effect of a mixed catholyte solution of alkali cations in the presence of an imidazolium-based ionic liquid has not been well-explored. Herein, 1-ethyl-3-methylimidazolium tetrafluoroborate, [EMIM][BF 4 ], in water was investigated with alkali salts to unravel the interaction effects for CO 2 electroreduction on Ag. Although both [EMIM] + and alkali cations have individually improved CO 2 to CO conversion on Ag in water, electrochemical results showed that alkali cations hindered imidazolium-mediated CO 2 electroreduction in most conditions. Li + , in particular, was sharply inhibitory compared to other alkali cations and strongly redirected the selectivity to hydrogen evolution. The nature of the alkali cation inhibition was investigated with spectroscopic techniques, including in situ surface-enhanced Raman spectroscopy (SERS) and dynamic electrochemical impedance spectroscopy (DEIS). Along with computational insights from density functional theory (DFT), the electrochemical and spectroscopic data suggest that alkali cations inhibit [EMIM]-mediated CO 2 reduction by competing for surface adsorption sites, preventing the potential-dependent structural reorientation of imidazolium, and promoting hydrogen evolution by bringing solvated water to the cathode surface.

cations↗

In silico λ-dynamics predicts protein binding specificities to modified RNAs

Abstract RNA modifications shape gene expression through a variety of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Biochemistry & Molecular Biology↗

OpenTurbine [SWR-23-07]

OpenTurbine is an open-source code designed for the simulation of wind turbine structural dynamics. It is a flexible multi-body dynamics solver wherein the wind turbine is modeled as a collection of beam finite elements, rigid bodies, and constraints. OpenTurbine is equipped with an application programming interface primarily designed for coupling to computational fluid dynamics codes, but it is also suited for coupling to aerodynamics solvers based on blade element momentum theory.

Sprague, MichaelA.↗

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗

Generic, Sparse Tensor Core for Neural Networks

Sparse neural network attracts more attention for model compression, fast execution, and power reduction. The state-of-the-art designed sparse tensor core for structured and static sparsity, which did not support well for generic or dynamic sparsity. We design a sparse tensor core to support generic sparsity pruning with a novel hybrid and blocked sparse matrix storage format, HB-ELL, which saves computation and storage while keeping the most significant elements, as well as supporting dynamic sparsity for data flow in neural networks. We achieve better performance with preliminary results than the state-of- the-art on an NVIDIA GPU simulator.

Wu, Xiaolong↗

Chapter 21. Atomic-Level Dislocation Dynamics in Irradiated Metals

Primary damage and microstructure evolution in structural nuclear materials operating under conditions of a high flux of energetic atomic particles and high temperature and stress lead to the formation of a high concentration, non-homogeneous distribution of defect clusters in the form of dislocation loops, voids, gas-filled bubbles and radiation-induced precipitates of nanometer scale. They cause changes in many material properties. Being obstacles to dislocation glide, they strongly affect the mechanical properties, with an increase in yield and flow stresses and a reduction in ductility. Atomic-scale computer simulations can provide details of how these effects are influenced by the obstacle structure, applied stress, strain rate and temperature. Processes such as obstacle cutting, transformation, absorption and drag are observed. Some recent results for body-centered and face-centered cubic metals are described in this review and, where appropriate, comparisons are drawn with predictions based on the elasticity theory of crystal defects. Perspectives on how to use this information at higher scales and in particular in mesoscale, dislocation dynamics simulations are also discussed.

Osetskiy, Yury N.↗

A general quantum algorithm for open quantum dynamics demonstrated with the Fenna-Matthews-Olson complex

Using quantum algorithms to simulate complex physical processes and correlations in quantum matter has been a major direction of quantum computing research, towards the promise of a quantum advantage over classical approaches. In this work we develop a generalized quantum algorithm to simulate any dynamical process represented by either the operator sum representation or the Lindblad master equation. We then demonstrate the quantum algorithm by simulating the dynamics of the Fenna-Matthews-Olson (FMO) complex on the IBM QASM quantum simulator. This work represents a first demonstration of a quantum algorithm for open quantum dynamics with a moderately sophisticated dynamical process involving a realistic biological structure. We discuss the complexity of the quantum algorithm relative to the classical method for the same purpose, presenting a decisive query complexity advantage of the quantum approach based on the unique property of quantum measurement.

97 MATHEMATICS AND COMPUTING↗

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗

Self-inhibited State of Venezuelan Equine Encephalitis Virus (VEEV) nsP2 Cysteine Protease: A Crystallographic and Molecular Dynamics Analysis

The Venezuelan equine encephalitis virus (VEEV) belongs to the Togaviridae family and is pathogenic to both humans and equines. The VEEV non-structural protein 2 (nsP2) is a cysteine protease (nsP2pro) that processes the polyprotein and thus it is a drug target for inhibitor discovery. The atomic structure of the VEEV nsP2 catalytic domain was previously characterized by both X-ray crystallography and computational studies. A modified nsP2pro harboring a N475A mutation in the N terminus was observed to exhibit an unexpected conformation: the N-terminal residues bind to the active site, mimicking binding of a substrate. The large conformational change of the N terminus was assumed to be induced by the N475A mutation, as N475 has an important role in stabilization of the N terminus and the active site. This conformation was first observed in the N475A mutant, but we also found it while determining a crystal structure of the catalytically active nsP2pro containing the wild-type N475 active site residue and K741A/K767A surface entropy reduction mutations. This suggests that the N475A mutation is not a prerequisite for self-inhibition. Here, we describe a high resolution (1.46 Å) crystal structure of a truncated nsP2pro (residues 463–785, K741A/K767A) and analyze the structure further by molecular dynamics to study the active and self-inhibited conformations of nsP2pro and its N475A mutant. A comparison of the different conformations of the N-terminal residues sheds a light on the interactions that play an important role in the stabilization of the enzyme.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗

Cyclic Peptides for Lanthanide Binding

Lanthanide ions are difficult to separate from one another due to their similar chemical properties. The discovery of lanthanide-binding peptides and proteins in nature has led to an increased interest in the possibility of utilizing the strong binding of peptides to lanthanide ions for their separations; as such, there has been an effort to identify or design peptides with improved lanthanide binding and selectivity toward particular lanthanide ions. Here, in this study, we designed and characterized lanthanide-binding cyclic peptides (LBCPs) with molecular dynamics simulations, electronic structure calculations, and emission spectroscopy. Luminescent decay measurements were done to determine the number of water molecules coordinated to the Eu 3+ ion in Eu-LBCP complexes and compare to the predicted number of water molecules by computation to assess the lanthanide-binding affinity of LBCPs. Measured stability constants show binding of the LBCPs to the Eu 3+ ion with stronger than micromolar affinity. We were able to identify multiple peptides that selectively bind to middle lanthanides. We describe the structural basis of the lanthanide-binding selectivity trend with strongest binding to the middle lanthanides, followed by the heavier lanthanides, and finally to the lighter ions.

ions↗

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗