Halina Szatylowicz

Aromaticity from the viewpoint of molecular geometry: application to planar systems.

to Planar Systems Tadeusz M. Krygowski,*,† Halina Szatylowicz,*,‡ Olga A. Stasyuk,‡ Justyna Dominikowska, and Marcin Palusiak †Department of Chemistry, Warsaw University, Pasteura 1, 02-093 Warsaw, Poland ‡Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland Department of Theoretical and Structural Chemistry, Faculty of Chemistry, University of Łod́z,́ Pomorska 163/165, 90-236 Łod́z,́ Poland

A tautomeric equilibrium between functionalized 2-formylphenylboronic acids and corresponding 1,3-dihydro-1,3-dihydroxybenzo[c][2,1]oxaboroles

Functionalized 2-formylphenylboronic acids undergo an unprecedented tautomeric rearrangement in solution to form corresponding 1,3-dihydro-1,3-dihydroxybenzo[c][2,1]oxaboroles. X-Ray analyses of selected examples revealed diverse solid-state molecular structures from a planar open form with a hydrogen-bonded carbonyl group (X = 3-F) through a twisted conformer showing a weak carbonyl–boron interaction (X = 3,5-Br2) to a cyclic oxaborole derivative (X = 3-Br). Variable-temperature 1H NMR spectroscopy was used to determine equilibrium constants as well as enthalpies and entropies of tautomerization in a mixed solvent [D6]acetone–D2O (95 : 5). A computational approach to the process by DFT (B3LYP) and MP2 methods has also been performed.

H-Bonded Complexes of Aniline with HF/F− and Anilide with HF in Terms of Symmetry-Adapted Perturbation, Atoms in Molecules, and Natural Bond Orbitals Theories

The hydrogen-bonded isoelectronic complexes of aniline with HF/F- and an ionic form of aniline with HF were investigated by use of computational methods: Symmetry-Adapted Perturbation Theory (SAPT), Atoms in Molecules (AIM), and Natural Bond Orbitals (NBO) approaches. All computations were based on structural models previously generated at the B3LYP/6-311+(d,p) level. The differences between neutral (Ph-NH2...HF)and anionic (Ph-NH2...F- and Ph-NH-...HF) complexes were clearly outlined. The discussed charged complexes serve as Lewis acids and base, HF and F-, respectively. It was found that electrostatic and induction energy terms, obtained as a result of the SAPT method, are most dependent on the type of H-bonding (i.e.,charged or neutral). The electrostatic term is the most distinctive between the neutral and charge-assisted hydrogen bonds in the investigated two-body systems, whereas the latter is more significant in the case of weaker interactions (larger H...B distances). Application of Principal Component Analysis (PCA) to energy components obtained from the SAPT procedure indicated that all of them are relatively well intercorrelated.The above-mentioned terms together with the exchange energy terms are the most important contributions ofthe main principal component, which describes 95% of the total variance. Comparison of AIM parameters in bond critical points for modeled H-bond systems shows a good agreement with those from equilibrium complexes, both experimental and calculated ones. It was found that charged H-bonded complexes exhibit larger fluctuation of electron density and its Laplacian in bond critical points, in line with SAPT analysis. NBO results confirmed the effect of the strength of interaction on property changes both in the region of H-bonding and outside of it. The latter, more distant consequences follow the Bent-Walsh rule for all studied complexes.

Effect of the H-bonding on aromaticity of purine tautomers.

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.

How H-bonding affects aromaticity of the ring in variously substituted phenol complexes with bases. 4. Molecular geometry as a source of chemical information.

Aromaticity of the ring of variously substituted phenols in their H-bonded complexes with various bases was a subject of analysis based on 664 geometries retrieved from CSD and by use of the aromaticity index HOMA. GEO and EN, the components of the HOMA index, describing a decrease of aromaticity due to an increase of bond alternation (GEO term) and bond elongation (EN term), were also studied. There is an approximate monotonic dependence of HOMA and GEO on the H-bond strength estimated by the C-O bond length of the hydroxyl group in phenols.

Hydrogen bonding as a modulator of aromaticity and electronic structure of selected ortho-hydroxybenzaldehyde derivatives.

Properties of hydrogen bonds can induce changes in geometric or electronic structure parameters in the vicinity of the bridge. Here, we focused primarily on the influence of intramolecular H-bonding on the molecular properties in selected ortho-hydroxybenzaldehydes, with additional restricted insight into substituent effects. Static models were obtained in the framework of density functional theory at B3LYP/6-311+G(d,p) level. The electronic structure parameters evolution was analyzed on the basis of Atoms In Molecules (AIM) and Natural Bond Orbitals methods. The aromaticity changes related to the variable proton position and presence of substituents were studied using Harmonic Oscillator Model of Aromaticity (HOMA), Nucleus-Independent Chemical Shift (NICS) and AIM-based parameter of Matta and Hernández-Trujillo. Finally, Car-Parrinello molecular dynamics was applied to study variability of the hydrogen bridge dynamics. The interplay between effects of the substitution and variable position of the bridged proton was discussed. It was found that the hydrogen bond energies are ca. 9-10 kcal/mol, and the bridged proton exhibits some degree of penetration into the acceptor region. The covalent character of the studied hydrogen bond was most observable when the bridged proton reached the middle position between the donor and acceptor regions. The aromaticity indexes showed that the aromaticity of the central phenyl ring is strongly dependent on the bridged proton position. Correlations between these parameters were found and discussed. In the applied time-scale, the analysis of time evolution of geometric parameters showed that the resonance strengthening does not play a crucial role in the studied compounds.

The role of the long-range exchange corrections in the description of electron delocalization in aromatic species

