24小时热门版块排行榜    

查看: 1289  |  回复: 3

fish637

铜虫 (初入文坛)

[求助] prodrg生成小分子不同文件时的问题 已有1人参与

想用Gromacs模拟小分子,在prodrg中输入分子式,如CH3(CH2)25CH3,结果生成小分子坐标文件和itp时,不是以单个原子为单元的,而是以CH3或CH2为一组的,提供的itp如下。有什么办法能得到atom命名为C、H的,type类型为C、H的itp,而不是以原子组CH3或者CH2的形式。请大神们帮帮忙~~
[ moleculetype ]
; Name nrexcl
LIG      3

[ atoms ]
;   nr      type  resnr resid  atom  cgnr   charge     mass
     1       CH3     1  LIG     CAZ     1   -0.023  15.0350   
     2       CH2     1  LIG     CAX     1    0.001  14.0270   
     3       CH2     1  LIG     CAW     1    0.004  14.0270   
     4       CH2     1  LIG     CAT     1    0.004  14.0270   
     5       CH2     1  LIG     CAR     1    0.004  14.0270   
     6       CH2     1  LIG     CAP     1    0.005  14.0270   
     7       CH2     1  LIG     CAN     1    0.005  14.0270   
     8       CH2     1  LIG     CAM     2    0.000  14.0270   
     9       CH2     1  LIG     CAJ     2    0.000  14.0270   
    10       CH2     1  LIG     CAH     2   -0.001  14.0270   
    11       CH2     1  LIG     CAF     2    0.001  14.0270   
    12       CH2     1  LIG     CAD     3    0.000  14.0270   
    13       CH2     1  LIG     CAB     3    0.000  14.0270   
    14       CH2     1  LIG     CAA     3   -0.001  14.0270   
    15       CH2     1  LIG     CAC     3    0.001  14.0270   
    16       CH2     1  LIG     CAE     4    0.000  14.0270   
    17       CH2     1  LIG     CAG     4    0.000  14.0270   
    18       CH2     1  LIG     CAI     4    0.000  14.0270   
    19       CH2     1  LIG     CAK     5    0.000  14.0270   
    20       CH2     1  LIG     CAL     5    0.000  14.0270   
    21       CH2     1  LIG     CAO     6    0.004  14.0270   
    22       CH2     1  LIG     CAQ     6    0.003  14.0270   
    23       CH2     1  LIG     CAS     6    0.003  14.0270   
    24       CH2     1  LIG     CAU     6    0.003  14.0270   
    25       CH2     1  LIG     CAV     6    0.005  14.0270   
    26       CH2     1  LIG     CAY     6    0.004  14.0270   
    27       CH3     1  LIG     CBA     6   -0.022  15.0350   

