1、1.关于激发态偶极用 TD-DFT 计算化合物的垂直激发态，想知道它的激发态偶极矩，输出文件中给出的这部分内容包括Ground to excited state Transition electric dipole moments (Au): 请问这和用单个数字表示的激发态偶极矩的关系是怎样的，是根据矢量来加出来的还是？菜鸟级，还请大牛多多指教作者：yuxuan1976我也在算激发态偶极矩，Jasminer 说的没错，你这里的数据是跃迁偶极矩作者：loovfnd好像 TDDFT 不能用 root=n 那么算激发态偶极距，试试 CIS 作者：loovfnd跃迁偶极矩不能等同于基态激发态的偶极矩变
2、化。CIS 先选择一个激发态，比如第二个。此时基态的偶极距知道了。然后 CIS（root=2，） 算出一个偶极距，这就是第二个激发态的偶极距。最后二者相减。2.G09 计算基态跃迁到第一激发态的偶极矩改变Note that this will be done for the first structure in the optimization and the last. The final case will correspond to the optimized excited state structure. The first could be a vertical excitatio
3、n if the initial geometry was the ground state optimized structure.使用不同的方法，不同的基组会得到不同的结果这很正常。因为他们求解的方法不同，忽略的都不同，所以求出的解释不同的。只有采用同一种方法，同一个基组才具有可比性。（当然这是针对你这个计算，MC 模拟每次计算，计算出的结果都不一样，因为每次随机数不可能一样）如果没有外场的条件下，偶极矩应该是一个常量。这是没错的，那是对实验来说。理论计算不可能和实验数据十分吻合，都是有误差的。而且不同的软件，就是利用同样的方法，也不是完全一致3.如何求激发态偶极矩我需要 CS2 激发态偶
4、极矩，做法如下：Gview 构建输入文件，用 AM1 优化，Gview 转成输入文件，CIS 计算激发态，再用 AM1 优化，.偶极矩一直都是：0，0，0 我用 CIS 时输入命令为# CIS（NStates=10）/6-31G(d) 这里有误吗？另外，是不是因为CS2 是直线型分子，所以无论基态激发态，偶极距均为零啊？如果只是算基态结构下激发后的偶极矩，直接在 cis 时用density=current，选择你要的 root 作布居就可以了。比较一下基态和激发态的 out 文件中布居的其它部分吧如果不一样证明输出的是激发态的相关数据，那就是说确实激发态偶极矩也为零了请问 Gaussian 的
5、优化后输出文件里的的偶极矩是指什么极矩？它的大小和磁性（磁矩）有关吗？越大越好吗？就我所了解的，偶极矩跟磁性没有必然的关系。偶极矩是描述分子中正负电荷分离的程度，体现分子极性大小，磁性主要是跟自旋平行的电子数有关，一般体现了体系单电子的多少。如何计算各激发态之间的跃迁偶极呢？比如从 S1 到 S2 之间的跃迁？另外设置 Root 和 Nstates 是怎么回事呢？Root=1,计算结果会显示三个激发态，S1,S2,S3,到底是什么意思呢？nstates 是你要算几个激发态root 是你要对哪个态做布居加 density(transition=n 或(n,m)就可得到基态到第 n 激发态或第 n
6、 激发态到第 m 激发态的跃迁偶极加上 alltransiondensity,各个激发态间的跃迁偶极,振子强度啥的就都有了加上 Alltransitiondensities 关键词，只能得到不同态之间的振子强度，跃迁偶极等信息，跃迁能信息怎么得到啊，比如想知道 S2 到 S1 的跃迁能？4. gaussian 算分子偶极矩Quote:Originally posted by wxjbuilder at 2010-05-06 21:34:00:记得输出文件中默认的就有偶极距当然也可以加入关键词好像是 Polarity（记不太清楚）基态算法：# opt b3lyp/6-31g(d) scrf=(c
7、pcm,solvent=acetone,read) geom=connectivity跃迁偶极矩算法：# b3lyp/6-31g(d) scrf=(cpcm,solvent=acetone) geom=connectivity TD(sing lets,Nstates=10,alltransitiondensities,direct)现在请教激发态偶极矩的算法？还有跃迁的能量是不是就是算跃迁偶极矩时的能量？5.请问激发态个数算 CS2：cis=(nstates=32,root=1)/6-31+G*分别取 root1，2，3.7，到 root8 的时候 Dipole ment 就开始为零了。就是
8、说，如果取感兴趣的态，rootx,就只能算出 7 个激发态的偶极距（我的理解）；可是，如果用#ZINDO（alltransitiondensities,root=100）,就可以算出 32 个激发态，包括基态与激发态之间、激发态之间的偶极距等都可以得到（已毕业师兄留下的）。到底应该可以算出多少个激发态呢？请大家指点！万分感谢作者：S07111072CIS 计算激发态一般会高估轨道函数，你可以用 TDDFT 试试看作者：yuxuan1976QUOTE:Originally posted by S07111072 at 2009-5-13 11:09:CIS 计算激发态一般会高估轨道函数，你可以用
9、 TDDFT 试试看 请问如果用 TDDFT,命令该怎么写？一直没用过 TDDFT，和 CIS 命令的输入一样吗？关键字怎么选？多谢了！作者：yuxuan1976还有个问题，TDDFT 是不是不可以用来优化激发态？作者：mengfcQUOTE:Originally posted by yuxuan1976 at 2009-5-13 19:33:还有个问题，TDDFT 是不是不可以用来优化激发态？ TDDFT 好象不可以优化激发态吧？论坛上关于 TD 激发态结果分析 12009-04-19 19:36Excitation energies and oscillator strengths:Exc
10、ited State 1: Triplet-B1U 3.7768 eV 328.27 nm f=0.00008 - 11 0.529528 - 17 -0.45942This state for optimization and/or second-order correction.Copying the excited state density for this state as the 1-particle RhoCI density.问题 1。 为什么激发态 1，激发能 3.07768eV 有两个跃迁？8 - 11 8 - 17 而且还有一个系数是-0.45942，代表什么意思啊2。
11、最后一句话是什么意思？书上说这是“THE STATE OF INTEREST”我将此例更多的输入输出信息列于此，以方便说明。输入：#T RCIS(NStates=2,50-50)/6-31+G(d) Test输出：Excitation energies and oscillator strengths:Excited State 1: Triplet-B1U 3.7768 eV 328.27 nm f=0.00008 - 11 0.529528 - 17 -0.45942This state for optimization and/or second-order correction.Cop
12、ying the excited state density for this state as the 1-particle RhoCI density.Excited State 2: Triplet-B3U 7.4297 eV 166.88 nm f=0.00008 - 10 0.607418 - 15 -0.34842Excited State 3: Singlet-B3U 7.8294 eV 158.36 nm f=0.18848 - 10 0.642888 - 15 -0.272058 - 18 -0.10534Excited State 4: Singlet-B1U 7.9753
13、 eV 155.46 nm f=0.55328 - 11 0.660688 - 17 -0.20317*Population analysis using the SCF density.*Orbital symmetries:Occupied (AG) (B1U) (AG) (B1U) (B2U) (AG) (B3G) (B3U)Virtual (B1U) (AG) (B2G) (B2U) (B3U) (B3G) (AG) (B1U)(B2G) (AG) (B2U) (B1U) (B3G) (B1U) (AG) (B2U)(B3U) (B2G) (AG) (B1U) (B2U) (B3G)
14、(AG) (B1U)(B1U) (B3G) (B3U) (B1G) (AU) (AG) (B2U) (AG) (B1U)(B2G) (B3G) (B1U) (AG) (B1U)The electronic state is 1-AG.先来解释输入：我们初始输入的是单重态，为 1-AG。因为作的是 CIS，所以在所有的激发中，只允许单电子激发。“NStates=2,50-50”说明，我们要解的是两个单重态和两个三重态，这从最后的输出（4 个激发态，两三两单）也可以得到证实。再来详细看输出：当输入中没有指定 Root=N 的时候，默认的是 N=1，意思是说我们感兴趣的（也就是要做的）是第一激发态，
15、这个态被优化并计算其广义密度。这就是“This state for optimization and/or second-order correction. Copying the excited state density for this state as the 1-particle RhoCI density.”的来由。激发态 1，是轨道 8 上的一个电子跃迁到轨道 11 上，跃迁后两个电子自旋相同，因而产生的激发态 1 是三重态。轨道 8 是 B3U，轨道 11 是 B2G，耦合后为 B1U，所以激发态 1 为 Triplet-B1U。系数是该组态在该激发态的 CI 展开式中的最大系
16、数，此处也可由“IOp(9/40=N)”来控制输出更多的系数（所有大于 10(-N)的系数）。激发态 1 并非单组态，也包含了由“8 - 17”的单电子跃迁产生的组态。这里要注意，符合激发态 1 的跃迁必须是能产生 Triplet-B1U 态的跃迁。轨道 17 和 11 的对称性相同，所以“8 - 17”也在激发态 1 的组态之中。激发态 3 为单重态。例如 “8 - 10” ，是说轨道 8 上的一个电子跃迁到轨道 10 上，跃迁后两个电子自旋相反，因而产生的激发态 3 是单重态。轨道 8 是 B3U，轨道 10 是 AG，耦合后为 B3U，所以激发态 3 为 Singlet-B3U。激发态
17、2、4 可同理分析可是，正如楼上人所说，817 跃迁的系数是个负数 -0.45942，做何解释呢？一般系数的正负应该表示的是相因子。因为 CI 的展开式里是各组态函数的组合，所以这儿的正负跟我们平常想象的电子云的取向好像联系不起来。这儿的正负是个很模糊的东西，因为在各组态里已经包含了所有轨道的贡献，因而不大好说。一般关心的是其绝对值。另外，“计算结果中得到了二个单重态和二个三重态” ，这是由我们输入的 “NStates=2,50-50”来控制的。若改为“NStates=3,50-50”，按 Manual 中描述的便可得到三个单重态和三个三重态，但我照此运行 e9_01.gjf却得到两个单重态和
18、四个三重态?在书里 P213 写到cis=(root=n) specifies which excited state is to be studied (used for geometry optimizations, population analysis, and other single-state procedures). The default is the first(lowest excited state(n=1)在这个例子中，我写了 root=1,所以 OPT 后的就是第 1 激发态的优化结果，即我想要的第一激发态的激发能和键长等相关？如果 ROOT=0，是不是对象基态啊？
19、当输入中没有指定 Root=N 的时候，默认的是 N=1，意思是说我们感兴趣的（也就是要做的）是第一激发态，这个态被优化并计算其广义密度。这就是“This state for optimization and/or second-order correction. Copying the excited state density for this state as the 1-particle RhoCI density.”的来由。有时 50-50 确实有问题，所以建议还是分开计算 。免得出错。不过你提到的问题我没遇到过，有试算过小分子，root1 和 2 时结构不同！我不知道你是怎么输入文
20、件的。反正。 root=n 代表让程序只优化你想要的那个第 n 态。nstate 这个参数不加时，程序默认计算 3 个态。大于 3 就要用这个参数 。问题 3比如第一个激发态对应着 8-11 和 8-17 两个跃迁。为什么是两个呢？然后每一个激发对应的跃迁波动方程系数有正有负，为什么呢？第一个和第四个激发态对应相同的两个跃迁，而多重度不同，why？振子强度是不是代表了跃迁的可能行？一般不用管因子的正负号。真正有意义的是因子的平方和。要求 0.52952*0.52952 + 0.45942*0.45942 约等于 0.5，如果是 UCIS 计算的话，就要求约等于 1.0。之所以有多个跃迁，是因为这几个跃迁产生的电子组态都能组合出给定的电子态。应注意，这里虽然有多个跃迁，但并不是说这是个多参考类型的计算，它和多组态或多参考方法有本质的不同。同一个电子组态完全可以产生不同多重度的电子态，原因就是电子自旋的方向可以不同。f 与跃迁矩平方和激发能之乘积成正比，吸收光谱强度 I 也有这样的关系。因此，I 正比于 f，从 f 可预测吸收峰的强度。
