第三步：与第一步相似，将fchk重新载入Multiwfn。然后以此进入功能18 → 11 → 载入log文件 → 输入2考察第二个激发态（S1）→2 磁偶极矩→n。
第四步：将Multiwfn输出的新的trdipcontri.txt文件复制到VMD的根目录中，并且覆盖之前的电偶极矩的trdipcontri.txt文件。回到刚才未关闭的VMD窗口，在文本输入窗口输入source loaddip.tcl，然后输入draw color red更换了箭头颜色和跃迁电偶极矩之间容易区分。再然后输入了dip all。得到了分子整体的跃迁磁偶极矩，如图2所示。
(, 下载次数 Times of downloads: 45)
图 2
我有如下几个问题不明白，想请各位老师解答：
1.这样在同一个图中同时绘制跃迁电偶极矩和跃迁磁偶极矩的方式对不对？
2.如图2所示跃迁电偶极矩和跃迁磁偶极矩方向相反，它们之间的夹角看似接近180°。但是根据TD-DFT的输出文件log（如图3所示）里显示E-M Angle只有1.30°显然图2中的不符，图2中的夹角更似乎是180°-1.30°=178.7°。请问这是什么原因导致的？
(, 下载次数 Times of downloads: 60)
图 3
(, 下载次数 Times of downloads: 41)
图 4
(, 下载次数 Times of downloads: 51)
图 5
请各位老师能帮忙解答一下，非常感谢！