Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular dynamic simulations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Subdiffusive High-Pressure Hydrogen Gas Dynamics in Elastomers

We report elastomeric rubber materials serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. High-pressure exposure and subsequent rapid decompression often lead to cavitation and stress-induced damage of the elastomer due to localization of the hydrogen gas. Here, we use all-atom classical molecular dynamics simulations to assess the impact of compositional variations on gas diffusion within the commonly used elastomer ethylene–propylene–diene monomer (EPDM). With the aim to build a predictive understanding of precursors to cavitation and to motivate material formulations that are less sensitive to hydrogen-induced failure, we perform systematic simulations of gas dynamics in EPDM as a function of temperature, gas concentration, and cross-link density. Our simulations reveal anomalous, subdiffusive hydrogen motion at pressure and intermediate times. We identify two groups of gas with different mobilities: one group exhibiting high mobility and one group exhibiting low mobility due to their motion being impeded by the polymer. With decreasing temperatures, the low-mobility group shows increased gas localization, the necessary precursor for cavitation damage in these materials. At lower temperatures, increasing cross-link density led to greater hydrogen gas mobility and a lower fraction of caged hydrogen, indicating that increasing cross-link density may reduce precursors to cavitation. Finally, we use a two-state kinetic model to determine the energetics associated with transitions between these two mobility states.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electric-Field-Induced Assembly of an Ionic Liquid–Water Interphase Enables Efficient Heavy Metal Electrosorption

Controlling ion desolvation, transport, and charge transfer at the electrode-electrolyte interface (EEI) is critical to enabling the rational design of efficient and selective separation of targeted heavy metals and decontamination of industrial waste water. One challenge is to sufficiently resolve and interrogate the intermediate transformation steps between solvated metal cations and bare cations for electroreduction at the EEI and establish pathways to modulate these steps to achieve efficient energy transfer for targeted reactive separations. Herein, we obtained a predictive understanding of the effect of adding a hydrophobic ionic liquid, 1-ethyl-3-methylimidazolium chloride (EMIMCl), to an aqueous electrolyte on modulating the desolvation and electrosorption of Pb 2+ cations using a combination of experimental and theoretical techniques. Here, we reveal formation of a compact interphase layer consisting of EMIMCl-Pb complex clusters under an applied electric field using operando electrochemical Raman spectroscopy, atomic force microscopy, and electrochemical impedance spectroscopy measurements combined with classical molecular dynamics simulations. The application of a larger negative potential is shown to result in formation of a well oriented layer with the positive imidazolium ring of the EMIMCl lying on the electrode and the long hydrophobic alkyl chain extending into the bulk electrolyte. This oriented layer is demonstrated to facilitate desolvation of incoming solvated Pb 2+ cations and decrease the charge transfer resistance for Pb electrodeposition, which has important implications for selective removal of Pb from contaminated mixtures. Overall, our studies open up new opportunities to modulate ion desolvation using hydrophobic ionic liquids in aqueous electrolytes for efficient heavy metal separation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tuning the Functionalities of Porous Liquids for Emergent Gas-Capture Properties

Type 3 Porous Liquids (PLs) are a class of materials with the potential to revolutionize gas capture, storage, and utilization. These PLs are formed by suspending sorbent nanoparticles (e.g., metal−organic frameworks) in sterically excluded solvents, creating permanent porosity for gas capture in a processable, low-viscosity phase. Herein, a computational study revealed sorbent surface functionalization strategies to enhance CO 2 sorption, and the molecular structural signatures underpinning the enhancements in gas uptake. PLs composed of a ZIF-8 surface functionalized with 3-amino-1,2,4-triazole (Atz) in glyceryl triacetate were targeted for emergent CO 2 capture, exceeding that of the unfunctionalized ZIF-8 PL. ZIF-8 was surface functionalized with Atz at various surface coverage fractions (f), and classical molecular dynamics simulations predicted an increase in CO 2 sorption capacity with increasing f, up to f = 0.75. Additionally, detailed structural analyses revealed that solvent orientational order, derived from the solvent triplet-angle distribution, can identify the gas-capture potential of a PL without requiring computationally expensive direct modeling of the CO 2 sorption. Combined with experimental validation, initial computational screening of PL compositions promises to accelerate the discovery of PL compositions for novel gas separation materials platforms.

gas capture↗
DOE OSTI2021

Water or Anion? Uncovering the Zn 2+ Solvation Environment in Mixed Zn(TFSI) 2 and LiTFSI Water-in-Salt Electrolytes

