Impact of Surface Ligand Identity and Density on the Thermodynamics of H Atom Uptake at Polyoxovanadate-Alkoxide Surfaces

An understanding of how molecular structure influences the thermodynamics of H atom transfer is critical to designing efficient catalysts for reductive chemistries. Herein, we report experimental and theoretical investigations summarizing structure-function relationships of polyoxovanadate-alkoxides that influence bond dissociation free energies of hydroxide ligands located at the surface of the cluster. We evaluate the thermochemical descriptors of O-H bond strength for a series of clusters, namely [V6O13-x(OH)x(TRIOLR)2]-2 (x = 2, 4, 6; R = NO2, Me) and [V6O11-x(OMe)2(OH)x(TRIOLNO2)2]-2, via computational analysis and open circuit potential measurements. Our findings reveal that modifications to the TRIOL ligand (e.g., changing from the previously reported electron withdrawing nitro-backed ligand to the electron-donating methyl variant) have limited influence on the strength of surface O-H bonds as a result of near complete thermodynamic compensation in these systems (i.e., correlated changes in redox potential and cluster basicity). In contrast, changes in surface density of alkoxide ligands via direct alkoxylation of the polyoxovanadate-alkoxide surface result in measurable increases in bond dissociation free energies of surface O-H bonds for the mixed-valent derivatives. Our findings indicate that the extent of (de)localization of electron density across the cluster core has an impact on the bond dissociation free energies of surface O-H bonds across all oxidation states of the assembly.

Participation of electrochemically inserted protons in the hydrogen evolution reaction on tungsten oxides

Understanding the mechanisms by which electrodes undergo the hydrogen evolution reaction (HER) is necessary to design better materials for aqueous energy storage and conversion. Here, we investigate the HER mechanism on tungsten oxide electrodes, which are stable in acidic electrolytes and can undergo proton-insertion coupled electron transfer concomitant with the HER. Electrochemical characterization showed that anhydrous and hydrated tungsten oxides undergo changes in HER activity coincident with changes in proton composition, with activity in the order HxWO3·H2O > HxWO3 > HxWO3·2H2O. We used operando X-ray diffraction and density functional theory to understand the structural and electronic changes in the materials at high states of proton insertion, when the oxides are most active towards the HER. H0.69WO3·H2O and H0.65WO3 have similar proton composition, structural symmetry, and electronic properties at the onset of the HER, yet exhibit different activity. We hypothesize that the electrochemically inserted protons can diffuse in hydrogen bronzes and participate in the HER. This would render the oxide volume, and not just the surface, as a proton and electron reservoir at high overpotentials. HER activity is highest in HxWO3·H2O, which optimizes both the degree of proton insertion and solid-state proton transport kinetics. Our results highlight the interplay between the HER and proton insertion-coupled electron transfer on transition metal oxides, many of which are non-blocking electrodes towards protons.

Single Atom Alloys Segregation in the Presence of Ligands

Single atom alloys (SAAs) have gained remarkable attention due to their tunable properties leading to enhanced catalytic performance, such as high activity and selectivity. The stability of SAAs is dictated by surface segregation, which can be affected by the presence of surface adsorbates. Research efforts have primarily focused on the effect of commonly found catalytic reaction intermediates, such as CO and H, on the stability of SAAs. However, there is a knowledge gap in understanding the effect of ligands from colloidal nanoparticle (NP) synthesis on surface segregation. Herein, we combine density functional theory (DFT) and machine learning to investigate the effect of thiol and amine ligands on the stability of colloidal SAAs. DFT calculations revealed rich segregation energy (Eseg) data of SAAs with d8 (Pt, Pd, Ni) and d9 (Ag, Au, Cu) metals exposing (111) and (100) surfaces, in the presence and absence of ligands. Using these data, we developed an accurate four-feature neural network using a multilayer perceptron regression (NN MLP) model. The model captures the underlying physics behind surface segregation in the presence of adsorbed ligands by incorporating features representing the thermodynamic stability of metals through the bulk cohesive energy, structural effects using the coordination number of the dopant and the ligands, the binding strength of the adsorbate to the metals, strain effects using the Wigner-Seitz radius, and electronic effects through electron affinities. We found that the presence of ligands makes the thermodynamics of segregation milder compared to the bare (nonligated) SAA surfaces. Importantly, the adsorption configuration (e.g., top vs bridge) and the binding strength of the ligand to the SAA surface (e.g., amines vs thiols) play an important role in altering the Eseg trends compared to the bare surface. We also developed an accurate NN MLP model that predicts Eseg in the presence of ligands to find thermodynamically stable SAAs, leading to the rapid and efficient screening of colloidal SAAs. Our model captures several experimental observations and elucidates complex physics governing segregation at nanoscale interfaces.

The Role of Grain Boundary Sites for the Oxidation of Copper Catalysts during the CO Oxidation Reaction

The oxidation of transition metal surfaces is a process that takes place readily at ambient conditions and that, depending on the specific catalytic reaction at hand, can either boost or hamper activity and selectivity. Cu catalysts are no exception in this respect since they exhibit different oxidation states for which contradicting activities have been reported, as, for example, in the catalytic oxidation of CO. Here, we investigate the impact of low-coordination sites on nanofabricated Cu nanoparticles with engineered grain boundaries on the oxidation of the Cu surface under CO oxidation reaction conditions. Combining multiplexed in situ single particle plasmonic nanoimaging, ex situ transmission electron microscopy imaging, and density functional theory calculations reveals a distinct dependence of particle oxidation rate on grain boundary density. Additionally, we found that the oxide predominantly nucleates at grain boundary-surface intersections, which leads to nonuniform oxide growth that suppresses Kirkendall-void formation. The oxide nucleation rate on Cu metal catalysts was revealed to be an interplay of surface coordination and CO oxidation behavior, with low coordination favoring Cu oxidation and high coordination favoring CO oxidation. These findings explain the observed single particle-specific onset of Cu oxidation as being the consequence of the individual particle grain structure and provide an explanation for widely distributed activity states of particles in catalyst bed ensembles.

