 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 853 | 回復: 2 | |||
草田露銀蟲 (初入文壇)
|
[求助]
關(guān)于力場的編寫 求擅長towhee的大神指導！ 已有1人參與
|
|
打算模擬二氯甲烷的氣液相密度，選用了OPLS-1996中適合二氯甲烷的一些參數(shù)，編寫力場如下： towhee_ff Version 14 Number of Nonbonded Types 2 Potential Type Lennard-Jones Classical Mixrule Geometric Atom Type Number 1 Nonbond Coefficients 0.3800000000E+01 0.5937955880E+02 0.3800000000E+01 0.2968977940E+02 Mass 0.1201070000E+02 Element C Bond Pattern null Base Charge 0.5000000000E+00 Polarizability 0.1000000000E+01 Force Field Name OPLS-1996 Atom Names C C C C Atom Type Number 2 Nonbond Coefficients 0.3400000000E+01 0.1509649800E+03 0.3400000000E+01 0.7548249000E+02 Mass 0.3545300000E+02 Element Cl Bond Pattern null Base Charge -0.2500000000E+00 Polarizability 0.1000000000E+01 Force Field Name OPLS-1996 Atom Names Cl Cl Cl Cl Number of Bonded Terms 1 Bond Type Number 1 Bond Style 1 Bond Coefficients 0.1772000000E+01 Vibration Order single Force Field Name OPLS-1996 Number of Atoms with Same Parameters 1 Atom Names C Cl Number of Angle Terms 1 Angle Type Number 1 Angle Style 1 Angle Coefficients 0.1118000000E+03 0.9057898800E+05 Angle Order wild Force Field Name OPLS-1996 Number of Atoms with Same Parameters 1 Atom Names Cl C Cl Number of Torsion Terms 0 Number of Improper Terms 0 Number of Angle-Angle Terms 0 Number of One-Five Types 0 Number of Bond Increments 0 input文件如下： inputformat 'Towhee' ensemble 'nvt' temperature 324.0d0 nmolty 1 nmolectyp 300 numboxes 2 stepstyle 'cycles' nstep 30000 controlstyle 'equilibration' potentialstyle 'internal' ffnumber 1 ff_filename /home/user1/towhee/towhee-7.0.2/ForceFields/towhee_ff_OPLS-1996 classical_potential 'Lennard-Jones' classical_mixrule 'Geometric ' lshift .false. ltailc .true. rmin 0.8d0 rcut 10.0d0 rcutin 10.0d0 electrostatic_form 'coulomb' coulombstyle 'ewald_fixed_kmax' kalp 5.6d0 kmax 5 dielect 1.0d0 linit .ture. initboxtype 'dimensions' initstyle 'full cbmc' 'full cbmc' initlattice 'simple cubic' 'simple cubic' initmol 240 60 inix iniy iniz 8 8 7 4 4 4 hmatrix 29.90d0 0.0d0 0.0d0 0.0d0 29.90d0 0.0d0 0.0d0 0.0d0 29.90d0 121.82d0 0.0d0 0.0d0 0.0d0 121.82d0 0.0d0 0.0d0 0.0d0 121.82d0 pmvol 0.002d0 pmvlpr 1.0d0 rmvol 0.1d0 tavol 0.5d0 pmcb 0.354d0 pmcbmt 1.0d0 pmall 0.0d0 pmtracm 0.675d0 pmtcmt 1.0d0 rmtrac 0.5d0 tatrac 0.5d0 pmrotate 1.0d0 pmromt 1.0d0 rmrot 0.05d0 tarot 0.5d0 cbmc_formulation 'Martin and Frischknecht 2006' cbmc_setting_style 'Martin and Frischknecht' # input_style 'basic connectivity map' nunit 3 nmaxcbmc 3 lpdbnames F forcefield 'OPLS-1996' charge_assignment 'bond increment' unit ntype qqatom 1 'Cl' -0.250d0 vibration 1 2 improper torsion 0 unit ntype qqatom 2 'C' 0.500d0 vibration 2 1 3 improper torsion 0 unit ntype qqatom 3 'Cl' -0.250d0 vibration 1 2 improper torsion 0 在towhee軟件中運行時，outpu文件中出現(xiàn)這樣的錯誤， Building the input file for molecule type: 1 unit: 1 name:Cl unit: 2 name:C unit: 3 name:Cl parameters not found for bond atoms: 1 and: 2 vib: 1 names: Cl C 希望好心人幫忙指點一二  另附上力場及input文件 |
回復此樓» 猜你喜歡
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有4人回復
-
086502化學工程342求調(diào)劑
已經(jīng)有3人回復
-
求調(diào)劑
已經(jīng)有9人回復
-
352分 化工與材料
已經(jīng)有5人回復
-
085602 化工專碩 338分 求調(diào)劑
已經(jīng)有10人回復
-
0856材料化工調(diào)劑 總分330
已經(jīng)有10人回復
-
330一志愿中國海洋大學 化學工程 085602 有讀博意愿 求調(diào)劑
已經(jīng)有4人回復
-
一志愿哈爾濱工業(yè)大學材料與化工方向336分
已經(jīng)有6人回復
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有6人回復
-
085600 286分 材料求調(diào)劑
已經(jīng)有4人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 一些關(guān)于模擬汽液相平衡的問題（towhee新手）
已經(jīng)有4人回復
- ms分子模擬畫結(jié)構(gòu)時需要考慮單鍵雙鍵和三鍵嗎？
已經(jīng)有16人回復
- MS的力場文件函數(shù)如何調(diào)用
已經(jīng)有5人回復
- Towhee能模擬氣體吸附等溫線嗎？
已經(jīng)有10人回復
- 使用amber力場，一些參數(shù)無法獲得
已經(jīng)有5人回復
- 求用到AMBER force field的開源軟件
已經(jīng)有10人回復
- 怎么自己做力場-及氣體力場的選擇
已經(jīng)有6人回復
- 【討論】分子模擬中力場參數(shù)的調(diào)節(jié)討論
已經(jīng)有7人回復
043114076
木蟲 (正式寫手)
- 應助: 49 (小學生)
- 金幣: 3364.9
- 散金: 52
- 紅花: 23
- 帖子: 788
- 在線: 369.5小時
- 蟲號: 570592
- 注冊: 2008-06-07
- 性別: GG
- 專業(yè): 理論和計算化學
草田露2樓2015-06-23 10:24:23
草田露
銀蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227.1
- 散金: 24
- 帖子: 26
- 在線: 60.3小時
- 蟲號: 3080499
- 注冊: 2014-03-22
- 性別: MM
- 專業(yè): 化工熱力學和基礎數(shù)據(jù)
3樓2015-06-23 11:29:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 草田露 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +12 | minht 2026-03-26 | 12/600 |
|
|---|---|---|---|---|
|
|
[有機交流]
高溫高壓反應求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 070300求調(diào)劑306分 +3 | 26要上岸 2026-03-27 | 3/150 |
|
|
|
[考研] 07化學280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 生物學學碩，一志愿湖南大學，初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
|
[考研] 材料科學與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學學碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-24 | 3/150 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|