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

Quantum Chemical Calculation of p K a s of Environmentally Relevant Functional Groups: Carboxylic Acids, Amines, and Thiols in Aqueous Solution

Developing accurate quantum chemical approaches for calculating p K a s is of broad interest. Useful accuracy can be obtained by using density functional theory (DFT) in combination with a polarizable continuum solvent model. However, some classes of molecules present problems for this approach, yielding errors greater than 5 p K units. Various methods have been developed to improve the accuracy of the combined strategy. These methods perform well but either do not generalize or introduce additional degrees of freedom, increasing the computational cost. The Solvation Model based on Density (SMD) has emerged as one of the most commonly used continuum solvent models. Nevertheless, for some classes of organic compounds, e.g., thiols, the p K a s calculated with the original SMD model show errors of 6–10 p K units, and we traced these errors to inaccuracies in the solvation free energies of the anions. To improve the accuracy of p K a s calculated with DFT and the SMD model, we developed a scaled solvent-accessible surface approach for constructing the solute–solvent boundary. By using a “direct” approach, in which all quantities are computed in the presence of the continuum solvent, the use of thermodynamic cycles is avoided. Furthermore, no explicit water molecules are required. Three benchmark data sets of experimentally measured p K a values, including 28 carboxylic acids, 10 aliphatic amines, and 45 thiols, were used to assess the optimized SMD model, which we call SMD with a scaled solvent-accessible surface (SMD sSAS ). Of the methods tested, the M06-2X density functional approximation, 6-31+G(d,p) basis set, and SMD sSAS solvent model provided the most accurate p K a s for each set, yielding mean unsigned errors of 0.9, 0.4, and 0.5 p K units, respectively, for carboxylic acids, aliphatic amines, and thiols. As a result, this approach is therefore useful for efficiently calculating the p K a s of environmentally relevant functional groups.

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

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

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

Computer simulation of methanol exchange dynamics around cations and anions

In this paper, we present the first computer simulation of methanol exchange dynamics between the first and second solvation shells around different cations and anions. After water, methanol is the most frequently used solvent for ions. Methanol has different structural and dynamical properties than water, so its ion solvation process is different. To this end, we performed molecular dynamics simulations using polarizable potential models to describe methanol-methanol and ion-methanol interactions. In particular, we computed methanol exchange rates by employing the transition state theory, the Impey-Madden-McDonald method, the reactive flux approach, and the Grote-Hynes theory. We observed that methanol exchange occurs at a nanosecond time scale for Na+ and at a picosecond time scale for other ions. We also observed a trend in which, for like charges, the exchange rate is slower for smaller ions because they are more strongly bound to methanol. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. The calculations were carried out using computer resources provided by the Office of Basic Energy Sciences.

Roy, Santanu↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2016

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid–Base and Ligand Binding Equilibria of Aquacobalamin

In redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. Moreover, a major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co ligand binding equilibrium constants (Kon/off), pKas and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co III , Co II , and Co I species, respectively, and the second model features saturation of each vacant axial coordination site on Co II and Co I species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pK a s and 2.3 log units for two log K on/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co axial ligand binding, leading to substantial errors in predicted pK a s and K on/off values. Finally, these findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

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

The Influence of Water on the Alkyl Region Structure in Variable Chain Length Imidazolium-Based Ionic Liquid/Water Mixtures

In this work, solutions of room temperature ionic liquids (RTILs) and water were studied by observing the reorientational dynamics of the fluorescent probe perylene. Perylene is solvated in the alkyl regions of the RTILs. Its D 2h symmetry made it possible to extract dynamical information on both in-plane and out-of-plane reorientation from time resolved fluorescence anisotropy measurements. Perylene reorientation reports on its interactions with the alkyl chains. The RTILs were a series of 1-alkyl-3-methylimidazolium tetrafluoroborates (C n mimBF 4 , where n is the number of carbons in the alkyl chain), and the effects on perylene’s dynamics were observed when varying the alkyl chain length of the cation (n = 4, 6, 8, and 10; butyl, hexyl, octyl, decyl) and varying the water content from pure RTIL to roughly three water molecules per RTIL ion pair. Time correlated single photon counting was used to measure the fluorescence anisotropy decays to determine the orientational dynamics. The friction coefficients for both the in-plane and out-of-plane reorientation were determined to eliminate the influence of changes in viscosity caused by both the addition of water and the different alkyl chain lengths. The friction coefficients provided information on the interactions of the perylene with it alkyl environment and how these interactions changed with chain length and water content. As chain length increased, the addition of water had less of an effect on the local alkyl environment surrounding the perylene. The friction coefficients generally increased with higher water contents; the inplane orientational motion was hindered significantly more than the out-of-plane motion. The restructuring of the alkyl regions is likely a consequence of a rearrangement of the ionic imidazolium head groups to accommodate partial solvation by water, which results in a change in the arrangement of the alkyl chains. At very high water content, BmimBF4 broke this general trend, with both in-plane and out-of-plane rotational friction decreasing above a water content of one water per ion pair. This decrease indicates a major reorganization of the overall liquid structure in high water content mixtures. In contrast to BmimBF 4 , the longer chain length RTILs are not infinitely miscible with water, and do not show evidence of a major reorganization before reaching saturation and phase-separating. The results suggest that phase separation in longer chain length BF 4 RTILs is a consequence of their inability to undergo the reorganization of the alkyl regions necessary to accommodate high water concentrations.

