The Quantum chemical cluster approach has been shown to be quite powerful and efficient in the modeling of enzyme active sites and reaction mechanisms. In this thesis, the reaction mechanisms of several enzymes have been investigated using the hybrid density functional B3LYP. The enzymes studied include four dinuclear zinc enzymes, namely dihydroorotase, N-acyl-homoserine lactone hydrolase, RNase Z, and human renal dipeptidase, two trinuclear zinc enzymes, namely phospholipase C and nuclease P1, two tungstoenzymes, namely formaldehyde ferredoxin oxidoreductase and acetylene hydratase, aspartate α-decarboxylase, and mycolic acid cyclopropane synthase. The potential energy profiles for various mechanistic scenarios have been calculated and analyzed. The role of the metal ions as well as important active site residues has been discussed.   In the cluster approach, the effects of the parts of the enzyme that are not explicitly included in the model are taken into account using implicit solvation methods.   For all six zinc-dependent enzymes studied, the di-zinc bridging hydroxide has been shown to be capable of performing nucleophilic attack on the substrate. In addition, one, two, or even all three zinc ions participate in the stabilization of the negative charge in the transition states and intermediates, thereby lowering the barriers.   For the two tungstoenzymes, several different mechanistic scenarios have been considered to identify the energetically most feasible one. For both enzymes, new mechanisms are proposed.   Finally, the mechanism of mycolic acid cyclopropane synthase has been shown to be a direct methyl transfer to the substrate double bond, followed by proton transfer to the bicarbonate.   From the studies of these enzymes, we demonstrate that density functional calculations are able to solve mechanistic problems related to enzymatic reactions, and a wealth of new insight can be obtained.

