| 6 | 1/1 | 返回列表 |
| 查看: 1464 | 回復(fù): 5 | |||||
xiuhui木蟲(chóng) (正式寫(xiě)手)
|
[求助]
算頻率出錯(cuò)。
|
|
最近,我用G03計(jì)算一系列有機(jī)分子體系,當(dāng)烷基鏈增大時(shí),如C21,計(jì)算頻率就會(huì)出現(xiàn) CPHF failed to converge in LinEq2.錯(cuò)誤。不知道該如何解決,請(qǐng)高手指教。 計(jì)算方法如下。 ****************************************** %Nproc=4 Will use up to 4 processors via shared memory. %nproc=4 Will use up to 4 processors via shared memory. %chk=C11-xianan-C11.chk %mem=10GB ----------------------------------------------------- # b3lyp/6-311+g(d,p) opt(tight,nolinear) freq=noraman 錯(cuò)誤信息: 1 vectors were produced by pass991. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass992. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass993. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass994. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass995. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass996. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass997. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass998. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 13. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass999. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. PrismC: NFx= 2048 NFxT= 14 NFxU= 14. 1 vectors were produced by pass***. CPHF failed to converge in LinEq2. |
高斯專(zhuān)輯 |
專(zhuān)家顧問(wèn) (職業(yè)作家)
地溝油冶煉專(zhuān)家
![]() |
專(zhuān)家經(jīng)驗(yàn): +458 |
|
參見(jiàn)第九項(xiàng) http://people.smu.edu/wzou/qcsc/faq/gaussian.html |

木蟲(chóng) (著名寫(xiě)手)
歷史控+原著控

木蟲(chóng) (正式寫(xiě)手)
|
第九項(xiàng)是這樣的 9. 問(wèn): HF/DFT解析頻率計(jì)算不收斂,該如何解決? 答: 在原子數(shù)很少的情況下(通常5個(gè)原子以下),可以考慮改用數(shù)值頻率(Freq=Numerical或Freq=DoubleNumer)。或者修改控制CPHF/CPKS方程收斂的IOP選項(xiàng),有以下幾個(gè)。在單點(diǎn)能輸出之后,找到類(lèi)似于下面的輸出(如果計(jì)算加上了#P開(kāi)關(guān),位于"(Enter .../g09/l1002.exe)"的下一行): Minotr: Closed shell wavefunction. IDoAtm=111 Direct CPHF calculation. Differentiating once with respect to electric field. with respect to dipole field. Differentiating once with respect to nuclear coordinates. Using symmetry in CPHF. Requested convergence is 1.0D-10 RMS, and 1.0D-09 maximum. Secondary convergence is 1.0D-12 RMS, and 1.0D-12 maximum. 最后兩行定義了解CPHF/CPKS方程的四個(gè)收斂標(biāo)準(zhǔn)。 第一個(gè)標(biāo)準(zhǔn)由IOp(10/7=N)控制:10^(-N),默認(rèn)N=8 第二個(gè)標(biāo)準(zhǔn)由系統(tǒng)設(shè)置,等于第一個(gè)標(biāo)準(zhǔn)乘10 第三、四個(gè)標(biāo)準(zhǔn)通常相等,由IOp(10/16=N)控制:10^(-N),默認(rèn)N=12。用于溶劑計(jì)算。 如果是氣態(tài)分子計(jì)算,把IOp(10/7)改小一點(diǎn)就行了。其它可能用到的還有: IOp(10/73):可以增大CPHF方程的迭代次數(shù),默認(rèn)1000 IOp(10/18=3):CPHF方程僅迭代一次,不檢查收斂?梢杂脕(lái)檢查定性結(jié)果,通常會(huì)帶來(lái)幾個(gè)到幾十個(gè)波數(shù)的振動(dòng)頻率誤差,不能用于正式計(jì)算。 另外還需要注意,IOP不能在多步任務(wù)中傳遞,因此如果使用上面的IOP選項(xiàng),結(jié)構(gòu)優(yōu)化和頻率計(jì)算要分開(kāi)算。 我算的大分子體系,針對(duì)大分子有什么解決方案嗎? |
木蟲(chóng) (正式寫(xiě)手)
專(zhuān)家顧問(wèn) (職業(yè)作家)
地溝油冶煉專(zhuān)家
![]() |
專(zhuān)家經(jīng)驗(yàn): +458 |
|
頻率計(jì)算加上開(kāi)關(guān)iop(10/33=1),可以判斷屬于以下情況中的哪一種。 如果CPHF有收斂趨勢(shì),可以降低收斂閾值IOp(10/7),或者增加迭代次數(shù)IOp(10/73),也可以同時(shí)設(shè)置這兩個(gè)參數(shù)。 如果CPHF沒(méi)有收斂趨勢(shì)(通常迭代是發(fā)散的)。少數(shù)情況下是數(shù)值誤差造成的,用同樣的結(jié)構(gòu)重新做頻率計(jì)算,運(yùn)氣好的話也許能避免這個(gè)問(wèn)題。一般的解決辦法是換理論方法(泛函和基組)重新優(yōu)化結(jié)構(gòu)和計(jì)算頻率。如果原子數(shù)在20以下,也可以算半數(shù)值頻率(Freq=Numerical),但這通常比解析頻率耗時(shí)。 |

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +7 | 迷人的哈哈 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料專(zhuān)碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
[考研] 335分 | 材料與化工專(zhuān)碩 | GPA 4.07 | 有科研經(jīng)歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分,求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過(guò) +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|