How to determine whether intramolecular H...H interactions can be classified as dihydrogen bonds

Different types of intramolecular H...H interactions have been analyzed using the MP2/6-311++G(d,p) level of approximation. These are C-H...H-B, C-H...H-Al, C-H...H-C, C-H...H-O, O-H...H-Al and O-H...H-B contacts. Quantum theory of atoms in molecules and natural bond orbitals methods were applied to better understand the nature of these interactions. It was found that some of the species analyzed possess the characteristics of typical hydrogen bonds, such as the O-H...O ones. The electron charge transfer from the Lewis base to the antibonding X-H (for example O-H) orbital of the Lewis acid is one such characteristic. The NBO method may be considered decisive to classify any system as dihydrogen bonded.     

pi-H...O hydrogen bonds: multicenter covalent pi-H interaction acts as the proton-donating system

The MP2 method and the Pople-style basis sets 6-311++G(d,p), 6-311++G(2df,2pd), and 6-311++G(3df,3pd) were used to perform calculations on H3O+...C2H2 and C2H3+...C2H2 complexes and related species. Hydrogen bonds existing for the analyzed complexes were investigated as well as related pi-H...O --> pi...H-O and pi-H...pi --> pi...H-pi proton-transfer processes. For some of the complexes analyzed the multicenter pi-H interaction possessing the properties of a covalent bond acts as a proton donor; more generally it is classified as the Lewis acid. The quantum theory of "atoms in molecules" (QTAIM) was also applied to deepen the nature of these interactions in terms of characteristics of bond critical points. The pi-H...O, O-H...pi, and pi-H...pi interactions analyzed here may be classified as hydrogen bonds since their characteristics are the same as or at least similar to those of typical hydrogen bonds. H...pi interactions are common in crystal structures of organic and organometallic compounds. The analyses performed here show a continuum of such interactions since there are H...pi contacts possessing the characteristics of weak intermolecular interactions on the one hand and pi-H multicenter covalent bonds on the other. Ab initio and QTAIM results support the latter statements.     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

PH~3...H~2O体系分子间相互作用的从头计算MO研究

在MP2/6-311++G(3d,3p)水平,对PH~3...H~2O体系可能存在的氢键复合物进行了全自由度能量梯度优化,发现PH~3...H~2O体系存在三个能量极小结构A,B和C。其中结构A和B以H~2O作为质子授体,结构C以PH~3作为质子授体,结构A较结构B和C分别稳定6.52kJ/mol和8.18kJ/mol。结构A具有C~s对称性,其结构中P原子和O原子间距离为354.78nm,键角OHP为171.35ⅲ,属于接近于直线的传统型氢键结构。进一步在高级电子相关校正的MP4SDTQ下,用6-311++G(3df,3pd)基加上键函数{3s3p2d1f},通过BSSE校正,精确计算了结构A的结合能为-10.84kJ/mol。     

Heteronuclear intermolecular resonance-assisted hydrogen bonds. The structure of pyrrole-2-carboxamide (PyCa)

The crystal and molecular structure of pyrrole-2-carboxamide (PyCa) determined by single crystal X-ray diffraction is presented. Molecular conformations of PyCa are also analyzed by FT-IR and NMR techniques. Additionally DFT calculations at the B3LYP/6-311++G(d,p) level of approximation are performed for dimers of PyCa and for related species. The existence of two tautomeric forms for the analyzed dimers differing in H-bond motifs, N-H...O or O-H...N, is studied. The geometrical and energetic features of such H-bonds show that these interactions may be classified as intermolecular resonance-assisted hydrogen bonds. Additionally the Bader theory is applied to determine and to analyze bond critical points.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.     

Very strong C-H...O, N-H...O, and O-H...O hydrogen bonds involving a cyclic phosphate.

Very short C-H...O, N-H...O, and O-H...O hydrogen bonds have been generated utilizing the cyclic phosphate [CH2(6-t-Bu-4-Me-C6H2O)2]P(O)OH (1). X-ray structures of (i) 1 (unsolvated, two polymorphs), 1...EtOH, and 1...MeOH, (ii) [imidazolium](+)[CH2(6-t-Bu-4-Me-C6H2O)2PO2](-)...MeOH [2], (iii) [HNC5H4-N=N-C5H4NH](2+)[(CH2(6-t-Bu-4-Me-C6H2O)2PO2)2](2-)...4CH3CN...H2O [3], (v) [K, 18-crown-6](+)[(CH2(6-t-Bu-4-Me-C6H2O)2P(O)OH)(CH2(6-t-Bu-4-Me-C6H2O)2PO2)](-)...2THF [4], (vi) 1...cytosine...MeOH [5], (vii) 1...adenine...1/2MeOH [6], and (viii) 1...S-(-)-proline [7] have been determined. The phosphate 1 in both its forms is a hydrogen-bonded dimer with a short O-H...O distance of 2.481(2) [triclinic form] or 2.507(3) A [monoclinic form]. Compound 2 has a helical structure with a very short C-H...O hydrogen bond involving an imidazolyl C-H and methanol in addition to N-H...O hydrogen bonds. A helical motif is also seen in 5. In 3, an extremely short N-H...O hydrogen bond [N...O 2.558(4) A] is observed. Compounds 6 and 7 also exhibit short N-H...O hydrogen bonds. In 1...EtOH, a 12-membered hydrogen-bonded ring motif, with one of the shortest known O-H...O hydrogen bonds [O...O 2.368(4) A], is present. 1...MeOH is a similar dimer with a very short O(-H)...O bond [2.429(3) A]. In 4, the deprotonated phosphate (anion) and the parent acid are held together by a hydrogen bond on one side and a coordinate/covalent bond to potassium on the other; the O-H...O bond is symmetrical and very strong [O...O 2.397(3) A].     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.     

