Clusters e nanoclusters têm recebido grande atenção devido à suas propriedades físicas e químicas, as quais divergem bastante dos materiais na fase bulk. Essas propriedades podem variar de acordo com a composição e tamanho do cluster. Uma compreensão da evolução das propriedades em relação a estes parâmetros é de grande importância para potencializar diversas aplicações, entretanto, esse entendimento permanece insatisfatório. Este trabalho foi dividido em duas etapas, em que a primeira busca investigar parâmetros energéticos, por meio do cálculo da energia de excesso, e estruturais, analisando parâmetro de ordem química, função de distribuição radial central, comprimento médio de ligação e número de coordenação efetiva, dos sistemas ZrnCum-n, ZrnAlm-n, CunAlm-n e ZrnCunAlm-2n para n = 55 e 561 átomos com o incremento n tomado de 1 em 1 para o sistema de 55 átomos e de 20 em 20 para os sistemas de 561 átomos. A segunda etapa consiste de investigar como variam as propriedades energéticas e estruturais do sistema ZrCu em função da evolução do tamanho do sistema. Para alcançar os objetivos propostos, neste trabalho foi usado o algoritmo de otimização global de clusters e nanopartículas basin-hopping Monte Carlo revisado. O potencial de interação atômica utilizado é o método do átomo imerso, que é bastante utilizado na descrição de sistemas metálicos. Os resultados obtidos sugerem que: (i) os sistemas puros apresentaram energia de coesão mais alta que seu análogo material na fase bulk, sugerindo que estes tendem a aglomerar-se formando estruturas bulk. Para os sistemas binários e ternários, foi identificado que todas as composições são energeticamente estáveis devido aos valores negativos obtidos pelo excesso de energia e, para o sistema ZrCu verificou-se a presença de efeitos de tamanho. (ii) Com relação à estrutura, as composições puras estudadas apresentaram simetria icosaédrica. Para o estudo da evolução do tamanho do sistema, Zr e Cu apresentaram estrutura com simetria icosaédrica até a composição de 561 átomos, além deste tamanho a simetria icosaédrica é quebrada. Para os sistemas binários e ternários foi obtido que os átomos tendem a distribuir-se dentro do nanocluster além de apresentarem quebra da simetria icosaédrica apresentando ausência de camadas atômicas ordenadas acompanhada de redução da coordenação efetiva. Os sistemas ZrCu e ZrAl demonstraram seguir a lei de Vegard, enquanto que os sistemas CuAl e ZrCuAl apresentaram desvio da lei de Vegard providos por efeitos eletrônicos, além de apresentarem a presença de efeitos de tamanho. / Clusters and nanocluster have attracted great attention due to their physical and chemical properties, very different from their analogous bulk. These properties can vary with composition and size cluster. An understanding of the properties evolution with respect these parameters is essential to improve several applications. However, this understanding is not complete. This study was piecemeal in two stage, being the first the investigation of energetic properties, by excess energy analisys, and structural properties, by chemical order parameter, radial distribution function, effective coordination number and average bond length, from ZrnCum-n, ZrnAlm-n, CunAlm-n and ZrnCunAlm-2n systems, where n = 55, 561 atoms and the increment n vary in one unit for 55-atoms system and twenty unit for 561-atoms system. The second stage is the investigation of how vary the energetic and structural properties from the size evolution ZrCu system. To do this study, was employed the global optimization algorith for cluster and nanoparticle Revised basin-hopping Monte Carlo, were this method use the classical calculation to determine the total energy of the system. The interatomic potential used was the embedded atom method, that was very usefull to describe metallic systems. Our results suggest: (i) the unary systems present cohesive energy higher than their analogous bulk, that indicate the trend of clusters to form bulk. To the binary and ternary systems, we had that all systems are favorable to form nanoalloys by negative value of excess energy. From ZrCu system, the stability decrease when increase the size of system. With respect the structure, the unary compounds present icosahedral symmetry. From the size-evolution study, the unary compounds present icosahedral symmetry until 561-atoms composition, after this size the icosahedral symmetry is broken. To binary and ternary systems, the atoms trend form mixture into the nanocluster, the icosahedral symmetry is broken with respect the unary compounds and presenting absence of ordered layers followed by effective coordination reduction. The ZrCu and ZrAl systems follow the Vegard law, while the CuAl and ZrCuAl systems present deviation from Vegard law, because electronic effects.

