Search

Based on the results of periodic density functional theory calculations, we have recently proposed that the height of self-diffusion barriers can serve as a descriptor for dendrite growth in batteries [M. Jackle et al., Energy Environ. Sci. 11, 3400 (2018)]. However, in the determination of the self-diffusion barriers, the electrochemical environment has not been taken into account. Still, due to the presence of electrical double layers at electrode/electrolyte interfaces, strong electric fields can be present close to the interfacial region. In a first step toward including the electrochemical environment, we have calculated barriers for terrace-diffusion on lithium, magnesium, and silver surfaces and across-step self-diffusion on lithium in the presence of electric fields. Whereas the electric field effect is more pronounced on a stepped surface than on flat terraces, overall we find a negligible influence of electric fields on self-diffusion barriers which we explain by the good screening properties of metals.

Ionic liquids (ILs) are considered as attractive electrolyte solvents in modern battery concepts such as Li-ion batteries. Here we present a comprehensive review of the results of previous model studies on the interaction of the battery relevant IL 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP]

The description of electrode-electrolyte interfaces is based on the concept of the formation of an electric double layer. This concept was derived from continuum theories extended by introducing point charge distributions. Based on ab initio molecular dynamics simulations, we analyze the electric double layer in an approach beyond the point charge scheme by instead assessing charge polarizations at electrochemical metal-water interfaces from first principles. We show that the atomic structure of water layers at room temperature leads to an oscillatory behavior of the averaged electrostatic potential. We address the relation between the polarization distribution at the interface and the extent of the electric double layer and subsequently derive the electrode potential from the charge polarization.

The structure of a liquid water layer on Pt(111) has been studied by ab initio molecular dynamics simulations based on periodic density functional theory calculations. First the reliability of the chosen exchange-correlation function has been validated by considering water clusters, bulk ice structures, and bulk liquid water, confirming that the dispersion corrected RPBE-D3/zero functional is a suitable choice. The simulations at room temperature yield that a water layer that is six layers thick is sufficient to yield liquid water properties in the interior of the water film. Performing a statistical average along the trajectory, a mean work function of 5.01 V is derived, giving a potential of zero charge of Pt(111) of 0.57 V vs. standard hydrogen electrode, in good agreement with experiments. Therefore we propose the RPBE-D3/zero functional as the appropriate choice for first-principles calculations addressing electrochemical aqueous electrolyte/metal electrode interfaces.

Ionic liquids (ILs) are considered as attractive electrolyte solvents in modern battery concepts such as Li-ion batteries. Here we present a comprehensive review of the results of previous model studies on the interaction of the battery relevant IL 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP]+[TFSI]−) with a series of structurally and chemically well-defined model electrode surfaces, which are increasingly complex and relevant for battery applications [Ag(111), Au(111), Cu(111), pristine and lithiated highly oriented pyrolytic graphite (HOPG), and rutile TiO2(110)]. Combining surface science techniques such as high resolution scanning tunneling microscopy and X-ray photoelectron spectroscopy for characterizing surface structure and chemical composition in deposited (sub-)monolayer adlayers with dispersion corrected density functional theory based calculations, this work aims at a molecular scale understanding of the fundamental processes at the electrode | electrolyte interface, ...

Lithium and magnesium exhibit rather different properties as battery anode materials with respect to the phenomenon of dendrite formation which can lead to short-circuits in batteries. Diffusion processes are the key to understanding structure forming processes on surfaces. Therefore, we have determined adsorption energies and barriers for the self-diffusion on Li and Mg using periodic density functional theory calculations and contrasted the results to Na which is also regarded as a promising electrode material in batteries. According to our calculations, magnesium exhibits a tendency towards the growth of smooth surfaces as it exhibits lower diffusion barriers than lithium and sodium, and as an hcp metal it favors higher-coordinated configurations in contrast to the bcc metals Li and Na. These characteristic differences are expected to contribute to the unequal tendencies of these metals with respect to dendrite growth.

We have treated the ionization probability of iodine molecules scattered from diamond by a semiclassical surface hopping scheme, namely Tully’s fewest-switches algorithm [J. Chem. Phys. 93, 1061 (1990)]. The interaction is described by a model potential that has been adjusted to empirical data. We start with a one-dimensional two-state model in which just the molecular distance from the surface and the neutral and negatively charged state of I2 are considered. We determine the ionization probability within the adiabatic and diabatic representation and compare it with exact quantum calculations. For this particular problem we find that the diabatic picture shows too little coherence, while the adiabatic representation yields satisfactory results. In the second part we have successively increased the complexity of the simulation by additionally taking a surface oscillator coordinate, the molecular rotation and vibration into account. Including more degrees of freedom damps out the Stuckelberg oscillations p...

