| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1217 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
chembetsey木蟲 (小有名氣)
|
[求助]
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é)果也略有差異。這是否正常? |
木蟲 (著名寫手)
|
波函數(shù)在交叉點(diǎn)附近不穩(wěn)定。建議限制波函數(shù)對(duì)稱性,分別指定電子態(tài)的對(duì)稱性為A'和A''進(jìn)行兩次優(yōu)化。限制對(duì)稱性后的波函數(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ì)都不正確。 |
木蟲 (著名寫手)
|
這個(gè)分子的計(jì)算并不是很容易!它的基態(tài)在c3v下是簡(jiǎn)并的,因?yàn)閖ahn teller效應(yīng)而分裂,產(chǎn)生三個(gè)對(duì)稱性破缺的極小點(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ì)算了 |
木蟲 (小有名氣)
|
忘記說(shuō)了,我也算了CCSD(T),沒(mé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ì)稱性下做的。我也曾試著打破對(duì)稱性,還是有同樣的問(wèn)題。 有一點(diǎn)我覺(jué)得奇怪的是,頻率計(jì)算完之后給出的Maximum Force, RMS Force等一系列值都無(wú)比大,盡管我構(gòu)型優(yōu)化的時(shí)候還用了opt=tight。 |
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)279求調(diào)劑 +4 | 哈哈哈^_^ 2026-03-31 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-30 | 4/200 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 287求調(diào)劑 +14 | land xuxu 2026-03-26 | 14/700 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。。 +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|