Desenvolupament computacional de la semblança molecular quàntica

La present Tesi Doctoral, titulada desenvolupament computacional de la semblança molecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblança basades 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 llurs propietats.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 les distribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenat Atomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC, parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions o capes 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 a representacions 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òniques hidrogenoides, 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 els algorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació més acurada 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, en consonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinades interaccions.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.

Development, implementation and application of electronic structural descriptors to the analysis of the chemical bonding, aromaticity and chemical reactivity

En la literatura sobre mecànica quàntica és freqüent trobar descriptors basats en la densitat de parells o la densitat electrònica, amb un èxit divers segons les aplicacions que atenyin. Per tal de que tingui sentit químic un descriptor ha de donar la definició d'un àtom en una molècula, o ésser capaç d'identificar regions de l'espai molecular associades amb algun concepte químic (com pot ser un parell solitari o zona d'enllaç, entre d'altres). En aquesta línia, s'han proposat diversos esquemes de partició: la teoria d'àtoms en molècules (AIM), la funció de localització electrònica (ELF), les cel·les de Voroni, els àtoms de Hirshfeld, els àtoms difusos, etc. L'objectiu d'aquesta tesi és explorar descriptors de la densitat basats en particions de l'espai molecular del tipus AIM, ELF o àtoms difusos, analitzar els descriptors existents amb diferents nivells de teoria, proposar nous descriptors d'aromaticitat, així com estudiar l'habilitat de totes aquestes eines per discernir entre diferents mecanismes de reacció.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

A theoretical study on the XeF(2) molecule

A relativistic four-component study was performed for the XeF(2) molecule by using the Dirac-Coulomb (DC) Hamiltonian and the relativistic adapted Gaussian basis sets (RAGBSs). The comparison of bond lengths obtained showed that relativistic effects on this property are small (increase of only 0.01 angstrom) while the contribution of electron correlation, obtained at CCSD(T) or CCSD-T levels, is more important (increase of 0.05 angstrom). Electron correlation is also dominant over relativistic effects for dissociation energies. Moreover, the correlation-relativity interaction is shown to be negligible for these properties. The electron affinity, the first ionization potential and the double ionization potential are obtained by means of the Fock-space coupled cluster (FSCC) method, resulting in DC-CCSD-T values of 0.3 eV, 12.5 eV and 32.3 eV, respectively. Vibrational frequencies and some anharmonicity constants were also calculated under the four-component formalism by means of standard perturbation equations. All these molecular properties are, in general, ill satisfactory agreement with available experimental results. Finally, a partition in terms of charge-charge flux-dipole flux (CCFDF) contributions derived by means of the quantum theory of atoms in molecules (QTAIM) in non-relativistic QCISD(FC)/3-21G* calculations was carried out for XeF(2) and KrF(2). This analysis showed that the most remarkable difference between both molecules lies on the charge flux contribution to the asymmetric stretching mode, which is negligible in KrF(2) but important in XeF(2). (c) 2008 Elsevier B.V. All rights reserved.

A theoretical study on the gas phase reactions of the anions NbO3-, NbO5-, and NbO2(0H)(2)(-) with H2O and O-2

The potential energy surfaces at the singlet (s) and the triplet (t) electronic states associated with the gas-phase ion/molecule reactions of NbO3-, NbO5-, and NbO2(OH)(2)(-) with H2O and O-2 have been investigated by means of DFT calculations at the B3LYP level. An analysis of the results points out that the most favorable reactive channel comprises s-NbO3- reacting with H2O to give an ion-molecule complex s-NbO3(H2O)without a barrier. From this minima, an intramolecular hydrogen transfer takes place between the incoming water molecule and an oxygen atom of the NbO3- fragment to render the most stable minimum, s-NbO2(OH)(2)(-). This oxyhydroxide system reacts with O-2 along a barrierless process to obtain the triplet t-NbO4(OH)(2)(-)-A intermediate, and the crossing point, CP1, between s and t electronic states has been characterized. The next step is the hydrogen-transfer process between the oxygen atom of a hydroxyl group and the one adjacent oxygen atom to render a minimum with the two OH groups near each other, t-NbO4(OH)(2)(-)-B. From this point, the last hydrogen migration takes place, to obtain the product complex, t-NbO5(H2O)(-), that can be connected with the singlet separated products, s-NbO5- and H2O. Therefore, a second crossing point, CP2, has been localized. The nature of the chemical bonding of the key minima (NbO3-, NbO2(OH)(2)(-), NbO4(OH)(2)(-)-B, and NbO5-) in both electronic states of the reaction and an interaction with O-2 has been studied by topological analysis of Becke-Edgecombe electron-localization function (ELF) and atoms-in-molecules (AIM) methodology. The niobium-oxygen interactions are characterized as unshared-electron (ionic) interactions and some oxygen-oxygen interactions as protocovalent bonds.