Thesis (Ph.D. (Physics)) --University of Limpopo, 2011 / The importance of gold for scientific uses is of fundamental importance to research and technology developments. The bulk gold shows reluctance to participate in chemical reactions, the effect which has been corrected by the change in the size towards nanoclusters. It is therefore imperative that the structure of gold nanomaterials is understood for better applications in catalysis and other developments. Molecular dynamics and the density functional theory have proven to be good tools in computational material science and have thus been used to greater lengths. Molecular dynamics simulations on different gold nanoclusters and nanotubes were successfully carried out at different thermodynamic conditions. The effect of size on the melting of materials was duly tested and our results to some extend agree with what has already been reported. Gold nanoclusters show melting below the bulk and the melting temperatures increase with cluster size. However, the Au55 cluster shows different results in that it melts above the bulk due to structural reconstruction. The structure of the clusters changes from spherical shapes to tetragonal or face centred cubic (fcc) structures. Gold nanotubes show no resistance to temperature and different configurations are obtained in different ensembles. Single wall nanotubes form spherical clusters in the NVT while the NPT conditions give patches of clusters at elevated temperatures. The multi wall nanotubes also form spherical clusters in the NVT but fcc structures are obtained in the NPT Berendsen ensemble towards melting. Ab initio calculations in DMOL3 code on different gold nanoclusters show the stability of the clusters to increase with size and the Au3 and Au8 clusters contain the most stable structures. The Au-Au bond length in the dimer was obtained to within reasonable agreement with experiments and other theoretical works. Doping of the clusters further improved their stability although different impurities give different observations. The QMERA code calculations show that a gold atom on top of the surface causes slanting of the outer MD layers. The morphology of the quantum atoms also changes as compared to the neutral surface and the results are compared by the DMOL3 code which confirms the QMERA results. / Mintek, and the National Research Foundation

Low-dimensional magnetic systems are of interest due to several new effects and modifications that occur at sizes below the average domain grain boundary within the bulk material. Molecule-like magnetite (Fe3O4) nanoparticles, with sizes ranging from one to two nm were synthesized and characterized in order to investigate new properties arising from quantum size effects. These small systems will provide opportunities to investigate magnetism of zero-dimension systems. A zero-dimensional object is usually called a quantum dot or artificial atom because its electronic states are few and sharply separated in energy, resembling those within an atom. Since the surface to volume ratio is the highest for zero-dimensional systems, most of the changes to magnetic behavior will be observed in ultra-fine magnetic particles. Chemically functional magnetic nanoparticles, comprised of a Fe3O4 magnetite core encased in a thin aliphatic carboxylate, have been prepared by sequential high temperature decomposition of organometallic compounds in a coordinating solvent. In this work, aliphatic carboxylic acid chain length, reaction temperature and duration were varied to produce small core diameters. In order to correlate size effects with changes in particle formation, it is important to have a through understanding of the structural components. This includes studies of the core size, surface effects, decomposition, electronic properties and magnetic behavior. Quantum size effects were observed in the (Fe3O4)X(carboxylate)Y monolayer protected clusters (MPCs) when the average core diameter was ≤ 2.0 nm, evidenced by a blue shifted absorbance band maxima, suggesting the onset of quantum confinement. These (Fe3O4)X(carboxylate)Y MPCs also posses a complex interplay between surface and finite size effects, which govern the magnetic properties of these zero-dimensional systems. These MPCs are all superparamagnetic above their blocking temperatures with total magnetic anisotropy values greater than the bulk value due to an increase in surface and magnetocrystalline anisotropy. A non-linear decrease in saturation magnetization (MS) [Bohr Magneton] per cluster) as a function of the reciprocal of core radius have been attributed to surface effects such as a magnetically inactive layer or an increase in spin disorder as core diameter decreases. The reduced core dimensions of these MPCs make them ideal candidates for further investigation of quantum magnetic systems.