La present Tesi Doctoral, titulada desenvolupament computacional de la semblançamolecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblançabasades en la comparació de funcions de densitat electrònica.El primer capítol, Semblança quàntica, és introductori. S'hi descriuen les funcions de densitat de probabilitat electrònica i llur significança en el marc de la mecànica quàntica. Se n'expliciten els aspectes essencials i les condicions matemàtiques a satisfer, cara a una millor comprensió dels models de densitat electrònica que es proposen. Hom presenta les densitats electròniques, mencionant els teoremes de Hohenberg i Kohn i esquematitzant la teoria de Bader, com magnituds fonamentals en la descripció de les molècules i en la comprensió de llurspropietats.En el capítol Models de densitats electròniques moleculars es presenten procediments computacionals originals per l'ajust de funcions densitat a models expandits en termes de gaussianes 1s centrades en els nuclis. Les restriccions físico-matemàtiques associades a lesdistribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenatAtomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC,parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions ocapes de l'expansió, d'acord als requisits de no negativitat. La qualitat d'aquestes densitats i de les mesures de semblança derivades es verifica abastament. Aquest model ASA s'estén arepresentacions dinàmiques, físicament més acurades, en quant que afectades per les vibracions nuclears, cara a una exploració de l'efecte de l'esmorteïment dels pics nuclears en les mesures de semblança molecular. La comparació de les densitats dinàmiques respecte les estàtiques evidencia un reordenament en les densitats dinàmiques, d'acord al que constituiria una manifestació del Principi quàntic de Le Chatelier. El procediment ASA, explícitament consistent amb les condicions de N-representabilitat, s'aplica també a la determinació directe de densitats electròniqueshidrogenoides, en un context de teoria del funcional de la densitat.El capítol Maximització global de la funció de semblança presenta algorismes originals per la determinació de la màxima sobreposició de les densitats electròniques moleculars. Les mesures de semblança molecular quàntica s'identifiquen amb el màxim solapament, de manera es mesuri la distància entre les molècules, independentment dels sistemes de referència on es defineixen les densitats electròniques. Partint de la solució global en el límit de densitats infinitament compactades en els nuclis, es proposen tres nivells de aproximació per l'exploraciósistemàtica, no estocàstica, de la funció de semblança, possibilitant la identificació eficient del màxim global, així com també dels diferents màxims locals. Es proposa també una parametrització original de les integrals de recobriment a través d'ajustos a funcions lorentzianes, en quant que tècnica d'acceleració computacional. En la pràctica de les relacions estructura-activitat, aquests avenços possibiliten la implementació eficient de mesures de semblança quantitatives, i, paral·lelament, proporcionen una metodologia totalment automàtica d'alineació molecular. El capítol Semblances d'àtoms en molècules descriu un algorisme de comparació dels àtoms de Bader, o regions tridimensionals delimitades per superfícies de flux zero de la funció de densitat electrònica. El caràcter quantitatiu d'aquestes semblances possibilita la mesura rigorosa de la noció química de transferibilitat d'àtoms i grups funcionals. Les superfícies de flux zero i elsalgorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació mésacurada pel càlcul de les propietats atòmiques.Finalment, en el capítol Semblances en estructures cristal·lines hom proposa una definició original de semblança, específica per la comparació dels conceptes de suavitat o softness en la distribució de fonons associats a l'estructura cristal·lina. Aquests conceptes apareixen en estudis de superconductivitat a causa de la influència de les interaccions electró-fonó en les temperatures de transició a l'estat superconductor. En aplicar-se aquesta metodologia a l'anàlisi de sals de BEDT-TTF,s'evidencien correlacions estructurals entre sals superconductores i no superconductores, enconsonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinadesinteraccions.Conclouen aquesta tesi un apèndix que conté el programa ASAC, implementació de l'algorisme ASA, i un capítol final amb referències bibliogràfiques. / The present Doctoral Thesis, entitled Computational Development of Quantum Molecular Similarity, fundamentally deals on the calculation of similar-ity measures arising from the comparison of electron density functions.The first chapter, Quantum Similarity, is introductory. Electron probability functions are described, emphasizing their significance in Quantum Mechanics, and their mathematical constrains.In the chapter Models of molecular electron densities, original procedures to fit electron densities to 1s Gaussian expansions are presented. Mathematical constrains attached to probability distribution functions are explicitly considered,in the procedure named Atomic Shell Approximation (ASA). This procedure, im-plementedin the computer program ASAC, uses an initial, nearly complete func-tional space, from where functions or shells are variationally selected, according to the non-negativity requirements. The quality of these model densities and theaccuracy of the derived similarity measures are extensively verified. The ASA model is also extended to dynamic distributions, presumably a more physical rep-resentation of free molecule and ligand electron densities. The ASA procedure,explicitly consistent with the N-representability conditions, is adapted to the direct determination of hydrogenoid electron densities, in a context of the Density Func-tional Theory.The chapter Global Maximization of the Similarity Function describes original algorithms to determine the maximum overlap of two molecular electron densities. Similarity measures are identified with the maximum overlap in order tomeasure the distances among molecules, independently on the reference frame-.work where they are defined. Starting from the known global solution attached to hypothetical, infinitely compacted molecular electron densities, one proposes three levels of approach for an efficient scanning and global maximization of the non-deformed similarity function. Parametrazing overlap integrals through Lorentzian-likefunctions is also proposed to speed up computations. In the practice of structure-activity relationships, the presented advances provide an efficient implemen-tation of quantitative similarity measures, and, moreover, provide a new, com-pletely automatic methodology for molecular superposition and alignments.The chapter Similarities of atoms in molecules describes an algorithm for the comparison of Bader atoms. The accurate similarity measures obtained pro-vide a rigorous quantification of the degree of transferability of atoms and func-tional groups.Finally, in the chapter Similarities among crystalline structures, it is pro-posed a similarity definition for the comparison of crystalline structures regarding the concept of softness. This concept emerges from the BCS theory of supercon-ductivity.It appears related to the influence of electron-phonon interactions in the transition temperatures to the superconducting state. The application of this methodology in analyzing BEDT-TTF salts reveals a structural correlation among su-perconductors and non-superconductors, according to pointed hypothesis regard-ing the influence of some intermolecular interactions.The present Thesis concludes listing the ASAC code, implementation of the ASA algorithm, together with a chapter containing bibliographic references.

