版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

yangxinyeah

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 0.7
散金: 43
帖子: 62
在線: 114.6小時
蟲號: 1127193
注冊: 2010-10-20
專業(yè): 化學(xué)

[求助] 在Towhee中用TraPPE-EH力場模擬benzene的VLE曲線已有1人參與

大家好：

最近我用Towhee軟件在模擬苯的氣液平衡數(shù)據(jù)時發(fā)現(xiàn)只能選擇TraPPE-UA力場。TraPPE-EH力場中只定義了CO2和N2。所以我在參數(shù)文件中把Siepmann的苯的TraPPE-EH力場數(shù)據(jù)增加進(jìn)去。并且在input文件中也定義了苯的鏈接信息。但是程序報錯，請問論壇里哪位高手知道這個怎么解決嗎？小女子在此拜謝了。還請大家?guī)臀野伞?br />
報錯信息為：
input_style: advanced connectivity map
nunit:          12
nmaxcbmc:          12
lpdbnames:  F
using the TraPPE-EH  force field
charge_assignment: bond increment
improper_assignment: none
match_style: standard
Building the input file for molecule type:    1
unit: 1 name:Caro
Mismatch of strings in buildmolec
Expected String: unit
Found String: impr
Fatal error, quitting.

我的輸入文件：
inputformat
'Towhee'
ensemble
'nvt'
temperature
520.0d0
nmolty
1
nmolectyp
250
numboxes
2
stepstyle
'cycles'
nstep
20000
printfreq
2000
blocksize
4000
moviefreq
4000
backupfreq
4000
runoutput
'full'
pdb_output_freq
4000
loutdft
.false.
loutlammps
.false.
loutdlpoly
.false.
pressurefreq
10
trmaxdispfreq
200
volmaxdispfreq
200
potentialstyle
'internal'
ffnumber
1
ff_filename
/home/liuxd/towhee-7.1.0/ForceFields/towhee_ff_TraPPE-EH
classical_potential
'Lennard-Jones'
classical_mixrule
'Lorentz-Berthelot'
lshift
.false.
ltailc
.true.
rmin
0.6d0
rcut
14.0d0
rcutin
10.0d0
electrostatic_form
'coulomb'
coulombstyle
'ewald_fixed_kmax'
kalp
5.6d0
kmax
8
dielect
1.0d0
linit
.true.
initboxtype
'dimensions'
initstyle
'full cbmc'
'full cbmc'
initlattice
'simple cubic'
'simple cubic'
initmol
210
40
inix iniy iniz
8 8 7
4 4 4
hmatrix
30.0d0 0.0d0 0.0d0
0.0d0 30.0d0 0.0d0
0.0d0 0.0d0 30.0d0
60.d0 0.0d0 0.0d0
0.0d0  60.d0 0.0d0
0.0d0 0.0d0  60.d0
pmvol
0.002d0
      pmvlpr
      1.0d0
      rmvol
      0.1d0
      tavol
      0.5d0
pmcell
0.0d0
      pmcellpr
      1.0d0
      pmcellpt
      0.5d0
      rmcell
      1.0d0
      tacell
      0.5d0
pm2boxrbswap
0.0d0
      pm2rbswmt
      1.0d0
      pm2rbswpr
      1.0d0
pm2boxcbswap
0.01d0
      pm2cbswmt
      1.0d0
      pm2cbswpr
      1.0d0
pm1boxcbswap
0.0d0
      pm1cbswmt
      1.0d0
pmcb
0.4d0
      pmcbmt
      1.0d0
      pmall
      0.0d0
pmtracm
0.7d0
      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'