Global and Local Connectivities Describe Hydrogen Intercalation in Metal Oxides

We introduce global connectivity, Ω, which quantifies the detailed connectivity of a material and accurately captures reactivity behavior. We demonstrate that combining global and local connectivity describes how metal oxides interact with hydrogen. Blending density functional theory, graph theory, and machine learning we built a reactivity model which accurately predicts hydrogen intercalation potentials of different metal oxides experimentally measured in the lab. The use of global connectivity can accelerate materials design through the development of novel structure-property relationships.

Confined Ru Sites in a 13X Zeolite for Ultrahigh H2 Production from NH3 Decomposition

Catalytic NH₃ synthesis and decomposition offer a new promising way to store and transport renewable energy in the form of NH₃ from remote or offshore sites to industrial plants. To use NH₃ as a hydrogen carrier, it is important to understand the catalytic functionality of NH₃ decomposition reactions at an atomic level. Here, we report for the first time that Ru species confined in a 13X zeolite cavity display the highest specific catalytic activity of over 4000 h^-1 for the NH₃ decomposition with a lower activation barrier, compared to most reported catalytic materials in the literature. Mechanistic and modeling studies clearly indicate that the N-H bond of NH₃ is ruptured heterolytically by the frustrated Lewis pair of Ru^δ+-O^δ- in the zeolite identified by synchrotron X-rays and neutron powder diffraction with Rietveld refinement as well as other characterization techniques including solid-state nuclear magnetic resonance spectroscopy, in situ diffuse reflectance infrared transform spectroscopy, and temperature-programmed analysis. This contrasts with the homolytic cleavage of N-H displayed by metal nanoparticles. Our work reveals the unprecedented unique behavior of cooperative frustrated Lewis pairs created by the metal species on the internal zeolite surface, resulting in a dynamic hydrogen shuttling from NH₃ to regenerate framework Brønsted acid sites that eventually are converted to molecular hydrogen.

Ab Initio Thermochemistry of Highly Flexible Molecules for Thermal Decomposition Analysis

Pyrolysis is a promising technology for chemical recycling of waste plastics, since it enables the generation of high-value chemicals with low capital and operating cost. The calculation of thermodynamic equilibrium composition using the Gibbs free energy minimization approach can determine pyrolysis operating conditions that produce desired products. However, the availability of thermochemical data can limit the application of equilibrium calculations. While density functional theory (DFT) calculations have been commonly used to produce accurate thermochemical data (e.g., enthalpies of formation) of small molecules, the accuracy and computational cost of these calculations are both challenging to handle for large, flexible molecules, exhibiting multiple conformations at elevated (i.e., pyrolysis) temperatures. In this work, we develop a computational framework to calculate accurate, temperature-dependent thermochemistry of large and flexible molecules by combining force field based conformational search, DFT calculations, thermochemical corrections, and Boltzmann statistics. Our framework produces accurately calculated thermochemistry that is used to predict equilibrium thermal decomposition profiles of octadecane, a model compound of polyethylene. Our thermochemistry results are compared against literature data demonstrating a great agreement, and the predicted decomposition profiles rationalize a series of pyrolysis experimental observations. Our work systematically addresses entropic contributions of large molecules and suggests paths for accurate and yet computationally feasible calculations of Gibbs free energies. The first-principles-based thermodynamic equilibrium analysis proposed in this work can be a significant step toward predicting temperature-dependent product distributions from plastic pyrolysis and guide experimentation on chemical plastic recycling.

Understanding formation thermodynamics of structurally diverse zeolite oligomers with first principles calculations

The mechanisms of many zeolitic processes, including nucleation and interzeolite transformation, are not fully understood owing to complex growth mixtures that obfuscate in situ monitoring of molecular events. In this work, we provide insights into zeolite chemistry by investigating the formation thermodynamics of small zeolitic species using first principles calculations. We systematically study how formation energies of pure-silicate and aluminosilicate species differ by structure type and size, temperature, and the presence of alkali or alkaline earth metal cations (Na⁺, K⁺, and Ca²⁺). Highly condensed (cage-like) species are found to be strongly preferred to simple rings in the pure-silicate system, and this thermodynamic preference increases with temperature. Introducing aluminum leads to more favorable formation thermodynamics for all species. Moreover, for species with a low Si/Al ratio (≤2), a thermodynamic preference does not exist among structure types; instead, a pool of diverse aluminosilicate structures compete in formation. Metal cation effects strongly depend on the presence of aluminum, cage size, cation type, and location, since each of these factors can alter electrostatic interactions between cations and zeolitic species. We reveal that confined metal cations may destabilize pure-silicate cages due to localized interactions; conversely, they stabilize aluminosilicates due to strong cation-framework attractions in sufficiently large cages. Importantly, this work rationalizes a series of experimental observations and can potentially guide efforts for controlling zeolite nucleation/crystallization processes.

Demystifying the Chemical Ordering of Multimetallic Nanoparticles

