Des del descobriment del buckminster ful.lerè el 1985, s'ha despertat un interés enorme per entendre la reactivitat química així com les propietats d'aquests compostos. La funcionalització exoèdrica del ful.lerè més abundant, el C60, està força ben establerta. Tanmateix, la investigació en aquest camp encara continua oberta ja que s'han sintetitzat una gran varietat de derivats molt prometedors donades les seves futures aplicacions. La tesi comprèn quinze capítols que contenen set publicacions relacionades. Els primers dos estudis es basen en la reacció Diels-Alder sobre els anomenats metal.loful.lerens endoèdrics TNT X3N@C78, X= Sc, Y. Aquest projecte de investigació està motivat pel desconeixament existent sobre les possibles conseqüències de l'encapsulació del grup X3N. El tercer estudi descriu minuciosament els canvis detectats en la funcionalització exoèdrica un cop s'ha produït l'encapsulació dels diferents gasos nobles. En aquesta tesi s'estudia en detall l'ús de l'aproximació ONIOM per a estudiar reaccions de cicloaddició en compostos ful.lerènics. Els resultats d'aquest projecte són d'alt interès per a la realització dels estudis posteriors sobre la reacció de Diels-Alder i la 1,3-dipolar en ful.lerens i derivats. Finalment, l'última part d'aquesta tesi es basa en les propietats antioxidants de determinats ful.lerens. A l'últim treball inclòs en aquesta tesi s'estudia en detall el mecanismo de reacció per a la eliminació del ió superòxid en presència de ful.lerens. / Since the buckminster fullerene discovery in 1985, a huge interest for understanding the chemical reactivity as well as the chemical properties of fullerene compounds has been awakened. The exohedral functionalization of the archetypal compound C60 is nowadays considered to be quite well-established. Still, the research in this field is open as a wide variety of derivatives with intriguing potential applications have been synthesized. The thesis is divided into fifteen chapters that contain seven related publications. The first two studies are based on the Diels-Alder reaction involving the Trimetallic Nitride Template (TNT) endohedral metallofullerenes X3N@C78, X = Sc, Y . This investigation project was basically motivated by the unclear evidence about the possible consequences of the X3N. The third study thoroughly describes the change on the exohedral functionalization upon noble gas encapsulation. In the fourth study included in this thesis, the performance of the ONIOM approach for studying cycloaddition reactions involving fullerene compounds is studied in detail. Results from the latter project are of interest for the following studies involving the 1,3-dipolar and the Diels-Alder cycloaddition reactions where the ONIOM strategy is employed. Finally, the last part of this thesis is based on the antioxidant properties of fullerene compounds, where the mechanism of action for the superoxide removal involving fullerene compounds is unraveled. The understanding of the SOD removal mechanism could represent a big improvement to design new fullerene derivatives with higher antioxidant properties.

