Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation structure”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ion Dissociation Dynamics in an Aqueous Premelting Layer

Using molecular dynamics simulations and methods of importance sampling, we study the thermodynamics and dynamics of sodium chloride in the aqueous premelting layer formed spontaneously at the interface between ice and its vapor. We uncover a hierarchy of time scales that characterize the relaxation dynamics of this system, spanning the picoseconds of ionic motion to the tens or hundreds of nanoseconds associated with fluctuations of the liquid-crystal interface in their presence. We find that ions distort both local interfaces, incurring restoring forces that result in the ions preferentially residing in the middle of the layer. While ion pair dissociation is thermodynamically favorable, these structural and dynamic effects cause its rate to vary by over an order of magnitude through the layer, with a maximum rate significantly depressed from the corresponding bulk value. The solvation environment of ions in the premelting layer is distinct from that in a bulk liquid, being dominated by slow reorganization of water molecules and a water structure intermediate between ice and its melt.

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

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Molecular-Level Insights into the Influence of Ionic Liquids on the Structure and Dynamics of Neutral and Charged Polyimides

We use all-atom molecular dynamics simulations to investigate how increasing ionic functionalization influences polyimide (PI) behavior in ionic liquid (IL) solvation environments. Here, we examine three polymer systems with varying charge densities: a neutral polyimide [N−PI] containing imidazole rings and progressively introduce charge through quaternization to create singly [C−PI] + and doubly charged [C−2PI] 2+ variants across a wide range of IL concentration (0−90 wt %). Through comprehensive structural, mechanical, and electrostatic analyses, we reveal that polymer charge density plays a central role in shaping IL organization and interaction with the polymer matrix. At low IL content, charged systems exhibit strong electrostatic complexation, leading to chain compaction, localized ESP environments, and elevated dielectric constants. As IL concentration increases, the effect of the different polymer charge states becomes less significant. Notably, an intermediate composition regime at approximately 50 wt % IL is associated with changes in IL-rich domain connectivity and overall system behavior across all three PI systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Compact Peptoid Molecular Brushes for Nanoparticle Stabilization

Controlling the interfaces and interactions of colloidal nanoparticles (NPs) via tethered molecular moieties is crucial for NP applications in engineered nanomaterials, optics, catalysis, and nanomedicine. Despite a broad range of molecular types explored, there is a need for a flexible approach to rationally vary the chemistry and structure of these interfacial molecules for controlling NP stability in diverse environments, while maintaining a small size of the NP molecular shell. Here, we demonstrate that low-molecular-weight, bifunctional comb-shaped, and sequence-defined peptoids can effectively stabilize gold NPs (AuNPs). The generality of this robust functionalization strategy was also demonstrated by coating of silver, platinum, and iron oxide NPs with designed peptoids. Each peptoid (PE) is designed with varied arrangements of a multivalent AuNP-binding domain and a solvation domain consisting of oligo-ethylene glycol (EG) branches. Among designs, a peptoid (PE5) with a diblock structure is demonstrated to provide a superior nanocolloidal stability in diverse aqueous solutions while forming a compact shell (similar to 1.5 nm) on the AuNP surface. We demonstrate by experiments and molecular dynamics simulations that PE5-coated AuNPs (PE5/AuNPs) are stable in select organic solvents owing to the strong PE5 (amine)-Au binding and solubility of the oligo-EG motifs. At the vapor-aqueous interface, we show that PE5/AuNPs remain stable and can self-assemble into ordered 2D lattices. The NP films exhibit strong near-field plasmonic coupling when transferred to solid substrates.

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

Luminescent Bimetallic Two-Coordinate Gold(I) Complexes Utilizing Janus Carbenes

