No description available.

A scheme, analogous to the linear combination of atomic orbitals (LCAO), is used to calculate rates of reactions for the fusion of nuclei con¯ned in molecules. As an example, the possibility of nuclear fusion in rotationally excited H2O molecules of angular momentum 1¡ is estimated for the p + p + 16O ! 18Ne¤(4:522; 1¡) nuclear transition. Due to a practically exact agreement of the energy of the Ne resonance and of the p + p + 16O threshold, the possibility of an enhanced transition probability is investigated. / Physics / M.Sc.

Orientador: Rogério Custodio / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Química / Made available in DSpace on 2018-08-21T05:01:18Z (GMT). No. of bitstreams: 1 Streit_Livia_D.pdf: 1852639 bytes, checksum: 1d6bd6190d21b329d992dfe5f74f8993 (MD5) Previous issue date: 2012 / Resumo: O estudo de processos fotoeletrônicos, como energias de dupla ionização é de grante interesse teórico e experimental na área da física molecular. Cálculos e medidas experimentais de energias de dupla ionização e de energias de ligação de elétron envolvem o entendimento de correlação e relaxação eletrônica. A consideração de tais efeitos é indispensável para a obtenção de resultados precisos. No presente trabalho, tais processos fotoeletrônicos são estudados com métodos Pós-Hartree-Fock e Monte Carlo Quântico (MCQ). No último caso, o método foi adaptado de maneira a utilizar um tipo diferente de função de onda. A primeira parte desta tese envolve o estudo de energias de dupla ionização. O objetivo prinicipal é a análise do tipo de função de onda tentativa no método MCQ. Três diferentes formas foram estudadas: Uma função Hartree-Fock monodeterminante com uma função de correlação eletrônica explícita de Boys-Handy otimizada para a molécula; uma função CI multideterminante; e ainda, uma função de onda multideterminante construída a partir de orbitais de Dyson. Para obter a função de onda de Dyson, um extenso trabalho envolvendo estudo e programação do programa de cálculo de estrutura eletrônica Gaussian 09 foi necessário. A primeira parte da tese ainda envolve o estudo de energias de dupla ionização por métodos Interação de Configurações Multireferência MRCI e Gaussian-3 G3. Os resultados obtidos são comparados aos resultados obtidos por MCQ. A segunda parte da tese compreende o estudo de energias de ligação de elétron de complexos de gases nobres com o ânion Uracila. A coexistência de ânions em que o orbital singularmente ocupado apresenta características de um orbital de valência e de um orbital muito difuso é investigada através de análise de geometria e do cálculo de energias de afinidade eletrônica e de ligação de elétron dos complexos neutros e aniônicos / Abstract: The study of photoelectron processes such as double ionization energies has become a very interesting field in experimental and theoretical molecular physics. Calculation and experimental measure of double ionization energies and electron binding energies involve understanding of electron correlation and relaxation energies. Consideration of those effects are a requisite in order to predict precise values. In this work, those photoelectron processes are studied with post-Hartree-Fock and Quantum Monte Carlo (QMC) methods. In the latter case, the method was also modified to take in account a different type of wavefunction. The first part of this thesis is concerned with the study of double ionization energies. The main objective of this work is the analysis of the type of trial wavefunction in QMC method. Three different forms were studied: A single determinant Hartree-Fock function, with the Boys-Handy explicit electron correlation function optimized for the molecule; a multi-determinant CI wavefunction; and finally, a multi-determinant wavefunction built from Dyson orbitals. In order to obtain the Dyson wavefunction, an extensive work on studying and programming the electron structure package Gaussian 09 was necessary. The first part still comprises the study of double ionization energies with MRCI and G3 methods. The results thus obtained were compared to the QMC results. The second part of this thesis embraces the study of electron binding energies of uracil¿noble-gas complexes. The coexistence of valence and diffuse¿bound anions of U(Xe), suggested by experimental studies is investigated with geometry analysis and calculation of adiabatic electron affinities and vertical electron detachment energies of the neutral an anionic complexes / Doutorado / Físico-Química / Doutor em Ciências