Three-dimensional quantum and classical dynamical calculations of the dissociative adsorption of hydrogen have been performed, in which, besides one reaction path coordinate, the lateral degrees of freedom of the hydrogen center of mass were taken into account. These calculations were compared to results obtained by classical and quantum sudden approximations, which assessed the importance of tunneling, zero-point effects, and also steering in the dissociative adsorption of hydrogen. For energies below the minimum barrier height, tunneling is of course the relevant mechanism for dissociation, but above the minimum barrier height quantization and zero-point effects become more prominent. Zero-point effects suppress the dissociation probability; however, for energies slightly above the minimum barrier height, steering of the particles is only operative in the quantum dynamics and can thereby almost compensate the suppression of the quantum sticking probabilities due to zero-point effects, compared to the cl...

The adsorption dynamics of water on Pt(111) was studied using ab initio molecular dynamics simulations based on density functional theory calculations including dispersion corrections. Sticking probabilities were derived as a function of initial kinetic energy and water coverage. In addition, the energy distribution upon adsorption was monitored in order to analyze the energy dissipation process. We find that on the water pre-covered surface the sticking probability is enhanced because of the attractive water-water interaction and the additional effective energy dissipation channels to the adsorbed water molecules. The water structures forming directly after the adsorption on the pre-covered surfaces do not necessarily correspond to energy minimum structures.

The influence of one‐dimensional lateral surface corrugation on the dissociative adsorption and associative desorption has been investigated by quantum dynamical model calculations. The study includes potential energy surfaces where either the adsorption is hindered by high barriers as in the system H2/Cu(111), or where activated as well as nonactivated paths to adsorption exist, which should be relevant for the description of, e.g., H2/Pd(100). In the high‐barrier system the dynamics for normal energies below the minimum barrier height is determined by the minimum energy path, whereas at larger energies the specific distribution in the barrier height and location has a strong influence on the adsorption and desorption dynamics. At surfaces where also nonactivated paths to adsorption are present, strong quantum effects due to the wave nature of the hydrogen beam are found although sticking at such surfaces is classically possible at all incident energies. It is shown that effects usually attributed to the...

The interaction of hydrogen with palladium surfaces represents one of the model systems for the study of the adsorption and absorption at metal surfaces. Theoretical gas-surface dynamics studies have usually concentrated on the adsorption dynamics on clean surfaces. Only recently it has become possible, based on advances in the electronic structure codes and improvements in the computer power, to address the much more complex problem of the adsorption dynamics on precovered surfaces. Here, I present ab initio molecular dynamics (AIMD) simulations based on periodic density functional theory (DFT) calculations of the adsorption of H(2) on hydrogen-precovered Pd(100) for a broad variety of different hydrogen coverage structures. The stability of the adsorbate structures and the adsorption dynamics are analyzed in detail. Calculated sticking probabilities are larger than expected for pure site-blocking consistent with experimental results. It turns out that the adsorption dynamics on the strongly corrugated surfaces depends sensitively on the dynamic response of the substrate atoms upon the impact of the impinging H(2) molecules. In addition, for some structures the adsorption probability was evaluated as a function of the kinetic energy. Adsorbate structures corresponding to the same coverage but with different arrangements of the adsorbed atoms can lead to a qualitatively different dependence of the adsorption probability on the kinetic energy changing also the order of the preferred structures, as far as the adsorption is concerned, as a function of the kinetic energy. This indicates that dynamical effects such as steering and dynamical trapping play an important role in the adsorption on these precovered substrates.

The adsorption of benzene, thiophene, and pyridine on the (111) surface of gold and copper have been studied using density functional theory (DFT). Adsorption geometries and energies as well as the nature of bonding have been analyzed and compared to experimental results. Dispersion effects between neighboring molecules and between molecules and the surface have been taken into account via a semiempirical C(6)R(-6) approach. The C(6) coefficients for metal atoms have been deduced using both atomic properties and a hybrid QM:QM approach. Whereas the pure DFT calculations underestimate the adsorption energies significantly, a good agreement with experimental results is obtained using the DFT-D method based on the QM:QM hybrid approach.

