The interaction of oxygen with small gold clusters

Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal-organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule-MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.

Various Bond Interactions between NO and Anionic Gold Clusters: A Theoretical Calculation

We studied the electronic and geometrical structures of AunNO- (n = 1-20) using the B3LYP method with relatively large basis sets to understand the size dependent reactivities of Aun- with...

Exploring the Effects of a Doping Silver Atom on Anionic Gold Clusters' Reactivity with O2

Reactivities of AgAu_n-1^- (n = 3-10) with O₂ at a low temperature were studied using an instrument combining a magnetron sputter cluster source, a microflow reactor, and a time-of-flight mass spectrometer. Their reaction products as well as size-dependent kinetic rates were nearly identical to those of corresponding Au_n^- (n = 3-10). Previous experiments showed that the Ag atom in AgAu_n-1^- (n = 3-10) was fully or partially enclosed by the gold atoms. We studied the adsorption of O₂ on these reported structures using the B3LYP theory with relatively large basis sets. The theoretical results indicate that the adsorption sites as well as the adsorption energies of O₂ on AgAu_n-1^- (n = 3-10) are nearly identical to those on the corresponding Au_n^- (n = 3-10). The O₂ adsorption on a series of proposed isomers of AgAu_n-1^- (denoted as Au_n-1Ag^-), in which the silver atom was on the protruding site, was explored using the same theoretical methods. The O₂ tends to bond with the protruding Ag atoms, and the binding energies are apparently higher than those on the corresponding Au_n^- and AgAu_n-1^-. The adsorption and activation of O₂ on Au_n^-, AgAu_n-1^-, and Au_n-1Ag^- were correlated with their global electron detachment energies (VDEs) as well as the element types of the adsorption sites. Generally, low VDE values and silver sites facilitate the O₂ adsorption, and these two factors separately dominate in various cluster species. The revealed effects of a doping silver atom in small gold clusters are helpful to understand the role of the residual silver components in many nano gold catalysts.

Catalytic Oxidation of Benzyl Alcohol to Benzaldehyde on Au8 and Au6Pd2 Clusters: A DFT Study on the Reaction Mechanism

Density functional theory calculations were performed to investigate the reaction mechanism of the aerobic oxidation of benzyl alcohol to benzaldehyde catalyzed by Au and Au–Pd clusters. Two consecutive reaction mechanisms were examined with Au8 and Au6Pd2 clusters: (1) the oxidation of benzyl alcohol with dissociated O atoms on metal clusters generating benzaldehyde and H2O; and (2) oxidation with adsorbed oxygen molecules generating benzaldehyde and H2O2. The calculations show that the aerobic oxidation of benzyl alcohol energetically prefers to proceed in the former mechanism, which agrees with the experimental observation. We demonstrate that the role of Au centers around the activation of molecular oxygen to peroxide-like species, which are capable of the H–abstraction of benzyl alcohol. The roles of Pd in the Au6Pd2 cluster are: (1) increasing the electron distribution to neighboring Au atoms, which facilitates the activation of O2; and (2) stabilizing the adsorption complex and transition states by the interaction between positively charged Pd atoms and the π-bond of benzyl alcohol, both of which are the origin of the lower energy barriers than those of Au8.