Orientador: Rogerio Custodio / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Quimica / Made available in DSpace on 2018-08-14T14:18:23Z (GMT). No. of bitstreams: 1 Paschoal_JulianadeLima_M.pdf: 914971 bytes, checksum: 5a3a529f3c0006227e418c9e1629e6a5 (MD5) Previous issue date: 2006 / Resumo: Uma estratégia recente denominada Monte Carlo Quântico (MCQ) permite acessar a função de onda exata de um sistema resolvendo a equação de Schrödinger. Dentre as alternativas de MCQ destacam-se o Monte Carlo Quântico Variacional (MCQV) e o Monte Carlo Quântico de Difusão (MCQD). O MCQV determina o valor médio de qualquer propriedade atômica ou molecular associada a uma função de onda arbitrária empregando o algoritmo de Metropolis. O MCQD, por sua vez, baseia-se na solução da equação de Schrödinger dependente do tempo através de um processo de difusão em equilíbrio com um processo cinético de primeira ordem. Neste trabalho os objetivos são: a) comparar os efeitos de funções de base de Slater com diferentes expoentes nos níveis de teoria do MCQV e MCQD; b) testar funções de onda baseadas no modelo Hartree e Hartree-Fock no MCQV e MCQD e c) avaliar o efeito da localização de orbitais nestes métodos. Esses objetivos foram avaliados em átomos, moléculas diatômicas e alguns hidretos poliatômicos contendo elementos do segundo período da tabela periódica. Inicialmente, usou-se de uma função de onda representada por um determinante de Slater com orbitais obtidos através da combinação linear de funções de Slater através do método Hartree-Fock. Os expoentes do conjunto de base utilizado foram determinados através das Regras de Slater, otimização Hartree-Fock em ambiente molecular e otimizaçãoHartree-Fock em ambiente atômico. O MCQV e o MCQD foram empregados para a obtenção da energia média do sistema. Posteriormente, substituíram-se as funções de Slater por funções STO-6G. Os mesmos expoentes do conjunto de base utilizados nos cálculos com funções de Slater foram empregados para os cálculos STO-6G. Finalmente, utilizou-se o produto de Hartree como função de onda para os cálculos MCQV e MCQD com as funções de base já mencionadas. As principais conclusões desse trabalho são: a) o MCQD, conforme esperado, apresenta menores energias quando comparado ao MCQV; b) Cálculos MCQD usando determinante de Slater, conjunto de base com otimização de expoente para a molécula e átomo e nehum fator de correção forneceu energias comparadas a Gaussianas do tipo double-zeta no método coupled cluster incluindo excitações simples e duplas; c) Funções de base STO-6G devem ser utilizadas com cautela para representar funções STO; d) as energias calculadas através do produto de Hartree apresentam um comportamento que se distancia das funções Hartree-Fock quando orbitais localizados são usados; e) resultados melhores são esperados quando orbitais são auto-consistentes com respeito ao método de Hartree. / Abstract: A recent strategy called Quantum Monte Carlo (QMC) allows to access the exact wave function of a system solving Schrödinger¿s equation. Among the alternatives of QMC, Variational Quantum Monte Carlo (VQMC) and Diffusion Quantum Monte Carlo (DQMC) are distinguished. VQMC determines the average value of any atomic or molecular property associated to an arbitrary wave function using Metropolis algorithm. DQMC, on the other hand, is based on the solution of the time-dependent Schrödinger equation from a diffusion process in equilibrium with a first-order kinetic process. In this work the objectives were: a) to compare the effect of the Slater basis set with exponents adjusted in different environments at the VQMC and DQMC levels of theory; b) to test wave functions based on the Hartree and Hartree-Fock models along with VQMC and DQMC; c) to evaluate the effect of orbital localization within these methods. These objectives are evaluated in atoms, diatomic molecules and some polyatomic hydrates containing elements from the second period of the Periodic Table. Initially, a conventional wave function represented by a single Slater determinant is used with orbitals from the linear combination of Slater¿s functions from the Hartree- Fock method. The basis set exponents are determined from the Slater rules, Hartree-Fock atom optimized and Hartree-Fock molecule optimized. VQMC and DQMC yielded the average energy of each system. Later, Slater¿s functions are changed to the STO-6G basis functions. The same basis set exponents are applied for the STO-6G calculations. Finally, the Hartree product is used as a wave function for the VQMC and DQMC calculations with the same basis functions already mentioned. The main conclusions fro this work are: a) DQMC, as expected, presents lower energies when compared to VQMC; b) DQMC calculations using single Slater determinant and basis set with molecule and atom optimized exponents and no correlation factor provided energies compared to a Gaussian double zeta basis set at the coupled cluster including singles and doubles excitations level of theory; c) STO-6G must be used with caution in order to represent STO functions; d) the energies calculated with the Hartree product presented a behavior not far from the Hartree-Fock wave functions when localized orbitals were used; e) better results are expected if orbitals are self-consistent with respect to the Hartree method. / Mestrado / Físico-Química / Mestre em Química

