Structural and electronic properties of small copper clusters: a first principles study

Molecular Adsorption of H2 on Small Neutral Silver-Copper Bimetallic Nanoparticles: A Search for Novel Hydrogen Storage Materials

In the search for novel hydrogen storage materials, neutral silver-copper bimetallic nanoparticles up to the size of eight atoms (Cu _m Ag _n : m + n ≤ 8) have been computationally studied. Density functional theory with the B3LYP exchange-correlation functional and the combined basis sets of LanL2DZ and aug-cc-pVQZ were used in all of the calculations. H₂ adsorption studies on the most stable cluster geometries of all of the neat and heterogeneous entities found that 12 potential candidates, CuAg₄, Cu₆, Cu₅Ag, Cu₄Ag₂, Cu₃Ag₃, Cu₂Ag₄, CuAg₆, Cu₅Ag₃, Cu₄Ag₄, Cu₃Ag₅, Cu₂Ag₆, and CuAg₇, fall within the recommended physisorption range of -18 to -6 kJ mol^-1. A correlation in the behavior of binding energy, vibrational frequency, average bond distance, highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap, and chemical hardness with H₂ adsorption was observed. This analysis further revealed that the H₂ adsorption to the cluster was either a parallel or a perpendicular alignment. The analysis of the electron configuration of each atom in the cluster and the H₂ molecule and the charge transfer analysis of these 12 clusters also showed that the physisorption in the perpendicular mechanism is due to an induced dipole interaction, while that in the parallel mechanism is due to a weak ionic interaction. The clusters identified with perpendicular adsorption, CuAg₄H₂, Cu₆H₂, Cu₃Ag₃H₂, and Cu₂Ag₄H₂, polarized the H₂ molecule but had no charge transfer with the H₂ molecule and those identified with parallel adsorption, Cu₅AgH₂, Cu₄Ag₂H₂, CuAg₆H₂, Cu₅Ag₃H₂, Cu₄Ag₄H₂, Cu₃Ag₅H₂, Cu₂Ag₆H₂, and CuAg₇H₂, pulled the electrons from the H₂ molecule and had charge transfer with the H₂ molecule. The shapes of the frontier molecular orbital diagrams of the HOMO and LUMO also followed this observation.

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.

A density functional study of Rh13

A density functional study was performed for the Rh13 cluster using the linear combination of Gaussian-type orbitals density functional theory (LCGTO-DFT) approach. The calculations employed both the local density approximation (LDA) as well as the generalized gradient approximation (GGA) in combination with a quasi-relativistic effective core potential (QECP). Initial structures for the geometry optimization were taken along Born–Oppenheimer molecular dynamics (BOMD) trajectories. The BOMD trajectories were performed at different temperatures and considered different potential energy surfaces (PES). As a result, several hundred isomers of the Rh13 cluster in different spin multiplicities were optimized with the aim to determine the lowest energy structures. All geometry optimizations were performed without any symmetry restriction. A vibrational analysis was performed to characterize these isomers. Structural parameters, relative stability energy, harmonic frequencies, binding energy, and most relevant Kohn–Sham (KS) molecular orbitals are reported. The obtained results are compared with available data from the literature. This study predicts a low symmetry biplanarlike structure as the ground-state structure of Rh13 with 11 unpaired electrons. This isomer was first noticed by inspection of first-principle Born–Oppenheimer molecular dynamics (BOMD) simulations between 300 and 600 K. This represents the most extensive theoretical study on the ground-state structure of the Rh13 cluster and underlines the importance of BOMD simulations to fully explore the PES landscapes of complicated systems.

Assessment of density functional theory optimized basis sets for gradient corrected functionals to transition metal systems: The case of small Nin (n⩽5) clusters

Density functional calculations have been performed for small nickel clusters, Ni(n), Ni(n) (+), and Ni(n)(-) (n<or=5), using the linear combination of Gaussian-type orbital density functional theory approach. Newly developed nickel all-electron basis sets optimized for generalized gradient approximation (GGA) as well as an all-electron basis set optimized for the local density approximation were employed. For both neutral and charged systems, several isomers and different multiplicities were studied in order to determine the lowest energy structures. A vibrational analysis was performed in order to characterize these isomers. Structural parameters, harmonic frequencies, binding energies, ionization potentials, and electron affinities are reported. This work shows that the employed GGA basis sets for the nickel atom are important for the correct prediction of the ground state structures of small nickel clusters and that the structural assignment of these systems can be performed, with a good resolution, over the ionization potential.