Applications of aqueous zinc batteries for grid-scale energy storage are limited by their poor reversibility and the competing water splitting reaction. The recent invention of a water-in-salt (WIS) electrolyte concept provides a new route enabling a stable and highly reversible aqueous zinc battery chemistry. In the present work, a mixed zinc bis(trifluoromethane sulfonyl)imide (Zn(TFSI) 2 ) and LiTFSI WIS electrolyte was studied using X-ray total scattering, X-ray absorption, and Fourier transform infrared spectroscopy in conjunction with classical molecular dynamics simulations. It was found that, in the highly concentrated WIS electrolyte consisting of 1 m Zn(TFSI) 2 and 20 m LiTFSI, Zn 2+ cations are mainly solvated by six waters in their first solvation shell, while the TFSI – anions are completely excluded. This ion solvation picture is fundamentally different from the previous understandings. The results suggest that additional studies are needed to fully understand the unusual stability and reversibility of zinc-WIS electrolyte-based batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria that encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo , the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We used molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo . Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed the overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceu m model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

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

Revitalizing Iron Redox by Anion-Insertion-Assisted Ferro- and Ferri-Hydroxides Conversion at Low Alkalinity

Iron hydroxides are desirable alkaline battery electrodes for low cost and environmental beneficence. However, hydrogen evolution on charging and Fe 3 O 4 formation on discharging cause low storage capacity and poor cycling life. Here, we report that green rust (GR) (Fe 2+ 4 Fe 3+ 2 (HO – ) 12 SO 4 ), formed via sulfate insertion, promotes Fe(OH) 2 /FeOOH conversion and shows a discharge capacity of ~211 mAh g –1 in half-cells and Coulombic efficiency of 93% after 300 cycles in full-cells. Theoretical calculations show that Fe(OH) 2 /FeOOH conversion is facilitated by intercalated sulfate anions. Classical molecular dynamics simulations reveal that electrolyte alkalinity strongly impacts the energetics of sulfate solvation, and low alkalinity ensures fast transport of sulfate ions. Anion-insertion-assisted Fe(OH) 2 /FeOOH conversion, also achieved with Cl – ion, paves a pathway toward efficient utilization of Fe-based electrodes for sustainable applications.

25 ENERGY STORAGE↗
DOE OSTI2023

The Energetic Origins of Pi–Pi Contacts in Proteins

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein’s stability and structure. Pi–pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi–pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi–pi contact motifs. We find that neutral pi–pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter–Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi–pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi–pi contacts in polar media and like typical ion–pi interactions in nonpolar environments. Like-charged pairs such as arginine–arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi–pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

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

Insights into solvent and surface charge effects on Volmer step kinetics on Pt (111)

Abstract The mechanism of pH-dependent hydrogen oxidation and evolution kinetics is still a matter of significant debate. To make progress, we study the Volmer step kinetics on platinum (111) using classical molecular dynamics simulations with an embedded Anderson-Newns Hamiltonian for the redox process and constant potential electrodes. We investigate how negative electrode electrostatic potential affects Volmer step kinetics. We find that the redox solvent reorganization energy is insensitive to changes in interfacial field strength. The negatively charged surface attracts adsorbed H as well as H + , increasing hydrogen binding energy, but also trapping H + in the double layer. While more negative electrostatic potential in the double layer accelerates the oxidation charge transfer, it becomes difficult for the proton to move to the bulk. Conversely, reduction becomes more difficult because the transition state occurs farther from equilibrium solvation polarization. Our results help to clarify how the charged surface plays a role in hydrogen electrocatalysis kinetics.

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

High-throughput screening of hypothetical metal-organic frameworks for thermal conductivity

Thermal energy management in metal-organic frameworks (MOFs) is an important, yet often neglected, challenge for many adsorption-based applications such as gas storage and separations. Despite its importance, there is insufficient understanding of the structure-property relationships governing thermal transport in MOFs. To provide a data-driven perspective into these relationships, here we perform large-scale computational screening of thermal conductivity k in MOFs, leveraging classical molecular dynamics simulations and 10,194 hypothetical MOFs created using the ToBaCCo 3.0 code. We found that high thermal conductivity in MOFs is favored by high densities (> 1.0 g cm -3 ), small pores (< 10 Å), and four-connected metal nodes. We also found that 36 MOFs exhibit ultra-low thermal conductivity (< 0.02 W m -1 K -1 ), which is primarily due to having extremely large pores (~65 Å). Furthermore, we discovered six hypothetical MOFs with very high thermal conductivity (>10 Wm -1 K -1 ), the structures of which we describe in additional detail.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights on the dissolution of water in an albite melt at high pressures and temperatures from a direct structural analysis

