 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1803 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
yangxinyeah銅蟲 (小有名氣)
|
[求助]
在Towhee中用TraPPE-EH力場(chǎng)模擬benzene的VLE曲線 已有1人參與
|
|
|
大家好： 最近我用Towhee軟件在模擬苯的氣液平衡數(shù)據(jù)時(shí)發(fā)現(xiàn)只能選擇TraPPE-UA力場(chǎng)。TraPPE-EH力場(chǎng)中只定義了CO2和N2。所以我在參數(shù)文件中把Siepmann的苯的TraPPE-EH力場(chǎng)數(shù)據(jù)增加進(jìn)去。并且在input文件中也定義了苯的鏈接信息。但是程序報(bào)錯(cuò)，請(qǐng)問論壇里哪位高手知道這個(gè)怎么解決嗎？小女子在此拜謝了。還請(qǐng)大家?guī)臀野伞?br /> 報(bào)錯(cuò)信息為： 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 修改后的力場(chǎng)文件： 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ù)此樓» 猜你喜歡
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有9人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有31人回復(fù)
-
070300化學(xué)專業(yè)279調(diào)劑
已經(jīng)有3人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
085404 22408 315分
已經(jīng)有5人回復(fù)
-
一志愿a區(qū)211，085601-307分求調(diào)劑
已經(jīng)有10人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
290分調(diào)劑求助
已經(jīng)有7人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有3人回復(fù)
-
考研調(diào)劑求助
已經(jīng)有4人回復(fù)
零度不結(jié)冰
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 540.9
- 散金: 30
- 紅花: 4
- 帖子: 117
- 在線: 112.5小時(shí)
- 蟲號(hào): 1156355
- 注冊(cè): 2010-11-26
- 性別: GG
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
3樓2016-04-22 11:41:20
yangxinyeah
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 0.7
- 散金: 43
- 帖子: 62
- 在線: 114.6小時(shí)
- 蟲號(hào): 1127193
- 注冊(cè): 2010-10-20
- 專業(yè): 化學(xué)
2樓2016-04-22 07:05:40
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
把所有的 improper torsion 0 去掉。 好好看看教程中關(guān)于improper_assignment=none的情況 http://towhee.sourceforge.net/in ... improper_assignment |
4樓2016-04-23 16:17:33
zhoupe123
木蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 4301.5
- 散金: 50
- 帖子: 141
- 在線: 130.6小時(shí)
- 蟲號(hào): 2313837
- 注冊(cè): 2013-03-03
- 性別: GG
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
5樓2016-04-25 15:31:24
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 11ggg 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +9 | 柳溪邊 2026-03-26 | 9/450 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。。� +5 | Acaciad 2026-03-25 | 5/250 |
|
|
[考研] 315求調(diào)劑 +5 | akie... 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
|
[考研] 323分 食品與營(yíng)養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +12 | Youth_ 2026-03-30 | 12/600 |
|
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考博] 26申博 +3 | 加油沖�。� 2026-03-26 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|