We introduce a vector-based interpretation of the chemical bond within the quantum theory of atoms in molecules (QTAIM), the bond-path framework set B = {p, q, r}, to follow variations in the 3D morphology of all bonds for the four infrared active normal modes of benzene. The bond-path framework set B comprises three unique paths p, q, and r where r is the familiar QTAIM bond concept of bond-path (r) while the two new paths p and q are formulated from the least and most preferred directions of electron density accumulation, respectively. We find 3D distortions including bond stretching/compression, torsion, and curving. We introduce two fractional measures to quantify these variations away from linearity of the bond.benzene, infrared spectroscopy, normal mode, quantum theory of atoms in molecules 1 | INTRODUCTION Two different approaches are used for the interpretation of vibrational spectroscopies using Raman Infrared (IR) [1][2][3] and inelastic incoherent neutron scattering of molecular structures. The first is the use of group theory with mathematical models of the forms and frequencies of the molecular vibrations including using normal mode coordinate analysis and symmetry assignments. The second is the use of empirical characteristic frequencies for chemical functional groups. We shall briefly illustrate the disadvantages of the first approach with the useful prototype of water ice, because it contains a mix of weak and strong bonding types. Previously, one of the current authors correlated the symmetry assignments of the zone-center normal modes of the proton ordered ice VIII and the corresponding frequencies with Raman, infrared, and incoherent neutron scattering spectra. [3] The comparison of the calculated frequencies with the Raman experiment spectra for the ice VIII O H sigma bonds for the symmetric and antisymmetric stretching modes was reasonable with only a 2% error for both of these modes. In calculations on this prototype system the errors were much larger for the bending, rotational and translational modes: up to 15%, 27%, and 10%, respectively. Contributions to the high frequency normal modes of vibration are dominated by the stronger O H sigma bonds conversely; the lower frequency modes are dominated by the vibrations of the weaker hydrogen bonds and O O bonding interactions. A reason for the much greater accuracy of the O H sigma