The water dissolution mechanism in silicate melts under high pressures is not well understood. Here we present the first direct structure investigation of a water-saturated albite melt to monitor the interactions between water and the network structure of silicate melt at the molecular level. In situ high-energy X-ray diffraction was carried out on the NaAlSi 3 O 8 -H 2 O system at 800 °C and 300 MPa, at the Advanced Photon Source synchrotron facility. The analysis of the X-ray diffraction data was augmented with classical Molecular Dynamics simulations of a hydrous albite melt, incorporating accurate water-based interactions. The results show that metal–oxygen bond breaking at the bridging sites occurs overwhelmingly at the Si site upon reaction with H 2 O, with subsequent Si–OH bond formation and negligible Al–OH formation. Furthermore, we see no evidence for the dissociation of the Al 3+ ion from the network structure upon breaking of the Si–O bond in the hydrous albite melt. The results also indicate that the Na + ion is an active participant in the modifications of the silicate network structure of the albite melt upon water dissolution at high P–T conditions. We do not find evidence for the Na + ion dissociating from the network structure upon depolymerization and subsequent formation of NaOH complexes. Instead, our results show that the Na + ion persists as a structure modifier with a shift away from Na–BO bonding to an increase in the extent of Na-NBO bonding, in parallel with pronounced depolymerization of the network. Our MD simulations show that the Si–O and Al–O bond lengths are expanded by about 6% in the hydrous albite melt compared to those of the dry melt at high P–T conditions. The changes in the network silicate structure of a hydrous albite melt at high pressure and temperature, as revealed in this study, must be considered in the advancement of water dissolution models of hydrous granitic (or alkali aluminosilicate) melts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reversible assembly of silica nanoparticles at water–hydrocarbon interfaces controlled by SDS surfactant

Achieving reversible and tunable assembly of silica nanoparticles at liquid–liquid interfaces is vital for a wide range of scientific and technological applications including sustainable subsurface energy applications, catalysis, drug delivery and material synthesis. In this study, we report the mechanisms controlling the assembly of silica nanoparticles (dia. 50 nm and 100 nm) at water–heptane and water–toluene interfaces using sodium dodecyl sulfate (SDS) surfactant with concentrations ranging from 0.001–0.1 wt% using operando ultrasmall/small-angle X-ray scattering, cryogenic scanning electron microscopy imaging and classical molecular dynamics simulations. The results show that the assembly of silica nanoparticles at water–hydrocarbon interfaces can be tuned by controlling the concentrations of SDS. Silica nanoparticles are found to: (a) dominate the interfaces in the absence of interfacial SDS molecules, (b) coexist with SDS at the interfaces at low surfactant concentration of 0.001 wt% and (c) migrate toward the aqueous phase at a high SDS concentration of 0.1 wt%. Energetic analyses suggest that the van der Waals and electrostatic interactions between silica nanoparticles and SDS surfactants increase with SDS concentration. However, the favorable van der Waals and electrostatic interactions between the silica nanoparticles and toluene or heptane decrease with increasing SDS concentration. As a result, the silica nanoparticles migrate away from the water–hydrocarbon interface and towards bulk water at higher SDS concentrations. These calibrated investigations reveal the mechanistic basis for tuning silica nanoparticle assembly at complex interfaces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Exemption of lattice collapse in Ni–MnO 2 birnessite regulated by the structural water mobility

Lattice collapse and associated mechanical fracture frequently occur in Li-intercalated metal oxide cathodes at the deep charge state upon Li-ion removal, governed by chemical compositions and the resulting electron density of the oxygen atoms. However, similar lattice collapse for metal oxide electrodes in aqueous storage and its mitigation have not been well studied. In this work, we reported the lattice collapse of MnO 2 layered birnessite during the aqueous de-sodiation process at high voltage due to the structural water motion, as evidenced by in situ XRD. Unlike non-aqueous Li-intercalated electrodes, Ni-dopants mitigated the lattice collapse of birnessite at deep charge states. Moreover, density functional theory (DFT) calculations showed that Ni doping induces charge depletion of lattice oxygen due to its higher electronegativity than Mn. This charge reduction, in turn, yields a significant decrease in electrostatic repulsion between the oxygens belonging to lattice and structural water. Classical molecular dynamics simulations based on atomic charges obtained from DFT elucidate that Ni-doping facilitates immobilization of structural water in (Ni)MnO 2 and prevents lattice collapse upon Na-ion removal. (Ni)MnO 2 exempted from lattice collapse shows an improved storage capacity relative to MnO 2 while maintaining similar cycling stability.

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

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

