Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Atomic Edge-Guided Polyethylene Crystallization on Monolayer Two-Dimensional Materials

In this report we combine an advanced synthesis of two-dimensional (2D) materials (MoSe 2 ) having well-defined atomic edge configurations with ab initio and atomistic molecular dynamics (MD) simulations to study how atomic edges interact with polyethylene (HDPE) chains in a dilute solution assembly process. Our results reveal that Mo-terminated zigzag (Mo-ZZ) edges act as preferred nucleation sites and strongly interact with HDPE chains. The HDPE chains align in parallel with the Mo-ZZ edges and form arrays of lamellae that are perpendicular to the edges. Interestingly, atomic edge configurations are observed to dramatically change such interactions. The crystallization discrepancy at different edges was demonstrated on the same piece of MoSe 2 with different types of edges. The ab initio and MD simulations between n-alkane (n = 5 and 25), a segment of HDPE, and MoSe 2 suggest that the atomic structures of MoSe 2 can affect their interactions with n-alkane chains. Following the Mo-ZZ edge preferred nucleation principle, controlled long-range alignment of HDPE lamellae can be realized by creating multilayer MoSe 2 with parallel atomic steps. This research opens a pathway toward an atomic level understanding of polymer–2D nanomaterial interactions. It also bridges the gap between atomic-level and long-range mesoscopic structures and introduces a novel strategy for long-range structural control.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of Ionic Group Distribution on the Structure and Dynamics of Amorphous Polymer Melts

Ionizable groups tethered to a polymer backbone often associate to form ionic clusters, whose features are determined by the balance of the polymer backbone and electrostatics. These assemblies impact the dynamics of the macromolecules and their ability to transport ions. Here, using fully atomistic molecular dynamics (MD) simulations, we investigate the effects of the distribution of ionizable groups along the polymer backbone on cluster characteristics and the resulting impacts on the structure and dynamics of amorphous polymers. Particularly, we probe polystyrene sulfonates (PSS) with random, precise, and block configurations of the SO 3 – sulfonate groups along the backbone with Na + as the counterion. We find that the distribution of the ionic groups affects the shape and distribution of the clusters as well as the internal packing of the ionizable groups in the cluster and the number of unique chains that participate in each cluster, affecting the structure and the dynamics of the polymers. Furthermore, the signature of ionic clusters, observed in the static structure factor S(q) for all three distributions, is significantly more pronounced for the precise and blocky polymers compared to the random one. Remarkably, we find that the local mobility of the polymer segments is not only affected by the number and size of the clusters but also by the number of polymer chains associated with clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theory of Cation Solvation in the Helmholtz Layer of Li-Ion Battery Electrolytes

The solvation environments of Li + in conventional nonaqueous battery electrolytes, such as LiPF 6 in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize transport properties and solid electrolyte interphase (SEI) formation. Solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here, we develop and test a theory of cation solvation in the Helmholtz layer of nonaqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF 6 EC/EMC mixtures as a function of surface charge, where we find the theory can qualitatively capture the solvation environments. Next, we turn to the Helmholtz layer, where we find the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li + and the solvents, again where we find reasonable agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li + binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding battery electrolytes. Our developed formalism can be parametrized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer through reducing the number of available coordination sites, which can be used to determine what could react and form the SEI.

Helmholtz↗
DOE OSTI2025

Ionic Associations and Hydration in the Electrical Double Layer of Water-in-Salt Electrolytes

