Graph Theory Cracks the Hidden Geometry of Giant Carbon Molecules

Some of the most important molecules in modern materials science are also among the most beautiful: vast, sixfold-symmetric sheets of fused benzene rings that resemble miniature fragments of graphene. Now, a pair of mathematicians has produced an exact, closed-form topological fingerprint for three entire families of these giant carbon structures, work that could give chemists a sharper toolkit for predicting how polycromatic hydrocarbons behave before a single experiment is run.

Mohammed Alsharafi and Yusuf Zeren of Yildiz Technical University, publishing in the journal Discover Chemistry, focused on three benzenoid systems: the hexabenzocoronene series, the hexa-cata-hexabenzocoronene series, and the dodeca-benzo-circumcoronene series. Each family is built by adding successive sixfold-symmetric layers of hexagonal cells to a central core, producing molecules whose carbon skeletons grow in a perfectly predictable way. The familiar molecule C42H18, hexabenzocoronene itself, is the second member of the first family, while the circumcoronene-type molecule C84H24, containing thirty-one fused benzene rings, is the second member of the third.

The study rests on a simple but powerful idea from chemical graph theory: represent each carbon atom as a vertex and each carbon-carbon bond as an edge, then compress the resulting network into numbers called topological indices. These numerical invariants have been a staple of quantitative structure-property and structure-activity research since the Zagreb indices were introduced by Gutman and Trinajstić in 1972, originally to probe the relationship between total pi-electron energy and molecular structure. The forgotten index followed in 2015, the hyper-Zagreb index in 2013, and the Y-index more recently, each capturing a different aspect of how atoms branch and connect.

What makes the new work distinctive is not the indices themselves, which are all established invariants, but the unified treatment. Alsharafi and Zeren partitioned every edge in each molecular graph according to the degrees of its two endpoint atoms. In a benzenoid hydrocarbon, degree-two vertices sit on the molecular boundary and carry a hydrogen atom, while degree-three vertices are fused carbons buried in the interior network. Edges joining two boundary atoms, one boundary and one interior atom, or two interior atoms therefore encode different aspects of the boundary-to-interior balance, and every index can be evaluated exactly by summing or multiplying over these three edge classes.

From that partition, the authors derived explicit quadratic formulas for the first and second Zagreb indices, the forgotten index, and the Y-index across all three families, alongside multiplicative versions, modified reciprocal versions, degree-weighted graph entropies, generating functions, and coindices that measure pairs of non-adjacent vertices. For the hexabenzocoronene family, for example, the first Zagreb index works out to 162s squared minus 222s plus 84, where s counts the layers, and the Y-index to 1458s squared minus 2238s plus 876. Similar expressions cover the other two families, giving chemists exact structural descriptors for any member of any series without ever drawing the molecule.

The comparative analysis delivers perhaps the most striking result. For every principal additive descriptor, the hexabenzocoronene family attains the minimum value among the three systems for all sizes, while the maximum flips depending on scale: the dodeca-benzo-circumcoronene family leads for small molecules, up to layer count six, but the hexa-cata-hexabenzocoronene family overtakes it from layer count seven onward. The proof is disarmingly clean, reducing to a handful of quadratic differences whose signs change at a single root between six and seven. In the asymptotic limit of very large molecules, the hexa-cata descriptors converge to exactly four-thirds of their hexabenzocoronene counterparts, while the circumcoronene-type values converge to parity.

That crossover has a physical reading. The hexa-cata construction adds more boundary atoms and more boundary-sensitive edges per layer than the circumcoronene annulus, so its descriptors, which weight fused interior connectivity most heavily, eventually dominate as the molecules grow. The multiplicative indices tell the same story through their logarithms, with the hexabenzocoronene family proven minimal and the ratios of the other families’ logarithms approaching the same four-thirds and one limits. These are exact mathematical statements, verified numerically for the first three members of each series.

The entropy calculations add an information-theoretic layer. Defined from the weighted edge distributions underlying each index, these quantities summarize how heterogeneous the bond environment is across the molecule. The authors are careful to stress that these are graph-information measures, not thermodynamic entropies, and that larger index values primarily reflect compact fused regions rather than any direct experimental property. That candor matters, because the temptation to over-interpret topological numbers is a perennial hazard in mathematical chemistry.

Why should anyone outside graph theory care? Hexabenzocoronene derivatives are serious business in materials research. They have been selectively functionalized, incorporated into conjugated copolymers for organic field-effect transistors and polymer solar cells, studied as nanographene model compounds for aromatic character, and investigated in stacked complexes with fullerene. Graphene-like polycyclic aromatic systems more broadly are prized for their electronic properties, and the boundary-to-interior distinction that drives these indices is structurally meaningful in exactly those systems. Exact formulas for how descriptors scale with molecular size could eventually help screen candidate structures computationally.

The authors are equally explicit about the limits of the present contribution. No property-prediction model is fitted here; the formulas are offered as candidate graph-theoretical inputs for future quantitative structure-property and structure-activity investigations, where their predictive value must be established against measured or quantum-chemical data through statistical validation. Even so, the work sets a high bar for rigor: a complete, closed-form, cross-family characterization of three chemically important molecular series, with every extremal claim proved rather than asserted. For a field increasingly reliant on machine-learning descriptors of molecular structure, having exact ground truth for entire homologous families is a quietly valuable resource.

Subject of Research: Degree-based topological indices and graph entropies of large benzenoid hydrocarbon systems

Article Title: Topological characterization of large benzenoid hydrocarbon systems and their chemical significance

Article References: Alsharafi, M., & Zeren, Y. (2026). Topological characterization of large benzenoid hydrocarbon systems and their chemical significance. Discover Chemistry, 3(1), Article 574. https://doi.org/10.1007/s44371-026-01023-7

Image Credits: AI Generated

DOI: 10.1007/s44371-026-01023-7

Keywords: chemical graph theory, topological indices, benzenoid hydrocarbons, hexabenzocoronene, polycyclic aromatic hydrocarbons, Zagreb indices, graph entropy, QSPR, QSAR, nanographene, molecular descriptors, circumcoronene