Response to comment on "Superposition of waves or densities: Which is the nature of chemical resonance?"

Quantitative Resonance Theory Based on the Clar Sextet Model

Despite its great explanatory power in understanding the chemistry of polycyclic aromatic hydrocarbons (PAHs) and related systems, the Clar sextet rule still remains an intuitive and qualitative model with notable exceptions in some cases. Here we develop a quantitative theory of chemical resonance based on semilocalized Clar-type resonance structures (named the Clar resonators) consisting of variable numbers of Clar sextets and C═C bonds. The constructed wave functions of the Clar resonators are used to expand the actual wave function of the π-conjugated system obtained from a DFT or Hartree-Fock calculation. The resultant weights and one-electron energies of the Clar resonators can serve as a quantitative measure of the importance of these resonators. Implementing the theory in our open-source python code EzReson and applying it to over a thousand PAH molecules of different sizes and shapes, we show that the weight of the Clar resonators increases exponentially with increasing number of sextets and that their energy decreases linearly with the latter, thus confirming the general validity of the Clar rule. On the basis of such a large-scale resonance analysis, we propose three extended Clar rules, along with a unified quantitative model, that are able to evaluate the importance of all Clar resonators and the ring aromaticity for PAHs. Using the present theories, we uncover the essential role that the minor Clar resonators may play in correctly understanding the resonance stabilization and local aromaticity of rings, which was totally overlooked in the original Clar model.

Extension and Quantification of the Fries Rule and Its Connection to Aromaticity: Large-Scale Validation by Wave-Function-Based Resonance Analysis

The Fries rule is a simple, intuitive tool to predict the most dominant Kekulé structures of polycyclic aromatic hydrocarbons (PAHs), which is valuable for understanding the structure, stability, reactivity, and aromaticity of these conjugated compounds. However, it still remains an empirical hypothesis, with limited qualitative applications. Herein, we verify, generalize, and quantify the Fries rule based on the recently developed resonance analysis of the DFT wave functions of over 1500 PAH and fullerene molecules with over a billion Kekulé structures. The extended rules, counting the numbers of electrons within all rings (not just sextets), are able to rank the relative importance of all Kekulé structures for all considered systems. The statistically meaningful quantification also opens a way to evaluate ring aromaticity based on the resonance theory, which generally agrees well with conventional aromaticity descriptors. Furthermore, we propose a purely graph-based aromaticity indicator nicely applicable to PAHs and fullerenes, with no need of any quantum chemistry calculations, so that it can make valuable predictions for molecular properties that are related to local aromaticity.