Estudo teórico QTAIM e DFT dos compostos de coordenação: efeito quelato, titanocenos e ligação química

This work is a study of coordination compounds by quantum theory of atoms in molecules (QTAIM), based on the topological analysis of the electron density of molecular systems, both theoretically and experimentally obtained. The coordination chemistry topics which were studied are the chelate effect, bent titanocene and chemical bond in coordination complexes. The chelate effect was investigated according to topological and thermodynamic parameters. The exchange of monodentate ligands on polydentate ligands from same transition metal increases the stability of the complex both from entropy and enthalpy contributions. In some cases, the latter had a higher contribution to the stability of the complex in comparison with entropy. This enthalpic contribution is explained according to topological analysis of the M-ligand bonds where polidentate complex had higher values of electron density of bond critical point, Laplacian of electron density of bond critical point and delocalization index (number of shared electrons between two atoms). In the second chapter, was studied bent titanocenes with bulky cyclopentadienyl derivative π-ligand. The topological study showed the presence of secondary interactions between the atoms of π-ligands or between atoms of π-ligand and -ligand. It was found that, in the case of titanocenes with small difference in point group symmetry and with bulky ligands, there was an nearly linear relationship between stability and delocalization index involving the ring carbon atoms (Cp) and the titanium. However, the titanocene stability is not only related to the interaction between Ti and C atoms of Cp ring, but secondary interactions also play important role on the stability of voluminous titanocenes. The third chapter deals with the chemical bond in coordination compounds by means of QTAIM. The quantum theory of atoms in molecules so far classifies bonds and chemical interactions in two categories: closed shell interaction (ionic bond, hydrogen bond, van der Waals interaction, etc) and shared interaction (covalent bond). Based on topological parameters such as electron density, Laplacian of electron density, delocalization index, among others, was classified the chemical bond in coordination compounds as an intermediate between closed shell and shared interactions

Estudo teórico aplicado ao petróleo: estabilidade dos alcenos e recuperação terciária

The present work aims to study the theoretical level of some processes employed in the refining of petroleum fractions and tertiary recovery of this fluid. In the third chapter, we investigate a method of hydrogenation of oil fractions by QTAIM (Quantum Theory of Atoms in Molecules) and thermodynamic parameters. The study of hydrogenation reactions, and the stability of the products formed, is directly related to product improvement in the petrochemical refining. In the fourth chapter, we study the theoretical level of intermolecular interactions that occur in the process of tertiary oil recovery, or competitive interactions involving molecules of non-ionic surfactants, oil and quartz rock where oil is accumulated. Calculations were developed using the semiempirical PM3 method (Parametric Model 3). We studied a set of ten non-ionic surfactants, natural and synthetic origin. The study of rock-surfactant interactions was performed on the surface of the quartz (001) completely hydroxylated. Results were obtained energetic and geometric orientations of various surfactants on quartz. QTAIM was obtained through the analysis of the electron density of interactions, and thus, providing details about the formation of hydrogen bonds and hydrogen-hydrogen systems studied. The results show that the adsorption of ethoxylated surfactants in the rock surface occurs through the hydrogen bonding of the type CH---O, and surfactants derivatives of polyols occurs by OH---O bonds. For structures adsorption studied, the large distance of the surfactant to the surface together with the low values of charge density, indicate that there is a very low interaction, characterizing physical adsorption in all surfactants studied. We demonstrated that surfactants with polar group comprising oxyethylene units, showed the lowest adsorption onto the surface of quartz, unlike the derivatives of polyols