ConspectusMultimetallic nanoparticles (NPs) have highly tunable properties due to the synergy between the different metals and the wide variety of NP structural parameters such as size, shape, composition, and chemical ordering. The major problem with studying multimetallic NPs is that as the number of different metals increases, the number of possible chemical orderings (placements of different metals) for a NP of fixed size explodes. Thus, it becomes infeasible to explore NP energetic differences with highly accurate computational methods, such as density functional theory (DFT), which has a high computational cost and is typically applied to up to a couple of hundred metal atoms. Here, we demonstrate a methodology advancing NP simulations by effectively exploring the vast materials space of multimetallic NPs and accurately identifying the ones with the most thermodynamically preferred chemical orderings. With accuracies reaching that of DFT, our methodology is applicable to practically any NP size, shape, and metal composition. We achieve this by significantly advancing the bond-centric (BC) model, a physics-based model that has been previously shown to rapidly predict bimetallic NP cohesive energies (CEs). Specifically, the BC model is trained in a way to understand how the bimetallic bond strength changes under different coordination environments present on a NP and how the metal composition of every site affects the detailed coordination environment using fractional coordination numbers. This newly modified BC model leads to an improvement from 0.331 (original model) to 0.089 eV/atom in CE predictions when compared to DFT values on a robust data set of 90 different NPs consisting of PtPd, AuPt, and AuPd NPs with varying compositions and chemical orderings. By incorporating the modified BC model into an in-house-developed genetic algorithm (GA) we can effectively and accurately predict the most stable chemical orderings of large, realistic bimetallic NPs consisting of thousands of metal atoms. This is demonstrated on AuPd bimetallic NPs, a challenging system due to the similarity in the cohesion of the two metals. By training our BC model using a unique DFT calculation on a bimetallic NP (one calculation for two metals combining together), we expand to explore the chemical ordering of multimetallic NPs. We first demonstrate the application of our methodology on a AuPdPt NP and validate our stability predictions with literature data. Then, we effectively explore the vast materials space of multimetallic NPs consisting of combinations of Au, Pt, and Pd as a function of metal composition. Our thermodynamic stability trends are presented in a ternary diagram revealing detailed, and yet, unexpected chemical ordering trends. Our computational framework can aid both experimental and computational researchers toward effectively screening multimetallic NP stability. Moreover, we provide an outlook of how this framework can be applied to catalyst discovery, high-entropy alloys, and single-atom alloys.

Heteroatom Manipulation of Zeolite Crystallization: Stabilizing Zn-FAU against Interzeolite Transformation

The preparation of metastable zeolites is often restricted to a limited range of synthesis conditions, which is exemplified in commercial syntheses lacking organics to stabilize the crystal structure. In the absence of an organic structure-directing agent, interzeolite transformation is a common phenomenon that can lead to undesirable products or impurities. Many studies have investigated the substitution of Si and Al in zeolite frameworks with alternative elements (heteroatoms) as a means of tailoring the properties of zeolites; however, relatively few studies have systematically explored the impact of heteroatoms on interzeolite transformations and their concomitant effects on zeolite crystallization. In this study, we examine methods to prepare isostructures of faujasite (FAU), which is one of the most commercially relevant zeolites and also a thermodynamically metastable structure. A survey of multivalent elements revealed that zinc is capable of stabilizing FAU at high temperatures and inhibiting its frequent transformation to zeolite gismondine (GIS). Using combined experimental and computational studies, we show that zinc alters the chemical nature of growth mixtures by sequestering silicates. Zinc heteroatoms incorporate in the FAU framework with a loading-dependent coordination. Our collective findings provide an improved understanding of driving forces for the FAU-to-GIS interzeolite transformation where we observe that heteroatoms (e.g., zinc) can stabilize zeolite FAU over a broad range of synthesis conditions. Given the growing interest in heteroatom-substituted zeolites, this approach to preparing zinc-containing FAU may prove applicable to a broader range of zeolite structures.

Size-dependent shape distributions of platinum nanoparticles

While it is well established that nanoparticle shape can depend on equilibrium thermodynamics or growth kinetics, recent computational work has suggested the importance of thermal energy in controlling the distribution of shapes in populations of nanoparticles. Here, we used transmission electron microscopy to characterize the shapes of bare platinum nanoparticles and observed a strong dependence of shape distribution on particle size. Specifically, the smallest nanoparticles (<2.5 nm) had a truncated octahedral shape, bound by 〈111〉 and 〈100〉 facets, as predicted by lowest-energy thermodynamics. However, as particle size increased, the higher-energy 〈110〉 facets became increasingly common, leading to a large population of non-equilibrium truncated cuboctahedra. The observed trends were explained by combining atomistic simulations (both molecular dynamics and an empirical square-root bond-cutting model) with Boltzmann statistics. Overall, this study demonstrates experimentally how thermal energy leads to shape variation in populations of metal nanoparticles, and reveals the dependence of shape distributions on particle size. The prevalence of non-equilibrium facets has implications for metal nanoparticles applications from catalysis to solar energy.

The Sensitivity of Metal Oxide Electrocatalysis to Bulk Hydrogen Intercalation: Hydrogen Evolution on Tungsten Oxide

Metal oxides are attracting increased attention as electrocatalysts owing to their affordability, tunability, and reactivity. However, these materials can undergo significant chemical changes under reaction conditions, presenting challenges for characterization and optimization. Herein, we combine experimental and computational methods to demonstrate that bulk hydrogen intercalation governs the activity of tungsten trioxide (WO₃) toward the hydrogen evolution reaction (HER). In contrast to the focus on surface processes in heterogeneous catalysis, we demonstrate that bulk oxide modification is responsible for experimental HER activity. Density functional theory (DFT) calculations reveal that intercalation enables the HER by altering the acid-base character of surface sites and preventing site blocking by hydration. First-principles microkinetic modeling supports that the experimental HER rates can only be explained by intercalated H_xWO₃, whereas nonintercalated WO₃ does not catalyze the HER. Overall, this work underscores the critical influence of hydrogen intercalation on aqueous cathodic electrocatalysis at metal oxides.

