我的步骤如下
第一步：将Gaussian输出文件fchk载入了Multiwfn。然后依次进入功能18 → 11 → 载入log文件 → 输入1考察第一个激发态（S1）→1 电跃迁偶极矩→n，最后保存了.pdb文件。
第二步：将Multiwfn输出的trdipcontri.txt以及教程中给的loaddip.tcl文件一并复制到VMD的根目录中。将之前输出的.pdb文件载入了VMD的main窗口。然后在文本窗口输入了source loaddip.tcl → draw color red → dip all，得到了分子整体的跃迁电偶极矩矢量图，并且保持VMD窗口和程序不关闭。
第三步：在Multiwfn中输入0 return → 18 → 11 → 输入1考察第一个激发态（S1）→2 磁跃迁偶极矩→n
第四步：将Multiwfn输出的新的trdipcontri.txt文件复制到VMD的根目录中，并且覆盖之前的电偶极矩的trdipcontri.txt文件。回到刚才未关闭的VMD窗口，在文本输入窗口输入source loaddip.tcl，然后输入了dip all。得到了分子整体的跃迁磁偶极矩，如图1所示。