14 SOLAR ENERGY↗
DOE OSTI2017

Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube

Recently, a collective variable “proton indicator” was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. Finally, the results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.

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

Comprehensive Experimental and Computational Spectroscopic Study of Hexacyanoferrate Complexes in Water: From Infrared to X-ray Wavelengths

We present a joint experimental and computational study of the hexacyanoferrate aqueous complexes at equilibrium in the 250 meV to 7.15 keV regime. The experiments include the vibrational spectroscopy of the cyanide ligands, the valence electronic absorption spectra and Fe 1s core hole spectra using X-ray absorption and emission techniques. These experiments are accompanied by density functional theory (DFT) based quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations and spectroscopy calculations of the experiments spanning the IR to X-ray wavelengths. Our simulations, including explicit solvent effects, across this large energy window allows for a systematic comparison of vibrational, valence electronic and core-level spectra. The analysis reveals how the local solvation environment influences the geometric and electronic structure and impacts the line shape and intensity of the observed spectral features in the infrared, optical and X-ray wavelengths.

Ross, Matthew↗
DOE OSTI2018

Ab Initio Molecular Dynamics Reveal Spectroscopic Siblings and Ion Pairing as New Challenges for Elucidating Prenucleation Aluminum Speciation

We report the characterization of prenucleation species is essential to understand crystallization mechanisms across many chemical systems and often involves the use of vibrational spectroscopy. Nowhere is this more evident than in the development of “green” aluminum processing technologies, where detailed understanding of the speciation of aluminum and its polynuclear analogues in highly alkaline, low water solutions is elusive. The aluminate anion Al(OH) 4 – predominates in alkaline conditions, yet equilibrium with dimeric species, either μ-oxo Al 2 O(OH) 6 2– or di-μ-hydroxo Al 2 (OH) 8 2– , can be assumed. Using ab initio molecular dynamics with full solvation and the presence of counterions, this work reconciles previous contradictory studies that had concluded only a single species under relevant solution conditions. We reveal that the two dimers are energetically separated by 2 kcal/mol in pure water but that the stability of each can be reversed by ion pairing expected in saturated salt solutions. Simulated Raman and IR spectra for each species (accounting for anharmonicity and the fluctuating solvating environment) provide the first proof that the considered species are “spectroscopic siblings”, whose multiple overlapping bands prevent definitive assertions in terms of speciation when compared to the experimental spectra. Lastly, these observations are likely to hold in higher order aluminate oligomers and as such present a massive challenge toward understanding the crystallization mechanisms relevant to aluminum processing.

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

Rate versus Free Energy Change for Attaching Highly Mobile Electrons to Molecules in Nonpolar Liquids

The inverted region of the Marcus theory, usually absent for bimolecular electron transfer reactions, is clearly observed for electron attachment reactions to molecules in nonpolar fluids. Application of pressure increased the energies of the solvated electron reactants letting us continuously adjust the free energy change. Inverted behavior is enabled by the very high mobilities of the solvated electrons which raise the diffusion-controlled encounter rates so high that they do not limit the reaction rates. The nonpolar media used in these experiments reduce reorganization energies, enhancing inverted behavior. Still, for every case showing an inverted region, the presence of low-lying excited states in the product radical anions led to regions of increasing rate constants that began at the energies of excited states of –0.54 to –1.2 eV. Here, while continuum models predict no solvent reorganization energy in nonpolar liquids, fits to the data found solvent reorganization energies of 0.05–0.4 eV supporting ideas advanced in theories of Matyushov.

14 SOLAR ENERGY↗
DOE OSTI2020

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

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

Surface-Mediated Solvent Decomposition in Li–Air Batteries: Impact of Peroxide and Superoxide Surface Terminations

