| 3 | 1/1 | 返回列表 |
| 查看: 1516 | 回復(fù): 2 | ||
[求助]
【新人求教】結(jié)構(gòu)優(yōu)化使用分子力場方法是否合適? 已有2人參與
|
|
本行是有機,最近剛剛接觸計算,計算的都是小分子的芳香共軛分子。 在使用gaussian的時候,我一般是把結(jié)構(gòu)優(yōu)化好的.gjf文件載入程序。 稍微查閱了一下相關(guān)資料,別人貌似都是用RHF或者B3LYP之類的方法對結(jié)構(gòu)進(jìn)行優(yōu)化,但是我自己用的方法都是用gaussview的clean功能找到分子的最高對稱性構(gòu)型,然后用Chem3d的MM2做分子力場方法的結(jié)構(gòu)優(yōu)化。 我自己的考慮是使用半經(jīng)驗方法得到的幾何結(jié)構(gòu)應(yīng)該要比從頭算精確,尤其對只含C,H,O,N之類輕原子的小分子體系來說。 我的這種想法是否正確呢?方法上是否有需要改進(jìn)的呢? 希望各位大神不吝賜教! |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
即便是CHNO的體系,較好的半經(jīng)驗方法(如PM6)至多也只比HF這種意義的從頭算優(yōu)化的結(jié)構(gòu)更準(zhǔn)確,比起較好的GGA級別的DFT的結(jié)構(gòu)更準(zhǔn)確一般都是不可能的,就算更準(zhǔn)確也只是巧合,尤其是可靠性遠(yuǎn)不及從頭算(諸如有一次用PM7優(yōu)化磷脂,居然把尾巴上的雙鍵的平面優(yōu)化成扭曲的) chem3D的MMFF94比MM2要強,不建議用MM2。MMFF94優(yōu)化出的結(jié)構(gòu)算是不錯的,如果要求更準(zhǔn)確,可以以其作為初始結(jié)構(gòu)再用諸如M06-2X/6-31G**或更高級別優(yōu)化進(jìn)一步結(jié)構(gòu)。 |
木蟲 (著名寫手)
|
Sob說的已經(jīng)比較清楚了。 我現(xiàn)在將我的做法說一下吧,我只是用GV 和chemdraw/3D用來建模,對結(jié)構(gòu)不太把握的時候可以用3D/MM2來大概的優(yōu)化出結(jié)構(gòu)(個人認(rèn)為這個優(yōu)化結(jié)果的可信度不高),然后再放到Gaussian里面去優(yōu)化。 當(dāng)然有時也會發(fā)現(xiàn)MM2尋找的最小點結(jié)構(gòu)與實際情況相差太遠(yuǎn),尤其是研究弱相互作用的時候,MM2結(jié)果是不可信的。 另外,LZ也提到了相關(guān)文獻(xiàn)上利用HF和B3LYP進(jìn)行優(yōu)化,建議LZ也這樣做,并將相關(guān)文獻(xiàn)引上。既然你做的是小分子,我建議你不要用HF方法,可以用Sob提到的方法與基組或者文獻(xiàn)中的B3LYP方法,基組參考最新文獻(xiàn)。個人更推薦Sob提到的基組與方法。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|