Structural study of gold clusters

Adsorption of Multiple NO Molecules on Anionic Gold Clusters: Electron Transfer and Disproportionation Enhanced by Molecular Aggregation

The reactions of Aun- clusters with multiple nitric oxide (NO) molecules are explored at 150 K by utilizing a mini-flow-tube reactor and a time-of-flight mass spectrometer. Adsorption of multiple NO molecules is observed on most Aun-, while disproportionation reactions only occur on even-sized Aun- with n = 4, 6, 8, 20 and odd-sized ones with n = 5 and 7. Theoretical calculations reveal the geometric structures and electronic states of the products containing bimolecular and trimolecular NO units, where two NO molecules typically form dimers. Different from NO monomers that weakly interact with odd-sized Aun- and form electron-sharing covalent bonds with Au10-(D3h) and Au16-, NO dimers can extract significant charge from parent Aun-. Regarding the three NO molecules, a predilection toward condensation into trimers on even-sized Aun- is observed, while the tendency is more toward an adsorption pattern of a dimer plus a monomer on odd-sized Aun-. The NO trimers register even higher charge gain from Aun- as compared with the NO dimers, which leads to an elevated degree of activation and induces the progression of disproportionation reactions. Therefore, when considering the reaction between NO and Aun-, it appears that NO has a propensity to form dimers or trimers on Aun-. This behavior of aggregate formation substantially enhances the ability of NO to absorb negative charges from Aun- although the occurrence of disproportionate dissociation reactions is initiated only for specific sizes.

Mechanism and dynamics of Baeyer-Villiger oxidation of furfural to maleic anhydride in presence of H2O2 and Au clusters

CONTEXT

The increasing demand for fuels and chemicals in the world has prompted the exploration of various forms of renewable energy resources. Using C5-based furfural as the platform to replace the fossil energy resources is greatly attractive because of its abundance and environmental friendliness. Here we study the activity, selectivity, and possible reaction pathways for the Baeyer-Villiger oxidation of furfural over small Au clusters using hydrogen peroxide as oxidant. Furfural reacts with hydrogen peroxide in the presence of the catalysts with 93% selectivity towards maleic anhydride. Natural population analysis, frontier molecular orbital analysis, and spectroscopic analysis are used to illustrate the interaction mechanism between C5H4O2, H2O2, and Au. Reaction pathways leading to the formation of maleic anhydride are also explored. The reaction of C5H4O2 with H2O2 in the absence of a catalyst bears a relatively high transition state energy barrier of 2.98 eV for the first step involving absorption of H atom of H2O2 on the -OH group of C5H4O2. This is in agreement with the blank experiment where there were rare oxidation products observed in the absence of the metal cluster catalysts. On the other hand, transition state energies in the presence of the Au metal clusters are lower and the most feasible pathway is where the substrate and H2O2 co-bind on the Au catalyst and H2O2 molecule transfers an oxygen to the substrate, leading to the cleavage of the O-O bond.

METHODS

DFT calculations were done with B3PW91 functional. 6-311G(df, p) basis set was used for C, O, and H and aug-cc-pVDZ-PP was used for gold atoms. Gaussian 09 software was used for the calculations. Multiwfn 3.7 dev was used for the quantum theory of atoms-in-molecules (QTAIM) investigations.

Isolable Tetragold(0) Clusters with Polarity-Tunable exo-Au-Au Bond via Intramolecular σ-Aromatization

Intramolecular π-aromatization is a trait of many organic compounds that enhances the stability of their structures and polarizes related C-C π bonds. In contrast, rare study is focused on this phenomenon in metal clusters. Many existing homometallic clusters exhibit aromaticity, often characterized by nonpolar metal-metal bonds and a high degree of symmetry. However, synthesizing low-symmetric homometallic clusters with high-polar metal-metal bonds is challenging due to their limited thermodynamic stability. Herein, we report a facile strategy for the synthesis of [Au(μ2 -ER2 )]3 -AuPMe3 (E=Ge, Sn; R2 =1,1,4,4-tetrakis(trimethylsilyl)butane-1,4-diyl) clusters and reveal a novel stabilization mode, intramolecular σ-aromatization. Our electronic structure analyses show that these low-symmetric clusters possess a ten-electron σ-aromatic system, which is achieved via intramolecular σ-aromatization. Moreover, the strength of σ-aromaticity gives rise to a polarity-tunable exo-Au-Au bond.

Molecular adsorption on coinage metal subnanoclusters: A DFT+D3 investigation

