Structure and shape variations in intermediate-size copper clusters

Theoretical study of NH3 , H2 S, and HCN adsorption enhancement on defective graphene-supported Cu19 clusters

Recent studies have shown that graphene-supported metal clusters can enhance catalytic reactivity compared with corresponding metal clusters. In this study, the adsorptions of NH3 , H2 S, and HCN on Cu19 and defective graphene-supported Cu19 clusters are investigated using plane-wave density functional theory. The results reveal the three gas molecules can be adsorbed on three types of top sites of Cu atoms, respectively. The adsorption energies of the corresponding adsorption sites on the defective graphene-supported Cu19 clusters are all increased compared with those on the Cu19 clusters. The orbital-resolved, crystal orbital Hamilton population analysis demonstrates that the larger the integrated crystal orbital Hamilton population, the stronger the adsorption between the gas molecule and the bonded Cu atom. The center of antibonding states on the defective graphene-supported Cu19 is shifted upward relative to Fermi level compared to the corresponding one on pure Cu19 , which explains the enhanced adsorption energy on defective graphene-supported Cu19. In addition, the closer d-band center to the Fermi level on the defective graphene-supported Cu19 indicates a stronger adsorption capacity than on pure Cu19 .

Enhanced CO adsorption on α-graphyne-supported and defective graphene-supported Cu19 clusters and a modified induction energy model

Ab initio molecular dynamics calculations were carried out to study the adsorption of CO on Cu19, α-graphyne-supported Cu19 and defective graphene-supported Cu19 clusters. The average adsorption energies on the three clusters are significantly increased by 68%, 104%, and 123% compared to the experimental value on the pristine Cu(110) surface. Furthermore, the α-graphyne-supported and defective graphene-supported Cu19 clusters exhibit greater adsorption strength than the pure Cu19 cluster, with 22% and 33% higher adsorption energies, respectively. The crystal orbital Hamilton population analysis shows that for the same type of adsorption site, the adsorption energy is linearly related to the bond interaction strength between the adsorbate and the substrate. We propose a modified induction energy model to predict the increase of chemisorption energy on α-graphyne-supported and defective graphene-supported Cu19 clusters based on the bare Cu19 cluster. The chemisorption energy enhancement predicted by the improved induction energy model has very good agreement with that calculated based on the ab initio molecular dynamics method and is more accurate than that predicted by the original induction energy model.

Low-energy configurations of Pt6Cu6 clusters and their physical-chemical characterization: a high-accuracy DFT study

Based on a combination of many-body potentials, an analysis of the inertia tensors and a Density Functional Theory framework, we use a method to harvest the lowest energy states of any set of cluster systems. Then, this methodology is applied to the Pt₆Cu₆ cluster case and the structural, chemical, electronic, anisotropy, magnetic and vibrational properties of the lowest energy isomers are studied. Unexpectedly, some tens of isomers with much lower energy than the precedent believed ground state [J. Chem. Phys., 131(4):044701] are found, which indicates the goodness of this methodology. Some of the isomers obtained present the point groups C_s, C_2v according to Schoenflies notation, while others do not exhibit specific symmetry operations. The global chemical descriptors as the ionization potential, the electron affinity and the chemical hardness have oscillating behaviors with overall decreasing trends as the energy of the isomer grows up, indicating a higher rate of deactivation by sintering processes and a higher strength of the adsorption of small molecules on these systems. We present interesting results of the electronic, magnetic, anisotropy, vibrational and thermal properties of these clusters and discuss them; what can be useful information for future experiments and technical applications in varied fields as catalysis, spintronics, molecular magnetism or magnetic storage information.

Computational Modeling of the Interaction of Silver Clusters with Carbohydrates

