应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 菜鸟求助：Fe-ZSM-5模型优化中出现问题
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 703  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yp_03

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[求助] 菜鸟求助：Fe-ZSM-5模型优化中出现问题

刚刚接触量化计算，建模过程中遇到许多问题，希望各位大侠可以指点下。
我构建了一个Fe-ZSM-5分子筛模型，模型是从ZSM分子筛中截取一个10元环。优化中放松Fe-O，Al及Al周围的Al-O-Si原子，其他原子固定。优化过程中出现问题，下面是输入文件及输出文件，请各位高人指点！在下不胜感激！

输入文件：
%nprocshared=4
%mem=4000MB
%chk=C:\Users\YangPeng\Desktop\yangpeng Fe-ZSM-5.chk
#p opt b3lyp/gen geom=connectivity optcyc=100 pseudo=read scfcon=5 scfcyc=250

Title Card Required

0 1
O                0   14.24320133   32.16273799   11.75215457
O               -1   10.02089287   30.98668454   17.06094357
O               -1   12.30512533   31.07353174   15.75988594
O               -1   13.48039018   30.36691140   13.50734713
O               -1    9.86619874   32.47295617    9.42149914
O               -1    7.39512880   32.94074672   10.17584974
O               -1    7.43129112   31.14853608   16.65091252
O               -1    5.62319103   32.19859775   15.05547383
O                0   14.88669005   32.40584760    9.21254193
O               -1    6.46295318   32.34860732   12.56506531
O                0   12.40723066   32.82428573    9.97904383
O                0   14.29604392   34.54149818   10.62793512
Si               0   14.48087236   30.66495516   12.28119831
Si              -1   11.62407447   31.01037010   17.21864788
Si              -1   13.54266934   30.27019518   15.11329960
Si               0   11.36491313   31.98740200    9.08506453
Al               0   13.95652286   32.98614460   10.39269346
Si              -1    8.42775493   33.01377767    8.94313010
Si              -1    6.33236807   32.04464229   10.98802533
Si              -1    5.49059688   33.00588234   13.67030760
Si              -1    8.59651109   30.31756638   17.39080765
Si              -1    5.95869403   30.96299898   16.02929620
H               -1    5.87107928   34.40964109   13.88383633
H               -1    4.09991603   32.94232744   13.19821912
H               -1    5.90504680   29.70879281   15.26445056
H               -1    4.97102986   30.91525796   17.11701872
H               -1    8.58568952   28.92705054   16.91410166
H               -1    8.36860568   30.34001216   18.84285976
H               -1   12.06631465   29.79298089   17.91383645
H               -1   12.02266256   32.19176673   17.99731953
H               -1   14.80514273   30.85428217   15.58859248
H               -1   13.47428157   28.85924769   15.52005541
H               -1   15.87104238   30.52042541   12.73665244
H               -1   14.22002329   29.71178057   11.19293549
H               -1   15.84262964   32.44476694    9.29169228
H               -1   14.48510494   32.02108696    8.43005318
H               -1   11.64642570   32.21536893    7.66039554
H               -1   11.49635804   30.55501055    9.38817668
H               -1    8.55298930   34.40644083    8.48958220
H               -1    7.93268155   32.18031384    7.83807789
H               -1    6.60178256   30.62098149   10.74002176
H               -1    4.97162010   32.37095489   10.53770871
H               -1   15.24828067   34.66196886   10.64615588
Fe               0   12.16019654   31.62827685   11.67234455
O                0   10.93923202   30.65908231   12.61165685

Fe
Lanl2DZ
****
O Si H 0
6-31G*
****

Fe
Lanl2DZ

输出文件：
前面省略……
Cycle  16  Pass 1  IDiag  1:
RMSU=  4.31D-06    CP:  9.59D-01  9.78D-02  8.08D-01  1.52D-01  3.52D-01
                    CP:  8.29D-01  6.01D-01  7.98D-01  1.04D+00  1.27D+00
                    CP:  1.25D+00  1.46D+00  2.41D+00  2.00D+00  1.67D+00
E= -3764.35001229196     Delta-E=       -0.000004328110 Rises=F Damp=F
DIIS: error= 1.42D-05 at cycle  16 NSaved=  16.
NSaved=16 IEnMin=16 EnMin= -3764.35001229196     IErMin=16 ErrMin= 1.42D-05
ErrMax= 1.42D-05 EMaxC= 1.00D+00 BMatC= 1.69D-07 BMatP= 2.58D-06
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com:  0.212D-05-0.480D-05 0.116D-03-0.324D-05 0.279D-03-0.188D-02
Coeff-Com: -0.199D-03-0.262D-03 0.325D-02 0.431D-02 0.158D-02-0.460D-01
Coeff-Com: -0.773D-01 0.554D-01 0.437D+00 0.623D+00
Coeff:      0.212D-05-0.480D-05 0.116D-03-0.324D-05 0.279D-03-0.188D-02
Coeff:     -0.199D-03-0.262D-03 0.325D-02 0.431D-02 0.158D-02-0.460D-01
Coeff:     -0.773D-01 0.554D-01 0.437D+00 0.623D+00
Gap=     0.041 Goal=   None    Shift=    0.000
RMSDP=7.28D-06 MaxDP=5.31D-04 DE=-4.33D-06 OVMax= 1.19D-03