Gold and silver subnanoclusters with few atoms are prominent candidates for catalysis-related applications, primarily because of the large fraction of lower-coordinated atoms exposed and ready to interact with external chemical species. However, an in-depth energetic analysis is necessary to characterize the relevant terms within the molecular adsorption process that can frame the interactions within the Sabatier principle. Herein, we investigate the interaction between Ag_n and Au_n subnanoclusters (clu, n = 2-7) and N₂ , NO, CO, and O₂ molecules, using scalar-relativistic density functional theory calculations within van der Waals D3 corrections. The onefold top site is preferred for all chemisorption cases, with a predominance of linear (≈180°) and bent (≈120°) molecular geometries. A larger magnitude of adsorption energy is correlated with smaller distances between molecules and clusters and with the weakening of the adsorbates bond strength represented by the increase of the equilibrium distances and decrease of molecular stretching frequencies. From the energetic decomposition, the interaction energy term was established as an excellent descriptor to classify subnanoclusters in the adsorption/desorption process concomitant with the Sabatier principle. The limiting cases: (i) weak molecular adsorption on the subnanoclusters, which may compromise the reaction activation, where an interaction energy magnitude close to 0 eV is observed (e.g., physisorption in N₂ /Ag₆ ); and (ii) strong molecular interactions with the subnanoclusters, given the interaction energy magnitude is larger than at least one of the individual fragment binding energies (e.g., strong chemisorption in CO/Au₄ and NO/Au₄ ), conferring a decrease in the desorption rate and an increase in the possible poisoning rate. However, the intermediate cases are promising by involving interaction energy magnitudes between zero and fragment binding energies. Following the molecular closed-shell (open-shell) electronic configuration, we find a predominant electrostatic (covalent) nature of the physical interactions for N₂ ⋯clu and CO ⋯clu (O₂ ⋯clu and NO⋯clu), except in the physisorption case (N₂ /Ag₆ ) where dispersive interaction is dominant. Our results clarify questions about the molecular adsorption on subnanoclusters as a relevant mechanistic step present in nanocatalytic reactions.

Efficient and accurate density-based prediction of macromolecular polarizabilities

Accurately and efficiently predicting macromolecules' polarizabilities is an open problem. In this work, we employ a few simple density-based quantities from the information-theoretic approach (ITA) to predict polarizability of proteins. We first build quantitative structure/property relationships between molecular polarizabilities and ITA quantities. We then verify the broad applicability of ITA quantities for polarizability prediction for inorganic, organic, and biological systems with both localized and delocalized electronic structure. As a proof-of-concept application, we predict the molecular polarizabilities of complex proteins. Based on the linear regression equations for 20 natural amino acid residues, 400 dipeptides, and 8000 tripeptides, one then predicts the molecular polarizability of a larger peptide or even a protein once the molecular wavefunction is obtained. Because it is extremely costly to determine the wavefunction for a macromolecule like a protein, we propose to combine the ITA with the linear-scaling generalized energy-based fragmentation (GEBF) method to predict the macromolecular polarizability. In GEBF, the total molecular polarizability is obtained as a linear combination of the corresponding quantities from a series of small subsystems. We can predict them based on the subsystem wavefunction and linear regression equations rather than compute them from the nearly-intractable coupled-perturbed Hartree-Fock or Kohn-Sham equations for the whole macromolecule. Computational results showcase that the GEBF-ITA protocol should be an inexpensive yet accurate theoretical tool for predicting macromolecular polarizabilities.

The adsorption of single Au atom and nucleation on γ-Al2O3 surfaces

Single-atom catalysts (SACs) in heterogeneous catalysts have attracted increasing attention and the adsorption and nucleation of single atom on the surface are closely related to the performance of the catalyst. The present work employed density functional theory calculations to examine the adsorption of single Au atom and nucleation on γ-Al₂O₃ surfaces at the atomic level. The effect of surface hydroxyls group on the adsorption and nucleation of single Au atom on γ-Al₂O₃ surfaces is explored. It was found that the spillover reactions of surface hydroxyls H atoms with the deposited Au_- are not available on the hydroxylated surface. The interaction of Au to the clean surface is the stronger than to the hydroxylated surface. The even-odd alternations of Au_x and weak binding of single Au atoms to γ-Al₂O₃ leads to large even-numbered Au cluster on the surface. Density of states and electron density difference analysis show that the electronic structure of Au/γ-Al₂O₃ is quite different from the reported Cu and Pd on Al₂O₃.

Phosphine-Stabilized Hidden Ground States in Gold Clusters Investigated via a Aun(PH3)m Database

Nanoclusters are promising materials for catalysis and sensing due to their large surface areas and unique electronic structures which can be tailored through composition, geometry, and chemistry. However, relationships correlating synthesis parameters directly to outcomes are limited. While previous computational studies have mapped the potential energy surface of specific systems of bare nanoclusters by generating and calculating the energies of reasonable structures, it is known that environmental ions and ligands crucially impact the final shape and size. In this work, phosphine-stabilized gold is considered as a test system and DFT calculations are performed for clusters with and without ligands, producing a database containing >10000 structures for Aun(PH3)m (n ≤ 12). We find that the ligation of phosphines affects the thermodynamic stability, bonding, and electronic structure of Au nanoclusters, specifically such that "hidden" ground state cluster geometries are stabilized that are dynamically unstable in the pure gold system. Further, the addition of phosphine introduces steric effects that induce a transition from planar to nonplanar structures at 4-5 Au atoms rather than up to 13-14 Au atoms, as previously predicted for bare clusters. This work highlights the importance of considering the ligand environment in the prediction of nanocluster morphology and functionality, which adds complexity as well as a rich opportunity for tunability.

