用高斯优化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


谢谢各位老师的指导

刘杭航

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

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

wzkchem5

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

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

刘杭航

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

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

刘杭航

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

嗯嗯，我之前都是用vasp的

刘杭航

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

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

niobium

你的cluster构建得不好

喵星大佬

用QE，VASP或者CASTEP，周期性处理这个体系

biogon

甭用高斯处理这种体系