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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular Design of Functional Polymers for Silica Scale Inhibition

Silica polymerization, which involves the condensation reaction of silicic acid, is a fundamental process with wide-ranging implications in biological systems, material synthesis, and scale formation. The formation of a silica-based scale poses significant technological challenges to energy-efficient operations in various industrial processes, including heat exchangers and water treatment membranes. Despite the common strategy of applying functional polymers for inhibiting silica polymerization, the underlying mechanisms of inhibition remain elusive. In this study, we synthesized a series of nitrogen-containing polymers as silica inhibitors and elucidated the role of their molecular structures in stabilizing silicic acids. Polymers with both charged amine and uncharged amide groups in their backbones exhibit superior inhibition performance, retaining up to 430 ppm of reactive silica intact for 8 h under neutral pH conditions. In contrast, monomers of these amine/amide-containing polymers as well as polymers containing only amine or amide functionalities present insignificant inhibition. Molecular dynamics simulations reveal strong binding between the deprotonated silicic acid and a polymer when the amine groups in the polymer are protonated. Notably, an extended chain conformation of the polymer is crucial to prevent proximity between the interacting monomeric silica species, thereby facilitating effective silica inhibition. Furthermore, the hydrophobic nature of alkyl segments in polymer chains disrupts the hydration shell around the polymer, resulting in enhanced binding with ionized silicic acid precursors compared to monomers. Finally, our findings provide novel mechanistic insights into the stabilization of silicic acids with functional polymers, highlighting the molecular design principles of effective inhibitors for silica polymerization.

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

Effect of Mixed-Solvent Environments on the Selectivity of Acid-Catalyzed Dehydration Reactions

The composition of the liquid phase can alter the rates of individual reaction steps and thus alter the selectivity of acid-catalyzed reactions, but these solvent effects are difficult to anticipate for design purposes. Herein, we report the kinetics and selectivity of Brønsted acid-catalyzed 1,2-propanediol dehydration in pure water and in aqueous mixtures of the polar aprotic cosolvents γ-valerolactone, 1,4-dioxane, tetrahydrofuran, N-methyl-2-pyrrolidone, tetramethylene sulfoxide, and dimethyl sulfoxide at 433 K. We find that the major product of 1,2-propanediol dehydration is propanal in most mixed-solvent environments with selectivities between 1 and 68 mol %. In contrast, 1,2-propanediol dehydration in aqueous mixtures of dimethyl sulfoxide affords acetone as the major product with up to 48% selectivity with minimal propanal formation. We use classical molecular dynamics simulations to probe these solvent effects by computing the difference between the solvation free energies of 1,2-propanediol and propanal in aqueous mixtures of polar aprotic cosolvents and in pure water. We find that the difference in the solvation free energies is correlated with the rates of propanal formation in all mixed-solvent environments, indicating that the solvent-mediated stabilization of the product state relative to the reactant state translates to increased selectivity toward the same product. Similar agreement between simulated solvation free energies and experimental reaction rates/selectivities is established for the acid-catalyzed dehydration of cis- and trans-1,2-cyclohexanediol and 1,3-cyclohexanediol. Finally, analysis of the solvation environment around 1,2-propanediol shows that dimethyl sulfoxide uniquely competes against water to solvate reactive hydroxyl groups, which causes a change in reaction mechanism in this solvent system that leads to the formation of acetone rather than propanal. Here, these results represent a step toward the computationally efficient screening of solvent systems for acid-catalyzed, liquid-phase processes.

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

Mechanistic Insights on Permeation of Water over Iron Cations in Nanoporous Silicon Oxide Films for Selective H 2 and O 2 Evolution

