计算化学公社

标题: 求助，NWchem对Au2分子进行BSSE校正计算报错 [打印本页]

作者
Author: lzeng 时间: 2025-10-23 12:53
标题: 求助，NWchem对Au2分子进行BSSE校正计算报错
我用的是pbe0/dhf-qzvpp-2c方法基组对Au2分子进行计算，以下是我的输入文件，基组定义部分肯定是没错的，成功优化出来了，但是加上BSSE就会报错，想请教各位老师我该怎么修改
echo

start pbe0-bsse
PRINT low
geometry noautosym
Au 0.000000 0.000000 0.000000
Au 0.000000 0.000000 2.5000000
end

basis "ao basis" spherical
Au  s
  24.200000000    1.4474134140
  22.000000000    -3.9770944675
  20.000000000    2.7597109130
Au  s
  5.4043553001    1.0000000000
Au  s
  1.4532546773    1.0000000000
Au  s
0.68294298420    1.0000000000
Au  s
0.19002755303    1.0000000000
Au  s
0.88155981740E-01 1.0000000000
Au  s
0.37739555187E-01 1.0000000000
Au  p
  15.500000000    0.13980086947
  14.000000000    -0.22761055697
  6.4227368205    0.33769361626
Au  p
  1.6595601681    -0.63892663191
Au  p
0.79402913993    1.0000000000
Au  p
0.35125155397    1.0000000000
Au  p
0.81000000000E-01 1.0000000000
Au  p
0.25000000000E-01 1.0000000000
Au  d
  8.6000000000    0.15056335030
  7.8531899379    -0.20276601819
  1.7650856798    0.32865434018
Au  d
0.79482263260    1.0000000000
Au  d
0.33085396360    1.0000000000
Au  d
0.12398535801    1.0000000000
Au  f
1.76074 1.0
Au  f
0.72482 1.0
Au  f
0.29838 1.0
Au  f
0.12283 1.0
Au  g
1.69917 1.0
Au  g
0.65941 1.0
Au  p
  17.6710       -0.0748
8.0323          0.5285
Au  p
3.6510       -0.8456
Au  d
3.8000000000       1.0000000000
END

ECP
Au nelec 60
Au ul
2    1.0000000    0.0000000
Au s
2    13.5232180 426.6418670
2    6.2643840 36.8006680
Au p
2    11.4138670 87.0020910
2    10.3292150 174.0043700
2    5.7074240    8.8706100
2    4.8281650 17.9024380
Au d
2    7.4309630 49.8836550
2    8.3219900 74.6845490
2    4.6096420    6.4862270
2    3.5115070    9.5468210
Au f
2    3.0846390    8.7916400
2    3.0247430 11.6584560
Au g
2    3.9784420 -5.2343370
2    4.0114910 -6.7381420
END
SO
Au p
2 11.413867 -174.004182
2 10.329215    174.004370
2    5.707424    -17.741219
2    4.828165    17.902438
Au d
2    7.430963 -49.883655
2    8.321990    49.789700
2    4.609642    -6.486227
2    3.511507    6.364547
Au f
2    3.084639 -5.861094
2    3.024743    5.829228
Au g
2    3.978442    2.617169
2    4.011491 -2.695257
END

RELATIVISTIC
ZORA
END

charge 2
scf; direct; end

dft
iterations 50
xc pbe0
mult 1
decomp
disp vdw 4
end

bsse
mon first 1
input_wghost "scf\; maxiter 200\; end"
mon second 2
END

task sodft optimize

作者
Author: sobereva 时间: 2025-10-24 02:49
算化学键问题根本不应当用counterpoise，下文说了
计算化学键键能时以counterpoise方式考虑BSSE不仅是多余的甚至是有害的
http://sobereva.com/381（http://bbs.keinsci.com/thread-6047-1-1.html）

作者
Author: lzeng 时间: 2025-10-24 09:00

sobereva 发表于 2025-10-24 02:49
算化学键问题根本不应当用counterpoise，下文说了
计算化学键键能时以counterpoise方式考虑BSSE不仅是多余 ...

社长您好，您说的情况我了解了，但是审稿人认为弱相互作用考虑BSSE校正很重要，叫我们补充上，还是想请老师指导一下用NWchem计算时应该怎么写输入文件，可能以后的学习当中也能用到，非常感谢

作者
Author: sobereva 时间: 2025-10-25 01:35

lzeng 发表于 2025-10-24 09:00
社长您好，您说的情况我了解了，但是审稿人认为弱相互作用考虑BSSE校正很重要，叫我们补充上，还是想请老 ...

直接怼回去，别被外行审稿人牵着走

北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）这页幻灯片还专门说了这事
(, 下载次数 Times of downloads: 1)

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