Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Synergistic Enhancement of Lead and Selenate Uptake at the Barite (001)–Water Interface

The interactions of heavy metals with minerals influence the mobility and bioavailability of toxic elements in natural aqueous environments. The sorption of heavy metals on covalently bonded minerals is generally well described by surface complexation models (SCMs). However, understanding sorption on sparingly soluble minerals is challenging because of the dynamically evolving chemistry of sorbent surfaces. The interpretation can be even more complicated when multiple metal ions compete for sorption. In the present study, we observed synergistically enhanced uptake of lead and selenate on the barite (001) surface through two sorption mechanisms: lattice incorporation that dominates at lower coverages and two-dimensional monolayer growth that dominates at higher coverages. Further, we also observed a systematic increase in the sorption affinity with increasing co-sorbed ion coverages, different from the assumption of invariant binding constants for individual adsorption processes in classical SCMs. Computational simulations showed thermodynamically favorable co-incorporation of lead and selenate by simultaneously substituting for barium and sulfate in neighboring sites, resulting in the formation of molecular clusters that locally match the net dimension of the substrate lattice. These results emphasize the importance of ion–ion interactions at mineral–water interfaces that control the fate and transport of contaminants in the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Finally, library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Solvation Dynamics of Wet Ethaline: Water is the Magic Component

The past two decades witnessed the development of a new type of solvent system, named deep eutectic solvents, which have become increasingly investigated because they offer new and potentially favorable properties, such as wide tunability in electrochemical, mechanical, and transport properties. Deep eutectic solvent (DES) systems are composed of at least one main solvent and an additional component that is meant to interrupt the original solvent/solvent interactions, thereby introducing lower melting points relative to each individual component. Ethaline (a 1:2 mole% mixture of choline chloride and ethylene glycol) is one of the most promising DES systems. However, it is also known to be very hygroscopic, which is a constant concern because water absorption during the use of ethaline alters its properties. Within this work, we demonstrate that modest amounts of water addition (1-10%) to ethaline are of little concern for practical use and can even lead to performance improvements, such as accelerated relaxation and solvation. In contrast, very small amounts of <1% of water lead to additional slowing of the solvent response. Thus, we suggest that the attempt to dry ethaline below 1% moisture is rather counter-productive if one attempts to achieve effective solvation and charge transport properties from DESs. This study investigates the effect of water content on the diffusional relaxation dynamics of ethaline. A set of independent spectroscopic experiments and computational simulations are aimed to provide insight into the solvent response of the DES system using femtosecond time-resolved absorption spectroscopy (fs-TA), broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR) diffusometry and broadband relaxometry, and molecular dynamics simulations (MDS) on ethaline with 0, 0.1, 1, 10, and 28.5wt% added water. For dry ethaline, we identify choline chloride as the rate-limiting solvation component in ethaline. However, the role of the solvent components changes gradually as water is added. Further, we provide quantitative solvent relaxation rates using the different presented time-resolved spectroscopic techniques and find remarkable agreement between them. Based on the solvent relaxation rates and combined with MDS, we develop a molecular understanding of the individual solvent components and their interactions in dry and wet ethaline with varying amounts of water content.

25 ENERGY STORAGE↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

Adsorption, Orientation, and Speciation of Amino Acids at Air–Aqueous Interfaces for the Direct Air Capture of CO 2

Amino acids make up a promising family of molecules capable of direct air capture (DAC) of CO 2 from the atmosphere. Under alkaline conditions, CO 2 reacts with the anionic form of an amino acid to produce carbamates and deactivated zwitterionic amino acids. The presence of the various species of amino acids and reactive intermediates can have a significant effect on DAC chemistry, the role of which is poorly understood. In this study, all-atom molecular dynamics (MD) based computational simulations and vibrational sum frequency generation (vSFG) spectroscopy studies were conducted to understand the role of competitive interactions at the air–aqueous interface in the context of DAC. We find that the presence of potassium bicarbonate ions, in combination with the anionic and zwitterionic forms of amino acids, induces concentration and charge gradients at the interface, generating a layered molecular arrangement that changes under pre- and post-DAC conditions. In parallel, an enhancement in the surface activity of both anionic and zwitterionic forms of amino acids is observed, which is attributed to enhanced interfacial stability and favorable intermolecular interactions between the adsorbed amino acids in their anionic and zwitterionic forms. The collective influence of these competitive interactions, along with the resulting interfacial heterogeneity, may in turn affect subsequent capture reactions and associated rates. Finally, these effects underscore the need to consider dynamic changes in interfacial chemical makeup to enhance DAC efficiency and to develop successful negative emission and storage technologies.

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

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tunable Nanostructures from Inverse Surfactants