Electrocatalysts encapsulated by an ultrathin and semipermeable oxide layer offer a promising avenue for efficient, selective, and cost-effective production of hydrogen through photoelectrochemical water splitting. This architecture is especially attractive for Z-scheme water splitting, for which a nanoporous oxide film can be leveraged to mitigate undesired, yet kinetically facile, reactions involving redox shuttles, such as aqueous iron cations, by limiting transport of these species to catalytically active sites. In this work, molecular dynamics simulations were combined with electrochemical measurements to provide a mechanistic understanding of permeation of water and Fe(III)/Fe(II) redox shuttles through nanoporous SiO 2 films. It is shown that even for SiO 2 pores with a width as small as 0.8 nm, water does not experience any energy barrier for permeating into the pores due to a favorable interaction with hydrophilic silanol groups on the oxide surface. In contrast, permeation of Fe(III) and Fe(II) into microporous SiO 2 pores is limited due to high energy barriers, which stem from a combination of distortion and dehydration of the second and third ion solvation shells. Furthermore, our simulations and experimental results show that SiO 2 coatings can effectively mitigate undesired Fe(III)/Fe(II) redox reactions at underlying electrodes by attenuating permeation of iron cations, while allowing water to permeate and thus participate in water splitting reactions. In a broader context, our study demonstrates that selectivity of solvated cations can be manipulated by controlling the pore size and surface chemistry of oxide films.

14 SOLAR ENERGY↗
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Membrane hydrophobicity determines the activation free energy of passive lipid transport

The collective behavior of lipids with diverse chemical and physical features determines a membrane's thermodynamic properties. Yet, the influence of lipid physicochemical properties on lipid dynamics, in particular interbilayer transport, remains underexplored. Here, we systematically investigate how the activation free energy of passive lipid transport depends on lipid chemistry and membrane phase. Through all-atom molecular dynamics simulations of 11 chemically distinct glycerophospholipids, we determine how lipid acyl chain length, unsaturation, and headgroup influence the free energy barriers for two elementary steps of lipid transport: lipid desorption, which is rate limiting, and lipid insertion into a membrane. Consistent with previous experimental measurements, we find that lipids with longer, saturated acyl chains have increased activation free energies compared to lipids with shorter, unsaturated chains. Lipids with different headgroups exhibit a range of activation free energies; however, no clear trend based solely on chemical structure can be identified, mirroring difficulties in the interpretation of previous experimental results. Compared to liquid-crystalline phase membranes, gel phase membranes exhibit substantially increased free energy barriers. Overall, we find that the activation free energy depends on a lipid's local hydrophobic environment in a membrane and that the free energy barrier for lipid insertion depends on a membrane's interfacial hydrophobicity. Both of these properties can be altered through changes in lipid acyl chain length, lipid headgroup, and membrane phase. Thus, the rate of lipid transport can be tuned through subtle changes in local membrane composition and order, suggesting an unappreciated role for nanoscale membrane domains in regulating cellular lipid dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Synthesis of metallic glass nanoparticles by inert gas condensation

Inert Gas Condensation (IGC) is a versatile method for the synthesis of ultrafine nanoparticles. While it has been mostly used in the synthesis of nanocrystalline particles it has also been employed in the synthesis of metallic glass nanoparticles, which can be further used in the generation of nanostructured metallic glasses, known as nanoglasses. In this work, we use molecular dynamics simulations to investigate the IGC processes leading to the synthesis of Cu 64 Zr 36 amorphous nanoparticles. We investigate the growth of nanoparticles up to 7.3 nm diameter using long simulations of up to 130 ns, under relatively low pressures in the range of 4 to 10 bars. Results indicate a significant effect of pressure on the final structure of the amorphous nanoparticles, yet no clear influence on the segregation of Cu atoms to its surface. Cu segregation occurring at all pressures results in a thin Cu rich layer, about 2.5 Å thick, that coats all nanoparticles. Here, statistics of atomic Voronoi polyhedra shows that the fraction of full icosahedra in the nanoparticles produced at 4 bars is equivalent to that of bulk metallic glasses produced by melt and quench at 10 9 K/s. Higher pressures result in a lower fraction of full icosahedra.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phosphorylation barcodes direct biased chemokine signaling at CXCR3

