Sequential Selective Dissolution of Coinage Metals in Recyclable Ionic Media.

Coinage metals Cu, Ag, and Au are essential for modern electronics and their recycling from waste materials is becoming increasingly important to guarantee the security of their supply. Designing new sustainable and selective procedures that would substitute currently used processes is crucial. Here, we describe an unprecedented approach for the sequential dissolution of single metals from Cu, Ag, and Au mixtures using biomass-derived ionic solvents and green oxidants. First, Cu can be selectively dissolved in the presence of Ag and Au with a choline chloride/urea/H2O2 mixture, followed by the dissolution of Ag in lactic acid/H2O2. Finally, the metallic Au, which is not soluble in either solution above, is dissolved in choline chloride/urea/Oxone. Subsequently, the metals were simply and quantitatively recovered from dissolutions, and the solvents were recycled and reused. The applicability of the developed approach was demonstrated by recovering metals from electronic waste substrates such as printed circuit boards, gold fingers, and solar panels. The dissolution reactions and selectivity were explored with different analytical techniques and DFT calculations. We anticipate our approach will pave a new way for the contemporary and sustainable recycling of multi-metal waste substrates.

Cation-induced changes in the inner- and outer-sphere mechanisms of electrocatalytic CO2 reduction.

The underlying mechanism of cation effects on CO2RR remains debated. Herein, we study cation effects by simulating both outer-sphere electron transfer (OS-ET) and inner-sphere electron transfer (IS-ET) pathways during CO2RR via constrained density functional theory molecular dynamics (cDFT-MD) and slow-growth DFT-MD (SG-DFT-MD), respectively. Our results show without any cations, only OS-ET is feasible with a barrier of 1.21 eV. In the presence of K+ (Li+), OS-ET shows a very high barrier of 2.93 eV (4.15 eV) thus being prohibited. However, cations promote CO2 activation through IS-ET with the barrier of only 0.61 eV (K+) and 0.91 eV (Li+), generating the key intermediate (adsorbed CO[Formula: see text]). Without cations, CO2-to-CO[Formula: see text](ads) conversion cannot proceed. Our findings reveal cation effects arise from short-range Coulomb interactions with reaction intermediates. These results disclose that cations modulate the inner- and outer-sphere pathways of CO2RR, offering substantial insights on the cation specificity in the initial CO2RR steps.

Cations Determine the Mechanism and Selectivity of Alkaline Oxygen Reduction Reaction on Pt(111).

The proton-coupled electron transfer (PCET) mechanism of the oxygen reduction reaction (ORR) is a long-standing enigma in electrocatalysis. Despite decades of research, the factors determining the microscopic mechanism of ORR-PCET as a function of pH, electrolyte, and electrode potential remain unresolved, even on the prototypical Pt(111) surface. Herein, we integrate advanced experiments, simulations, and theory to uncover the mechanism of the cation effects on alkaline ORR on well-defined Pt(111). We unveil a dual-cation effect where cations simultaneously determine i) the active electrode surface by controlling the formation of Pt-O and Pt-OH overlayers and ii) the competition between inner- and outer-sphere PCET steps. The cation-dependent transition from Pt-O to Pt-OH determines the ORR mechanism, activity, and selectivity. These findings provide direct evidence that the electrolyte affects the ORR mechanism and performance, with important consequences for the practical design of electrochemical systems and computational catalyst screening studies. Our work highlights the importance of complementary insight from experiments and simulations to understand how different components of the electrochemical interface contribute to electrocatalytic processes.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats.

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

Influence of a Cu-zirconia interface structure on CO2 adsorption and activation.

CO2 adsorption and activation on a catalyst are key elementary steps for CO2 conversion to various valuable products. In the present computational study, we screened different Cu-ZrO2 interface structures and analyzed the influence of the interface structure on CO2 binding strength using density functional theory calculations. Our results demonstrate that a Cu nanorod favors one position on both tetragonal and monoclinic ZrO2 surfaces, where the bottom Cu atoms are placed close to the lattice oxygens. In agreement with previous calculations, we find that CO2 prefers a bent bidentate configuration at the Cu-ZrO2 interface and the molecule is clearly activated being negatively charged. Straining of the Cu nanorod influences CO2 adsorption energy but does not change the preferred nanorod position on zirconia. Altogether, our results highlight that CO2 adsorption and activation depend sensitively on the chemical composition and atomic structure of the interface used in the calculations. This structure sensitivity may potentially impact further catalytic steps and the overall computed reactivity profile.

