求助：用ORCA的SOS和Polar计算的极化率差别很大

yaoli2

我用SOS和Polar计算得到极化率，两者差别很大。我是想重复文章(https://www.sciencedirect.com/sc ... i/S0009261416000130)的结果，所以Polar用的基组跟文章中一样。SOS输入：%tddft NRoots 400 end
                ! cam-b3lyp 6-311+G** VeryTightSCF
SOS结果（使用Multiwfn得到的）：
Polarizability tensor:              1              2              3
  1      84.180658       0.000000      -0.000000
  2      -0.000000     204.065713       0.000000
  3      -0.000000       0.000000      26.087063
Isotropic average polarizability:     104.777811
Isotropic average polarizability volume:      15.526469 Angstrom^3
Polarizability anisotropy (definition 1):     157.200019
Eigenvalues:      26.087063      84.180658     204.065713
Polarizability anisotropy (definition 2):     148.931852

Polar：%elprop Polar 3 end
         ! M062X 6-311+G** VeryTightSCF
Polar结果（log输出文件）：
-------------------------
THE POLARIZABILITY TENSOR
-------------------------


The raw cartesian tensor (atomic units):
  227.69034     79.90205      1.07136
   79.90205    312.63257     21.13822
    1.07136     21.13822    144.33364
diagonalized tensor:
  140.23204    182.17119    362.25332


    0.16720     -0.84426     -0.50919
   -0.19596      0.47770     -0.85639
    0.96625      0.24297     -0.08557


Isotropic polarizability : 228.21885

sobereva

连泛函都不统一

SOS的结果对考虑的态数很敏感，需要做随态数的收敛性测试

基于Multiwfn的SOS功能计算，必须让量化程序输出足够多的组态，下文明确说了。ORCA里%tddft中必须写tprint 1E-8或更小的值
使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率
http://sobereva.com/232

标题必须说清楚，避免任何歧义。SOS是Multiwfn提供的功能，而不是ORCA直接提供的

yaoli2

好的，谢谢。