| 查看: 655 | 回复: 2 | ||
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
- 附件 2 : Fe-ZSM-5.out
2012-07-17 16:32:26, 3.64 K
2012-07-17 16:32:48, 261.99 K
» 猜你喜欢
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有155人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
七嗪类物质合成求助
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】MS的界面模型建立问题，旋转
已经有14人回复
- 菜鸟求助液相色谱简单问题！！！求关注
已经有12人回复
- 菜鸟求助液相色谱简单问题！！！求关注！！！求围观
已经有13人回复
- 菜鸟新手求助神经网络问题
已经有6人回复
- 【求助】菜鸟求助四氢呋喃除水问题
已经有24人回复
- 【求助/交流】一只菜鸟关于细胞培养的若干问题，求助老鸟~~
已经有6人回复
- 【求助】菜鸟问有关VASP作业提交的问题！
已经有21人回复
- 【求助】菜鸟请教大虾关于DSC的几个问题
已经有12人回复
- 【求助】菜鸟求教一些关于塑料加工实验的问题
已经有5人回复
- 【求助】菜鸟的问题
已经有44人回复
- 【求助】菜鸟求助回收率计算的问题
已经有7人回复
- 【求助】2011年国考菜鸟疑问求助
已经有9人回复
meteoric30
木虫 (正式写手)
- 应助: 43 (小学生)
- 金币: 1719.5
- 散金: 12
- 红花: 2
- 帖子: 436
- 在线: 320.7小时
- 虫号: 1046846
- 注册: 2010-06-24
- 性别: GG
- 专业: 催化化学
2楼2012-07-17 19:22:11
yp_03
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3378.7
- 散金: 237
- 红花: 4
- 帖子: 395
- 在线: 223小时
- 虫号: 656490
- 注册: 2008-11-17
- 性别: GG
- 专业: 其他无机非金属材料
|
你好，错误就是优化到中途出错。只优化了两步就优化不下去了。由于我刚刚接触高斯，所以具体也不知道哪里错了，请指点一下。以下是出错的地方，具体不知道什么意思。 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