Synthesis of Self-Assembled Single Atomic Layer Gold Crystals-Goldene

We report, for the first time, a technique to synthesize free-standing, one-atom thick 2D gold crystals (namely, goldene) and self-assembled 2D periodic arrays of goldene. High-resolution transmission electron microscopy (HRTEM) imaging of goldene revealed herringbone and honeycomb lattices, which are primarily gold surface features due to its reconstruction. Imaging of these surface-only features by a nonsurface characterization technique such as HRTEM is an unequivocal proof of the absence of three-dimensionality in goldene. Atomic force microscopy confirmed 1-2 Å thickness of goldene. High-resolution X-ray photoelectron spectroscopy (HR-XPS), selective area electron diffraction, and energy-dispersive X-ray spectroscopy confirmed the chemical identity of goldene. We discovered the phenomenon of electric field-induced self-assembly of goldene supracrystals with a herringbone structure and developed an electric field printing (e-print) technique for goldene arrays. Goldene showed a semiconductor response with a knee voltage of ∼3.2 V, and I/V spectroscopy revealed periodic room temperature Coulomb blockade oscillations. These observations are consistent with the theoretical calculations reported in the literature predicting enhanced Coulombic interactions between gold valence electrons and the nucleus in stable 2D gold. Goldene exhibited multiple, intense, and well-resolved optical absorption peaks and several fine bands across the UV-vis region, and we calculated its optical band gap to be 3.59 eV. Magnetic force microscopy measurements of goldene periodic arrays showed a ∼5 mV peak amplitude confirming its ferromagnetism. Optical and magnetic properties of goldene are consistent with those reported in the literature for 2D planar gold clusters with less than 12 atoms.

The effect of different energy portions on the 2D/3D stability swapping for 13-atom metal clusters

The complexity of Cu₁₃, Ag₁₃, and Au₁₃ coinage-metal clusters was investigated through their energy contributions via a density functional theory study, considering improvements in the PBE functional, such as van der Waals (vdW) corrections, spin-orbit coupling (SOC), Hubbard term (+U), and their combinations. Investigating two-dimensional (planar 2D) and three-dimensional (distorted 3D, CUB - cuboctahedral, and ICO - icosahedral) configurations, we found that vdW corrections are dominant in modulating the stability swapping between 2D and ICO (3D) for Ag₁₃ (Au₁₃), whereas for Cu₁₃ its role is increasing the relative stability between 2D (least stable) and 3D (most stable), setting ICO as the reference. Among the energy portions that constitute the relative total energy, the dimensionality difference correlates with the magnitude of the relative dispersion energy (large for 2D/ICO and small for 3D/ICO) as the causal factor responsible for an eventual stability swapping. For instance, empirical vdW corrections may favor Ag₁₃ as ICO, while semi empirical ones tend to swap the stability by favoring 2D. The same tendency is observed for Au₁₃, except when SOC is included, which enlarges the stability of 3D over 2D. Energy decomposition analysis combined with the natural orbitals for the chemical valence approach confirmed the correlations between the dimensionality difference and the magnitude of the relative dispersion energies. Our structural analysis protocol was able to capture the local distortion effects (or even their absence) through the quantification of the Hausdorff chirality measure. Here, ICO, CUB, and 2D are achiral configurations for all coinage-metal clusters, whereas Cu₁₃ as 3D presents a slight chirality when vdW correction based on many body dispersion is used, at the same time Ag₁₃ as 3D turned out to be chiral for all calculation protocols as evidence of the role of the chemical composition.

Rapid electrochemical recognition of trimethoprim in human urine samples using new modified electrodes (CPE/Ag/Au NPs) analysing tunable electrode properties: experimental and theoretical studies

