Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOLVATION”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Challenges in Modeling Electrochemical Reaction Energetics with Polarizable Continuum Models

A major challenge in the modeling of electrochemical phenomena is the accurate description of the interface between an electrolyte and a charged conductor. Polarizable continuum models (PCM) have been gaining popularity because they offer a computationally inexpensive method of modeling the electrolyte. In this Perspective, we discuss challenges from using one such model which treats the ions using a linearized Poisson–Boltzmann (LPB) distribution. From a physical perspective, this model places charge unphysically close to the surface and adsorbates, and it includes excessively steep ramping of the dielectric constant from the surface to the bulk solvent. Both of these issues can be somewhat mitigated by adjusting parameters built into the model, but in doing so, the resultant capacitance deviates from experimental values. Likewise, hybrid explicit-implicit approaches to the solvent may offer a more realistic description of hydrogen bonding and solvation to reaction intermediates, but the corresponding capacitances also deviate from experimental values. These deviations highlight the need for a careful adjustment of parameters in order to reproduce not only solvation energies but also other physical properties of solid–liquid interfaces. Continuum approaches alone also necessarily do not capture local variations in the electric field from cations at the interface, which can affect the energetics of intermediates with substantial dipoles or polarizability. Finally, since the double-layer charge can be varied continuously, LPB/PCM models provide a way to determine electrochemical barriers at constant potential. However, double-layer charging and the atomic motion associated with reaction events occur on significantly different timescales. As a result, we suggest that more detailed approaches, such as the modified Poisson–Boltzmann model and/or the addition of a Stern layer, may be able to mitigate some but not all of the challenges discussed.

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

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

Trends in Aqueous Hydration Across the 4f Period Assessed by Reliable Computational Methods

The geometric and electronic structures, as well as the thermodynamic properties of trivalent lanthanide hydrates {Ln(H₂O)₈,₉ 3+ and Ln(H₂O)₈,₉(H₂O)₁₂,₁₄ 3+ , Ln = La – Lu} have been examined using unrestricted density functional theory (UDFT), unrestricted M€oller-Plesset perturbation theory (UMP2), and multiconfigurational self-consistent field methods (MCSCF). While Ln-hydrates with 2-5 unpaired f-electrons have some multiconfigurational character, the correlation energy lies within 5-7 kcal/mol across the period and for varying coordination numbers. As such DFT yields structural parameters and thermodynamic data quite close to experimental values. Both UDFT and UMP2 predict free energies of water addition to the Ln(H₂O)₈ 3+ species to become less favorable across the period; however, it is a non-linear function of the surface charge density of the ion. UDFT further predicts that the symmetry of the metal-water bond lengths is sensitive to the specific f-electron configuration, presumably because of repulsive interactions between filled f-orbitals and water lone-pairs. Within the Ln(H₂O)₈,₉(H₂O)₁₂,₁₄ 3+ clusters, interactions between solvation shells overrides this orbital effect, increasing the accuracy of the geometric parameters and calculated vibrational frequencies. Calculated atomic charges indicate that the water ligands each donate 0.1 to 0.2 electrons to the Ln(III) metals, with increasing electron donation across the period. Significant polarization and charge transfer between solvation shells is also observed. The relationship between empirical effective charges and calculated atomic charges is discussed with suggestions for reconciling the trends across the period.

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

Thermochemistry of Lewis Adducts of BH 3 and Nucleophilic Substitution of Triethylamine on NH 3 BH 3 in Tetrahydrofuran

The thermochemistry of the formation of Lewis base adducts of BH 3 in tetrahydrofuran (THF) solution and the gas phase and the kinetics of substitution on ammonia borane by triethylamine are reported. The dative bond energy of Lewis adducts were predicted using density functional theory at the B3LYP/DZVP2 and B3LYP/6-311+G** levels and correlated ab initio molecular orbital theories, including MP2, G3(MP2), and G3(MP2)B3LYP, and compared with available experimental data and accurate CCSD(T)/CBS theory results. The analysis showed that the G3 methods using either the MP2 or the B3LYP geometries reproduce the benchmark results usually to within ~1 kcal/mol. Energies calculated at the MP2/aug-cc-pVTZ level for geometries optimized at the B3LYP/DZVP2 or B3LYP/6-311+G** levels give dative bond energies 2-4 kcal/mol larger than benchmark values. The enthalpies for forming adducts in THF were determined by calorimetry and compared with the calculated energies for the gas phase reaction: THFBH 3 + L → LBH 3 + THF. The formation of NH 3 BH 3 in THF was observed to yield significantly more heat than gas phase dative bond energies predict, consistent with strong solvation of NH 3 BH 3 . Substitution of NEt 3 on NH 3 BH 3 is an equilibrium process in THF solution (K ≈ 0.2 at 25 °C). The reaction obeys a reversible bimolecular kinetic rate law with the Arrhenius parameters: log A = 14.7 ± 1.1 and E a = 28.1 ± 1.5 kcal/mol. Finally, simulation of the mechanism using the SM8 continuum solvation model shows the reaction most likely proceeds primarily by a classical S N 2 mechanism.

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

