Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics 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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Virtual texture analysis to understand microstructure effects on deformation twinning and detwinning behavior in BCC metals

Understanding and predicting deformation twinning contributions to plastic deformation in BCC metals has been a long-standing challenge due to the interplay with dislocation slip and non-Schmid effects that render an asymmetry under tension and compression. This paper uses molecular dynamics simulations to understand the effect of unloading and grain orientation on deformation twinning in a nanocrystalline Fe (nc-Fe) system as a model BCC metal. The nc-Fe system is loaded under uniaxial stress tension and compression to understand the effect of grain orientation (Schmid effects) on deformation twinning behavior and the tension/compression asymmetry (non-Schmid effects). A new virtual texture analysis “VirTex” tool is used to understand the role of grain orientation on the nucleation of twins and their contributions to the observed stress–strain response. For certain grain orientations, the twinnability is observed to be different in tension and compression. In addition, the flow stress accommodation from twins in certain grains is observed to be different in tension and compression and different from that for the grains. Subsequent unloading leads to detwinning in the deformed microstructures, where the extent of detwinning depends on the strain from which the system is unloaded and on the morphology of the twin. Lastly, the simulations are carried out to analyze the role of the Schmid factor on the twinnability and asymmetry in tension and compression.

Kannan, Aadhithyan [University of Connecticut, Sto↗
DOE OSTI2025

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗
DOE OSTI2023

Plastic homogeneity in nanoscale heterostructured binary and multicomponent metallic eutectics: An overview

Heterostructured materials comprised of relatively soft/hard disparate phases typically exhibit composite strengthening but lack plastic deformability at ambient temperatures. However, heterostructured systems comprised of nanoscale phases can simultaneously enhance yield strength and strain hardening, thereby promoting uniform distribution of plastic flow. In this review, the atomic-scale deformation mechanisms in model systems of eutectic alloys, Al-Al 2 Cu and Al-Si, refined to nanoscales via laser rapid solidification are discussed, and compared with literature on multi-component (high entropy) eutectics such as Ni-Al-Fe-based with Cr and/or Co additions. The nano-lamellar Al-Al 2 Cu structures exhibit unit defect mechanisms not reported in monolithic Al 2 Cu intermetallic: localized shear on {0 1 1} and shear-induced faults on {1 2 1} planes, constrained by closely-spaced dislocation arrays in Al confined by Al/Al 2 Cu interfaces. The unexpected plasticity mechanisms are enabled by slip continuity in nanoscale Al-Al 2 Cu eutectics associated with the orientation relationship and interface habit planes. In nano-fibrous Al-Si eutectic, tensile ductility at strength approaching 600 MPa is observed resulting from dislocation plasticity in the nano-Al channels and cracking in Si nanofibers. Molecular dynamics simulations show that Al dislocations easily cross-slip (screw) or climb (edge) along Al-Si interfaces, making slip transmission difficult. The propagation of nano-cracks is suppressed by surrounding strain hardening Al, retaining good ductility of the sample, in spite of lack of direct slip transmission. Finally, the critical unit mechanisms of slip transmission and interface-enabled plasticity observed in nanoscale eutectic binary systems can also explain the strength-ductility relationship in multi-component eutectics and homogeneously distributed plastic flow with increasing microstructural heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗
DOE OSTI2022

Opening twisted polymer chains for simultaneously high printability and battery fast-charge

Fast-charging is key to the widespread adoption of battery-based electric vehicles. However, improving fast-charging through architecture optimization is expensive. To reduce costs and expand to commercialization, we applied roll-to-roll screen printing technology to create channels and decrease the tortuosity of electrodes. For the first time, this work successfully opened twisted polymer chains within high-solid-content inks to improve their screen printability and battery performance of as-printed electrodes. With LiNi 0.6 Mn 0.2 Co 0.2 O 2 as active materials, the 60% solid content ink presents superior screen printability after opening the twisted binder chains. As-printed electrode exhibits 33% higher charge capacity at 6 C than printed electrode with chains twisted ink at mass loading of 6.5 mg/cm 2 . Furthermore, coarse-grained molecular dynamics simulations are performed to study the underlying mechanism systematically. In conclusion, the new ink preparation procedure provides a scalable, effective strategy for manufacturing screen-printable battery ink and promotes screen-printed electrode technology.

25 ENERGY STORAGE↗
DOE OSTI2022

Regimes of evaporation and mixing behaviors of nanodroplets at transcritical conditions

The objective of this paper is to examine the fundamental mechanisms responsible for the transition between subcritical evaporation and supercritical dense-fluid-mixing in the absence of convection effects, specifically focusing on the liquid–vapor interfacial dynamics. To isolate the dynamics of this transition process, we characterize the different physical behaviors exhibited by an -dodecane nanoscale droplet placed in different nitrogen ambient conditions across the fuel’s critical point. We employ a continuum-based interface-resolving diffuse-interface method to explore the underlying phase-exchange mechanisms that bring about such distinct dynamics. Following the comparison against molecular dynamics simulations and experiments of evaporating droplets and experimental data for vapor–liquid equilibria, a parametric study at various ambient conditions and droplet sizing is performed to identify four regimes of evaporation/mixing behaviors: sub- and supercritical droplet evaporation, and sub- and supercritical dense-fluid-mixing. It is shown that the distinction in the phase-exchange mechanisms in these four regimes are brought about by the different thermodynamic phases the droplet center can exhibit during the evaporation/mixing process: subcritical liquid, supercritical liquid-like, subcritical gaseous, and supercritical gas-like, respectively. It is shown that the subcritical dense-fluid-mixing behavior is a direct result of nanoconfinement of the liquid–vapor interfacial structure and thus is not present for large droplet sizes. Finally, the present study also shows that the supercritical phase-exchange dynamics can follow two different pathways: supercritical droplet-like evaporation and supercritical dense-fluid-mixing. Furthermore, promoting the early transition to supercritical dense-fluid-mixing can significantly expedite the phase-exchange process through the disintegration of the liquid-like droplet core.

42 ENGINEERING↗
DOE OSTI2025

Tailoring Cu-Zr gradient nanoglass structures: Influence of nanoparticle size and cooling rates on glass-glass interfaces

The study of gradient nanoglasses (GNGs) has gained attention due to their unique mechanical properties and potential applications in advanced materials. This study employs molecular dynamics simulations to synthesize a GNG using Cu-Zr metallic glass nanoparticles (NPs) sized from 3 to 15 nm. The NPs were produced by melting and quenching metallic clusters at a relatively slow quench rate of 10 9 K/s. The synthesis of GNG is elucidated along with the characterization of its heterogeneous metallic glass nanostructure. A seamless GNG structure is formed through cold compression of Cu 64 Zr 36 amorphous NPs of varying sizes. The influence of NP size on the GNG structure is investigated, utilizing deeply relaxed NPs, which exhibit a characteristic Cu segregation pattern on their surfaces. The results highlight an increase in structural heterogeneity due to heterogeneous mass transport and the development of local composition and density variations caused by Cu segregation at glass-glass interfaces (GGIs). A reduction in NP size is correlated with decreased Cu atomic displacements and local density at GGIs, suggesting that larger NPs may produce stronger GGIs. This research presents a novel methodology for synthesizing heterogeneous metallic glasses, demonstrating the capacity to control and customize nanostructure heterogeneity through the manipulation of NP sizes and cooling rates. Furthermore, these findings enhance our understanding of structural evolution during nanoglass synthesis and lay the foundation for further exploration in nanomaterial synthesis and characterization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Insights into the dual cleavage activity of the GH16 laminarinase enzyme class on β-1,3 and β-1,4 glycosidic bonds

Glycoside hydrolases (GH) are involved in the degradation of a wide diversity of carbohydrates and present several biotechnological applications. Many GH families are composed of enzymes with a single well-defined specificity. In contrast, enzymes from the GH16 family can act on a range of different polysaccharides, including β-glucans and galactans. SCLam, a GH16 member derived from a soil metagenome, an endo-β-1,3(4)-glucanase (EC 3.2.1.6), can cleave both β-1,3 and β-1,4 glycosidic bonds in glucans, such as laminarin, barley β-glucan, and cello-oligosaccharides. A similar cleavage pattern was previously reported for other GH16 family members. However, the molecular mechanisms for this dual cleavage activity on (1,3)- and (1,4)-β-D-glycosidic bonds by laminarinases have not been elucidated. In this sense, we determined the X-ray structure of a presumably inactive form of SCLam co-crystallized with different oligosaccharides. The solved structures revealed general bound products that are formed due to residual activities of hydrolysis and transglycosylation. Biochemical and biophysical analyses and molecular dynamics simulations help to rationalize differences in activity towards different substrates. Our results depicted a bulky aromatic residue near the catalytic site critical to select the preferable configuration of glycosidic bonds in the binding cleft. Altogether, these data contribute to understanding the structural basis of recognition and hydrolysis of β-1,3 and β-1,4 glycosidic linkages of the laminarinase enzyme class, which is valuable for future studies on the GH16 family members and applications related to biomass conversion into feedstocks and bioproducts.

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

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗
DOE OSTI2021

Binding mechanism and binding free energy of amino acids and citrate to hydroxyapatite surfaces as a function of crystallographic facet, pH, and electrolytes

Hydroxyapatite (HAP) is the major mineral phase in bone and teeth. The interaction of individual amino acids and citrate ions with different crystallographic HAP surfaces has remained uncertain for decades, creating a knowledge gap to rationally design interactions with peptides, proteins, and drugs. In this contribution, we quantify the binding mechanisms and binding free energies of the 20 end-capped natural amino acids and citrate ions on the basal (001) and prismatic (010)/(020) planes of hydroxyapatite at pH values of 7 and 5 for the first time at the molecular scale. We utilized over 1500 steered molecular dynamics simulations with highly accurate potentials that reproduce surface and hydration energies of (hkl) hydroxyapatite surfaces at different pH values. Charged residues demonstrate a much higher affinity to HAP than charge-neutral species due to the formation of superficial ion pairs and ease of penetration into layers of water molecules on the mineral surface. Binding free energies range from 0 to -60 kJ/mol and were determined with ~10% uncertainty. The highest affinity was found for citrate, followed by Asp(-) and Glu(-), and followed after a gap by Arg(+), Lys(+), as well as by His(+) at pH 5. The (hkl)-specific area density of calcium ions, the protonation state of phosphate ions, and subsurface directional order of the ions in HAP lead to surface-specific binding patterns. Amino acids without ionic side groups exhibit weak binding, between -3 and 0 kJ/mol, due to difficulties to penetrate the first layer of water molecules on the apatite surfaces. We explain recognition processes that remained elusive in experiments, in prior simulations, discuss agreement with available data, and reconcile conflicting interpretations. In conclusion, the findings can serve as useful input for the design of peptides, proteins, and drug molecules for the modification of bone and teeth-related materials, as well as control of apatite mineralization.

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

Designing sequence-defined peptoids for fibrillar self-assembly and silicification

In the biological environment, mineral crystals exquisitely controlled by biomacromolecules often show intricate hierarchical structures and superior mechanical properties. Among these biominerals, spicules, hybrid silica/protein superstructures serving as skeletal elements in demosponges, represent an excellent example for motivating the synthesis of silica materials. Herein, by designing sequence-defined peptoids containing side chains with a strong binding to silica, we demonstrated that self-assembly of these peptoids into fiber structures enables the mimicking of both biocatalytic and templating functions of silicatein filaments for the formation of silica fibers at near-neutral pH and ambient temperature. We further showed that the presence of amino groups is significant for the nucleation of silica on self-assembled peptoid nanofibers. Molecular dynamics simulation further confirmed that having silica-binding of amino side chains is critical for self-assembled peptoid fibers in triggering silica formation. Here, we demonstrated that tuning inter-peptoid interactions by varying carboxyl and amino side chains significantly influences the assembly kinetics and final morphologies of peptoid assemblies as scaffolds for directing silica mineralization to form silica spheres, fibers, and sheets. The formation of silica shell on peptoid fibers increased the mechanical property of peptoid hydrogel materials by nearly 1000-fold, highlighting the great potential of using silicification to enhance the mechanical property of hydrogel materials for applications including tissue engineering. Since peptoids are highly robust and programmable, we expect that self-assembly of peptoids containing solid-binding side chains into hierarchical materials opens new opportunities in the design and synthesis of highly tunable scaffolds that direct the formation of composite nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Physicochemical control of solvation and molecular assembly of charged amphiphilic oligomers at air-aqueous interfaces

Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM + ) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies.

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