Pharmaceutical effluents are a serious environmental issue, which require to be treated by a suitable technique; thus, the electrochemical process is actively considered as a viable method for the treatment. In this work, new carbon paste electrodes (CPEs) were fabricated by compressing gold and silver nanoparticles (NPs), namely, CPE/Ag NPs, CPE/Au NPs, and CPE/Ag/Au NPs and then completely characterized by different analytical methods. The performance of the electrodes was studied after determining their surface area (×10^-6 cm²) as 4.17, 5.05, 5.27, and 5.12, producing high anodic currents for K₄[Fe(CN)₆] compared to the commercial electrode. This agrees with the results of impedance study, where the electron transfer rate constants (k_app, ×10^-3 cm s^-1) were determined to be 28.7, 42.6, 41.0, and 101.4 for CPE, CPE/Ag NPs, CPE/Au NPs, and CPE/Ag/Au NPs, respectively, through the Bode plot-phase shifts. This is consistent with the charge transfer resistance (R_CT, Ω), resulting as 171 for CPE/Ag/Au NPs < 395 for CPE/Ag NPs < 427 for CPE/Au NPs and < 742 for CPE. Therefore, these electrodes were employed to detect trimethoprim (TMP) since metallic NPs contribute good crystallinity, stability, conduciveness, and surface plasmon resonance to the CPE, convalescing the sensitivity; comprehensively, they were applied for its detection in real water and human urine samples, and the limit of detection (LOD) was as low as 0.026, 0.032, and 0.026 μmol L^-1 for CPE/Ag NPs, CPE/Au NPs, and CPE/Ag/Au NPs, respectively. In contrast, unmodified CPE was unable to detect TMP due to the lack of efficiency. The developed technique shows excellent electrochemical recovery of 92.3 and 97.1% in the urine sample. Density functional theory (DFT) was used to explain the impact of the metallic center in graphite through density of states (DOS).

Assessment of the van der Waals, Hubbard U parameter and spin-orbit coupling corrections on the 2D/3D structures from metal gold congeners clusters

The coinage-metal clusters possess a natural complexity in their theoretical treatment that may be accompanied by inherent shortcomings in the methodological approach. Herein, we performed a scalar-relativistic density functional theory study, considering Perdew, Burke, and Ernzerhof (PBE) with (empirical and semi empirical) van der Waals (vdW), spin-orbit coupling (SOC), +U (Hubbard term), and their combinations, to treat the Cu 13 , Ag 13 , and Au 13 clusters in different structural motifs. The energetic scenario is given by the confirmation of the 3D lowest energy configurations for Cu 13 and Ag 13 within all approaches, while for Au 13 there is a 2D/3D competition, depending on the applied correction. The 2D geometry is 0.43 eV more stable with plain PBE than the 3D one, the SOC, +U, and/or vdW inclusion decreases the overestimated stability of the planar configurations, where the most surprising result is found by the D3 and D3BJ vdW corrections, for which the 3D configuration is 0.29 and 0.11 eV, respectively, more stable than the 2D geometry (with even higher values when SOC and/or +U are added). The D3 dispersion correction represents 7.9% (4.4%) of the total binding energy for the 3D (2D) configuration, (not) being enough to change the sd hybridization and the position of the occupied d -states. Our predictions are in agreement with experimental results and in line with the best results obtained for bulk systems, as well as with hybrid functionals within D3 corrections. The properties description undergoes small corrections with the different approaches, where general trends are maintained, that is, the average bond length is smaller (larger) for lower (higher)-coordinated structures, since a same number of electrons are shared by a smaller (larger) number of bonds, consequently, the bonds are stronger (weaker) and shorter (longer) and the sd hybridization index is larger (smaller). Thus, Au has a distinct behavior in relation to its lighter congeners, with a complex potential energy surface, where in addition to the relevant relativistic effects, correlation and dispersion effects must also be considered.

Systematic Investigation of the Structure, Stability, and Spin Magnetic Moment of CrM n Clusters (M = Cu, Ag, Au, and n = 2-20) by DFT Calculations

Binary clusters of transition-metal and noble-metal elements have been gathering momentum for not only advanced fundamental understanding but also potential as elementary blocks of novel nanostructured materials. In this regard, the geometries, electronic structures, stability, and magnetic properties of Cr-doped Cu _n , Ag _n , and Au _n clusters (n = 2-20) have been systematically studied by means of density functional theory calculations. It is found that the structural evolutions of CrCu _n and CrAg _n clusters are identical. The icosahedral CrCu₁₂ and CrAg₁₂ are crucial sizes for doped copper and silver species. Small CrAu _n clusters prefer the planar geometries, while the larger ones appear as on the way to establish the tetrahedral CrAu₁₉. Our results show that while each noble atom contributes one s valence electron to the cluster shell, the number of chromium delocalized electrons is strongly size-dependent. The localization and delocalization behavior of 3d orbitals of the chromium decide how they participate in metallic bonding, stabilize the cluster, and give rise to and eventually quench the spin magnetic moment. Moreover, molecular orbital analysis in combination with a qualitative interpretation using the phenomenological shell model is applied to reveal the complex interplay between geometric structure, electronic structure, and magnetic moment of clusters. The finding results are expected to provide greater insight into how a host material electronic structure influences the geometry, stability, and formation of spin magnetic moments in doped systems.

Quantifying electron-correlation effects in small coinage-metal clusters via ab initio calculations