[ bonds ]
; ai  aj  fu    c0, c1, ...
   2   1   2    0.153   7150000.0    0.153   7150000.0 ;   CAX  CAZ   
   2   3   2    0.153   7150000.0    0.153   7150000.0 ;   CAX  CAW   
   3   4   2    0.153   7150000.0    0.153   7150000.0 ;   CAW  CAT   
   4   5   2    0.153   7150000.0    0.153   7150000.0 ;   CAT  CAR   
   5   6   2    0.153   7150000.0    0.153   7150000.0 ;   CAR  CAP   
   6   7   2    0.153   7150000.0    0.153   7150000.0 ;   CAP  CAN   
   7   8   2    0.153   7150000.0    0.153   7150000.0 ;   CAN  CAM   
   8   9   2    0.153   7150000.0    0.153   7150000.0 ;   CAM  CAJ   
   9  10   2    0.153   7150000.0    0.153   7150000.0 ;   CAJ  CAH   
  10  11   2    0.153   7150000.0    0.153   7150000.0 ;   CAH  CAF   
  11  12   2    0.153   7150000.0    0.153   7150000.0 ;   CAF  CAD   
  12  13   2    0.153   7150000.0    0.153   7150000.0 ;   CAD  CAB   
  13  14   2    0.153   7150000.0    0.153   7150000.0 ;   CAB  CAA   
  14  15   2    0.153   7150000.0    0.153   7150000.0 ;   CAA  CAC   
  15  16   2    0.153   7150000.0    0.153   7150000.0 ;   CAC  CAE   
  16  17   2    0.153   7150000.0    0.153   7150000.0 ;   CAE  CAG   
  17  18   2    0.153   7150000.0    0.153   7150000.0 ;   CAG  CAI   
  18  19   2    0.153   7150000.0    0.153   7150000.0 ;   CAI  CAK   
  19  20   2    0.153   7150000.0    0.153   7150000.0 ;   CAK  CAL   
  20  21   2    0.153   7150000.0    0.153   7150000.0 ;   CAL  CAO   
  21  22   2    0.153   7150000.0    0.153   7150000.0 ;   CAO  CAQ   
  22  23   2    0.153   7150000.0    0.153   7150000.0 ;   CAQ  CAS   
  23  24   2    0.153   7150000.0    0.153   7150000.0 ;   CAS  CAU   
  24  25   2    0.153   7150000.0    0.153   7150000.0 ;   CAU  CAV   
  25  26   2    0.153   7150000.0    0.153   7150000.0 ;   CAV  CAY   
  26  27   2    0.153   7150000.0    0.153   7150000.0 ;   CAY  CBA   

[ pairs ]
; ai  aj  fu    c0, c1, ...
   1   4   1                                           ;   CAZ  CAT   
   2   5   1                                           ;   CAX  CAR   
   3   6   1                                           ;   CAW  CAP   
   4   7   1                                           ;   CAT  CAN   
   5   8   1                                           ;   CAR  CAM   
   6   9   1                                           ;   CAP  CAJ   
   7  10   1                                           ;   CAN  CAH   
   8  11   1                                           ;   CAM  CAF   
   9  12   1                                           ;   CAJ  CAD   
  10  13   1                                           ;   CAH  CAB   
  11  14   1                                           ;   CAF  CAA   
  12  15   1                                           ;   CAD  CAC   
  13  16   1                                           ;   CAB  CAE   
  14  17   1                                           ;   CAA  CAG   
  15  18   1                                           ;   CAC  CAI   
  16  19   1                                           ;   CAE  CAK   
  17  20   1                                           ;   CAG  CAL   
  18  21   1                                           ;   CAI  CAO   
  19  22   1                                           ;   CAK  CAQ   
  20  23   1                                           ;   CAL  CAS   
  21  24   1                                           ;   CAO  CAU   
  22  25   1                                           ;   CAQ  CAV   
  23  26   1                                           ;   CAS  CAY   
  24  27   1                                           ;   CAU  CBA   

[ angles ]
; ai  aj  ak  fu    c0, c1, ...
   1   2   3   2    109.5       520.0    109.5       520.0 ;   CAZ  CAX  CAW   
   2   3   4   2    109.5       520.0    109.5       520.0 ;   CAX  CAW  CAT   
   3   4   5   2    109.5       520.0    109.5       520.0 ;   CAW  CAT  CAR   
   4   5   6   2    109.5       520.0    109.5       520.0 ;   CAT  CAR  CAP   
   5   6   7   2    109.5       520.0    109.5       520.0 ;   CAR  CAP  CAN   
   6   7   8   2    109.5       520.0    109.5       520.0 ;   CAP  CAN  CAM   
   7   8   9   2    109.5       520.0    109.5       520.0 ;   CAN  CAM  CAJ   
   8   9  10   2    109.5       520.0    109.5       520.0 ;   CAM  CAJ  CAH   
   9  10  11   2    109.5       520.0    109.5       520.0 ;   CAJ  CAH  CAF   
  10  11  12   2    109.5       520.0    109.5       520.0 ;   CAH  CAF  CAD   
  11  12  13   2    109.5       520.0    109.5       520.0 ;   CAF  CAD  CAB   
  12  13  14   2    109.5       520.0    109.5       520.0 ;   CAD  CAB  CAA   
  13  14  15   2    109.5       520.0    109.5       520.0 ;   CAB  CAA  CAC   
  14  15  16   2    109.5       520.0    109.5       520.0 ;   CAA  CAC  CAE   
  15  16  17   2    109.5       520.0    109.5       520.0 ;   CAC  CAE  CAG   
  16  17  18   2    109.5       520.0    109.5       520.0 ;   CAE  CAG  CAI   
  17  18  19   2    109.5       520.0    109.5       520.0 ;   CAG  CAI  CAK   
  18  19  20   2    109.5       520.0    109.5       520.0 ;   CAI  CAK  CAL   
  19  20  21   2    109.5       520.0    109.5       520.0 ;   CAK  CAL  CAO   
  20  21  22   2    109.5       520.0    109.5       520.0 ;   CAL  CAO  CAQ   
  21  22  23   2    109.5       520.0    109.5       520.0 ;   CAO  CAQ  CAS   
  22  23  24   2    109.5       520.0    109.5       520.0 ;   CAQ  CAS  CAU   
  23  24  25   2    109.5       520.0    109.5       520.0 ;   CAS  CAU  CAV   
  24  25  26   2    109.5       520.0    109.5       520.0 ;   CAU  CAV  CAY   
  25  26  27   2    109.5       520.0    109.5       520.0 ;   CAV  CAY  CBA   