Lattice-scale variations in viscosity are correlated with solution structure at mineral-water interfaces

At solid-liquid interfaces, the viscosity increases markedly from the bulk due to the collective interactions of ions and water molecules, influencing phenomena relevant to nanofluidics, colloidal dynamics, and electrochemistry. Here, in this study, we investigated dissipative forces at the boehmite-water interface using 3D atomic force microscopy. We observed an increase in interfacial solution viscosity, η, by 10-100-fold as the nanoprobe approached the surface in normal direction, with up to four oscillatory features showing average peaks of η/η bulk = 44-71. Moreover, the viscosity showed sub-nanometer variations within 0.5 nm from the interface, templated by the underlying crystal lattice and correlated with the interfacial solution structure. Beyond a near-wall region of approximately 1.2 nm, the dissipative response was comparable to that in bulk solution. Molecular dynamics simulations, along with statistical mechanical analyses, provided details on hydrodynamic structures near the interface. Specifically, the lattice-dependent dissipative responses are correlated with extensive hydrogen bonding by interfacial water molecules, which increased friction, particularly along the [001] direction. These results demonstrate how solution viscosity at mineral-water interfaces is anisotropic and correlated with the local solution structure, providing insights into the dynamics of nanocrystal attachment.

Viscosity↗
DOE OSTI2024

