“…1, with the parameters U and V ij representing the on-site and long-range Coulomb interactions, respectively. The one-electron hopping matrix elements t ij are restricted between nearest-neighbor carbon atoms i and j, with the value t 0 = -2.4 eV corresponding to uniform carbon-carbon bond-length r 0 = 1.4 Å, in accordance with our earlier calculations on conjugated polymers, polyaromatic hydrocarbons, and graphene quantum dots [18][19][20][21][22][23][24] . For the non-uniform bond-lengths, the values of corresponding t ij are determined from the exponential formula t ij = t 0 e (r 0 −r ij )/δ , extensively used by us earlier 25 , in which r ij is the distance between ith and jth carbon atoms (in Å), t 0 = -2.4 eV r 0 = 1.4 Å, and δ = 0.73 Å is a parameter depicting electron-phonon coupling.…”