Predicting Segregation Energy in Single Atom Alloys Using Physics and Machine Learning

Single atom alloys (SAAs) show great promise as catalysts for a wide variety of reactions due to their tunable properties, which can enhance the catalytic activity and selectivity. To design SAAs, it is imperative for the heterometal dopant to be stable on the surface as an active catalytic site. One main approach to probe SAA stability is to calculate surface segregation energy. Density functional theory (DFT) can be applied to investigate the surface segregation energy in SAAs. However, DFT is computationally expensive and time-consuming; hence, there is a need for accelerated frameworks to screen metal segregation for new SAA catalysts across combinations of metal hosts and dopants. To this end, we developed a model that predicts surface segregation energy using machine learning for a series of SAA periodic slabs. The model leverages elemental descriptors and features inspired by the previously developed bond-centric model. The initial model accurately captures surface segregation energy across a diverse series of FCC-based SAAs with various surface facets and metal-host pairs. Following our machine learning methodology, we expanded our analysis to develop a new model for SAAs formed from FCC hosts with FCC, BCC, and HCP dopants. Our final, five-feature model utilizes second-order polynomial kernel ridge regression. The model is able to predict segregation energies with a high degree of accuracy, which is due to its physically motivated features. We then expanded our data set to test the accuracy of the five features used. We find that the retrained model can accurately capture E _seg trends across different metal hosts and facets, confirming the significance of the features used in our final model. Finally, we apply our pretrained model to a series of Ir- and Pd-based SAA cuboctahedron nanoparticles (NPs), ranging in size and FCC dopants. Remarkably, our model (trained on periodic slabs) accurately predicts the DFT segregation energies of the SAA NPs. The results provide further evidence supporting the use of our model as a general tool for the rapid prediction of SAA segregation energies. By creating a framework to predict the metal segregation from bulk surfaces to NPs, we can accelerate the SAA catalyst design while simultaneously unraveling key physicochemical properties driving thermodynamic stabilization of SAAs.

A Classification Model to Identify Direct-Acting Mutagenic Polycyclic Aromatic Hydrocarbon Transformation Products

Polycyclic aromatic hydrocarbons (PAHs) are a complex group of environmental contaminants, many having long environmental half-lives. As these compounds degrade, the changes in their structure can result in a substantial increase in mutagenicity compared to the parent compound. Over time, each individual PAH can potentially degrade into several thousand unique transformation products, creating a complex, constantly evolving set of intermediates. Microbial degradation is the primary mechanism of their transformation and ultimate removal from the environment, and this process can result in mutagenic activation similar to the metabolic activation that can occur in multicellular organisms. The diversity of the potential intermediate structures in PAH-contaminated environments renders hazard assessment difficult for both remediation professionals and regulators. A mixture of structural and energetic descriptors has proven effective in existing studies for classifying which PAH transformation products will be mutagenic. However, most existing studies of environmental PAH mutagens primarily focus on nitrogenated derivatives, which are prevalent in the atmosphere and not as relevant in soil. Additionally, PAH products commonly found in the environment can range from as large as five rings to as small as a single ring, requiring a broadly inclusive methodology to comprehensively evaluate mutagenic potential. We developed a combination of supervised and unsupervised machine learning methods to predict environmentally induced PAH mutagenicity with improved performance over currently available tools. K-means clustering with principal component analysis allows us to identify molecular clusters that we hypothesize to have similar mechanisms of action. Recursive feature elimination identifies the most influential descriptors. The cluster-specific regression outperforms available classifiers in predicting direct-acting mutagens resulting from the microbial biodegradation of PAHs and provides direction for future studies evaluating the environmental hazards resulting from PAH biodegradation.

CuZrO3: If it exists it should be a sandwich

CuZrO₃ has been hypothesized to be a catalytic material with potential applications for CO₂ reduction. Unfortunately, this material has received limited attention in the literature, and to the best of our knowledge the exact crystal structure is still unknown. To address this challenge, we utilize several different structural prediction techniques in concert, including the Universal Structure Predictor: Evolutionary Xtallography (USPEX), the Materials Project Structure Predictor, and the Open Quantum Materials Database (OQMD). Leveraging these structural prediction techniques in conjunction with Density-Functional Theory (DFT) calculations, we determine a possible structure for CuZrO₃, which resembles a "sandwich" morphology. Our calculations reveal that this new structure is significantly lower in energy than a previously hypothesized perovskite structure, albeit it still has a thermodynamic preference to decompose into CuO and ZrO₂. In addition, we experimentally tried to synthesize CuZrO₃ based on literature reports and compared computational to experimental X-ray Diffraction (XRD) patterns confirming that the final product is a mixture of CuO and ZrO₂. Finally, we conducted a computational surface energetics and CO₂ adsorption study on our discovered sandwich morphology, demonstrating that CO₂ can adsorb and activate on the material. However, these CO₂ adsorption results deviate from previously reported results further confirming that the CuZrO₃ is a metastable form and may not be experimentally accessible as a well-mixed oxide, since phase segregation to CuO and ZrO₂ is preferred. Taken together, our combined computational and experimental study provides evidence that the synthesis of CuZrO₃ is extremely difficult and if this oxide exists, it should have a sandwich-like morphology.