Des del seu descobriment, a la molècula C60 se li coneixen una varietat de derivats segons el tipus de funcionalització amb propietats fisicoquímiques específiques de gran interès científic. Una sel·lecció de derivats corresponents a addicions simple o múltiple al C60 s'ha considerat en aquest treball d'investigació. L'estudi a nivell de química computacional de diversos tipus d'addició al C60 s'han portat a terme per tal de poder donar resposta a aspectes que experimentalment no s'entenen o són poc clars.Els sistemes estudiats en referència a l'addició simple al C60 han estat en primer lloc els monoiminoful·lerens, C60NR, (de les dues vies proposades per la seva síntesi, anàlisis cinètic i termodinàmic han ajudat a explicar els mecanismes de formació i justificar l'addició a enllaços tipus [5,6]), i en segon lloc els metanoful·lerens i els hidroful·lerens substituits, C60CHR i C60HR, (raons geomètriques, electròniques, energètiques i magnètiques justifiquen el diferent caràcter àcid ente ambdós derivats tenint en compte una sèrie de substituents R amb diferent caràcter electrònic donor/acceptor). Els fluoroful·lerens, C60Fn, i els epoxid ful·lerens, C60On, (anàlisi sistemàtic dels seus patrons d'addició en base a poder justificar la força que els governa han aportat dades complementàries a les poques que existeixen experimentalment al respecte). / Since the discovery of C60 molecule a large number of derivatives molecules have been described with a great scientific interest of their specific physical and chemical properties. A selection of single and multiple addition products has been considered in this investigation. Study at Computational Chemistry level for this selected derivatives have been carried out in order to give answer to several points that experimentally are not understandable or not enough clear.As single addition derivatives, firstly were studied the monoiminofullerenes, C60NR, (two routes of synthesis have been considered, kinetic and thermodynamic analysis have help to explain formation mechanisms and justify the possible addition at [5,6]-type bonds), and secondly the methanofullerenes and the substituted hydrofullerenes, C60CHR and C60HR, (geometric, electronic, energetic and magnetic reasons justify the different acid character between both series of derivatives taking care on the influence of R substituents with different donor/acceptor character). The fluorofullerenes, C60Fn, and the epoxide fullerenes, C60On, have been studied as multiple addition derivatives (systematic analysis of the addition pattern have been performed in the way to find reasons to justify the driving force of the multiple addition process).

In this thesis, six small molecules (S-allylcysteine, S-allyl mercaptocysteine, allicin, methyl propyl disulfide, allyl methyl sulfide and dipropylsulfide) that are found in garlic and onion, and are known to be beneficial for human health were studied using molecular mechanics, semi-empirical methods, ab-initio (Restricted Hartree Fock), and density functional theory. Using the same methods, a synthetic pyrethroid pesticide molecule, called cyfluthrin, was also studied. Structural, vibrational and electronic properties of these molecules were found. These theoretical studies could clarify the role of these molecules on human health before they are commercially developed and used. In addition, unfolding dynamics of small peptide sequences (DDATKTFT and its variants) in immunoglobulin G-binding protein G was investigated. Protein folding and unfolding is one of the most important unsolved problems in molecular biology. Because of the large number of atoms involved in protein folding, it is a massive computational problem. The hope is that, one could understand this mechanism with the help of molecular dynamics simulation on small peptides. One of our findings is that the location of the hydrogen bonds is important for the stability of the peptide.

Computational chemists are concerned about two aspects when choosing between the myriad of theoretical methodologies: the accuracy (the "truth") and the computational cost (the tractability). Among the least expensive methods are the Hartree-Fock (HF), density functional theory (DFT), and second-order Moller-Plesset perturbation theory (MP2) methods. While each of these methods yield excellent results in many cases, the inadequate inclusion of certain types of electron correlation (either high-orders or nondynamical) can produce erroneous results. The compromise for the computation of noncovalent interactions often comes from empirically scaling DFT and/or MP2 methods to fit benchmark data sets. The DFT method with an empirically fit dispersion term (DFT-D) often yields semi-quantitative results. The spin-component scaled MP2 (SCS-MP2) method parameterizes the same- and opposite-spin correlation energies and often yields less than 20% error for prototype noncovalent systems compared to chemically accurate CCSD(T) results. There is no simple fix for cases with a large degree of nondynamical correlation (such as transition metal-salen complexes). While testing standard and new DFT functionals on the spin-state energy gaps of transition metal-salen complexes, no DFT method produced reliable results compared to very robust CASPT3 results. Therefore each metal-salen complex must be evaluated on a case-by-case basis to determine which methods are the most reliable. Utilizing a combination of DFT-D and SCS-MP2 methods, the reaction mechanism for the addition of cyanide to unsaturated imides catalyzed by the Al(Cl)-salen complex was performed. Various experimental observations are rationalized through this mechanism.

A molecule can crystallise in more than one crystal structure, a common phenomenon in organic compounds known as polymorphism. Different polymorphic forms may have significantly different physical properties, and a reliable prediction would be beneficial to the pharmaceutical industry. However, crystal structure prediction (CSP) based on the knowledge of the chemical structure had long been considered impossible. Previous failures of some CSP attempts led to speculation that the thermodynamic calculations in CSP methodologies failed to predict the kinetically favoured structures. Similarly, regarding the stabilities of co-crystals relative to their pure components, the results from lattice energy calculations and full CSP studies were inconclusive. In this thesis, these problems are addressed using the state-of-the-art CSP methodology implemented in the GRACE software. Firstly, it is shown that the low-energy predicted structures of four organic molecules, which have previously been considered difficult for CSP, correspond to their experimental structures. The possible outcomes of crystallisation can be reliably predicted by sufficiently accurate thermodynamic calculations. Then, the polymorphism of 5- chloroaspirin is investigated theoretically. The order of polymorph stability is predicted correctly and the isostructural relationships between a number of predicted structures and the experimental structures of other aspirin derivatives are established. Regarding the stabilities of co-crystals, 99 out of 102 co-crystals and salts of nicotinamide, isonicotinamide and picolinamide reported in the Cambridge Structural Database (CSD) are found to be more stable than their corresponding co-formers. Finally, full CSP studies of two co-crystal systems are conducted to explain why the co-crystals are not easily obtained experimentally.

The understanding of molecular effects in nanoscale environments is becoming increasingly relevant for various emerging fields. These include spectroscopy for molecular identification as well as in finding molecules for energy harvesting. Theoretical quantum chemistry has been increasingly useful to address these phenomena to yield an understanding of these effects. In the first part of this dissertation, we study the chemical effect of surface-enhanced Raman scattering (SERS). We use quantum chemistry simulations to study the metal-molecule interactions present in these systems. We find that the excitations that provide a chemical enhancement contain a mixed contribution from the metal and the molecule. Moreover, using atomistic studies we propose an additional source of enhancement, where a transition metal dopant surface could provide an additional enhancement. We also develop methods to study the electrostatic effects of molecules in metallic environments. We study the importance of image-charge effects, as well as field-bias to molecules interacting with perfect conductors. The atomistic modeling and the electrostatic approximation enable us to study the effects of the metal interacting with the molecule in a complementary fashion, which provides a better understanding of the complex effects present in SERS. In the second part of this dissertation, we present the Harvard Clean Energy project, a high-throughput approach for a large-scale computational screening and design of organic photovoltaic materials. We create molecular libraries to search for candidates structures and use quantum chemistry, machine learning and cheminformatics methods to characterize these systems and find structure-property relations. The scale of this study requires an equally large computational resource. We rely on distributed volunteer computing to obtain these properties. In the third part of this dissertation we present our work related to the acceleration of electronic structure methods using graphics processing units. This hardware represents a change of paradigm with respect to the typical CPU device architectures. We accelerate the resolution-of-the-identity Moller-Plesset second-order perturbation theory algorithm using graphics cards. We also provide detailed tools to address memory and single-precision issues that these cards often present.

No description available.

The present state of modeling radio-induced effects at the cellular level neglects to account for the microscopic inhomogeneity of the nucleus from the non-aqueous contents by approximating the entire cellular nucleus as a homogenous medium of water. Charged particle track-structure calculations utilizing this principle of superposition are thereby neglecting to account for approximately 30% of the molecular variation within the nucleus. To truly understand what happens when biological matter is irradiated, charged particle track-structure calculations need detailed knowledge of the secondary electron cascade, resulting from interactions with not only the primary biological component – water – but also the non-aqueous contents, down to very low energies. This paper presents developments for a novel approach, which to our knowledge has never been done before, to reducing the homogenous water approximation. The purpose of our work is to develop of a completely self-consistent computational method for predicting molecule-specific ionization, excitation, and scattering cross sections in the very low energy regime that can be applied in a condensed history Monte Carlo track-structure code. The present methodology begins with the calculation of a solution to the many-body Schrödinger equation and proceeds to use Monte Carlo methods to calculate the perturbations in the internal electron field to determine the aforementioned processes. Results are computed for molecular water in the form of linear energy loss, secondary electron energies, and ionization-to-excitation ratios and compared against the low energy predictions of the GEANT4-DNA physics package of the Geant4 simulation toolkit.

A summary of research conducted since September 2007 at the University of Toronto in the laboratory of Professor Robert A. Batey is presented in this thesis, which is divided into four chapters. The first chapter contains a two-part introduction, where aryl- and aliphatic-Claisen rearrangements are discussed in part 1, and the nucleophilic addition of organoboron reagents to unsaturated C–N functionalities is described in part 2. Chapter 2 contains research involving synthetic and theoretical studies of aryl-Claisen rearrangements and other sigmatropic reactions. The work towards developing the lanthanide-catalyzed domino aryl-Claisen rearrangement for the synthesis of contiguous aryl–C(sp³) moieties is presented first. This is followed by computational studies involving E/Z-selectivity differences for the aryl-Claisen rearrangement, which was an issue noted for the domino aryl-Claisen reaction of a linear substrate. The mechanistic origins of E/Z-selectivity differences for the mono aryl-Claisen rearrangement, which was experimentally ambiguous for over 40 years, is resolved through computational methods. A theoretical analysis of selectivity differences for the allylic azide rearrangement is also described. The third section contains a discussion of Eu(fod)3-catalyzed aryl-Claisen rearrangements on vinyl bromide systems and preliminary studies involving application of the substrates in cross-coupling reactions, and other attempted mono- and domino sigmatropic rearrangements are presented in the fourth section. In chapter 3, the search for computational methods that can accurately predict experimental free energy of activation barriers for the aliphatic-Claisen rearrangement through benchmarking studies with a priori kinetic barrier and kinetic isotope effect data is described. Methods were found to predict new valid transition states and predict ΔG‡ values with a mean unsigned error of 0.3 kcal/mol relative to experimental values. In chapter 4, the development of new allylboration reaction is outlined, involving the double allylboration of nitriles and anhydrides, and initial studies towards the first aminoallylboration reactions of N-aluminoaldimines to form 1,2-diamines.

As the demand for clean and renewable energy sources continues to grow, much attention has been given to solid oxide fuel cells (SOFCs) due to their efficiency and low operating temperature. However, the components of SOFCs must still be improved before commercialization can be reached. Of particular interest is the solid electrolyte, which conducts oxygen ions from the cathode to the anode. Samarium-doped ceria (SDC) is the electrolyte of choice in most SOFCs today, due mostly to its high ionic conductivity at low temperatures. However, the underlying principles that contribute to high ionic conductivity in doped ceria remain unknown, and so it is difficult to improve upon the design of SOFCs. This thesis focuses on identifying the atomistic interactions in SDC which contribute to its favourable performance in the fuel cell. Unfortunately, information as basic as the structure of SDC has not yet been found due to the difficulty in experimentally characterizing and computationally modelling the system. For instance, to evaluate 10.3% SDC, which is close to the 11.1% concentration used in fuel cells, one must investigate 194 trillion configurations, due to the numerous ways of arranging the Sm ions and oxygen vacancies in the simulation cell. As an exhaustive search method is clearly unfeasible, we develop a genetic algorithm (GA) to search the vast potential energy surface for the low-energy configurations, which will be most prevalent in the real material. With the GA, we investigate the structure of SDC for the first time at the DFT+U level of theory. Importantly, we find key differences in our results from prior calculations of this system which used less accurate methods, which demonstrate the importance of accurately modelling the system. Overall, our simulation results of the structure of SDCagree with experimental measurements. We identify the structural significance of defects in the doped ceria lattice which contribute to oxygen ion conductivity. Thus, the structure of SDC found in this work provides a basis for developing better solid electrolytes, which is of significant scientific and technological interest. Following the structure search, we perform an investigation of the electronic properties of SDC, to understand more about the material. Notably, we compare our calculated density of states plot to XPS measurements of pure and reduced SDC. This allows us to parameterize the Hubbard (U) term for Sm, which had not yet been done. Importantly, the DFT+U treatment of the Sm ions also allowed us to observe in our simulations the magnetization of SDC, which was found by experiment. Finally, we also study the SDC surface, with an emphasis on its structural similarities to the bulk. Knowledge of the surface structure is important to be able to understand how fuel oxidation occurs in the fuel cell, as many reaction mechanisms occur on the surface of this porous material. The groundwork for such mechanistic studies is provided in this thesis.