| 24小時(shí)熱門(mén)版塊排行榜 |
| 8 | 1/1 | 返回列表 |
| 查看: 1215 | 回復(fù): 7 | |||
chembetsey木蟲(chóng) (小有名氣)
|
[求助]
QCISD(T) 頻率計(jì)算:穩(wěn)定點(diǎn)有多個(gè)虛頻問(wèn)題
|
|
在計(jì)算CH3O這樣一個(gè)簡(jiǎn)單的分子的時(shí)候,QCISD(T)的優(yōu)化成功了,但是頻率計(jì)算卻得出數(shù)個(gè)很大的虛頻。有沒(méi)有人也碰到同樣的事情過(guò)? 一個(gè)可能性是,優(yōu)化雖然成功,但是得到的并不是一個(gè)穩(wěn)定的結(jié)構(gòu)。我用QCISD算基于QCISD(T)優(yōu)化結(jié)構(gòu)的頻率卻沒(méi)有虛頻。(當(dāng)然這樣做出來(lái)的結(jié)果是meanless的)這是否說(shuō)明QCISD(T)下優(yōu)化的結(jié)構(gòu)并沒(méi)有錯(cuò)? 另,QCISD(T)的結(jié)構(gòu)優(yōu)化和頻率計(jì)算都只能是numerical 的,初始結(jié)構(gòu)輸入似乎只能是用z-matrix。根據(jù)輸入的z-matrix的不同,似乎優(yōu)化出來(lái)的結(jié)果也略有差異。這是否正常? |
木蟲(chóng) (著名寫(xiě)手)
|
這個(gè)分子的計(jì)算并不是很容易!它的基態(tài)在c3v下是簡(jiǎn)并的,因?yàn)閖ahn teller效應(yīng)而分裂,產(chǎn)生三個(gè)對(duì)稱(chēng)性破缺的極小點(diǎn)和三個(gè)鞍點(diǎn)。 有兩種可能性,一個(gè)是這里所有c3v結(jié)構(gòu)的點(diǎn)都是圓錐交叉點(diǎn),導(dǎo)數(shù)不連續(xù),而且會(huì)有兩個(gè)負(fù)無(wú)窮的虛頻,如果的jahn teller distortion比較小,計(jì)算數(shù)值導(dǎo)數(shù)時(shí)可能會(huì)在圓錐交叉附近取了點(diǎn),那么得到的結(jié)果就會(huì)不準(zhǔn)確,可以觀察一下目前的結(jié)構(gòu),考慮降低數(shù)值導(dǎo)數(shù)的步長(zhǎng),步長(zhǎng)一定要遠(yuǎn)小于到c3v結(jié)構(gòu)的距離 還有可能就是你可能取到了那三個(gè)鞍點(diǎn)中的一個(gè) 另外對(duì)于這個(gè)問(wèn)題,由于接近簡(jiǎn)并,最好還是用多組態(tài)方法。在有,為什么選擇用qcisd而不是ccsd呢?現(xiàn)在很少見(jiàn)到qcisd計(jì)算了 |
木蟲(chóng) (小有名氣)
|
忘記說(shuō)了,我也算了CCSD(T),沒(méi)有類(lèi)似的問(wèn)題。得到的構(gòu)型沒(méi)有虛頻。原本是想都做出來(lái)作對(duì)比用。 因?yàn)閖ahn teller效應(yīng),C3v的構(gòu)型不如Cs的穩(wěn)定。所以我的計(jì)算是在Cs的對(duì)稱(chēng)性下做的。我也曾試著打破對(duì)稱(chēng)性,還是有同樣的問(wèn)題。 有一點(diǎn)我覺(jué)得奇怪的是,頻率計(jì)算完之后給出的Maximum Force, RMS Force等一系列值都無(wú)比大,盡管我構(gòu)型優(yōu)化的時(shí)候還用了opt=tight。 |
木蟲(chóng) (著名寫(xiě)手)
|
波函數(shù)在交叉點(diǎn)附近不穩(wěn)定。建議限制波函數(shù)對(duì)稱(chēng)性,分別指定電子態(tài)的對(duì)稱(chēng)性為A'和A''進(jìn)行兩次優(yōu)化。限制對(duì)稱(chēng)性后的波函數(shù)是穩(wěn)定的。 最好用有cs可約表示的內(nèi)坐標(biāo),然后凍結(jié)所有A''內(nèi)坐標(biāo)(計(jì)算數(shù)值導(dǎo)數(shù)時(shí)也不考慮這些內(nèi)坐標(biāo)),這樣應(yīng)該可約比較順暢地得到這兩個(gè)限制優(yōu)化的極小點(diǎn)。a'和a''中一個(gè)能量更高的應(yīng)該是鞍點(diǎn),低的那個(gè)是極小點(diǎn)。 在姜泰勒系統(tǒng)使用qcisd和ccsd需要額外小心。ccsd和qcisd的哈密頓矩陣都不是赫米矩陣,所以ccsd和qcisd勢(shì)能面間的交叉點(diǎn)的拓?fù)湫再|(zhì)都不正確。 |
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
專(zhuān)家顧問(wèn) (職業(yè)作家)
地溝油冶煉專(zhuān)家
![]() |
專(zhuān)家經(jīng)驗(yàn): +458 |

木蟲(chóng) (著名寫(xiě)手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +10 | Li李魚(yú) 2026-03-26 | 10/500 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
[考研] 309求調(diào)劑 +15 | 誰(shuí)不是少年 2026-03-29 | 15/750 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +5 | 小小機(jī)靈蟲(chóng) 2026-03-29 | 11/550 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 安徽大學(xué)專(zhuān)碩生物與醫(yī)藥專(zhuān)業(yè)(086000)324分,英語(yǔ)已過(guò)四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂(lè)雞翅 2026-03-26 | 4/200 |
|
|
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|