Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

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 451 records · Page 25

Multi-scale Rule-of-Mixtures Model of Carbon Nanotube/Carbon Fiber/Epoxy Lamina

A unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes is modeled with a multi-scale method, the atomistically informed rule-of-mixtures. This multi-scale model is designed to include the effect of the carbon nanotubes on the constitutive properties of the lamina. It included concepts from the molecular dynamics/equivalent continuum methods, micromechanics, and the strength of materials. Within the model both the nanotube volume fraction and nanotube distribution were varied. It was found that for a lamina with 60% carbon fiber volume fraction, the Young's modulus in the fiber direction varied with changes in the nanotube distribution, from 138.8 to 140 GPa with nanotube volume fractions ranging from 0.0001 to 0.0125. The presence of nanotube near the surface of the carbon fiber is therefore expected to have a small, but positive, effect on the constitutive properties of the lamina.

Frankland, Sarah-Jane V.↗

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Hot Structure Concepts and Materials for Space Exploration Project

Hot structures are multifunctional material systems that integrate the functions of thermal protection systems (TPS) and primary structure to increase mass efficiency in aerospace vehicles. The goal of this project was to improve our capabilities for developing lightweight, reusable, and reliable hot structures concepts to enable advanced space exploration missions. Specifically, current hot structure carbon/carbon (C/C) composite material systems have low interlaminar strength, which can lead to material failures at low load levels, and large material property variability in as-fabricated parts. Two parallel paths were followed to address these shortcomings. First, analysis models were developed to identify the conditions under which excessive interlaminar stresses occur and to evaluate the potential of stitching the C/C laminae to improve the interlaminar thermal-structural response. This effort included the development of high-fidelity, fracture-based thermal-structural finite element analysis models to evaluate critical features needed to produce efficient hot structure concepts. The second path focused on developing atomistic and mesoscale simulation tools to model the fabrication process in an effort to reduce the time and cost of production while improving the quality of the resulting parts. Aspects of the process currently being investigated include the infusion of the matrix precursor into the fiber preform, the cure of the precursor molecules into their polymeric form, and the carbonization and graphitization of the polymer to form the carbon matrix material. Results of the 6-month study are presented herein.

thermal↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

Highly transferable atomistic machine-learning potentials from curated and compact datasets across the periodic table

Machine learning atomistic potentials trained using density functional theory (DFT) datasets allow for the modeling of complex material properties with near-DFT accuracy while imposing a fraction of its computational cost. The curation of the DFT datasets can be extensive in size and time-consuming to train and refine. In this study, we focus on addressing these barriers by developing minimalistic and flexible datasets for many elements in the periodic table regardless of their mass, electronic configuration, and ground state lattice. These DFT datasets have on average, ~4000 different structures and 27 atoms per structure, which we found sufficient to maintain the predictive accuracy of DFT properties and notably with high transferability. We envision these highly curated training sets as starting points for the community to expand, modify, or use with other machine learning atomistic potential models, whatever may suit individual needs, further accelerating the utilization of machine learning as a tool for material design and discovery.

42 ENGINEERING↗

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

59 BASIC BIOLOGICAL SCIENCES↗

hippynn Python Package

hippynn is a python package for defining, training, and applying neural networks to atomistic systems. In particular, it focuses on Hierarchical Interacting Particle Neural Networks (HIP-NNs), a deep learning architecture for atomistic systems. HIP-NNs take input data describing the properties of atomistic systems (often generated using ab-initio quantum mechanics) to learn fast and accurate models. A trained HIP-NN can predict potential energy surfaces, atomic charges, and more. hippynn uses PyTorch for portable high-performance code, including both CPU and GPU support. hippynn allows for extensive customization, including user-defined models and loss functions, to facilitate future research into extensions of HIP-NN as well as other atomistic deep learning models.

Lubbers, Nicolas↗

Experimental characterization and atomistic simulation of grain boundary segregation in Mg-Y alloys