La present tesi està centrada en l'ús de la Teoria de Semblança Quàntica per a calcular descriptors moleculars. Aquests descriptors s'utilitzen com a paràmetres estructurals per a derivar correlacions entre l'estructura i la funció o activitat experimental per a un conjunt de compostos. Els estudis de Relacions Quantitatives Estructura-Activitat són d'especial interès per al disseny racional de molècules assistit per ordinador i, en particular, per al disseny de fàrmacs. Aquesta memòria consta de quatre parts diferenciades. En els dos primers blocs es revisen els fonaments de la teoria de semblança quàntica, així com l'aproximació topològica basada en la teoria de grafs. Ambdues teories es fan servir per a calcular els descriptors moleculars. En el segon bloc, s'ha de remarcar la programació i implementació de programari per a calcular els anomenats índexs topològics de semblança quàntica. La tercera secció detalla les bases de les Relacions Quantitatives Estructura-Activitat i, finalment, el darrer apartat recull els resultats d'aplicació obtinguts per a diferents sistemes biològics. / The present thesis is centred in the use of the Quantum Similarity Theory to calculate molecular descriptors. These molecular descriptors are used as structural parameters to derive correlations between the structure and the function or experimental activity for a set of compounds. Quantitative Structure-Activity Relationship studies are of special interest for the rational Computer-Aided Molecular Design and, in particular, for Computer-Aided Drug Design. The memory has been structured in four differenced parts. The two first blocks revise the foundations of quantum similarity theory, as well as the topological approximation, based in classical graph theory. These theories are used to calculate the molecular descriptors. In the second block, the programming and implementation of Topological Quantum Similarity Indices must be remarked. The third section details the basis for Quantitative Structure-Activity Relationships and, finally, the last section gathers the application results obtained for different biological systems.

Noncovalent interactions are of pivotal importance in many areas of chemistry, biology, and materials science, and the intermolecular interactions involving aromatic rings in particular, are fundamental to molecular organization and recognition processes. The work detailed in this thesis involves the application of state-of-the-art ab initio electronic structure theory methods to elucidate the nature of pi-pi interactions. The binding energies, and geometrical and orientational preferences of the simplest prototype of aromatic pi-pi interactions, the benzene dimer, are explored. We obtain the first converged values of the binding energies using highly accurate methods and large basis sets. Results from this study predict the T-shaped and parallel-displaced configurations of benzene dimer to be nearly isoenergetic. The role of substituents in tuning pi-pi interaction is investigated. By studying dimers of benzene with various monosubstituted benzenes (in the sandwich and two T-shaped configurations), we surprisingly find that all of the substituted sandwich dimers considered bind more strongly than benzene dimer. We also find that these interactions can be tuned by a modest degree of substitution. Energy decomposition analysis using symmetry-adapted perturbation theory (SAPT) reveals that models based solely on electrostatic effects will have difficulty in reliably predicting substituent effects in pi-pi interactions.

The field of organic electronics, in which -conjugated, organic molecules and polymers are used as the active components (e.g., semiconductor, light emitter/harvester, etc.), has lead to a number a number of key technological developments that have been founded within fundamental research disciplines. In the Dissertation that follows, the research involves the use of quantum-chemical techniques to elucidate fundamental aspects of both intermolecular and intramolecular electron-transfer processes in organic, -conjugated molecules. The Dissertation begins with an introduction and brief review of organic molecular systems used as electron-transport semiconducting materials in device applications and/or in the fundamental studies of intramolecular mixed-valence processes. This introductory material is then followed by a brief review of the electronic-structure methods (e.g., Hartree-Fock theory and Density Functional Theory) and electron-transfer theory (i.e., semiclassical Marcus theory) employed throughout the investigations. The next three Chapters deal with investigations related to the characterization of non-rigid, -conjugated molecular systems that have amorphous solid-state properties used as the electron-transport layer in organic electronic and optoelectronic devices. Chapters 3 and 4 involve studies of silole- (silacyclopentadiene)-based materials that possess attractive electronic and optical properties in the solid state. Chapter 5 offers a preliminary study of dioxaborine-based molecular structures as electron-transport systems. In Chapters 6 8, the focus of the work shifts to investigations of organic mixed-valence systems. Chapter 6 centers on the examination of tetraanisylarylenediamine systems where the inter-redox site distances are approximately equal throughout the series. Chapter 7 examines the bridge-length dependence of the geometric structure, charge-(de)localization, and electronic coupling for a series of vinylene- and phenylene-vinylene-bridged bis-dianisylamines. In Chapter 8, the role of symmetric vibrations in the delocalization of the excess charge is studied in a dioxaborine radical-anion and a series of radical-cation bridged-bisdimethylamines. Finally, Chapter 9 provides a synopsis of the work and goals for future consideration.