We investigate many-electron correlation effects in neutral and charged coinage-metal clusters Cun, Agn, and Aun (n = 1-4) via ab initio calculations using fixed-node diffusion Monte Carlo (FN-DMC) simulations, density functional theory (DFT), and the Hartree-Fock (HF) method. From very accurate FN-DMC total energies of the clusters and the HF results in the infinity large complete-basis-set limit, we obtain correlation energies in these strongly correlated many-electron clusters involving d orbitals. The obtained bond lengths of the clusters, atomic binding and dissociation energies, ionization potentials, and electron affinities are in satisfactory agreement with the available experiments. In the analysis, the electron correlation effects on these observable physical quantities are quantified by relative correlation contributions determined by the difference between the calculated FN-DMC and HF results. We show that the correlation contribution is not only significant for the quantities related to electronic structures of the coinage-metal clusters, such as electron affinity, but it is also essential for the stability of the atomic structures of these clusters. For example, the electron correlation contribution is responsible for more than 90% of the atomic binding energies of the small neutral copper clusters. We also demonstrate the orbital-occupation dependence of the correlation energy and electron pairing of the valence electrons in these coinage-metal clusters from the electron correlation-energy gain and spin-multiplicity change in the electron addition processes, which are reflected in their ionization potentials and electron affinities.

Gold nanostructures on iron oxide surfaces and their interaction with CO

We review results of density functional theory calculations of the adsorption of single gold atoms and formation of sub-nanometer Aunstructures (n= 2 to 5) on most stable iron oxide surfaces: hematite (0001), and magnetite (111) and (001). Structural, energetic, and electronic properties of Aunstructures on both Fe- and O-rich oxide terminations are discussed. Different chemical character of the two oxide terminations is reflected in distinctly stronger binding of gold at the oxygen- than at the iron-terminated surface, and in different changes of the adsorption binding energy with the size of the Auncluster. On the iron-terminated oxide surface the binding energy increases whereas on the oxygen-rich termination it decreases with the number of Au atoms in the structure. Upon CO adsorption on magnetite surface all Aunstructures have a net positive charge and CO binds to the most cationic Au atom of a cluster. Interactions of Aunand CO with magnetite (111) show many similarities with those on hematite (0001) surface. The influence of the substrate relaxation effects on adsorption energy is also discussed.

Effect of Gold Electronic State on the Catalytic Performance of Nano Gold Catalysts in n-Octanol Oxidation

This study aims to identify the role of the various electronic states of gold in the catalytic behavior of Au/M_xO_y/TiO₂ (where M_xO_y are Fe₂O₃ or MgO) for the liquid phase oxidation of n-octanol, under mild conditions. For this purpose, Au/M_xO_y/TiO₂ catalysts were prepared by deposition-precipitation with urea, varying the gold content (0.5 or 4 wt.%) and pretreatment conditions (H₂ or O₂), and characterized by low temperature nitrogen adsorption-desorption, X-ray powder diffraction (XRD), energy dispersive spectroscopy (EDX), scanning transmission electron microscopy-high angle annular dark field (STEM HAADF), diffuse reflectance Fourier transform infrared (DRIFT) spectroscopy of CO adsorption, temperature-programmable desorption (TPD) of ammonia and carbon dioxide, and X-ray photoelectron spectroscopy (XPS). Three states of gold were identified on the surface of the catalysts, Au⁰, Au¹⁺ and Au³⁺, and their ratio determined the catalysts performance. Based on a comparison of catalytic and spectroscopic results, it may be concluded that Au⁺ was the active site state, while Au⁰ had negative effect, due to a partial blocking of Au⁰ by solvent. Au³⁺ also inhibited the oxidation process, due to the strong adsorption of the solvent and/or water formed during the reaction. Density functional theory (DFT) simulations confirmed these suggestions. The dependence of selectivity on the ratio of Brønsted acid centers to Brønsted basic centers was revealed.

Prediction of Structures and Atomization Energies of Coinage Metals, (M)n, n < 20: Extrapolation of Normalized Clustering Energies to Predict the Cohesive Energy

The geometries of the group 11 coinage metals (n = 2-20) were optimized to determine the lowest energy isomers for each cluster size, singlets for even numbers and doublets for odd numbers. For copper and silver, 2-D (planar) geometries were favored up to n = 6. For gold, 2D (planar) geometries were favored up to n = 13. Normalized clustering energies were plotted as a function of cluster size (n^-1/3, for n = 4-20) with various DFT functionals and the CCSD(T)-F12b method and were extrapolated to predict the bulk cohesive energy. In the case of copper and silver, there is excellent agreement between the cohesive energies predicted at the CCSD(T)-F12b level of theory and the experimental values. For gold, the CCSD(T)-F12b values needed to be corrected for spin-orbit relativistic effects to obtain good agreement with experiment. Electronic properties including the HOMO-LUMO gaps for the even clusters and the spin densities for the odd clusters were calculated. The lowest gap is predicted to occur for n = 16 where the HOMO and LUMO are very similar in shape.

Energy landscape of Au13: a global view of structure transformation

