Contrasting motif preferences of platinum and gold nanoclusters between 55 and 309 atoms

Characterisation of the morphology of surface-assembled Au nanoclusters on amorphous carbon

In this study, aberration-corrected scanning transmission electron microscopy is employed to investigate the morphology of Au clusters formed from the aggregation of single atoms sputtered onto an amorphous carbon surface. The morphologies of surface-assembled clusters of N > 100 atoms are referenced against the morphologies of size-selected clusters determined from previously published results. We observe that surface-assembled clusters (at the conditions employed here) are approximately spherical in shape. The structural isomers of the imaged clusters have also been identified, and the distribution of structural types is broadly in agreement with those from size-selected cluster deposition sources. For clusters of approximately 147 atoms, we find a preference for icosahedra over decahedra and truncated octahedra, but at this size there is a high proportion of unidentified/amorphous structures. At around 309 atoms, we find a preference for decahedra over icosahedra and truncated octahedra, but over half the structures remain unidentifiable/amorphous. For sizes above approximately 561 atoms we are able to identify most of the structures, and find decahedra are still the most favoured, although in competition with single-crystal fcc morphologies. The similarity in structure between surface-assembled and size-selected clusters from a cluster source provides evidence of the relevance of size-selected cluster studies to clusters synthesised by other, industrially relevant, methodologies.

Calculations of the effect of catalyst size and structure on the electrocatalytic reduction of CO2 on Cu nanoclusters

The structure and catalytic properties of Cu nanoclusters of sizes between 55 and 147 atoms were examined to understand if small Cu clusters could provide enhancement over traditional catalysts for the electrocatalysis of CO2 to CO and carbon-based fuels, such as CH4 and CH3OH, compared to bulk Cu surfaces and large Cu nanoparticles. Clusters studied included Cu55, Cu78, Cu101, Cu124, and Cu147, the structures of which were determined using global optimisation. The majority of Cu clusters examined were icosahedral, including the perfect closed-shell, partial-shell, elongated and distorted icosahedral clusters. Free energy diagrams for the reduction of CO2 showed the potential required for the formation of CO is notably smaller for all cluster sizes considered, relative to Cu(111). Less variation is observed for the limiting potential for the formation of CH4 and CH3OH. However, it was found that clusters that are either a distorted motif or contain vacancy defects yielded the best activity and provide an interesting synthesis target for future experiments.

Comparison of Taboo Search Methods for Atomic Cluster Global Optimization with a Basin-Hopping Algorithm

The basin-hopping algorithm (BHA) allows for the efficient exploration of atomic cluster potential energy surfaces by random perturbations in configuration space, followed by energy minimizations. Here, the taboo search method is incorporated to prevent the search from revisiting recently visited regions of the search space. Two taboo search modes are implemented, one mode resets the search to random coordinates upon encountering the taboo region, while the other simply rejects any proposed move into the taboo region. These two modes are tested and compared on a variety of potential energy surfaces─several clusters where atomic interactions are described by the Lennard-Jones potential, and Au55 where a semi-empirical tight binding potential is used to describe atomic interactions. Some differences in performance between the two taboo search modes were noted for LJ38 and Au55, with the mode that rejects all hops into the taboo region performing better, offering a means to improve the efficiency of the BHA for multifunnel systems. However, both taboo search modes failed to significantly improve performance on multifunnel systems where more than two funnels were present in the system.

Theory of Anisotropic Metal Nanostructures

A significant challenge in the development of functional materials is understanding the growth and transformations of anisotropic colloidal metal nanocrystals. Theory and simulations can aid in the development and understanding of anisotropic nanocrystal syntheses. The focus of this review is on how results from first-principles calculations and classical techniques, such as Monte Carlo and molecular dynamics simulations, have been integrated into multiscale theoretical predictions useful in understanding shape-selective nanocrystal syntheses. Also, examples are discussed in which machine learning has been useful in this field. There are many areas at the frontier in condensed matter theory and simulation that are or could be beneficial in this area and these prospects for future progress are discussed.

New Magic Au24 Cluster Stabilized by PVP: Selective Formation, Atomic Structure, and Oxidation Catalysis

An unprecedented magic number cluster, Au24Cl x (x = 0-3), was selectively synthesized by the kinetically controlled reduction of the Au precursor ions in a microfluidic mixer in the presence of a large excess of poly(N-vinyl-2-pyrrolidone) (PVP). The atomic structure of the PVP-stabilized Au24Cl x was investigated by means of aberration-corrected transmission electron microscopy (ACTEM) and density functional theory (DFT) calculations. ACTEM video imaging revealed that the Au24Cl x clusters were stable against dissociation but fluctuated during the observation period. Some of the high-resolution ACTEM snapshots were explained by DFT-optimized isomeric structures in which all the constituent atoms were located on the surface. This observation suggests that the featureless optical spectrum of Au24Cl x is associated with the coexistence of distinctive isomers. X-ray photoelectron spectroscopy and Fourier-transform infrared spectroscopy of CO adsorbates revealed the electron-rich nature of Au24Cl x clusters due to the interaction with PVP. The Au24Cl x :PVP clusters catalyzed the aerobic oxidation of benzyl alcohol derivatives without degradation. Hammett analysis and the kinetic isotope effect indicated that the hydride elimination by Au24Cl x was the rate-limiting step with an apparent activation energy of 56 ± 3 kJ/mol, whereas the oxygen pressure dependence of the reaction kinetics suggested the involvement of hydrogen abstraction by coadsorbed O2 as a faster process.

