Previous work has shown that commonly-used GGA and hybrid functionals overstabilize the red and orange polymorphs relative to the yellow ones. This can be understood as an SIE effect stemming from over-delocalization of the π electrons in near-planar geometries, where the system is better conjugated. As shown in Fig. , both GGAs and hybrid functionals over-stabilize the more planar conformations relative to benchmark CCSD(T) calculations. Hybrid functionals reduce the SIE and improve the conformational energy profile, but errors remain significant and are large enough to erroneously alter the polymorph stability rankings. Torsional profiles computed using (DC-)B3LYP+D3 and (DC-)PBE0+D3 are shown in Fig. and are similar to the (DC-)SCAN0+D3 results in Fig. .