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 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 recycled and reused. The applicability of 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 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.

Ligand assisted hydrogenation of levulinic acid on Pt(111) from first principles calculations

In this study, we investigate the hydrogenation reaction of levulinic acid to 4-hydroxypentanovic acid on a ligand-modified Pt(111) using DFT. Modifying nanoparticle surfaces with ligands can have beneficial effects on...

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

CO₂ adsorption and activation on a catalyst are key elementary steps for CO₂ conversion to various valuable products. In the present computational study, we screened different Cu-ZrO₂ interface structures and analyzed the influence of the interface structure on CO₂ binding strength using density functional theory calculations. Our results demonstrate that a Cu nanorod favors one position on both tetragonal and monoclinic ZrO₂ surfaces, where the bottom Cu atoms are placed close to the lattice oxygens. In agreement with previous calculations, we find that CO₂ prefers a bent bidentate configuration at the Cu-ZrO₂ interface and the molecule is clearly activated being negatively charged. Straining of the Cu nanorod influences CO₂ adsorption energy but does not change the preferred nanorod position on zirconia. Altogether, our results highlight that CO₂ 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.

Chemoselective heterogeneous iridium catalyzed hydrogenation of cinnamalaniline

A selective atmospheric hydrogenation of unsaturated imines over heterogeneous iridium catalyst is described, in addition the selectivity is elucidated by DFT-calculations.

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.

Binding Behavior of Carbonmonoxide to Gold Atoms on Ag(001)

The adsorption behavior of single CO molecules at 4 K bound to Au adatoms on a Ag(001) metal surface is studied with scanning tunneling microscopy (STM) and inelastic electron tunneling spectroscopy (IETS). In contrast to earlier observations two different binding configurations are observed—one on top of a Au adatom and the other one adsorbed laterally to Au on Ag(001). Moreover, IETS reveals different low-energy vibrational energies for the two binding sites as compared to the one for a single CO molecule bound to Ag(001). Density functional theory (DFT) calculations of the adsorption energies, the diffusion barriers, and the vibrational frequencies of the CO molecule on the different binding sites rationalize the experimental findings.

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 Au⁰ 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 Au_n 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.

First-principles insight into CO hindered agglomeration of Rh and Pt single atoms on m-ZrO2

Kinetic and thermodynamic stability of single-atom and nanocluster catalysts is addressed under reaction conditions within a DFT-parametrised multi-scale thermodynamic framework combining atomistic, non-equilibrium, and nanothermodynamics.

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 [Cu₂₅H₁₀(SPhCl₂)₁₈]^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 [Cu₂₅H₁₀(SPhCl₂)₁₈]^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/ZrO₂ with different size of the rhodium species, we find that a Rh adatom and a Rh₄ 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.

Gold assisted oxygen dissociation on a molybdenum-doped CaO(001) surface

Using density functional theory (DFT) calculations, we address the adsorption of O₂ and the coadsorption of gold species and oxygen molecules on a Mo-doped CaO(001) surface with 1.25% impurity concentration.

Review: monoclinic zirconia, its surface sites and their interaction with carbon monoxide

This review is focused on the surface sites of monoclinic zirconia, especially surface hydroxyls and their manipulation, and the interaction of CO with zirconia to form linearly adsorbed CO and formates.

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

Formates form via gas-phase or linearly adsorbed CO and decompose either reversibly to CO or reductively to CO₂ and H₂.

Quantum well states in two-dimensional gold clusters on MgO thin films

The electronic structure of ultrasmall Au clusters on thin MgO/Ag(001) films has been analyzed by scanning tunneling spectroscopy and density functional theory. The clusters exhibit two-dimensional quantum well states, whose shapes resemble the eigenstates of a 2D electron gas confined in a parabolic potential. From the symmetry of the highest occupied (HOMO) and lowest unoccupied molecular orbital (LUMO) of a particular cluster, its electron filling and charge state is determined. In accordance with a Bader charge analysis, aggregates containing up to 20 atoms accumulate one to four extra electrons due to a charge transfer from the MgO/Ag interface. The HOMO-LUMO gap is found to close for clusters containing between 70 and 100 atoms.

Electronic structure of MgO-supported Au clusters: quantum dots probed by scanning tunneling microscopy

We investigate via density functional theory (DFT) the appearance of small MgO-supported gold clusters with 8 to 20 atoms in a scanning tunneling microscope (STM) experiment. Comparison of simulations of ultrathin films on a metal support with a bulk MgO leads to similar results for the cluster properties relevant for STM. Simulated STM pictures show the delocalized states of the cluster rather than the atomic structure. This finding is due to the presence of s- derived delocalized states of the cluster near the Fermi energy. The properties of theses states can be understood from a jellium model for monovalent gold.