We are interested in complex electronic structures of various atomic and ionics systems. We use an ab initioapproach, the multiconfigurational Dirac-Hartree-Fock (MCDHF), to compute atomic structures and properties.We contribute in three main ways to the already existent literature: by developing and implementing originalcomputer programs, by investigating possibilities of alternative computational methodologies and strategies, andfinally by performing accurate atomic structure calculations to support other research fields, i.e. nuclear physics,astrophysics or experimental physics, through the theoretical estimation of relevant atomic data.We raise questions about the choice of the optimal orbital basis by considering finite basis sets, MCDHF orbitalbases and natural-orbital bases. We demonstrate the promising potential of the latter in the context of hyperfinestructures and hope that others will find interest in pursuing our analysis. Ultimately, our work put forward someweaknesses of the traditional optimization strategy based on the layer-by-layer optimization strategy.We also perform large-scale calculations to determine accurate atomic properties such as energy levels, hyperfinestructures, isotope shifts, transition parameters, radiative lifetimes and Landé g factors. We show through thevariety of atomic properties and atomic systems studied, the difficulty of describing, in the relativistic framework,the correlation between the spatial position of electrons due to their Coulomb repulsion.This thesis is organized in two main parts. The first one is dedicated to the theoretical and computationalbackgrounds that are needed to understand the theoretical models and the interpretation of our results. Thesecond part presents and summarizes our articles and manuscripts. They are separated in four groups, A, B, C,and D, around the themes of the atomic orbital bases, the applications to nuclear physics, the applications toastrophysics, and investigations of negative ions. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

The central problem in the modern electronic structure theory is the calculation of cor- relation energy, possibly by an approach that would account for both static and dynamic correlation in an efficient, balanced and accurate way. In this thesis, I present a collection of methods that combine the effective treatment of dynamic correlation by the coupled cluster theory with density matrix renormalization group, a well-established technique for calculations of strongly correlated systems. The connection between them is achieved via the tailored coupled clusters (TCC) ansatz, which conveniently does not impose any ad- ditional computational costs. After the successful initial assessment, we developed more efficient implementations of these methods by employing the local approaches based on pair natural orbitals. This way, we extended the range of possible applications to larger systems with thousands of basis functions. To assess the accuracy of TCC as well as its local counterparts, we performed a variety of benchmark calculations ranging from small, yet challenging systems such as the nitrogen molecule or tetramethyleneethane diradical, to larger molecules like oxo-Mn(Salen) or Fe(II)-porphyrin model. 1

No description available.

No description available.

