Maurício J. Piotrowski

Transition-metal 13-atom clusters assessed with solid and surface-biased functionals

First-principles density-functional theory studies have reported open structures based on the formation of double simple-cubic (DSC) arrangements for Ru(13), Rh(13), Os(13), and Ir(13), which can be considered an unexpected result as those elements crystallize in compact bulk structures such as the face-centered cubic and hexagonal close-packed lattices. In this work, we investigated with the projected augmented wave method the dependence of the lowest-energy structure on the local and semilocal exchange-correlation (xc) energy functionals employed in density-functional theory. We found that the local-density approximation (LDA) and generalized-gradient formulations with different treatment of the electronic inhomogeneities (PBE, PBEsol, and AM05) confirm the DSC configuration as the lowest-energy structure for the studied TM(13) clusters. A good agreement in the relative total energies are obtained even for structures with small energy differences, e.g., 0.10 eV. The employed xc functionals yield the same total magnetic moment for a given structure, i.e., the differences in the bond lengths do not affect the moments, which can be attributed to the atomic character of those clusters. Thus, at least for those systems, the differences among the LDA, PBE, PBEsol, and AM05 functionals are not large enough to yield qualitatively different results.

Theoretical Investigation of the Adsorption Properties of CO, NO, and OH on Monometallic and Bimetallic 13-Atom Clusters: The Example of Cu13, Pt7Cu6, and Pt13.

We report a density functional theory investigation of the adsorption properties of CO, NO, and OH on the Cu13, Pt7Cu6, and Pt13 clusters in the cationic, neutral, and anionic states with the aim to improve our atomistic understanding of the adsorption properties on bimetallic clusters compared with monometallic clusters. The adsorption energy of CO and NO are substantially stronger on Pt13 than on Cu13, and hence, CO and NO bind preferentially on Pt sites on Pt7Cu6. Thus, it can contribute to drive the migration of the Pt atoms from the core to the surface region in large PtCu nanoalloys. The CO and NO adsorption energies on the bimetallic cluster are enhanced by a few percent compared with the energies of the monometallic clusters, which shows that the Pt-Cu interaction can contribute to an increase in the adsorption energy. In contrast with CO and NO trends, the OH adsorption energies on Cu13, Pt7Cu6, and Pt13 deviates only up to 0.31 eV, and hence, there is no clear preference for Cu or Pt sites on Pt7Cu6 or an enhancement of the adsorption energy on the bimetallic systems. We found a reduction of the CO and NO vibrational frequencies upon adsorption, which indicates a weakening of the CO and NO binding energies, and it is supported by a slight increase in the bond lengths. However, the OH vibrational frequency increases upon adsorption, which indicates an enhancement of the OH binding energy, which is supported by a slight decrease in the bond length by about 0.01 Å. It can be explained by the large charge transfer from the clusters to the O atom, which enhances the electrostatic interaction in the O-H bonding.

A theoretical investigation of the structural and electronic properties of 55-atom nanoclusters: The examples of Y–Tc and Pt

Several studies have found that the Pt55 nanocluster adopts a distorted reduced core structure, DRC55, in which there are 8-11 atoms in the core and 47-44 atoms in the surface, instead of the compact and high-symmetry icosahedron structure, ICO55, with 13 and 42 atoms in the core and surface, respectively. The DRC structure has also been obtained as the putative global minimum configuration (GMC) for the Zn55 (3d), Cd55 (4d), and Au55 (5d) systems. Thus, the DRC55 structure has been reported only for systems with a large occupation of the d-states, where the effects of the occupation of the valence anti-bonding d-states might play an important role. Can we observe the DRC structure for 55-atom transition-metal systems with non-occupation of the anti-bonding d-states? To address this question, we performed a theoretical investigation of the Y 55, Zr55, Nb55, Mo55, Tc55, and Pt55 nanoclusters, employing density functional theory calculations. For the putative GMCs, we found that the Y 55 adopts the ICO55 structure, while Nb55 and Mo55 adopt a bulk-like fragment based on the hexagonal close-packed structure and Tc55 adopts a face-centered cubic fragment; however, Zr55 adopts a DRC55 structure, like Zn55, Cd55, Pt55, and Au55. Thus we can conclude that the preference for DRC55 structure is not related to the occupation of the anti-bonding d-states, but to a different effect, in fact, a combination of structural and electronic effects. Furthermore, we obtained that the binding energy per atom follows the occupation of the bonding and anti-bonding model, i.e., the stability of the studied systems increases from Y to Tc with a small oscillation for Mo, which also explains the equilibrium bond lengths. We obtained a larger magnetic moment for Y 55 (31 μB) which can be explained by the localization of the d-states in Y at nanoscale, which is not observed for the remaining systems (0-1 μB).

Density functional investigation of the adsorption effects of PH3 and SH2 on the structure stability of the Au55 and Pt55 nanoclusters

Although several studies have been reported for Pt55 and Au55 nanoclusters, our atomistic understanding of the interplay between the adsorbate-surface interactions and the mechanisms that lead to the formation of the distorted reduced core (DRC) structures, instead of the icosahedron (ICO) structure in gas phase, is still far from satisfactory. Here, we report a density functional theory (DFT) investigation of the role of the adsorption effects of PH3 (one lone pair of electrons) and SH2 (two lone pairs) on the relative stability of the Pt55 and Au55 nanoclusters. In gas phase, we found that the DRC structures with 7 and 9 atoms in the core region are about 5.34 eV (Pt55) and 2.20 eV (Au55) lower in energy than the ICO model with Ih symmetry and 13 atoms in the core region. However, the stability of the ICO structure increases by increasing the number of adsorbed molecules from 1 to 18, in which both DRC and ICO structures are nearly degenerate in energy at the limit of 18 ligands, which can be explained as follows. In gas phase, there is a strong compression of the cationic core region by the anionic surface atoms induced by the attractive Coulomb interactions (core+-surface-), and hence, the strain release is obtained by reducing the number of atoms in the cationic core region, which leads to the 55 atoms distorted reduced core structures. Thus, the Coulomb interactions between the core+ and surface- contribute to break the symmetry in the ICO55 structure. On the other hand, the addition of ligands on the anionic surface reduces the charge transfer between the core and surface, which contributes to decrease the Coulomb interactions and the strain on the core region of the ICO structure, and hence, it stabilizes a compact ICO structure. The same conclusion is obtained by adding van der Waals corrections to the plain DFT calculations. Similar results are obtained by the addition of steric effects, which are considered through the adsorption of triphenylphosphine (PPh3) molecules on Au55, in which the relative stability between ICO and DRC is the same as for PH3 and SH2. However, for Pt55, we found an inversion of stability due to the PPh3 ligand effects, where ICO has higher stability than DRC by 2.40 eV. Our insights are supported by several structural, electronic, and energetic analyses.