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.

<div>The global demand for energy of mankind, the ever-increasing cost of energy, and the expected depletion of fossil energy carriers within the next centuries urge the exploration of alternative and more sustainable ways to provide energy. The current quest for energy-efficient technologies for the replacement of existing cooling devices has made the magnetocaloric effect a field of current scientific interest. Cooling technologies based on magnetic refrigerants are expected to have a better environmental impact compared with those based on the gas compression-expansion cycle. This technology provides an alternative for refrigeration applications with advantages, such as high energy efficiency, environmental friendliness, and low power consumption. In search of promising magnetocaloric materials, several rare earth-depleted transition metal-based materials were designed and investigated.</div><div>In this work, RCrxAl2-x and RZnAl (R = Gd, Tb, Dy, Ho) belonging to the ternary rare-earth transition-metal Laves phases, were chosen as the starting point to establish the effect of valence electron concentration (VEC) on the magnetic behavior and magnetocaloric effect. Our result and the results from the previously studied RTAl phases (T = Cu, Ni, Co, Fe, Mn) shows that the perturbation of the valence electron concentration at the Fermi level is found to be the driving force that dictates the crystal structure, magnetocaloric and magnetic properties of these systems. Most notably, the decrease in the valence electron concentration at the Fermi level leads to an increase in the curie temperature.</div><div>In addition, we have further extended this theory to GdNiAl2 systems. GdNiAl2 is a known magnetocaloric material which exhibits an isothermal magnetic entropy change of ΔSM = 16.0 Jkg-1K-1 at TC = 28K under a magnetic field change from 0-5T. However, the low TC limits its application as a room temperature refrigerant. We, therefore, substituted Co for (Ni/Al) in the structure of GdNiAl2, intending to substantially perturb the position of the Fermi level of Ni since that will lead to a decrease in the VEC and hence elevate the TC. The study was also extended to another Icosagenides (Ga,), which saw the substitution of Ga for Al in GdNiAl2 and its Co substituted analogs. The Ga analogs exhibit complex magnetic behavior with a cascade (multiple) of magnetic transitions, as opposed to the rather simple magnetism of their Al congeners.</div>