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 361 records · Page 20

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

The Solvation of the E. coli CheY Phosphorylation Site Mapped by XFMS

The Escherichia coli CheY protein belongs to a large bacterial response regulator superfamily. X-ray hydroxy radical foot-printing with mass spectroscopy (XFMS) has shown that allosteric activation of CheY by its motor target triggers a concerted internalization of aromatic sidechains. We reanalyzed the XFMS data to compare polar versus non-polar CheY residue positions. The polar residues around and including the 57D phosphorylated site had an elevated hydroxy radical reactivity. Bioinformatic measures revealed that a water-mediated hydrogen bond network connected this ring of residues with the central 57D. These residues solvated 57D to energetically stabilize the apo-CheY fold. The abundance of these reactive residues was reduced upon activation. This result was supported by the bioinformatics and consistent with the previously reported activation-induced increase in core hydrophobicity. It further illustrated XFMS detection of structural waters. Direct contacts between the ring residues and the phosphorylation site would stabilize the aspartyl phosphate. In addition, we report that the ring residue, 18R, is a constant central node in the 57D solvation network and that 18R non-polar substitutions determine CheY diversity as assessed by its evolutionary trace in bacteria with well-studied chemotaxis. These results showcase the importance of structured water dynamics for phosphorylation-mediated signal transduction.

59 BASIC BIOLOGICAL SCIENCES↗

ReaxFF molecular dynamics simulations of electrolyte–water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal–chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of Peripheral Substituents on Fe(II) Spin State in Complexes with Tridentate Schiff‐Base Ligands

Abstract Six Fe(II) complexes with substituted 6‐((quinoline‐8‐ylimino)methyl)phenolates (R 2 qsal) were synthesized. X‐ray crystal structure determination for [Fe(R 2 qsal) 2 ], where R=F ( 1 ), Cl ( 2 ), Br ( 3 ), I ( 4 ), CF 3 /H ( 5 ), or CN/H ( 6 ), revealed mononuclear complexes that crystallize either solvent‐free in case of 1 , 5 , and 6 or as CH 2 Cl 2 solvates in the case of 2 , 3 , and 4 . Complexes 1 – 4 , which contain π‐donating substituents, exhibit only high‐spin (HS) state, as shown by crystal structure analysis and magnetic measurements. In contrast, complexes 5 and 6 , which contain electron‐withdrawing substituents, exist in the low‐spin (LS) state at lower temperatures but exhibit crossover to the HS state above 300 K. Electronic structure calculations show a good correlation between the relative energies of the LS and HS states and the experimentally observed behavior.

Gakiya‐Teruya, Miguel↗

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