计算化学公社

 找回密码 Forget password
 注册 Register
Views: 236|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 模仿文献方法(pbe/sbjkc)优化过程失败

[复制链接 Copy URL]

41

帖子

0

威望

200

eV
积分
241

Level 3 能力者

本帖最后由 zyd 于 2025-4-11 08:53 编辑

老师们好。最近我在重复文献的内容:
这是文献中采用的方法

然后我首先重复团簇结构的及优化,团簇结构初始文件是从文件中得到的。我的输入文件如下:
%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文件出现这样的信息

log文件出现这样的信息,请问这里那里出现了问题,请各位老师指教

Cd9Se7H22_03.gjf (4.6 KB, 下载次数 Times of downloads: 1)


Cd9Se7H22_03.out (1.27 KB, 下载次数 Times of downloads: 2)



41

帖子

0

威望

200

eV
积分
241

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2025-4-11 08:54:45 | 只看该作者 Only view this author
log文件比较大,上传不了,所以只截取去了最后出现的错误信息

4104

帖子

4

威望

8871

eV
积分
13055

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2025-4-11 10:54:23 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2025-4-11 10:55 编辑

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

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2025-4-11 13:59:20 | 只看该作者 Only view this author
zyd 发表于 2025-4-11 08:54
log文件比较大,上传不了,所以只截取去了最后出现的错误信息

如置顶的新社员必读贴和论坛首页的公告栏非常明确所示,超过500KB的文本型文件上传前必须先压缩再上传,以节约论坛空间、节约下载耗时。认真把社员必读贴好好看了再发帖
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

5#
发表于 Post on 2025-4-11 14:01:07 | 只看该作者 Only view this author
文献明明用的是PBE而不是PBE0

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

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

41

帖子

0

威望

200

eV
积分
241

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-4-11 15:42:49 | 只看该作者 Only view this author
zjxitcc 发表于 2025-4-11 10:54
您上传的out文件没有 有效信息,需要log文件。您可以先按下文自行排查原因(尤其注意结构优化步数是否超过 ...

好嘞,谢谢

41

帖子

0

威望

200

eV
积分
241

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-4-11 15:43:53 | 只看该作者 Only view this author
sobereva 发表于 2025-4-11 13:59
如置顶的新社员必读贴和论坛首页的公告栏非常明确所示,超过500KB的文本型文件上传前必须先压缩再上传, ...

嗯嗯,知道啦,谢谢老师

41

帖子

0

威望

200

eV
积分
241

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-4-11 15:44:23 | 只看该作者 Only view this author
sobereva 发表于 2025-4-11 14:01
文献明明用的是PBE而不是PBE0

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

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 15:31 , Processed in 0.193697 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list