It has long been a challenge in physics and chemistry to acquire a global picture of the energy landscape of a specific material, as well as the kinetic transformation process between configurations of interest. Here we have presented a comprehensive approach to deal with the structure transformation problem, along with the illustration of the energy landscape, as exemplified with the case of Au13. A configuration space based on interatomic distances was proposed and demonstrated to have a strong correlation between structure and energy, with application in structure analysis to screen for trial transition pathways. As several representative configurations and their transition pathways ascertained and by projecting on a plane, a visual two-dimensional contour map was sketched revealing the unique energy landscape of Au13. It shows that the 2D and 3D clusters form two funnels in the high-dimensional configuration space, with a transition pathway with a 0.976 eV barrier bridging them.

C. Ribeiro A, Kolobova E, M. D. R. S. Martins L, A. C. Carabineiro S, German D, Pichugina D, Jiang C, J. L. Pombeiro A, Bogdanchikova N, Cortés Corberán V, Pestryakov A. Supported Gold Nanoparticles as Catalysts in Peroxidative and Aerobic Oxidation of 1-Phenylethanol under Mild Conditions

The efficiency of Au/TiO₂ based catalysts in 1-phenylethanol oxidation was investigated. The role of support modifiers (La₂O₃ or CeO₂), influence of gold loading (0.5% or 4%) and redox pretreatment atmosphere, catalyst recyclability, effect of oxidant: tert-butyl hydroperoxide (TBHP) or O₂, as well as the optimization of experimental parameters of the reaction conditions in the oxidation of this alcohol were studied and compared with previous studies on 1-octanol oxidation. Samples were characterized by temperature-programmed oxygen desorption (O₂-TPD) method. X-ray photoelectron spectroscopy (XPS) measurements were carried out for used catalysts to find out the reason for deactivation in 1-phenylethanol oxidation. The best catalytic characteristics were shown by catalysts modified with La₂O₃, regardless of the alcohol and the type of oxidant. When O₂ was used, the catalysts with 0.5% Au, after oxidative pretreatment, showed the highest activity in both reactions. The most active catalysts in 1-phenylethanol oxidation with TBHP were those with 4% Au and the H₂ treatment, while under the same reaction conditions, 0.5% Au and O₂ treatment were beneficial in 1-octanol oxidation. Despite the different chemical nature of the substrates, it seems likely that Au⁺(Au^δ+) act as the active sites in both oxidative reactions. Density functional theory (DFT) simulations confirmed that the gold cationic sites play an essential role in 1-phenylethanol adsorption.

Physico-Chemical Insights into Gas-Phase and Oxide-Supported Sub-Nanometre AuCu Clusters

Catalysis by AuCu nanoclusters is a promising scientific field. However, our fundamental understanding of the underlying mechanisms of mixing in AuCu clusters at the sub-nanometre scale and their physico-chemical properties in both the gas-phase and on oxide supports is limited. We have identified the global minima of gas-phase and MgO(100)-supported AuCu clusters with 3–10 atoms using the Mexican Enhanced Genetic Algorithm coupled with density functional theory. Au and Cu adatoms and supported dimers have been also simulated at the same level of theory. The most stable composition, as calculated from mixing and binding energies, is obtained when the Cu proportion is close to 50%. The structures of the most stable free AuCu clusters exhibit Cu-core/Au-shell segregation. On the MgO surface however, there is a preference for Cu atoms to lie at the cluster-substrate interface. Due to the interplay between the number of interfacial Cu atoms and surface-induced cluster rearrangement, on the MgO surface 3D structures become more stable than 2D structures. The O-site of MgO surface is found to be the most favourable adsorption site for both metals. All dimers favour vertical (V) configurations on the surface and their adsorption energies are in the order: AuCu < CuCu < AuAu < AuCu (where the underlined atom is bound to the O-site). For both adatoms and AuCu dimers, adsorption via Cu is more favourable than Au-adsorbed configurations, but, this disagrees with the ordering for the pure dimers due to a combination of electron transfer and the metal-on-top effect. Binding energy (and second difference) and HOMO-LUMO gap calculations show that even-atom (even-electron) clusters are more stable than the neighbouring odd-atom (odd- electron) clusters, which is expected for closed- and open-shell systems. Supporting AuCu clusters on the MgO(100) surface decreases the charge transfer between Au and Cu atoms calculated in free clusters. The results of this study may serve as a foundation for designing better AuCu catalysts.

The influence of support materials on the structural and electronic properties of gold nanoparticles – a DFT study

This study investigates the effect of commonly used support materials (MgO, C, CeO₂) on small gold particles using dispersion corrected density functional theory (DFT-D).

Aromatic character of [Au13]5+ and [MAu12]4+/6+ (M = Pd, Pt) cores in ligand protected gold nanoclusters – interplay between spherical and planar σ-aromatics

Ligand-protected superatoms are able to behave as both spherical and planar aromatic species, providing a strong link between spherical and planar σ-aromatics, which can be controlled selectively by tuning their redox charge states.