Structure and shape variations in intermediate-size copper clusters

Using extensive, unbiased searches based on density-functional theory, we explore the structural evolution of Cu(n) clusters over the size range n=8-20. For n=8-16, the optimal structures are plateletlike, consisting of two layers, with the atoms in each layer forming a trigonal bonding network similar to that found in smaller, planar clusters (n<or=6). For n=17 and beyond, there is a transition to compact structures containing an icosahedral 13-atom core. The calculated ground-state structures are significantly different from those predicted earlier in studies based on empirical and semiempirical potentials. The evolution of the structure and shape of the preferred configuration of Cu(n), n<or=20, is shown to be nearly identical to that found for Na clusters, indicating a shell-model-type behavior in this size range.

Homonuclear transition-metal trimers

Density-functional theory has been used to determine the ground-state geometries and electronic states for homonuclear transition-metal trimers constrained to equilateral triangle geometries. This represents the first application of consistent theoretical methods to all of the ten 3d block transition-metal trimers, from scandium to zinc. A search of the potential surfaces yields the following electronic ground states and bond lengths: Sc3(2A1',2.83 A), Ti3(7E',2.32 A), V3(2E",2.06 A), Cr3(17E',2.92 A), Mn3(16A2',2.73 A), Fe3(11E",2.24 A), Co3(6E",2.18 A), Ni3(3A2",2.23 A), Cu3(2E',2.37 A), and Zn3(1A1',2.93 A). Vibrational frequencies, several low-lying electronic states, and trends in bond lengths and atomization energies are discussed. The predicted dissociation energies DeltaE(M3-->M2+M) are 49.4 kcal mol(-1)(Sc3), 64.3 kcal mol(-1)(Ti3), 60.7 kcal mol(-1)(V3), 11.5 kcal mol(-1)(Cr3), 32.4 kcal mol(-1)(Mn3), 61.5 kcal mol(-1)(Fe3), 78.0 kcal mol(-1)(Co3), 86.1 kcal mol(-1)(Ni3), 26.8 kcal mol(-1)(Cu3), and 4.5 kcal mol(-1)(Zn3).

High accuracy ab initio studies of Li6+, Li6−, and three isomers of Li6

The structures and energetics of Li(6) (+), Li(6) (-) and three isomers of Li(6) are investigated using the coupled-cluster singles, doubles and perturbative triples [CCSD(T)] method with valence and core-valence correlation consistent basis sets of double- to quadruple-zeta quality (cc-pVXZ and cc-pCVXZ, where X=D-Q). These results are compared with qualitatively different predictions by less reliable methods. Our results conclusively show that the D(4h) isomer is the global minimum structure for Li(6). It is energetically favored over the C(5v) and D(3h) structures by about 5.1 and 7.1 kcal mol(-1), respectively, after the inclusion of the zero-point vibrational energy (ZPVE) correction. Our most accurate total atomization energies are 123.2, 117.6, and 115.7 kcal mol(-1) for the D(4h), C(5v), and D(3h) isomers, respectively. Comparison of experimental optical absorption spectra with our computed electronic spectra also indicate that the D(4h) isomer is indeed the most stable structure. The cation, anion, and some higher spin states are investigated using the less expensive cc-pCVDZ basis set. Adiabatic ionization energies and electron affinities are reported and compared with experimental values. Predictions of molecular properties are found to be sensitive to the basis set used and to the treatment of electron correlation.

Electric dipole polarizabilities of copper clusters

The static electric dipole polarizabilities of Cu(9)-Cu(61) have been measured via a molecular beam deflection method. The clusters display per-atom polarizabilities that decrease monotonically with size, from approximately 16 A(3) per atom Cu(9-10) to approximately 5 A(3) (Cu(45-61)). Absent are any discernible discontinuities or odd-even alternations due to electronic shell filling or electron pairing effects. For the smallest clusters, the experimental polarizabilities are approximately 3 times larger than those predicted classically for conducting ellipsoids, and approach the classical values only for clusters containing more than approximately 45 atoms.