We present a computational study of the interface of a Pt electrode and an aqueous electrolyte employing semi-empirical dispersion corrections and an implicit solvent model within first-principles calculations. The electrode potential is parametrized within the computational hydrogen electrode scheme. Using one explicit layer, we find that the most realistic interface configuration is a water bilayer in the H-up configuration. Furthermore, we focus on the contribution of the dispersion interaction and the presence of water on H, O, and OH adsorption energies. This study demonstrates that the implicit water scheme represents a computationally efficient method to take the presence of an aqueous electrolyte interface with a metal electrode into account.

We report the first six-dimensional quantum dynamical study of the dissociative adsorption of H(2) on a (110) surface. We have performed quantum coupled-channel calculations for the system H(2)/Pd(110) based on a potential energy surface (PES) that was derived from ab initio electronic structure calculations. In particular, we have focused on the effects of the corrugation and anisotropy of the PES on the H(2) dissociation probability. Our results agree well with the available experimental data for the sticking probability as a function of the initial kinetic energy and the angle of incidence. Because of the coupling between the anisotropy and corrugation of the potential energy surface our calculations predict an unusual rotational heating and a rather small rotational alignment in desorption.

The structure of liquid water has been addressed by ab initio molecular dynamics simulations based on density functional theory. Exchange-correlation effects have been described by the popular PBE and RPBE functionals within the generalized gradient approximation as these functionals also yield satisfactory results for metals which is important to model electrochemical interfaces from first principles. In addition, dispersive interactions are included by using dispersion-corrected schemes. It turns out that the dispersion-corrected RPBE functional reproduces liquid water properties quite well in contrast to the PBE functional. This is caused by the replacement of the over-estimated directional hydrogen-bonding in the PBE functional by non-directional dispersive interactions.

We have studied the interaction of CH(4) with Al(2) and Al(3) neutral and charged clusters in the two lowest lying spin states using density functional theory. These calculations, via extended search, are used to determine the stable positions of H and CH(3) near the cluster, and the transition state to break the H-CH(3) bond. In all cases, stable methyl-aluminum-hydrides are possible. The H desorption is studied by means of vibration analysis and application of transition state theory. A common observed trend is that, in breaking the H-CH(3) bond, the interacting H atom is attached to the "surface" of the clusters attracting some negative charge of approximately 0.2e. The charge transfer is illustrated using the corresponding orbitals near the transition state in conjunction with the computed Mulliken population analysis. Thermal vibrations, generally, do not enhance the reaction. In all exothermic cases, the binding energy toward CH(3)+HAl(n) (charge) increases with increasing charge of the original Al(n) ((q=-1,0,1)) cluster. Although Al lacks occupied d-orbitals, the small Al clusters reduce the (free methane) CH(3)-H dissociation barrier except for Al(3) ((q=-1,0)). The relevant reactions in desorption require approximately 400-700 degrees C.

The adsorption, dissociation, diffusion, and desorption of oxygen interacting with the Pt(111) surface have been studied using kinetic Monte Carlo simulations. This study has been motivated by uncertainties in the theoretical and the experimental derivations of O(2)Pt(111) reaction barriers. The simulations reproduce all known experimental data within basically one set of parameters, thus yielding microscopic insights into the elementary reaction steps occurring in the interaction of oxygen with Pt(111) and providing reliable estimates for adsorption energies and diffusion and desorption barriers. In particular, we confirm that the distance of oxygen atoms directly after dissociation is caused by ballistic hot atom motion rather than by diffusive motion. We address the equilibrium structure of oxygen atoms at high coverages. At low temperatures, chains of oxygen pairs are formed. We show that this mechanism can be explained by a lowered dissociation in the vicinity of already adsorbed atoms. Finally we discuss the role of the lateral interaction between the oxygen atoms in the oxygen desorption process.

The sticking and scattering of O(2)Pt(111) has been studied by tight-binding molecular dynamics simulations based on an ab initio potential energy surface. We focus, in particular, on the sticking probability as a function of the angle of incidence and the energy and angular distributions in scattering. Our simulations provide an explanation for the seemingly paradox experimental findings that adsorption experiments suggest that the O(2)Pt(111) interaction potential should be strongly corrugated while scattering experiments indicate a rather small corrugation. The potential energy surface is indeed strongly corrugated which leads to a pronounced dependence of the sticking probability on the angle of incidence. The scattered O(2) molecules, however, experience a rather flat surface due to the fact that they are predominantly scattered at the repulsive tail of the potential.