The structural, electronic, and optical properties of a series of organic pi-conjugated polymer, oligomer, or molecular materials of interest for applications in organic electronics are described. For this purpose, quantum-chemical techniques ranging from Density Functional Theory to Hartree-Fock ab initio and semiempirical methods are used to evaluate the charge-transport, charge-transfer, and electroluminescence properties of pi-conjugated organic materials. First, the effect of electronic polarization on the charge-transport parameters of organic semiconductors is discussed. A generalized methodology based on a basis set orthogonalization procedure is developed to determine reliable charge-transport characteristics. The charge-transport parameters of a number of organic semiconductors such as oligoacenes and derivatives are studied with this methodology. Then, triplet emitters, in particular iridium complexes, that achieve high efficiency electroluminescence in organic light-emitting diodes are discussed. The effects of ligand substitution and orientation on the luminescent properties of iridium compounds are investigated in order to develop structure-properties relationships. The emission properties of these complexes are found to be governed by an interplay between metal-to-ligand charge transfer excitations and ligand-centered and/or interligand excitations. The extent of mixing of these various excitations turns out to be highly dependent on the nature of the substituents. Design strategies to shift the emission color towards deep blue are proposed. Finally, several classes of materials acting as hosts for phosphorescent emitters are studied. It is shown that restricting the conjugation length leads to high energy gap host materials suitable for blue phosphorescent emitters.

This thesis is devoted to first-principles simulations of bio- and nano-materials,focusing on various soft x-ray spectra, ground-state energies and structures of isolated largemolecules, bulk materials, and small molecules in ambient solutions. K-edge near-edge x-ray absorption fine structure (NEXAFS) spectra, x-ray emission spectra, andresonant inelastic x-ray scattering spectra of DNA duplexes have been studied by means oftheoretical calculations at the density functional theory level. By comparing a sequence of DNAduplexes with increasing length, we have found that the stacking effect of base pairs has verysmall influence on all kinds of spectra, and suggested that the spectra of a general DNA can bewell reproduced by linear combinations of composed base pairs weighted by their ratio. The NEXAFS spectra study has been extended to other realistic systems. We have used cluster modelswith increasing sizes to represent the infinite crystals of nucleobases and nucleosides, infinitegraphene sheet, as well as a short peptide in water solution. And the equivalent core holeapproximation has been extensively adopted, which provides an efficient access to these largesystems. We have investigated the influence of external perturbations on the nitrogen NEXAFSspectra of guanine, cytosine, and guanosine crystals, and clarified early discrepancies betweenexperimental and calculated spectra. The effects of size, stacking, edge, and defects to theabsorption spectra of graphene have been systematically analyzed, and the debate on theinterpretation of the new feature has been resolved. We have illustrated the influence of watersolvent to a blocked alanine molecule by using the snapshots generated from molecular dynamics. Multi-scale computational study on four short peptides in a self-assembled cage is presented. It isshown that the conformation of a peptide within the cage does not corresponds to its lowest-energyconformation in vacuum, due to the Zn-O bond formed between the peptide and the cage, and theconfinement effect of the cage. Special emphasis has been paid on a linear-scaling method, the generalized energy basedfragmentation energy (GEBF) approach. We have derived the GEBF energy equation at the Hartree-Focklevel with the Born approximation of the electrostatic potential. Numerical calculations for amodel system have explained the accuracy of the GEBF equation and provides a starting point forfurther refinements. We have also presented an automatic and efficient implementation of the GEBFapproach which is applicable for general large molecules. / QC 20110404

