| 查看: 1795 | 回復(fù): 12 | ||
dxyan木蟲 (正式寫手)
|
[求助]
做單分子結(jié)構(gòu)優(yōu)化的困惑,求達(dá)人解答。 已有1人參與
|
|
以前做過小分子構(gòu)型優(yōu)化,是通過調(diào)節(jié)鍵長和鍵角,設(shè)計一系列初始結(jié)構(gòu),然后分別計算構(gòu)型優(yōu)化和頻率計算,然后通過能量+頻率來判斷最優(yōu)結(jié)構(gòu)。而對于雙分子體系,根據(jù)可能的作用位點來擺放不同的初始結(jié)構(gòu),也是通過能量+頻率來判斷最優(yōu)結(jié)構(gòu)。然而現(xiàn)在要算結(jié)構(gòu)較為復(fù)雜的大分子,對于這個單分子的優(yōu)化有些束手無策。不知道該如何設(shè)置初始的構(gòu)型。有以下幾個問題還請達(dá)人多多指教。 1 。對于結(jié)構(gòu)優(yōu)化是不是計算后四個收斂標(biāo)準(zhǔn)都是YES,就肯定是最優(yōu)結(jié)構(gòu). 2 。對于結(jié)構(gòu)復(fù)雜的分子從哪些角度設(shè)計初始構(gòu)型,由于是學(xué)化工的,化學(xué)功底薄弱,這方面有些無從下手。 3。另外看到有些帖子講確定最優(yōu)結(jié)構(gòu)可以進(jìn)行勢能面掃面,這個用剛性掃描好還是柔性掃描好?具體通過gaussview怎么設(shè)置? |
木蟲 (正式寫手)
|
1.四個Yes也可能只是極小值,而不是最小值,只是相對優(yōu)化的結(jié)構(gòu),從來沒有最優(yōu)化的說法,只有在相應(yīng)的方法下得到的優(yōu)化結(jié)構(gòu) 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
|
2.設(shè)計分子的初始結(jié)構(gòu)可以通過參考文獻(xiàn)獲得相應(yīng)的鍵長,角度信息,使得初始結(jié)構(gòu)更接近于實驗值,這樣,優(yōu)化的時間會縮短,優(yōu)化的結(jié)果可能也會更理想 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
|
3.聽你的意思你不是專門做計算的,那通常就沒有必要做勢能面掃描,要知道,勢能面掃描所需要的工作量和得到的結(jié)果的意義不成正比,甚至,很多時候,外行人在不知道該算什么的情況下,一味的就勢能面掃描,一副工作量很大的樣子,這是有些滑稽的,實際上,只要你的方法,基組恰當(dāng),優(yōu)化得到的結(jié)果就已經(jīng)夠用了,我不知道樓主是否知道什么是勢能面掃描,但我可以簡單的告訴你,越復(fù)雜的分子,做這個就越困難,結(jié)果的意義甚至越小,一句話,直接優(yōu)化就可以了 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
|
那是不是也就是說4個YES就表示是在勢能面上是能量比較低的點,如果是極小值得話是否應(yīng)該也達(dá)到優(yōu)化要求了?是要在在相應(yīng)的方法下得到一個優(yōu)化的結(jié)構(gòu),可是這個優(yōu)化結(jié)構(gòu)如何確定就是可以接受的構(gòu)型?因為不同的初始結(jié)構(gòu)算出來可能都是一個可以收斂但是構(gòu)型不同的結(jié)構(gòu)。我的主要目的是想得到不同分子的優(yōu)化結(jié)構(gòu)(比如3個不同的分子),再分別和第四個分子算相互作用能,NOB等,來比較相互作用的強(qiáng)弱。如果前面3個單分子的結(jié)構(gòu)都不是優(yōu)化結(jié)構(gòu)的話,這樣和后面第四個分子算出來的相互作用的結(jié)果不是可能不正確,這樣的話不是無法比較? |
木蟲 (正式寫手)
|
如果是文獻(xiàn)沒有算過的分子那該如何設(shè)計初始結(jié)構(gòu)呢?有些文獻(xiàn)上也沒給出鍵長和角度信息。我這次算的分子別人做過相似的結(jié)構(gòu),就是母核是2個環(huán)狀的結(jié)構(gòu),連了一個長的脂肪鏈。別人優(yōu)化的結(jié)構(gòu)中長的脂肪鏈?zhǔn)且粋偏向環(huán)的蜷曲結(jié)構(gòu)。我要算的分子不同處在于這個長脂肪鏈上多了3個雙鍵,其他都是一樣的。結(jié)果按照他的優(yōu)化構(gòu)型加上雙鍵優(yōu)化后,脂肪鏈完全展開。而且第一次優(yōu)化結(jié)果全部是4個YES,不知道為什么?其他不同的初始構(gòu)型我也不知道怎么擺;仡^看看自己以前算的小分子,設(shè)計了不同初始構(gòu)型,然后分別計算構(gòu)型優(yōu)化和頻率計算,然后通過能量+頻率來判斷最優(yōu)結(jié)構(gòu),F(xiàn)在發(fā)現(xiàn)當(dāng)時我所選的那些優(yōu)化構(gòu)型,收斂結(jié)構(gòu)居然都是4個NO(當(dāng)時還不懂這點),所以現(xiàn)在非常困惑。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
通常的做法是每個分子單獨要先優(yōu)化,然后再算多個分子,因為你在計算相互作用能的時候也要用到單個分子的能量,另外,你在計算多個分子的時候,實際上也會對每個分子進(jìn)行優(yōu)化。也有的作者是這樣做的,將多個分子中的單個分子提取出來,固定左邊,算單點能獲得單獨分子的能量 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
送紅花一朵 |
你所說的通過能量和頻率來查看是否會優(yōu)化完畢,我不知道你說的是不是查看虛頻,你之前的算的結(jié)構(gòu)可能有虛頻,如果你算的不是過度態(tài),那這樣的結(jié)果是不合理的,必須消除虛頻 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 5/250 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|