Density-functional tight-binding approach for metal clusters, nanoparticles, surfaces and bulk: application to silver and gold

Density-functional based tight-binding (DFTB) is an efficient quantum mechanical method that can describe a variety of systems, going from organic and inorganic compounds to metallic and hybrid materials. The present topical review addresses the ability and performance of DFTB to investigate energetic, structural, spectroscopic and dynamical properties of gold and silver materials. After a brief overview of the theoretical basis of DFTB, its parametrization and its transferability, we report its past and recent applications to gold and silver systems, including small clusters, nanoparticles, bulk and surfaces, bare and interacting with various organic and inorganic compounds. The range of applications covered by those studies goes from plasmonics and molecular electronics, to energy conversion and surface chemistry. Finally, perspectives of DFTB in the field of gold and silver surfaces and NPs are outlined.

Au147 nanoparticles: Ordered or amorphous?

Structural aspects of the Au₁₄₇ cluster have been investigated through a density functional based tight binding global optimization involving a parallel tempering molecular dynamics scheme with quenching followed by geometries relaxation at the Density Functional Theory (DFT) level. The focus is put on the competition between relaxed ordered regular geometries and disordered (or amorphous) structures. The present work shows that Au₁₄₇ amorphous geometries are relevant low energy candidates and are likely to contribute in finite temperature dynamics and thermodynamics. The structure of the amorphous-like isomers is discussed from the anisotropy parameters, the atomic coordinations, the radial and pair distribution functions, the IR spectra, and the vibrational DOS. With respect to the regular structures, the amorphous geometries are shown to be characterized by a larger number of surface atoms, a less dense volume with reduced coordination number per atom, a propensity to increase the dimension of flat facets at the surface, and a stronger anisotropy. Moreover, all amorphous clusters have similar IR spectra, almost continuous with active frequencies over the whole spectral range, while symmetric clusters are characterized by a few lines with large intensities.

Optimum Particle Size for Gold-Catalyzed CO Oxidation

The structure sensitivity of gold-catalyzed CO oxidation is presented by analyzing in detail the dependence of CO oxidation rate on particle size. Clusters with less than 14 gold atoms adopt a planar structure, whereas larger ones adopt a three-dimensional structure. The CO and O₂ adsorption properties depend strongly on particle structure and size. All of the reaction barriers relevant to CO oxidation display linear scaling relationships with CO and O₂ binding strengths as main reactivity descriptors. Planar and three-dimensional gold clusters exhibit different linear scaling relationship due to different surface topologies and different coordination numbers of the surface atoms. On the basis of these linear scaling relationships, first-principles microkinetics simulations were conducted to determine CO oxidation rates and possible rate-determining step of Au particles. Planar Au₉ and three-dimensional Au₇₉ clusters present the highest CO oxidation rates for planar and three-dimensional clusters, respectively. The planar Au₉ cluster is much more active than the optimum Au₇₉ cluster. A common feature of optimum CO oxidation performance is the intermediate binding strengths of CO and O₂, resulting in intermediate coverages of CO, O₂, and O. Both these optimum particles present lower performance than maximum Sabatier performance, indicating that there is sufficient room for improvement of gold catalysts for CO oxidation.

Benchmark Study of Density Functional Theory for Neutral Gold Clusters, Aun (n = 2-8)

Neutral gold clusters, Au_n (n = 2-8), were optimized using coupled cluster singles and doubles with perturbative triples (CCSD(T)) with a triple-ζ-level basis set to develop reliable reference values for their structural and energy parameters in order to assess the performance of density functionals. The performance of 44 density functional theory (DFT) methods for calculating molecular structures and relative energies is assessed with respect to CCSD(T). In addition, their performance when calculating vertical ionization potentials (vIPs) of Au_n (n = 2-8) is also assessed by comparison with experimental values. The revTPSS functional shows good performance for calculating both the structural and energy properties of Au_n (n = 2-8), whereas B3P86 shows a remarkable performance in calculating the vIPs. The quadruple-ζ-level valence basis set is necessary for obtaining accurate energy values in CCSD(T) calculations.

Small gold clusters catalyzing the conversion of glycerol to epichlorohydrin

The conversion of glycerol to epichlorohydrin (GTE) is of great interest because the product is widely used in plastics, rubbers and adhesives, and also contributes to the disposal of the reactant glycerol, a major by-product in biodiesel production.

Unraveling the Planar-Globular Transition in Gold Nanoclusters through Evolutionary Search