A series of bimetallic carbene-metal-amide (cMa) complexes have been prepared with bridging biscarbene ligands to serve as a model for the design of luminescent materials with large oscillator strengths and small energy differences between the singlet and triplet states (ΔE ST ). The complexes have a general structure (R 2 N)Au(:carbene—carbene:)Au(NR 2 ). The bimetallic complexes show solvation-dependent absorption and emission that is analyzed in detail. It is found that the molar absorptivity of the bimetallic complexes is correlated with the energy barrier to rotation of the metal–ligand bond. The bimetallic cMa complexes also exhibit short emission lifetimes (τ = 200–300 ns) with high photoluminescence efficiencies (Φ PL > 95%). The radiative rates of bimetallic cMa complexes are 3–4 times faster than that of the corresponding monometallic complexes. Analysis of temperature-dependent luminescence data indicates that the lifetime for the singlet state (τS 1 ) of bimetallic cMa complexes is near 12 ns with a ΔE ST of 40–50 meV. Here the presented compounds provide a general design for cMa complexes to achieve small values for ΔE ST while retaining high radiative rates. Solution-processed organic light-emitting devices (OLEDs) made using two of the complexes as luminescent dopants show high efficiency and low roll-off at high luminance.

14 SOLAR ENERGY↗
DOE OSTI2026

Discerning Influences from Enthalpy and Entropy at Aqueous Interfaces Involved in Biomass Conversions in Porous Catalysts

Project Summary: The goal in this proposal is to learn how solvent influences the enthalpies and entropies of catalytic species in zeolite pores. Specifically, the physical, chemical, and structural features of solvent, catalytic species, and zeolite pores that determine enthalpies and entropies of solvation in solution phase biomass conversions will be interrogated using multiscale simulations and machine learning. Multiscale simulations are based off of existing strategies and employ quantum mechanics and classical molecular dynamics, providing an excellent balance between chemical accuracy and computational expense. They are capable of calculating enthalpies and entropies of solvation separately and have been validated in prior work to achieve high accuracy compared to their parent methods.

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

Charge accumulation and solvation in $β$-NiOOH: Surface chemistry of an OER catalyst from ML-aided simulations

Electrochemical water splitting is a key technology for a sustainable energy transition, providing a route to store surplus electricity from renewable sources. A central bottleneck is the sluggish oxygen evolution reaction (OER), which drives the search for catalysts that are active, stable, and inexpensive enough for large-scale deployment. Within this context, pure and doped NiO x H y combine high activity with low cost, making them prime candidates for alkaline OER. Yet, despite extensive study, the atomistic structure of NiOOH under operando conditions and the associated reaction mechanisms remain debated. Here, we investigate the structural complexity of pure β-NiOOH, the scaffold for its doped derivatives. We systematically investigate the oxidation of the surface adsorbates via proton-coupled electron transfer steps across relevant facets and sites, identifying the most probable sequence of deprotonation events. Our results reveal asymmetric charge accumulation on Wulff-relevant surfaces and show how applied potential can promote morphological restructuring. Explicit solvation is included through machine-learning interatomic potential molecular dynamics of the NiOOH/water interface, which allows us to resolve the hydrophobic and hydrophilic character of different surfaces and the associated interfacial water structure. Together, these insights demonstrate how surface chemistry and solvation jointly govern the stability of NiOOH and the accumulation of surface charge, with possible implications for catalytic performance.

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

An Aluminum-Based Metal–Organic Cage for Cesium Capture

