Aromaticity changes along the reaction coordinate connecting the cyclobutadiene dimer to cubane and the benzene dimer to hexaprismane

Aromaticity and reactivity are two deeply connected concepts. Most of the thermally allowed cycloadditions take place through aromatic transition states, while transition states of thermally forbidden reactions are usually less aromatic, if at all. In this work, we perform a numerical experiment to discuss the change of aromaticity that occurs along the reaction paths that connect two antiaromatic units of cyclobutadiene to form cubane and two aromatic rings of benzene to yield hexaprismane. It is found that the aromaticity profile along the reaction coordinate of the [4+4] cycloaddition of two antiaromatic cyclobutadiene molecules goes through an aromatic highest energy point and finishes to an antiaromatic cubane species. Up to our knowledge, this represents the first example of a theoretically and thermally forbidden reaction path that goes through an intermediate aromatic region. In contrast, the aromaticity profile in the [6+6] cycloaddition of two aromatic benzene rings show a slow steady decrease of aromaticity from reactants to the highest energy point and from this to the final hexaprismane molecule a plunge of aromaticity is observed. In both systems, the main change of aromaticity occurs abruptly near the highest energy point, when the distance between the centers of the two rings is about 2.2 angstrom.