Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

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

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, in this work, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Effect of Molecular Structure of Quinones and Carbon Electrode Surfaces on the Interfacial Electron Transfer Process

Quinones can undergo thermodynamically reversible proton-coupled electron transfer reactions and are being applied as electroactive compounds in aqueous organic batteries. However, the electrochemical reversibility of these compounds is affected not only by their molecular structure but also by the properties of a carbon-based electrode surface. This study combines experimental and theoretical approaches to understand this dependence. Specifically, we study the electron transfer kinetics of two synthesized quinone derivatives and two commercially available ones with a glassy carbon, a highly ordered pyrolytic graphite, and a high-edge-density graphite electrode (HEDGE). The electrochemical reversibility is notably improved on the HEDGE, which shows a higher density of defects and presents oxygenated functional groups at its surface. The electron transfer kinetics are controlled by adsorbed species onto the HEDGE. Molecular dynamics simulation and quantum mechanics calculations suggest defects with oxygen-containing functional groups, such as C-O and C$=$O, on HEDGE surfaces drive the interaction with the functional groups of the molecules, during physisorption from van der Waals forces. The presence of sulfonic acid side groups and a greater number of aromatic rings in the molecular structure may contribute to a higher stabilization of quinone derivatives on HEDGEs. We propose that high-performance carbon-based electrodes can be obtained without catalysts for organic batteries, by the engineering of carbon-based surfaces with edge-like defects and oxygenated functional groups.

25 ENERGY STORAGE↗
DOE OSTI2024

The dynamic adsorption affinity of ligands is a surrogate for the passivation of surface defects

Abstract Surface defects in semiconducting materials, though they have been widely studied, remain a prominent source of loss in optoelectronic devices; here we sought a new angle of approach, looking into the dynamic roles played by surface defects under atmospheric stressors and their chemical passivants in the lifetime of optoelectronic materials. We find that surface defects possess properties distinct from those of bulk defects. ab initio molecular dynamics simulations reveal a previously overlooked reversible degradation mechanism mediated by hydrogen vacancies. We find that dynamic surface adsorption affinity (DAA) relative to surface treatment ligands is a surrogate for passivation efficacy, a more strongly-correlated feature than is the static binding strength emphasized in prior reports. This guides us to design targeted passivator ligands with high molecular polarity: for example, 4-aminobutylphosphonic acid exhibits strong DAA and provides defect passivation applicable to a range of perovskite compositions, including suppressed hydrogen vacancy formation, enhanced photovoltaic performances and operational stability in perovskite solar cells.

Science & Technology - Other Topics↗
DOE OSTI2023

Large-scale computational polymer solubility predictions and applications to dissolution-based plastic recycling

Dissolution-based plastic recycling is a promising approach to separate and recover high quality pure polymer resins from multicomponent plastic waste by exploiting differences in polymer solubility. The design of a dissolution-based polymer recycling process requires the selection of appropriate solvent systems and operating temperatures to dissolve only target polymers. Determining these parameters experimentally is challenging due to the wide range of solvents and temperatures possible for a given set of target polymers. In this work, we report a computational scheme that employs molecular dynamics simulations and the Conductor-like Screening Model for Realistic Solvents to predict polymer solubilities. Using this scheme, we established a computational solubility database for 8 common polymers and 1007 solvents at multiple temperatures and measured selected solubilities experimentally to validate computational predictions. Analysis of functional groups within this large database then provides chemical heuristics relating the molecular structures of good and non-solvents for selected polymers. We further developed a tool that automates the selection of solvents for all possible sequences in which target polymers can be selectively dissolved to guide the design of dissolution-based plastic recycling processes. Here, we demonstrate the application of these methods via multiple experimental case studies of representative dissolution-based polymer recycling processes in which pure polymer resins were successfully recovered from physical mixtures of polymers.

42 ENGINEERING↗
DOE OSTI2022

Processes controlling helium bubble dynamics at varying temperatures in simulated radioactive materials

