| 5 | 1/1 | 返回列表 |
| 查看: 1423 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
樂石||木蟲 (正式寫手)
|
[交流]
LAMMPS, CHARMM-gen force field, topology file builder, PDB to DATA 已有1人參與
|
||
|
可能有時候大家會覺得有機(jī)物的data結(jié)構(gòu)文件搭建比較費時,所以會用一些軟件進(jìn)行從pdb到data的自動轉(zhuǎn)換。但是,MS也好,其他軟件也好,總有找不到信息的時候,因為ff數(shù)據(jù)庫不可能涵蓋所有原子種類。Amber會稍好一些, 因為它只關(guān)心原子和幾個原子相連,不太關(guān)心和誰相連,但是這樣可能會降低準(zhǔn)確度。如果你用charmm ff的話,可以參照我寫的python代碼。力場參數(shù)文件為par_all36_cgenff.prm。 https://github.com/yangl9/topobuilder/ 該代碼也只能夠?qū)崿F(xiàn)從pdb到data轉(zhuǎn)換的半自動化。因為任何時候,par_all36_cgenff.prm是不能給出所有原子種類信息的,這種參數(shù)測試工作太耗時了。這時我們需要做一些比較靠譜的近似,cgenff就可以給出和你想要的結(jié)構(gòu)差不多的而且存在于par_all36_cgenff.prm文件中的立場參數(shù)。但是問題在于cgenff有時給出的信息也非常不靠譜,比如penalty超過50.( https://cgenff.paramchem.org/commonFiles/newsEvents.php)這時我們需要手動更正一些信息,這就是半自動化。 我的代碼能夠自動從par_all36_cgenff.prm搜取已存在的信息,而對于那些不存在的,允許手動輸入。當(dāng)你輸入鍵的信息之后,代碼會根據(jù)所輸入鍵(1-2)的信息去找其所對應(yīng)的角的力場參數(shù)(1-2-3);如果你輸入了角(1-2-3)的信息,代碼會根據(jù)所輸入角的信息去找相關(guān)dihedral的力場參數(shù)(1-2-3-4)。(如果存在一個dihedral對應(yīng)多行信息的情況,也需要去手動更正。) 有任何問題歡迎隨時反饋,郵箱請見代碼首行。 18.03.2019 |
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-08 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
[考研] 301求調(diào)劑 +5 | 一二LV 2026-03-07 | 5/250 |
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
[考研] 物理學(xué)求調(diào)劑 +3 | sx宋 2026-03-05 | 3/150 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 能動297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研] 325求調(diào)劑 +5 | 學(xué)家科 2026-03-04 | 5/250 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 一志愿東北大學(xué)材料專碩328,求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|