 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 2114 | 回復: 6 | ||||
趙紅霞鐵桿木蟲 (著名寫手)
|
[求助]
聚合物的pdb文件生成psf文件用于NAMD動力學計算 已有1人參與
|
|
想用全原子做聚合物（聚電解質(zhì)）在溶液里的構(gòu)象，因為體系大，materials studio算起來慢，想在NAMD里算。在生成psf文件時遇到問題，請各位支招。 1. 先在materials studio構(gòu)建好了聚合物分子，導出為pdb格式。 2. 下一步生成psf文件，不知道哪種方法好用，我暫且先選擇在VMD里借助psfgen轉(zhuǎn)化。pgn文件如下： package require psfgen topology par_all27_prot_lipid_na.inp segment U {pdb PCE1_period30nm11.pdb} coordpdb PCE1_period30nm11.pdb U guesscoord writepdb PCE1_period30nm11.pdb writepsf PCE1_period30nm11.psf 問題來了，因為charmm或者amber力場主要針對蛋白和核酸體系，里面規(guī)定了標準的殘基名，根據(jù)殘基來規(guī)定殘基里每個原子的力場參數(shù)。 我的polymer的殘基名無法識別。請問如何解決？ PS：雖然charmm或者amber力場主要針對蛋白和核酸體系，但是我看到文獻里仍然有人拿它們做別的polymer。那么psf如何生成呢？ |
回復此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有16人回復
-
304求調(diào)劑
已經(jīng)有6人回復
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復
-
求材料調(diào)劑
已經(jīng)有8人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復
-
一志愿華中科技大學，080502，354分求調(diào)劑
已經(jīng)有4人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有6人回復
-
085410人工智能專碩317求調(diào)劑（0854都可以）
已經(jīng)有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
趙紅霞
鐵桿木蟲 (著名寫手)
- 應助: 73 (初中生)
- 金幣: 6197.9
- 散金: 141
- 紅花: 33
- 帖子: 1541
- 在線: 700小時
- 蟲號: 1544094
- 注冊: 2011-12-20
- 性別: MM
- 專業(yè): 理論和計算化學
2樓2015-12-17 22:57:38
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗: +49 - 模擬EPI: 2
- 應助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 分子模擬
3樓2015-12-18 01:35:44
趙紅霞
鐵桿木蟲 (著名寫手)
- 應助: 73 (初中生)
- 金幣: 6197.9
- 散金: 141
- 紅花: 33
- 帖子: 1541
- 在線: 700小時
- 蟲號: 1544094
- 注冊: 2011-12-20
- 性別: MM
- 專業(yè): 理論和計算化學
4樓2015-12-18 07:29:07
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗: +49 - 模擬EPI: 2
- 應助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計算化學
- 管轄: 分子模擬
|
我不太明白你的意思。lammps里當然可以用charmm力場。可是你的體系里分子的類型如果和charmm力場不匹配，你怎么調(diào)用charmm力場參數(shù)來做模擬? 我不清楚你是怎么生成lammps輸入文件的。如果是用charmm力場，一般通過charmm2lammps.pl這個腳本來生成輸入文件。即便使用這個腳本，你也需要charmm格式的topology和parameter文件，以及psf文件。如果你的pdb文件里原子類型和charmm力場topology文件里的原子類型根本不匹配，是根本不可能生成輸入文件的。 你現(xiàn)在的問題就是，你pdb文件里原子的類型和charmm力場的topology里原子還有分子類型不匹配，所以無法生成psf文件。 首先，你要搞清楚你要用什么力場，如果你要用charmm力場，那就把你的pdb文件里原子類型改成charmm的原子類型。 其次，你要確定charmm的topology文件里是否有你想模擬的分子。如果沒有，你應該參照charmm topology文件的格式，自己建一個新的topology文件，文件中應該包含你要模擬的高分子的原子質(zhì)量，成鍵，鍵角，二面角等信息。并注意topology文件里的原子名和residue name與你的pdb文件相匹配。至于有沒有腳本可以自動生成topology文件，我不太清楚，請自行查找。 |
5樓2015-12-18 08:57:34
趙紅霞
鐵桿木蟲 (著名寫手)
- 應助: 73 (初中生)
- 金幣: 6197.9
- 散金: 141
- 紅花: 33
- 帖子: 1541
- 在線: 700小時
- 蟲號: 1544094
- 注冊: 2011-12-20
- 性別: MM
- 專業(yè): 理論和計算化學
6樓2015-12-18 20:08:48
王飯飯_
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 129
- 散金: 10
- 帖子: 41
- 在線: 6.5小時
- 蟲號: 28174891
- 注冊: 2022-01-16
- 性別: MM
- 專業(yè): 計算機軟件
7樓2023-02-06 20:27:49
相關版塊跳轉(zhuǎn)
我要訂閱樓主 Andsharia 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研]
 271材料工程求調(diào)劑
+4
|
.6lL 2026-03-18 | 4/200 |
|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研]
0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|