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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular understanding of ion transport in a zwitterionic electrolyte

Zwitterions (ZIs) are unique molecules that carry both positive and negative charges, resulting in overall charge neutrality and high dielectric constants. These distinctive properties have enabled broad applications of zwitterionic functionality, including the emerging use of ZIs in lithium-ion battery electrolytes. As a contribution to this developing field, we use all-atom molecular dynamics simulations to investigate the ion transport mechanisms in amorphous mixtures of a zwitterionic liquid containing a range of LiTFSI salt concentrations. Furthermore, the local coordination environment around the Li + ions plays a strong role in governing ionic conductivity, as well as the enhancement of Li + transport numbers with increasing salt concentration. Addition of small amounts of water leads to increased conductivity and ion mobilities due to the water coordinating with the Li + ions, which reduces direct interactions with larger charged species.

Classical molecular dynamic simulations↗
DOE OSTI2022

Strategies for Controlling the Spatial Orientation of Single Molecules Tethered on DNA Origami Templates Physisorbed on Glass Substrates: Intercalation and Stretching

Nanoarchitectural control of matter is crucial for next-generation technologies. DNA origami templates are harnessed to accurately position single molecules; however, direct single molecule evidence is lacking regarding how well DNA origami can control the orientation of such molecules in three-dimensional space, as well as the factors affecting control. Here, we present two strategies for controlling the polar (θ) and in-plane azimuthal (Φ) angular orientations of cyanine Cy5 single molecules tethered on rationally-designed DNA origami templates that are physically adsorbed (physisorbed) on glass substrates. By using dipolar imaging to evaluate Cy5's orientation and super-resolution microscopy, the absolute spatial orientation of Cy5 is calculated relative to the DNA template. The sequence-dependent partial intercalation of Cy5 is discovered and supported theoretically using density functional theory and molecular dynamics simulations, and it is harnessed as our first strategy to achieve θ control for a full revolution with dispersion as small as ±4.5°. In our second strategy, Φ control is achieved by mechanically stretching the Cy5 from its two tethers, being the dispersion ±10.3° for full stretching. These results can in principle be applied to any single molecule, expanding in this way the capabilities of DNA as a functional templating material for single-molecule orientation control. The experimental and modeling insights provided herein will help engineer similar self-assembling molecular systems based on polymers, such as RNA and proteins.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Molecular dynamics studies of grain boundary mobility and anisotropy in BCC γ-uranium

Grain morphologies such as grain size and aspect ratio in uranium-based metallic fuels are important microstructural features that can impact various fuel performance properties such as fission-gas-induced swelling, thermal transport, high burnup structure formation, and radiation resistance. Accurate prediction of the fuel grain morphologies requires knowledge of critical grain growth parameters such as grain boundary (GB) mobility and anisotropy. In this work, molecular dynamics simulations were performed to study the GB mobility and its anisotropy in pure body-centered-cubic (BCC) γ uranium. Nine GBs with different combinations of misorientation angles (20°, 30°, 45°) and rotation axes (<100>, <110>, <111>), as well as an additional <111> 38.2° GB were studied using three interatomic potentials. Here, it is found that the GB mobility anisotropy has complex trends, depending on both rotation axis and misorientation. However, in general the <110> rotation axis has the fastest GB mobility at the same misorientation. The results of this work can be used as not only baselines for future studies of GB mobility in uranium-based alloys such as uranium-molybdenum (U-Mo) fuels, but also input for mesoscale modeling of grain growth in uranium-based alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Morphology–Diffusivity Relationships in Fluorine-Free Random Terpolymers for Proton-Exchange Membranes

Here, using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new family of fluorine-free terpolymers designed as proton-exchange membranes. Simulated random terpolymers consist of three monomers with a 5-carbon backbone with a phenylsulfonate, phenyl, or no pendant group and have ion exchange capacities (IECs) ranging from 1.06–4.14 mmol/g. At a hydration level of 9, cluster analysis reveals macrophase separation between water and terpolymers with IEC < 2.1 mmol/g and continuous, percolated hydrophilic and hydrophobic nanoscale domains at higher IECs. Channel width distribution analysis of the percolated morphologies revealed that more hydrophobic units produce less uniform channels. Decreasing the surface area per sulfonate group and increasing the fractal dimension of the hydrophilic domains correlate with increased water diffusivity, due to a more acidic interface and more isotropic water channels. Relative to the previously studied phenylsulfonate homopolymer, these terpolymers with lower IECs have only modestly lower water diffusion, and we anticipate other advantages related to processability.

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

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Elucidate the molecular basis of ampholytic chitosan as a high-performance cryoprotectant to myosin denaturation: The importance of saccharide charges

