明白了，谢谢各位老师指导

先按普通单重态算，然后对TS结构做个stable波函数稳定性检测，如果不稳定，用stable=opt看看得到什么波函数，然后guess=read在新优化出来的波函数对应的势能面上进一步用opt=TS找过渡态

反应物、产物是单重态。过渡态单重态、三重态都要算（如果都能找得到过渡态的话）。若单、三态下能垒比较接近，还要优化出自旋交叉点的结构。