The main focus of this thesis is ab initio modeling of core level spectra with a high-level quantum chemical description both of the chemical interactions and of local atomic multiplet effects. In particular, the combination of calculations and synchrotron-based core-level spectroscopy aims at understanding the local structure of the electronic valence in transition metal complexes, and the details of the solvation mechanisms of electrolyte solutions, systems relevant for the dye-sensitized solar cell. Configurational sampling in solution is included through molecular dynamics simulations. Transition metal complexes are studied with x-ray absorption (XA) and resonant inelastic scattering (RIXS) spectroscopy, characterizing excited states with atomic site specificity. The theoretical multiconfigurational method, applying an active-space partitioning of the molecular orbitals (RASSCF), is used to assign the transitions observed in spectra of hydrated Ni2+ explicitly, including charge transfer and multiplet effects. Furthermore, the solvent-induced binding energy properties of the I- and I3- anions in aqueous, ethanol, and acetonitrile solutions are analyzed using photoelectron spectroscopy (XPS). The study shows that specific ion–solvent interactions are important for the core-level binding energy shifts in solution. The special case with I3- dissolved in water, where hydrogen bonding causes breaking of the molecular symmetry, is treated and proves that the geometry changes influence the photoelectron spectrum of aqueous I3- directly.

At the atomic scale, interactive physically-based modeling tools are more and more in demand. Unfortunately, solving the underlying physics equations at interactive rates is computationally challenging. In this dissertation, we propose new algorithms that allow for interactive modeling of chemical structures. We first present a modeling tool to construct structural models of hydrocarbon systems. The physically-based feedbacks are based on the Brenner potential. In order to be able to interactively edit systems containing numerous atoms, we introduce a new adaptive simulation algorithm. Then, we introduce what we believe to be the first interactive quantum chemistry simulation algorithm at the Atom Superposition and Electron Delocalization Molecular Orbital (ASED-MO) level of theory. This method is based on the divide-and-conquer (D&C) approach, which we show is accurate and efficient for this non-self-consistent semi-empirical theory. We then propose a novel Block-Adaptive Quantum Mechanics (BAQM) approach to interactive quantum chemistry. BAQM constrains some nuclei positions and some electronic degrees of freedom on the fly to simplify the simulation. Finally, we demonstrate several applications, including one study of graphane formation, interactive simulation for education purposes, and virtual prototyping at the atomic scale, both on desktop computers and in virtual reality environments.

The nature of the ground state and the first (lowest) singlet excited state geometrical conformations and electronic transitions in the aromatic five-membered heterocyclic oligomers &ndash / oligothiophenes (nT), oligofurans (nF), and oligopyrroles (nP)- containing up to six monomer units (total of 18 molecules) were explored using several computational methodologies. Geometry optimizations were carried out at Austin Model 1 (AM1), Restricted Hartree-Fock (RHF/6-31G*), and Density Functional Theory (DFT, B3LYP/6-31G*) levels for the ground-state conformations of these structurally well-defined heterocyclic oligomers. The Configuration Interaction Singles (CIS) method with the 6-31G* basis set was chosen in computation of the optimal geometry of the lowest singlet excited state. Lowest singlet excitation S1&szlig / S0 energies were calculated using the Zerner&rsquo / s Intermediate Neglect of Differential Overlap for Spectroscopy (ZINDO/S), CIS (CIS/6-31G*), and Time-Dependent DFT (TDDFT/6-31G* and TDDFT/6-31+G*) methods. In computation of the emission S1&agrave / S0 energies, we have employed all methods above except ZINDO/S. In investigation of geometries of the ground and lowest singlet excited state, we compared the bond length alternation (BLA) parameters, Dri in the conjugated backbone of the oligomers. Saturation of the geometrical parameters at the center of oligomers was observed after a certain chain length. Among all methodologies used in computation of excitation (S1&szlig / S0) and emission (S1&agrave / S0) energies, TDDFT results showed the best agreement with experimental data. Fits of computed and experimental excitation energies to an exponential function using the least squares method enabled us to predict Effective Conjugation Length (ECL) values. We obtained the ECLs of 17 (17), 16 (15), and 14 (13) monomer units for polythiophene (PTh), polyfuran (PFu), and polypyrrole (PPr), which have very good agreement with the results obtained from the fits of experimental data (the values in parentheses).