The uses of charged poly/oligo-saccharides enable the retarding of protein denaturation against various environmental stresses during food storage and manufacturing. However, at subzero temperatures, the molecular basis of such stabilization behaviors, i.e., cryoprotections, remain less explored. Here, in this study, we introduced an ampholytic saccharide, carboxymethyl chitooligosaccharide (CMCO) that effectively inhibited the freezing-induced myosin denaturation. The in-depth cryoprotective mechanism was systematically investigated by using molecular dynamic simulation and multispectral characterizations. Results showed that CMCO may interact with myosin through hydrogen bonding and electrostatic interactions, which caused the expelling of water at protein surfaces and the reduced conformational flexibility of myosin molecules. Due to this water replacement event, both secondary and ternary structures of myosin became freezing-resistant, leading to the inhibited protein aggregations and retained functionalities, such as solubility, Ca 2+ -ATPase activity, and gelling properties. Moreover, cryoprotective behaviors of CMCO were charge-dependent. CMCO with a higher degree of carboxymethyl substitution (DS: 1.2) was inclined to bind and stabilize myosin molecules better than the low-DS (DS: 0.8) one, even though both outperformed other cryoprotective saccharides. Therefore, this investigation not only introducing a high-performance myosin cryoprotectant, but also elaborated the cryoprotective mechanism of ampholytic saccharides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular Understanding of Charge Storage in MoS 2 Supercapacitors with Ionic Liquids

Owing to high electrical conductivity and ability to reversibly host a variety of inserted ions, 2D metallic molybdenum disulfide (1T-MoS 2 ) has demonstrated promising energy storage performance when used as a supercapacitor electrode. However, its charge storage mechanism is still not fully understood, in particular, how the interlayer spacing of 1T-MoS 2 would affect its capacitive performance. In this work, molecular dynamics simulations of 1T-MoS 2 with interlayer spacing ranging from 0.615 to 1.615 nm have been performed to investigate the resulting charge storage capacity in ionic liquids. Simulations reveal a camel-like capacitance-potential relation, and MoS 2 with an interlayer spacing of 1.115 nm has the highest volumetric and gravimetric capacitance of 118 F cm –3 and 42 F g –1 , respectively. Although ions in MoS 2 with an interlayer spacing of 1.115 nm diffuse much faster than with interlayer spacings of 1.365 and 1.615 nm, the MoS 2 with larger interlayer spacing has a much faster-charging process. Our analyses reveal that the ion number density and its charging speed, as well as ion motion paths, have significant impacts on the charging response. Here, this work helps to understand how the interlayer spacing affects the interlayer ion structures and the capacitive performance of MoS 2 , which is important for revealing the charge storage mechanism and designing MoS 2 supercapacitor.

25 ENERGY STORAGE↗
DOE OSTI2021

New interaction potentials for alkaline earth silicate and borate glasses

New interaction potentials were developed for molecular dynamics simulations to study the role of Mg and Ca in modifying the structure and properties of alkaline earth silicates and borates. Competition between the depolymerization of the silica network and the formation of new bonds between oxygen atoms and modifiers leads to the enhancement of the elastic moduli with increasing modifier content in alkaline earth silicate glasses. Compared with calcium silicate, the higher elastic moduli of magnesium silicate result from a higher connectivity of the overall glass network due to the incorporation of fourfold coordinated magnesium and a more rigid connection between oxygen atoms and modifiers. In contrast to the silicates, the effect of modifier on the elastic moduli of alkaline earth borates is dominated by the formation of fourfold coordinated boron (N 4 ). Calcium borate with higher N 4 shows a more rigid network structure and higher elastic moduli.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CryoFold: Determining protein structures and data-guided ensembles from cryo-EM density maps

Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mechanical properties of heterogeneous metallic glasses: Insights from brick-and-mortar designs

The brick-and-mortar architectural paradigm is commonly applied in structural design to achieve an exceptional combination of strength and plasticity due to its highly tunable mechanical attributes. This study applies molecular dynamics simulations of tensile loading to investigate the mechanical properties and failure mechanisms of brick-and-mortar metallic glasses (BMMGs). The focus is primarily on the implications of the aspect ratio of the bricks and interlayer thickness on the strength, plasticity, and deformation mechanisms. Results indicate that the failure mode is typified by multiple shear bands localized within the softer mortar regions, creating a staggered network under smaller aspect ratios that significantly enhances plastic deformation. However, a transition to a single dominant SB occurs when the aspect ratio exceeds 4.1. In addition, results indicate that a concurrent improvement in strength and plasticity can be achieved by modulating the distance between brick layers in BMMGs. This enhancement originates from the extensive generation and interaction of shear transformation zones. Furthermore, this study highlights that a desirable balance between strength and plasticity can be obtained in BMMGs with appropriate brick aspect ratios and interlayer thicknesses, providing a potential design strategy for the advancement of novel metallic glasses with superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of iron redox ratio on the structures of boroaluminosilicate glasses