A viable Li/O 2 battery will require the development of stable electrolytes that do not continuously decompose during cell operation. In some recent experiments it is suggested that reactions occurring at the interface between the liquid electrolyte and the solid lithium peroxide (Li 2 O 2 ) discharge phase are a major contributor to these instabilities. To clarify the mechanisms associated with these reactions, a variety of atomistic simulation techniques, classical Monte Carlo, van der Waals-augmented density functional theory, ab initio molecular dynamics, and various solvation models, are used to study the initial decomposition of the common electrolyte solvent, dimethoxyethane (DME), on surfaces of Li 2 O 2 . Comparisons are made between the two predominant Li 2 O 2 surface charge states by calculating decomposition pathways on peroxide-terminated (O 2 2– ) and superoxide-terminated (O 2 1– ) facets. For both terminations, DME decomposition proceeds exothermically via a two-step process comprised of hydrogen abstraction (H-abstraction) followed by nucleophilic attack. In the first step, abstracted H dissociates a surface O 2 dimer, and combines with a dissociated oxygen to form a hydroxide ion (OH – ). In the remaining surface oxygen then attacks the DME, resulting in a DME fragment that is strongly bound to the Li 2 O 2 surface. DME decomposition is predicted to be more exothermic on the peroxide facet; nevertheless, the rate of DME decomposition is faster on the superoxide termination. The impact of solvation (explicit vs implicit) and an applied electric field on the reaction energetics are investigated. Finally, our calculations suggest that surface-mediated electrolyte decomposition should out-pace liquid-phase processes such as solvent auto-oxidation by dissolved O 2 .

25 ENERGY STORAGE↗
DOE OSTI2016

Effects of High and Low Salt Concentration in Electrolytes at Lithium–Metal Anode Surfaces

The use of high concentration salts in electrolyte solutions of lithium-sulfur (Li-S) batteries has been shown beneficial for mitigating some effects such as polysulfide shuttle and dendrite growth at the Li metal anode. Such complex solutions have structural, dynamical, and reactivity associated issues that need to be analyzed for a better understanding of the reasons behind such beneficial effects. A passivation interfacial layer known as solid-electrolyte interphase (SEI) is generated during battery cycling as a result of electron transfer from the metal anode causing electrolyte decomposition. Here in this work, we investigate using density functional theory and ab initio molecular dynamics simulations the salt decomposition, solvation effects, interactions among intermediate products and other species, and potential components of the SEI layer as a function of chemical nature and concentration of the salt, for lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and lithium bis(fluorosulfonyl)imide (LiFSI) at 1M and 4M concentrations in dimethoxyethane. It is found that LiTFSI undergoes a less complete reduction and facilitates charge transfer from the anode, whereas LiFSI shows a more complete decomposition forming LiF as one of the main SEI products. In addition, the specific decomposition mechanisms of each salt clearly point to the initial SEI components and the potential main products derived from them. Finally, very complex networks are found among the salt and solvent molecules in their attempt to maximize Li ion solvation that is quantified through the determination of coordination numbers.

25 ENERGY STORAGE↗
DOE OSTI2017

First-Principles Investigation of Lithium Polysulfide Structure and Behavior in Solution

We present the Lithium-Sulfur battery is a promising next generation energy storage technology that could meet the demands of modern society with a theoretical specific energy near 2500 W h kg -1 . However, this battery chemistry faces unique problems such as the parasitic polysulfide shuttle reaction which hinders battery performance severely. This shuttle phenomenon is caused by solubilities of intermediate reaction products in the electrolyte during the reduction chemistry of the battery. With molecular simulation and computational chemistry tools, we studied the thermodynamics, solvation structure, and dynamics of the long-chain lithium polysulfide species Li 2 S 6 and Li 2 S 8 in dimethoxyethane and 1,3-dioxolane to gain a deeper fundamental understanding of this process. We determined the structure of the 1st solvation shell for Li + as well as those of Li 2 S 6 , Li 2 S 8 closed and Li 2 S 8 linear in pure solvents and solvents with extra Li + added. Finally, the lithium polysulfide species were found not to favor dissociation and would most likely exist as fully lithiated species in solution.

25 ENERGY STORAGE↗
DOE OSTI2017

Atomic-Scale Quantification of Interfacial Binding between Peptides and Inorganic Crystals: The Case of Calcium Carbonate Binding Peptide on Aragonite