In this article, we address the role of the long‐range exchange corrections in description of the cyclic delocalization of electrons in aromatic systems at the density functional theory level. A test set of diversified monocyclic and polycyclic aromatics is used in benchmark calculations involving various exchange‐correlation functionals. A special emphasis is given to the problem of local aromaticity in acenes, which has been a subject of long‐standing debate in the literature. The presented results indicate that the noncorrected exchange‐correlation functionals significantly overestimate cyclic delocalization of electrons in heteroaromatics and aromatic systems with fused rings, which in the case of acenes leads to conflicting local aromaticity predictions from different criteria.

How amino and nitro substituents direct electrophilic aromatic substitution in benzene: an explanation with Kohn–Sham molecular orbital theory and Voronoi deformation density analysis

The substituent effect of the amino and nitro groups on the electronic system of benzene has been investigated quantum chemically using quantitative Kohn-Sham molecular orbital theory and a corresponding energy decomposition analysis (EDA). The directionality of electrophilic substitution in aniline can accurately be explained with the amount of contribution of the 2pz orbitals on the unsubstituted carbon atoms to the highest occupied π orbital. For nitrobenzene, the molecular π orbitals cannot explain the regioselectivity of electrophilic substitution as there are two almost degenerate π orbitals with nearly the same 2pz contributions on the unsubstituted carbon atoms. The Voronoi deformation density analysis has been applied to aniline and nitrobenzene to obtain an insight into the charge rearrangements due to the substituent. This analysis method identified the orbitals involved in the C-N bond formation of the π system as the cause for the π charge accumulation at the ortho and para positions in the case of the NH2 group and the largest charge depletion at these same positions for the NO2 substituent. Furthermore, we showed that it is the repulsive interaction between the πHOMO of the phenyl radical and the πHOMO of the NH2 radical that is responsible for pushing up the πHOMO of aniline and therefore activating this π orbital of the phenyl ring towards electrophilic substitution.

Effect of the substituent and hydrogen bond on the geometry and electronic properties of OH and O(-) groups in para-substituted phenol and phenolate derivatives.

Interrelations between intra- and intermolecular interactions were analyzed by using computational modeling of the para-X-substituted derivatives of phenol and phenolate (where X = NO, NO(2), CHO, COMe, COOH, CONH(2), Cl, F, H, Me, OMe, and OH) and their equilibrium H-bonded complexes with HB and B(-) (where HB = HF and HCN and B(-) = F(-) and CN(-)). B3LYP/6-311++G** computation was applied. Both the substituent effect and H-bonding changed the electronic properties of the -O(-) and -OH groups and geometric parameters of phenol and phenolate derivatives and their H-bonded complexes. C-O bond lengths and aromaticity indices of the ring were found to depend linearly on σ(p)(-) of the substituents. In the first case the greatest sensitivity on the substituent effect was for 4-X-C(6)H(4)OH···CN(-) and 4-X-C(6)H(4)O(-)···HF complexes, whereas for 4-X-C(6)H(4)O(-)···HCN systems it was comparable with that for phenol derivatives and a little smaller than that for 4-X-C(6)H(4)O(-) derivatives. This means that the strength of H-bonding may considerably change the sensitivity of the C-O bond length to the substituent effect. The greatest sensitivity of the aromaticity indices, both HOMA and NICS(1)zz, to σ(p)(-) was found for phenolate and then for phenolate H-bonded complexes, followed by phenol complexes, and the lowest sensitivity was observed for phenol derivatives. The interatomic proton-acceptor distance, being a measure of the H-bond strength, was found to depend linearly on σ(p)(-) of the substituents with a positive slope for O···HB (HF or HCN) interactions and a negative slope for OH···B(-) interactions. NBO charges on the oxygen and hydrogen atoms also depend on σ(p)(-) of the substituents. In the latter case for strong H-bonded complexes (energy less than ∼-20 kcal/mol) the substituent effect works oppositely for 4-X-C(6)H(4)OH···B(-) in comparison with the 4-X-C(6)H(4)O(-)···HB systems. Moreover, following the Espinoza et al. [J. Chem. Phys. 2002, 117, 5529] and Grabowski et al. [J. Phys. Chem. B 2006, 110, 6444] classifications, the above and q(H) vs proton-acceptor distance relationships suggest a partially covalent character of the hydrogen bond for these complexes and the degree of its covalent nature depending on the substituent.

Interference of H-bonding and substituent effects in nitro- and hydroxy-substituted salicylaldehydes

AbstractTwo intramolecular interactions, i.e., (1) hydrogen bond and (2) substituent effect, were analyzed and compared. For this purpose, the geometry of 4- and 5-X-substituted salicylaldehyde derivatives (X = NO2, H or OH) was optimized by means of B3LYP/6-311 + G(d,p) and MP2/aug-cc-pVDZ methods. The results obtained allowed us to show that substituents (NO2 or OH) in the para or meta position with respect to either OH or CHO in H-bonded systems interact more strongly than in the case of di-substituted species: 4- and 3-nitrophenol or 4- and 3-hydroxybenzaldehyde by ∼31%. The substituent effect due to the intramolecular charge transfer from the para-counter substituent (NO2) to the proton-donating group (OH) is ∼35% greater than for the interaction of para-OH with the proton-accepting group (CHO). The total energy of H-bonding for salicylaldehyde, and its derivatives, is composed of two contributions: ∼80% from the energy of H-bond formation and ∼20% from the energy associated with reorganization of the electron structure of the systems in question. FigureSubstituent effect stabilization energy (SESE) estimation for the salicylaldehyde and its 4- and 5-X-substituted derivatives