#trappe-EH benzene
input_style
'advanced connectivity map'
nunit
12
nmaxcbmc
12
lpdbnames
F
forcefield
'TraPPE-EH'
charge_assignment
'bond increment'
improper_assignment
'none'
match_style
'standard'
unit ntype
1 'Caro'
vibration
3
2 'aromatic'
6 'aromatic'
7 'single'
improper torsion
0
unit ntype
2 'Caro'
vibration
3
1 'aromatic'
3 'aromatic'
8 'single'
improper torsion
0
unit ntype
3 'Caro'
vibration
3
2 'aromatic'
4 'aromatic'
9 'single'
improper torsion
0
unit ntype
4 'Caro'
vibration
3
3 'aromatic'
5 'aromatic'
10 'single'
improper torsion
0
unit ntype
5 'Caro'
vibration
3
4 'aromatic'
6 'aromatic'
11 'single'
improper torsion
0
unit ntype
6 'Caro'
vibration
3
5 'aromatic'
1 'aromatic'
12 'single'
improper torsion
0
unit ntype
7 'Haro'
vibration
1
1 'single'
improper torsion
0
unit ntype
8 'Haro'
vibration
1
2 'single'
improper torsion
0
unit ntype
9 'Haro'
vibration
1
3 'single'
improper torsion
0
unit ntype
10 'Haro'
vibration
1
4 'single'
improper torsion
0
unit ntype
11 'Haro'
vibration
1
5 'single'
improper torsion
0
unit ntype
12 'Haro'
vibration
1
6 'single'
improper torsion
0

修改后的力場文件：
towhee_ff Version
      15
Number of Nonbonded Types
         6
Potential Type
Lennard-Jones
Classical Mixrule
Lorentz-Berthelot
Atom Type Number
         1
Nonbond Coefficients
0.2800000000E+01
0.2700000000E+02
0.0000000000E+00
0.0000000000E+00
Mass
0.1201100000E+02
Element
C
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
C_co2
C_co2
C_co2
none
Atom Type Number
         2
Nonbond Coefficients
0.3050000000E+01
0.7900000000E+02
0.0000000000E+00
0.0000000000E+00
Mass
0.1599900000E+02
Element
O
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
O_co2
O_co2
O_co2
none
Atom Type Number
         3
Nonbond Coefficients
0.3310000000E+01
0.3600000000E+02
0.0000000000E+00
0.0000000000E+00
Mass
0.1400700000E+02
Element
N
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
N_n2
N_n2
N_n2
none
Atom Type Number
         4
Nonbond Coefficients
0.0000000000E+00
0.0000000000E+00
0.0000000000E+00
0.0000000000E+00
Mass
0.0000000000E+00
Element
-
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
COM_n2
COM_n2
COM_n2
none
Atom Type Number
         5
Nonbond Coefficients
0.3600000000E+01
0.3070000000E+02
0.0000000000E+00
0.0000000000E+00
Mass
0.1201100000E+02
Element
C
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
Caro
Caro
Caro
none
Atom Type Number
         6
Nonbond Coefficients
0.2360000000E+01
0.2545000000E+02
0.0000000000E+00
0.0000000000E+00
Mass
0.1007900000E+01
Element
H
Bond Pattern
null
Base Charge
0.0000000000E+00
Polarizability
0.1000000000E+01
Force Field Name
TraPPE-EH
Atom Names
Haro
Haro
Haro
none
Number of Bonded Terms
         4
Bond Type Number
         1
Bond Style
         1
Bond Coefficients
0.1160000000E+01
Vibration Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         1
Atom Names
C_co2    O_co2
Bond Type Number
         2
Bond Style
         1
Bond Coefficients
0.5500000000E+00
Vibration Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         1
Atom Names
N_n2    COM_n2
Bond Type Number
         3
Bond Style
         1
Bond Coefficients
0.1080000000E+01
Vibration Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         1
Atom Names
Caro    Haro
Bond Type Number
         4
Bond Style
         1
Bond Coefficients
0.1392000000E+01
Vibration Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         1
Atom Names
Caro    Caro
Number of Angle Terms
         2
Angle Type Number
         1
Angle Style
         0
