应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 关于力场的编写 求擅长towhee的大神指导！
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 858  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

草田露

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[求助] 关于力场的编写 求擅长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文件
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-06-22 11:31:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

【答案】应助回帖

感谢参与，应助指数 +1
我感觉应该还是你的ｆｆ文件内容有错误。这个文件对于空格也是敏感的，不能多也不能少。建议使用程序自带的ｆｆ文件，如果要修改也要先在文件中做一步一步的修改。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

草田露
2楼2015-06-23 10:24:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

草田露

主管区长

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
送红花一朵
引用回帖:
2楼: Originally posted by 043114076 at 2015-06-23 10:24:23
我感觉应该还是你的ｆｆ文件内容有错误。这个文件对于空格也是敏感的，不能多也不能少。建议使用程序自带的ｆｆ文件，如果要修改也要先在文件中做一步一步的修改。

谢谢师兄！
回复此楼
3楼2015-06-23 11:29:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 草田露 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +12 墨墨漠 2026-03-25 12/600 2026-04-01 01:42 by 1018329917
[考研] 一志愿：西北大学，英一数一408-284分求调剂 +7 12.27 2026-03-27 7/350 2026-03-31 21:59 by lbsjt
[考研] 0856材料化工调剂 总分330 +17 zhubinhao 2026-03-27 17/850 2026-03-31 18:56 by 小张做实验
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 304求调剂 +8 素年祭语 2026-03-31 11/550 2026-03-31 18:20 by 无际的草原
[考研] 085600，321分求调剂 +7 大馋小子 2026-03-31 7/350 2026-03-31 18:02 by 无际的草原
[考研] 307分求调剂 +6 (o~o) 2026-03-31 6/300 2026-03-31 17:22 by 唐沐儿
[考研] 化学工程085602 305分求调剂 +28 RichLi_ 2026-03-25 36/1800 2026-03-31 14:56 by JourneyLucky
[考研] 土木304求调剂 +3 顶级擦擦 2026-03-31 3/150 2026-03-31 13:59 by 西京学院招办
[考研] 085600 295分求调剂 +13 W55j 2026-03-30 15/750 2026-03-31 13:29 by 王亮_大连医科大
[考研] 334求调剂 +7 Trying] 2026-03-31 7/350 2026-03-31 12:33 by 无际的草原
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 一志愿郑大材料工程290求调剂 +12 Youth_ 2026-03-30 12/600 2026-03-31 03:34 by 蒙奇奇521
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 324求调剂 +9 hanamiko 2026-03-26 11/550 2026-03-30 14:27 by JourneyLucky
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 0856调剂 +5 求求让我有书读� 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
信息提示
关闭
请填处理意见
关闭