版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

kangsgo

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 93.6
帖子: 19
在線: 15.9小時(shí)
蟲(chóng)號(hào): 2389350
注冊(cè): 2013-03-30
專(zhuān)業(yè): 蛋白質(zhì)組學(xué)

[求助] gromacs配體力場(chǎng)電荷的問(wèn)題

我先教程在線服務(wù)器PRODRG 2.5的方法建立的力場(chǎng)電荷文件如下：

引用回帖:

; nr    type  resnr resid  atom  cgnr charge    mass

   1    CH2    1  EVO    CAF    1 0.000  13.0190

   2       C    1  EVO    CAE    2 0.023  12.0110

   3       NR    1  EVO    NAC    2 0.043  14.0067

   4       H    1  EVO    HAC    2 -0.066 1.0080

   5       C    1  EVO    CAI    3 0.116  12.0110

   6    CR1    1  EVO    CAP    3 -0.022  12.0110

   7       HC    1  EVO    HAR    3 0.011 1.0080

   8    CR1    1  EVO    CAU    3 -0.022  12.0110

   9       HC    1  EVO    HAY    3 0.011 1.0080

10    CR1    1  EVO    CAT    3 -0.022  12.0110

11       HC    1  EVO    HAX    3 0.011 1.0080

12    CR1    1  EVO    CAO    3 -0.022  12.0110

13       HC    1  EVO    HAQ    3 0.011 1.0080

14       C    1  EVO    CAH    3 -0.036  12.0110

15       C    1  EVO    CAG    3 -0.036  12.0110

16    CH2    1  EVO    CAM    4 0.113  14.0270

17    CH1    1  EVO    CAN    4 0.112  14.0270

18       NL    1  EVO    NAB    4 0.773  14.0067

19       OA    1  EVO    OAA    4 -0.160  15.9994

20    CH2    1  EVO    CAJ    4 0.162  13.0190

21       C    1  EVO    CAL    5 -0.016  12.0110

22    CR1    1  EVO    CAR    5 -0.009  12.0110

23       HC    1  EVO    HAT    5 0.025 1.0080

24    CR1    1  EVO    CAW    6 -0.015  12.0110

25       HC    1  EVO    HA0    6 0.015 1.0080

26    CR1    1  EVO    CAV    6 -0.015  12.0110

27       HC    1  EVO    HAZ    6 0.015 1.0080

28    CR1    1  EVO    CAQ    7 -0.016  12.0110

29       HC    1  EVO    HAS    7 0.013 1.0080

30       C    1  EVO    CAK    7 0.145  12.0110

31       NL    1  EVO    NAD    7 -0.236  14.0067

32    CH3    1  EVO    CAS    7 0.094  15.0350

群里和教程說(shuō)這個(gè)做出來(lái)不準(zhǔn)，用ATB較好，于是我用它做了一個(gè)，如下：

引用回帖:

[ atoms ]
;  nr  type  resnr  resid  atom  cgnr  charge mass total_charge
1 HC 1 SDSG H13 1 0.089 1.0080
2    C 1 SDSG C15 1 -0.127  12.0110
3 HC 1 SDSG H11 1 0.089 1.0080
4 HC 1 SDSG H12 1 0.089 1.0080    ;  0.140
5 NT 1 SDSG    N3 2 -0.530  14.0067
6    C 1 SDSG    C2 2 0.710  12.0110
7 HC 1 SDSG    H2 2 0.035 1.0080
8    C 1 SDSG    C1 2 -0.239  12.0110
9 NT 1 SDSG    N2 2 -0.437  14.0067
10  HS14 1 SDSG    H1 2 0.364 1.0080
11    C 1 SDSG    C3 2 0.097  12.0110    ; -0.000
12    C 1 SDSG    C5 3 0.309  12.0110
13    C 1 SDSG    C4 3 -0.122  12.0110
14    C 1 SDSG C12 3 -0.372  12.0110
15 HC 1 SDSG    H8 3 0.185 1.0080    ;  0.000
16    C 1 SDSG C11 4 -0.139  12.0110
17 HC 1 SDSG    H7 4 0.139 1.0080    ;  0.000
18    C 1 SDSG C16 5 -0.223  12.0110
19 HC 1 SDSG H14 5 0.171 1.0080
20    C 1 SDSG C17 5 -0.068  12.0110
21 HC 1 SDSG H15 5 0.120 1.0080    ;  0.000
22    C 1 SDSG    C7 6 0.552  12.0110
23    C 1 SDSG    C8 6 -0.355  12.0110
24    C 1 SDSG C13 6 -0.383  12.0110
25 HC 1 SDSG    H9 6 0.186 1.0080    ;  0.000
26    C 1 SDSG C14 7 -0.023  12.0110
27 HC 1 SDSG H10 7 0.109 1.0080
28    C 1 SDSG C19 7 -0.216  12.0110
29 HC 1 SDSG H17 7 0.130 1.0080    ;  0.000
30    C 1 SDSG C18 8 -0.012  12.0110
31 HC 1 SDSG H16 8 0.149 1.0080    ;  0.137
32 NT 1 SDSG    N1 9 -0.560  14.0067
33    C 1 SDSG    C6 9 0.704  12.0110
34    O 1 SDSG    O1 9 -0.589  15.9994
35    C 1 SDSG C10 9 -0.019  12.0110
36 HC 1 SDSG    H5 9 0.079 1.0080
37 HC 1 SDSG    H6 9 0.108 1.0080    ; -0.277
38    C 1 SDSG    C9 10 -0.188  12.0110
39 HC 1 SDSG    H3 10 0.096 1.0080
40 HC 1 SDSG    H4 10 0.092 1.0080    ;  0.000
; total charge of the molecule:  -0.000