Iron oxide is commonly found in natural or industrial glass compositions and can exist as ferrous (Fe 2+ ) or ferric (Fe 3+ ) species, with their ratios depending on glass composition, temperature, pressure and the redox reactions during the glass forming process. The iron redox ratio plays an important role on silicate glass structures and consequently various properties. This work aims to study the effect of iron oxide, and particularly the iron redox ratio, on the structures of borosilicate and boroaluminosilicate glasses using molecular dynamics simulations with newly developed iron potential parameters that are compatible with the borosilicate potentials. The results provide detailed cation coordination states of both iron species and the effect of redox ratio on boron coordination and other structural features. Particularly, competition for charge compensating modifier cations (such as Na + ) among the fourfold-coordinated cations such as B 3+ , Al 3+ , and Fe 3+ is investigated by calculating the cation–cation pair distribution functions and coordination preferential ratios. The results show that the trivalent ferric ions, with a shorter Fe–O bond distance and better defined first coordiation shell with mainly four-fold coordination, act as a glass former whereas the divalent ferrous ions mainly play the role of glass modifier. The ferrous/ferric ratio (Fe 2+ /Fe 3+ ) was found to affect the glass chemistry and hence glass properties by regulating the amount of four-coordinated boron, the fraction of non-briding oxygen and other features. Finally, the results are compared with available experimental data to gain insights of the complex structures and charge compensation schemes of the glass system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Concentric core-shell tracks and spectroscopic properties of SrTiO 3 under intense electronic excitation

A deeper understanding of the intrinsic link between the irradiation-induced microstructures and the corresponding spectroscopic properties is becoming increasingly attractive for the research fields of materials science, physics, and information technology. In this work, the structural damage response, spectroscopic features, and associated physical mechanisms of SrTiO 3 single crystals under swift heavy ion irradiation were comparatively analyzed using a combination of experimental and theoretical approaches. Corresponding to 0.11–5.00 MeV/u ion irradiation with electronic energy loss ranging from 4.0 to 29.3 keV/nm, the inelastic thermal spike calculations combined with molecular dynamics simulations are compared with the experimental observations, revealing the track fine structures (individual spherical defects with a disordered region, and discontinuous and continuous tracks consisting of an amorphous core and a disordered outer shell) and demonstrating the dominant effects of the deposition energy and lattice temperature on track damage formation and evolution; thus, two essential thresholds for defect formation of ~0.60 eV/atom and amorphous region formation of ~1.81 eV/atom were identified to better describe the concentric core-shell track structure. The enhanced nanohillock formation in the surface region is attributed to the combined action of kinetic and potential energy depositions, and the fluence dependence in regulating the hillock dimensions is also presented. The measured refractive index profiles along different crystal axes further indicate the anisotropy of irradiation-induced lattice expansion. With increasing ion fluence and damage level, the intrinsic bandgap gradually decreased, and the concentration of additional radiative recombination centers (1.81 eV, 2.31 eV, and 2.42 eV) accordingly increased in the SrTiO 3 crystal, providing the possibility of regulating related defects to achieve tunable/selective photon emission for more critical technological applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Molecular design and fabrication of PIM-1/polyphosphazene blend membranes with high performance for CO 2 /N 2 separation

New polymeric blend membranes for CO 2 separation were synthesized based on insights from molecular dynamics simulation. A molecular-level structure-property relationship in polymers of intrinsic microporosity (PIM) based blend membranes, was investigated in detail computationally. Calculated local density profiles and energy of interaction of the blend membranes, composed of PIM-1 and various polyphosphazenes, showed that using the polyphosphazene with a higher concentration of ether side chains can improve the compatibility with PIM-1. Furthermore, based on the findings of computational studies, blend membranes were experimentally fabricated from PIM-1 and polyphosphazenes with various polyether side chain concentrations. Polyether concentration in polyphosphazenes was correlated with the film properties and gas transport performance of the blend membranes. Blend membranes showed very high CO 2 permeability (3100-5300 barrer) and improved CO 2 /N 2 selectivity (24-28), outperforming all other PIM-based blend membranes reported to date. Moreover, the CO 2 permeability performance of the blend membranes was tested 566 hours under real post-combustion flue gas from a coal-fired power plant, including CO 2 , N 2 , H 2 O, O 2 , SO x and NO x .

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