Au nanoclusters are of technological relevance for catalysis, photonics, sensors, and of fundamental scientific interest owing to planar to globular structural transformation at an anomalously high number of atoms i.e. in the range 12–14. The nature and causes of this transition remain a mystery. In order to unravel this conundrum, high throughput density functional theory (DFT) calculations, coupled with a global structural optimization scheme based on a modified genetic algorithm (GA) are conducted. More than 20,000 Au₁₂, Au₁₃, and Au₁₄ nanoclusters are evaluated. With any DFT functional, globular and planar structures coexist across the size range of interest. The planar-globular transition is gradual at room temperature rather than a sharp transition as previously believed. The effects of anionicity, s-d band hybridization and long range interactions on the dimensional transition are quantified by using the structures adjacent to the minima. Anionicity marginally changes the relative stability of the clusters. The degree of s-d hybridization is varied via changing the Hubbard U value which corroborate that s-d hybridization alone does not stabilize planar structures. van der Waals interactions, on the other hand, stabilize globular structures. These results elucidate the balance between the different reasons of the dimensional transition in gold nanoclusters.

Analysis of Different Fragmentation Strategies on a Variety of Large Peptides: Implementation of a Low Level of Theory in Fragment-Based Methods Can Be a Crucial Factor

We have investigated the performance of two classes of fragmentation methods developed in our group (Molecules-in-Molecules (MIM) and Many-Overlapping-Body (MOB) expansion), to reproduce the unfragmented MP2 energies on a test set composed of 10 small to large biomolecules. They have also been assessed to recover the relative energies of different motifs of the acetyl(ala)18NH2 system. Performance of different bond-cutting environments and the use of Hartree-Fock and different density functionals (as a low level of theory) in conjunction with the fragmentation strategies have been analyzed. Our investigation shows that while a low level of theory (for recovering long-range interactions) may not be necessary for small peptides, it provides a very effective strategy to accurately reproduce the total and relative energies of larger peptides such as the different motifs of the acetyl(ala)18NH2 system. Employing M06-2X as the low level of theory, the calculated mean total energy deviation (maximum deviation) in the total MP2 energies for the 10 molecules in the test set at MIM(d=3.5Å), MIM(η=9), and MOB(d=5Å) are 1.16 (2.31), 0.72 (1.87), and 0.43 (2.02) kcal/mol, respectively. The excellent performance suggests that such fragment-based methods should be of general use for the computation of accurate energies of large biomolecular systems.

Modeling the atomistic growth behavior of gold nanoparticles in solution

The properties of gold nanoparticles strongly depend on their three-dimensional atomic structure, leading to an increased emphasis on controlling and predicting nanoparticle structural evolution during the synthesis process. In order to provide this atomistic-level insight and establish a link to the experimentally-observed growth behavior, a kinetic Monte Carlo simulation (KMC) approach is developed for capturing Au nanoparticle growth characteristics. The advantage of this approach is that, compared to traditional molecular dynamics simulations, the atomistic nanoparticle structural evolution can be tracked on time scales that approach the actual experiments. This has enabled several different comparisons against experimental benchmarks, and it has helped transition the KMC simulations from a hypothetical toy model into a more experimentally-relevant test-bed. The model is initially parameterized by performing a series of automated comparisons of Au nanoparticle growth curves versus the experimental observations, and then the refined model allows for detailed structural analysis of the nanoparticle growth behavior. Although the Au nanoparticles are roughly spherical, the maximum/minimum dimensions deviate from the average by approximately 12.5%, which is consistent with the corresponding experiments. Also, a surface texture analysis highlights the changes in the surface structure as a function of time. While the nanoparticles show similar surface structures throughout the growth process, there can be some significant differences during the initial growth at different synthesis conditions.

CO adsorption on small Aun (n = 1-4) structures supported on hematite. I. Adsorption on iron terminated α-Fe2O3 (0001) surface

This is the first of two papers dealing with the adsorption of Au and formation of Aun nanostructures (n = 1-4) on hematite (0001) surface and adsorption of CO thereon. The stoichiometric Fe-terminated (0001) surface of hematite was investigated using density functional theory in the generalized gradient approximation of Perdew-Burke-Ernzerhof (PBE) form with Hubbard correction U, accounting for strong electron correlations (PBE+U). The structural, energetic, and electronic properties of the systems studied were examined for vertical and flattened configurations of Aun nanostructures adsorbed on the hematite surfaces. The flattened ones, which can be viewed as bilayer-like structures, were found energetically more favored than vertical ones. For both classes of structures the adsorption binding energy increases with the number of Au atoms in a structure. The adsorption of Aun induces charge rearrangement at the Aun/oxide contact which is reflected in work function changes. In most considered cases Aun adsorption increases the work function. A detailed analysis of the bonding electron charge is presented and the corresponding electron charge rearrangements at the contacts were quantified by a Bader charge analyses. The interaction of a CO molecule with the Aun nanostructures supported on α-Fe2O3 (0001) and the oxide support was studied. It is found that the CO adsorption binding to the hematite supported Aun structures is more than twice as strong as to the bare hematite surface. Analysis of the Bader charges on the atoms showed that in each case CO binds to the most positively charged (cationic) atom of the Aun structure. Changes in the electronic structure of the Aun species and of the oxide support, and their consequences for the interactions with CO, are discussed.