计算化学公社

标题: 用高斯优化Cu2O上吸附的CO2分子第一步就不收敛出错 [打印本页]

作者
Author: 刘杭航 时间: 2020-9-24 16:22
标题: 用高斯优化Cu2O上吸附的CO2分子第一步就不收敛出错
新手小白想请教一下各位老师，老师让我优化一下Cu2O上吸附一个CO2分子，我选了一个比较小的重复单元，将边缘饱和了一下，进行一个简单的优化计算，但是第一步SCF就不收敛，之前搜过sob老师的那篇解决scf不收敛的方法帖子，选用了不少方法，还是在第一步就不收敛提示l502。最后使用了scf=xqc的方法，第一步就出错，信息是l101,照着提示改了也不对，想请教一下各位老师是我的输入文件哪里出错了吗？谢谢谢谢。出错信息
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Thu Sep 24 13:24:56 2020, MaxMem=  134217728 cpu:    4.7
(Enter /usr/local/g03/l107.exe)
Allocation:
   ICR=    1 ICP= 112 IAnSav= 223
   IRI= 260 IWt= 6476 IDP=  12692  ICScr=  12803
   IRC=  12914    IC=  19130    IX=  25457 IDF=  25568
IDDF=  25679  IXInc=  38000  IEVal=  38111  IEVec=  38222
IScr1=  50543  IScr2=  50654 IEnd=  50764
Recover reactant geometry.  Energy = -0.39337242D+04
Leave Link  107 at Thu Sep 24 13:24:57 2020, MaxMem=  134217728 cpu:    1.4
(Enter /usr/local/g03/l101.exe)

Symbolic Z-matrix:
Charge =  0 Multiplicity = 1

There are no atoms!
Error termination via Lnk1e in /usr/local/g03/l101.exe at Thu Sep 24 13:24:57 2020.
Job cpu time:  4 days 14 hours 16 minutes 56.1 seconds.
File lengths (MBytes):  RWF= 663 Int=    0 D2E=    0 Chk=    54 Scr=    1

下面是我的输入文件：

%NPROCSHARED=8
%mem=1GB
%chk=Cu2O112CO2_opt.chk
#p B3LYP/gen pseudo=read opt  opt=cartesian scfcon=5 scf=(maxcycle=80,xqc)

Cu2O112CO2_opt

-5 1
Cu             -1 0.32448400 -2.72691300 -0.72921900
Cu             -1 -2.70149400 -2.53573200 -1.02417500
Cu             -1 1.85154200 -0.54544600 0.75047200
Cu             -1 -4.49650700 -2.92750300 1.40575100
Cu             0 2.14685800 2.12832000 -0.55558600
Cu             0 3.40750900 1.53183200 2.23479600
Cu             0 0.30878800 1.72718600 1.93245200
Cu             -1 -4.20041500 -0.16308400 0.16055800
Cu             -1 1.55545000 -3.30986500 1.99566400
Cu             0 -3.90509900 2.51068200 -1.14549900
Cu             0 -2.64444800 1.91419500 1.64488200
Cu             0 -5.74316900 2.10954800 1.34253900
Cu             -1 3.35046300 -2.91809400 -0.43426200
O             -1 4.81929400 1.89565800 3.33635500
O             -1 1.77053000 -2.37513500 0.39563100
O             0 1.93255400 1.28424300 1.10531200
O             -1 -4.28142700 -1.99277300 -0.19428200
O             0 -4.11940300 1.66660500 0.51539900
O             -1 -1.23266300 2.27802000 2.74644200
O             -1 -4.70749800 -3.84444400 2.97533700
H             -1 -5.52059600 -3.56498300 3.40238800
O             -1 4.90041300 -3.45074900 -1.24840600
H             -1 5.33183700 -4.11807800 -0.70975400
O             -1 -7.32766000 2.54177500 2.14967500
H             -1 -7.39151900 2.08762200 2.99304200
O             -1 -3.69597800 3.33434400 -2.76622500
H             -1 -4.46493900 3.87887600 -2.95002800
O             -1 1.34445900 -4.22680600 3.56525100
H             -1 0.53136100 -3.94734500 3.99230100
Cu             -1 -5.69987300 -2.33797000 -1.29779900
O             -1 -7.11831900 -2.68316700 -2.40131600
H             -1 -7.89968400 -2.85993200 -1.87233200
H             -1 4.53321300 1.84297300 4.25122300
O             -1 -1.09404600 -3.07199600 -1.83266300
O             0 2.29231482 3.48259659 -0.36089177
C             0 2.51438267 4.68425245 -0.06042593
O             0 2.73644942 5.88590858 0.24003906

1-13 30 0
sdd
****
14-29 31-37 0
6-31G*
****

Cu 0
ECP10MWB 3 10
F-Komponente
1
2 1.000000 0.000000
S-F
2
2 30.220000 355.770158
2 13.190000 70.865357
P-F
2
2 33.130000 233.891976
2 13.220000 53.947299
D-F
2
2 38.420000 -31.272165
2 13.260000 -2.741104

谢谢各位老师的指导

作者
Author: biogon 时间: 2020-9-24 16:42
甭用高斯处理这种体系

作者
Author: 喵星大佬 时间: 2020-9-24 16:51
用QE，VASP或者CASTEP，周期性处理这个体系

作者
Author: niobium 时间: 2020-9-25 11:24
你的cluster构建得不好

作者
Author: 刘杭航 时间: 2020-9-25 19:19

biogon 发表于 2020-9-24 16:42
甭用高斯处理这种体系

嗯嗯，我之前都是用vasp,老师让我用高斯试一下

作者
Author: 刘杭航 时间: 2020-9-25 19:20

喵星大佬发表于 2020-9-24 16:51
用QE，VASP或者CASTEP，周期性处理这个体系

嗯嗯，我之前都是用vasp的

作者
Author: 刘杭航 时间: 2020-9-25 19:21

niobium 发表于 2020-9-25 11:24
你的cluster构建得不好

我这给是直接切出来一个面选了最小的重复单元，请问老师这样构建不合理嘛

作者
Author: wzkchem5 时间: 2020-9-25 19:39

刘杭航发表于 2020-9-25 19:21
我这给是直接切出来一个面选了最小的重复单元，请问老师这样构建不合理嘛

所有的铜都应该和两个晶格氧配位（不算CO2），也就是和CO2配位的铜应该是三配位。你相当于是让CO2吸附在一配位铜上，金属原子即使在表面也很难采取一配位。
此外，吸附CO2的铜原子应该位于团簇的一个表面的中心。团簇模型是在软件不支持周期性边界条件或支持得不好的时候，不得已的一种方法，通过包括吸附位点附近有限多个原子来近似模拟无穷体系，那么必然希望表面的原子在吸附位点的各个方向上都有分布，像你现在的构建方式相当于在模拟台阶位，不是在模拟bulk表面。
最后关于计算方法：B3LYP算这种体系估计不靠谱，建议用PBE或PBE0。另外应该加色散校正。

作者
Author: 刘杭航 时间: 2020-9-26 16:27

wzkchem5 发表于 2020-9-25 19:39
所有的铜都应该和两个晶格氧配位（不算CO2），也就是和CO2配位的铜应该是三配位。你相当于是让CO2吸附在 ...

嗯嗯，我再去试一下，谢谢您！

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