Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

Surface collective oscillations of metal clusters and spheres: Random-phase-approximation sum-rules approach

Using the once and thrice energy-weighted moments of the random-phase-approximation strength function, we have derived compact expressions for the average energy of surface collective oscillations of clusters and spheres of metal atoms. The L=0 volume mode has also been studied. We have carried out quantal and semiclassical calculations for Na and Ag systems in the spherical-jellium approximation. We present a rather thorough discussion of surface diffuseness and quantal size effects on the resonance energies.

Static dipole polarizability of alkali-metal clusters: Electronic exchange and correlation effects

Nonlocal approximations for the electronic exchange and correlation effects are used to compute, within density-functional theory, the polarizability and surface-plasma frequencies of small jelliumlike alkali-metal clusters. The results are compared with those obtained using the local-density approximation and with available experimental data, showing the relevance of these effects in obtaining an accurate description of the surface response of metallic clusters.

Barrier for the reaction X20+ + X20+ -> X402+ in alkali-metal clusters related to electron density at the bond midpoint of the supermolecule (X20+)2

Using the extended Thomas-Fermi version of density-functional theory (DFT), calculations are presented for the barrier for the reaction Na20++Na20+¿Na402+. The deviation from the simple Coulomb barrier is shown to be proportional to the electron density at the bond midpoint of the supermolecule (Na20+)2. An extension of conventional quantum-chemical studies of homonuclear diatomic molecular ions is then effected to apply to the supermolecular ions of the alkali metals. This then allows the Na results to be utilized to make semiquantitative predictions of position and height of the maximum of the fusion barrier for other alkali clusters. These predictions are confirmed by means of similar DFT calculations for the K clusters.

Deformed-jellium model for the fission of multiply charged metal clusters

A deformed-jellium model is used to calculate the fission barrier height of positive doubly charged sodium clusters within an extended Thomas-Fermi approximation. The fissioning cluster is continuously deformed from the parent configuration until it splits into two fragments. Although the shape of the fission barrier obviously depends on the parametrization of the fission path, we have found that remarkably, the maximum of the barrier corresponds to a configuration in which the emerging fragments are already formed and rather well apart. The implication of this finding in the calculation of critical numbers for fission is illustrated in the case of multiply charged Na clusters.

Collective spin excitations of alkali-metal clusters

The response function of alkali-metal clusters, modeled as jellium spheres, to dipole (L=1) and quadrupole (L=2) spin-dependent fields is obtained within the time-dependent local-spin-density approximation of density-functional theory. We predict the existence of low-energy spin modes of surface type, which are identified from the strength function. Their collectivity and evolution with size are discussed.

Essential factors for control of the equilibrium in the reversible rearrangement of M3N@Ih-C80 fulleropyrrolidines: Exohedral functional groups versus endohedral metal clusters

The effects of exohedral moieties and endohedral metal clusters on the isomerization of M3N@Ih-C80 products from the Prato reaction through [1,5]-sigmatropic rearrangement were systematically investigated by using three types of fulleropyrrolidine derivatives and four different endohedral metal clusters. As a result, all types of derivatives provided the same ratios of the isomers for a given trimetallic nitride template (TNT) as the thermodynamic products, thus indicating that the size of the endohedral metal clusters inside C80 was the single essential factor in determining the equilibrium between the [6,6]-isomer (kinetic product) and the [5,6]-isomer. In all the derivatives, the [6,6]- and [5,6]-Prato adducts with larger metal clusters, such as Y3N and Gd3N, were equally stable, which is in good agreement with DFT calculations. The reaction rate of the rearrangement was dependent on both the substituent of exohedral functional groups and the endohedral metal-cluster size. Further DFT calculations and 13C NMR spectroscopic studies were employed to rationalize the equilibrium in the rearrangement between the [6,6]- and [5,6]-fulleropyrrolidines

Deformed-jellium model for the fission of multiply charged metal clusters

A deformed-jellium model is used to calculate the fission barrier height of positive doubly charged sodium clusters within an extended Thomas-Fermi approximation. The fissioning cluster is continuously deformed from the parent configuration until it splits into two fragments. Although the shape of the fission barrier obviously depends on the parametrization of the fission path, we have found that remarkably, the maximum of the barrier corresponds to a configuration in which the emerging fragments are already formed and rather well apart. The implication of this finding in the calculation of critical numbers for fission is illustrated in the case of multiply charged Na clusters.

Calculation of the electronic properties of neutral and ionized divalent-metal clusters

The electronic properties of neutral and ionized divalent-metal clusters have been studied using a microscopic theory, which takes into account the interplay between van der Waals (vdW) and covalent bonding in the neutral clusters, and the competition between hole delocalization and polarization energy in the ionized clusters. By calculating the ground-state energies of neutral and ionized. Hg_n clusters, we determine the size dependence of the bond character and the ionization potential I_p(n). For neutral Hg_n clusters we obtain a transition from van del Waals to covalent behaviour at the critical size n_c ~ 10-20 atoms. Results for I_p(Hg_n) with n \le 20 are in good agreement with experiments, and suggest that small Hg_n^+ clusters can be viewed as consisting of a positive trimer core Hg_3^+ surrounded by n - 3 polarized neutral atoms.

On the transition from localized to delocalized electronic states in divalent-metal clusters

The transition from van der Waals to covalent bonding, which is expected to occur in divalent-metal clusters with increasing cluster size, is discussed. We propose a model which takes into account, within the same electronic theory, the three main competing contributions, namely the kinetic energy of the electrons, the Coulomb interactions between electrons, and the s \gdw p intraatomic transitions responsible for van der Waals like bonding. The model is solved by taking into account electron correlations using a generalized Gutzwiller approximation (slave boson method). The occurrence of electron localization is studied as a function of the interaction parameters and cluster size.

Theory for the change of bond character in divalent-metal clusters

To determine the size dependence of the bonding in divalent-metal clusters we use a many-electron Hamiltonian describing the interplay between van der Waals (vdW) and covalent interactions. Using a saddle-point slave-boson method and taking into account the size-dependent screening of charge fluctuations, we obtain for Hg_n a sharp transition from vdW to covalent bonding for increasing n. We show also, by solving the model Hamiltonian exactly, that for divalent metals vdW and covalent bonding coexist already in the dimers.

Transition metal clusters on graphite

The topic of this work is 3d transition metals deposited on graphite. Spin-polarised density-functional calculations are used to obtain the magnetic moments of deposited adatoms and dimers. Interatomic potentials are also deduced. These are used in molecular dynamics simulations to study cluster formation and to investigate cluster morphology.

Interação de átomos leves com clusters de metais de transição

Density Functional Theory (DFT) calculations on the interactions of small atoms (H, C, O, and S) on first-row transition metal clusters were performed. The results show that the adsorption site may vary between the metal surface and the edge of the cluster. The adsorption energies, adatom-nearest neighbor and adatom-metal plane distances were also determined. Finally, the authors present a discussion about the performance of these metals as anodes on solid oxide fuel cells. The results obtained agree with empirical data, indicating that the theoretical model used is adequate