Water-induced surface ordering facilitating the microcutting of ductile metals

Metal cutting is a crucial process in modern manufacturing. Enhancing the machinability of metals can significantly improve their production efficiency and surface integrity. Coating surface-active media (SAM) on the free surface of the metals before cutting is an easy method to improve machinability, which usually pertains to the category of the renowned Rehbinder effect. However, the existing SAM is usually hazardous and complex materials. Besides, the effect of SAM on the local structure of the metal surface remains unclear. In this study, water is employed as a simple yet often overlooked SAM in the microcutting of copper. Using water as SAM also allows the employment of X-ray absorption fine structure spectroscopy (XAFS) to study the local structure of copper with and without water coating. Results show that water coating on the free surface of copper can significantly reduce the cutting force and chip thickness, and improve the surface finish. Interestingly, removing the water coating enables the recovery of the cutting force, demonstrating a reversible effect. Based on the XAFS results and molecular dynamics simulation, a water-induced surface ordering mechanism is proposed to explain the findings from the microcutting experiments. This mechanism suggests that water molecules can induce surface ordering in copper, resulting in reduced surface energy and fracture toughness of copper, thus enhancing machinability. In conclusion, this work provides valuable insights into the comprehension of the Rehbinder effect and shows that picometer-scale modifications of the surface atom arrangement can considerably alter the deformation mode of metals, paving the way for the development of new manufacturing processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep eutectic solvent-based polymer electrolyte for solid-state lithium metal batteries

