 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 848 | 回復: 2 | ||
草田露銀蟲 (初入文壇)
|
[求助]
關于力場的編寫 求擅長towhee的大神指導！ 已有1人參與
|
|
打算模擬二氯甲烷的氣液相密度，選用了OPLS-1996中適合二氯甲烷的一些參數，編寫力場如下： 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文件中出現這樣的錯誤， 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文件 |
回復此樓» 猜你喜歡
一志愿北京化工大學 070300 學碩 336分 求調劑
已經有4人回復
-
352求調劑
已經有3人回復
-
一志愿東華大學化學070300，求調劑
已經有8人回復
-
277材料科學與工程080500求調劑
已經有7人回復
-
317求調劑
已經有18人回復
-
293求調劑
已經有5人回復
-
280分求調劑 一志愿085802
已經有7人回復
-
0854電子信息求調劑
已經有3人回復
-
263求調劑
已經有4人回復
-
石河子大學（211、雙一流）碩博研究生長期招生公告
已經有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 一些關于模擬汽液相平衡的問題（towhee新手）
已經有4人回復
- ms分子模擬畫結構時需要考慮單鍵雙鍵和三鍵嗎？
已經有16人回復
- MS的力場文件函數如何調用
已經有5人回復
- Towhee能模擬氣體吸附等溫線嗎？
已經有10人回復
- 使用amber力場，一些參數無法獲得
已經有5人回復
- 求用到AMBER force field的開源軟件
已經有10人回復
- 怎么自己做力場-及氣體力場的選擇
已經有6人回復
- 【討論】分子模擬中力場參數的調節(jié)討論
已經有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è): 化工熱力學和基礎數據
3樓2015-06-23 11:29:33
相關版塊跳轉
我要訂閱樓主 草田露 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 315分，誠求調劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 289材料與化工（085600）B區(qū)求調劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 287求調劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 319求調劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 291 求調劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學071000，求調劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 280求調劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 求調劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
|
[考研] 0703化學297求調劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 336求調劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 0805材料320求調劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 22 350 本科985求調劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 346求調劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 301求調劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 085410人工智能專碩317求調劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 312求調劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 085601求調劑 +4 | Du.11 2026-03-16 | 4/200 |
|