获得两个偶极矩之差，指定第几激发态指令哪个是正确的？

jxystudio

大致思路：需要使用双能级公式，所以涉及到两个偶极矩之差

第一步：基态已经优化完毕，得到了基态的偶极矩

第二步：计算了前50个激发态,查看.log文件，得到S21的f是最大的

但是双能级公式的那个偶极矩，好像是需要第21激发态的偶极矩，并非跃迁偶极矩（可以按照sob老师的方法，得到），综合论坛以及手册，于是准备按照，以下指令计算第21激发态的偶极矩

# b3lyp/gen td=(root=21) density
# b3lyp/gen td=(read,root=21) density
# b3lyp/gen td=(root=21，nstates=50) density=rhoci
# b3lyp/gen td=(root=21，nstates=50) density

不知道哪个是正确的，求助大家

UW_0728.

只考察第21个激发态的性质的话，nstates可以设25-30，50偏大了，会浪费一些时间
其余问题看http://sobereva.com/227，按sob老师这篇博文的说法算特定激发态的偶极矩我感觉是应该加density=rhoci的

sobereva

b3lyp/gen td(root=21,nstates=30) density=rhoci
另外，不要盲目看哪个f是最大的，下文说的关键态是指激发能最低的一批态里面明显f不是很小的。即便很高的态的f是最大的，但由于激发能大，对beta的贡献可能也不大
谈谈计算第一超极化率的双能级公式
http://sobereva.com/361（http://bbs.keinsci.com/thread-5102-1-1.html）

下文介绍了Multiwfn可以计算各个激发态自身对beta的贡献，可以容易地识别出哪个态作为关键态
使用Multiwfn对第一超极化率做双能级和三能级模型分析
http://sobereva.com/512（http://bbs.keinsci.com/thread-14501-1-1.html）