Índice de aromaticidade baseado na deslocalização, degenerescência e densidade (D3BIA) para acenos e moléculas homoaromáticas

The aromaticity index is an important tool for the investigation of aromatic molecules. This work consists on new applications of the aromaticity index developed by teacher Caio Lima Firme, so-called D3BIA (density, delocalization, degeneracy-based index of aromaticity). It was investigated its correlation with other well-known aromaticity indexes, such as HOMA (harmonic oscillator model of aromaticity), NICS (nucleus independent chemical shielding), PDI (para-delocalization index), magnetic susceptibility (), and energetic factor in the study of aromaticity of acenes and homoaromatic species based on bisnoradamantanyl cage. The density functional theory (DFT) was used for optimization calculations and for obtaining energetic factors associated with aromaticity and indexes HOMA and NICS. From quantum theory of atoms in molecules (QTAIM) it was obtained the indexes D3BIA, PDI and . For acenes, when the over-mentioned indexes were applied it was observed no correlation except for D3BIA and HOMA (R2=0.752). For bisnoradamantenyl dication and its derivatives, it was obtained a good correlation between D3BIA and NICS. Moreover, it was evaluated solely one of the factors used on D3BIA calculation, the delocalization index uniformity (DIU), so as to investigate its possible influence on stability of chemical species. Then, the DIU was compared with the formation Gibbs free energy of some pairs of carbocations, isomers or not, which each pair had small difference in point group symmetry and no difference among other well-known stability factors. The obtained results indicate that DIU is a new stability factor related to carbocations, that is, the more uniform the electron density delocalization, the more stable the is carbocation. The results of this work validate D3BIA and show its importance on the concept of aromaticity, indicating that it can be understood from degeneracy of atoms belonging the aromatic site, the electronic density in the aromatic site and the degree of uniformity of electron delocalization

Electronic structure of polyacetylene with bonded oxygen

In this work we have studied the electronic structure of finite polyacetylene chains with structural oxygen-bonding models following data from 13C-NMR experiments. We have used a combination of Austin Method One and Hydrogenic Atoms in Molecules version 3 methods to perform geometric and spectroscopic calculations. Our results show that the electronically-active states are generally unaffected by the incorporation of oxygen. © 1992.

Structural, electronic and optical properties of Fe(III) complex with pyridine-2,6-dicarboxylic acid: A combined experimental and theoretical study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A combined experimental and theoretical study on the formation of ag filaments on beta-Ag2MoO4 induced by electron irradiation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Isatin-Schiff base copper(II) complexes-A DFT study of the metal-ligand bonding situation

Herein, we report results of calculations based on density functional theory (BP86/TZVP) of a set of isatin-Schiff base copper(II) and related complexes, 1-12, that have shown significant pro-apoptotic activity toward diverse tumor cells. The interaction of the copper(II) cation with different ligands has been investigated at the same level of theory. The strength and character of the Cu(II)-L bonding was characterized by metal-ligand bond lengths, vibrational frequencies, binding energies, ligand deformation energies, and natural population analysis. The metal-ligand bonding situation was also characterized by using two complementary topological approaches, the quantum theory of atoms-in-molecules (QTAIM) and the electron localization function (ELF). The calculated electronic g-tensor and hyperfine coupling constants present significant agreement with the EPR experimental data. The calculated parameters pointed to complex 10 as the most stable among the isatin-Schiff base copper(II) species, in good agreement with experimental data that indicate this complex as the most reactive in the series. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

The infrared fundamental intensities of some cyanopolyynes