Co-Solvent Electrolyte Engineering for Stable Anode-Free Zinc Metal Batteries

We report anode-free metal batteries can in principle offer higher energy density, but this requires them to have extraordinary Coulombic efficiency (>99.7%). Although Zn-based metal batteries are promising for stationary storage, the parasitic side reactions make anode-free batteries difficult to achieve in practice. In this work, a salting-in-effect-induced hybrid electrolyte is proposed as an effective strategy that enables both a highly reversible Zn anode and good stability and compatibility toward various cathodes. The as-prepared electrolyte can also work well under a wide temperature range (i.e., from -20 to 50 °C). It is demonstrated that in the presence of propylene carbonate, triflate anions are involved in the Zn 2+ solvation sheath structure, even at a low salt concentration (2.14 M). The unique solvation structure results in the reduction of anions, thus forming a hydrophobic solid electrolyte interphase. The waterproof interphase along with the decreased water activity in the hybrid electrolyte effectively prevents side reactions, thus ensuring a stable Zn anode with an unprecedented Coulombic efficiency (99.93% over 500 cycles at 1 mA cm –2 ). More importantly, we design an anode-free Zn metal battery that exhibits excellent cycling stability (80% capacity retention after 275 cycles at 0.5 mA cm –2 ). This work provides a universal strategy to design co-solvent electrolytes for anode-free Zn metal batteries.

25 ENERGY STORAGE↗
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 OSTI2011

Cluster-Models for Uranyl(VI) Adsorption on alpha-Alumina

Aqueous complexation, adsorption and redox chemistry of actinide species at mineral surfaces have significant impact on their transport and reactive behaviour in chemically and physically heterogeneous environments. The adsorption configurations and energies of micro-solvated uranyl dication, UO2 2+, on fully hydroxylated and H-deficient a-alumina (0001)-like finite cluster models were studied. The significant size of the models provides faithful representations of features that have emerged from periodic calculations, but most importantly, they afford us a systematic study of the preferred adsorption configurations, effect of secondary solvation shells and explicit treatment of the total charge. The energetics computed from the difference between the optimized structures and the appropriate reference states, point at a preference for inner sphere type complex.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Anion Effects on Interfacial Absorption of Gases in Ionic Liquids. A Molecular Dynamics Study

Molecular dynamics simulations with many-body interactions were carried out to systematic study the effect of anion type, tetrafluoroborate [BF4] or hexafluorophosphate [PF6], paired with the cation 1-butyl-3-methylimidazolium [bmim], on the interfacial absorption of gases in room temperature ionic liquids (RTILs). The potentials of mean force (PMF) of CO2 and H2O at 350 K were calculated across the air-liquid interfaces of [bmim][BF4] and [bmim][PF6]. We found that the PMFs for H2O exhibited no interfacial minima at both interfaces, while the corresponding PMFs for CO2 had significant free energy minima there. However, the PMFs for H2O showed a much higher interfacial free energy than in the bulk for [bmim][BF4], but only a slightly higher interfacial free energy for [bmim][PF6] than in bulk. The reason for this was due to the more hydrophilic nature of the [BF4] anion, and the fact that [BF4] was found to have little propensity for the interface. Our results show that H2O is much more likely to be found at the air-[bmim][PF6] interface than at the air-[bmim][BF4] interface. The free energies of solvation were found to be more negative for [bmim][BF4] than [bmim][PF6] for water and similar for CO2. This observation is consistent with experimental Henry’s law coefficients. Our results show that anion type, in addition to affecting the free energy of solvation into RTILs, should also significantly influence the uptake mechanism. This work was supported by the US Department of Energy Basic Energy Sciences' Chemical Sciences, Geosciences & Biosciences Division. Pacific Northwest National Laboratory is operated by Battelle for the US Department of Energy.

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

Hybrid Quantum Mechanical and Molecular Mechanics Study of the SN2 Reaction of CCl4 + OH- in Aqueous Solution: The Potential of Mean Force, Reaction Energetics, and Rate Constants