Silver nanoparticles are recognized for their numerous physical, biological, and pharmaceutical applications. In the present study, the interaction of silver clusters with monosaccharide molecules is examined to identify which molecule works better as a reducing agent in the application of a green synthesis approach. Geometry optimization of clusters containing one, three, and five silver atoms is performed along with the optimization of α-d-glucose, α-d-ribose, d-erythrose, and glyceraldehyde using density functional theory. Optimized geometries allow identifying the interaction formed in the silver cluster and monosaccharide complexes. An electron localization function analysis is performed to further analyze the interaction found and explain the reduction process in the formation of silver nanoparticles. The overall results indicate that glyceraldehyde presents the best characteristics to serve as the most efficient reducing agent.

Effects of Temperature on Enantiomerization Energy and Distribution of Isomers in the Chiral Cu13 Cluster

In this study, we report the lowest energy structure of bare Cu₁₃ nanoclusters as a pair of enantiomers at room temperature. Moreover, we compute the enantiomerization energy for the interconversion from minus to plus structures in the chiral putative global minimum for temperatures ranging from 20 to 1300 K. Additionally, employing nanothermodynamics, we compute the probabilities of occurrence for each particular isomer as a function of temperature. To achieve that, we explore the free energy surface of the Cu₁₃ cluster, employing a genetic algorithm coupled with density functional theory. Moreover, we discuss the energetic ordering of isomers computed with various density functionals. Based on the computed thermal population, our results show that the chiral putative global minimum strongly dominates at room temperature.

Systematic cluster growth: a structure search method for transition metal clusters

The systematic cluster growth (SCG) method is a biased structure search strategy based on a seeding process for investigating the structural evolution and growth pattern of transition metal clusters.

Reliability of semiempirical and DFTB methods for the global optimization of the structures of nanoclusters

In this work, we explore the possibility of using computationally inexpensive electronic structure methods, such as semiempirical and DFTB calculations, for the search of the global minimum (GM) structure of chemical systems. The basic prerequisite that these inexpensive methods will need to fulfill is that their lowest energy structures can be used as starting point for a subsequent local optimization at a benchmark level that will yield its GM. If this is possible, one could bypass the global optimization at the expensive method, which is currently impossible except for very small molecules. Specifically, we test our methods with clusters of second row elements including systems of several bonding types, such as alkali, metal, and covalent clusters. The results reveal that the DFTB3 method yields reasonable results and is a potential candidate for this type of applications. Even though the DFTB2 approach using standard parameters is proven to yield poor results, we show that a re-parametrization of only its repulsive part is enough to achieve excellent results, even when applied to larger systems outside the training set.

Adsorption-Induced Liquid-to-Solid Phase Transition of Cu Clusters in Catalytic Dissociation of CO2

Sub-nanometer metal clusters widely existing in catalysts have a large ensemble of metastable isomers that can interconvert during catalytic reactions, exhibiting complex dynamical catalytic effects. In this work, we systematically investigate the temperature dependent structural dynamics of the Cu₁₃ cluster in CO₂ dissociation using ab initio molecular dynamics and the free energy calculation method. We find an abnormal entropic effect due to adsorption-induced liquid-to-solid phase transition of the cluster during the course of the elementary dissociation step at transition temperatures. In the dissociation product, the formation of a rigid Cu₃O unit decreases the dynamical fluidity of the cluster and increases the melting temperature, causing a decrease in the entropy of the dissociation product. Our work demonstrates the nontrivial effects of surface adsorption on phase transition behaviors of dynamic clusters and offers a new perspective to dynamic catalysis.

A method for predicting basins in the global optimization of nanoclusters with applications to AlxCuy alloys

The problem of obtaining the geometrical configuration of a molecule that minimizes its potential energy is a very complicated one for a series of applications, ranging from determining the structure of biological macromolecules to nanoclusters of atoms. Global optimization tools are available for this task, and many of them are based in performing successive local optimizations, where the starting geometries for these steps are determined by an intelligent algorithm. Here we develop a method to save computing time in the optimization of nanoclusters by predicting if a given minimum has been previously visited during local optimization steps. Our application to Cu-Al nanoalloys indicates that it is possible to save a substantial amount of computational cost. The application also reveals new promising Al_xCu_y clusters and explain their stabilities in terms of the jellium model.