[ dihedrals ]
; ai  aj  ak  al  fu    c0, c1, m, ...
   4   3   2   1   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAT  CAW  CAX  CAZ   
   5   4   3   2   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAR  CAT  CAW  CAX   
   6   5   4   3   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAP  CAR  CAT  CAW   
   7   6   5   4   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAN  CAP  CAR  CAT   
   8   7   6   5   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAM  CAN  CAP  CAR   
   9   8   7   6   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAJ  CAM  CAN  CAP   
  10   9   8   7   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAH  CAJ  CAM  CAN   
  11  10   9   8   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAF  CAH  CAJ  CAM   
  12  11  10   9   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAD  CAF  CAH  CAJ   
  13  12  11  10   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAB  CAD  CAF  CAH   
  14  13  12  11   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAA  CAB  CAD  CAF   
  15  14  13  12   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAC  CAA  CAB  CAD   
  16  15  14  13   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAE  CAC  CAA  CAB   
  17  16  15  14   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAG  CAE  CAC  CAA   
  18  17  16  15   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAI  CAG  CAE  CAC   
  19  18  17  16   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAK  CAI  CAG  CAE   
  20  19  18  17   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAL  CAK  CAI  CAG   
  21  20  19  18   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAO  CAL  CAK  CAI   
  22  21  20  19   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAQ  CAO  CAL  CAK   
  23  22  21  20   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAS  CAQ  CAO  CAL   
  24  23  22  21   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAU  CAS  CAQ  CAO   
  25  24  23  22   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAV  CAU  CAS  CAQ   
  26  25  24  23   1      0.0    5.9 3      0.0    5.9 3 ; dih   CAY  CAV  CAU  CAS   
  27  26  25  24   1      0.0    5.9 3      0.0    5.9 3 ; dih   CBA  CAY  CAV  CAU
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhzj2011

新虫 (初入文坛)

应该跟力场有关吧,你看看prodrg能不能换成全原子力场
2楼2014-10-07 16:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish637

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by zhzj2011 at 2014-10-07 16:13:17
应该跟力场有关吧,你看看prodrg能不能换成全原子力场

prodrg使用时,没有选力场的项,并且使用说明有指出它是按组生成的,不知道楼主有么有好的办法
3楼2014-10-07 20:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wgz_tju

金虫 (小有名气)

【答案】应助回帖

★ ★
fish637: 金币+2, 有帮助 2014-10-21 10:55:46
引用回帖:
3楼: Originally posted by fish637 at 2014-10-07 20:44:42
prodrg使用时,没有选力场的项,并且使用说明有指出它是按组生成的,不知道楼主有么有好的办法...

它默认用gromos力场
4楼2014-10-20 15:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fish637 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见