Dynamics of weak interactions in the ligand layer of meta-mercaptobenzoic acid protected gold nanoclusters Au68(m-MBA)32 and Au144(m-MBA)40.

Atomically precise metal nanoclusters, stabilized and functionalized by organic ligands, are emerging nanomaterials with potential applications in plasmonics, nano-electronics, bio-imaging, nanocatalysis, and as therapeutic agents or drug carriers in nanomedicine. The ligand layer has an important role in modifying the physico-chemical properties of the clusters and in defining the interactions between the clusters and the environment. While this role is well recognized from a great deal of experimental studies, there is very little theoretical information on dynamical processes within the layer itself. Here, we have performed extensive molecular dynamics simulations, with forces calculated from the density functional theory, to investigate thermal stability and dynamics of the ligand layer of the meta-mercaptobenzoic acid (m-MBA) protected Au68 and Au144 nanoclusters, which are the first two gold nanoclusters structurally solved to atomic precision by electron microscopy [Azubel et al., Science, 2014, 345, 909 and ACS Nano, 2017, 11, 11866]. We visualize and analyze dynamics of three distinct non-covalent interactions, viz., ligand-ligand hydrogen bonding, metal-ligand O[double bond, length as m-dash]C-OHAu interaction, and metal-ligand Ph(π)Au interaction. We discuss their relevance for defining, at the same time, the dynamic stability and reactivity of the cluster. These interactions promote the possibility of ligand addition reactions for bio-functionalization or allow the protected cluster to act as a catalyst where active sites are dynamically accessible inside the ligand layer.

Towards Atomically Precise Supported Catalysts from Monolayer-Protected Clusters: The Critical Role of the Support.

Controlling the size and uniformity of metal clusters with atomic precision is essential for fine-tuning their catalytic properties, however for clusters deposited on supports, such control is challenging. Here, by combining X-ray absorption spectroscopy and density functional theory calculations, it is shown that supports play a crucial role in the evolution of monolayer-protected clusters into catalysts. Based on the acidic nature of the support, cluster-support interactions lead either to fragmentation of the cluster into isolated Au-ligand species or ligand-free metallic Au0 clusters. On Lewis acidic supports that bind metals strongly, the latter transformation occurs while preserving the original size of the metal cluster, as demonstrated for various Aun sizes. These findings underline the role of the support in the design of supported catalysts and represent an important step toward the synthesis of atomically precise supported nanomaterials with tailored physico-chemical properties.

Escaping scaling relationships for water dissociation at interfacial sites of zirconia-supported Rh and Pt clusters.

Water dissociation is an important reaction involved in many industrial processes. In this computational study, the dissociation of water is used as a model reaction for probing the activity of interfacial sites of globally optimized ZrO2 supported Pt and Rh clusters under the framework of density functional theory. Our findings demonstrate that the perimeter sites of these small clusters can activate water, but the dissociation behavior varies considerably between sites. It is shown that the studied clusters break scaling relationships for water dissociation, suggesting that these catalysts may achieve activities beyond the maximum imposed by such relations. Furthermore, we observed large differences in the thermodynamics of the water dissociation reaction between global minimum and near-global minimum isomers of the clusters. Overall, our results highlight the uniqueness of interfacial sites in catalytic reactions and the need for developing new concepts and tools to deal with the associated complexity.

Atomically Precise, Thiolated Copper-Hydride Nanoclusters as Single-Site Hydrogenation Catalysts for Ketones in Mild Conditions.

Copper-hydrides are known catalysts for several technologically important reactions such as hydrogenation of CO, hydroamination of alkenes and alkynes, and chemoselective hydrogenation of unsaturated ketones to unsaturated alcohols. Stabilizing copper-based particles by ligand chemistry to nanometer scale is an appealing route to make active catalysts with optimized material economy; however, it has been long believed that the ligand-metal interface, particularly if sulfur-containing thiols are used as stabilizing agent, may poison the catalyst. We report here a discovery of an ambient-stable thiolate-protected copper-hydride nanocluster [Cu25H10(SPhCl2)18]3- that readily catalyzes hydrogenation of ketones to alcohols in mild conditions. A full experimental and theoretical characterization of its atomic and electronic structure shows that the 10 hydrides are instrumental for the stability of the nanocluster and are in an active role being continuously consumed and replenished in the hydrogenation reaction. Density functional theory computations suggest, backed up by the experimental evidence, that the hydrogenation takes place only around a single site of the 10 hydride locations, rendering the [Cu25H10(SPhCl2)18]3- one of the first nanocatalysts whose structure and catalytic functions are characterized fully to atomic precision. Understanding of a working catalyst at the atomistic level helps to optimize its properties and provides fundamental insights into the controversial issue of how a stable, ligand-passivated, metal-containing nanocluster can be at the same time an active catalyst.

Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials.

