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

Probing Hydrogen-Bond Anharmonicity in Cocrystals with Colossal Thermal Expansion

Organic cocrystals with sizable anisotropic thermal expansion (TE) are attractive building blocks for thermally responsive soft actuators. Such materials often display intermolecular H-bonding, whose thermal strength fluctuations are frequently overlooked if displacements and/or conformational changes are visible during TE. In this work, we employed variable-temperature single-crystal synchrotron diffraction to investigate the linear expansion coefficients of two organic cocrystals, 1 and 2, which exhibit colossal thermal expansivity exceeding 140 MK –1 . Despite such pronounced TE, solid-state 2 H and 13 C NMR experiments show no evidence of fast internal reorientations or large-amplitude molecular displacement during expansion. Although nondirectional π–π stacking dominates the colossal TE, our results indicate that hydrogen bonds play a non-negligible role in the anisotropic response. Molecular dynamics simulations reveal a temperature-dependent broadening of H-bond length distributions, which arises from the increased thermal sampling of an anharmonic potential energy surface. Our data support an expansion mechanism based on cooperativity in π-π interactions, accompanied by thermal softening and fluctuation-driven stretching of directional hydrogen bonds in these molecular solids. This mechanism underscores the significance of directional interactions in organic cocrystals, which is a key factor in the design of organic thermoresponsive crystalline materials.

Hernández-Morales, Ernesto A. [Univ. Nacional Auto↗
DOE OSTI2026

Atomic-Scale Imaging Reveals Polar-π Interactions in Two-Dimensional Molecular Superlattices

Controlling coassembly of synthetic oligomers into binary superlattices at the atomic level is challenging. Here, we report a strategy for programming polar-π interactions in oligomeric peptoids, a class of sequence-defined peptidomimetics, facilitating the formation of homogeneous two-dimensional (2D) superlattices. N-2-phenylethyl and N-(2-perfluorophenyl)ethyl side chains, similar in size, but with contrasting electrostatic characteristics, were introduced at defined sequence positions to generate favorable dipolar aromatic interactions. The resulting nanosheets exhibit different crystal motifs depending on the side chain interactions: systems containing only one type of aromatic side chain form a parallel V-shaped motif driven by π-π interactions, whereas a combination of both types of aromatic side chains, either within one backbone or through the coassembly of two distinct peptoids, adopt an antiparallel V-shaped superlattice with higher thermal stability, driven by polar-π interactions. Cryogenic transmission electron microscopy directly resolved the packing arrangement of perfluorophenyl and phenyl rings in individual nanosheet superlattices, confirming that intermolecular polar-π interaction dominates the superlattice motifs and increases lattice stability. Molecular dynamics simulations and density functional theory calculations further substantiate the energetic favorability of polar-π interactions over π-π interactions, rationalizing the formation of homogeneous superlattices with enhanced thermal stability. Our discoveries establish a design principle for binary coassembly using sequence-defined oligomers, which enables control over unit cell geometry, lattice stability, and molecular registration through aromatic side chain polarization and sequence control. This ability to program atomic-scale binary superlattices opens new avenues for designing functional 2D soft materials.

Lee, Yen Jea [Lawrence Berkeley National Laborator↗
DOE OSTI2026

Conformational Isomerization of Imide Anions Governs Solvation and Transport in Water-in-Salt Electrolytes

The behavior of highly concentrated electrolytes departs radically from the dilute-solution theory, yet the molecular origin of this transformation remains unresolved. Here, we identify the conformational isomerization of molecular ions as a decisive, previously unrecognized control parameter governing structure and transport in crowded aqueous electrolytes. Across a series of fluorosulfonimide anions, we show that increasing concentration drives a collective shift from extended transoid to compact cisoid conformers, revealed by small-angle X-ray scattering, vibrational spectroscopy, pulsed-field gradient NMR, and molecular dynamics simulations. This conformational transition triggers a collapse of the hydrogen-bonded water network and the emergence of densely packed ionic domains with confined water, producing a qualitative change in Li+ transport from solvent-mediated diffusion to network-confined hopping. Anion size and asymmetry systematically tune the onset of this transition, demonstrating that molecular geometry dictates mesoscale organization and dynamics in the ion-rich regime. Our results establish ion conformation, not merely composition or coordination, as a fundamental thermodynamic variable in concentrated solutions, providing a chemical framework that unifies solvation structure and transport in water-in-salt electrolytes and suggesting new principles for designing dense ionic media.

Nguyen, Huong TD↗
DOE OSTI2021

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

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

Water, Ion, and Proton Dynamics and Interactions in Nanoconfined Systems

The research was directed at understanding the role that nanoconfinement has on the dynamics and structure of liquids confined in nanometer size structures. Experiments were conducted using ultrafast infrared methods. The methods were two dimensional infrared (2D IR) spectroscopy and polarization pump-probe (PSPP) experiments. In addition FT-IR spectroscopy was used. To augment the experimental work, high level molecular dynamics simulations of the systems were conducted. The system used for confinement were nanoporous silica with essentially monodispersed pores. Pores of different diameters were studies, i.e., 2.5 nm to 8 nm. Water and organic liquids were investigated. The experiments were performed on a vibrational probe molecule, i.e., the CN stretch of the anion SeCN-. The 2D IR experiments measured spectral diffusion, which directly reported on the structural dynamics of the confined liquids. 2D IR was also used to perform chemical exchange measurements, which observed the formation and dissociation of molecular complexes with surface silanols. The PSPP experiments measured orientational relaxation of the probe molecule as well as it vibrational lifetime. All of the observables depended on the diameter of the pores and the type of liquid. The simulations of the observables provided detailed structural information on the structure of the confined liquids and there interactions with the surfaces. In addition to the study of the confinement in nanopores, water in hydrogels was also studied. The same experimental methods were used along with MD simulations. It was found the influence of hydrogels on water did not come from confinement of the water in the hydrogel polymer network, but rather from the interactions of the water with the individual monomers that compose the polymer. In addition, the temperature dependence of a polymer/water system was tracked from below to above the lower critical solution temperature to observe the nature of the phase separated polymer and how water responded to the phase separation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
DOE OSTI2020

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

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

Understanding Gas Absorption in Multivalent Ionic Liquids via Solute-Solvent Interaction Analyses

Solute-solvent interactions are often evaluated to understand and explain gas solubilities. This information is regularly used to interpret gas solubility and separation within ionic liquids (ILs). Here, we use several different quantum chemical and previous molecular dynamics simulation results to show that strong solute-solvent interactions neither lead to high solubility nor to high selectivity, with respect to other solutes. Instead, the solute-solvent interaction strength with different isolated IL ions is affected by the electrostatic extremes. Also, the solute-solvent interaction with different IL molecules is strongly correlated to the molecular volumes or polarizabilities (within a series of same-anion based ILs).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