Some cyanopolyynes, HCnN (n = 1, 3, ... , 17), are investigated by means of calculations at the MP2/cc-pVTZ and CCSD/cc-pVDZ levels. Although the MP2/cc-pVTZ results for geometries and molecular dipole moments are encouraging, the CCSD/cc-pVDZ level was superior for the study of infrared fundamental intensities. The main bands are also analyzed with a charge-charge flux-dipole flux (CCFDF) partition model based on quantities given by the Quantum Theory of Atoms in Molecules (QTAIM). The intensity of vibrations corresponding to the stretching of CH bonds (3471-3473 cm(-1)) increases in line with the number of carbon atoms (from 61 to 146 km mol(-1) between HCN and HC13N). This increase is due to the charge flux contribution while the other contributions remain roughly unaltered except for HCN. Moreover, the hydrogen atom loses an almost constant amount of electronic charge during the CH bond enlargement and a small fraction of this charge spreads to atoms farther and farther away from hydrogen as the molecule size increases. The band associated with the doubly degenerate CH bending vibrations (643-732 cm(-1)) presents approximately the same intensity in all the studied cyanopolyynes (from 67 to 76 km mol(-1)). The CCFDF/QTAIM contributions are also nearly the same for these bending modes in HC5N and larger systems. The intensity of the mode mostly identified as CN stretching (around 2378-2399 cm(-1) except for HCN) increases from HCN up to HC7N (from 0.3 to 83 km mol(-1)) and nearly stabilizes around 80-90 km mol(-1) for larger systems. The CCFDF/QTAIM contributions for this mode also change significantly up to HC7N and remain almost constant in larger systems. We also observed the appearing of a very relevant band between 2283 and 2342 cm(-1). This mode is mainly associated with the symmetric stretching of CC triple bonds near the molecule center and exhibits large charge fluxes while the other contributions are almost negligible in the largest cyanopolyynes. The two vibrational bands associated with the smallest frequencies are also studied and extrapolation equations are suggested to predict their positions in larger cyanopolyynes. (C) 2012 Elsevier B.V. All rights reserved.

Gas-phase reactivity of 2-hydroxy-1,4-naphthoquinones: a computational and mass spectrometry study of lapachol congeners

In order to understand the influence of alkyl side chains on the gas-phase reactivity of 1,4-naphthoquinone derivatives, some 2-hydroxy-1,4-naphthoquinone derivatives have been prepared and studied by electrospray ionization tandem mass spectrometry in combination with computational quantum chemistry calculations. Protonation and deprotonation sites were suggested on the basis of gas-phase basicity, proton affinity, gas-phase acidity (?Gacid), atomic charges and frontier orbital analyses. The nature of the intramolecular interaction as well as of the hydrogen bond in the systems was investigated by the atoms-in-molecules theory and the natural bond orbital analysis. The results were compared with data published for lapachol (2-hydroxy-3-(3-methyl-2-butenyl)-1,4-naphthoquinone). For the protonated molecules, water elimination was verified to occur at lower proportion when compared with side chain elimination, as evidenced in earlier studies on lapachol. The side chain at position C(3) was found to play important roles in the fragmentation mechanisms of these compounds. Copyright (c) 2012 John Wiley & Sons, Ltd.

Modern charge density studies: the entanglement of experiment and theory

This tutorial review article is intended to provide a general guidance to a reader interested to learn about the methodologies to obtain accurate electron density mapping in molecules and crystalline solids, from theory or from experiment, and to carry out a sensible interpretation of the results, for chemical, biochemical or materials science applications. The review mainly focuses on X-ray diffraction techniques and refinement of experimental models, in particular multipolar models. Neutron diffraction, which was widely used in the past to fix accurate positions of atoms, is now used for more specific purposes. The review illustrates three principal analyses of the experimental or theoretical electron density, based on quantum chemical, semi-empirical or empirical interpretation schemes, such as the quantum theory of atoms in molecules, the semi-classical evaluation of interaction energies and the Hirshfeld analysis. In particular, it is shown that a simple topological analysis based on a partition of the electron density cannot alone reveal the whole nature of chemical bonding. More information based on the pair density is necessary. A connection between quantum mechanics and observable quantities is given in order to provide the physical grounds to explain the observations and to justify the interpretations.