计算化学公社

标题: 模仿文献方法（pbe/sbjkc）优化过程失败 [打印本页]

作者
Author: zyd 时间: 2025-4-11 08:52
标题: 模仿文献方法（pbe/sbjkc）优化过程失败
本帖最后由 zyd 于 2025-4-11 08:53 编辑

老师们好。最近我在重复文献的内容：
这是文献中采用的方法
(, 下载次数 Times of downloads: 17)
然后我首先重复团簇结构的及优化，团簇结构初始文件是从文件中得到的。我的输入文件如下：
%chk=Cd9Se7H22_03.chk
#p opt freq PBE1PBE/genecp SCRF=(CPCM,Solvent=Water)

Title Card Required

0 1
H             0 4.27007600 4.19441900 -0.64813600
H             0 -4.27007600 4.19441900 -0.64813600
H             0 -5.76751100 1.60078500 -0.64813600
H             0 -1.49743500 -5.79520300 -0.64813600
H             0 1.49743500 -5.79520300 -0.64813600
H             0 5.76751100 1.60078500 -0.64813600
H             0 1.21906800 5.63451900 -0.88504800
H             0 -1.21906800 5.63451900 -0.88504800
H             0 -5.48917100 -1.76151600 -0.88504800
H             0 -4.27010300 -3.87300400 -0.88504800
H             0 4.27010300 -3.87300400 -0.88504800
H             0 5.48917100 -1.76151600 -0.88504800
H             0 0.00000000 4.93097400 1.10493700
H             0 -4.27034900 -2.46548700 1.10493700
H             0 4.27034900 -2.46548700 1.10493700
H             0 2.13526700 2.96216500 2.83891600
H             0 -2.13526700 2.96216500 2.83891600
H             0 -3.63294400 0.36811200 2.83891600
H             0 -1.49767700 -3.33027800 2.83891600
H             0 1.49767700 -3.33027800 2.83891600
H             0 3.63294400 0.36811200 2.83891600
H             0 0.00000000 0.00000000 4.58598800
Se             -1 0.00000000 4.93095900 -0.38806300
Se             -1 -4.27033500 -2.46547900 -0.38806300
Se             -1 4.27033500 -2.46547900 -0.38806300
Se             -1 -2.13567300 1.23303100 -0.38806500
Se             -1 0.00000000 -2.46606300 -0.38806500
Se             -1 2.13567300 1.23303100 -0.38806500
Se             -1 0.00000000 0.00000000 3.09898800
Cd             -1 -4.27040900 2.46552200 -1.25911700
Cd             -1 0.00000000 -4.93104400 -1.25911700
Cd             -1 4.27040900 2.46552200 -1.25911700
Cd             -1 0.00000000 2.46548500 -1.25903800
Cd             -1 -2.13572100 1.23305900 2.22693500
Cd             -1 -2.13517300 -1.23274300 -1.25903800
Cd             -1 0.00000000 -2.46611800 2.22693500
Cd             -1 2.13517300 -1.23274300 -1.25903800
Cd             -1 2.13572100 1.23305900 2.22693500

Se    0
SP 4 1.00
   3.711000             0.055744             -.006014
   1.586000             -.510520             -.121447
   0.533900             0.480755             0.452607
   0.208500             0.810292             0.669751
SP 1 1.00
   0.078210             1.000000             1.000000
D 2 1.00
   0.4754120          1.000000
   0.2077760          1.000000
****
Cd    0
SP 4 1.00
   63.380000             -.000464             -.004874
   6.7140000             2.292877             0.189630
   5.6020000          -3.902884             -.491872
   1.9710000             2.084604             1.213501
SP 1 1.00
   0.768900             1.000000             1.000000
SP 2 1.00
   2.696000             -.064866             -.053338
   0.175800             1.019661             1.007102
SP 1 1.00
   0.056470             1.000000             1.000000
D 3 1.00
   15.510000             0.012065
   2.9410000             0.370168
   1.3790000             0.680760
D 1 1.00
   0.578200             1.000000
D 1 1.00
   0.230000             1.000000
****
H 0
3-21G
****

SE    0
SE-ECP    3    28
f potential
  1
1    2.0935600          -5.5849800
s-f potential
  2
0    20.5564400          13.2808100
2    3.1549000          51.0001100
p-f potential
  2
0    2.5849500             1.7964000
2    2.5852600          30.4378700
d-f potential
  2
0    1.3829000             2.4000100
2    1.4267100             6.5994800
CD    0
CD-ECP    3    28
f potential
  1
1    9.8132500          -8.6986400
s-f potential
  3
0    1.6047600             6.5693400
2    3.3772700          -59.3639500
2    5.0887600          93.5961000
p-f potential
  3
0    1.2517100             3.6786300
2    3.4259200          -104.5069400
2    4.2141500          145.6731600
d-f potential
  2
0    22.3533700             3.3309400
2    6.7540800          62.0661600

然后计算完之后，显示的out文件出现这样的信息
(, 下载次数 Times of downloads: 19)
log文件出现这样的信息，请问这里那里出现了问题，请各位老师指教
(, 下载次数 Times of downloads: 15)
(, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 2)

作者
Author: zyd 时间: 2025-4-11 08:54
log文件比较大，上传不了，所以只截取去了最后出现的错误信息

作者
Author: zjxitcc 时间: 2025-4-11 10:54
本帖最后由 zjxitcc 于 2025-4-11 10:55 编辑

您上传的out文件没有有效信息，需要log文件。您可以先按下文自行排查原因（尤其注意结构优化步数是否超过程序上限），几乎所有Gaussian报错及解决措施都在以下帖子里，不一定要让其他人帮你找原因
Gaussian 常见报错及解决方法；新手求助报错时的注意事项
http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: sobereva 时间: 2025-4-11 13:59

zyd 发表于 2025-4-11 08:54
log文件比较大，上传不了，所以只截取去了最后出现的错误信息

如置顶的新社员必读贴和论坛首页的公告栏非常明确所示，超过500KB的文本型文件上传前必须先压缩再上传，以节约论坛空间、节约下载耗时。认真把社员必读贴好好看了再发帖

作者
Author: sobereva 时间: 2025-4-11 14:01
文献明明用的是PBE而不是PBE0

如果几何优化不收敛，仔细看此文
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: zyd 时间: 2025-4-11 15:42

zjxitcc 发表于 2025-4-11 10:54
您上传的out文件没有有效信息，需要log文件。您可以先按下文自行排查原因（尤其注意结构优化步数是否超过 ...

好嘞，谢谢

作者
Author: zyd 时间: 2025-4-11 15:43

sobereva 发表于 2025-4-11 13:59
如置顶的新社员必读贴和论坛首页的公告栏非常明确所示，超过500KB的文本型文件上传前必须先压缩再上传， ...

嗯嗯，知道啦，谢谢老师

作者
Author: zyd 时间: 2025-4-11 15:44

sobereva 发表于 2025-4-11 14:01
文献明明用的是PBE而不是PBE0

如果几何优化不收敛，仔细看此文

谢谢老师指点，我仔细看一下

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