PSI4 如何做激发态计算

小桥流水1102

最近在学习用PSI4 软件进行EOM-CC2（EOM-CCSD）激发态计算，计算中需要使用关键词roots_per_irrep，同过查看说明书和软件示例文件，仍不知其意义和使用。诚心请教PSI4大神，这个关键词如何使用？

以下是PSI4对 roots_per_irrep的解释和相关示例：

roots_per_irrep：

Number of excitedstates per irreducible representation for EOM-CC and CC-LR calculations. Irrepsdenote the final state symmetry, not the symmetry of the transition.

(1)  Type:array

(2)  Default:No Default

PSI4 示例1（psi4/samples/cc47/input.dat）：

#! EOM-CCSD/cc-pVDZ on H2O2 with two excited states in each irrep

set {

  basis cc-pVDZ

  roots_per_irrep [2, 2]

}

示例2（psi4/samples/cc12/input.dat）：

#! Single point energies of multiple excited states with EOM-CCSD

set {

  basis cc-pVDZ

  roots_per_irrep [2, 2, 2, 2]

}

cloudfire

你好，请问没有学过点群的不可约表示，请问怎么学会计算激发态呀，然后怎么输出波函数信息你知道吗？

sobereva

小桥流水1102 发表于 2018-10-26 10:00
感谢 sob 老师指点。我还想问下：PSI4计算中设置ROOTS_PER_IRREP [2,2,2,2] （每个不可约表示计算2个根） ...

Nstates=2与之对应，手册里写着

小桥流水1102

sobereva 发表于 2018-10-25 22:43
例：
使用EOM-CC2/6-31G*计算水分子，每个不可约表示计算2个根（当然，实际研究中不能用这么点的基组）
m ...

感谢 sob 老师指点。我还想问下：PSI4计算中设置ROOTS_PER_IRREP [2,2,2,2] （每个不可约表示计算2个根） 与Gaussian进行EOM-CCSD计算中（Nstates=10）是怎样的关联？

sobereva

例：
使用EOM-CC2/6-31G*计算水分子，每个不可约表示计算2个根（当然，实际研究中不能用这么点的基组）
molecule h2o {
   O     0.0000000000   0.0000000000  -0.0647163165
   H     0.0000000000  -0.7490459647   0.5135474533
   H     0.0000000000   0.7490459647   0.5135474533
}

set basis 6-31G*
set ROOTS_PER_IRREP [2,2,2,2]
energy ('eom-cc2')
还可以换成eom-ccsd。