| 3 | 1/1 | 返回列表 |
| 查看: 692 | 回復(fù): 2 | ||
yp_03木蟲 (正式寫手)
|
[求助]
菜鳥求助:Fe-ZSM-5模型優(yōu)化中出現(xiàn)問題
|
|
剛剛接觸量化計(jì)算,建模過程中遇到許多問題,希望各位大俠可以指點(diǎn)下。 我構(gòu)建了一個(gè)Fe-ZSM-5分子篩模型,模型是從ZSM分子篩中截取一個(gè)10元環(huán)。優(yōu)化中放松Fe-O,Al及Al周圍的Al-O-Si原子,其他原子固定。優(yōu)化過程中出現(xiàn)問題,下面是輸入文件及輸出文件,請各位高人指點(diǎn)!在下不勝感激! 輸入文件: %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. |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
你好,錯(cuò)誤就是優(yōu)化到中途出錯(cuò)。只優(yōu)化了兩步就優(yōu)化不下去了。由于我剛剛接觸高斯,所以具體也不知道哪里錯(cuò)了,請指點(diǎn)一下。以下是出錯(cuò)的地方,具體不知道什么意思。 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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|