Water-in-Salt-Electrolytes (WiSEs) are an exciting class of concentrated electrolytes finding applications in energy storage devices because of their expanded electrochemical stability window, good conductivity and cation transference number, and fire-extinguishing properties. These distinct properties are thought to originate from the presence of an anion-dominated ionic network and interpenetrating water channels for cation transport, which indicates that associations in WiSEs are crucial to understanding their properties. Currently, associations have mainly been investigated in the bulk, while little attention has been given to the electrolyte structure near electrified interfaces. Here, we develop a theory for the electrical double layer (EDL) of WiSEs, where we consistently account for the thermoreversible associations of species into Cayley tree aggregates. The theory predicts an asymmetric structure of the EDL. At negative voltages, hydrated Li + dominates, and cluster aggregation is initially slightly enhanced before disintegration at larger voltages. At positive voltages, when compared to the bulk, clusters are strictly diminished. Performing atomistic molecular dynamics (MD) simulations of the EDL of WiSE provides EDL data for validation and bulk data for parametrization of our theory. Validating the predictions of our theory against MD showed good qualitative agreement. Furthermore, we performed electrochemical impedance measurements to determine the differential capacitance of the studied LiTFSI WiSE and also found reasonable agreement with our theory. Overall, the developed approach can be used to investigate ionic aggregation and solvation effects in the EDL, which, among other properties, can be used to understand the precursors for solid-electrolyte interphase formation.

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

Morphology and Dynamics in Hydroxide-Conducting Polysulfones

In alkaline zinc–manganese dioxide batteries, there is a need for selective polymeric separators that have good hydroxide ion conductivity but that prevent the transport of zincate (Zn(OH) 4 ) 2- . Here we investigate the nanoscale structure and hydroxide transport in two cationic polysulfones that are promising for these separators. We present the synthesis and characterization for a tetraethylammonium-functionalized polysulfone (TEA-PSU) and compare it to our previous work on an N-butylimidazolium-functionalized polysulfone (NBI-PSU). We perform atomistic molecular dynamics (MD) simulations of both polymers at experimentally relevant water contents. The MD simulations show that both polymers develop well phase separated nanoscale water domains that percolate through the polymer. Calculation of the total scattering intensity from the MD simulations reveal weak or nonexistent ionomer peaks at low wave vectors. The lack of an ionomer peak is due to a loss of contrast in the scattering. The small water domains in both polymers, with median diameters on the order of 0.5–0.7 nm, lead to hydroxide and water diffusion constants that are 1–2 orders of magnitude smaller than their values in bulk water. This confinement lowers the conductivity but also may explain the strong exclusion of zincate from the PSU membranes seen experimentally.

25 ENERGY STORAGE↗
DOE OSTI2021

Molecular Level Differences in Ionic Solvation and Transport Behavior in Ethylene Oxide-Based Homopolymer and Block Copolymer Electrolytes

Block copolymer electrolytes (BCE) such as polystyrene-block-poly(ethylene oxide) (SEO) blended with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and composed of mechanically robust insulating and rubbery conducting nanodomains are promising solid-state electrolytes for Li batteries. Here, we compare ionic solvation, association, distribution, and conductivity in SEO-LiTFSI BCEs and their homopolymer PEO-LiTFSI analogs toward a fundamental understanding of the maximum in conductivity and transport mechanisms as a function of salt concentration. Ionic conductivity measurements reveal that SEO-LiTFSI and PEO-LiTFSI exhibit similar behaviors up to a Li/EO ratio of 1/12, where roughly half of the available solvation sites in the system are filled, and conductivity is maximized. As the Li/EO ratios increase to 1/5 the conductivity, of the PEO-LiTFSI drops nearly 3-fold, while the conductivity of SEO-LiTFSI remains constant. FTIR spectroscopy reveals that additional Li cations in the homopolymer electrolyte are complexed by additional EO units when the Li/EO ratio exceeds 1/12, while in the BCE, the proportion of complexed and uncomplexed EO units remains constant; Raman spectroscopy data at the same concentrations show that Li cations in the SEO-LiTFSI samples tend to coordinate more to their counteranions. Atomistic-scale molecular dynamics simulations corroborate these results and further show that associated ions tend to segregate to the SEO-LiTFSI domain interfaces. The opportunity for "excess" salt to be sequestered at BCE interfaces results in the retention of an optimum ratio of uncompleted and complexed PEO solvation sites in the middle of the conductive nanodomains of the BCE and maximized conductivity over a broad range of salt concentrations.

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

Evaluating the Role of Metastable Surfaces in Mechanochemical Reduction of Molybdenum Oxide

