Theoretical Investigation of the Potential-Dependent CO Adsorption on Copper Electrodes.

Despite the importance of CO adsorption in many electrocatalytic reaction mechanisms, there has been little investigation of the dependence of the free energy of CO adsorption on the applied potential. Herein, we report on the potential-dependent adsorption of CO on Cu electrodes using a grand-canonical density functional theory approach. We demonstrate that, within the working potential range of electrocatalytic CO2 reduction on Cu(111) and Cu(100), the CO adsorption strength can change by over 0.1 eV. Our analyses explain the potential dependence through an interfacial capacitance loss upon CO adsorption as well as orbital relaxation induced by the electrode potential. Via sensitivity analysis with respect to two electrolyte model parameters (solvent dielectric constant and Debye screening length), we find that the surface excess charge density is a useful descriptor of the CO adsorption free energy.

Identifying the Local Atomic Environment of the "Cu3+" State in Alkaline Electrochemical Systems.

CuO-based catalysts are active for the oxygen evolution reaction (OER), although the active form of copper for the OER is still unknown. We combine operando Raman experiments and density functional theory (DFT) electronic structure calculations to determine the form of Cu(O)xOHy present under OER conditions. Raman spectra show a distinct feature related to the active "Cu3+" species, which is only present under highly oxidizing conditions. DFT is used to produce theoretical Raman standards and match the unique Raman feature of copper under OER potentials. This method identifies a range of Cu3+-containing compounds which match the distinct Raman feature. We then integrate experimental electrochemistry to progressively eliminate possible structures and determine the stoichiometry of the active form as CuOOH, which likely takes the form of a surface-adsorbed hydroxide on a CuO surface. Bader charge analysis, site-projected wavefunctions, and density of states analysis show that electron density is removed from O 2p orbitals upon hydroxide adsorption, suggesting that the electronic structure exhibits d9L Cu2+ behavior rather than the local electronic structure of a formal Cu3+.

Negative Dielectric Constant of Water at a Metal Interface.

Water polarizability at a metal interface plays an essential role in electrochemistry. We devise a classical molecular dynamics approach with an efficient description of metal polarization and a novel ac field method to measure the local dielectric response of interfacial water. Water adlayers next to the metal surface exhibit higher-than-bulk in-plane and negative out-of-plane dielectric constants, the latter corresponding physically to overscreening of the applied field. If we account for the gap region at the interface, the average out-of-plane dielectric constant is quite low (ε_{⊥}≈2), in agreement with reported measurements on confined thin films.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach.

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H2O-···-H2O), the free energy to transform an H3O+ cation into an H5O2+ cation, and the hydrogen bond free energy of protonated water clusters (H3O+-···-H2O and H5O2+-···-H2O). The computed equilibrium hydrogen bond free energy of H2O-···-H2O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H3O+ by roughly 1-2 kBT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H+ cation and H2O in solution were found to be stronger than between two H2O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.