The graphene moiré superstructure formed on Ru(0001) (g/Ru(0001)) has shown the potential as a template to self-assemble super-lattices of metal nanoparticles as model catalysts. To explore the possibility of rational catalyst design on g/Ru(0001), detailed density functional theory (DFT) calculations have been performed to investigate the adsorption and diffusion of Rh and Au adatoms on g/Ru(0001). The consequences of different hopping rates for cluster nucleation have been explored by performing Monte Carlo-based statistical analysis, which suggests that diffusing species other than adatoms need to be taken into account to develop an accurate description of cluster nucleation and growth on this surface. DFT calculations have also been carried out to investigate the adsorption and diffusion of 18 4d (Y-Ag) and 5d (La-Au) transition metal adatoms on g/Ru(0001). Given the necessity to study larger diffusing species than adatoms, DFT calculations have been performed to study the adsorption and diffusion of Rh and Au dimers and trimers on g/Ru(0001). It was shown that the mobility of Rh clusters decreases with the increase of cluster size, while for Au, dimers diffuse faster than monomers and trimers on the moiré surface. We then used a genetic algorithm combined with DFT calculations to predict the lowest energy structure of a Au8 cluster on g/Ru(0001). Our prediction leads us to propose that Au clusters aggregates through Oswald ripening with Au dimer being the major diffusing species. Finally, we examined the morphology of a Cu19 cluster on g/Cu(111) using MD simulations with COMB3 potential. We also studied the mobility of Cu clusters on g/Cu(111) at elevated temperatures. The analysis suggests that g/Cu(111) may not be a suitable substrate for the formation and growth of isolated Cu clusters. All these calculation results have provided us a better understanding and useful insights into the nucleation and growth mechanism of metal clusters on graphene moiré.

Doctor of Philosophy / Department of Physics / Talat S. Rahman / In this work we present a study of the geometric, electronic, vibrational and magnetic properties of several nanostructured systems for which experimental data call for a theoretical understanding. In order to investigate the effect of magnetic dipolar interactions on the magnetization of nanomagnets arranged in finite lattices, we utilize a phenomenological classical approach, which is based on the Landau-Lifshitz equation. Dipolar interactions lead to hysteretic behavior of the magnetization curves and established that the external field sweep rate, sample temperature, and shape anisotropy play a role in determining the specifics. Our results (derived from a classical approach) for magnets arranged in a square lattice suggest that stepped hysteresis curves do not have necessarily a quantum origin (quantum tunneling of the magnetization). We also find that in the square lattice small changes in the dipolar strength introduce sudden transitions in the magnetic hysteresis. For the examination of geometric vibrational and electronic structure of systems of interest, we turn to density functional theory (DFT), which is the leading technique for modeling nanoscale systems from first principles. We have applied DFT to either address some old queries of surface science, such as the dynamics of the CO-chemisorbed Cu(001) surface, or to contribute to the forefront of hydrogen-based economy through the comprehension of the growth and diffusion of Pt islets on Ru(0001), or to predict the geometric and electronic properties of materials to-be-created, as in the case of core-shell bimetallic nanoclusters. In the case of CO on Cu(001), although the bond has been considered to be weak enough so as to treat the adsorbate and substrate separately, our calculations are able to reproduce measurements and provide evidence that the dynamics of the molecule is influenced by the substrate and vice versa, as well as by intermolecular interactions. Taking into account the adsorbate-substrate interplay, has furthermore clarified issues that were pending for the clean surface and led to the correct interpretation of some features in the phonon dispersion of the chemisorbed surface. DFT has also directed us to the conclusion that the catalytic properties of few-atom Pt islets on Ru nanoclusters are preserved by the low probability of these islets to diffuse through the edges of the Ru nanoclusters. Moreover, the analysis of the Ag_{27}Cu_7 nanoalloy from ab initio methods has opened a wide panorama in terms of the geometry, coordination, energetics, and electronic structure of alloyed phases, in general,that may aid in the assembling on new materials.

