求助gaussian结构优化分子筛不收敛

sommnus

求助各位大佬，再进行一个分子筛的优化的时候调整结构难以收敛，但是我也看不出来哪里有问题。
%chk=D:\ZDJ\1.chk
# b3lyp/3-21g* em=gd3bj

Title Card Required

0 1
Al                -4.87630298    1.65528604    0.46596085
O                 -3.01621167    1.84850694    0.46596085
P                 -2.09730156    0.31689355    0.46596085
O                 -3.53288773   -0.29652791    0.46596085
Al                -4.21113672   -1.89726542    0.46596085
O                 -5.86730544   -2.20600355    0.46596085
P                 -5.67929413   -4.25900761    0.46596085
O                 -7.45700936   -4.98594154    0.46596085
O                -11.57738885   -1.54430510    0.46596085
Si               -13.02102036   -2.15400215    0.46596085
O                -13.32763274   -0.63152099    0.46596085
Si               -12.15461329    0.41305637    0.46596085
O                -11.62083580    2.10296005    0.46596085
Al                -9.99213653    2.25304684    0.46596085
O                 -8.97569495    3.66480272    0.46596085
P                 -7.34822359    2.88561721    0.46596085
O                 -5.70647716    3.19640936    0.46596085
Al                -7.79346333   -2.98239536    0.46596085
O                 -9.65597343   -2.90542174    0.46596085
Si                -9.92575067   -1.36171754    0.46596085
H                 -4.07637916    1.55449793    1.78976734
H                 -3.68440805   -2.62267192    1.73041421
H                 -7.48126704   -3.88420649    1.68734128
O                 -3.62246227    1.67552467   -0.97472437
H                 -3.62737842    2.54022483   -1.39170093
O                 -3.56207107   -2.79115342   -1.09217201
H                 -2.71507320   -3.19921633   -0.89806051
O                 -7.39053688   -3.98578844   -1.10850814
H                 -8.16296386   -4.49067634   -1.37316784
O                 -1.14716679    0.04867005   -0.93024834
H                 -0.22058462    0.18269676   -0.71792236
O                 -4.91075047   -4.78781429    1.89907080
H                 -4.74562071   -5.73193061    1.84451315
O                 -5.00791955   -4.72095581   -0.78595707
O                 -1.26184973    0.08104509    1.69364513
O                 -6.55475262    3.76139090    1.70189138
H                 -6.92233462    4.64670210    1.75392564
O                 -6.77658358    3.43960692   -0.85000538
O                -13.15483470    0.45629717   -0.84710762
H                -13.88833535    1.05029885   -0.67178274
O                -13.04495924    0.44069920    1.71667273
O                -13.47280686   -3.00876686    2.01971960
O                -13.47280686   -3.00876686   -1.08779790
H                -14.32448287   -3.43871409   -0.98105833
O                 -9.31847705   -0.55259210   -1.05897678
H                 -9.71176544   -0.97668759   -1.82517992
O                 -9.31847705   -0.55259210    1.99089847
H                 -9.62751515   -1.04078074    2.75755831
O                -10.33965501    3.04477851    1.70304547
H                -10.66010960    3.94971434    1.70304547
H                 -9.29273920    1.02795874    0.87865270
Ni               -15.00886585   -3.93447995    1.77542869
Mo               -14.59051616    1.58404406    1.39030436

sobereva

你当前用的模型根本什么都不是，完全离谱、违背基本结构化学常识，典型的反面教材。用3-21G*算这个还极其糟糕，甚至对O连极化函数都没有。这根本不是适合Gaussian跑的体系。在缺乏基本常识的情况下做计算纯粹白浪费时间和计算资源。



用免费的第一性原理程序CP2K算这种体系最为理想，速度极快，还完美考虑周期性，结合Multiwfn创建输入文件使用还超级容易（http://sobereva.com/587），还有专门的北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）对CP2K的使用讲得极其全面透彻并给了巨量例子，很容易上手。

zjxitcc

建议立即请师兄师姐喝饮料，请他们带你入门。如果组里没有做计算稍微专业一些的人士，建议让你导师立即派你到专业的课题组学习（至少1年为宜）