Hydrogen Evolution Electrocatalyst Design: Turning Inert Gold into Active Catalyst by Atomically Precise Nanochemistry

Electrocatalytic hydrogen evolution reaction (HER) holds promise in the renewable clean energy scheme. Crystalline Au and Ag are, however, poor in catalyzing HER, and the ligands on colloidal nanoparticles are generally another disadvantage. Herein, we report a thiolate (SR)-protected Au₃₆Ag₂(SR)₁₈ nanocluster with low coverage of ligands and a core composed of three icosahedral (I_h) units for catalyzing HER efficiently. This trimeric structure, together with the monomeric I_h Au₂₅(SR)₁₈^- and dimeric I_h Au₃₈(SR)₂₄, constitutes a unique series, providing an opportunity for revealing the correlation between the catalytic properties and the catalyst's structure. The Au₃₆Ag₂(SR)₁₈ surprisingly exhibits high catalytic activity at lower overpotentials for HER due to its low ligand-to-metal ratio, low-coordinated Au atoms and unfilled superatomic orbitals. The current density of Au₃₆Ag₂(SR)₁₈ at -0.3 V vs RHE is 3.8 and 5.1 times that of Au₂₅(SR)₁₈^- and Au₃₈(SR)₂₄, respectively. Density functional theory (DFT) calculations reveal lower hydrogen binding energy and higher electron affinity of Au₃₆Ag₂(SR)₁₈ for an energetically feasible HER pathway. Our findings provide a new strategy for constructing highly active catalysts from inert metals by pursuing atomically precise nanoclusters and controlling their geometrical and electronic structures.

Correlating structural rules with electronic properties of ligand-protected alloy nanoclusters

Thiolate protected gold nanoclusters (TPNCs) are a unique class of nanomaterials finding applications in various fields, such as biomedicine, optics, and catalysis. The atomic precision of their structure, characterized through single crystal x-ray diffraction, enables the accurate investigation of their physicochemical properties through electronic structure calculations. Recent experimental efforts have led to the successful heterometal doping of TPNCs, potentially unlocking a large domain of bimetallic TPNCs for targeted applications. However, how TPNC size, bimetallic composition, and location of dopants influence electronic structure is unknown. To this end, we introduce novel structure-property relationships (SPRs) that predict electronic properties such as ionization potential (IP) and electron affinity (EA) of AgAu TPNCs based on physically relevant descriptors. The models are constructed by first generating a hypothetical AgAu TPNC dataset of 368 structures with sizes varying from 36 to 279 metal atoms. Using our dataset calculated with density functional theory (DFT), we employed systematic analyses to unravel size, composition, and, importantly, core-shell effects on TPNC EA and IP behavior. We develop generalized SPRs that are able to predict electronic properties across the AgAu TPNC materials space. The models leverage the same three fundamental descriptors (i.e., size, composition, and core-shell makeup) that do not require DFT calculations and rely only on simple atom counting, opening avenues for high throughput bimetallic TPNC screening for targeted applications. This work is a first step toward finely controlling TPNC electronic properties through heterometal doping using high throughput computational means.

The role of ligands in atomically precise nanocluster-catalyzed CO2 electrochemical reduction

Ligand effects are of major interest in catalytic reactions owing to their potential critical role in determining the reaction activity and selectivity. Herein, we report ligand effects in the CO2 electrochemical reduction reaction at the atomic level with three unique Au25 nanoclusters comprising the same kernel but different protecting ligands (-XR, where X = S or Se, and R represents the carbon tail). It is observed that a change in the carbon tail shows no obvious impact on the catalytic selectivity and activity, but the anchoring atom (X = S or Se) strongly affects the electrocatalytic selectivity. Specifically, the S site acts as the active site and sustains CO selectivity, while the Se site shows a higher tendency of hydrogen evolution. Density functional theory (DFT) calculations reveal that the energy penalty associated with the *COOH formation is lower on the S site by 0.26 eV compared to that on the Se site. Additionally, the formation energy of the product (*CO) is lower on the sulfur-based Au nanocluster by 0.43 eV. We attribute these energetic differences to the higher electron density on the sulfur sites of the Au nanocluster, resulting in a modified bonding character of the reaction intermediates that reduce the energetic penalty for the *COOH and *CO formation. Overall, this work demonstrates that S/Se atoms at the metal-ligand interface can play an important role in determining the overall electrocatalytic performance of Au nanoclusters.

Predicting ligand removal energetics in thiolate-protected nanoclusters from molecular complexes

Thiolate-protected metal nanoclusters (TPNCs) have attracted great interest in the last few decades due to their high stability, atomically precise structure, and compelling physicochemical properties. Among their various applications, TPNCs exhibit excellent catalytic activity for numerous reactions; however, recent work revealed that these systems must undergo partial ligand removal in order to generate active sites. Despite the importance of ligand removal in both catalysis and stability of TPNCs, the role of ligands and metal type in the process is not well understood. Herein, we utilize Density Functional Theory to understand the energetic interplay between metal-sulfur and sulfur-ligand bond dissociation in metal-thiolate systems. We first probe 66 metal-thiolate molecular complexes across combinations of M = Ag, Au, and Cu with twenty-two different ligands (R). Our results reveal that the energetics to break the metal-sulfur and sulfur-ligand bonds are strongly correlated and can be connected across all complexes through metal atomic ionization potentials. We then extend our work to the experimentally relevant [M25(SR)18]- TPNC, revealing the same correlations at the nanocluster level. Importantly, we unify our work by introducing a simple methodology to predict TPNC ligand removal energetics solely from calculations performed on metal-ligand molecular complexes. Finally, a computational mechanistic study was performed to investigate the hydrogenation pathways for SCH3-based complexes. The energy barriers for these systems revealed, in addition to thermodynamics, that kinetics favor the break of S-R over the M-S bond in the case of the Au complex. Our computational results rationalize several experimental observations pertinent to ligand effects on TPNCs. Overall, our introduced model provides an accelerated path to predict TPNC ligand removal energies, thus aiding towards targeted design of TPNC catalysts.

