| 5 | 1/1 | 返回列表 |
| 查看: 2016 | 回復(fù): 8 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】請問哪種方法的模擬結(jié)果更加可信? 已有6人參與
|
||||
|
在模擬氨基酸分子在水溶液中最穩(wěn)定的分子構(gòu)象時,文獻中常用的方法有:DFT,ab initio, B3LYP/6-31+G, HF, MP2, MS, Discover3等,不考慮計算時間的話,請問這些方法中哪些方法的計算結(jié)果最可信(最準確)? DFT 和MS, Discover3 相比,誰更準確? 本人對模擬不太懂,多謝指點啊。 ![]() [ Last edited by wuli8 on 2010-4-26 at 13:22 ] |
第一性原理 |
木蟲 (正式寫手)
|
hi,做分子模擬,比如模擬氨基酸的穩(wěn)定構(gòu)象,從熱力學(xué)上講,就是要找到自由能最低的構(gòu)象。 對于一個真空中的單分子,能量就是自由能,因此只要搜索它各種可能的構(gòu)象(靜態(tài)),然后看看哪個能量最低即可。 在溶液中,溶質(zhì)構(gòu)象受到溶劑的影響,穩(wěn)定構(gòu)象可能不一樣,自由能需要通過動態(tài)模擬確定,這就是分子動力學(xué)MD(Molecular dynamics)。分子動力學(xué)模擬的關(guān)鍵也是計算能量,相比較靜態(tài)的方法,動態(tài)的,對溶液體相(Bulk solution)的模擬,需要精確計算氨基酸溶質(zhì)和溶劑的分子間相互作用能。 因此,無論哪種方法,能量計算的精度決定了計算結(jié)果的好壞。 能量計算分為考慮電子結(jié)構(gòu)的(從頭算HF,MP2,密度泛函DFT, 統(tǒng)稱第一性原理理論),不考慮電子結(jié)構(gòu)的(純經(jīng)典分子力學(xué)MM,, Molecular Mechanics: 比如CHARMM, AMBER); 還有一些半經(jīng)驗的(AM1, PM3),介于二者之間,但是也是考慮電子結(jié)構(gòu)的。 HF和DFT是兩大類的從頭算能量計算方法,它們通過計算電子的分布得出能量。他們的區(qū)別在于基本的思路不一樣,HF從分子軌道(原子軌道的線性組合)導(dǎo)出能量信息,而DFT通過電子云密度導(dǎo)出能量信息。 前面有朋友說MP2最準確,的確是這樣的。它是一種后HF(Post HF)方法。因為HF計算沒有充分考慮電子相關(guān)(correlation, 簡單說是電子電子相互作用),因此MP2在HF的結(jié)果上進行修正(MP2是二級微擾方法,還有MP4等高級微擾方法)。MP2是上述方法中最為耗時的計算。而DFT是計算精度與速度都比較適中的的方法。 關(guān)于B3LYP/631+G*, 這是一種基于DFT的方法,B3LYP指的是電子交換與相關(guān)的一種處理方法。631+G*是基組,即求解Schrodinger方程所用的基函數(shù),它決定了計算的精度。還有6311++G(d,p)等等,一般來說,基組越高越好。 所有基于電子結(jié)構(gòu)的能量計算都很慢,不適合于大規(guī)模的分子動力學(xué)模擬(也有從頭算分子動力學(xué),大多基于DFT),因此,一般分子動力學(xué)模擬都采用分子力學(xué)的方法計算能量。相比較HF, DFT, 分子力學(xué)(又叫分子力場force field)把分子看成是用彈簧連起來的帶電小球,描述小球間相互作用的函數(shù)叫勢函數(shù),又叫分子力場force field,一個分子力場一般包括一套原子類型的定義,一套勢函數(shù),以及不同原子類型的參數(shù)。適用于生物體系模擬的分子力場有CHARMM, AMBER等。 GROMACS是一個進行分子動力學(xué)模擬的軟件,你需要用它來運行模擬,包括搭建你的模擬所用的盒子,選定所用的力場,設(shè)定模擬時間,分析模擬結(jié)果等等。 Gaussian是進行電子結(jié)構(gòu)計算所用的軟件,可以進行HF, DFT, MP2等各種電子結(jié)構(gòu)計算,也能進行一些基于分子力學(xué)的計算,但是gaussian不適合于進行分子動力學(xué)模擬計算。 如果想做基于電子結(jié)構(gòu)的分子動力學(xué)模擬(第一性原理分子動力學(xué)),你需要使用VASP, ABINIT, DMOL, CPMD, cp2k 等軟件,它們一般基于DFT理論。為了加速計算,它們一般把內(nèi)層電子簡化為一個球,只計算價層電子的結(jié)構(gòu)。即使這樣,第一性原理分子動力學(xué)往往需要大規(guī)模并行計算,非常昂貴。 以上是我的一些理解,供初學(xué)者參考。不對的地方,請大家指教。 [ Last edited by ChemiAndy on 2010-4-27 at 03:38 ] |
木蟲 (正式寫手)
Dr.

木蟲 (正式寫手)
Dr.

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招08考數(shù)學(xué) +5 | laoshidan 2026-03-20 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 6/300 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 一志愿南大,0703化學(xué),分數(shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|