Alpha-radiation damage is known to degrade materials exposed to harsh radiation environments via swelling and embrittlement by nano-scale helium bubbles. How these bubbles grow under varying temperature is essential to understanding and predicting material failure and optimization. Using in situ X-ray scattering, we probed ≈1012 bubbles before, during and after annealing at different temperatures in He-implanted aluminum to reveal a critical growth activation temperature that depends strongly on the helium implantation conditions. An analysis of the data aided by our molecular dynamics simulations identifies multiple temperature-dependent bubble growth mechanisms above and below the critical growth activation temperature. Moreover, the irreversible bubble growth that occurs after annealing for ≈4 h was evaluated using dislocation dynamics and found to negligibly affect the yield strength but may affect subsequent bubble growth at temperature. These in situ results provide valuable insight into bubble growth never observed before and reveal, for the first time, bubble shrinkage after cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase-field model of deformation twin-grain boundary interactions in hexagonal systems

Here, a phase-field model for describing dynamic interactions between deformation twins and grain boundaries in a hexagonal close packed (HCP) metal is established. It is applied to simulating the coupled evolution mechanisms of grains and twins in magnesium (Mg) by probing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries. We analyze the effect of the strain relaxation near a grain boundary, which is related to the geometric compatibility arising from the misorientation angle between adjoining grains. Phase-field simulations demonstrate that twin transmission across a grain boundary into a neighboring grain leads to grain boundary migration towards the neighboring grain with a reduced grain boundary width. The preferred nucleation site of new twin variant in the neighboring grain is related to the elastic interaction energy distribution. These predicted transmission behaviors agree well with existing experimental observations and molecular dynamics simulations. By analyzing set of systematic phase-field simulation results, we establish a rotation-angle-related twin variant selection rule for analyzing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries in Mg and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

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

Directed Gas-Phase Formation of Azulene (C 10 H 8 ): Unraveling the Bottom-Up Chemistry of Saddle-Shaped Aromatics

The azulene (C 10 H 8 ) molecule, the simplest polycyclic aromatic hydrocarbon (PAH) carrying a fused seven- and five-membered ring, is regarded as a fundamental molecular building block of saddle-shaped carbonaceous nanostructures such as curved nanographenes in the interstellar medium. However, an understanding of the underlying gas-phase formation mechanisms of this nonbenzenoid 10π-Hückel aromatic molecule under low-temperature conditions is in its infancy. Here, by merging crossed molecular beam experiments with electronic structure calculations and molecular dynamics simulations, our investigations unravel an unconventional low-temperature, barrierless route to azulene via the reaction of the simplest organic radical, methylidyne (CH), with indene (C 9 H 8 ) through ring expansion. This reaction might represent the initial step toward to the formation of saddle-shaped PAHs with seven-membered ring moieties in hydrocarbon-rich cold molecular clouds such as the Taurus Molecular Cloud-1 (TMC-1). These findings challenge conventional wisdom that molecular mass growth processes to nonplanar PAHs, especially those containing seven-membered rings, operate only at elevated pressure and high-temperature conditions, thus affording a versatile low-temperature route to contorted aromatics in our galaxy.

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

Molecular-scale origins of solution nanostructure and excess thermodynamic properties in a water/amphiphile mixture

The molecular and nanoscale origins of nonideality in excess thermodynamic properties are essential to understanding cosolvent mixtures, yet they remain challenging to determine. Here, we consider a binary mixture of water and an amphiphile, N,N,N',N'-tetramethylmalonamide (TMMA), which is characterized by strong hydrogen bonding between the two components and no hydrogen bonding between amphiphiles. Using molecular dynamics simulation, validated with excess volume measurements and X-ray scattering, we identify three distinct solution regimes across the composition range of the binary mixture and find that the transition between two of these regimes, marked by the water percolation threshold, is closely correlated with minima in the excess volume and excess enthalpy. Structural analysis of the simulations reveals an interplay between local interactions and solution nanostructure, determined by the relative strength of the water–water and water–amphiphile hydrogen bonding interactions. By comparison with other amphiphiles, such as linear alcohols, the relative strength of like and unlike interactions between water and amphiphile affects the relationship between thermodynamics and structural regimes. This provides insight into how molecular forces of mutual solvation interact across length scales and how they manifest in excess thermodynamic properties.

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

Selective filling of n-hexane in a tight nanopore

