 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 654 | 回復: 1 | ||
[求助]
基組外推保存變量的問題
|
|
對mrci計算的能量做外推，i 是對鍵長做循環(huán)，遇到了奇怪的問題，對不同的 i 值外推計算都能完成，輸出能量在i=4時還能保存到我設定的變量中，而i=5時的外推已經(jīng)完成，但在保存結果是卻顯示 VARIABLE UNDEFINED，最后的列表也只顯示到i=4的結果，從i=5開始全是0，不知是哪里設置有問題，望高人指點！ 輸入： ........... ........... text,========= extrapolation 1 ========== text, extrapolate reference energy using EX1 and correlation energy using L3; extrapolate,basis=avtz:avqz:av5z,method_c=l3,method_r=ex1,npc=2 emc1(i)=energr(i) !eref ecorr1(i)=energy(i)-emc1(i) ecorrq1(i)=energd(i)-emc1(i) etot1(i)=emc1(i)+ecorr1(i) etotq1(i)=emc1(i)+ecorrq1(i) text,========= extrapolation 2 ========== text,Extrapolate reference energy by EX1 and correlation energy by LH3 extrapolate,basis=avtz:avqz:av5z,method_c=LH3,method_r=ex1,npc=2 emc2(i)=energr(i) ecorr2(i)=energy(i)-emc2(i) ecorrq2(i)=energd(i)-emc2(i) etot2(i)=emc2(i)+ecorr2(i) etotq2(i)=emc2(i)+ecorrq2(i) {table,emc1,ecorr1,ecorrq1,etot1,etotq1,emc2,ecorr2,ecorrq2,etot2,etotq2 head,emc1,ecorr1,ecorrq1,etot1,etotq1,emc2,ecorr2,ecorrq2,etot2,etotq2} 輸出： （ i=5） *** ========= extrapolation 1 ==========text,Compute energies, extrapolate reference energy using EX1 and correlation energy using L3; Extrapolating MRCI energy with basissets AVTZ AVQZ AV5Z Running default procedure: HF-SCF000 MULTI000 MRCI Evaluating energy using basis set AVTZ Evaluating energy using basis set AV5Z Extrapolating reference energy using E(n)=ECBS+A*exp(-C*n) L=3 E= -173.71689968 EFIT= -173.71689962 L=4 E= -173.73229278 EFIT= -173.73229276 L=5 E= -173.73572689 EFIT= -173.73572688 E(CBS)= -173.73671302 Coefficients: 1.7844 1.5002 Extrapolating correlation energy using function n^{-X} with X= 3.0 L=4 E= -0.39469968 EFIT= -0.39469968 L=5 E= -0.40534419 EFIT= -0.40534419 E(CBS)= -0.41651219 Coefficients: 1.3960 Extrapolating correlation energy using function n^{-X} with X= 3.0 L=4 E= -0.42010195 EFIT= -0.42010195 L=5 E= -0.43126779 EFIT= -0.43126779 E(CBS)= -0.44298277 Coefficients: 1.4644 BASIS EREF ECORR ECORR(+Q) ETOT ETOT(+Q) AVTZ -173.71689968 -0.36617376 -0.38982303 -174.08307344 -174.10672271 AVQZ -173.73229278 -0.39469968 -0.42010195 -174.12699246 -174.15239473 AV5Z -173.73572689 -0.40534419 -0.43126779 -174.14107108 -174.16699468 CBS -173.73671302 -0.41651219 -0.44298277 -174.15322521 -174.17969579 SETTING EMC1(5) = ENERGR(I) VARIABLE ENERGY UNDEFINED, ASSUMING 0 VARIABLE EMC1 UNDEFINED, ASSUMING 0 SETTING ECORR1(5) = 0.00000000 AU VARIABLE ENERGD UNDEFINED, ASSUMING 0 VARIABLE EMC1 UNDEFINED, ASSUMING 0 SETTING ECORRQ1(5) = 0.00000000 AU VARIABLE EMC1 UNDEFINED, ASSUMING 0 SETTING ETOT1(5) = 0.00000000 AU VARIABLE EMC1 UNDEFINED, ASSUMING 0 SETTING ETOTQ1(5) = 0.00000000 AU *** ========= extrapolation 2 ==========text,Extrapolate reference energy by EX1 and correlation energy by LH3 Extrapolating MRCI energy with basissets AVTZ AVQZ AV5Z Running default procedure: HF-SCF000 MULTI000 MRCI Evaluating energy using basis set AVTZ Evaluating energy using basis set AVQZ Extrapolating reference energy using E(n)=ECBS+A*exp(-C*n) L=3 E= -173.71689969 EFIT= -173.71752755 L=4 E= -173.73229278 EFIT= -173.72401600 L=5 E= -173.71689968 EFIT= -173.72454860 E(CBS)= -173.72459623 Coefficients: 12.7804 2.5000 Extrapolating correlation energy using function (n+ 0.5)^{-X} with X= 3.0 L=4 E= -0.39469968 EFIT= -0.39469968 L=5 E= -0.36617376 EFIT= -0.36617376 E(CBS)= -0.33162993 Coefficients: -5.7472 Extrapolating correlation energy using function (n+ 0.5)^{-X} with X= 3.0 L=4 E= -0.42010194 EFIT= -0.42010194 L=5 E= -0.38982303 EFIT= -0.38982303 E(CBS)= -0.35315638 Coefficients: -6.1004 BASIS EREF ECORR ECORR(+Q) ETOT ETOT(+Q) AVTZ -173.71689969 -0.36617375 -0.38982302 -174.08307344 -174.10672271 AVQZ -173.73229278 -0.39469968 -0.42010194 -174.12699246 -174.15239472 AV5Z -173.71689968 -0.36617376 -0.38982303 -174.08307344 -174.10672271 CBS -173.72459623 -0.33162993 -0.35315638 -174.05622616 -174.07775261 SETTING EMC2(5) = ENERGR(I) VARIABLE ENERGY UNDEFINED, ASSUMING 0 VARIABLE EMC2 UNDEFINED, ASSUMING 0 SETTING ECORR2(5) = 0.00000000 AU VARIABLE ENERGD UNDEFINED, ASSUMING 0 VARIABLE EMC2 UNDEFINED, ASSUMING 0 SETTING ECORRQ2(5) = 0.00000000 AU VARIABLE EMC2 UNDEFINED, ASSUMING 0 SETTING ETOT2(5) = 0.00000000 AU VARIABLE EMC2 UNDEFINED, ASSUMING 0 SETTING ETOTQ2(5) = 0.00000000 AU EMC1 ECORR1 ECORRQ1 ETOT1 ETOTQ1 EMC2 ECORR2 ECORRQ2 ETOT2 ETOTQ2 -173.0982253 -0.36524581 -0.38631198 -173.4634711 -173.4845373 -173.0982253 -0.36524581 -0.38631198 -173.4634711 -173.4845373 -173.3650561 -0.39562373 -0.41910869 -173.7606799 -173.7841648 -173.3650561 -0.39562373 -0.41910869 -173.7606799 -173.7841648 -173.5400886 -0.40634896 -0.43101916 -173.9464376 -173.9711078 -173.5193235 -0.36572507 -0.38811021 -173.8850485 -173.9074337 -173.6581719 -0.41721537 -0.44308127 -174.0753872 -174.1012531 -173.6455559 -0.33088267 -0.35179713 -173.9764386 -173.9973531 0.0000000 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.0000000 0.0000000 |
回復此樓» 猜你喜歡
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復
-
哈爾濱工業(yè)大學（深圳）綠色化學團隊招收2026年秋季入學博士生
已經(jīng)有11人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有111人回復
-
澳大利亞昆士蘭科技大學PhD博士獎學金
已經(jīng)有0人回復
-
【211博士招生】環(huán)境化學、地學、毒理方向，擅長污染物遷移轉(zhuǎn)化降解、計算模擬者優(yōu)先
已經(jīng)有17人回復
-
求助如何利用vasp計算磁光克爾效應
已經(jīng)有1人回復
-
070300總分298求調(diào)劑
已經(jīng)有2人回復
-
一志愿廈門大學化學學碩307求調(diào)劑
已經(jīng)有10人回復
-
理論與計算化學方向調(diào)劑——初試388分（數(shù)學144）
已經(jīng)有0人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 自相關函數(shù)為什么最終會接近于0
已經(jīng)有8人回復
- 談談能量的基組外推
已經(jīng)有24人回復
- 突然感覺DFT已經(jīng)完勝MP2了
已經(jīng)有7人回復
- 從He2等的量子化學圖形解讀范德華作用
已經(jīng)有43人回復
- matlab 計算過程中工作空間變量保存及讀取 求解決內(nèi)存不足的方法
已經(jīng)有5人回復
- 【sobereva個人文集】亂彈DFT-D（從sobereva的博客中轉(zhuǎn)載過來的）
已經(jīng)有27人回復
2樓2014-04-16 14:21:06
相關版塊跳轉(zhuǎn)
我要訂閱樓主 kent1022 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[有機交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應機理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩求調(diào)劑 +8 | hyl3153942 2026-03-29 | 8/400 |
|
|
|
[考研] 一志愿華中師范化學332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +12 | 辰沐5211314 2026-03-26 | 12/600 |
|
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 085602 化學工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[碩博家園] 招收生物學/細胞生物學調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 化學調(diào)劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|