Neste trabalho foram estudados computacionalmente os [2.2]ciclofanos ([2.2]paraciclofano (1), anti-[2.2]metaciclofano (2a), sin-[2.2]metaciclofano (2b) e [2.2]metaparaciclofano (3)), que são os [2n]ciclofanos mais simples, contendo dois anéis fenílicos conectados por duas pontes etilênicas. Os ciclofanos têm apresentado inúmeras aplicações importantes, podendo atuar como auxiliares em sínteses assimétricas e como catalisadores que simulam funções enzimáticas, apresentando seletividade em relação aos substratos. Eles são empregados tanto em químicab supramolecular quanto em áreas biomédicas. Estudos que empregam ressonância de spin eletrônico ou que investigam propriedades ópticas não-lineares dos [2.2]ciclofanos indicam que os mesmos apresentam interações transanulares, que ocorrem através de recobrimento direto entre orbitais pertencentes a anéis diferentes, through-space, ou através de recobrimento entre orbitais dos anéis e das pontes, through-bond. As interações transanulares possuem um papel fundamental na química dos ciclofanos, alterando o comportamento reacional destes compostos e as transições espectroscópicas. Apesar dos métodos de preparação de ciclofanos, desde os mais simples aos mais complexos, serem intensamente investigados, estudos computacionais, que busquem compreender as correlações entre tensão e aromaticidade, estrutura eletrônica e o mecanismo de ocorrência das interações transanulares, são raramente encontrados na literatura. Desse modo, o objetivo deste trabalho foi estudar as interações transanulares, bem como correlacionar as diferenças estruturais, a tensão sobre anéis e pontes, cargas atômicas, aromaticidade e os deslocamentos químicos, não apenas para os isômeros dos [2.2]ciclofanos, mas também seus derivados fluorados (perfluoração de um dos anéis dos [2.2]ciclofanos), bem como avaliar os efeitos de diversos substituintes (CN, Cl, C=O, NH2 e NO2) e da protonação na estrutura eletrônica do isômero [2.2]paraciclofano. As otimizações de geometria de 1, realizadas com diferentes métodos e conjuntos de funções de base, mostraram que os modelos MP2/6-31+G(d,p) e B3PW91/6-31+G(d,p) fornecem os melhores resultados em comparação com os dados de raios-x. Buscas conformacionais mostraram que 2a e 2b são confôrmeros com energias diferentes e que 3 possui dois confôrmeros degenerados. As energias relativas e de tensão das pontes, seguiram a mesma ordem, indicando que a tensão sobre as pontes e a repulsão entre as nuvens ? dos anéis aromáticos são determinantes para a estabilidade dos [2.2]ciclofanos. As reações isodésmicas indicaram que os anéis comportam-se como absorvedores de tensão. NICS e HOMA mostraram que apesar das perdas de planaridade dos anéis a aromaticidade é mantida. O método NBO confirmou que todos os [2.2]ciclofanos apresentam interações through-bond, mas apenas 2a e 2b apresentaram interações through-space significantes. A análise AIM mostrou que as interações transanulares observadas são do tipo camada fechada (iônica ou ligação de hidrogênio) e que estabilizam os [2.2]ciclofanos. Para os derivados fluorados as principais alterações geométricas observadas foram para os diedros das pontes. As reações isodésmicas revelaram que as tensões das pontes e as energias relativas são afetadas pela fluoração. Além disso, os anéis dos isômeros fluorados absorvem mais tensão que os anéis dos isômeros não fluorados. NICS e HOMA mostraram que a substituição por flúor aumenta a aromaticidade dos [2.2]ciclofanos. A análise NBO indicou que a perfluoração aumentou o número e a intensidade das interações through-space, mas as mesmas ficaram restritas principalmente aos derivados fluorados de 2a e 2b. A mesma análise evidenciou que há uma conjugação dos pares de elétrons dos átomos de flúor com o sistema ?. Por outro lado, a análise AIM sugeriu que a substituição não aumenta o número de interações through-space, mas confirmou a conjugação dos pares de elétrons dos átomos de flúor. Os demais substituintes empregados afetam os parâmetros geométricos do [2.2]paraciclofano (1) de maneira diferenciada. A análise particionada das reações isodésmicas mostrou que as tensões nos anéis e nas pontes dependem não apenas do substituinte empregado, mas também da posição da substituição. NICS e HOMA indicaram que a aromaticidade no anel não-substituído dos derivados substituídos é maior que em 1. A análise NBO revelou que a substituição e a protonação aumentam a ocorrência de interações transanulares through-space. O método AIM indicou a presença de interações transanulares apenas para o derivado substituído com NH2 e CN e para a espécie protonada. No entanto, tais interações apresentaram características de interações de camada fechada. com pequenas estabilizações. As cargas atômicas e os deslocamentos químicos confirmaram as mudanças na densidade eletrônica, observadas através do método AIM. / In this work, the [2.2]cyclophanes ([2.2]paracyclophane (1), anti-[2.2]metacyclophane (2a), syn-[2.2]metacyclophane (2b) e [2.2]metaparacyclophane (3)), which are the simplest [2n] cyclophanes that contain two phenyl rings connected by two ethanediyl linkages, were studied computationally. Cyclophanes have presented several important applications, such as auxiliary in asymmetric synthesis, catalysts that simulate enzymatic functions, presenting selectivity in relation to the substrates. They are employed either in supramolecular chemistry or in biomedical areas. Studies that apply electron spin resonance or that investigate the non-linear optical properties of [2.2]cyclophanes, indicate that these compounds present transannular interactions, which occur through direct overlap of orbitals lying in different rings, throughspace, or through overlap between orbitals from rings and bridges, through-bond. The transannular interactions have a fundamental role in cyclophane chemistry, changing the reactional behavior of these compounds, and the spectroscopic transitions. Despite the fact that the well known methods of preparation, from the simplest to the most complex cyclophanes, have been studied intensively, computational studies that intent to comprehend the correlations between tension and aromaticity, electronic structure, and the mechanism of the transannular interactions are rarely found in the literature. Therefore, the aim of this work was not only to study the transannular interactions, correlating the structural differences, tension in rings and bridges, atomic charges, aromaticity, and chemical shifts of the [2.2]cyclophanes isomers but also to extent a similar treatment to the fluorinated derivatives. In addition, the effects of substituents such as (CN, Cl, C=O, NH2, and NO2) and the protonation on the electronic structure of [2.2]paracyclophane were also evaluated. The geometry optimizations of 1, carried out by using different methods and basis set, showed that the models MP2/6-31+G(d,p) and B3PW91/6-31+G(d,p) provide the best results in comparison with the x-ray data. Conformational searches showed that 2a and 2b are the conformers that present the same energy and the isomer 3 has two degenerated conformers. The strain energies of the bridges followed the same tendency as the relative energies, indicating that the tension on the bridges and the repulsions between the ? clouds of the aromatic rings are the key factors that determine the [2.2]cyclophane stabilities. The isodesmic reactions indicated that the rings are absorbents of tension. NICS and HOMA showed that the aromaticity of the rings is preserved despite the changes on the planarity. The NBO method confirmed that all [2.2]cyclophanes present through-bond interactions, but only 2a and 2b exhibit noteworthy through-space interactions. The AIM analysis pointed out that the transannular interactions behave as closed shell interactions (ionic or hydrogen bond), stabilizing the [2.2]cyclophanes. The main geometric changes, observed to the fluorinated derivatives, were those related with the dihedral angle of bridges. The isodesmic reactions pointed out that the tensions of bridges and the relative energies are affected by the fluorination. In addition, the fluorinated rings absorb more tension than the non-fluorinated rings. NICS and HOMA showed that the substitution by fluorine increases the aromaticity of the [2.2]cyclophanes. The NBO analysis indicated that the number of through-space interactions increase with the fluorination, but it is restrict to the derivatives of 2a and 2b. In addition, the same analysis pointed out a conjugation of the fluorine lone pairs with the ? system. On the other hand, the AIM analysis suggested that the substitution do not increase the number of through-space interactions, but confirmed the conjugation of the fluorine lone pairs. The other substituents can affect the geometric parameter of 1 noticeably. The partitioned analysis of isodesmic reactions showed that the tensions in bridges and rings not only depend on the substituents employed but also on the position of substitution. NICS and HOMA pointed out that the aromaticity is bigger in the non-substituted rings of [2.2]paracyclophane derivatives than in 1. The NBO analysis showed that the substitution and protonation increase the number of through-space interactions. AIM method indicated the transannular interactions occur only to the derivate substituted by NH2 and CN, and to the protonated specie. However, these interactions presented features of closed shell interactions with small stabilizations. The atomic charges and the chemical shifts confirmed the changes of the electronic density, observed through the AIM method.

