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.

Transition-metal-centered monocyclic boron wheels are important candidates in the family of planar hypercoordinate species that show intriguing structure, stability and bonding situation. Through the detailed potential energy surface explorations of MB9 − (M Fe, Ru, Os) clusters, we introduce herein OsB9 − to be a new member in the transition-metalcentered borometallic molecular wheel gallery. Previously, FeB9 − and RuB9 − clusters were detected by photoelectron spectroscopy and the structures were reported to have singlet D9h symmetry. Our present results show that the global minimum for FeB9 − has a molecular wheel-like structure in triplet spin state with Cs symmetry, whereas its heavier homologues are singlet molecular wheels with D9h symmetry. Chemical bonding analyses show that RuB9 − and OsB9 − display a similar type of electronic structure, where the dual σ + π aromaticity, originated from three delocalized σ bonds and three delocalized π bonds, accounts for highly stable borometallic molecular wheels.

No description available.

No description available.

In quantum chemistry, scientists aim to solve the time-independent Schrödinger equation by employing a variety of approximation techniques whose accuracy are typically inversely proportional to their computational cost. This problem is amplified when it comes to chiral molecules, whose stereochemical assignments and associated chiroptical properties can be incredibly sensitive to small changes in their three-dimensional structure, requiring highly accurate theoretical methods. On the other hand, due to the polynomial scaling with system size, it is sometimes impractical to apply such methods to chemical compounds of broad scientific interest, especially when a multitude of low-energy conformations have to be accounted for as well. As a result, the assignment of absolute configurations to chiral compounds remains a tedious task. However, the characterization of these compounds is something that many different scientists are significantly invested in. The ultimate goal, then, is twofold: to gain useful insight by utilizing the electronic structure methods at your disposal while simultaneously developing new approximation techniques that can be used to push the boundaries on what is currently capable in computational chemistry. Therefore, we start by applying widely accepted density functional theory methods to predict optical rotations and electronic circular dichroism for naturally occurring antiplasmodial and anticancer drug candidates. We find that by comparing the computational results directly with those obtained through experimental measurement, we can provide reliable absolute config- uraitonal assignments to a variety of chiral compounds with numerous stereogenic centers. We also present the first ever prediction of vibrational circular dichroism with second-order Møller-Plesset perturbation theory. This extension opens the door to systematically improvable correlated wave function methods that can be employed when density functional theory fails or when higher accuracy results are required. / Doctor of Philosophy / Theoretical chemistry aims to draw a line from a molecule's three-dimensional structure to a set of physical observables, allowing for the efficient prediction of such properties. One family of chemical compounds for which this task becomes increasingly difficult is known as chiral molecules. A chiral compound is defined as one that has a non-superimposable mirror image. The concept of chirality is most tangibly seen with a pair of human hands, which demonstrate this same mirror-like behavior. In the same way that a person has left and right hands, a three-dimensonal handedness can be used to characterize many compounds that are essential to life including enzymes, sugars, and proteins. Although procedures have been developed to consistently isolate pure samples of such compounds, the correct identification of each hand poses a much larger task and costs the global pharmaceutical industry tens to hundreds of millions of dollars every year. As such, gaining insight about these incredibly valuable compounds and their associated properties is a never ending goal for many scientists. One such way to gain insight is through the direct comparison of experimental and calculated properties, namely chiroptical properties. These unique properties define how chiral compounds interact with light. While experimental scientists are limited in the degree to which they can probe a molecule's structure, theoretical chemists have the advantage of knowing the exact three-dimensional structure for the compound they are studying. On the other hand, theoretical chemists rely on comparison to experimental results to develop new methods or apply the available techniques to predict molecular properties. This work begins by attempting to match calculated properties to experimentally measured ones in order to confirm the detailed molecular structure of natural product drug candidates. Through multiple such computational studies, it is shown that the current methods are sometimes limited in the knowledge that they can provide. As a result, it is absolutely necessary to continue to improve on the existing methods. We go on to provide a first-of-its-kind implementation that allows for theoretical chemists to compare their results to experiment in a way that was not previously possible.