Using a specific explicitly solvated interface model between a calcium carbonate binding peptide and crystalline aragonite, we investigate the electronic structure, atomic bonding, solvation effect, and the role of hydrogen bonding on the cohesion, stability, and functionality of this complex hybrid system using density functional calculation. The large interface model is strategically constructed using a stepwise procedure followed by ab initio molecular dynamics to obtain the optimal conformation. The calculated data on the electronic structure and bonding are analyzed in terms of three structural parts: aragonite, peptide, and water. Next, we focus on the binding between aragonite (001) surface and the peptide mediated by water. Finally, specific interatomic bonding between the amino acids in peptide and the (001) surface of aragonite is quantified. A single quantum mechanical metric, the total bond order density (TBOD), infers the dynamic interplay of different competing interactions. Four amino acids HIS1, ARG6, MET7, and TRP11 in the peptide sequence have strong interfacial Ca–O bonding and O···H hydrogen bonding between aragonite and peptide. The calculated Young’s modulus 33.37 GPa is in line with the measured value for nacre. Our approach for interfacial study between aragonite and a calcium carbonate binding peptide offers a broad perspective for probing complex interactions between the biomimetic interfaces. Finally, TBOD can be used as an effective parameter in ranking the efficacy of peptide–surface interactions and in providing a programmable design for bio-inspired material interfaces based on computational means.

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

Effect of Solvent and Substrate on the Surface Binding Mode of Carboxylate-Functionalized Aromatic Molecules

Here, the efficiency of dye-sensitized solar cells (DSSCs) is strongly influenced by dye molecule orientation and interactions with the substrate. Understanding the factors controlling the surface orientation of sensitizing organic molecules will aid in the improvement of both traditional DSSCs and other devices that integrate molecular linkers at interfaces. Here, we describe a general approach to understand relative dye–substrate orientation and provide analytical expressions predicting orientation. We consider the effects of substrate, solvent, and protonation state on dye molecule orientation. In the absence of solvent, our model predicts that most carboxylic acid-functionalized molecules prefer to lie flat (parallel) on the surface, due to van der Waals interactions, as opposed to a tilted orientation with respect to the surface that is favored by covalent bonding of the carboxylic acid group to the substrate. When solvation effects are considered, however, the molecules are predicted to orient perpendicular to the surface. We extend this approach to help understand and guide the orientation of metal–organic framework (MOF) thin-film growth on various metal–oxide substrates. A two-part analytical model is developed on the basis of the results of DFT calculations and ab initio MD simulations that predicts the binding energy of a molecule by chemical and dispersion forces on rutile and anatase TiO 2 surfaces, and quantifies the dye solvation energy for two solvents. The model is in good agreement with the DFT calculations and enables rapid prediction of dye molecule and MOF linker binding preference on the basis of the size of the adsorbing molecule, identity of the surface, and the solvent environment. We establish the threshold molecular size, governing dye molecule orientation, for each condition.

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

Theoretical Approaches to Describing the Oxygen Reduction Reaction Activity of Single-Atom Catalysts

Single-atom catalysts have recently emerged as promising low-cost alternatives to Pt for the oxygen reduction reaction (ORR). Given the unique properties that distinguish these systems from traditional transition-metal electrocatalysts, it is essential to benchmark and establish appropriate computational approaches to study these novel materials. Herein, we employ multiple levels of theory, including wave function methods, density functional theory (DFT), and classical simulations, to investigate Cu-modified covalent triazine framework catalysts (Cu/CTF). We consider three major aspects of treating this system computationally. First, we present a step-wise approach to predict the ORR mechanism and adsorbate coverages on Cu/CTF. We then benchmark various DFT methods to coupled-cluster theory with the domain-based local pair natural orbital approximation, which indicates that HSE06 and PBE0 hybrid functionals most accurately describe the adsorption energies of ORR adsorbates on Cu/CTF. We finally employ thermodynamic integration and other techniques to consider solvation effects, which play significant roles in predicting the energies of reaction intermediates and the overall ORR pathway. Lastly, our findings indicate that accurate descriptions of both the electronic structure and solvation are necessary to understand the ORR activity of Cu/CTF.

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

Solvent–Adsorbate Interactions and Adsorbate-Specific Solvent Structure in Carbon Dioxide Reduction on a Stepped Cu Surface

In this work, the structure of water and its interactions with various carbon dioxide reduction intermediates adsorbed on a Cu(211) surface is investigated using density functional theory. We find that the presence of adsorbates has a significant and adsorbate-specific effect on the local water structure and that solvation can stabilize adsorbate conformations different from those found in vacuum. We describe relationships between the hydrogen bonding capability of an adsorbate, the dipole moment of the adsorbate, the energetic strength of water–adsorbate interactions, and the change induced in the local water orientation by the adsorbate. Mechanistic implications are discussed. Here, we investigate and quantify the error associated with using arbitrary locally optimized solvent structures in calculations of relevant physical quantities, such as solvated binding energies and work functions. Possible effects of thermal motion on calculations of the work function are investigated using ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