Bacterial Membranes Are More Perturbed by the Asymmetric Versus Symmetric Loading of Amphiphilic Molecules

Characterizing the biophysical properties of bacterial membranes is critical for understanding the protective nature of the microbial envelope, interaction of biological membranes with exogenous materials, and designing new antibacterial agents. Presented here are molecular dynamics simulations for two cationic quaternary ammonium compounds, and the anionic and nonionic form of a fatty acid molecule interacting with a Staphylococcus aureus bacterial inner membrane. The effect of the tested materials on the properties of the model membranes are evaluated with respect to various structural properties such as the lateral pressure profile, lipid tail order parameter, and the bilayer’s electrostatic potential. Conducting asymmetric loading of molecules in only one leaflet, it was observed that anionic and cationic amphiphiles have a large impact on the Staphylococcus aureus membrane’s electrostatic potential and lateral pressure profile as compared to a symmetric distribution. Nonintuitively, we find that the cationic and anionic molecules induce a similar change in the electrostatic potential, which points to the complexity of membrane interfaces, and how asymmetry can induce biophysical consequences. Finally, we link changes in membrane structure to the rate of electroporation for the membranes, and again find a crucial impact of introducing asymmetry to the system. Understanding these physical mechanisms provides critical insights and viable pathways for the rational design of membrane-active molecules, where controlling the localization is key.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Columnar Grain-driven Plasticity and Cracking in Nanotwinned FCC Metals

The mechanisms of strengthening and plasticity in columnar-grained metals with preferentially oriented nano-sized twins have been examined traditionally by considering dislocation processes, but rarely from the perspective of grain boundary (GB) deformation. Here, the effects of GB strain accommodation on plastic deformation in four different columnar-grained nanocrystalline nanotwinned (nt) face-centered-cubic metals (Cu, Ag, Al, and Ni) were studied by large-scale molecular dynamics simulations. It is observed that in tensile deformation parallel to coherent twin boundaries (CTBs), the dislocation mechanisms in each metal are identical and associated with GB emissions of jog and threading dislocations at small and large CTB spacings, respectively. However, CTB strengthening effects are increasingly more pronounced in columnar-grained nt metals as their shear modulus increases, which is rationalized by the dependence of GB stress concentrations on twin size, metal type and strain rate. Also, while flow stresses in nt-Cu, nt-Ag, and nt-Al metals increase linearly with decreasing CTB spacing, a maximum strength limit is reached in nt-Ni below a critical CTB spacing of 6 nm. The strength limit in nt-Ni results from columnar GB cracking induced by prominent GB sliding. For columnar-grained microstructures, GB sliding is equivalent in nt-Ag and nt-Al and slightly lower in nt-Cu but markedly higher in nt-Ni. These findings underscore the importance of new GB deformation mechanisms on plasticity and fracture in columnar-grained nt metals and enrich our understanding of CTB strengthening in fcc metals synthesized in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intrinsic and extrinsic effects on the fracture toughness of ductile metallic glasses

The potential application of metallic glasses (MGs) as structural materials demands accurate measurements of their toughness and the understanding of the underlying factors affecting it. Currently, it is challenging to precisely measure the toughness of MGs, especially ductile MGs. Measured toughness values are widely scattered, even for MGs with identical compositions. That is attributed to the combined effect of intrinsic and extrinsic factors including processing, sample geometries, and loading conditions. A fundamental understanding of the influences of these elements is thus of great significance. In the present study, molecular dynamics simulations are performed to investigate the influence of intrinsic and extrinsic effects on the fracture toughness of CuZr MGs. In particular, focused is placed on the effects of cooling rate, notch shape and size, strain rate, and temperature. The results indicate that the fracture toughness of a MG scales with the cooling rate used to prepare it. This is attributed to increased free volume content generated at high cooling rates, which enhances plastic deformation and amplifies the associated energy dissipation during plastic deformation events. Fracture toughness also increases with the strain rate, arguably due to strain rate-induced crack extension delay. Overall, the results demonstrate that the largest fracture toughness are achieved when MG samples are fabricated at high cooling rates and subjected to high strain rate deformation. In addition, results suggest that the fracture toughness decreases with increasing temperature, due to the significant decrease in strength. Here, the correlations revealed between these crucial intrinsic and extrinsic parameters and the calculated MG fracture toughness support the development of a framework to understand the root of the discrepancies in the measurement of the toughness of MGs and provide insights into the design of tough MGs for structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and post-simulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. Here, the program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn’s simulation capability and analysis module, we report complete active space self-consistent field–based AIMS simulations of the 1,2- dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