Magnitudes and chemical consequences of R(3)N(+)-C-H...O[double bond]C hydrogen bonding

The magnitude of the stabilizing interaction between an aliphatic C[bond]H bond attached to an ammonium nitrogen and a carbonyl oxygen was evaluated by ab initio calculations at the MP2/6-311++G** level of theory. Attractive R(3)N(+)-C-H...O[double bond]C interactions play an important role in supramolecular recognition and various types of stereoselective catalysis. Our calculations show that R(3)N(+)-C-H...O[double bond]C is the strongest hydrogen bond of the C-H...O type known to date. Such hydrogen bonds remain as stabilizing interactions even in water for amide acceptors.     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.     

Influence of N-H...O and C-H...O hydrogen bonds on the 17O NMR tensors in crystalline uracil: computational study

We report a computational study for the 17O NMR tensors (electric field gradient and chemical shielding tensors) in crystalline uracil. We found that N-H...O and C-H...O hydrogen bonds around the uracil molecule in the crystal lattice have quite different influences on the 17O NMR tensors for the two C=O groups. The computed 17O NMR tensors on O4, which is involved in two strong N-H...O hydrogen bonds, show remarkable sensitivity toward the choice of cluster model, whereas the 17O NMR tensors on O2, which is involved in two weak C-H...O hydrogen bonds, show much smaller improvement when the cluster model includes the C-H...O hydrogen bonds. Our results demonstrate that it is important to have accurate hydrogen atom positions in the molecular models used for 17O NMR tensor calculations. In the absence of low-temperature neutron diffraction data, an effective way to generate reliable hydrogen atom positions in the molecular cluster model is to employ partial geometry optimization for hydrogen atom positions using a cluster model that includes all neighboring hydrogen-bonded molecules. Using an optimized seven-molecule model (a total of 84 atoms), we were able to reproduce the experimental 17O NMR tensors to a reasonably good degree of accuracy. However, we also found that the accuracy for the calculated 17O NMR tensors at O2 is not as good as that found for the corresponding tensors at O4. In particular, at the B3LYP/6-311++G(d,p) level of theory, the individual 17O chemical shielding tensor components differ by less than 10 and 30 ppm from the experimental values for O4 and O2, respectively. For the 17O quadrupole coupling constant, the calculated values differ by 0.30 and 0.87 MHz from the experimental values for O4 and O2, respectively.     

The magnitude of [C-H...O] hydrogen bonding in molecular and supramolecular assemblies

Ab initio calculations at the MP2/6-311++G** level on model systems (N-methylpyridinium complexes of dimethyl ether and dimethyl phosphate anion) provide quantitative measures of the large stabilization energies that arise from [C-H...O] contacts in charged systems. These attractive interactions control (i) the self-assembly of bipyridinium-based catenanes and rotaxanes in solution, (ii) the self-organization of left-handed Z-DNA with alternating [dC-dG] sequences in the solid state, and (iii) the binding of pyridinium derivatives with single- and double-stranded DNA. Slightly attractive interactions occur between the donor ether and phosphate moieties and a neutral pyridine molecule in the gas phase. Electrostatic potential and solvation calculations demonstrate that [C-H...O] interactions which involve a cationic [C-H] donor are dominated by electrostatic terms.     

O-H...O interactions involving doubly charged anions: charge compression in carbonate-bicarbonate crystals

The O-H...O interaction formed by the anions HCO(3)(-) and CO(3)(2-) has been investigated on the basis of data retrieved from the Inorganic Crystal Structure Database (ICSD) and by means of ab initio computations. It has been shown that the O-H...O separations associated with HCO(3)(-)...(3)(2-) interactions are shorter than those found in crystals containing hydrogen carbonate monoanions such as HCO(3)(-)...HCO(3)(-). Ab initio MP2/6-311G++(2d,2p) computations on the crystal Na(3)(HCO(3))(CO(3)).2H(2)O have shown that the interaction between the monoanion donor and the dianion acceptor, for example HCO(3)(-)...CO(3)(2-), is more repulsive than that between singly charged ions, for example HCO(3)(-)...HCO(3)(-), but is largely overcompensated for by anion-cation electrostatic attractions. The shortening of the (-)O-H...O(2-) interaction relative to the (-)O-H...O(-) interaction has been explained as a consequence of the increased charge compression, that is of the stronger cation-anion interactions established by the CO(3)(2-) dianions with respect to those established by monoanions, and does not reflect an increase in the strength of the (-)O-H ...O(-) interaction. To expand the structural sample in the crystal packing analysis, the structure of the novel mixed salt K(2)Na(HCO(3))(CO(3)).2H(2)O has been determined by single-crystal X-ray diffraction and compared with the structure of the salt Na(3)(HCO(3))(CO(3)).2H(2)O used in the computations.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

