First principles studies of the potential-induced lifting of the Au(100) surface reconstruction

Thermodynamic Cyclic Voltammograms Based on Ab Initio Calculations: Ag(111) in Halide-Containing Solutions

Cyclic voltammograms (CVs) are a central experimental tool for assessing the structure and activity of electrochemical interfaces. Based on a mean-field ansatz for the interface energetics under applied potential conditions, we here derive an ab initio thermodynamics approach to efficiently simulate thermodynamic CVs. All unknown parameters are determined from density functional theory (DFT) calculations coupled to an implicit solvent model. For the showcased CVs of Ag(111) electrodes in halide-anion-containing solutions, these simulations demonstrate the relevance of double-layer contributions to explain experimentally observed differences in peak shapes over the halide series. Only the appropriate account of interfacial charging allows us to capture the differences in equilibrium coverage and total electronic surface charge that cause the varying peak shapes. As a case in point, this analysis demonstrates that prominent features in CVs do not only derive from changes in adsorbate structure or coverage but can also be related to variations of the electrosorption valency. Such double-layer effects are proportional to adsorbate-induced changes in the work function and/or interfacial capacitance. They are thus especially pronounced for electronegative halides and other adsorbates that affect these interface properties. In addition, the analysis allows us to draw conclusions on how the possible inaccuracy of implicit solvation models can indirectly affect the accuracy of other predicted quantities such as CVs.

An ab initio study of electrochemical vs. electromechanical properties: the case of CO adsorbed on a Pt(111) surface

We have studied electrochemical vibrational and energy properties of CO/Pt(111) in the framework of periodic density functional theory (DFT) calculations. We have used a modified version of the previously developed Filhol-Neurock method to correct the unphysical contributions arising from homogeneous background countercharge in the case of thick metallic slabs. The stability of different CO adsorption sites on Pt(111) (Top, Bridge, Hcp, Fcc) has been studied at constant electric field. The energies are dominated by the surface dipole interaction with the external electric field: a strong positive electric field favors the surfaces with the lower dipole moment (that correspond to the ones with the lower coordination). The Stark tuning slope of the CO stretching frequency for a Top site was calculated for different surface coverages in very good agreement with both experimental and other theoretical results. Finally, we have performed an analysis of the origin of Stark shifts showing that the total Stark effect can be split into two competing components. The first one corresponds to the direct effect of charging on the C-O chemical bond: it is referred as an electrochemical effect. The second is the consequence of the surface dipole interaction with the applied electric field that modifies the C-O distance, inducing a change of the C-O force constant because of C-O bond anharmonicity: it is referred as the electromechanical effect. In the CO/Pt(111) case, the dominant contribution is electromechanical. The electrochemical contribution is very small because the electronic system involved in the surface charging is mostly non-bonding as analyzed by looking at the surface Fukui function.

Multiscale Modeling of Au-Island Ripening on Au(100)

We describe a multiscale modeling hierarchy for the particular case of Au-island ripening on Au(100). Starting at the microscopic scale, density functional theory was used to investigate a limited number of self-diffusion processes on perfect and imperfect Au(100) surfaces. The obtained structural and energetic information served as basis for optimizing a reactive forcefield (here ReaxFF), which afterwards was used to address the mesoscopic scale. Reactive force field simulations were performed to investigate more diffusion possibilities at a lower computational cost but with similar accuracy. Finally, we reached the macroscale by means of kinetic Monte Carlo (kMC) simulations. The reaction rates for the reaction process database used in the kMC simulations were generated using the reactive force field. Using this strategy, we simulated nucleation, aggregation, and fluctuation processes for monoatomic high islands on Au(100) and modeled their equilibrium shape structures. Finally, by calculating the step line tension at different temperatures, we were able to make a direct comparison with available experimental data.