| 5 | 1/1 | 返回列表 |
| 查看: 2029 | 回復: 8 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】請問哪種方法的模擬結果更加可信? 已有6人參與
|
||||
|
在模擬氨基酸分子在水溶液中最穩(wěn)定的分子構象時,文獻中常用的方法有:DFT,ab initio, B3LYP/6-31+G, HF, MP2, MS, Discover3等,不考慮計算時間的話,請問這些方法中哪些方法的計算結果最可信(最準確)? DFT 和MS, Discover3 相比,誰更準確? 本人對模擬不太懂,多謝指點啊。 ![]() [ Last edited by wuli8 on 2010-4-26 at 13:22 ] |
第一性原理 |
木蟲 (正式寫手)
Dr.

木蟲 (正式寫手)
Dr.

木蟲 (正式寫手)
|
hi,做分子模擬,比如模擬氨基酸的穩(wěn)定構象,從熱力學上講,就是要找到自由能最低的構象。 對于一個真空中的單分子,能量就是自由能,因此只要搜索它各種可能的構象(靜態(tài)),然后看看哪個能量最低即可。 在溶液中,溶質構象受到溶劑的影響,穩(wěn)定構象可能不一樣,自由能需要通過動態(tài)模擬確定,這就是分子動力學MD(Molecular dynamics)。分子動力學模擬的關鍵也是計算能量,相比較靜態(tài)的方法,動態(tài)的,對溶液體相(Bulk solution)的模擬,需要精確計算氨基酸溶質和溶劑的分子間相互作用能。 因此,無論哪種方法,能量計算的精度決定了計算結果的好壞。 能量計算分為考慮電子結構的(從頭算HF,MP2,密度泛函DFT, 統(tǒng)稱第一性原理理論),不考慮電子結構的(純經(jīng)典分子力學MM,, Molecular Mechanics: 比如CHARMM, AMBER); 還有一些半經(jīng)驗的(AM1, PM3),介于二者之間,但是也是考慮電子結構的。 HF和DFT是兩大類的從頭算能量計算方法,它們通過計算電子的分布得出能量。他們的區(qū)別在于基本的思路不一樣,HF從分子軌道(原子軌道的線性組合)導出能量信息,而DFT通過電子云密度導出能量信息。 前面有朋友說MP2最準確,的確是這樣的。它是一種后HF(Post HF)方法。因為HF計算沒有充分考慮電子相關(correlation, 簡單說是電子電子相互作用),因此MP2在HF的結果上進行修正(MP2是二級微擾方法,還有MP4等高級微擾方法)。MP2是上述方法中最為耗時的計算。而DFT是計算精度與速度都比較適中的的方法。 關于B3LYP/631+G*, 這是一種基于DFT的方法,B3LYP指的是電子交換與相關的一種處理方法。631+G*是基組,即求解Schrodinger方程所用的基函數(shù),它決定了計算的精度。還有6311++G(d,p)等等,一般來說,基組越高越好。 所有基于電子結構的能量計算都很慢,不適合于大規(guī)模的分子動力學模擬(也有從頭算分子動力學,大多基于DFT),因此,一般分子動力學模擬都采用分子力學的方法計算能量。相比較HF, DFT, 分子力學(又叫分子力場force field)把分子看成是用彈簧連起來的帶電小球,描述小球間相互作用的函數(shù)叫勢函數(shù),又叫分子力場force field,一個分子力場一般包括一套原子類型的定義,一套勢函數(shù),以及不同原子類型的參數(shù)。適用于生物體系模擬的分子力場有CHARMM, AMBER等。 GROMACS是一個進行分子動力學模擬的軟件,你需要用它來運行模擬,包括搭建你的模擬所用的盒子,選定所用的力場,設定模擬時間,分析模擬結果等等。 Gaussian是進行電子結構計算所用的軟件,可以進行HF, DFT, MP2等各種電子結構計算,也能進行一些基于分子力學的計算,但是gaussian不適合于進行分子動力學模擬計算。 如果想做基于電子結構的分子動力學模擬(第一性原理分子動力學),你需要使用VASP, ABINIT, DMOL, CPMD, cp2k 等軟件,它們一般基于DFT理論。為了加速計算,它們一般把內層電子簡化為一個球,只計算價層電子的結構。即使這樣,第一性原理分子動力學往往需要大規(guī)模并行計算,非常昂貴。 以上是我的一些理解,供初學者參考。不對的地方,請大家指教。 [ Last edited by ChemiAndy on 2010-4-27 at 03:38 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 292求調劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 材料求調劑 +8 | @taotao 2026-03-21 | 8/400 |
|
|
[考研]
材料學碩,求調劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 求調劑,一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 324求調劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 284求調劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 341求調劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 材料科學與工程 317求調劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 297求調劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 尋找調劑 +5 | 倔強芒? 2026-03-21 | 8/400 |
|
|
[考研] 化學調劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
[考研] 考研調劑 +6 | 來好運來來來 2026-03-21 | 7/350 |
|
|
[考研] 機械學碩總分317求調劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 網(wǎng)絡空間安全0839招調劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 0854電子信息求調劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 一志愿北化315 求調劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 276求調劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 求調劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|