Solvent manipulation of the pre-reduction metal-ligand complex and particle-ligand binding for controlled synthesis of Pd nanoparticles

Understanding how to control the nucleation and growth rates is crucial for designing nanoparticles with specific sizes and shapes. In this study, we show that the nucleation and growth rates are correlated with the thermodynamics of metal-ligand/solvent binding for the pre-reduction complex and the surface of the nanoparticle, respectively. To obtain these correlations, we measured the nucleation and growth rates by in situ small angle X-ray scattering during the synthesis of colloidal Pd nanoparticles in the presence of trioctylphosphine in solvents of varying coordinating ability. The results show that the nucleation rate decreased, while the growth rate increased in the following order, toluene, piperidine, 3,4-lutidine and pyridine, leading to a large increase in the final nanoparticle size (from 1.4 nm in toluene to 5.0 nm in pyridine). Using density functional theory (DFT), complemented by 31P nuclear magnetic resonance and X-ray absorption spectroscopy, we calculated the reduction Gibbs free energies of the solvent-dependent dominant pre-reduction complex and the solvent-nanoparticle binding energy. The results indicate that lower nucleation rates originate from solvent coordination which stabilizes the pre-reduction complex and increases its reduction free energy. At the same time, DFT calculations suggest that the solvent coordination affects the effective capping of the surface where stronger binding solvents slow the nanoparticle growth by lowering the number of active sites (not already bound by trioctylphosphine). The findings represent a promising advancement towards understanding the microscopic connection between the metal-ligand thermodynamic interactions and the kinetics of nucleation and growth to control the size of colloidal metal nanoparticles.

Mechanistic understanding of methane-to-methanol conversion on graphene-stabilized single-atom iron centers

DFT calculations and kinetic modeling elucidate solvent effects and complex mechanisms for the room-temperature methane-to-methanol conversion on an FeN4/graphene catalyst.

Atom-by-Atom Evolution of the Same Ligand-Protected Au21, Au22, Au22Cd1, and Au24 Nanocluster Series

Atom-by-atom manipulation on metal nanoclusters (NCs) has long been desired, as the resulting series of NCs can provide insightful understanding of how a single atom affects the structure and properties as well as the evolution with size. Here, we report crystallizations of Au22(SAdm)16 and Au22Cd1(SAdm)16 (SAdm = adamantanethiolate) which link up with Au21(SAdm)15 and Au24(SAdm)16 NCs and form an atom-by-atom evolving series protected by the same ligand. Structurally, Au22(SAdm)16 has an Au3(SAdm)4 surface motif which is longer than the Au2(SAdm)3 on Au21(SAdm)15, whereas Au22Cd1(SAdm)16 lacks one staple Au atom compared to Au24(SAdm)16 and thus the surface structure is reconstructed. A single Cd atom triggers the structural transition from Au22 with a 10-atom bioctahedral kernel to Au22Cd1 with a 13-atom cuboctahedral kernel, and correspondingly, the optical properties are dramatically changed. The photoexcited carrier lifetime demonstrates that the optical properties and excited state relaxation are highly sensitive at the single atom level. By contrast, little change in both ionization potential and electron affinity is found in this series of NCs by theoretical calculations, indicating the electronic properties are independent of adding a single atom in this series. The work provides a paradigm that the NCs with continuous metal atom numbers are accessible and crystallizable when meticulously designed, and the optical properties are more affected at the single atom level than the electronic properties.

Predicting the Energetics of Hydrogen Intercalation in Metal Oxides Using Acid-Base Properties

The ability to predict intercalation energetics from first principles is attractive for identifying candidate materials for energy storage, chemical sensing, and catalysis. In this work, we introduce a computational framework that can be used to predict the thermodynamics of hydrogen intercalation in tungsten trioxide (WO₃). Specifically, using density functional theory (DFT), we investigated intercalation energetics as a function of adsorption site and hydrogen stoichiometry. Site-specific acid-base properties determined using DFT were used to develop linear structure screening models that informed a kernel ridge energy prediction model. These regressions provided a series of hydrogen binding energy predictions across stoichiometries ranging from WO₃ to H_0.625WO₃, which were then converted to equilibrium potentials for hydrogen intercalation. Experimental validation using cyclic voltammetry measurements yielded good agreement with the predicted intercalation potentials. This methodology enables fast exploration of a large geometric configuration space and reveals an intuitive physical relationship between acidity, basicity, and the thermodynamics of hydrogen intercalation. Furthermore, the combination of theoretical and experimental results suggests H_0.500WO₃ as a maximum stable stoichiometry for the bronzes that arises from competition with hydrogen evolution rather than the inability of WO₃ to accommodate additional hydrogen. Our experimental results further indicate hydrogen insertion in WO₃ is highly irreversible for low H-stoichiometries, which we propose to be a consequence of the semiconductor-to-metal transition that occurs upon initial H-intercalation. Overall, the agreement between theory and experiment suggests that local acid-base characteristics govern hydrogen intercalation in tungsten trioxide, and this insight can aid the accelerated discovery of redox-active metal oxides for catalytic hydrogenations.

