关于激发态的不可约表示疑问

万卷书万里路

我在计算一篇文章里的算例时候发现，文章中提到的激发态地不可约表示和我用文中提到的同样的泛函基组、同样的程序计算出来的都不一样，是输出文件的问题，还是我做激发态单点计算的输入文件的问题？参考文献：J. Phys. Chem. C 2013, 117, 21466−21475。
烦请各位老师指教~
文献中的数据截图：

以红色框起来的3元环为例，他原文说到，对于α和β峰对应的跃迁是B3u对称性，p-band是B2u对称性（p-band主要来自HOMO-LUMO跃迁时），我开始的时候是自己用Gaussian计算的文件，关键词分别是 # b3lyp/6-31g*(优化结构)， # td=(nstates=20) PBEPBE/tzvp IOp(9/40=4) IOp(3/76=1000003750) IOp(3/7=0625006250) (激发态单点计算，PBE38)。之后用Multiwfn分析各个跃迁，发现他在Table1里表示的峰能量、振子强度大小、跃迁归属、权重比例和跃迁偶极矩的大小数值虽有出入，因为泛函基组和他原文不同，但趋势是一样的，在看到激发态的对称性时，是完全不一样的，下面是我的输出文件截图：
我这里找到他文章中对应的几个，我算出来的α  β 是B2u, 文献里提到的是B3u, 我算出来p-band是B1u, 文献里是B2u，这对称性跟他前文Itroduction里提到的、Table1里表示的都不一样。
我去确认了他Introduction里引用的参考文献，确实是这样说的，也就是说，应该是我计算错了？
之后，我用文献里的方法LB94/TZP重新计算激发态，显示的结果跟用PBE38在Gaussian里计算的对称性一样，还是和原文不一样，因为是同样的方法计算，所以表1里的数值都对应上了，只有对称性，我查了一下，也是不能理解这个对称性表示的问题，为什么对不上。

不知道是哪一步出错了，什么原因，请教各位老师~