Anionic Copper Clusters Reacting with NO: An Open-Shell Superatom Cu18. J Phys Chem Lett 2020;11:5807-5814. [PMID: 32597656 DOI: 10.1021/acs.jpclett.0c01643]

Gas-phase metal clusters have been a subject of research interest for allowing reliable strategies to explore the stability and reactivity of materials at reduced sizes with atomic precision. Here we have prepared well-resolved copper cluster anions Cu_n^- (n = 7-37) and systematically studied their reactivity with O₂, NO, and CO. We found remarkable stability of an open-shell cluster Cu₁₈^-, which is comparable with the closed-shell clusters Cu₁₇^- and Cu₁₉^- within the picture of an electronic shell model. Even without having a magic number of valence electrons, intriguingly, the unpaired electron on the singly occupied molecular orbital of Cu₁₈^- is mainly contributed by the central copper atom, while the other 18 delocalized valence electrons occupy the lower-energy superatomic orbitals of the cluster. The finding of such an open-shell superatom Cu₁₈^-, with an electron configuration of 1S²1P⁶1D¹⁰2S¹||1F⁰, is interesting in the sense that an elementary cluster of coinage metal atoms could still behave as a superatom mimicking coinage metals like silver or gold atoms with an empty f orbital. The superatomic stability of this Cu₁₈^- cluster is reinforced by the unique electrostatic interaction between the Cu^- core and Cu₁₇ shell, which provides new insights into the chemistry of metal clusters.

Effect of atomicity on the oxidation of cationic copper clusters studied using thermal desorption spectrometry

The resistivity to oxidation of small copper clusters, Cun+ (n ≤ 5), in the gas phase with a precise atomicity at the molecular level was investigated using a combination of thermal desorption spectrometry and mass spectrometry. Oxide clusters, CunOm+, with more O atoms than those present with a stoichiometry of n : m = 1 : 1 were produced at room temperature in the presence of O2, and the weakly bound excess oxygen atoms involved in the clusters were removed by post heating. Non-oxidized Cu2+ and Cu3+ clusters were formed in the range of 323-923 K, whereas partially oxidized clusters, Cu4O2+ and Cu5O2+, were generated for n = 4 and 5. Considering the fact that CunOm+ (m = n/2 + 1) tends to be generated for n ≥ 6, the small copper clusters were concluded to be resistive to oxidation. The possible reaction paths for the oxidation of Cu2+ and Cu4+ clusters were obtained by density functional calculations, which were consistent with the experimental findings. The oxidation states of the Cu atoms in the clusters were discussed based on the natural charges of the atoms.

The stability of Cu clusters and their adsorption for CH4 and CH3 by first principle calculations

Two-dimensional (2D) and three-dimensional (3D) Cu _n clusters (n indicates the atom number) and their adsorption behaviors for both methane (CH₄) and methyl (CH₃) are studied in this work using the density functional theory method, where n ranges from 6 to 20. In these small clusters, it is found that the CH₄ molecule is always adsorbed on the top site with the adsorption energy between -0.05 eV and -0.21 eV. Considering methane dehydrogenation, stronger adsorption for CH₄ is required, so 2D clusters with n = 7, 14, 15, and 16 and 3D clusters with n = 6, 10, 12, and 17 are found to have relatively stronger adsorption. However, for the adsorption of CH₃, there is an obvious even-odd oscillation change in the size of 3D clusters, while it is not clear in 2D clusters since one cannot find an even-odd change as n > 14. The weaker adsorption for CH₃ occurs on 3D clusters when n is even except 6 and also on 2D clusters when n = 6, 7, 10, and 12 with higher carbon poisoning resistance. Based on these calculated results, some Cu clusters which show good potential ability for methane dehydrogenation are provided, especially when n = 10 and 12 for 3D structures, and n = 7 for the 2D ones.

Uniaxial magnetic anisotropy energy of bimetallic Co-Ni clusters from a first-principles perspective