The bimolecular nucleophilic substitution reaction of CCl{sub 4} and OH{sup -} in aqueous solution was investigated on the basis of a combined quantum mechanical and molecular mechanics method. A multilayered representation approach is employed to achieve high accuracy results at the CCSD(T) level of theory. The potential of mean force calculations at the DFT level and CCSD(T) level of theory yield reaction barrier heights of 22.7 and 27.9 kcal/mol, respectively. Both the solvation effects and the solvent-induced polarization effect have significant contributions to the reaction energetics, for example, the solvation effect raises the saddle point by 10.6 kcal/mol. The calculated rate constant coefficient is 8.6 x 10{sup -28} cm{sup 3} molecule{sup -1} s{sup -1} at the standard state condition, which is about 17 orders magnitude smaller than that in the gas phase. Among the four chloromethanes (CH{sub 3}Cl, CH{sub 2}Cl{sub 2}, CHCl{sub 3}, and CCl{sub 4}), CCl{sub 4} has the lowest free energy activation barrier for the reaction with OH{sup -1} in aqueous solution, confirming the trend that substitution of Cl by H in chloromethanes diminishes the reactivity.

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

Ethylene glycol reforming on Pt(111): first-principles microkinetic modeling in vapor and aqueous phases

First-principles, periodic density functional theory (DFT) calculations and mean-field microkinetic modeling have been used to investigate the decomposition of ethylene glycol for hydrogen production on Pt(111) in vapor and aqueous phases. All dehydrogenated species derived from ethylene glycol (C 2 H x O 2 , x = 0–6) and methanol (CH y O, y = 0–4), and all elementary C–C, C–H, and O–H bond breaking steps are included in the microkinetic model. Reaction path analysis in vapor phase indicates that both initial C–H and O–H dehydrogenation steps are kinetically relevant at all temperatures (470–530 K). Initial O–H bond cleavage is reversible at low temperatures but accounts for an increasingly dominant fraction of the total reaction flux at higher temperatures. C–C bond scission is observed only after significant dehydrogenation (x ≤ 3) and is more likely to happen in surface intermediates where one of the cleavage products is CO. Here, the process is highly selective to the production of H2 compared to methanol. For aqueous-phase model development, free energies of solvation were computed for all surface intermediates and transition states using a continuum solvation approach. Our aqueous-phase microkinetic model predicts a 0.4 eV lower apparent activation energy and an order of magnitude larger turnover frequencies. Initial C–H bond cleavage becomes more important but the general trends are similar to the vapor phase, suggesting that the reaction chemistry is similar in both vapor and aqueous phases.

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

Understanding extreme fast charge limitations in carbonate mixtures

Charging large format lithium ion batteries within ten to fifteen minutes requires changes to the electrolyte composition in addition to modification of electrode and cell architectures. Several approaches to address this need have been proposed; but there is not a lot of clarity on understanding which factors limit the performance of existing electrolytes. Furthermore, this work takes a closer look at the solvated components of a mixture of LiPF 6 in ethylene carbonate (EC) and ethyl methyl carbonate (EMC) in the context of extreme fast charging and relates these findings to cell-level requirements. Molecular dynamics studies of the Gen-2 electrolyte compositions with increasing salt concentrations, have been performed to estimate transport properties like diffusivity, transference number and conductivity. Molecular-level differences in the structure of solvation shells under extreme LiPF 6 concentrations are probed here and some key aspects on solvent structure that help overcome barriers to Li + transport under extreme fast charge are discussed.

25 ENERGY STORAGE↗
DOE OSTI2011

Treecode-based generalized Born method

We have developed a treecode-based O(N log N) algorithm for the generalized Born (GB) implicit solvation model. Our treecode-based GB (tGB) is based on the GBr6 [J. Phys. Chem. B 111, 3055 (2007)], an analytical GB method with a pairwise descreening approximation for the R6 volume integral expression. The algorithm is composed of a cutoff scheme for the effective Born radii calculation, and a treecode implementation of the GB charge–charge pair interactions. Test results demonstrate that the tGB algorithm can reproduce the vdW surface based Poisson solvation energy with an average relative error less than 0.6% while providing an almost linear-scaling calculation for a representative set of 25 proteins with different sizes (from 2815 atoms to 65456 atoms). For a typical system of 10k atoms, the tGB calculation is three times faster than the direct summation as implemented in the original GBr6 model. Thus, our tGB method provides an efficient way for performing implicit solvent GB simulations of larger biomolecular systems at longer time scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

The oxidation of tyrosine and tryptophan studied by a molecular dynamics normal hydrogen electrode