Polarizability of atomic Pt, Pt+, and Pt. J Chem Phys 2021;154:174302. [PMID: 34241047 DOI: 10.1063/5.0044996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Electrostatic properties are important for understanding and modeling many phenomena, such as the adsorption of a catalytic metal upon an oxide support. The charge transfer between the metal and the support can lead to positive or negative charges on the metal. Here, the static dipole polarizability is computed for atomic platinum in charge states 0, +1, and -1 in several low-lying electronic terms and levels. Core pseudopotentials are used along with coupled-cluster theory. The best results are estimates for the coupled-cluster CCSDTQ/q-aug-cc-pwCV∞Z-PP values for atomic terms, combined with compositional data from spin-orbit configuration interaction. The polarizability of the anion Pt^- is especially challenging for the theory with wildly varying results from different coupled-cluster perturbative approximations such as CCSD(T). For atomic mercury (Hg), selected as a nearby experimental value, our polarizability volume is larger than experiment by 0.8 bohrs³ (or 0.12 × 10^-30 m³). For the ground level of neutral platinum, Pt(³D₃), we find α₀ = (41.2 ± 1.1) bohrs³ or (6.10 ± 0.16) × 10^-30 m³. A handful of density functional theory methods are tested and found generally within 10% of our best values.

How O2-Binding Affects Structural Evolution of Medium Even-Sized Gold Clusters Aun- (n = 20-34)

We report the first joint anion photoelectron spectroscopy and theoretical study on how O₂-binding affects the structures of medium even-sized gold clusters, Au_n^- (n = 20-34), a special size region that entails a variety of distinct structures. Under the temperature conditions in the current photoelectron spectroscopy experiment, O₂-bound gold clusters were observed only for n = 22-24 and 34. Nevertheless, O₂ binding with the clusters in the size range of n = 20-34 can be still predicted based on the obtained global-minimum structures. Consequently, a series of structural transitions, from the pyramidal to fused-planar to core-shell structures, are either identified or predicted for the Au_nO₂^- clusters, where the O₂-binding is in either superoxo or peroxo fashion. The identified global-minimum structures of Au_nO₂^- (n = 20-34) also allow us to gain improved understanding of why the clusters Au_n^- (n = 26-32) are less reactive with O₂ in comparison to others.

Probing the structural, electronic, and adsorptive properties of Au16O2- clusters

Great progress has been made in O₂ adsorption on gold clusters. However, systematic investigations of O₂ adsorption on [Formula: see text] clusters have not been reported. Here, we present a systematic study of the structural, electronic, and adsorptive properties of [Formula: see text] clusters by density functional theory (DFT) calculations coupled with stochastic kicking method. Global minimum searches for [Formula: see text] reveal that exohedral derivatives are more favored. Furthermore, the obtained ground-state structure exhibits significant stability, as judged by its larger adsorption energy (1.16 eV) and a larger HOMO-LUMO gap (0.57 eV). The simulated photoelectron spectra (PES) of [Formula: see text] isomers will be instructive to identify the structures in future experiments. There are three interesting discoveries in the present paper: (1) O₂ undergoes chemical adsorption onto the parent [Formula: see text] clusters, but the amount of the adsorption energy is related to the parent [Formula: see text] clusters; (2) the process that O₂ undergoes dissociative adsorption onto the parent [Formula: see text] clusters is exothermic; (3) [Formula: see text] isomers show smaller X-A energy gaps than those of parent [Formula: see text] clusters, reflecting that their geometric and electronic structures are distorted remarkably due to dissociative adsorption of O₂.

Co-Adsorption of O2 and CO on Au2- : Infrared Photodissociation Spectroscopy and Theoretical Study of [Au2 O2 (CO)n ]- (n=2-6)

The co-adsorption of O₂ and CO on anionic sites of gold species is considered as a crucial step in the catalytic CO oxidation on gold catalysts. In this regard, the [Au₂ O₂ (CO)_n ]^- (n=2-6) complexes were prepared by using a laser vaporization supersonic ion source and were studied by using infrared photodissociation spectroscopy in the gas phase. All the [Au₂ O₂ (CO)_n ]^- (n=2-6) complexes were characterized to have a core structure involving one CO and one O₂ molecule co-adsorbed on Au₂^- with the other CO molecules physically tagged around. The CO stretching frequency of the [Au₂ O₂ (CO)]^- core ion is observed around ν˜ =2032-2042 cm^-1 , which is about 200 cm^-1 higher than that in [Au₂ (CO)₂ ]^- . This frequency difference and the analyses based on density functional calculations provide direct evidence for the synergy effect of the chemically adsorbed O₂ and CO. The low lying structures with carbonate group were not observed experimentally because of high formation barriers. The structures and the stability (i.e., the inertness in a sense) of the co-adsorbed O₂ and CO on Au₂^- may have relevance to the elementary reaction steps on real gold catalysts.

Nucleobases tagged to gold nanoclusters cause a mechanistic crossover in the oxidation of CO

A mechanistic crossover is observed upon using nucleobase tagged gold clusters as catalysts favoring the Eley–Rideal mechanism, over the conventional Langmuir–Hinshelwood pathway followed using pristine gold clusters during CO oxidation.

Aerobic oxidation of methanol to formic acid on Au8-: benchmark analysis based on completely renormalized coupled-cluster and density functional theory calculations

The left-eigenstate completely renormalized coupled-cluster (CC) method with singles, doubles, and noniterative triples [CR-CC(2,3)] and a few representative density functional theory (DFT) approaches have been applied to methanol oxidation to formic acid on a Au8(-) cluster, which is a model for aerobic oxidations on gold nanoparticles. It is demonstrated that CR-CC(2,3) supports the previous exothermic reaction mechanism, placing the initial rate-determining transition state, which corresponds to hydrogen transfer from the methoxy species to the molecular oxygen, at about 20 kcal/mol above the reactants, less than 40 kcal/mol above the O2 and CH3O(-) species coadsorbed on Au8(-), and considerably above the remaining two transition states along the reaction pathway. The DFT calculations using the previously exploited M06 hybrid functional show reasonable agreement with CR-CC(2,3), but B3LYP offers additional improvements in the description of the relevant activation energies. Pure functionals, including M06-L, BP86, and TPSS, do not work well, significantly underestimating the activation barriers, but dispersion corrections, as in B97-D, bring the results closer to the M06 accuracy level.

The shape of Au8: gold leaf or gold nugget?

The size at which nonplanar isomers of neutral, pristine gold nanoclusters become energetically favored over planar ones is still debated amongst theoreticians and experimentalists. Spectroscopy confirms planarity is preferred at sizes up to Au7, however, starting with Au8, the uncertainty remains for larger nanoclusters. Au8 computational studies have had different outcomes: the planar D4h "cloverleaf" isomer competes with the nonplanar Td, C2v and D2d "nugget" isomers for greatest energetic stability. We here examine the 2D vs. 3D preference in Au8 by presenting our own B2PLYP, MP2 and CCSD(T) calculations on these isomers: these methods afford a better treatment of long-range correlation, which is at the root of gold's characteristic aurophilicity. We then use findings from these high-accuracy computations to evaluate two less expensive DFT approaches, applicable to much larger nanoclusters: alongside the standard functional PBE, we consider M06-L (highly parametrized to incorporate long-range dispersive interactions). We find that increasing basis set size within the B2PLYP framework has a greater destabilizing effect on the nuggets than it has on the Au8 cloverleaf. Our CCSD(T) and B2PLYP predictions, replicated by DFT-PBE, all identify the cloverleaf as the most stable isomer; MP2 and DFT-M06-L show overestimation of aurophilicity, and favor, respectively, the nonplanar D2d and Td nuggets in its stead. We conclude that PBE, which more closely reproduces CCSD(T) findings, may be a better candidate density functional for the simulation of gold nanoclusters in this context.

Charge-dependent trends in structures and vibrational frequencies of [CO-Au-O2](q) (q = -1, 0, +1) complexes: evidence for cooperative interactions

Charge-dependent trends in the structures, vibrational frequencies, and binding energies of binary [AuCO](q) and [AuO(2)](q) and ternary [O(2)AuCO](q) complexes (q = -1,0,+1), have been investigated using density functional theory calculations. Three different geometrical motifs, given descriptive names of "separated", "pre-reactive", and "long-range", are identified for the ternary complexes. For the binary systems, the general trend is that the complexes become more diffuse as the charge becomes more negative, having longer intermolecular bond distances and weaker binding energies. The trends shown by the ternary complexes are more complicated, and are different for the various geometrical motifs. However, a general trend is that there is a cooperative interaction involving both the CO and O(2) with the Au center, which becomes more pronounced as the negative charge on the complexes increases from cationic to neutral to anionic. This cooperative interaction leads to increased electron density on the O(2) moiety, as is reflected in the bond-lengths and vibrational frequencies. Furthermore, it is found that for the pre-reactive complexes, the role of the Au center is to stabilize the formation of a conjugated π-system between the CO and O(2) molecules.

The binding of Ag+ and Au+ to ethene

For the reaction M(+)(C(2)H(4))(n-1) + C(2)H(4) --> M(+)(C(2)H(4))(n), where M = Ag, Au, the binding energies are predicted at the second order perturbation (MP2) and coupled cluster (CCSD(T)) levels of theory. As the basis set is systematically improved, the predicted M = Ag binding energies steadily improve, as compared to the experimental values. In fact, the complete basis set limit (CBS) predicted CCSD(T) binding energy for Ag(C(2)H(4))(+) is within experimental error. For MP2, as the basis set is improved, the agreement with experiment worsens. Gold ions are predicted to bind more strongly than silver ions to ethene ligands. Mulliken population analyses of the silver and gold systems exhibit delocalization of the positive charges of the metal ions onto the ethene ligands. Reduced variational space analysis indicates that electrostatic interactions are the principal contributor to the bonding in these systems. Multiconfigurational self-consistent field calculations do not support the Dewar-Chatt-Duncanson model of transition metal-alkene bonding in Au(C(2)H(4))(+).

Theory and simulation in heterogeneous gold catalysis

This critical review covers the application of quantum chemistry to the burgeoning area of the heterogeneous oxidation by Au. We focus on the most established reaction, the oxidation of CO at low temperature. The review begins with an overview of the methods available comparing the treatment of the electron-electron interaction and relativistic effects. The structure of Au particles and their interaction with oxide reviews is then discussed in detail. Calculations of the adsorption and reaction of CO and O2 are then considered and results from isolated and supported Au clusters compared (155 references).

CO oxidation catalyzed by single-walled helical gold nanotube

We study the catalytic capability of unsupported single-walled helical gold nanotubes Au(5,3) by using density functional theory. We use the CO oxidation as a benchmark probe to gain insights into high catalytic activity of the gold nanotubes. The CO oxidation, catalyzed by the Au(5,3) nanotube, proceeds via a two-step mechanism, CO + O2 --> CO2 +O and CO + O --> CO2. The CO oxidation is initiated by the CO + O2 --> OOCO --> CO2 + O reaction with an activation barrier of 0.29 eV. On the reaction path, a peroxo-type O-O-CO intermediate forms. Thereafter, the CO + O --> CO2 reaction proceeds along the reaction pathway with a very low barrier (0.03 eV). Note that the second reaction cannot be the starting point for the CO oxidation due to the energetically disfavored adsorption of free O2 on the gold nanotube. The high catalytic activity of the Au(5,3) nanotube can be attributed to the electronic resonance between electronic states of adsorbed intermediate species and Au atoms at the reaction site, particularly among the d states of Au atom and the antibonding 2pi* states of C-O and O1-O2, concomitant with a partial charge transfer. The presence of undercoordinated Au sites and the strain inherent in the helical gold nanotube also play important roles. Our study suggests that the CO oxidation catalyzed by the helical gold nanotubes is likely to occur at the room temperature.

Isomers of Au8

Using newly developed correlation consistent basis sets for gold, the relative energies for the neutral Au8 geometric isomers have been re-evaluated and the vertical ionization potentials calculated. The results using the correlation consistent basis sets show that second-order Moller-Plesset perturbation theory calculations strongly favor nonplanar Au8 structures for all basis sets that were employed. However, the general trend at the coupled cluster singles and doubles with perturbative triples level of theory is to increasingly favor planar structures as the basis set is improved. The effects of basis set and the effects of core-valence correlation are discussed.

Density functional study of the charge on Aun clusters (n=1–7) supported on a partially reduced rutile TiO2(110): Are all clusters negatively charged?

It is widely believed that small gold clusters supported on an oxide surface and adsorbed at the site of an oxygen vacancy are negatively charged. It has been suggested that this negative charge helps a gold cluster adsorb oxygen and weakens the O-O bond to make oxidation reactions more efficient. Given the fact that an oxygen vacancy is electron rich and that Au is a very electronegative element, the assumption that the Au cluster will take electron density from the vacancy is plausible. However, the density functional calculations presented here show that the situation is more complicated. The authors have used the Bader method to examine the charge redistribution when a Aun cluster (n=1-7) binds next to or at an oxygen vacancy on rutile TiO2(110). For the lowest energy isomers they find that Au1 and Au3 are negatively charged, Au5 and Au7 are positively charged, and Au2, Au4, and Au6 exchange practically no charge. The behavior of the Aun isomers having the second-lowest energy is also unexpected. Au2, Au3, Au5, and Au7 are negatively charged upon adsorption and very little charge is transferred when Au4 and Au6 are adsorbed. These observations can be explained in terms of the overlap between the frontier molecular orbitals of the gold cluster and the eigenstates of the support. Aun with even n becomes negatively charged when the lowest unoccupied molecular orbital has a lobe pointing in the direction of the oxygen vacancy or towards a fivefold coordinated Ti (5c-Ti) located in the surface layer; otherwise it stays neutral. Aun with odd n becomes negatively charged when the singly occupied molecular orbital has a lobe pointing in the direction of a 5c-Ti located at the vacancy site or in the surface layer, otherwise it donates electron density into the conduction band of rutile TiO2(110) becoming positively charged.

Origin of Oxide Sensitivity in Gold-Based Catalysts: A First Principle Study of CO Oxidation over Au Supported on Monoclinic and Tetragonal ZrO2

The catalytic performance of Au/oxide catalysts can vary significantly upon the change of oxide species or under different catalyst preparation conditions. Due to its complex nature, the physical origin of this phenomenon remains largely unknown. By extensive density functional theory calculations on a model system, CO oxidation on Au/ZrO2, this work demonstrates that the oxidation reaction is very sensitive to the oxide structure. The surface structure variation due to the transformation of the oxide phase or the creation of structural defects (e.g., steps) can greatly enhance the activity. We show that CO oxidation on typical Au/ZrO2 catalysts could be dominated by minority sites, such as monoclinic steps and tetragonal surfaces, the concentration of which is closely related to the size of oxide particle. Importantly, this variation in activity is difficult to understand following the traditional rules based on the O2 adsorption ability and the oxide reducibility. Instead, electronic structure analyses allow us to rationalize the results and point toward a general measure for CO + O2 activity, namely the p-bandwidth of O2, with important implications for Au/oxide catalysis.

CO adsorption on pure and binary-alloy gold clusters: a quantum chemical study

We performed density-functional theory analysis of nondissociative CO adsorption on 22 binary Au-alloy (Au(n)M(m)) clusters: n=0-3, m=0-3, and m+n=2 (dimers) or 3 (trimers), M=Cu/Ag/Pd/Pt. We report basis-set superposition error corrections to adsorption energies and include both internal energy of adsorption (DeltaU(ads)) and Gibbs free energy of adsorption (DeltaG(ads)) at standard conditions (298.15 K and 1 atm). We found onefold (atop) CO binding on all the clusters except Pd2 (twofold/bridged), Pt2 (twofold/bridged), and Pd3 (threefold). In agreement with the experimental results, we found that CO adsorption is thermodynamically favorable on pure Au/Cu clusters but not on pure Ag clusters and also observed the following adsorption affinity trend: Pd>Pt>Au>Cu>Ag. For alloy dimers we found the following patterns: Au2>M Au>M2 (M=Ag/Cu) and M2>M Au>Au2 (M=Pd/Pt). Alloying Ag/Cu dimers with (more reactive) Au enhanced adsorption and the opposite effect was observed for PdPt dimers. The Ag-Au, Cu-Au, and Pd-Au trimers followed the trends observed on dimers: Au3>M Au2>M2Au>M3 (M=Ag/Cu) and Pd3>Pd2Au>PdAu2>Au3. Interestingly, Pt-Au trimers reacted differently and alloying with Au systematically increased the adsorption affinity: PtAu2>Pt2Au>Pt3>Au3. A strikingly different behavior of Pt is also manifested by the triplet spin state and onefold (atop) binding in Pt3-CO which is in contradiction with the singlet spin state and threefold binding in Pd3-CO. We found a linear correlation between CO binding energy (BE) and elongation of the CO bond. For Ag-Au and Cu-Au clusters, the increase in CO BE (and elongation of the C-O bond which is probably due to the back donation) is accompanied by the decrease in the cluster-CO distance suggesting that the donation (from 5sigma highest occupied molecular orbital in CO to cluster lowest unoccupied molecular orbital) mechanism also contributes to the BE. For Pd-Au clusters, the cluster-CO distance (and CO bond length) increases with increase in the BE, suggesting that the donation mechanism may not be important for those clusters. No clear trend was observed for Pt-Au clusters.