Along with the growing precision in the control of matter at increasingly smaller size scales, a field of research, based onto magnetic materials of technical interest, such as bimetallic clusters, has been developed in very recent years. Thereby, here, we report on a complete study of bimetallic clusters composed of cobalt and nickel with up to 7 atoms using ab initio methods in the GGA approach. We applied an unrestricted search method based on the tensor of inertia eigenvalues to find the most stable configurations of the clusters, obtaining a diverse set of structures with different geometric properties. We explored the effect of composition on the structural properties, the chemical stability, the magnetization and the magnetic anisotropy energy (MAE) of the so-obtained systems. Our results indicate that the behavior of the clusters is mainly governed by the Co-Co interaction and to a lesser extent by the Co-Ni and Ni-Ni interactions. Furthermore, for a given cluster size the magnetic moment increases by 2 μB/Co-substitution plus 1 μB/Ni-substitution coming from the cobalt and nickel core d-states, while in some cases unpaired hybrid s-electrons can also give rise to itinerant magnetism. These features have been analyzed with the help of a Jellium model and have important consequences for the magnetism and the magnetic anisotropy of the clusters. The magnetic behavior and MAE present complex and intriguing landscapes, which suggests the possibility of finely controlling the magnetic states, by tuning the cluster composition, aiming at technical implementation in fields such as molecular magnetism or quantum computation. In particular, cases such as Co6Ni, Co4Ni3, CoNi6 Co2Ni5 and Co3Ni3 present high relative stability and enhanced magnetic moments (around 10 μB), what makes them promising candidates for applications such as subnanometer magnetic information storage.

Evolution of the structural, energetic, and electronic properties of the 3d, 4d, and 5d transition-metal clusters (30 TMn systems for n = 2–15): a density functional theory investigation

Subnanometric transition-metal (TM) clusters have attracted great attention due to their unexpected physical and chemical properties, leastwise compared to their bulk counterparts.

DFT simulations and microkinetic modelling of 1-pentyne hydrogenation on Cu20 model catalysts

Adsorption and dissociation of H2 and hydrogenation of 1-pentyne on neutral and anionic Cu20 clusters have been investigated using the density functional theory and microkinetic modelling. Molecular adsorption of H2 is found to occur strictly at atop sites. The H2 dimer is activated upon adsorption, and the dissociation occurs with moderate energy barriers. The dissociated H atoms reside preferentially on 3-fold face and 2-fold edge sites. Based on these results, the reaction paths leading to the partial and total hydrogenation of 1-pentyne have been studied step-by-step. The results suggest that copper clusters can display selective activity on the hydrogenation of alkyne and alkene molecules. The hydrogenated products are more stable than the corresponding initial reactants following an energetic staircase with the number of added H atoms. Stable semi-hydrogenated intermediates are formed before the partial (1-pentene) and total (pentane) hydrogenation stages of 1-pentyne. The microkinetic model analysis shows that C5H10 is the dominant product. Increasing the reactants (C5H8/H2) ratio enhances the formation of products (C5H10 and C5H12).

CO oxidation catalyzed by neutral and anionic Cu20 clusters: relationship between charge and activity

Reactions of CO and O2 on neutral and anionic Cu20 clusters have been investigated by spin-polarized density functional theory. Three reaction mechanisms of CO oxidation are explored: reactions with atomic oxygen (dissociated O2) as well as reactions with molecular oxygen, including Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) mechanisms. The adsorption energies, reaction pathways, and reaction barriers for CO oxidation are calculated systematically. The anionic Cu20(-) cluster can adsorb CO and O2 more strongly than the neutral counterpart due to the superatomic shell closing of 20 valence electrons which leaves one electron above the band gap. The activation of O2 molecule upon adsorption is crucial to determine the rate of CO oxidation. The CO oxidation proceeds efficiently on both Cu20 and Cu20(-) clusters, when O2 is pre-adsorbed dissociatively. The ER mechanism has a lower reaction barrier than the LH mechanism on the neutral Cu20. In general, CO oxidation occurs more readily on the anionic Cu20(-) (effective reaction barriers 0.1-0.3 eV) than on the neutral Cu20 cluster (0.3-0.5 eV). Moreover, Cu20(-) exhibits enhanced binding for CO2. From the analysis of the reverse direction of CO oxidation, it is observed that the transition of CO2 to CO + O can occur on the Cu20(-) cluster, which demonstrates that Cu clusters may serve as good catalyst for CO2 chemistry.

