Recently, TAO-DFT (i.e., thermally-assisted-occupation density functional theory) has been formulated to tackle such challenging problems.  Furthermore, the latter is shown to be related to the increasing concentration of active orbitals on the zigzag edges of the larger n-pyrenes.
THERM