Hierarchical materials in the natural world are often made through the self-assembly of amphiphilic molecules. Achieving similar structural complexity in synthetic materials requires understanding how various molecular parameters affect assembly behavior. In recent years, inverse surfactants─molecules with hydrophobic head groups and hydrophilic macromolecular tails─have been shown to self-assemble into supramolecular assemblies in aqueous solutions that show promise for a number of applications, including drug delivery. Here, we build an understanding of the morphological phase diagram of inverse surfactants using insights from scattering experiments, computer simulations, and statistical mechanics. The scattering and simulation results reveal that changing the headgroup size is an important molecular knob in controlling morphological transitions. The molecular size ratio of the hydrophobic group to the hydrophilic group emerges as a crucial dimensionless quantity in our theory and plays a determining role in setting the micelle structure and the transition from mesoscale to macroscale aggregates. Our minimal theory is able to qualitatively explain the key features of the morphological phase diagram, including the prevalence of fiber-like structures in comparison to spherical and planar micelles. Together, these findings provide a more complete picture of the molecular dependencies of assemblies of inverse surfactants, which we hope may aid in the de novo design of supramolecular structures.

Christakopoulos, Panagiotis [ORNL] (ORCID:00090004↗
DOE OSTI2022

Pinned Electrode/Electrolyte Interphase and Its Formation Origin for Sulfurized Polyacrylonitrile Cathode in Stable Lithium Batteries

Sulfurized polyacrylonitrile (SPAN) represents one of the most promising directions for high-energy-density lithium (Li)-sulfur batteries. However, the practical application of Li||SPAN is currently limited by the insufficient chemical/electrochemical stability of electrode/electrolyte interphase (EEI). Here, a pinned EEI layer is designed for stabilizing a SPAN cathode by regulating the EEI formation mechanism in an advanced LiFSI/ether/fluorinated-ether electrolyte. Computational simulations and experimental investigations reveal that, benefiting from the nonsolvating nature, the fluorinated-ether can not only act as a protective shield to prevent the Li polysulfides dissolution but also, more importantly, endow a diffusion-controlled EEI formation process. It promotes the formation of a uniform, protective, and conductive EEI layer pinning into SPAN surface region, enabling the high loading Li||SPAN batteries with superior cycling stability, wide temperature performance, and high-rate capability. Finally, this design strategy opens an avenue for exploring advanced electrolytes for Li||SPAN batteries and guides the interface design for broad types of battery systems.

25 ENERGY STORAGE↗
DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗
DOE OSTI2021

Origin of Proteolytic Stability of Peptide-Brush Polymers as Globular Proteomimetics

Peptide-brush polymers (PBPs), wherein every side-chain of the polymers is peptidic, represent a new class of proteomimetic with unusually high proteolytic resistance while maintaining bioactivity. Here, we sought to determine the origin of this behavior and to assess its generality via a combined theory and experimental approach. A series of PBPs with various polymer backbone structures were prepared and examined for their proteolytic stability and bioactivity. We discovered that an increase in the hydrophobicity of the polymer backbones is predictive of an elevation in proteolytic stability of the side-chain peptides. Computer simulations, together with small-angle X-ray scattering (SAXS) analysis, revealed globular morphologies for these polymers, in which pendant peptides condense around hydrophobic synthetic polymer backbones driven by the hydrophobic effect. As the hydrophobicity of the polymer backbones increases, the extent of solvent exposure of peptide cleavage sites decreases, reducing their accessibility to proteolytic enzymes. This study provides insight into the important factors driving PBP aqueous-phase structures to behave as globular, synthetic polymer-based proteomimetics.

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

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗
DOE OSTI2020

Elucidation of Active Oxygen Sites upon Delithiation of Li 3 IrO 4

Transformational increases in the storage capacity of battery cathodes could be achieved by tapping into the redox activity at oxide ligands in addition to conventional transition metal couples. However, the key signatures that govern such lattice oxygen redox (LOR) have not been ascertained. Li 3 IrO 4 has the largest reversible LOR, rendering it a unique model system. Here, X-ray spectroscopy and computational simulations reveal that LOR in Li 3 IrO 4 is selectively compensated via O sites with three lone pairs, which are activated by Li/Ir disorder. The two-electron LOR can be reversed to regenerate the initial state without unlocking competing bulk reactions observed in many other compounds. We uncover an intricate interplay between stoichiometry, O coordination, and nonbonding states in LOR and pinpoint spectroscopic signatures. This interplay is indispensable for designing materials with 3d metals that fulfill the promise of LOR to overcome the bottlenecks of current cathodes for future implementation in practical batteries.

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

Control of monodomain polymer-stabilized cuboidal nanocrystals of chiral nematics by confinement.

Liquid crystals are important components of optical technologies. Cuboidal crystals consisting of chiral liquid crystals-the so-called blue phases (BPs), are of particular interest due to their crystalline structures and fast response times, but it is critical that control be gained over their phase behavior as well as the underlying dislocations and grain boundaries that arise in such systems. Blue phases exhibit cubic crystalline symmetries with lattice parameters in the 100 nm range and a network of disclination lines that can be polymerized to widen the range of temperatures over which they occur. Here, we introduce the concept of strain-controlled polymerization of BPs under confinement, which enables formation of strain-correlated stabilized morphologies that, under some circumstances, can adopt perfect single-crystal monodomain structures and undergo reversible crystal-to-crystal transformations, even if their disclination lines are polymerized. We have used super-resolution laser confocal microscopy to reveal the periodic structure and the lattice planes of the strain and polymerization stabilized BPs in 3D real space. Our experimental observations are supported and interpreted by relying on theory and computational simulations in terms of a free energy functional for a tensorial order parameter. Simulations are used to determine the orientation of the lattice planes unambiguously. The findings presented here offer opportunities for engineering optical devices based on single-crystal, polymer-stabilized BPs whose inherent liquid nature, fast dynamics, and long-range crystalline order can be fully exploited.

confinement↗
DOE OSTI2023

DeCOIL: Optimization of Degenerate Codon Libraries for Machine Learning-Assisted Protein Engineering

With advances in machine learning (ML)-assisted protein engineering, models based on data, biophysics, and natural evolution are being used to propose informed libraries of protein variants to explore. Synthesizing these libraries for experimental screens is a major bottleneck, as the cost of obtaining large numbers of exact gene sequences is often prohibitive. Degenerate codon (DC) libraries are a cost-effective alternative for generating combinatorial mutagenesis libraries where mutations are targeted to a handful of amino acid sites. However, existing computational methods to optimize DC libraries to include desired protein variants are not well suited to design libraries for ML-assisted protein engineering. To address these drawbacks, we present DEgenerate Codon Optimization for Informed Libraries (DeCOIL), a generalized method which directly optimizes DC libraries to be useful for protein engineering: to sample protein variants that are likely to have both high fitness and high diversity in the sequence search space. Using computational simulations and wet-lab experiments, we demonstrate that DeCOIL is effective across two specific case studies, with potential to be applied to many other use cases. DeCOIL offers several advantages over existing methods, as it is direct, easy-to-use, generalizable, and scalable. With accompanying software, DeCOIL can be readily implemented to generate desired informed libraries.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Single Iridium Atom Doped Ni 2 P Catalyst for Optimal Oxygen Evolution

Single-atom catalysts (SACs) with 100% active sites have excellent prospects for application in the oxygen evolution reaction (OER). However, further enhancement of the catalytic activity for OER is quite challenging, particularly for the development of stable SACs with overpotentials <180 mV. Here, we report an iridium single atom on Ni 2 P catalyst (Ir SA -Ni 2 P) with a record low overpotential of 149 mV at a current density of 10 mA·cm –2 in 1.0 M KOH. The Ir SA -Ni 2 P catalyst delivers a current density up to ∼28-fold higher than that of the widely used IrO 2 at 1.53 V vs RHE. Both the experimental results and computational simulations indicate that Ir single atoms preferentially occupy Ni sites on the top surface. The reconstructed Ir–O–P/Ni–O–P bonding environment plays a vital role for optimal adsorption and desorption of the OER intermediate species, which leads to marked enhancement of the OER activity. Additionally, the dynamic “top-down” evolution of the specific structure of the Ni@Ir particles is responsible for the robust single-atom structure and, thus, the stability property. In conclusion, this Ir SA -Ni 2 P catalyst offers novel prospects for simplifying decoration strategies and further enhancing OER performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