Nanoclusters synthesized by synchrotron radiolysis in concert with wet chemistry

Wet chemical reduction of metal ions, a common strategy for synthesizing metal nanoparticles, strongly depends on the electric potential of the metal, and its applications to late transition metal clusters have been limited to special cases. Here, we describe copper nanoclusters grown by synchrotron radiolysis in concert with wet chemistry. The local structure of copper aggregates grown by reducing Cu(II) pentanedionate using synchrotron x-ray beam was studied in situ by x-ray absorption spectroscopy. A detailed analysis of the XANES and EXAFS spectra, compared with DFT calculations and full-potential non-muffin-tin multiple scattering calculations, identified the nanocluster as Cu13 with icosahedral symmetry. The novel "charged" nanoclusters tightly bound to electron-donating amido molecules, which formed as a result of photo-induced deprotonation of ligand amines, were stabilized by irradiation. Monodispersive deposition of nanoclusters was enabled by controlling the type and density of "monomers", in remarkable contrast to the conventional growth of metallic nanoparticles.

Size-selective catalytic growth of nearly 100% pure carbon nanocoils with copper nanoparticles produced by atomic layer deposition

In this paper, Cu nanoparticles with narrow size distribution are synthesized by reduction of CuO films produced by atomic layer deposition (ALD), which are used as catalysts for the catalytic growth of carbon nanostructures. By properly adjusting the ALD cycle numbers, the size of produced Cu nanoparticles can be well controlled. Uniform carbon nanocoils with near 100% purity can be obtained by using 50-80 nm Cu nanoparticles, while thin straight fibers and thick straight fibers are produced by applying 5-35 and 100-200 nm Cu nanoparticles, respectively. The mechanism of the particle size-dependent growth of the carbon nanostructure was analyzed based on the experimental results and theoretical simulation. Our results can provide important information for the preparation of helical carbon nanostructures with high purity. Moreover, this work also demonstrates that ALD is a viable technique for synthesizing nanoparticles with highly controllable size and narrow size distribution suitable for studying particle size-dependent catalytic behavior and other applications.

Structural and electronic properties of tungsten nanoclusters by DFT and basin-hopping calculations

The structural and electronic properties of small tungsten nanoclusters W_n (n = 2–16) were investigated by density functional theory (DFT) calculations.

Gas-induced formation of Cu nanoparticle as catalyst for high-purity straight and helical carbon nanofibers

The facile preparation of high-purity carbon nanofibers (CNFs) remains challenging due to the high complexity and low controllability in reaction. A novel approach using gas-induced formation of Cu crystals to control the growth of CNFs is developed in this study. By adjusting the atmospheric composition, controllable preparation of Cu nanoparticles (NPs) with specific size and shape is achieved, and they are further used as a catalyst for the growth of straight or helical CNFs with good selectivity and high yield. The preparation of Cu NPs and the formation of CNFs are completed by a one-step process. The inducing effect of N(2), Ar, H(2), and C(2)H(2) on the formation of Cu NPs is systematically investigated through a combined experimental and computational approach. The morphology of CNFs obtained under different conditions is rationalized in terms of Cu NP and CNF growth models. The results suggest that the shapes of CNFs, namely, straight or helical, depend closely on the size, shape, and facet activity of Cu NPs, while such a gas-inducing method offers a simple way to control the formation of Cu NPs.

Icosahedral to double-icosahedral shape transition of copper clusters

The lowest-energy isomers of Cu(N) clusters for N = 20-30 are identified using an unbiased search algorithm and density functional theory calculations. The low-energy structures over this size range are dominated by those based on a 13-atom icosahedral (I(h)) core and a 19-atom double icosahedron (DI(h)) core. A transition in the ground-state isomers from I(h)-based to DI(h)-based structures is predicted overt N = 21-23. We discuss this transition in the broader context of the growth pattern for Cu(N) over N = 2-30 that features regions of gradual evolution in which atoms successively add to the cluster surface, separated by sudden changes to a different structural organization and more compact shape. These transitions result from a competition between interatomic bonding energy and surface energy. The implications of this growth pattern for the further evolution of copper from microstructure to bulk are discussed.