Chemokine receptors, a group of G protein-coupled receptors (GPCRs), interact with transducers such as G proteins, ß-arrestins, and GPCR kinases (GRKs). In the chemokine system, many chemokine agonists act as “biased agonists” that preferentially activate distinct signaling effectors when binding to the same receptor, resulting in distinct physiological effects. Although one third of FDA-approved drugs target GPCRs, there has been limited success in targeting the chemokine system. Currently, there is little evidence that differential receptor phosphorylation, or “phosphorylation barcodes,” direct these biased responses at chemokine receptors. To address this knowledge gap, we used mass spectrometry to demonstrate that chemokines of CXCR3 promote different ensembles of phosphorylation barcodes that are associated with differential activation of G proteins, ß-arrestins and GRKs. Chemokine stimulation also resulted in distinct changes throughout the kinome in global phosphoproteomic studies. Mutation of specific CXCR3 phosphosites altered ß-arrestin conformation and impacted ß-arrestin activation in molecular dynamics simulations. T-cells expressing phosphorylation-deficient CXCR3 mutants resulted in distinct agonist- and receptor-specific chemotactic and signaling profiles that were not completely explained by engagement of G proteins, ß-arrestins, and GRKs alone. In conclusion, our results directly link distinct GPCR phosphorylation patterns with non-redundant chemokine signaling (Figure 1).

CXCR3↗
DOE OSTI2021

Solvation Structure and Dynamics of Mg(TFSI) 2 Aqueous Electrolyte

Using ab initio molecular dynamics (AIMD) simulations, classical molecular dynamics (CMD) simulations, small-angle X-ray scattering (SAXS), and pulsed-field gradient nuclear magnetic resonance (PFG-NMR), the solvation structure and ion dynamics of magnesium bis(trifluoromethanesulfonyl)imide (Mg(TFSI) 2 ) aqueous electrolyte at 1, 2, and 3 m concentrations are investigated. From AIMD and CMD simulations, the first solvation shell of an Mg 2+ ion is found to be composed of six water molecules in an octahedral configuration and the solvation shell is rather rigid. The TFSI - ions prefer to stay in the second solvation shell and beyond. Meanwhile, the comparable diffusion coefficients of positive and negative ions in Mg(TFSI) 2 aqueous electrolytes have been observed, which is mainly due to the formation of the stable [Mg(H 2 O) 6 ] 2+ complex, and, as a result, the increased effective Mg ion size. Finally, the calculated correlated transference numbers are lower than the uncorrelated ones even at the low concentration of 2 and 3 m, suggesting the enhanced correlations between ions in the multivalent electrolytes. This work provides a molecular-level understanding of how the solvation structure and multivalency of the ion affect the dynamics and transport properties of the multivalent electrolyte, providing insight for rational designs of electrolytes for improved ion transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Toward the Prediction and Control of Glass Transition Temperature for Donor–Acceptor Polymers

Semiconducting donor–acceptor (D–A) polymers have attracted considerable attention toward the application of organic electronic and optoelectronic devices. However, a rational design rule for making semiconducting polymers with desired thermal and mechanical properties is currently lacking, which greatly limits the development of new polymers for advanced applications. Here, polydiketopyrrolopyrrole (PDPP)-based D–A polymers with varied alkyl side-chain lengths and backbone moieties are systematically designed, followed by investigating their thermal and thin film mechanical responses. The experimental results show a reduction in both elastic modulus and glass transition temperature (Tg) with increasing side-chain length, which is further verified through coarse-grained molecular dynamics simulations. Informed from experimental results, a mass-per-flexible bond model is developed to capture such observation through a linear correlation between Tg and polymer chain flexibility. Using this model, a wide range of backbone Tg over 80 °C and elastic modulus over 400 MPa can be predicted for PDPP-based polymers. This study highlights the important role of side-chain structure in influencing the thermomechanical performance of conjugated polymers, and provides an effective strategy to design and predict Tg and elastic modulus of future new D–A polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interactions between Per- and Polyfluoroalkyl Substances (PFAS) at the Water–Air Interface