Metal–organic cages are a class of supramolecular structures that often require the careful selection of organic linkers and metal nodes. Of this class, few examples of metal–organic cages exist where the nodes are composed of main group metals. Herein, we have prepared an aluminum-based metal–organic cage, H 8 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Al-pdc-AA), using inexpensive and commercially available materials. The cage formation was achieved via solvothermal self-assembly of solvated aluminum and pyridine-dicarboxylic linkers in the presence of a capping agent, acetic acid. The obtained supramolecular structure was characterized by single-crystal X-ray diffraction (SCXRD), thermogravimetric analysis, and NMR spectroscopy. Based on crystal structure and computational analyses, the cage has a 3.7 Å diameter electron-rich cavity suitable for the binding of cations such as cesium (ionic radius of 1.69 Å). Here, the host–guest interactions were probed with 1 H and 133 Cs NMR spectroscopy in DMSO, where at low concentrations, Cs + binds to Al-pdc-AA in a 1:1 ratio. The binding site was identified from the crystal structure of CsH 7 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Cs + Al-pdc-AA), and a binding affinity of ~10 6 –10 7 M –1 was determined from NMR titration experiments. The Al-pdc-AA showed improved selectivity for cesium binding over alkali metal cations (Cs + > Rb + > K + >> Na + ~ Li + ). Collectively, the study reports a novel aluminum cage that can serve as a promising host for efficient and selective cesium removal.

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

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Infrared Spectroscopy of Li + Solvation in EmimBF 4 and in Propylene Carbonate: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF 4 in the ionic liquid 1-ethyl-3-methylimidazolium tetrafluoroborate (EmimBF 4 ) and in the organic solvent propylene carbonate (PC) are studied via infrared spectroscopy and ab initio molecular dynamics (AIMD) simulations. The measurements show that the major effects of LiBF 4 in both solutions, compared to their neat counterparts, are the appearance of a new broad band in the 300-450 cm -1 frequency region and a broadening of the IR structure in the 900-1200 cm -1 region with the development of a new peak at 980 cm -1 . Computational analysis indicates that hindered translational motions of Li + in its solvation cage are mainly responsible for the former, while the latter is due to Li + -induced structural changes and accompanying vibrational frequency shifts of constituent ions and molecules of the solutions. In addition, molecular motions in these and lower frequency regions are generally correlated, disclosing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, detailed analysis of these features via AIMD simulations of the spectrum and its components arising from auto- and cross-correlations of motions of constituent molecular species, combined with generalized normal modes of the solutions and normal modes of small Li + -containing clusters, is presented. Furthermore, other minor spectral changes caused by the lithium salt as well as the interaction-induced effect on IR spectra are also discussed.

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

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2024

Influence of Ti‐incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H 2 O 2

Abstract Vapor‐phase propylene (C 3 H 6 ) epoxidation kinetics with hydrogen peroxide (H 2 O 2 ) strongly reflects the physical properties of Ti‐incorporated zeolite catalysts and the presence of spectating molecules (“solvent”) near active sites even without a bulk liquid phase. Steady‐state turnover rates of C 3 H 6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH) x ) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH 3 CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H 2 O 2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring‐opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C 3 H 6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH 3 CN. Thus, higher PO selectivities achieved in hydrophilic Ti‐MFI at steady‐state reflect the preferential stabilization of transition states for C 3 H 6 epoxidation (a primary reaction) relative to PO ring‐opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH 3 CN within pores and the extended pore structure.

Kwon, Ohsung↗
DOE OSTI2024

Influence of Ti‐incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H 2 O 2

Vapor-phase propylene (C 3 H 6 ) epoxidation kinetics with hydrogen peroxide (H 2 O 2 ) strongly reflects the physical properties of Ti-incorporated zeolite catalysts and the presence of spectating molecules (“solvent”) near active sites even without a bulk liquid phase. Steady-state turnover rates of C 3 H 6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH) x ) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH 3 CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H 2 O 2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring-opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C 3 H 6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH 3 CN. Thus, higher PO selectivities achieved in hydrophilic Ti-MFI at steady-state reflect the preferential stabilization of transition states for C 3 H 6 epoxidation (a primary reaction) relative to PO ring-opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH 3 CN within pores and the extended pore structure.

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

Imaging a solvent-swollen polymer gel network by open liquid transmission electron microscopy

Visualizing the network of a solvent-swollen polymer gel remains problematic. To address this challenge, open transmission electron microscopy (TEM) was applied to thin gel films permeated by a nonvolatile ionic liquid. The targeted physical gels were prepared by cooling concentrated solutions of poly(ethylene glycol) in 1-ethyl-3-methyl imidazolium ethyl sulfate [EMIM][EtSO 4 ]. During the cooling, gelation occurred by a frustrated crystallization of the dissolved polymer, leading to a percolated, solvent-permeated semicrystalline network in which nanoscale polymer crystals acted as crosslinks. Crystalline features ranging from ~5 to ~200 nm were observed, with the visible network strands dominantly consisting of long curvilinear crystallites of ~15–20 nm diameter. Nascent spherulites irregularly decorated the network, creating a complex structural hierarchy that complicated analyses. Lacking diffraction contrast, TEM did not visualize the many disordered, fully solvated PEG chains present in the voids between crystals. Finally, recognizing that a network's three dimensionality is ambiguous when assessed through two-dimensional microscopy projections, a small gel region was studied by TEM tomography, revealing a nearly isotropic three-dimensional arrangement of the curvilinear crystallites, which displayed remarkably uniform cylindrical cross sections.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗
DOE OSTI2023