Angle Coefficients
0.1800000000E+03
0.1000000000E-04
Angle Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         2
Atom Names
O_co2    C_co2    O_co2
N_n2    COM_n2    N_n2
Angle Type Number
         2
Angle Style
         0
Angle Coefficients
0.1200000000E+03
0.1000000000E-04
Angle Order
wild
Force Field Name
TraPPE-EH
Number of Atoms with Same Parameters
         3
Atom Names
Caro    Caro    Caro
Haro    Caro    Haro
Caro    Caro    Haro
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
         4
Bond Increment Type Number
         1
Bond Increment Value
-0.3500000000E+00
Bond Increment Order
wild
Force Field Name
TraPPE-EH
Atom Names
O_co2    C_co2
Bond Increment Type Number
         2
Bond Increment Value
-0.4820000000E+00
Bond Increment Order
wild
Force Field Name
TraPPE-EH
Atom Names
N_n2    COM_n2
Bond Increment Type Number
         3
Bond Increment Value
0.0000000000E+00
Bond Increment Order
wild
Force Field Name
TraPPE-UA
Atom Names
Caro    Caro
Bond Increment Type Number
         4
Bond Increment Value
-0.9500000000E+00
Bond Increment Order
wild
Force Field Name
TraPPE-UA
Atom Names
Caro    Haro

回復(fù)此樓

» 猜你喜歡

本人考085602 化學(xué)工程專碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2016-04-21 13:34:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yangxinyeah

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 0.7
散金: 43
帖子: 62
在線: 114.6小時
蟲號: 1127193
注冊: 2010-10-20
專業(yè): 化學(xué)

求幫助啊。自己頂一個

贊一下

回復(fù)此樓

2樓2016-04-22 07:05:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

零度不結(jié)冰

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 540.9
散金: 30
紅花: 4
帖子: 117
在線: 112.5小時
蟲號: 1156355
注冊: 2010-11-26
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

towhee等很多程序都是串行讀入，根據(jù)報錯的信息，在對應(yīng)的輸入文件的字段進(jìn)行查找修改

贊一下

回復(fù)此樓

我的態(tài)度決定我的高度。

3樓2016-04-22 11:41:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

992648352

銀蟲 (小有名氣)

應(yīng)助: 10 (幼兒園)
金幣: 38
散金: 61
紅花: 6
帖子: 256
在線: 191.7小時
蟲號: 1249767
注冊: 2011-03-30
性別: GG
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

把所有的
improper torsion
0
去掉。
好好看看教程中關(guān)于improper_assignment=none的情況
http://towhee.sourceforge.net/in ... improper_assignment

贊一下

回復(fù)此樓

4樓2016-04-23 16:17:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhoupe123

木蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 4301.5
散金: 50
帖子: 141
在線: 130.6小時
蟲號: 2313837
注冊: 2013-03-03
性別: GG
專業(yè): 化學(xué)工程及工業(yè)化學(xué)

同用towhee模擬汽液相平衡，可以加個好友交流一下

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

5樓2016-04-25 15:31:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

992648352

銀蟲 (小有名氣)

月只藍(lán): 屏蔽內(nèi)容, 不版塊禁止發(fā)布QQ群信息。 2016-09-26 08:31:28

本帖內(nèi)容被屏蔽

6樓2016-04-26 11:42:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

匿名

用戶注銷 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1014.3
散金: 482
紅花: 13
帖子: 249
在線: 118.8小時
蟲號: 0
注冊: 2014-10-29
性別: MM
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

本帖僅樓主可見

贊一下

7樓2016-07-07 16:46:59

已閱申請模擬EPI 回復(fù)此樓編輯查看我的主頁

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 yangxinyeah 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy

24小時熱門版塊排行榜

[求助] 在Towhee中用TraPPE-EH力場模擬benzene的VLE曲線 已有1人參與

» 猜你喜歡

【答案】應(yīng)助回帖

[求助] 在Towhee中用TraPPE-EH力場模擬benzene的VLE曲線已有1人參與