Poly(ethylene) oxide (PEO)-based electrolytes have been widely studied for solid-state lithium batteries while their ionic conductivity and lithium-ion transference number still need to be further improved. Herein, using the combined experimental and theoretical approach, we demonstrate a novel, solid-state PEO-deep eutectic solvent (DES) electrolyte for the first time. We found that the in situ formation of DES can reduce the crystallinity of PEO matrix and more Li + ions can move freely owing to the weakened coordination between ether oxygens and Li-ions. Besides, we show that more Li + ions can be dissociated from Li salts in PEO-DES electrolyte using the molecular dynamics simulations. Such liquid-free PEO-DES electrolytes showed good ionic conductivity (2.1 × 10 −4 S cm −1 ) which is 160% higher than that of conventional PEO-LiTFSI (8.1 × 10 −5 S cm −1 ) electrolyte at 60 °C. Additionally, the PEO-DES electrolyte showed 136% increase of Li-ion transference number (0.33) compared with ionic liquid-doped PEO-LiTFSI (0.14) at 60 °C. Moreover, the PEO-DES exhibited good compatibility with Li metal and stable Li plating/stripping behavior with little morphology change of Li metal. In conclusion, this research also provides new insights into the enhancement mechanisms of novel polymer electrolytes, improving our fundamental understanding of critical challenges that have impeded the adoption of solid-state lithium metal batteries.

Deep eutectic solvent↗