Silicon nanowires (SiNWs) and the copper chalcogenides, namely copper sulfide (Cu2S) and selenide Cu2Se, have diverse applications in renewable energy technology. For example, SiNWs which have direct band gaps unlike bulk Si, have the potential to radically reduce the cost of Si based photovoltaic cells. However, they degrade quickly under ambient conditions. Various surface passivations have therefore been investigated for enhancing their stability but it is not yet well understood how they affect the electronic structure of SiNWs at a fundamental level. Here, we will explore, from first-principles simulation, how fluorine, methyl and hydrogen surface passivations alter the electronic structures of [111] and [110] SiNWs via strain and quantum confinement. We also show how electronic charge states in [111] and [110] SiNWs can be effectively modelled by simple quantum wells. In addition, we address the issue of why [111] SiNWs are less influenced by their surface passivation than [110] SiNWs. Like SiNWs, Cu2S and Cu2Se also make excellent photovoltaic cells. However, they are most well known for their exceptional thermoelectric performance. This is by virtue of their even more unique solid-liquid hybrid nature which combines the low thermal conductivity and good electrical characteristics required for a high thermoelectric efficiency. We use first-principles molecular dynamics simulations to show that Cu diffusion rates in Cu2S and Cu2Se can be as high as 10-5cm2s-1. We also relate their phonon power spectra to their low thermal conductivities. Furthermore, we evaluate the thermoelectric properties of Cu2S and Cu2Se using a combination of Boltzmann transport theory and first-principles electronic structure calculations. Our results show that both Cu2S and Cu2Se are capable of maintaining high Seebeck coefficients in excess of 200μVK-1 for hole concentrations as high as 3x1020cm-3.

Ce travail est une contribution ab initio à la caractérisation d'halogénures diatomiques de métaux de transition. Nous avons choisi de caractériser la structure électronique des chlorures de métaux de transition du groupe Vb (NbCl et TaCl) et du fluorure de nickel car une série de spectres les concernant ont été enregistrés mais aucune donnée théorique fiable n'était disponible dans la littérature. Pour étudier ces molécules, nous avons appliqué une procédure de calcul à deux étapes qui permet de tenir compte des effets de corrélation électronique et des effets relativistes. Dans la première étape, nous avons réalisé des calculs CASSCF/ICMRCI+Q de grande taille qui tiennent compte de l'énergie de corrélation et introduisent des effets relativistes scalaires. Dans la deuxième étape, le couplage spin-orbite est traité par la "state interacting method" implémentée dans le logiciel MOLPRO. Nous avons développé des stratégies de calcul basées sur ces méthodes de calcul et adaptées aux différentes molécules ciblées. Ainsi, pour les molécules NbCl et TaCl, nous avons utilisé des pseudopotentiels relativistes scalaires et spin-orbite, tandis que pour la molécule NiF, nous avons réalisé des calculs tous électrons. Nous avons d'abord testé la stratégie de calcul sur les cations Nb+ et Ta+. Ensuite, nous avons calculé pour la première fois les structures électroniques relativiste scalaire et spin-orbite des molécules NbCl (de 0 à 17000 cm-1) et TaCl (de 0 à 23000 cm-1). A l'aide de ces données théoriques, nous avons interprété les spectres expérimentaux en collaboration avec Bernath et al. Nous avons proposé plusieurs attributions de transitions électroniques en accord avec l'expérience mais nos résultats théoriques ne nous ont pas permis de les attribuer toutes. Néanmoins, nous avons mis en évidence une série d'autres transitions électroniques probables qui pourraient, à l'avenir, servir à l'interprétation de nouveaux spectres mieux résolus. Outre son intérêt expérimental, cette étude a permis de comparer les structures électroniques des molécules isovalencielles VCl, NbCl et TaCl, mettant en évidence des différences importantes. L'élaboration d'une nouvelle stratégie de calcul pour décrire les systèmes contenant l'atome de nickel représentait un véritable défi en raison de la complexité des effets de corrélation électronique. Notre stratégie de calcul a consisté à introduire ces effets en veillant à réduire au maximum la taille des calculs qui devenait considérable. Nous l'avons testée sur l'atome Ni et appliquée ensuite au calcul des structures électroniques relativiste scalaire et spin-orbite de la molécule NiF entre 0 à 2500 cm-1. Nous avons obtenus des résultats qui corroborent l'expérience.

No description available.

No description available.

No description available.