Per- and polyfluoroalkyl substances (PFAS)--so-called “forever chemicals”--contaminate the drinking water of about 100 million people in the U.S. alone and are inefficiently removed by standard treatment techniques. A key property of these compounds that underlies their fate and transport and the efficacy of several promising remediation approaches is that they accumulate at the water−air interface. This phenomenon remains incompletely understood, particularly under conditions relevant to natural and treatment systems where water−air interfaces often carry significant loads of other organic contaminants or natural organic matter. To understand the impact of organic loading on PFAS adsorption, we carried out molecular dynamics simulations of PFAS at varying interfacial densities. We find that adsorbed PFAS form strong mutual interactions (attraction between perfluoroalkyl chains and electrostatic interactions among charged head groups) that give rise to ordered interfacial coatings. Furthermore, these interactions often involve near-cancellation of hydrophobic attraction and Coulomb repulsion. Our findings explain an apparent paradox whereby PFAS adsorption isotherms often suggest minimal mutual interactions while simultaneously displaying a high sensitivity to the composition and density of interfacial coatings. Consideration of the compounds present with PFAS at the interface has the potential to allow for more accurate predictions of fate and transport and the design of more efficient remediation approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Molecular Origins of Temperature-Dependent X-ray Absorption in YAG:Dy

The absorption and emission of X-rays in dysprosium-doped yttrium aluminum garnet (YAG:Dy) has produced unexpected thermographic behavior, which is investigated using a combination of finite temperature ab initio molecular dynamic simulations, structural characterization, and electronic structure calculations of X-ray characteristics. Calculated average peak X-ray absorption spectra (XAS) from simulations between 300 and 600 K result in peak intensity loss due to thermalization effects, matching experimentally measured behavior of YAG:Dy. Here, investigation of atomic snapshots indicates structural factors that correlated with the X-ray behavior, with the first Y–O coordination sphere identified as the primary structural feature unique to high XAS intensity as calculated by radial and pair distribution functions.

Ab initio molecular dynamics↗
DOE OSTI2024

A strain density function to analyze particle size effects during high velocity impacts of yttria

Here, the impact behavior of a single yttria (Y 2 O 3 ) nanoparticle onto a Y 2 O 3 substrate is studied as a function of particle velocity (300–1200 m/s) and diameter (12–50 nm) using molecular dynamics simulations. To analyze the results, a strain density function is developed to provide both quantitative and qualitative information about the deformation mechanisms that contribute to the final state of the system. This function provides both clear evidence of shear localization and insight into the conditions for which localization occurs during the impact of Y 2 O 3 particles as they approach experimentally relevant sizes. Further analysis shows that localization and fragmentation only occur for Y 2 O 3 impacts with diameters of at least 50 nm, while particles with diameters of 25 nm and smaller deform primarily by amorphization and viscous flow. Implications for the deposition of films of Y 2 O 3 and other rare-earth oxides via the micro cold spray process are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Revealing the structure and electronic characteristics of Te-rich threshold switching materials for high-density integration

Ovonic threshold switching (OTS) selectors play an important role in the integration of advanced three-dimensional memory. Selectors based on tellurium (Te)-containing materials exhibit significant promise due to their low threshold voltages and superior consistency. Here we have theoretically studied the structure and electronic properties of a typical OTS material, amorphous GeTe 6 , to explore the switching mechanisms using ab initio molecular dynamics simulations. The results indicate that Ge atoms tend to bond with Te atoms, forming stable chemical bonds. The Te-centered clusters are predominantly in the form of distorted octahedrons, while the Ge-centered clusters are in the form of both octahedrons and tetrahedrons. Notably, the proportion of tetrahedrons within the 4-coordinated Ge-centered clusters reaches an impressive 66.9%. These tetrahedrons are randomly dispersed throughout the simulated cell, leading to a stable amorphous configuration. The mid-gap state observed in the mobility bandgap originates from the atomic chain composed of both over-coordinated Ge and Te atoms. It is the inherent stability of the chemical environment within amorphous GeTe 6 that enables it to maintain its amorphous phase under a repeated threshold voltage, a characteristic that distinguishes it from non-OTS materials, such as amorphous tellurium. Furthermore, our findings provide an in-depth understanding of the structure and electronic characteristics of amorphous GeTe 6 , which can promote the design and application of the Te-rich threshold switching materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