Properties of solid-liquid interfaces are of immense importance for electrocatalytic and electrochemical systems, but modeling such interfaces at the atomic level presents a serious challenge and approaches beyond standard methodologies are needed. An atomistic computational scheme needs to treat at least part of the system quantum mechanically to describe adsorption and reactions, while the entire system is in thermal equilibrium. The experimentally relevant macroscopic control variables are temperature, electrode potential, and the choice of the solvent and ions, and these need to be explicitly included in the computational model as well; this calls for a thermodynamic ensemble with fixed ion and electrode potentials. In this work, a general framework within density functional theory (DFT) with fixed electron and ion chemical potentials in the grand canonical (GC) ensemble is established for modeling electrocatalytic and electrochemical interfaces. Starting from a fully quantum mechanical description of multi-component GC-DFT for nuclei and electrons, a systematic coarse-graining is employed to establish various computational schemes including (i) the combination of classical and electronic DFTs within the GC ensemble and (ii) on the simplest level a chemically and physically sound way to obtain various (modified) Poisson-Boltzmann (mPB) implicit solvent models. The detailed and rigorous derivation clearly establishes which approximations are needed for coarse-graining as well as highlights which details and interactions are omitted in vein of computational feasibility. The transparent approximations also allow removing some of the constraints and coarse-graining if needed. We implement various mPB models within a linear dielectric continuum in the GPAW code and test their capabilities to model capacitance of electrochemical interfaces as well as study different approaches for modeling partly periodic charged systems. Our rigorous and well-defined DFT coarse-graining scheme to continuum electrolytes highlights the inadequacy of current linear dielectric models for treating properties of the electrochemical interface.

DFT Prediction of Enhanced Reducibility of Monoclinic Zirconia upon Rhodium Deposition.

Oxides are an important class of materials and are widely used, for example, as supports in heterogeneous catalysis. In a number of industrial catalytic processes, oxide supports actively participate in chemical transformations by releasing lattice oxygen anions. While this is intuitively understood for reducible oxides, the reducibility of irreducible oxides may be modified via nanoengineering or upon inclusion of foreign species. Our calculations predict that the ability of irreducible monoclinic zirconia to release oxygen improves substantially upon deposition of rhodium. Through a comprehensive screening of Rh/ZrO2 with different size of the rhodium species, we find that a Rh adatom and a Rh4 nanocluster have the largest impact on the reducibility of zirconia. With increasing size the effect of rhodium decays. Our findings demonstrate that the phenomenon of enhanced reducibility of irreducible oxides in the presence of metals should be considered when interpreting experimental and computational results, as reactions that involve release of oxygen from an oxide support might be possible for irreducible oxides.

First Principles Calculations for Hydrogenation of Acrolein on Pd and Pt: Chemoselectivity Depends on Steric Effects on the Surface.

The chemoselective hydrogenation of acrolein on Pt(111) and Pd(111) surfaces is investigated employing density functional theory calculations. The computed potential energy surfaces together with the analysis of reaction mechanisms demonstrate that steric effects are an important factor that governs chemoselectivity. The reactions at the C=O functionality require more space than the reactions at the C=C functionality. Therefore the formation of allyl alcohol is more favorable at low coverage, while the reduction of the C=C bond and the formation of propanal becomes kinetically more favorable at higher coverage. The elementary reaction steps are found to follow different reaction mechanisms, which are identified according to terminology typically used in organometallic catalysis. The transition state scaling (TSS) relationship is demonstrated and the origin of multiple TSS lines is linked to variation of an internal electronic structure of a carbon skeleton.

Carbon dioxide activation and reaction induced by electron transfer at an oxide-metal interface.

A model system has been created to shuttle electrons through a metal-insulator-metal (MIM) structure to induce the formation of a CO2 anion radical from adsorbed gas-phase carbon dioxide that subsequently reacts to form an oxalate species. The process is completely reversible, and thus allows the elementary steps involved to be studied at the atomic level. The oxalate species at the MIM interface have been identified locally by scanning tunneling microscopy, chemically by IR spectroscopy, and their formation verified by density functional calculations.

Water and carbon oxides on monoclinic zirconia: experimental and computational insights.