Ce projet concerne la synthèse et caractérisation de nouveaux matériaux composites à base de nanoclusteurs de thiolates d’or Au25(SR)18 (tGNCs), supportés sur divers polymères de coordination (MOFs), ainsi que sur ZrO2. L’activité catalytique de ces matériaux a été évaluée sur la transformation de différents substrats. Les tGNCs sont des matériaux atomiquement bien définis et connus pour être actifs dans des réactions d’oxydation. Les nanoparticules de MOFs sont des matériaux pouvant servir de support pour des tGNCs avec de bonnes dispersions. Certains MOFs sont connus pour avoir des propriétés acides et peuvent être actifs en catalyse. Parmi eux, MIL-101 (Cr), UiO-66 (Zr) et ZIF-8 (Zn) on été choisis en raison de leur propriétés acides et/ou de stabilité thermique. La synergie entre les tGNCs et les MOFs a été évaluée à travers la conversion catalytique de différents substrats tels le glucose, le fructose, l’alcool benzylique et le furfural, impliquant des étapes nécessitant un caractère acide et/ou oxydant. Globalement, il n’a pas été observé d’impact de la présence d’or sur la réactivité de ces substrats, et les tendences catalytiques sont celles obtenues avec les MOFs seuls. Cela est certainement dû à la stabilité thermique non suffisante des MOFs qui prévient une calcination efficace des tNGCs. Lorsque ces clusters sont déposés sur ZrO2, il a été possible de les calciner à différentes températures pour étudier l’effet du ligand et de la taille de particules, pour des réactions d’oxydation en phase liquide. Ainsi, il a été montré par exemple que la température de calcination a un impact significatif sur le comportement catalytique de ces composites, qui ont donné de bonnes activités pour l’oxydation de l’alcool benzylique en benzaldéhyde dans le toluène et en conditions douces, et pour l’esterification oxydante du furfural en furoate de méthyle / This research project reports the synthesis and characterization of new composite materials based on Au25(SR)18 thiolate gold nanoclusters (tGNCs), supported over a range of metal organic frameworks (MOFs), and ZrO2. The synthesized composite materials were tested for catalytic transformations of various substrates. tGNCs are atomically well defined materials known to be active in oxidation reactions. MOFs nanoparticles are materials suitable for high dispersion of tGNCs. Some MOFs are known to have acidity and can be active as catalysts. Among them, MIL-101 (Cr), UiO-66 (Zr) and ZIF-8 (Zn) were chosen due to their acidic and/or thermal stability properties. The synergy between tGNCs and MOFs has been tested through catalytic conversions of different substrates like glucose, fructose, benzylalcohol and furfural, involving steps requiring acidic and oxidative features. Globally, no impact of the presence of Au clusters was observed, and the composite materials showed the same catalytic trends as those obtained with the MOFs alone. This is mainly due to the not sufficient thermal stability of the MOFs that prevents efficient calcination of the tGNCs. In contrast, when deposited on ZrO2 it was possible to calcine Au25(SG)18 nanoclusters at different temperatures to study the ligand and particle size effects in liquid phase oxidation reactions. For example, the calcination temperature had a significant impact on the catalytic behaviour of this composite materials, which showed good activity for the oxidation of benzyl alcohol into benzaldehyde in toluene under mild conditions, and of furfural oxidative esterification into methyl-2-furoate