DEHP− extractant binding to trivalent lanthanide Er3+: Fast binding accompanied by concerted angular motions of hydration water

Solvent extraction of trivalent rare earth metal ions by organophosphorus extractants proceeds via binding of phosphoric acid headgroups to the metal ion. Water molecules in the tightly bound first hydration shell of the metal ions must be displaced by oxygen atoms from phosphoric acid headgroups. Here, we use classical molecular dynamics simulations to explore the event in which a fully hydrated Er3+ binds to its first phosphoric acid headgroup. Approach of the headgroup into the region between the first and second hydration shells leads to a fast ejection of a water molecule that is accompanied by reordering of the hydration water molecules, including discretization of their angular positions and collective rotation about the metal ion. The water molecule ejected from the first shell is located diametrically opposite from the binding oxygen. Headgroup binding places a headgroup oxygen closer to Er3+ than its first hydration shell and creates a loosely bound water that subsequently exchanges between the first shell and its environment. This second exchange of water also occurs at discrete angular positions. This geometrical aspect of binding may be of relevance to understanding the binding and transport of ion–extractant complexes that are expected to occur at the organic–aqueous liquid–liquid interface used in solvent extraction processes.

Chemistry↗
DOE OSTI2023

Molecular dynamics study on the interfacial properties of mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders on graphene and graphite surfaces

This study investigates the behavior of two different mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders, polyvinylpyrrolidone:polyvinylidene difluoride (PVP:PVDF) and polyvinylpyrrolidone:polyacrylic acid (PVP:PAA), at graphene and graphite interfaces using classical molecular dynamics simulations. The aim is to identify the best performing monomer binder blend and carbon-based material for the design of battery-optimized energy devices. The PVP:PAA monomer binder blend and graphite are found to have the best interaction energies, densification upon adsorption, and more ordered structure. The adsorption of both monomer binder blends is strongly guided by the higher affinity of PVP and PAA monomeric molecules for the surfaces compared to PVDF. The structure of adsorbed layers of PVP:PVDF monomer binder blend on graphene and graphite develops more quickly than PVP:PAA, indicating faster kinetics. This study complements a previous density functional theory study recently reported by our group and contributes to a better understanding of the nanoscopic features of relevant interfacial regions involving mixtures of monomers of PVP-based battery binders and different carbon-based materials. In conclusion, the effect of a blend of commonly used monomer binders on carbon-based materials is essential for obtaining tightly bound anode and cathode active materials in lithium-ion batteries, which is crucial for designing battery-optimized energy devices.

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

What do far-infrared spectra of solitary water in “water-in-solvent” systems reveal about water’s solvation and dynamics?

Classical molecular dynamics simulations of water in ionic and dipolar solvents were used to interpret the far-infrared (FIR) rotation/libration spectra of “solitary water” in terms of water’s rotational dynamics and interactions with solvents. Seven solvents represented by nonpolarizable all-atom force fields and a series of idealized variable-charge solvents were used to span the range of solvent polarities (hydrogen bonding) studied experimentally. Simulated spectra capture the solvent dependence observed, as well as the relationship between the frequencies of water libration (νL) and OH stretching bands (νOH). In more strongly interacting solvents, simulated νL are ∼20% higher than those of experiment. In all solvents, the simulated spectra are composites of rotational motions about the two axes perpendicular to water’s dipole moment, and the different frequencies of these two motions are responsible for the breadth of the libration band and the bimodal shape observed in halide ionic liquids. Simulations overestimate the separation of these two components in most solvents. The character of water rotational motions changes markedly with solvent polarity, from quasi-free rotation in nonpolar and weakly polar solvents to highly constrained libration in strongly hydrogen bonding environments. The changeover to librational motions dominating the spectrum occurs between solvents such as benzene (νL ∼ 250 cm−1) and acetonitrile (νL ∼ 400 cm−1). For solvents in the latter category, the mean frequency of the experimental FIR band provides a direct measure of mean-squared torques and, therefore, force constants associated with interactions constraining water’s librational motion.

Chemistry↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