Heteroatom Tracing Reveals the 30-Atom Au-Ag Bimetallic Nanocluster as a Dimeric Structure

Understanding the formation of face-centered cubic (fcc) nanostructures at the atomic level remains a major task. With atomically precise nanoclusters (NCs) as model systems, herein we devised an atom-tracing strategy by heteroatom doping into Au₃₀(SR)₁₈ (SR = S-^tC₄H₉) to label the specific positions in M₃₀(SR)₁₈ NCs (M = Au/Ag), which clearly reveals the dimeric nature of M₃₀. Interestingly, the specific position is also consistent with the Ag-doping site in M₂₁(SR)₁₅. Electronic orbital analysis shows intrinsic orbital localization at the two specific positions in M₃₀, which are decisive to the electronic structure of M₃₀, regardless of Au or Ag occupancy. The fcc dimeric NC, which would not be discovered without Ag tracing, provides a possible explanation for the wide accessibility of nonsuperatomic Au-SR NCs.

Rapid Prediction of Bimetallic Mixing Behavior at the Nanoscale

The nanoparticle (NP) design space allows for variations in size, shape, composition, and chemical ordering. In the search for low-energy structures, this results in an extremely large search space which cannot be screened by brute force methods. In this work, we develop a genetic algorithm to predict stable bimetallic NPs of any size, shape, and metal composition. Our method predicts nanostructures in agreement with experimental trends and it captures the detailed chemical ordering of an experimental 23,196-atom FePt NP with nearly atom-by-atom accuracy. Our developed screening process is extremely fast, allowing us to generate and analyze a database of 5454 low-energy bimetallic NPs. By identifying thermodynamically stable NPs, we rationalize bimetallic mixing at the nanoscale and reveal metal-, size-, and temperature-dependent mixing behavior. Importantly, our method is applicable to any bimetallic NP size, bridging the materials gap in nanoscale simulations, and guides experimentation in the lab by elucidating stability, mixing, and detailed chemical ordering behavior of bimetallic NPs.

Heterometal-Doped M23 (M = Au/Ag/Cd) Nanoclusters with Large Dipole Moments

Dipole moment (μ) is a critical parameter for molecules and nanomaterials as it affects many properties. In metal-thiolate (SR) nanoclusters (NCs), μ is commonly low (0-5 D) compared to quantum dots. Herein, we report a doping strategy to give giant dipoles (∼18 D) in M₂₃ (M = Au/Ag/Cd) NCs, falling in the experimental trend for II-VI quantum dots. In M₂₃ NCs, high μ is caused by the Cd-Br bond and the arrangement of heteroatoms along the C₃ axis. Strong dipole-dipole interactions are observed in crystalline state, with energy exceeding 5 kJ/mol, directing a "head-to-tail" alignment of Au_22-nAg_nCd₁(SR)₁₅X (SR = adamantanethiolate) dipoles. The alignment can be controlled by μ via doping. The optical absorption peaks of M₂₃ show solvent polarity-dependent shifts (∼25 meV) with negative solvatochromism. Detailed electronic structures of M₂₃ are revealed by density functional theory and time-dependent DFT calculations. Overall, the doping strategy for obtaining large dipole moments demonstrates an atomic-level design of clusters with useful properties.

Ligand exchange on Au38(SR)24: substituent site effects of aromatic thiols