This dissertation focuses on computational study of the geometry and energetics small molecules and nanoclusters involving transition metals (TM). These clusters may be used for various industrial applications including catalysis and photonics. Specifically, in this work we have studied hydrides and carbides of 3d-transition metal systems (Sc through Cu), small nickel and gold clusters. Qualitatively correct description of the bond dissociation is ensured by allowing the spatial and spin symmetry to break. We have tested applicability of new exchange-correlation functional and alternative theoretical descriptions (spin-contamination correction in broken symmetry DFT and ensemble Kohn-Sham (EKS)) as well. We studies TM hydrides and carbides systems to understand the importance of underlying phenomenon of bond breaking in catalytic processes. We have tested several exchange-correlation functionals including explicit dependence on kinetic energy density for the description of hydrides (both neutral and cationic) and carbides formed by 3d-transition metals. We find M05-2x and BMK dissociation energies are in better agreement with experiment (where available) than those obtained with high level wavefunction theory methods, published previously. This agreement with experiment deteriorates quickly for other functionals when the fraction of the Hartree-Fock exchange in DFT functional is decreased. Higher fraction of HF exchange is also essential in EKS formalism, but it does not help when spin-adapted unrestricted approach is employed. We analyze the electron spin densities using Natural Bond Orbital population analysis and find that simple description of 3d electrons as non-bonding in character is rarely correct. Unrestricted formalism results in appreciable spin-contamination for some of the systems at equilibrium, which motivated us to investigate it further in details. In order to correct the spin contamination effect on the energies, we propose a new scheme to correct for spin contamination arising in broken-symmetry DFT approach. Unlike conventional schemes, our spin correction is introduced for each spin-polarized electron pair individually and therefore is expected to yield more accurate energy values. We derive an expression to extract the energy of the pure singlet state from the energy of the broken-symmetry DFT description of the low spin state and the energies of the high spin states (pentuplet and two spin-contaminated triplets in the case of two spin-polarized electron pairs). We validate our spin-contamination correction approach by a simple example of H2 and applied to more complex MnH system. Ensemble KS formalism is also applied to investigate the dissociation of C2 molecule. We find that high fraction of HF exchange is essential to reproduce the results of EKS treatment with exact exchange-correlation functional. We analyze the geometry and energetics of small nickel clusters (Ni2-Ni5) for several lowest energy isomers. We also study all possible spin states of small nickel cluster isomers and report observed trends in energetics. Finally we determine the geometry and energetics of ten lowest energy isomers of four small gold clusters (Au2, Au4, Au6, and Au8). We have also investigated the influence of cluster geometry, ligation, solvation and relativistic effects on electronic structure of these gold clusters. The effect of one-by-one ligand attachment in vacuum and solvent environment is also studied. Performance of five DFT functionals are tested as well; Local Spin Density Approximation (SVWN5), Generalized Gradient Approximation (PBE), kinetic energy density-dependent functional (TPSS), hybrid DFT (B3LYP), and CAM-B3LYP which accounts for long-range exchange effects believed to be important in the analysis of metal bonding in gold complexes and clusters. Our results exhibit the ligand induced stability enhancement of otherwise less stable isomers of Au4, Au6 and Au8. Ligands are found to play a crucial role in determining the 2D to 3D transition realized in small gold clusters. In order to select an appropriate theory level to use in this study, we investigate the effect of attachment of four different ligands (NH3, NMe3, PH3, PMe3) on cluster geometry and energetics of Au2 and Au4 in vacuum and in solution. Our results benchmark the applicability of DFT functional model and polarization functions in the basis set for calculations of ligated gold cluster systems. We employ five different basis sets with increasing amount of polarization and diffuse functions; LANL2DZ, LANL2DZ-P, def2-SVP, def2-TZVP, and def2-QZVP. We obtain NMe3 = NH3 > PH3 > PMe3 order of ligand binding energies and observe shallow potential energy surfaces in all molecules. Our results suggest appropriate quantum-chemical methodologies to model small noble metal clusters in realistic ligand environment to provide reliable theoretical analysis in order to complement experiments.