Mechanochemistry and mechanocatalysis are gaining increasing attention as environmentally friendly chemical processes because of their solvent-free nature and scalability. Significant effort has been devoted for studying continuum-scale phenomena in mechanochemistry, such as temperature and pressure gradients, but the atomic-scale mechanisms remain relatively unexplored. In this work, we focus on the mechanochemical reduction of MoO3 as a case study. We use experimental techniques to determine the mechanochemical reduction conditions and density functional theory (DFT) simulations to establish an atomistic framework for identifying the metastable surfaces that are most likely to enable this process. Our results show that metastable surfaces can significantly lower or remove thermodynamic barriers for surface reduction and that kinetic energy from milling can facilitate the formation of metastable surfaces that have high surface fracture energies and are not thermally accessible. These findings indicate that metastable surfaces are an important aspect of mechanochemistry along with hot spots and other continuum-scale phenomena.

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

Stability and molecular pathways to the formation of spin defects in silicon carbide

Abstract Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate—the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon (V Si ) to carbon monovacancies (V C ), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of V Si relative to V C favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Engineering of Human Apolipoprotein for Enhanced Interaction with Bacterial Toxins

Sepsis is caused by an over-activation of the immune system in response to the presence of bacterial endotoxins in the bloodstream. LPS exemplifies a potent immune stimulator which can lead to an excessive activation leading to an overly pro-inflammatory reaction and if left untreated results in multi-organ failure and death. Human lipoproteins (HLP) have been shown to neutralize LPS and facilitate its clearance from circulation via the liver. However, host-derived lipoproteins are inhomogeneous in content and size, presenting a challenge for therapeutic applications. In this work, we have used IC facilities in order to computationally engineer HLP for an enhanced interaction with LPS. This work is a continuation of our previous efforts to derive an optimal Nanodisc (ND) lipid mixture for enhanced association with the endotoxin. We have performed fully atomistic Molecular dynamic (MD) simulations in order to carefully identify the most optimal residues in HLP, which can potentially lead to a stronger association to LPS.

42 ENGINEERING↗
DOE OSTI2023

Allosterically coupled conformational dynamics in solution prepare the sterol transfer protein StarD4 to release its cargo upon interaction with target membranes

Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1984

An atomistic model of slip formation

The results of an atomistic model for the simulation of the early stages of crack initiation in a two-dimensional triangular lattice are presented. In the current model, each particle in the system is treated discretely and assumed to be interacting with the surrounding particles via Lennard-Jones potentials. A uniaxial load (in incremental elongations) is applied to the rectangular two-dimensional slab in either the x or the y direction. After each incremental elongation the system is equilibrated using a static method. Initially, elastic behavior in the x and y directions is observed. Continued elongation results in plastic deformation. In lattices with point defects, the defects first move to the surface, creating vacancies which trigger plastic deformation.

Halicioglu, T.↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

First-principles effective Hamiltonian for finite-temperature modeling of nonperovskite ferroelectrics

First-principles-based effective Hamiltonian techniques have been widely employed for over three decades to investigate ferroelectricity and related phenomena in perovskite materials. These techniques offer high accuracy, transferability, compatibility with various finite-temperature algorithms, computational efficiency, and ease in incorporating interactions with external fields. They have been adapted to study diverse phenomena, ranging from topological dipole patterns in ferroelectric nanostructures to multicaloric effects. In this work, we develop an effective Hamiltonian for the nonperovskite ferroelectric HfO 2 (hafnia). Applying this methodology to explore the finite-temperature and finite-electric-field properties of ferroelectric hafnia revealed (1) exceptionally large intrinsic coercive fields, an order of magnitude higher than those observed in perovskite ferroelectrics; (2) their atomistic origin; and (3) the existence of a regime where the relationship between the coercive field and the energy barrier for polarization reversal is counterintuitive. Here, these developments could accelerate progress both in methodological advancements for simulating ferroics and in the atomistic understanding of a broad range of ferroelectrics.

Electric polarization↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
DOE OSTI2021

1D EDIPIC

EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.

Sydorenko, Dmytro↗