| 查看: 933 | 回复: 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文件 |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有80人回复
-
感觉是下周放榜了
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一些关于模拟汽液相平衡的问题（towhee新手）
已经有4人回复
- ms分子模拟画结构时需要考虑单键双键和三键吗？
已经有16人回复
- MS的力场文件函数如何调用
已经有5人回复
- Towhee能模拟气体吸附等温线吗？
已经有10人回复
- 使用amber力场，一些参数无法获得
已经有5人回复
- 求用到AMBER force field的开源软件
已经有10人回复
- 怎么自己做力场-及气体力场的选择
已经有6人回复
- 【讨论】分子模拟中力场参数的调节讨论
已经有7人回复
草田露
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 227.1
- 散金: 24
- 帖子: 26
- 在线: 60.3小时
- 虫号: 3080499
- 注册: 2014-03-22
- 性别: MM
- 专业: 化工热力学和基础数据
送红花一朵 3楼2015-06-23 11:29:33
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
2楼2015-06-23 10:24:23
25