The development of new storage materials will facilitate the use of hydrogen as a major energy carrier in near future. Among the chemical hydrides used as hydrogen storage materials for supplying hydrogen at ambient temperature, sodium borohydride seems to be an ideal one because it is stable under ordinary conditions and liberates hydrogen gas in a safe and controllable way in aqueous solutions. However, self hydrolysis of sodium borohydride is so slow that requires a suitable catalyst. This work aims the use of water dispersible cobalt(0) nanoclusters having large portion of atoms on the surface as catalyst for the hydrolysis of sodium borohydride. In-situ formation of cobalt(0) nanoclusters and catalytic hydrolysis of sodium borohydride were performed starting with a cobalt(II) chloride as precursor and sodium borohydride as reducing agent and substrate in the presence of a water soluble stabilizer. As stabilizer, water soluble polyacrylic acid as well as hydrogen phosphate ion were tested. Cobalt(0) nanoclusters were characterized by using all the available analytical methods including FT-IR, TEM, XPS, UV-visible electronic absorption spectroscopy. The kinetics of cobalt(0) nanoclusters catalyzed hydrolysis of sodium borohydride were studied depending on the catalyst concentration, substrate concentration, stabilizing agent concentration and temperature.

Computational approaches play an important role in today's materials science owing to the remarkable advances in modern supercomputing architecture and algorithms. Ab initio simulations solely based on a quantum description of matter are now very able to tackle materials problems in which the system contains up to a few thousands atoms. This dissertation aims to address the modern electronic structure calculation methods applied to a range of various materials such as liquid and amorphous phase materials, nanostructures, and small organic molecules. Our simulations were performed within the density functional theory framework, emphasizing the use of real-space ab initio pseudopotentials. On the first part of our study, we performed liquid and amorphous phase simulations by employing a molecular dynamics technique accelerated by a Chebyshev-subspace filtering algorithm. We applied this technique to find l- and a- SiO₂ structural properties that were in a good agreement with experiments. On the second part, we studied nanostructured semiconducting oxide materials, i.e., SnO₂ and TiO₂, focusing on the electronic structures and optical properties. Lastly, we developed an efficient simulation method for non-contact atomic force microscopy. This fast and simple method was found to be a very powerful tool for predicting AFM images for many surface and molecular systems. / text

The development of new storage materials will facilitate the use of hydrogen as a major energy carrier in near future. In hydrogen economy, chemical hydrides such as NaBH4, KBH4, LiH, NaH have been tested as hydrogen storage materials for supplying hydrogen at ambient temperature. Among these chemical hydrides, sodium borohydride seems to be an ideal hydrogen storage material because it is stable under ordinary conditions and liberates hydrogen gas in a safe and controllable way in aqueous solutions. However, self hydrolysis of sodium borohydride is so slow that it requires a suitable catalyst. All of the prior catalysts tested for the hydrolysis of sodium borohydride are heterogeneous and, therefore, have limited activity because of the small surface area. Here, we report for the first time the employment of water dispersible metal(0) nanoclusters having a large portion of atoms on the surface as a catalyst for the hydrolysis of sodium borohydride. In-situ formation of nickel(0) nanoclusters and catalytic hydrolysis of sodium borohydride were performed in the same medium. Nickel(0) nanoclusters are prepared from the reduction of nickel(II) acetylacetonate by sodium borohydride in aqueous solution and stabilized with hydrogenphosphate anions. The nickel(0) nanoclusters were characterized by using XPS, Powder XRD, FT-IR, UV-Vis and NMR spectroscopic methods. The kinetics of the nickel(0) nanoclusters catalyzed hydrolysis of sodium borohydride was studied depending on the catalyst concentration, substrate concentration, stabilizing agent concentration and temperature. Tha kinetic study shows that the nickel(0) nanocluster-catalyzed hydrolysis of sodium borohydride is first order with respect to catalyst concentration and zero order with respect to substrate concentration The activation parameters of this reaction were also determined from the evaluation of the kinetic data. The hydrogenphosphate-stabilized nickel(0) nanoclusters provide a lower activation energy (Ea= 55 kJ/mol) than bulk nickel (Ea=73 kJ/mol) for the hydrolysis of sodium borohydride.