No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação. / In the present work, we describe our results concerning our molecular dynamics investigation of the mechanical solvation dynamics within the linear response regime in model systems composed by liquid argon with a monoatomic or diatomic solute. The effect of molecular parameters (size, polarizability) and density has been elucidated for various solvation models. Time Correlation Functions for the solvation energy were calculated and separated into n-body (n = 2; 3) contributions distinguishing repulsive and attractive interactions in both liquid systems. In addition, we computed second time derivatives of these functions in order to describe translational, rotational, and roto-translational portions in the solutions containing the diatomics. We found that collective time correlation functions are well described by binary correlations at low liquid densities and, at high densities, ternary correlations become more important producing faster decaying collective time correlation functions due to partial cancellation effects. The repulsive and attractive time correlation functions exhibit a dynamic behavior that is independent on the solvation model due to linear scaling factors that only affect the absolute amplitudes of these functions. In general, the systems involving a rotational degree of freedom furnish smaller correlation times for the collective solvation dynamics, but stronger correlated two-body and three-body terms. Finally, this study shows that the solvation dynamics for the solution containing the diatomics relaxes predominatly by binary translational mechanisms when solvation models involving changes only in the polarizability parameter are considered. Binary attractive rotational mechanism become important in models with changes in the bond length.

Neste trabalho, sondas solvatocrômicas foram empregadas no estudo da solvatação em solventes puros e misturas aquosas de solventes polares práticos e apráticos. A partir da polaridade de misturas aquosas de solventes orgânicos, determinou-se a concentração de água interfacial de micelas catiônicas. Finalmente, avaliou-se o efeito da polaridade e força iônica interfaciais de micelas catiônicas sobre a velocidade da reação entre p-nitrofenildifenilfosfato e o íon fluoreto. No estudo de misturas aquosas de solventes orgânicos, um novo modelo de solvatação preferencial foi elaborado, segundo o qual três espécies presentes em solução competem pela camada de solvatação da sonda: água, solvente orgânico e o \"complexo\" formado por uma molécula de água e outra de solvente orgânico (Solv-Água). Através deste modelo, analisou-se o efeito da temperatura e das propriedades das sondas e dos solventes orgânicos sobre o fenômeno da solvatação. Em soluções aquosas de micelas catiônicas, observou-se uma desidratação da interface em função do aumento do grupo hidrofílico do tensotativo e da mudança de geometria micelar de esférica para cilíndrica. Verificou-se que a velocidade da reação entre pnitrofenildifenilfosfato e fluoreto é pouco afetada pela interface micelar e pode ser reproduzida, na ausência de micelas, em soluções com mesma concentração de água e força iônica que as micelas. / Solvatochromic probes have been employed to study the polarity of pure solvents and binary mixtures of water with protic and aprotic polar solvents. From polarity data of aqueous organic mixtures, the concentration of interfacial water of cationic micelles has been determined. In aqueous solutions of cationic micelles, the dehydration of interfacial region has been observed as a result of increasing the volume of the surfactant head group and changing micellar geometry from spherical to cylindrical. Finally, the effect of interfacial polarity and ionic force on the rate of the reaction between 4-nitrophenyldiphenylphosphate and fluoride ion has been evaluated. In studying aqueous organic mixtures, a new preferential solvation model has been developed, that considers the competition between three species in solution for the probe micro-solvation shell: water, organic solvent and a 1:1 \"complex\" formed by water and organic solvent (Solv-Água). Based on this new model, the effect of temperature, as well as probe and organic solvent properties on solvation has been analyzed. It has been observed that the rate of the reaction between p-nitrophenyldiphenylphosphate and fluoride ion shows little dependence on the properties of interfacial region of cationic micelles and can be reproduced, in the absence of micelles, in solutions containing the same water concentration and ionic force as the micellar pseudo-phase.