Noncovalent interactions between cytosine and SWCNT: curvature dependence of complexes via pi...pi stacking and cooperative CH...pi/NH...pi

Fragments of C24H12, adapted from a variety of armchair [(n,n), (n = 5, 7, and 8)] and zigzag [(m,0) (m = 8, 10, and 12)] single-walled carbon nanotube (SWCNT), are used to model corresponding SWCNTs with different diameters and electronic structures. The parallel binding mainly through pi...pi stacking interaction, as well as the perpendicular binding via cooperative NH...pi and CH...pi between cytosine and the fragments of SWCNT have been extensively investigated with a GGA type of DFT, PW91LYP/6-311++G(d,p). The eclipsed tangential (ET) conformation with respect to the six-membered ring of cytosine and the central ring of SWCNT fragments is less stable than the slipped tangential (ST) conformation for the given fragment; perpendicular conformations with NH2 and CH ends have higher negative binding energy than those with NH and CH ends. At PW91LYP/6-311++G(d,p) level, two tangential complexes are less bound than perpendicular complexes. However, as electron correlation is treated with MP2/6-311G(d,p) for PW91LYP/6-311++G(d,p) optimized complexes, it turns out there is an opposite trend that two tangential complexes become more stable than three perpendicular complexes. This result implies that electron correlation, a primary source to dispersion energy, has more significant contributions to the pi...pi stacking complexes than to the complexes via cooperative NH...pi and CH...pi interactions. In addition, it was found for the first time that binding energies for two tangential complexes become more negative with increasing nanotube diameter, while those for three perpendicular complexes have a weaker dependence on the curvature; i.e., binding energies are slightly less and less negative. The performance of a novel hybrid DFT, MPWB1K, was also discussed.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

C-H...N and C-H...S hydrogen bonds--influence of hybridization on their strength

Simple complexes connected through C-H...S and C-H...N interactions are investigated: CH4...NH3, C2H4...NH3, C2H2...NH3, CH4...SH2, C2H4...SH2, and C2H2...SH2. Ab initio and DFT calculations are performed (SCF, MP2, B3LYP) using different basis sets up to the MP2/aug-cc-pVQZ//MP2/aug-cc-pVDZ level of approximation. The Bader theory is applied since MP2/6-311++G(d,p) wave functions are used to find and to characterize bond critical points in terms of electron densities and their Laplacians. The influence of hybridization on the properties of C-H...S and C-H...N systems is also studied showing that the strength of such interactions increases in the following order: C(sp3)-H...Y, C(sp2)-H...Y, C(sp)-H...Y, where Y = S, N--it is in line with the previous findings on C-H...O hydrogen bonds. The results also show that CH4...SH2 and C2H4...SH2 complexes should be rather classified as van der Waals interactions and not as hydrogen bonds. The frequency associated with the C-H stretch of C(sp3)-H...S is blue-shifted.     

A comparative computational study of FKrCCH...Y, FCCKrH...Y, and FCCH...Y complexes (Y = BF, CO, N2, OH2, OH(CH3), O(CH3)2)

The structural and spectroscopic changes in complexes of FCCKrH...Y and FKrCCH...Y (Y = BF, CO, N(2), OH(2), OH(CH(3)), O(CH(3))(2)) were computed at the MP2∕6-31++G(d,p) level of theory and compared with the corresponding properties for FCCH...Y. The computed bond length changes and frequency shifts on complexation were rationalized by comparing with a perturbation model, which gives quantitative agreement with the standard ab initio results. A recently proposed model also gives a reasonable qualitative account of the observed trends in these complexes.     

Intramolecular H-bonds: DFT and QTAIM studies on 3-(aminomethylene)pyran-2,4-dione and its derivatives

Intramolecular N-H...O hydrogen bonds in 3-(aminomethylene)pyran-2,4-dione and its simple derivatives (F, Li, and BeH substituents) were analyzed theoretically. The systems were optimized at the B3LYP/6-311++G(d,p) level of approximation. For some fluorine derivatives the corresponding tautomers with O-H...N intramolecular H-bonds were investigated, and for such pairs of tautomers, the calculations on transition states of the N-H...O <--> N...H-O proton-transfer reaction were carried out. The geometrical and energetic parameters for these species were characterized. The topological parameters derived from Bader theory were also analyzed; these are characteristics of H-bond critical points and also of ring critical points. Besides N-H...O and O-H...N intramolecular hydrogen bonds, there are the other intramolecular interactions, mostly ionic such as Be(+delta)...(-delta)O, Li(+delta)...(-delta)O, and Li(+delta)...(-delta)F. The F...O interactions also exist for some of species investigated. They may be classified as energetically stabilizing ones since the corresponding bond paths and critical points exist. The numerous correlations and dependencies between geometrical, topological, and energetic parameters were detected and described.