| 8 | 1/1 | 返回列表 |
| 查看: 1213 | 回復(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) (著名寫手)
|
這個(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ì)算了 |
木蟲(chóng) (小有名氣)
|
忘記說(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。 |
木蟲(chóng) (著名寫手)
|
波函數(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ì)都不正確。 |
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
專家顧問(wèn) (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

木蟲(chóng) (著名寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +14 | 李李不服輸 2026-03-26 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +6 | 真名有冰 2026-03-29 | 7/350 |
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲(chóng)085600 2026-03-27 | 8/400 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲(chóng)不懶了 2026-03-27 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|