但是它為何有兩個(gè)的原子數(shù)都不相同·····所以我又用gaussian和ambertools做了一下resp，電荷如下：

引用回帖:

1  DUMM  DU M 0  -1  -2    0.000    .0       .0    .00000
2  DUMM  DU M 1 0  -1    1.449    .0       .0    .00000
3  DUMM  DU M 2 1 0    1.523 111.21    .0    .00000
4  C18 ca M 3 2 1    1.540 111.208  -180.000 -0.055157
5  C19 ca B 4 3 2    1.386 85.241 142.707 -0.250863
6  C14 ca B 5 4 3    1.382 118.685 127.363 -0.013132
7  C8 ca S 6 5 4    1.388 120.843    0.144 -0.295655
8  C6 c    S 7 6 5    1.484 118.878  -177.881  0.775848
9  O1 o    E 8 7 6    1.202 122.075    1.441 -0.589377
10  H10 ha E 6 5 4    1.073 121.147  -179.133  0.138140
11  H17 ha E 5 4 3    1.075 120.678 -52.361  0.151555
12  H16 ha E 4 3 2    1.076 61.050 14.700  0.143398
13  C13 ca M 4 3 2    1.383 124.432 -92.009 -0.348872
14  H9 ha E 13 4 3    1.072 119.445 73.243  0.169709
15  C7 ca M 13 4 3    1.397 119.914  -107.315  0.407232
16  N3 nh M 15  13 4    1.380 122.608 178.980 -0.314454
17  C15 c3 3 16  15  13    1.443 121.147 -10.183 -0.164667
18  H11 h1 E 17  16  15    1.081 109.161  -168.263  0.088033
19  H12 h1 E 17  16  15    1.083 110.363 -49.702  0.088033
20  H13 h1 E 17  16  15    1.089 113.402 71.436  0.088033
21  C2 c3 M 16  15  13    1.440 119.011 159.869  0.137290
22  H2 h2 E 21  16  15    1.083 106.205 157.152  0.153836
23  N1 n    M 21  16  15    1.454 111.980 41.122 -0.499094
24  C10 c3 M 23  21  16    1.462 112.625 176.392  0.087076
25  H5 h1 E 24  23  21    1.085 107.153 -53.056  0.061913
26  H6 h1 E 24  23  21    1.077 108.409  -169.658  0.061913
27  C9 c3 M 24  23  21    1.537 111.448 67.415 -0.224694
28  H3 hc E 27  24  23    1.087 109.239  -162.305  0.111802
29  H4 hc E 27  24  23    1.086 109.717 80.840  0.111802
30  C3 cd M 27  24  23    1.505 109.472 -40.598 -0.034776
31  C1 cc M 30  27  24    1.347 123.077 10.717  0.145019
32  N2 na M 31  30  27    1.378 110.755  -177.270 -0.671235
33  H1 hn E 32  31  30    0.993 124.025  -160.962  0.437444
34  C5 ca M 32  31  30    1.377 107.985 -3.570  0.344809
35  C4 ca M 34  32  31    1.401 108.072    3.182 -0.053032
36  C11 ca M 35  34  32    1.397 119.313 178.515 -0.134838
37  H7 ha E 36  35  34    1.075 120.594  -179.623  0.143002
38  C16 ca M 36  35  34    1.377 118.987    0.732 -0.233578
39  H14 ha E 38  36  35    1.076 119.885 179.756  0.152162
40  C17 ca M 38  36  35    1.402 120.836 -0.091 -0.098966
41  H15 ha E 40  38  36    1.075 119.352 179.829  0.139658
42  C12 ca M 40  38  36    1.377 121.238 -0.370 -0.326357
43  H8 ha E 42  40  38    1.076 121.084  -179.649  0.171039

所以我有點(diǎn)迷茫了···我應(yīng)該如何將其后面生成的電荷給轉(zhuǎn)入進(jìn)去呀？

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有6人回復(fù)
083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
初始318分求調(diào)劑（有工作經(jīng)驗(yàn)）已經(jīng)有3人回復(fù)
二本跨考鄭大材料306英一數(shù)二已經(jīng)有3人回復(fù)
化學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) 已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有4人回復(fù)
265求調(diào)劑已經(jīng)有9人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)

1樓 2015-11-11 23:11:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 kangsgo 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽(yáng)CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 317求調(diào)劑 +8	申子申申 2026-03-19	13/650	2026-03-21 00:09 by 劉國(guó)森
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門(mén)版塊排行榜

kangsgo

[求助] gromacs配體力場(chǎng)電荷的問(wèn)題

» 猜你喜歡

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频