| 24小時(shí)熱門(mén)版塊排行榜 |
| 9 | 1/1 | 返回列表 |
| 查看: 3817 | 回復(fù): 8 | ||
lujunyan1118金蟲(chóng) (小有名氣)
|
[求助]
用高斯QM/MM算頻率和QST3 已有2人參與
|
|
各位大蝦,小弟有個(gè)問(wèn)題。 我想用高斯的QM/MM算一個(gè)反應(yīng),高層大約100個(gè)原子,低層約500個(gè)殘基。 高層用B3LYP/6-31G*。 我優(yōu)化了反應(yīng)物和產(chǎn)物,并猜測(cè)了一個(gè)過(guò)渡態(tài),想用QST3找過(guò)渡態(tài)。但是算的時(shí)候總是說(shuō)內(nèi)存不夠,需要3506MW。可電腦沒(méi)這么多內(nèi)存。為何QST3要如此之大的內(nèi)存?我減少低層的殘基有效嗎? 另外,用QM/MM算頻率的時(shí)候也經(jīng)常說(shuō)內(nèi)存不夠,但是只用QM時(shí)不會(huì)出現(xiàn)這個(gè)問(wèn)題。 是不是QM/MM的時(shí)候,低層的頻率也計(jì)算了呢?能否讓高斯只計(jì)算高層的頻率? 謝謝各位! |
金蟲(chóng) (小有名氣)
|
一般也就用QM/MM算大分子體系核心體系構(gòu)型; 不知道你的計(jì)算命令具體設(shè)置,如果用dft/AMBER肯定算不了,計(jì)算量過(guò)大。 可以參考紐約大學(xué)Yingkai Zhang,他們組做地比較有心得 |

金蟲(chóng) (小有名氣)
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
|
低層設(shè)得太大了,尤其是算頻率,這極其困難。頻率是極其耗時(shí)也是相當(dāng)耗內(nèi)存的。 高層也偏大,只需要把關(guān)鍵性的原子包含就夠了,而低層把外面與之相接的一層殘基包含進(jìn)去也就夠了,不要把整個(gè)蛋白都當(dāng)成低層。 |
金蟲(chóng) (小有名氣)
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
金蟲(chóng) (小有名氣)
|
謝謝,我明白了。 還有一個(gè)問(wèn)題想請(qǐng)教下,我看文獻(xiàn)中有很多用QM/MM算大體系的催化機(jī)理,他們好像說(shuō)的是用Full enzyme,沒(méi)說(shuō)把中心區(qū)域摳出來(lái)。而且他們也算了反應(yīng)物,TS,產(chǎn)物的頻率。他們其實(shí)是只取了一部分蛋白算頻率,而沒(méi)有說(shuō)嗎? 我要通過(guò)比較反應(yīng)物,TS,產(chǎn)物能量來(lái)算活化能,那這三者的能量也必須用摳出來(lái)的那一部分結(jié)構(gòu)來(lái)算能量? |
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
|
QMMM在具體實(shí)現(xiàn)上,不同的人、不同程序做法可能很不同。ONIOM算是高斯(也包括Q-Chem)的一種特殊的實(shí)現(xiàn)方式,并不是常規(guī)的QMMM做法。 全酶的振動(dòng)分析高斯做不動(dòng),但不代表其它程序做不動(dòng)(比如專(zhuān)門(mén)做生物體系的amber就可以做得不錯(cuò))。 也可能很多文獻(xiàn)在優(yōu)化的時(shí)候用全酶,僅在做振動(dòng)分析獲得熱力學(xué)校正能時(shí)只考慮QM區(qū)域(雖然沒(méi)明確說(shuō)出來(lái)),這樣處理完全是合理的。 |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|
[考研] 生物學(xué)調(diào)劑招人!! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 材料專(zhuān)碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|