As a rare earth solute element in Mg alloys, Y has the beneficial effects of increasing both the strength and the ductility as well as weakening the crystallographic texture. To achieve a more fundamental understanding on how Y addition affects the microstructural evolution and mechanical properties, the Y segregation behavior at grain boundaries was investigated in Mg-1wt.%Y and Mg-7wt.%Y alloys at different conditions. The segregation intensity and its dependence on the grain boundary misorientation angle were experimentally characterized and computationally predicted. Strong segregation at grain boundaries was observed in both low and high Y-containing alloys. Y segregation was found to remain in alloy Mg-7Y after high-temperature annealing heat treatment at 540 °C. No direct correlation between the Y segregation intensity and the grain boundary misorientation angle could be established based on either the experimental characterization or the atomistic simulation with a spectral model. We thus conclude that grain boundary segregation of Y is independent of grain boundary misorientation angle.

Grain boundary↗

Cation–Ligand Interactions Dictate Salt Partitioning and Diffusivity in Ligand-Functionalized Polymer Membranes

Membranes are an attractive alternative to current thermal separations due to their scalability and energy efficiency in desalinating water. Unfortunately, many of the conventional membrane materials available today are unable to differentiate between ionic solutes, especially alkali cations, compromising their use in ion–ion separations. Inspired by the ion-specific interactions exhibited by biological ion channels, recent research efforts have focused on synthesizing and characterizing new polymeric materials that incorporate ligands into polymer networks to bias solubility and/or diffusivity of one cationic species over another. Despite these efforts, little is known about the influence of incorporating ligands into polymer membranes on solubility and diffusivity of the complexing species. In this study, we first build a qualitative model of salt partitioning, diffusivity, and permeability in generic cation-complexing ligand-functionalized polymer membranes. Next, to validate our model and hypotheses, we perform atomistic molecular dynamics simulations of a 12-crown-4-functionalized membrane in the presence of alkali halide salts at low concentration. Generally, cation complexation enhances cation solubility but decreases diffusivity. Interestingly, the reduction in diffusivity is predicted to be larger than the enhancement in solubility for materials which operate by the mechanisms proposed in our physical picture, ultimately resulting in a reduction in the permeability of the selectively complexing ion.

36 MATERIALS SCIENCE↗

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗

1 / ω electric-field noise in surface ion traps from correlated adsorbate dynamics

Ion traps are promising architectures for implementing quantum computers, but they suffer from excessive “anomalous” ion motional heating that limit their overall coherence and practicality for scalable quantum computing. The exact microscopic origins of anomalous heating remain an open question, but experiments point to adsorbates on trap electrodes as one likely source. Many different models of anomalous heating have been proposed, but these models have yet to pinpoint the atomistic origin of the experimentally observed 1/ω electric-field noise scaling seen in ion traps at frequencies between 0.1–10 MHz. In this work, we show that a model based on previously proposed surface-induced dipole fluctuations on adsorbates, but which also incorporates interparticle interaction dynamics through molecular dynamics simulations of up to multiple monolayers of adsorbates, gives rise to 1/ω frequency scaling at the MHz frequencies typically employed in ion traps. Overall, these results demonstrate that moderate-to-high densities of adsorbates can give rise to a set of activated motions that produce the 1/ω noise observed in ion traps and that collective adsorbate motions produce the observed noise spectra that a noninteracting model does not capture.

36 MATERIALS SCIENCE↗

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theories of homogeneous and electrochemical electron transfer in complex media and interfaces (Final Technical Report)

This project makes the next step in establishing practical theories of charge transfer in complex media. The development of formal models is supported by extensive atomistic simulations, quantum calculations of force-field parameters, and direct measurements of charge-transfer spectra. All theory development is supported by experiment, extensive numerical simulations, and through external collaborations.

14 SOLAR ENERGY↗

Fused Deposition Modeling Additive Manufacturing of Carbonized Structures via Waste-Enhanced Filaments

This report details the development, characterization, and use of coal-enhanced composite materials in additive manufacturing applications. High coal loading formulations—containing up to 70 wt.% coal—were successfully extruded and processed using commercially available 3D printers. Extensive experimental testing was conducted to assess the mechanical, thermal, and microstructural properties of the composites. In parallel, multi-scale computational modeling was employed to elucidate atomistic interactions and evaluate the effects of printing-induced defects on structural performance. Large-scale printability trials demonstrated the feasibility of fabricating complex components for tooling and construction applications, including wind turbine blade molds and modular wall sections. Techno-economic analyses demonstrated the cost-effectiveness and scalability of coal-enhanced composites for large-scale additive manufacturing applications such as wind turbine blade tooling.

01 COAL, LIGNITE, AND PEAT↗