计算化学公社

标题: 求助：用ORCA的SOS和Polar计算的极化率差别很大 [打印本页]

作者
Author: yaoli2 时间: 2023-2-1 09:07
标题: 求助：用ORCA的SOS和Polar计算的极化率差别很大
本帖最后由 yaoli2 于 2023-2-1 09:07 编辑

我用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

作者
Author: sobereva 时间: 2023-2-1 13:51
连泛函都不统一

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

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

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

作者
Author: yaoli2 时间: 2023-2-1 16:13
好的，谢谢。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)