Investigation of intramolecular hydrogen bonding in naphthoquinone derivatives by quantum chemical calculations

Intramolecular hydrogen-bonded naphthoquinone derivatives are investigated using both DFT (density functional theory) and MP2 (Moller-Plesset perturbation theory) method. Three different sets of naphthoquinone derivatives, each set consisting of syn and anti conformers with respect to the substitution by (-OH/-NH2) at peri and para position of the ring, are considered for investigation. In all cases, the syn conformer is found to be energetically more stable compared with the anti conformer. However, the hydrogen bond strength of the anti conformer is found to be more as obtained from NBO (natural bond orbital) analysis and QTAIM (quantum theory of atoms in molecules) theory. The O center dot center dot center dot H-O interactions are associated with covalent character while O center dot center dot center dot H-N interactions are noncovalent in nature. The NCI (noncovalent interaction) isosurface clearly indicates the formation of O center dot center dot center dot H-O and O center dot center dot center dot H-N intramolecular hydrogen bond. The energy partition analysis indicates electrostatic interaction as the dominant force in stabilizing these systems. The aromaticity calculation by HOMA (harmonic oscillator model of aromaticity) index and Bird index indicate the syn conformer with higher aromaticity compared to the anti conformer. The overall stability of the syn conformers is thus dominated by the aromatic stabilization energy instead of hydrogen bonding energy. The global reactivity descriptors further support the higher stability of the syn conformers according to maximum hardness principle (MHP) and minimum electrophilicity principle (MEP).