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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 菜鸟求助：Fe-ZSM-5模型优化中出现问题
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 700  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 324分 085600材料与化工 +7 呆鹅oor 2026-03-27 7/350 2026-03-31 21:18 by yuq
[考研] 材料科学与工程339求调剂 +3 hyz0119 2026-03-31 3/150 2026-03-31 18:31 by 1939136013狗壮
[考研] 285求调剂 +3 FZAC123 2026-03-30 3/150 2026-03-31 17:49 by 热情沙漠
[考研] 求调剂 +9 akdhjs 2026-03-31 10/500 2026-03-31 17:14 by zhyzzh
[考研] 0856 335分 | 封装or激光加工 老师看看我！！ +9 cccchenso 2026-03-29 9/450 2026-03-31 16:37 by lishahe
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 270求调剂 +3 小杰pp 2026-03-31 4/200 2026-03-31 12:59 by wxiongid
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +5 pipiver 2026-03-30 5/250 2026-03-31 07:37 by JourneyLucky
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
[考研] 085602化工求调剂（331分） +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 21:07 by peike
[考研] 105500药学求调剂，一志愿山东大学药学，348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 一志愿北京工业大学，324分求调剂 +6 零八# 2026-03-28 6/300 2026-03-29 21:20 by nanaliuyun
[考研] 085600，专业课化工原理，321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 070300求调剂306分 +4 26要上岸 2026-03-27 4/200 2026-03-28 13:06 by 唐沐儿
[考研] 292求调剂 +14 鹅鹅鹅额额额额� 2026-03-25 15/750 2026-03-28 08:45 by WYUMater
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 打过很多竞赛，085406控制工程300分，求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