Ab initio molecular dynamics simulations of the O2/Pt(1 1 1) interaction

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

Misconceptions in the Exploding Flask Demonstration Resolved through Students' Critical Thinking

As it connects to a large set of important fundamental ideas in chemistry and analytical techniques discussed in high school chemistry curricula, we review the exploding flask demonstration. In this demonstration, methanol vapor is catalytically oxidized by a Pt wire catalyst in an open container. The exothermicity of reactions occurring at the catalytic surface heats the metal to the extent that it glows. When restricting reactant and product gas flow, conditions may favor repetitive occurrence of a small explosion. We show how mass spectrometry and infrared spectroscopy allow for unravelling the chemical background of this demonstration and discuss various ideas on how to use it in a classroom setting to engage students' critical thinking about chemical research. Along the way, we show that two commonly published ideas about the chemical background of this demonstration are incorrect, and we suggest simple tests that may be performed in a high school setting either as an addition to the demonstration or as a student research project.

Water Adsorption and Insertion in MOF-5

The high surface areas and tunable properties of metal-organic frameworks (MOFs) make them attractive materials for applications in catalysis and the capture, storage, and separation of gases. Nevertheless, the limited stability of some MOFs under humid conditions remains a point of concern. Understanding the atomic-scale mechanisms associated with MOF hydrolysis will aid in the design of new compounds that are stable against water and other reactive species. Toward revealing these mechanisms, the present study employs van der Waals-augmented density functional theory, transition-state finding techniques, and thermodynamic integration to predict the thermodynamics and kinetics of water adsorption/insertion into the prototype compound, MOF-5. Adsorption and insertion energetics were evaluated as a function of water coverage, while accounting for the full periodicity of the MOF-5 crystal structure, that is, without resorting to cluster approximations or structural simplifications. The calculations suggest that the thermodynamics of MOF hydrolysis are coverage-dependent: water insertion into the framework becomes exothermic only after a sufficient number of H₂O molecules are coadsorbed in close proximity on a Zn-O cluster. Above this coverage threshold, the adsorbed water clusters facilitate facile water insertion via breaking of Zn-O bonds: the calculated free-energy barrier for insertion is very low, 0.17 eV at 0 K and 0.04 eV at 300 K. Our calculations provide a highly realistic description of the mechanisms underlying the hydrolysis of MOFs under humid working conditions.

Dynamics of O2 Chemisorption on a Flat Platinum Surface Probed by an Alignment-Controlled O2 Beam

O₂ adsorption on Pt surfaces is of great technological importance owing to its relevance to reactions for the purification of car exhaust gas and the oxygen reduction on fuel-cell electrodes. Although the O₂ /Pt(111) system has been investigated intensively, questions still remain concerning the origin of the low O₂ sticking probability and its unusual energy dependence. We herein clarify the alignment dependence of the initial sticking probability (S₀ ) using the single spin-rotational state-selected [(J,M)=(2,2)] O₂ beam. The results indicate that, at low translational energy (E₀ ) conditions, direct activated chemisorption occurs only when the O₂ axis is nearly parallel to the surface. At high energy conditions (E₀ >0.5 eV), however, S₀ for the parallel O₂ decreases with increasing E₀ while that of the perpendicular O₂ increases, accounting for the nearly energy-independent O₂ sticking probability determined previously by a non-state-resolved experiment.

Toward full simulation of the electrochemical oxygen reduction reaction on Pt using first-principles and kinetic calculations

First-principles molecular dynamics gave the kinetic and redox parameters of the oxygen reduction reaction in a fuel cell using a bias control scheme, and gave the current–voltage relationship.

Quantum and classical dynamics of reactive scattering of H2 from metal surfaces

State-of-the-art theoretical models allow nowadays an accurate description of H₂/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.