Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

Síntesi i estructura de lligands hexaazamacrocíclics i dels seus complexos de Zn. Estudis de catàlisi i de fenòmens de reconeixement molecular amb substrats aniònics

Es realitza una breu introducció a la Química Supramolecular, la qual estudia les estructures i funcions de les associacions que resulten de la unió d'espècies moleculars a través d'enllaços intermoleculars, no covalents. El reconeixement molecular, reactivitat i transport són les funcions bàsiques de les espècies supramoleculars. El present treball tracta de la química supramolecular d'anions i del potential efecte catalític de complexos de Zinc amb lligands hexaazamacrocilics en reaccions de hidròlisis d'esters de fosfat. Els lligands descrits en aquest treball són lligands macrocíclics hexadentats formats per amines secundàries o terciàries enllaçades a cadenes alquíliques de dos o tres àtoms de carbonis formant dos braços, que estan unides per un espaiador, El conjunt de lligands que s'estudien permeten estudiar els diferents factors, geomètrics, electrònics i estèrics que controlen els fenòmens de reconeixement molecular. El treball descriu en detall els procediments de síntesi i la caracterització d'aquest conjunt de lligands. S'ha realitzat un estudi sistemàtic dels diferents factors que afecten als fenòmens de reconeixement molecular entre lligands macrocíclics (L), descrits anteriorment i substrats aniònics (S): fosfats, polifosfats i dicarboxilats. Els estudis s'han realitzat tant en dissolució aquosa, a partir de mesures potenciomètriques, com en estat sòlid per difracció de Raig-X a partir de la obtenció de les estructures cristal·lines dels complexos. La força de l'enllaç en els complexos ternaris H:L:S es racionalitza en termes d'enllaç per pont d'hidrogen, interaccions electrostàtiques i interaccions per -Stacking. Així mateix es racionalitza la importància de les dimensions i forma de la cavitat en el grau d'interacció. Actualment hi ha un gran interès en l'estudi del paper que juguen els ions metàl·lics en el centre actiu dels metal·loenzims hidrolítics tal com la carboxipeptidasa, que catalitza la hidròlisi d'aminoàcids C-terminals de substrats polipeptídics, l'anhidrasa carbònica, que catalitza la reacció d'hidratació del CO2 per donar hidrogencarbonat, la fosfatasa alcalina, que actua catalitzant la hidròlisi no específica de monoesters de fosfats a pH alcalí. El Zn(II) és un metall que es troba sovint formant part dels centres actius dels metal·loenzims i és el segon element de transició més abundant en els organismes vius. Per aquest motiu són útils els models sintètics que puguin actuar com a mimetitzadors de reaccions bioinorgàniques en les que els metal·loenzims juguen un paper com a catalitzadors. Els models sintètics són complexos de metalls de transició formats per un lligand o lligands en el seu entorn de coordinació més immediat que tinguin efectes electrònics i estèrics els més similar possible als que es produeixen en el centre actiu dels enzims. Donat la importància dels complexos de Zn(II) com a models biomimètics de metal·loenzims hidrolítics, s'han sintetitzat i caracteritzat complexos dinuclears de Zn(II) amb lligands hexaazamacrocíclics descrits anteriorment. S'ha realitzat un estudi de la seva activitat catalítica en la hidròlisi d'un ester activat, l'acetat de p-nitrofenil, un substrat característic que s'utilitza per analitzar l'activitat catalítica dels models biomimètics.

A chemical Hamiltonian approach study of the basis set superposition error changes on electron densities and one- and two-center energy components

A comparision of the local effects of the basis set superposition error (BSSE) on the electron densities and energy components of three representative H-bonded complexes was carried out. The electron densities were obtained with Hartee-Fock and density functional theory versions of the chemical Hamiltonian approach (CHA) methodology. It was shown that the effects of the BSSE were common for all complexes studied. The electron density difference maps and the chemical energy component analysis (CECA) analysis confirmed that the local effects of the BSSE were different when diffuse functions were present in the calculations

A new all-round density functional based on spin states and SN2 barriers

We report here a new empirical density functional that is constructed based on the performance of OPBE and PBE for spin states and SN 2 reaction barriers and how these are affected by different regions of the reduced gradient expansion. In a previous study [Swart, Sol̀, and Bickelhaupt, J. Comput. Methods Sci. Eng. 9, 69 (2009)] we already reported how, by switching between OPBE and PBE, one could obtain both the good performance of OPBE for spin states and reaction barriers and that of PBE for weak interactions within one and the same (SSB-sw) functional. Here we fine tuned this functional and include a portion of the KT functional and Grimme's dispersion correction to account for π- π stacking. Our new SSB-D functional is found to be a clear improvement and functions very well for biological applications (hydrogen bonding, π -π stacking, spin-state splittings, accuracy of geometries, reaction barriers)

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Destabilizing Mutations Alter the Hydrogen Exchange Mechanism in Ribonuclease A

The effect of strongly destabilizing mutations, I106A and V108G of Ribonuclease A (RNase A), on its structure and stability has been determined by NMR. The solution structures of these variants are essentially equivalent to RNase A. The exchange rates of the most protected amide protons in RNase A (35ºC), the I106A variant (35ºC), and the V108G variant (10ºC) yield stability values of 9.9, 6.0, and 6.8 kcal/mol, respectively, when analyzed assuming an EX2 exchange mechanism. Thus, the destabilization induced by these mutations is propagated throughout the protein. Simulation of RNase A hydrogen exchange indicates that the most protected protons in RNase A and the V108G variant exchange via the EX2 regime, whereas those of I106A exchange through a mixed EX1 1 EX2 process. It is striking that a single point mutation can alter the overall exchange mechanism. Thus, destabilizing mutations joins high temperatures, high pH and the presence of denaturating agents as a factor that induces EX1 exchange in proteins. The calculations also indicate a shift from the EX2 to the EX1 mechanism for less protected groups within the same protein. This should be borne in mind when interpreting exchange data as a measure of local stability in less protected regions

Analysis of chemical bonding and aromaticity from electronic delocalization descriptors

Interactions between electrons determine the structure and properties of matter from molecules to solids. Therefore, the understanding of the electronic structure of molecules will enable us to extract relevant chemical information. In the first part of this thesis, we focus our attention on the analysis of chemical bonding by means of the Electron Localization Function (ELF) and the Domain-Averaged Fermi Hole analysis (DAFH). In the second part, we assess the performance of some indicators of aromaticity by analyzing their advantages and drawbacks. We propose a series of tests based on well-known aromaticity trends that can be applied to evaluate the aromaticity of current and future indicators of aromaticity in both organic and inorganic species. Moreover, we investigate the nature of electron delocalization in both aromatic and antiaromatic systems in the light of Hückel’s (4n + 2) rule. Finally, we analyze the phenomenon of multiple aromaticity in all-metal clusters.