SCF Done:  E(RB+HF-LYP) =  -3764.35001229     A.U. after   16 cycles
             Convg  =    0.7282D-05             -V/T =  2.0220
             S**2   =   0.0000
KE= 3.683266993604D+03 PE=-1.774405201516D+04 EE= 5.883110309083D+03
Leave Link  502 at Tue Jul 17 11:04:48 2012, MaxMem=  524288000 cpu:    1004.7
(Enter /opt/g03/l701.exe)
Compute integral first derivatives.
... and contract with generalized density number  0.
   1 Symmetry operations used in ECPInt.
ECPInt:  NShTT=   10878 NPrTT=   80013 LenC2=    8063 LenP2D=   29808.
LDataN:  DoStor=F MaxTD1= 5 Len=  102
LDataN:  DoStor=T MaxTD1= 5 Len=  102
Leave Link  701 at Tue Jul 17 11:04:49 2012, MaxMem=  524288000 cpu:       2.6
(Enter /opt/g03/l702.exe)
L702 exits ... SP integral derivatives will be done elsewhere.
Leave Link  702 at Tue Jul 17 11:04:49 2012, MaxMem=  524288000 cpu:       0.0
(Enter /opt/g03/l703.exe)
Compute integral first derivatives, UseDBF=F.
Integral derivatives from FoFDir, PRISM(SPDF).
Symmetry not used in FoFDir.
MinBra= 0 MaxBra= 2 Meth= 1.
IRaf=       0 NMat=   1 IRICut=       1 DoRegI=T DoRafI=F ISym2E= 0 JSym2E=0.
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Fe-ZSM-5zeolite.gjf
    • 2012-07-17 16:32:26, 3.64 K
  • 附件 2 : Fe-ZSM-5.out
    • 2012-07-17 16:32:48, 261.99 K

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-07-17 16:33:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

【答案】应助回帖

感谢参与，应助指数 +1
错误行贴出来啊~~实在懒得下载了~~遇到什么问题了？也没讲
赞一下 回复此楼
2楼2012-07-17 19:22:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yp_03

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
2楼: Originally posted by meteoric30 at 2012-07-17 19:22:11
错误行贴出来啊~~实在懒得下载了~~遇到什么问题了？也没讲

你好，错误就是优化到中途出错。只优化了两步就优化不下去了。由于我刚刚接触高斯，所以具体也不知道哪里错了，请指点一下。以下是出错的地方，具体不知道什么意思。
Cycle  16  Pass 1  IDiag  1:
RMSU=  4.31D-06    CP:  9.59D-01  9.78D-02  8.08D-01  1.52D-01  3.52D-01
                    CP:  8.29D-01  6.01D-01  7.98D-01  1.04D+00  1.27D+00
                    CP:  1.25D+00  1.46D+00  2.41D+00  2.00D+00  1.67D+00
E= -3764.35001229196     Delta-E=       -0.000004328110 Rises=F Damp=F
DIIS: error= 1.42D-05 at cycle  16 NSaved=  16.
NSaved=16 IEnMin=16 EnMin= -3764.35001229196     IErMin=16 ErrMin= 1.42D-05
ErrMax= 1.42D-05 EMaxC= 1.00D+00 BMatC= 1.69D-07 BMatP= 2.58D-06
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
赞一下 回复此楼
3楼2012-07-18 14:03:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yp_03 的主题更新
31/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大） +3 相信必会光芒万� 2026-03-31 3/150 2026-04-01 02:16 by 小青椒26
[考研] 调剂推荐 +7 清酒714 2026-03-26 8/400 2026-03-31 22:12 by 544594351
[考研] 一志愿：西北大学，英一数一408-284分求调剂 +7 12.27 2026-03-27 7/350 2026-03-31 21:59 by lbsjt
[考研] 086000生物与医药298调剂求助 +3 元元青青 2026-03-31 5/250 2026-03-31 19:58 by vgtyfty
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 7/350 2026-03-31 18:54 by 小张做实验
[考研] 086000调剂 +5 7901117076 2026-03-26 5/250 2026-03-31 17:45 by 544594351
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 0703化学 +20 妮妮ninicgb 2026-03-27 20/1000 2026-03-31 13:33 by 无际的草原
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135，现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 一志愿南开大学0710生物学359求调剂 +5 兔兔兔111223314 2026-03-29 7/350 2026-03-30 18:29 by 兔兔兔111223314
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +7 脚滑的守法公民 2026-03-27 7/350 2026-03-29 20:21 by dophin1985
[考研] 本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 打过很多竞赛，085406控制工程300分，求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
信息提示
关闭
请填处理意见
关闭