The thermochemical constants for the oxidation of tyrosine and tryptophan through proton coupled electron transfer in aqueous solution have been computed applying a recently developed density functional theory (DFT) based molecular dynamics method for reversible elimination of protons and electrons. This method enables us to estimate the solvation free energy of a proton (H{sup +}) in a periodic model system from the free energy for the deprotonation of an aqueous hydronium ion (H{sub 3}O{sup +}). Using the computed solvation free energy of H{sup +} as reference, the deprotonation and oxidation free energies of an aqueous species can be converted to pK{sub a} and normal hydrogen electrode (NHE) potentials. This conversion requires certain thermochemical corrections which were first presented in a similar study of the oxidation of hydrobenzoquinone [J. Cheng, M. Sulpizi, and M. Sprik, J. Chem. Phys. 131, 154504 (2009)]. Taking a different view of the thermodynamic status of the hydronium ion, these thermochemical corrections are revised in the present work. The key difference with the previous scheme is that the hydronium is now treated as an intermediate in the transfer of the proton from solution to the gas-phase. The accuracy of the method is assessed by a detailed comparison of the computed pK{sub a}, NHE potentials and dehydrogenation free energies to experiment. As a further application of the technique, we have analyzed the role of the solvent in the oxidation of tyrosine by the tryptophan radical. The free energy change computed for this hydrogen atom transfer reaction is very similar to the gas-phase value, in agreement with experiment. The molecular dynamics results however, show that the minimal solvent effect on the reaction free energy is accompanied by a significant reorganization of the solvent.

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

Communication: Thermodynamic analysis of critical conditions of polymer adsorption

Polymer adsorption to solid surfaces is a ubiquitous phenomenon, which has attracted long-lasting attention. Dependent on the competition between the polymer-solid adsorption and polymer-solvent solvation interactions, a chain may assume either 3d solvated conformation when adsorption is weak or 2d adsorbed conformation when adsorption is strong. The transition between these conformations occurring upon variation of adsorption strength is quite sharp, and in the limit of “infinite” chain length, can be treated as a critical phenomenon. We suggest a novel thermodynamic definition of the critical conditions of polymer adsorption from the equality of incremental chemical potentials of adsorbed and free chains. We show with the example of freely jointed Lennard-Jones chains tethered to an adsorbing surface that this new definition provides a link between thermodynamic and geometrical features of adsorbed chains and is in line with classical scaling relationships for the fraction of adsorbed monomers, chain radii of gyration, and free energy.

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

Thermodynamic properties and diffusion of water + methane binary mixtures

Thermodynamic and diffusion properties of water + methane mixtures in a single liquid phase are studied using NVT molecular dynamics. An extensive comparison is reported for the thermal pressure coefficient, compressibilities, expansion coefficients, heat capacities, Joule-Thomson coefficient, zero frequency speed of sound, and diffusion coefficient at methane concentrations up to 15% in the temperature range of 298–650 K. The simulations reveal a complex concentration dependence of the thermodynamic properties of water + methane mixtures. The compressibilities, heat capacities, and diffusion coefficients decrease with increasing methane concentration, whereas values of the thermal expansion coefficients and speed of sound increase. Increasing methane concentration considerably retards the self-diffusion of both water and methane in the mixture. These effects are caused by changes in hydrogen bond network, solvation shell structure, and dynamics of water molecules induced by the solvation of methane at constant volume conditions.

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

Communication: Nucleation of quantized vortex rings in {sup 4}He nanodroplets

Whereas most of the phenomena associated with superfluidity have been observed in finite-size helium systems, the nucleation of quantized vortices has proven elusive. Here we show using time-dependent density functional simulations that the solvation of a Ba{sup +} ion created by photoionization of neutral Ba at the surface of a {sup 4}He nanodroplet leads to the nucleation of a quantized ring vortex. The vortex is nucleated on a 10 ps timescale at the equator of a solid-like solvation structure that forms around the Ba{sup +} ion. The process is expected to be quite general and very efficient under standard experimental conditions.

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

Development of three-dimensional site-site Smoluchowski-Vlasov equation and application to electrolyte solutions

Site-site Smoluchowski-Vlasov (SSSV) equation enables us to directly calculate van Hove time correlation function, which describes diffusion process in molecular liquids. Recently, the theory had been extended to treat solute-solvent system by Iida and Sato [J. Chem. Phys. 137, 034506 (2012)]. Because the original framework of SSSV equation is based on conventional pair correlation function, time evolution of system is expressed in terms of one-dimensional solvation structure. Here, we propose a new SSSV equation to calculate time evolution of solvation structure in three-dimensional space. The proposed theory was applied to analyze diffusion processes in 1M NaCl aqueous solution and in lithium ion battery electrolyte solution. The results demonstrate that these processes are properly described with the theory, and the computed van Hove functions are in good agreement with those in previous works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