Uma das propriedades fundamentais de soluções de micelas aquosas é sua capacidade de solubililizar uma ampla faixa de solutos orgânicos com diferentes graus de polaridade e hidrofobicidade. Recentemente Quina et. aI., estudaram as interações químicas específicas que governam a solubilidade para um ampla faixa de solutos orgânicos em micelas usando relações lineares de energia livre de solvatação (LSERs), baseado nos parâmetros empíricos do soluto desenvolvido por Abraham. O presente trabalho trata a solubilização de solutos orgânicos em micelas baseado em parâmetros obtidos unicamente da estrutura molecular do soluto através de cálculos de química quântica. Nossos resultados mostram que valores experimentais de Ks para solubilização de solutos orgânicos em detergentes aniônicos (micelas de SDS, dodecil sulfato de sódio e agregados de SDS com PEG-polietilenoglicol) podem ser reproduzidos usando parâmetros teóricos. A relação linear teórica de energia livre de solvatação resultante é quimicamente condizente com LSERs baseados em parâmetros empíricos do soluto. No caso de detergentes catiônicos e não iônicos, a metodologia apresenta limitações, por que o método de cálculo de cargas formais usando a definição de análise populacional de Mulliken não descreve de maneira apropriada o caráter ácido e básico das pontes de hidrogênio entre os solutos e a água. Alternativas para resolver o problema da limitação imposta pelos parâmetros teóricos, mantendo a aproximação linear de energia livre, são discutidas. / One of the most fundamental properties of aqueous micelar solutions is their ability to solubilize a wide variety of organic solutes with quite distinct polarities and degree of hydrophobicity. Recently Quina et al., studied the specific chemical interactions governing the solubilities of a wide range of solutes in micelles using linear solvation energy relationships (LSERs) based on empirical solute parameters developed by Abraham. This work reports an investigation of the solubilization of organic solutes in micelles based on the use of theoretical solute parameters derived from quantum chemical calculations Our results show that the experimental Ks values for solubilization of organic solutes in anionic detergents (micelles of SDS-sodium dodecyl sulfate and aggregate formed between SDS and PEG-poly(ethylene glycol)) can be reproduced using theoretical parameters. The resultant theoretical linear solvation energy relationships are chemically consistent with LSERs based on empirical solute parameters. In the case of cationic and nonionic detergents, the theoretical parameters used to describe hydrogen bonding fail because methods of calculating formal charges using the concept of Mulliken populational analysis cannot describe hydrogen bonding donor and acceptor character in an appropriate way. Possible strategies to solve the theoretical descriptor problem in the context of the linear free energy approximation are discussed.

Absorção e Emissão de radiação eletromagnética por moléculas na região do UV-Vis fornece informações sobre os estados eletrônicos excitados, sendo propriedades de grande interesse devido a sua relação com processos biológicos, bem como suas possíveis aplicações em diagnósticos e tecnologia. Essas propriedades são sensíveis ao meio em que as moléculas se encontram, tornando-se assim natural a busca por métodos teóricos que possibilitam descrever essas interações. Neste trabalho, nós usamos a metodologia Sequential-Quantum Mechanics/Molecular Mechanics para estudar o espectro de absorção e de emissão de moléculas de relevância biológica, quando em solução. Simulações clássicas Monte Carlo foram usadas pra construir uma configuração eletrostática média do líquido para posterior cálculo das propriedades, feito através dos métodos multiconfiguracionais CASSCF e CASPT2. Cuidados especiais foram tomados para incluir a polarização eletrônica que o soluto sofre devido à presença do solvente. Nossa contribuição é a adaptação do método do Gradiente de Energia Livre para permitir a obtenção de geometrias de estados excitados do soluto em solução. A técnica foi implementada em um programa e aplicada com sucesso nos sistemas estudados. As perspectivas agora se abrem para a obtenção de intersecções cônicas em meio solvente, permitindo assim o estudo de decaimentos não-radiativos em sistemas solvatados. / Absorption and emission of electromagnetic radiation by molecules in the UV-Vis region yields information about the electronic excited states, being properties of great interest due its relation with biological process, as well as its possible applications in diagnostics and technology. These properties are sensible to the environment in which the molecules are, making natural the search for theoretical methods that describe these interactions. In this work, we have used the Sequential-Quantum Mechanics/Molecular Mechanics methodology to study the spectrum of absorption and emission of molecules of biological relevance, when in solution. Classical Monte Carlo simulations were carried out to construct an average solvent electrostatic configuration and used to represent the liquid in posterior quantum mechanics calculation of the properties, performed by using the multiconfigurational methods CASSCF and CASPT2. Special cares were taken to include the electronic polarization of the solute due the solvent. Our contribution is the adaptation of the Free Energy Gradient method to allow the determination of the solute excited states geometries in solution. The method was implemented in a computer program and successfully applied in the systems studied. The perspectives are now open to the calculation of conical intersection in solvent environment, enabling the study of nonradioactive decays in solvated systems.

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.