A natureza das ligações de hidrogênio e a influência destas interações na estrutura eletrônica de complexos neutros, catiônicos, aniônicos e radicalares foi estudada utilizando-se análises geométricas, energéticas, eletrônicas e topológicas. Inicialmente, verificaram-se alterações na aromaticidade do cátion pirílio após a complexação com uma a três moléculas de água. Tais complexos foram ainda estudados em meio reacional com constante dielétrica igual a da água com o emprego do modelo PCM (Polarizable Continuum Model). Adicionalmente, os efeitos da hidroxilação na estrutura eletrônica dos cátions benzopirílio e flavílio foram considerados. Posteriormente, analisaram-se os efeitos das fortes ligações de hidrogênio na ressonância do grupo carboxila em complexos formados entre o radical hidroperoxil e os ácidos fórmico, acético e trifluoroacético. Como extensão desse trabalho, estudos envolvendo complexos obtidos com e sem restrições na otimização de geometria possibilitaram obter informações a respeito da ressonância dos grupos carboxila e carboxilato quando o fluoreto de hidrogênio interage linear ou perpendicularmente com todos os átomos do ácido fórmico e do ânion formiato. O desenvolvimento das atividades supracitadas compreendeu a análise da função de onda pelos métodos NBO (Natural Bond Orbital), NSA (Natural Steric Analysis), NRT (Natural Resonance Theory) e AIM (Atoms in Molecules). As alterações em parâmetros geométricos e nas cargas atômicas foram consideradas. Uma análise energética foi realizada com o emprego do método de decomposição de energia proposto por Xantheas. As freqüências vibracionais e a intensidade das bandas do estiramento do grupo X-H, doador da ligação de hidrogênio, foram analisadas. As densidades de spin para os complexos radicalares também foram obtidas. A influência das ligações de hidrogênio e o efeito do substituinte na aromaticidade dos cátions foram verificados com o emprego dos métodos e índices NICS (Nucleus Independent Chemical Shifts), HOMA (Harmonic Oscillator Model of Aromaticity), HOSE (Harmonic Oscillator Stabilization Energy) e PDI (para-Delocalization Index). Os cálculos foram efetuados com os modelos B3LYP/6-31+G(d,p), B3LYP/6-311++G(3df,3pd) e UB3LYP/6-311++G(3df,3pd). Ocasionalmente, outras funções de base (EPR-III e cc-pVDZ), assim como o método MP2, foram utilizados para testar a confiabilidade dos resultados obtidos. As interações intermoleculares pouco alteraram a estrutura eletrônica e a aromaticidade do cátion pirílio. Analogamente, a substituição de um átomo de hidrogênio por um grupo hidroxila em diversas posições dos cátions benzopirílio e flavílio também não provocou modificações muito significativas na estrutura eletrônica desses cátions, embora tenha se verificado uma dependência da aromaticidade com a posição da hidroxila. Por outro lado, a distorção geométrica associada às ligações de hidrogênio foram responsáveis pelo incremento ou diminuição da ressonância do grupo carboxila nos ácidos fórmico, acético, trifluoroacético e do grupo carboxilato no ânion formiato. Os efeitos dos grupos doador e sacador de elétrons na estabilização dos complexos radicalares foram evidenciados. Adicionalmente, pode-se atribuir um caráter covalente parcial em algumas ligações de hidrogênio. / The nature of hydrogen bonds and their influence on electronic structure of neutral, cationic, anionic, and radical complexes was studied by using geometric, energetic, electronic, and topological analysis. The changes in aromaticity of the pyrylium cation upon complexation with one up to three water molecules were investigated. The PCM (Polarizable Continuum Model) model was employed to study the pyrylium-water complexes in a water reaction medium. In addition, the effects of hydroxylation on electronic structure of the benzopyrylium and flavilium cations were also considered. In addition, the effects of strong hydrogen bonds on carboxyl group resonance in the complexes formed between the hydroperoxyl radical and formic, acetic, and trifluoroacetic acids were analyzed. In extension of this work, studies including complexes, obtained with and without geometric restrictions, provided information about the resonance of the carboxyl and carboxylate groups when the hydrogen fluoride interacts, linear or perpendicularly, with all atoms of formic acid and formate anion. The analysis of the wavefunction by using NBO (Natural Bond Orbital), NSA (Natural Steric Analysis), NRT (Natural Resonance Theory), and AIM (Atoms in Molecules) methods was necessary to the development of the above mentioned activities. The changes in geometric parameters and atomic charges were also considered. An energetic analysis of complexes was done with the energy decomposition method proposed by Xantheas. The vibrational frequencies and the intensity of the X-H (hydrogen bond donor group) stretching bands were studied. The spin densities for the radical complexes were also obtained. The Nucleus Independent Chemical Shifts (NICS), Harmonic Oscillator Model of Aromaticity (HOMA), HOSE (Harmonic Oscillator Stabilization Energy), and PDI (para-Delocalization Index) aromaticity criteria were employed to verify the hydrogen bond influence and the effect of hydroxylation in the aromaticity of the cations. The calculations were carried out by using B3LYP/6-31+G(d,p), B3LYP/6-311++G(3df,3pd), and UB3LYP/6-311++G(3df,3pd) models. Occasionally, other basis set (EPR-III and cc-pVDZ), as well as the MP2 method, were applied to test the accuracy of the results. The intermolecular interactions lead to small alterations in the electronic structure and aromaticity of pyrylium cation. Similarly, the substitution at different positions of the benzopyrylium and flavilium cations by a hydroxyl group does not cause significant changes in the electronic structure of these cations. However, a dependence of the hydroxyl group position on aromaticity was observed. On the other hand, for formic, acetic, trifluoroacetic acids, as well as for the formate anion, the resonance of the carboxyl and carboxylate groups is affected not only by the geometric distortions but also by the hydrogen bonds. The effects of the electron-donating and electron-withdrawing groups in the stabilization of radical complexes were characterized. Furthermore, a partial covalent character can be attributed to some hydrogen bonds.