|
Gaussian 讲解56-Gaussian 09 关键词Dreiding，EOMCCSDDreiding 参见分子力学方法的说明 EOMCCSD 说明 请求用 EOM-CCSD 方法[Koch90, Stanton93, Koch94a, Kallay04]进行激发态计算。EOM-CCSD 是对 CCSD 的扩充，用于模拟激发态。它为激发态计算提供 CCSD 级别的精度，并需要大体相当的计算量（像 CCSD 一样，标度都是 N6 ）和附加磁盘空间。这个方法在一开始用 CIS 计算产生 激发态的初始猜测，之后进行 EOM-CCSD 分析。 注意：EOM-CCSD 计算使用阿贝尔群（不能使用更高的点群）。 选项 只能用 NState 和 NStPIR 其中的一个指定所需态数。如果指定了两者，则 NState 优先。 如果都没指定，默认为 NStPIR=2。 NStates=N 在 EOM 中，尝试求解最低的 N 个态。考虑到 CIS 和 EOM 之间可能存在激发态 重新排序，好的办法是把 N 设定为大于所需的态数。 NStPIR=K 在 EOM 中，对每一对称类型求解的态数。默认是 2。注意，对称类型对应于 最大的阿贝尔子群。 如果 K 小于零，则读取一个由空行终止的输入部分，在其中指定每个对称类 型（不可约表示）的态数。要快速确定对称性顺序，可以用链接 0 命令%KJob L301 运行一个预备计算。建议同时对 CIS 指定合理的态数 NCISState（见下）。 Singlets 求解单重激发态。该选项只影响闭壳层体系的计算，并且对这种情况是默认 的。 Triplets 求解三重激发态。该选项只影响闭壳层体系的计算。如果同时求解单、三态， 必须同时指定 Singlets。 NCISState=M 用 CIS 产生猜测的总态数。默认对 NState 是 N*Irr.Reps.，对是(K+2)* Irr.Reps.。 Root=N 指定“感兴趣的态”。默认为第一激发态（N=1）。 Convergence=N 设定计算的收敛，对能量是 10-N ，对波函是 10-(N-2) 。默认为 N=7。 CCConvergence=N 用 10-N 作为 CCSD 和基态 Z 矢量迭代的收敛。CCSDConvergence 是该选项的同 义词。默认为 N=8。 LRTransitionDensities 除了 EOM 型（非松弛）跃迁密度外，还请求计算线性响应跃迁密度[Koch90, Koch94a, Kallay04]。这个算法比默认的 EOM-CCSD 方法更严格，但计算量也更大。 注意，当 CCSD 提供了精确波函时（例如双电子体系），两种算法等价。只能 用于单重的闭壳层和开壳层体系。 EnergyOnly 为了节省时间，只算右本征矢。对于激发能计算足够，但不能算跃迁密度。 读写波幅有关的选项 默认保存波幅，用于随后的计算。可以选择从以前的计算中读入波幅。在随后的计算中 可以增加态数。用于产生猜测的 CIS 也读入矢量，如果需要更多的猜测，它会自动增加态数 （假设基组不变）。 SaveAmplitudes 把收敛的波幅保存到检查点文件，用于随后的计算（例如，使用更大的基组）。 使用该选项会导致极大的检查点文件，但是可以显著加快后面的计算。 ReadAmplitudes 从检查点文件读取收敛的波幅（如果存在的话）。注意新的计算可以使用与 原计算不同的基组，方法（如果适用），等。 ReadGroundStateAmplitudes 只读入基态的（和 Z 矢量的）振幅，不读激发态振幅。当从 EOM 单态计算变 为 EOM 三态计算时，会用到该选项。ReadGSAmplitudes 是该选项的同义词。 NewCIS 当读取 EOM 波幅时，从头开始新的 CIS 计算。该选项用于读入单态，但计算 单态和三态时。如果用 ReadAmplitudes 恢复的前一个波幅计算的基组与当 前基组不同，也需要该选项。 适用性 能量。 例子 使用 EOM-CCSD。经常使用的方法是，先执行初步的、小型的 EOM-CCSD 计算，其中求解 大量的态，接下来对感兴趣的态运行更准确的计算。下面是演示这个方法的执行路径部分：
这里是 EOM-CCSD 计算的一些示例输出。下面的输出行引入结果部分
接下来是跃迁电偶极矩，分为左右两部分。在每一行的最后给出偶极强度和振荡强度， 前后两部分的结果是相同的，因为偶极强度是跃迁电偶极矩的乘积。
对于每一个态，有一个单独部分列出了激发态的 CI 展开因子，以及相应的轨道阿贝尔 对称类型，分为左右两部分，接下来是激发类型：
有关的关键词 CCSD, CIS, SAC-CI |