The local ordering of polar solvents around crystalline carbon nitride nanosheets in solution

The crystalline graphitic carbon nitride, poly-triazine imide (PTI) is highly unusual among layered materials since it is spontaneously soluble in aprotic, polar solvents including dimethylformamide (DMF). The PTI material consists of layers of carbon nitride intercalated with LiBr. When dissolved, the resulting solutions consist of dissolved, luminescent single to multilayer nanosheets of around 60–125 nm in diameter and Li+ and Br- ions originating from the intercalating salt. To understand this unique solubility, the structure of these solutions has been investigated by high-energy X-ray and neutron diffraction. Although the diffraction patterns are dominated by inter-solvent correlations there are clear differences between the X-ray diffraction data of the PTI solution and the solvent in the 4–6 Å -1 range, with real space differences persisting to at least 10 Å. Structural modelling using both neutron and X-ray datasets as a constraint reveal the formation of distinct, dense solvation shells surrounding the nanoparticles with a layer of Br - close to the PTI-solvent interface. This solvent ordering provides a configuration that is energetically favourable underpinning thermodynamically driven PTI dissolution.

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

Dependency of solvation effects on metal identity in surface reactions

Solvent interactions with adsorbed moieties involved in surface reactions are often believed to be similar for different metal surfaces. However, solvents alter the electronic structures of surface atoms, which in turn affects their interaction with adsorbed moieties. To reveal the importance of metal identity on aqueous solvent effects in heterogeneous catalysis, we studied solvent effects on the activation free energies of the O–H and C–H bond cleavages of ethylene glycol over the (111) facet of six transition metals (Ni, Pd, Pt, Cu, Ag, Au) using an explicit solvation approach based on a hybrid quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface. A significant metal dependence on aqueous solvation effects was observed that suggests solvation effects must be studied in detail for every reaction system. The main reason for this dependence could be traced back to a different amount of charge-transfer between the adsorbed moieties and metals in the reactant and transition states for the different metal surfaces.

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

Investigating Lignin Aggregation and Interactions with Solvents during γ-Valerolactone (GVL) Pretreatment: A Combined Small Angle Neutron Scattering and Molecular Simulations Study

The strong tendency of lignin to aggregate in solution, coupled with limited understanding of how its molecular structure governs this behavior, hinders its effective utilization in biorefineries. Here, in this study, we investigated the solution behavior of lignin extracted from poplar using γ-valerolactone/water (GVL/H2O, 9:1 wt/wt) through combined small-angle neutron scattering (SANS) and molecular dynamics (MD) simulations. Lignin samples obtained at 100 °C (L100) and 120 °C (L120) differed in β–O–4 content, hydroxyl distribution, and S/G ratio, enabling direct assessment of how molecular composition governs solvation and aggregation. SANS showed that L120 formed rigid and elongated cylindrical aggregates at 25 °C that transitioned to more flexible spheroidal structures by 50 °C and remained stable up to 80 °C, whereas L100 adopted globular aggregates that progressively collapsed with increasing temperature. MD simulations reinforced these observations by showing that S-rich (L120-like) oligomers had larger radii of gyration, stronger solvent coordination driven by methoxy groups, and fewer lignin–lignin contacts. In contrast, G-rich (L100-like) oligomers displayed persistent aggregation and lower solubility. Collectively, these results reveal that increased aromatic methoxylation enhances lignin–solvent interactions and suppresses self-association, whereas reduced methoxylation and higher β–O–4 content promote persistent aggregation into colloid-like structures with restricted solvent penetration into the aggregate interior.

biorefinery pretreatment↗