Density functional investigations on structural and electronic properties of anionic and neutral sodium clusters Na(N) (N = 40-147): comparison with the experimental photoelectron spectra

Density functional calculations have been carried out to obtain low energy equilibrium geometries of anionic and neutral sodium clusters over a wide range of sizes 40 ≤ N ≤ 147, where N is the number of atoms. An exhaustive search for the low energy equilibrium geometries has been carried out. The density of states of the lowest energy geometries are compared with the experimental photoelectron spectra (Huber et al 2009 Phys. Rev. B 80 235425; Kostko et al 2007 Phys. Rev. Lett.98 043401) for N > 41. The agreement between theory and experiment is good for almost all the clusters and the changes in the spectrum with size correlate very well with the changes in the shapes as observed in the evolutionary trend of the ground state geometries.

Growth mechanism and chemical bonding in scandium-doped copper clusters: experimental and theoretical study in concert

Size matters! The electronic structure and size-dependent stability of neutral and cationic scandium-doped copper clusters have been investigated by mass spectrometric studies (for the cations) and also quantum chemical computations. The proposed reaction paths ultimately lead to the most stable Frank-Kasper-shaped Cu(16)Sc(+) cluster (shown here), which could be the germ of a new crystallization process.Electronic structure and size-dependent stability of scandium-doped copper cluster cations, Cu(n)Sc(+), were investigated by using a dual-target dual-laser vaporization production scheme followed by mass spectrometric studies and also quantum chemical computations in the density functional theory framework. The neutral species also were studied by using computational methods. Enhanced abundances and dissociation energies were measured in the case of Cu(n)Sc(+) for n=4, 6, 8, 10 and 16, the last of these identified as being extraordinary stable. Neutral clusters are stable with n=5, 7, 9 and 15, which are isoelectronic with respect to the number of the valence s electrons with the stable cationic clusters; hence a simple electron count determines cluster properties to a great extent. The Cu(17)Sc cluster was found to be a superatomic molecule, containing Cu(16)Sc(+) and Cu(-) units; however, the charge separation is not as pronounced as in the case of CuLi. Cu(15)Sc was found to be a stable cluster with a large dissociation energy and a closed electronic structure; hence this can be regarded as a superatom, analogous to the noble gases. The main factors determining the growth patterns of these clusters are the central position of the scandium atom and the successive filling of the shell orbitals. For smaller clusters, the reaction paths appear to diverge yielding various products; however all paths ultimately lead to the most stable Frank-Kasper shaped Cu(16)Sc cluster, which in turn can be the germ of the crystallization process.

Diversity driven unbiased search of minimum energy cluster configurations

The determination of the spatial distributions that atoms adopt to form condensed matter is a problem of crucial importance, since most physical properties depend on the atomic arrangement. This is especially relevant for clusters, where periodicity is nonexistent. Several optimization procedures have been implemented to tackle this problem, with ever increasing success. Here we put forward a search scheme which preserves as large a diversity as allowed by the use of phenomenological potentials, generating in an unbiased fashion a bank of configurations to be explored; a procedure we denominate diversity driven unbiased search (DDUS). It consists in the generation, using phenomenological potentials, of a data bank of putative minima rather than a single, or just a few, configurations which are based on the conformational space annealing method (CSA). All of the configurations in the bank are thereafter refined by means of DFT computations. Certainly, in spite of our efforts to generate a bank as diverse as possible, not all relevant structures might be included in it, since quantum effects are ignored. The procedure is applied to several examples of rhodium, palladium, silver, platinum and gold clusters, between 5 and 23 atoms in size. The main conclusion we reach is that unbiased search, among a significant number of candidates, quite often leads to rather unexpectedly low symmetry configurations, which turn out to be the lowest energy ones within our scheme.