Hexagonal flakes as fused parallelograms: A determinantal formula for Zhang‐Zhang polynomials of the O (2, m , n ) benzenoids

Energy Decomposition Scheme for Rectangular Graphene Flakes

We show-to our own surprise-that total electronic energies for a family of m × n rectangular graphene flakes can be very accurately represented by a simple function of the structural parameters m and n with errors not exceeding 1 kcal/mol. The energies of these flakes, usually referred to as multiple zigzag chains Z(m,n), are computed for m, n < 21 at their optimized geometries using the DFTB3 methodology. We have discovered that the structural parameters m and n (and their simple algebraic functions) provide a much better basis for the energy decomposition scheme than the various topological invariants usually used in this context. Most terms appearing in our energy decomposition scheme seem to have simple chemical interpretations. Our observation goes against the well-established knowledge stating that many-body energies are complicated functions of molecular parameters. Our observations might have far-reaching consequences for building accurate machine learning models.

Zhang-Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John-Sachs Theorem

Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John–Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.

Clar Covers of Overlapping Benzenoids: Case of Two Identically-Oriented Parallelograms

We present a complete set of closed-form formulas for the ZZ polynomials of five classes of composite Kekuléan benzenoids that can be obtained by overlapping two parallelograms: generalized ribbons Rb, parallelograms M, vertically overlapping parallelograms MvM, horizontally overlapping parallelograms MhM, and intersecting parallelograms MxM. All formulas have the form of multiple sums over binomial coefficients. Three of the formulas are given with a proof based on the interface theory of benzenoids, while the remaining two formulas are presented as conjectures verified via extensive numerical tests. Both of the conjectured formulas have the form of a 2×2 determinant bearing close structural resemblance to analogous formulas for the number of Kekulé structures derived from the John-Sachs theory of Kekulé structures.