Understanding the critical roles of ligands (e.g. thiolates, SR) in the formation of metal nanoclusters of specific sizes has long been an intriguing task since the report of ligand exchange-induced transformation of Au38(SR)24 into Au36(SR')24. Herein, we conduct a systematic study of ligand exchange on Au38(SC2H4Ph)24 with 21 incoming thiols and reveal that the size/structure preference is dependent on the substituent site. Specifically, ortho-substituted benzenethiols preserve the structure of Au38(SR)24, while para- or non-substituted benzenethiols cause its transformation into Au36(SR)24. Strong electron-donating or -withdrawing groups do not make a difference, but they will inhibit full ligand exchange. Moreover, the crystal structure of Au38(SR)24 (SR = 2,4-dimethylbenzenethiolate) exhibits distinctive ππ stacking and "anagostic" interactions (indicated by substantially short AuH distances). Theoretical calculations reveal the increased energies of frontier orbitals for aromatic ligand-protected Au38, indicating decreased electronic stability. However, this adverse effect could be compensated for by the AuH-C interactions, which improve the geometric stability when ortho-substituted benzenethiols are used. Overall, this work reveals the substituent site effects based on the Au38 model, and highlights the long-neglected "anagostic" interactions on the surface of Au-SR NCs which improve the structural stability.

Improving alkane dehydrogenation activity on γ-Al2O3 through Ga doping

Doping the surface of γ-Al₂O₃ with gallium enhances the alkane dehydrogenation catalytic activity.

Towards elucidating structure of ligand-protected nanoclusters

Developing a centralized database for ligand-protected nanoclusters can fuel machine learning and data-science-based approaches towards theoretical structure prediction.

The role of nanoparticle size and ligand coverage in size focusing of colloidal metal nanoparticles

Controlling the size distribution of nanoparticles is important for many applications and typically involves the use of ligands during synthesis. In this study, we show that the mechanism of size focusing involves a dependence of the growth rate on the size of the nanoparticles and the ligand coverage on the surface of the nanoparticles. To demonstrate these effects, we used in situ small angle X-ray scattering (SAXS) and population balance kinetic modeling (PBM) to investigate the evolution of size distribution during the synthesis of colloidal Pd metal nanoparticles. Despite temporal overlap of nucleation and growth, our in situ SAXS show size focusing of the distribution under different synthetic conditions (different concentrations of metal and ligand as well as solvent type). To understand the mechanism of size focusing using PBM, we systematically studied how the evolution of the nanoparticle size distribution is affected by nucleation rate, and dependence of the growth rate constant on ligand surface coverage, and size of the nanoparticles. We show that continuous nucleation contributes to size defocusing. However, continuous nucleation results in different reaction times for the nanoparticle population leading to time and size-dependent ligand surface coverage. Using density functional theory (DFT) calculations and Brønsted-Evans-Polanyi relations, we show that as the population grows, larger nanoparticles grow more slowly than smaller ones due to lower intrinsic activity and higher ligand coverage on the surface. Therefore, despite continuous nucleation, the faster growth of smaller nanoparticles in the population leads to size focusing. The size focusing behaviour (due to faster growth of smaller nanoparticles) was found to be model independent and similar results were demonstrated under different nucleation and growth pathways (e.g. growth via ion reduction on the surface and/or monomer addition). Our results provide a microscopic connection between kinetics and thermodynamics of nanoparticle growth and metal-ligand binding, and their effect on the size distribution of colloidal nanoparticles.

Unfolding adsorption on metal nanoparticles: Connecting stability with catalysis

Metal nanoparticles have received substantial attention in the past decades for their applications in numerous areas, including medicine, catalysis, energy, and the environment. Despite these applications, the fundamentals of adsorption on nanoparticle surfaces as a function of nanoparticle size, shape, metal composition, and type of adsorbate are yet to be found. Herein, we introduce the first universal adsorption model that accounts for detailed nanoparticle structural characteristics, metal composition, and different adsorbates by combining first principles calculations with machine learning. Our model fits a large number of data and accurately predicts adsorption trends on nanoparticles (both monometallic and alloy) that have not been previously seen. In addition to its application power, the model is simple and uses tabulated and rapidly calculated data for metals and adsorbates. We connect adsorption with stability behavior of nanoparticles, thus advancing the design of optimal nanoparticles for applications of interest, such as catalysis.

Structure-property relationships on thiolate-protected gold nanoclusters

Since their discovery, thiolate-protected gold nanoclusters (Au _n (SR) _m ) have garnered a lot of interest due to their fascinating properties and "magic-number" stability. However, models describing the thermodynamic stability and electronic properties of these nanostructures as a function of their size are missing in the literature. Herein, we employ first principles calculations to rationalize the stability of fifteen experimentally determined gold nanoclusters in conjunction with a recently developed thermodynamic stability theory on small Au nanoclusters (≤102 Au atoms). Our results demonstrate that the thermodynamic stability theory can capture the stability of large, atomically precise nanoclusters, Au₂₇₉(SR)₈₄, Au₂₄₆(SR)₈₀, and Au₁₄₆(SR)₅₇, suggesting its applicability over larger cluster size regimes than its original development. Importantly, we develop structure-property relationships on Au nanoclusters, connecting their ionization potential and electron affinity to the number of gold atoms within the nanocluster. Altogether, a computational scheme is described that can aid experimental efforts towards a property-specific, targeted synthesis of gold nanoclusters.

Elucidating the optical spectra of [Au25(SR)18]q nanoclusters

Density functional theory calculations reveal how optical spectra of [Au₂₅(SR)₁₈]^q nanoclusters (q = −1, 0, +1) change with different ligands.

Computational Study of Methane Activation on γ-Al2O3

The C-H activation of methane remains a longstanding challenge in the chemical industry. Metal oxides are attractive catalysts for the C-H activation of methane due to their surface Lewis acid-base properties. In this work, we applied density functional theory calculations to investigate the C-H activation mechanism of methane on various sites of low-index facets of γ-Al₂O₃. The feasibility of C-H activation on different metal-oxygen (acid-base) site pairs was assessed through two potential mechanisms, namely, the radical and polar. The effect of surface hydroxylation on C-H activation was also investigated to examine the activity of γ-Al₂O₃ under realistic catalytic surface conditions (hydration). On the basis of our calculations, it was demonstrated that the C-H activation barriers for polar pathways are significantly lower than those of the radical pathways on γ-Al₂O₃. We showed that the electronic structure (s- and p-band center) for unoccupied and occupied bands can be used to probe site-dependent Lewis acidity and basicity and the associated catalytic behavior. We identified the dissociated H₂ binding and final state energy as C-H activation energy descriptors for the preferred polar pathway. Finally, we developed structure-activity relationships for the C-H activation of methane on γ-Al₂O₃ that account for surface Lewis acid-base properties and can be utilized to accelerate the discovery of catalysts for methane (and shale gas) upgrade.

Rethinking Heterometal Doping in Ligand-Protected Metal Nanoclusters

Heterometal doping is a promising avenue toward tailoring properties of ligand-protected metal nanoclusters for specific applications. Though successful doping has been demonstrated in several structures, the underlying reasons for the dopant preference on occupying specific locations on the nanocluster with different concentrations remain unclear. In this study we apply our thermodynamic stability model, originally developed for ligand-protected monometallic nanoclusters, to rationalize the synthetic accessibility, dopant location, and concentrations of various heterometals on ligand-protected Au nanoclusters. Importantly, we demonstrate that the thermodynamic stability theory is a significant step forward in accurately describing doping effects on nanoclusters using first-principles calculations. With our computational predictions being in excellent agreement with a series of experiments, we introduce the thermodynamic stability theory as a new method for bimetallic nanocluster prediction.