 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2595 | 回復(fù): 12 | ||
[求助]
環(huán)己烷在gromacs中的gromos45a3力場中的itp文件怎么寫？ 已有1人參與
|
|
想用gromacs中的gromos45a3力場算一個環(huán)己烷的體系，但是自己寫了一個環(huán)己烷的itp文件，grompp時在報錯，不知道問題出在哪里？格式不對么？求指教。。。 [ moleculetype ] ;Name nrexcl CYLH 3 [ atoms ] CH2 CH2 0.0 0 CH2 CH2 0.0 0 CH2 CH2 0.0 0 CH2 CH2 0.0 0 CH2 CH2 0.0 0 CH2 CH2 0.0 0 [ bonds ] CH2 CH2 gb_26 CH2 CH2 gb_26 CH2 CH2 gb_26 CH2 CH2 gb_26 CH2 CH2 gb_26 CH2 CH2 gb_26 [ angles ] CH2 CH2 CH2 ga_14 CH2 CH2 CH2 ga_14 CH2 CH2 CH2 ga_14 CH2 CH2 CH2 ga_14 CH2 CH2 CH2 ga_14 CH2 CH2 CH2 ga_14 [ dihedrals ] CH2 CH2 CH2 CH2 gd_17 CH2 CH2 CH2 CH2 gd_17 CH2 CH2 CH2 CH2 gd_17 CH2 CH2 CH2 CH2 gd_17 CH2 CH2 CH2 CH2 gd_17 CH2 CH2 CH2 CH2 gd_17 圖片1.png |
回復(fù)此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- tinker使用oplsaa力場的原子類型定義問題
已經(jīng)有4人回復(fù)
- 請問gromacs導(dǎo)入EAM勢能函數(shù)的問題和lammps導(dǎo)入martini粗�；龅膯栴}
已經(jīng)有4人回復(fù)
- Gromacs怎么使用clayff力場？
已經(jīng)有6人回復(fù)
- Gromacs中添加CHARMM力場
已經(jīng)有41人回復(fù)
- Gromacs粗�；觯ㄗ约簩懀�
已經(jīng)有11人回復(fù)
- 離子液體模擬：如何從別的軟件移植力場到Gromacs（itp文件的含義以及如何寫itp文件）
已經(jīng)有4人回復(fù)
- gromacs 添加新的力場pdb2gmx報錯
已經(jīng)有5人回復(fù)
- Gromacs運行時itp文件語法出錯
已經(jīng)有5人回復(fù)
- gromacs計算時原系統(tǒng)帶電荷，如何添加離子中和電性？
已經(jīng)有8人回復(fù)
- gromacs，用pdb文件生成top文件時遇到的問題
已經(jīng)有4人回復(fù)
- gromacs模擬分子在二氧化硅表面的吸附
已經(jīng)有9人回復(fù)
- gromacs oplsaa 力場中的原子類型
已經(jīng)有6人回復(fù)
- 求Gromacs中Glycam力場數(shù)據(jù)包。。非常感謝
已經(jīng)有5人回復(fù)
- 【求助】gromacs出現(xiàn)了一個問題
已經(jīng)有10人回復(fù)
fangsteel
木蟲 (正式寫手)
- 應(yīng)助: 96 (初中生)
- 金幣: 1213.6
- 散金: 32
- 紅花: 22
- 帖子: 456
- 在線: 192.1小時
- 蟲號: 2248690
- 注冊: 2013-01-18
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2014-05-06 11:28:40
|
我按著格式修改后，mdrun后體系膨脹了，用的是npt系綜，密度只有5.4kg/m3，而環(huán)己烷的是密度是779呢？問題出在哪里呢？求指教 ~~~~(>_<  ~~~~ [moleculetype ] ;Name nrexcl CYLH 3 [ atoms ] ; nr at type res nr residu name at name cg nr charge 1 CH2 1 CYLH C2 1 0.0 14.027 2 CH2 1 CYLH C2 1 0.0 14.027 3 CH2 1 CYLH C2 1 0.0 14.027 4 CH2 1 CYLH C2 1 0.0 14.027 5 CH2 1 CYLH C2 1 0.0 14.027 6 CH2 1 CYLH C2 1 0.0 14.027 [ bonds ] ; i j funct length force.c. 1 2 1 0.1530 7.1500e+06 2 3 1 0.1530 7.1500e+06 3 4 1 0.1530 7.1500e+06 4 5 1 0.1530 7.1500e+06 5 6 1 0.1530 7.1500e+06 6 1 1 0.1530 7.1500e+06 [ angles ] ; i j k funct angle force.c. 1 2 3 1 111.00 530.00 2 3 4 1 111.00 530.00 3 4 5 1 111.00 530.00 4 5 6 1 111.00 530.00 5 6 1 1 111.00 530.00 6 1 2 1 111.00 530.00 [ dihedrals ] ; i j k l funct phi cp mult 1 2 3 4 1 0.000 5.92 3 2 3 4 5 1 0.000 5.92 3 3 4 5 6 1 0.000 5.92 3 4 5 6 1 1 0.000 5.92 3 5 6 1 3 1 0.000 5.92 3 6 1 2 3 1 0.000 5.92 3 |
3樓2014-05-08 08:32:21
fangsteel
木蟲 (正式寫手)
- 應(yīng)助: 96 (初中生)
- 金幣: 1213.6
- 散金: 32
- 紅花: 22
- 帖子: 456
- 在線: 192.1小時
- 蟲號: 2248690
- 注冊: 2013-01-18
- 性別: GG
- 專業(yè): 原子和分子物理
4樓2014-05-08 12:06:32
|
T= 300K, 分子數(shù)目512 這是mdp文件 include = define = integrator = md tinit = 0 dt = 0.001 nsteps = 10000000 comm-mode = Linear nstcomm = 10 nstlog = 1000 nstcalcenergy = -1 nstenergy = 1000 nstxout = 0 nstvout = 0 nstfout = 0 nstxtcout = 1000 xtc-precision = 1000 energygrps = CYLH nstlist = 10 ns_type = grid pbc = xyz rlist = 0.8 coulombtype = cutoff rcoulomb-switch = 0 rcoulomb = 0.8 vdw-type = cutoff rvdw-switch = 0 rvdw = 0.8 pme_order = 4 ewald_rtol = 1e-05 ewald_geometry = 3dc epsilon_surface = 0 optimize_fft = yes Tcoupl = nose-hoover nhchainlength = 1 tc-grps = system tau_t = 0.4 ref_t = 300 ;Pressure coupling Pcoupl = Parrinello-Rahman Pcoupltype = isotropic tau_p = 2.0 compressibility = 4.5e-5 ref_p = 1.01325 andersen_seed = 815131 gen_vel = no gen_temp = 300.0 gen_seed = 173529 constraints = all-bonds |
5樓2014-05-08 16:46:28
fangsteel
木蟲 (正式寫手)
- 應(yīng)助: 96 (初中生)
- 金幣: 1213.6
- 散金: 32
- 紅花: 22
- 帖子: 456
- 在線: 192.1小時
- 蟲號: 2248690
- 注冊: 2013-01-18
- 性別: GG
- 專業(yè): 原子和分子物理
6樓2014-05-08 21:06:36
7樓2014-05-13 14:55:06
fangsteel
木蟲 (正式寫手)
- 應(yīng)助: 96 (初中生)
- 金幣: 1213.6
- 散金: 32
- 紅花: 22
- 帖子: 456
- 在線: 192.1小時
- 蟲號: 2248690
- 注冊: 2013-01-18
- 性別: GG
- 專業(yè): 原子和分子物理
8樓2014-05-13 15:48:01
9樓2014-05-14 10:10:23
fangsteel
木蟲 (正式寫手)
- 應(yīng)助: 96 (初中生)
- 金幣: 1213.6
- 散金: 32
- 紅花: 22
- 帖子: 456
- 在線: 192.1小時
- 蟲號: 2248690
- 注冊: 2013-01-18
- 性別: GG
- 專業(yè): 原子和分子物理
10樓2014-05-14 10:41:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xwwjd 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|