Abstract Molecular sieving may occur when two molecules compete for a nanopore. In nearly all known examples, the nanopore is larger than the molecule that selectively enters the pore. Here, we experimentally demonstrate the ability of single-wall carbon nanotubes with a van der Waals pore size of 0.42 nm to separate n-hexane from cyclohexane—despite the fact that both molecules have kinetic diameters larger than the rigid nanopore. This unexpected finding challenges our current understanding of nanopore selectivity and how molecules may enter a tight channel. Ab initio molecular dynamics simulations reveal that n-hexane molecules stretch by nearly 11.2% inside the nanotube pore. Although at a relatively low probability (28.5% overall), the stretched state of n-hexane does exist in the bulk solution, allowing the molecule to enter the tight pore even at room temperature. These insights open up opportunities to engineer nanopore selectivity based on the molecular degrees of freedom.

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

The Role of the Solvent in the Condensed-Phase Dynamics and Identity of Chemical Bonds: The Case of the Sodium Dimer Cation in THF

When a solute molecule is placed in solution, is it acceptable to presume that its electronic structure is essentially the same as in the gas phase? In this paper, we address this question from a simulation perspective for the case of the sodium dimer cation (Na 2 + ) molecule in both liquid Ar and liquid tetrahydrofuran (THF). In previous work, we showed that when local specific interactions between a solute and solvent are energetically on the order of a hydrogen bond, the solvent can become part of the chemical identity of the solute. Here, using mixed quantum/classical molecular dynamics simulations, we see that for the Na 2 + molecule, solute--solvent interactions lead to two stable, chemically-distinct coordination states (Na(THF) 4 -Na(THF) 5 + and Na(THF) 5 -Na(THF) 5 + ) that are not only stable themselves as gas-phase molecules but that also have a completely new electronic structure with important implications for the excited-state photodissociation of this molecule in the condensed phase. Furthermore, we show through a set of comparative classical simulations that treating the solute's bonding electron explicitly quantum mechanically is necessary to understand both the ground-state dynamics and chemical identity of this simple diatomic molecule; even use of the quantum-derived potential of mean force is insufficient to describe the behavior of the molecule classically. Lastly, we calculate the results of a proposed transient hole-burning experiment that could be used to spectroscopically disentangle the presence of the different coordination states.

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

Quantitative Simulations of Siloxane Adsorption in Metal–Organic Frameworks

We present a transferable force field (FF) for simulating the bulk properties of linear and cyclic siloxanes and the adsorption of these species in metal–organic frameworks (MOFs). Unlike previous FFs for siloxanes, our FF accurately reproduces the vapor–liquid equilibria of each species in the bulk phase. The quality of our FF combined with the Universal Force Field using standard Lorentz–Berthelot combining rules for MOF atoms was assessed in a wide range of MOFs without open metal sites, showing good agreement with dispersion-corrected density functional theory calculations. Predictions with this FF show good agreement with the limited experimental data for siloxane adsorption in MOFs that is available. As an example of using the FF to predict adsorption properties in MOFs, we present simulations examining entropy effects in binary linear and cyclic siloxane mixture coadsorption in the large-pore MOF with structure code FOTNIN.

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

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

Herein, we present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.

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

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. Furthermore, this approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Autonomous reinforcement learning agent for stretchable kirigami design of 2D materials

Abstract Mechanical behavior of 2D materials such as MoS 2 can be tuned by the ancient art of kirigami. Experiments and atomistic simulations show that 2D materials can be stretched more than 50% by strategic insertion of cuts. However, designing kirigami structures with desired mechanical properties is highly sensitive to the pattern and location of kirigami cuts. We use reinforcement learning (RL) to generate a wide range of highly stretchable MoS 2 kirigami structures. The RL agent is trained by a small fraction (1.45%) of molecular dynamics simulation data, randomly sampled from a search space of over 4 million candidates for MoS 2 kirigami structures with 6 cuts. After training, the RL agent not only proposes 6-cut kirigami structures that have stretchability above 45%, but also gains mechanistic insight to propose highly stretchable (above 40%) kirigami structures consisting of 8 and 10 cuts from a search space of billion candidates as zero-shot predictions.

36 MATERIALS SCIENCE↗