Zirconium oxide (ZrO2, zirconia) is an interesting catalytic material to be used in biomass conversion, e.g., gasification and reforming. In this work, we show that reducing and hydrating pretreatments affect the surface sites on monoclinic zirconia. The multitechnique approach comprises temperature-programmed surface reactions (TPSR) under CO and CO2 at 100-550 °C, in situ DRIFTS investigations of the surface species and density functional theory (DFT) calculations. The key findings of the work are: (1) formates are formed either directly from gas-phase CO on terminal surface hydroxyls or via the linear CO surface species that are found exclusively on the reduced zirconia without water treatment; (2) formates are able to decompose at high temperature either reversibly to CO or reductively to CO2 and H2via surface reaction between formates and multicoordinated hydroxyls; and (3) a new weak reversible binding state of CO is found exclusively on ZrO2 that is first reduced and subsequently hydrated.

Electronic and vibrational properties of meso-tetraphenylporphyrin on silver substrates.

The electronic and vibrational properties of meso-tetraphenylporphyrin (mtpp) on silver substrates are investigated using UV-vis and surface-enhanced resonance Raman scattering (SERRS) spectroscopy. Whereas the vibrational signatures associated with the tetrapyrrole backbone exhibit minor variations throughout sequences of consecutively recorded SERRS spectra, the C═C stretching vibrational modes localized on the meso-phenyl moieties of mtpp exhibit noticeable intensity fluctuations, masked in the average SERRS response. We attribute the observed vibrational-state-specific blinking events to conformational changes in mtpp, namely, torsional flexibility which mediates the coupling between the π-framework of the meso-phenyls and the underlying metal substrate.

Orbiting orbitals: visualization of vibronic motion at a conical intersection.

The Jahn-Teller (JT) active unpaired electron of single metalloporphyrin radical anions is imaged through scanning tunneling microscopy. It is demonstrated that the electron is delocalized over the porphyrin macrocycle and its topographic image is determined by vibronic motion: the orbital of the electron adiabatically follows the zero-point pseudorotation of skeletal deformations. Transformation of the polar graphs of the observed images allows visualization of the adiabatic vibrational density to which the electron is coupled. The vibronic potential at the conical intersection is visualized and the half-integer angular momentum characteristic of the Berry phase is revealed in the radial function of the electron. The measurements underscore the economy of Jahn-Teller dynamics: small atomic displacements (∼10(-1) Å) determined by weak interactions (∼10 meV) control the motion of the electron on a 10 Å scale and determine the molecular shape and function.

Thiolate adsorption on Au(hkl) and equilibrium shape of large thiolate-covered gold nanoparticles.

The adsorption of thiolates on Au surfaces employing density-functional-theory calculations has been studied. The dissociative chemisorption of dimethyl disulfide (CH(3)S-SCH(3)) on 14 different Au(hkl) is used as a model system. We discuss trends on adsorption energies, bond lengths, and bond angles as the surface structure changes, considering every possible Au(hkl) with h, k, l ≤ 3 plus the kinked Au(421). Methanethiolate (CH(3)S-) prefers adsorption on bridge sites on all surfaces considered; hollow and on top sites are highly unfavourable. The interface tensions for Au(hkl)-thiolate interfaces is determined at low coverage. Using the interface tensions in a Wulff construction method, we construct atomistic models for the equilibrium shape of large thiolate-covered gold nanoparticles. Gold atoms in a nanoparticle change their equilibrium positions upon adsorption of thiolates towards shapes of higher sphericity and higher concentration of step-edge atoms.

First principles investigations of Pd-on-Au nanostructures for trichloroethene catalytic removal from groundwater.

Catalytic groundwater remediation from chlorinated organic solvents like trichloroethene (TCE) has been found to be more effective and sustainable than traditional non-destructive methods. Among the experimentally studied catalyst materials, Pd-decorated Au nanoparticles show the highest activity and selectivity combined with the best resistance towards poisoning by chemicals present in groundwater. In this study the thermochemistry and adsorption geometries of TCE and its hydrodechlorination products are investigated via density functional theory calculations. Various model systems for Pd-supported Au nanoparticles are addressed. The adsorption of TCE is endothermic on bare Au(111), almost thermoneutral or slightly exothermic on Pd-Au surface alloys and clearly exothermic on Pd overlayer structures on Au(111). The strongest chemisorption is on the di-σ configuration between Pd atoms over the smallest 2D Pd clusters containing only a few Pd atoms. These are not, however, the best catalysts as they are too small to co-adsorb hydrogen needed for hydrodechlorination reaction. We demonstrate good correlation between adsorption energies and the d-band center of the system. The variation of adsorption energy from the one Pd-Au composition to the other can be tentatively assigned to be due to the ligand and coordination effects. Also, the ensemble effects are important; without the right ensemble the adsorption is weak or endothermic.