Development of a Structural Comparison Method to Promote Exploration of the Potential Energy Surface in the Global Optimization of Nanoclusters

A structural comparison method (SCM) was created to quantify the structural diversity of nanoclusters and was implemented into a global optimization algorithm to evaluate structural diversity between generated clusters on the fly and promote exploration of the potential energy surface. The SCM evaluated topological differences between clusters using the common neighbor analysis and provided a numerical measure of similarity between the two clusters. The SCM was implemented into a genetic algorithm by integrating it into a new structure + energy fitness operator such that structural diversity of clusters in the population and their energies were used to assign fitness values to clusters. The efficiency of the genetic algorithm with this new fitness operator was benchmarked against several Lennard-Jones clusters (LJ₃₈, LJ₇₅, and LJ₉₈) known to be difficult cases for global optimization algorithms. For LJ₃₈ and LJ₇₅, this new structure + energy fitness operator performed equally well or better than the energy fitness operator. However, the efficiency of locating the global minimum of LJ₉₈ decreased using this new structure + energy fitness operator. Further analysis of the genetic algorithm with this fitness operator showed that the algorithm did indeed promote exploration of the PES of LJ₉₈ as desired but hindered refinement of clusters, preventing it from locating the global minimum even if the energy funnel of the global minimum had been located.

Born to be different: the formation process of Cu nanoparticles tunes the size trend of the activity for CO2 to CH4 conversion

We investigate the impact of the formation process of Cu nanoparticles on the distribution of adsorption sites and hence on their activity. Using molecular dynamics, we model formation pathways characteristic of physical synthesis routes as the annealing of a liquid droplet, the growth proceeding via the addition of single atoms, and the coalescence of individual nanoparticles. Each formation process leads to different and characteristic size-dependent distributions of their adsorption sites, catalogued and monitored on-the-fly by means of a suitable geometrical descriptor. Annealed or coalesced nanoparticles present a rather homogeneous distribution in the kind and relative abundance of non-equivalent adsorption sites. Atom-by-atom grown nanoparticles, instead, exhibit a more marked occurrence of adsorption sites corresponding to adatoms and small islands on (111) and (100) facets. Regardless of the formation process, highly coordinated sites are more likely in larger nanoparticles, while the abundance of low-coordination sites depends on the formation process and on the nanoparticle size. Furthermore, we show how each characteristic distribution of adsorption sites reflects in different size trends for the Cu-nanoparticle activity, taking as an example the electro-reduction of CO₂ into CH₄. To this end, we employ a multi-scale method and observe that the faceted but highly defected structures obtained during the atom-by-atom growth become more and more active with increasing size, with a mild dependence on the original seed. In contrast, the activity of Cu-nanoparticles obtained by annealing decreases with their size, while coalesced nanoparticles' activity shows a non-monotonic behaviour.

Exploring the Structure–Activity Relationship on Platinum Nanoparticles

The design of active and stable Pt-based nanoscale electrocatalysts for the oxygen reduction reaction (ORR) plays the central role in ameliorating the efficiency of proton exchange membrane fuel-cells towards future energy applications. On that front, theoretical studies have contributed significantly to this research area by gaining deeper insights and understanding of the ongoing processes. In this work, we present an approach capable of characterizing differently-shaped platinum nanoparticles undergoing thermally- and adsorbate-induced restructuring of the surface. Further, by performing ReaxFF-Grand Canonical Molecular Dynamics simulations we explored the water formation on these roughened (“realistic”) nanoparticles in a H2/O2 environment. Taking into consideration the coverage of oxygen-containing intermediates and occurring surface roughening the nanoparticles’ activities were explored. Hereby, we succeeded in locally resolving the water formation on the nanoparticles’ surfaces, allowing an allocation of the active sites for H2O production. We observed that exposed, low-coordinated sites as well as pit-shaped sites originating from roughening of vertices and edges are most active towards H2O formation.

Simulations of the Oxidation and Degradation of Platinum Electrocatalysts

Improved understanding of the fundamental processes leading to degradation of platinum nanoparticle electrocatalysts is essential to the continued advancement of their catalytic activity and stability. To this end, the oxidation of platinum nanoparticles is simulated using a ReaxFF reactive force field within a grand-canonical Monte Carlo scheme. 2-4 nm cuboctahedral particles serve as model systems, for which electrochemical potential-dependent phase diagrams are constructed from the thermodynamically most stable oxide structures, including solvation and thermochemical contributions. Calculations in this study suggest that surface oxide structures should become thermodynamically stable at voltages around 0.80-0.85 V versus standard hydrogen electrode, which corresponds to typical fuel cell operating conditions. The potential presence of a surface oxide during catalysis is usually not accounted for in theoretical studies of Pt electrocatalysts. Beyond 1.1 V, fragmentation of the catalyst particles into [Pt₆ O₈ ]^4- clusters is observed. Density functional theory calculations confirm that [Pt₆ O₈ ]^4- is indeed stable and hydrophilic. These results suggest that the formation of [Pt₆ O₈ ]^4- may play an important role in platinum catalyst degradation as well as the electromotoric transport of Pt^2+/4+ ions in fuel cells.