Made available in DSpace on 2017-07-20T12:40:20Z (GMT). No. of bitstreams: 1 Mariana_N_Levandowski.pdf: 3613595 bytes, checksum: 6d13f04a1872ade8572c9eae481770c6 (MD5) Previous issue date: 2017-02-13 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This study presents a theoretical and experimental analysis of the conformational equilibrium of phenylacetic acid and its derivatives, namely: o-fluorophenylacetic, ochlorophenylacetic, o-bromophenylacetic, o-nitrophenylacetic, hydroxyphenylacetic, m-fluorophenylacetic, p-fluorophenylacetic, p-nitrophenylacetic, p-hydroxyphenylacetic, p-mercaptophenylacetic and p-aminophenylacetic. Thetechniques used were computational calculations and nuclear magnetic resonance and infrared spectroscopy.The study made it possible to verify the conformational preference of each compound, determined by their energies and geometries, which were evaluated both for the isolated molecule, at distinct theory levels, and in different solvents, with the IEF-PCM solvation model. By using natural bond orbital calculations, it was possible to identify the main orbital interactions for each conformation. With the total deletion calculations, it was verified that there was a balance between the steric and hyperconjugative effects in the stabilization of the most stable structures. The majority conformers of halogenated compounds in ortho revealed the same spatial orientation between the carbonyl and the benzene ring, similar to that of the phenylacetic acid. In the derivatives, nitro and hydroxylated in ortho, there was the presence of an intramolecular hydrogen bond, but in none of these cases the interaction referred to one of the main conformers. The infrared experiments in solvents of different polarities, allowed the identification of more than one band, enabling the attribution of conformers present in the equilibrium as well as the observation of the same tendency as the one obtained theoretically. The coupling constant 1JCH was evaluated through nuclear magnetic resonance. No significant variation of this constant was observed, showing that the major contribution to the equilibrium is related to the majoritarian conformer, which prevails even with the change in the polarity of the medium, for each compound under evaluation. / Neste trabalho foi apresentado uma análise teórica e experimental do equilíbrio conformacional do ácido fenilacético e seus derivados, sendo eles: ofluorofenilacético, o-clorofenilacético, o-bromofenilacético, o-nitrofenilacético, ohidróxifenilacético, m-fluorofenilacético, p-fluorofenilacético, p-nitrofenilacético, phidróxifenilacético,p-mercaptofenilacético e p-aminofenilacético. As técnicas utilizadas foram os cálculos computacionais e as espectroscopias de ressonância magnética nuclear e infravermelho. Com este estudo foi possível verificar a preferência conformacional de cada composto, determinadas pelas suas energias e geometrias, as quais foram avaliadas tanto para a molécula isolada, em distintos níveis de teoria, quanto em diferentes solventes, com o modelo de solvatação IEFPCM. Através dos cálculos dos orbitais naturais de ligação foi possível identificar as principais interações orbitais para cada conformação. Com os cálculos de deleção total, verificou-se que houve um balanço entre os efeitos estérico e hiperconjugativo na estabilização das estruturas mais estáveis. Os confôrmeros majoritários dos compostos halogenados em orto mostraram a mesma orientação espacial entre a carbonila e o anel benzênico, semelhante ao ácido fenilacético. Nos derivados,nitrado e hidroxilado em orto, houve a presença de uma ligação de hidrogênio intramolecular, em nenhum dos casos esta interação era referente a um dos confôrmeros principais. Utilizando a técnica de infravermelho na região do fundamental da carbonila, em solventes de diferentes polaridades, permitiram identificar a presença de mais de uma banda, sendo possível a atribuição dos confôrmeros presentes no equilíbrio, bem como observar a mesma tendência obtida teoricamente. Com a espectroscopia de ressonância magnética nuclear foi avaliada a constante de acoplamento 1JCH. Não houve uma variação significativa desta constante, mostrando que a maior contribuição no equilíbrio está relacionada com o confôrmero majoritário, que prevalece mesmo com a mudança da polaridade do meio, de cada composto avaliado.

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.