 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1390 | 回復(fù): 4 | |||
[求助]
一些關(guān)于模擬汽液相平衡的問題（towhee新手） 已有1人參與
|
|
本人最近在學(xué)習(xí)Towhee，嘗試自己模擬簡單的純物質(zhì)。出現(xiàn)了不少問題，想在這里請教一下大家。希望高手們解答！ 1、自己寫的input文件，charge_assignment這塊我寫的是bond-increment,但是運(yùn)行的時(shí)候老是說找不到參數(shù) parameters not found for bond increment atoms 1 and 2 vib 1 names [CH ] [CH ] 下面是力場文件相應(yīng)的部分（格式應(yīng)該不會錯，我照著example寫的）。 Bond Increment Type Number 1 Bond Increment Value 0.0000000000E+00 Bond Increment Order aromatic Force Field Name TraPPE-UA Atom Names CH CH 2、后來這個問題得不到解決，因?yàn)樗肽M的物質(zhì)不考慮電荷，所以就把”charge_assignment“這塊改成了”none“。然后運(yùn)行，但是模擬結(jié)果 的汽液密度最終都差不多一樣了，這樣的話就是不分相了。后來嘗試改變盒子的大小，但是都得到的不分相。 “盒子里面的分子數(shù)目都是50/250，盒子的大小根據(jù)實(shí)際的文獻(xiàn)密度給的（大致）�！� 3、還有像請問下是不是各種move的概率問題，我大致給的概率是轉(zhuǎn)動和平動分別是0.35,0.3。體積給的是0.002、交換給的是0.003、pmcb給的是0.29。 |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有10人回復(fù)
-
289求調(diào)劑
已經(jīng)有15人回復(fù)
-
291 求調(diào)劑
已經(jīng)有7人回復(fù)
-
274求調(diào)劑
已經(jīng)有14人回復(fù)
-
北京林業(yè)大學(xué)碩導(dǎo)招生廣告
已經(jīng)有3人回復(fù)
-
309求調(diào)劑
已經(jīng)有5人回復(fù)
-
292求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 運(yùn)用Aspen Plus進(jìn)行分析汽液相平衡關(guān)系
已經(jīng)有3人回復(fù)
- 液相色譜平衡問題
已經(jīng)有17人回復(fù)
- 液相色譜平衡問題
已經(jīng)有11人回復(fù)
- 板式精餾塔設(shè)計(jì)中汽液相平衡曲線擬合
已經(jīng)有5人回復(fù)
- 請教島津高效液相儀平衡流動相以及換流動相時(shí)注意的問題
已經(jīng)有3人回復(fù)
- 年產(chǎn)5000噸的醋酸甲酯生產(chǎn)工藝過程中最小回流比 塔板數(shù) 再沸器熱負(fù)荷等的計(jì)算
已經(jīng)有9人回復(fù)
- fluent可以模擬氣液相反應(yīng)嗎？
已經(jīng)有9人回復(fù)
- Aspen Plus 軟件中氣液平衡問題
已經(jīng)有9人回復(fù)
- 液相色譜柱平衡時(shí)間過長
已經(jīng)有11人回復(fù)
- 工藝開發(fā)中蒸汽壓和氣液平衡數(shù)據(jù)的具體作用
已經(jīng)有3人回復(fù)
- aspen能否解決靜態(tài)非平衡關(guān)系的熱力學(xué)和PVT關(guān)系
已經(jīng)有10人回復(fù)
- monte carlo 計(jì)算汽液
已經(jīng)有6人回復(fù)
- 液相預(yù)熱、柱平衡問題
已經(jīng)有7人回復(fù)
- 液液相平衡中如何表示物料守恒
已經(jīng)有11人回復(fù)
- 汽液平衡數(shù)據(jù)怎么算
已經(jīng)有15人回復(fù)
- 液相鬼峰 每天平衡很久，求原因
已經(jīng)有18人回復(fù)
- Gibbs系綜和NPT系綜的區(qū)別
已經(jīng)有4人回復(fù)
- ASPEN 回歸含鹽汽液相平衡數(shù)據(jù)
已經(jīng)有7人回復(fù)
- 液相平衡時(shí)總是會出現(xiàn)一些小細(xì)峰……
已經(jīng)有4人回復(fù)
- 關(guān)于液相梯度洗脫的平衡的問題
已經(jīng)有19人回復(fù)
- 使用towhee進(jìn)行GEMC模擬相平衡的MC運(yùn)動定義問題..
已經(jīng)有15人回復(fù)
- 【求助】液相在平衡時(shí)，基線一直不平。 怎么辦？
已經(jīng)有13人回復(fù)
- 求助一些汽液相平衡數(shù)據(jù)？
已經(jīng)有1人回復(fù)
- 【討論】液相平衡時(shí)間
已經(jīng)有20人回復(fù)
- 【交流】液相色譜上樣前平衡問題
已經(jīng)有10人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
992648352(月只藍(lán)代發(fā)): 金幣+5, 鼓勵交流！ 2014-06-25 13:56:43
感謝參與，應(yīng)助指數(shù) +1
992648352(月只藍(lán)代發(fā)): 金幣+5, 鼓勵交流！ 2014-06-25 13:56:43
|
把力場文件貼上來吧，目測是你的格式有問題，比如 type number的1，前面的空格數(shù)可能不對。自己去原力場里面看看它的位置。 汽液相分不開，先試試減小汽液相盒子體積，你的汽液相盒子大小分別是多少？溫度？溫度有沒有接近臨界溫度。 另外你的cut值設(shè)的多少？可以嘗試減小截?cái)喟霃街怠?br /> 如果還不行，請?zhí)斎牒土鑫募纯础? |
2樓2014-06-22 20:21:42
|
恩， towhee_ff Version 14 Number of Nonbonded Types 2 Potential Type Lennard-Jones Classical Mixrule Lorentz-Berthelot Atom Type Number 1 Nonbond Coefficients 0.4680000000E+01 0.1000000000E+02 0.0000000000E+00 0.0000000000E+00 Mass 0.1301910000E+02 Element C Bond Pattern sp2 Base Charge 0.0000000000E+00 Polarizability 0.1000000000E+01 Force Field Name TraPPE-UA Atom Names CH CH CH CH Atom Type Number 2 Nonbond Coefficients 0.3600000000E+01 0.1800000000E+03 0.0000000000E+00 0.0000000000E+00 Mass 0.3206500000E+02 Element S Bond Pattern null Base Charge 0.0000000000E+00 Polarizability 0.1000000000E+01 Force Field Name TraPPE-UA Atom Names Scc(thio) Scc(thio) Scc(thio) Scc(thio) Number of Bonded Terms 2 Bond Type Number 1 Bond Style 1 Bond Coefficients 0.1710000000E+01 Vibration Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 1 Atom Names Scc(thio) CH Bond Type Number 2 Bond Style 1 Bond Coefficients 0.1400000000E+01 Vibration Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 1 Atom Names CH CH Number of Angle Terms 3 Angle Type Number 1 Angle Style 0 Angle Coefficients 0.1115000000E+03 0.1000000000E-04 Angle Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 1 Atom Names CH CH Scc(thio) Angle Type Number 2 Angle Style 0 Angle Coefficients 0.0922000000E+03 0.1000000000E-04 Angle Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 1 Atom Names CH Scc(thio) CH Angle Type Number 3 Angle Style 0 Angle Coefficients 0.1124000000E+03 0.1000000000E-04 Angle Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 1 Atom Names CH CH CH Number of Torsion Terms 1 Torsion Type Number 1 Torsion Style 8 One-Four Nonbond Logical T One-Four Coulombic Scaling 0.5000000000E+00 Torsion Coefficients Torsion Order wild Force Field Name TraPPE-UA Number of Atoms with Same Parameters 5 Atom Names Scc(thio) CH CH CH CH CH CH CH CH CH CH Scc(thio) CH CH Scc(thio) CH CH Scc(thio) CH CH Number of Improper Terms 0 Number of Angle-Angle Terms 0 Number of One-Five Types 0 Number of Bond Increments 2 Bond Increment Type Number 1 Bond Increment Value 0.0000000000E+00 Bond Increment Order aromatic Force Field Name TraPPE-UA Atom Names CH Scc(thio) Bond Increment Type Number 2 Bond Increment Value 0.0000000000E+00 Bond Increment Order aromatic Force Field Name TraPPE-UA Atom Names CH CH 輸入文件： inputformat 'Towhee' ensemble 'nvt' temperature 450.00d0 nmolty 1 nmolectyp 300 numboxes 2 stepstyle 'cycles' nstep 5000 printfreq 100 blocksize 1000 moviefreq 2000 backupfreq 2000 runoutput 'full' pdb_output_freq 1000 pressurefreq 10 trmaxdispfreq 100 volmaxdispfreq 100 potentialstyle 'internal' ffnumber 1 ff_filename /home/lc/towhee-7.0.4/ForceFields/towhee_ff_TraPPE-UA classical_potential 'Lennard-Jones' classical_mixrule 'Lorentz-Berthelot' lshift .false. ltailc .true. rmin 1.0d0 rcut 14.0d0 rcutin 10.0d0 electrostatic_form 'none' linit .TRUE. initboxtype 'dimensions' initstyle 'full cbmc' 'full cbmc' initlattice 'simple cubic' 'simple cubic' initmol 50 250 inix iniy iniz 3 5 4 7 6 6 hmatrix 50.0d0 0.0d0 0.0d0 0.0d0 50.0d0 0.0d0 0.0d0 0.0d0 50.0d0 30.0d0 0.0d0 0.0d0 0.0d0 30.0d0 0.0d0 0.0d0 0.0d0 30.0d0 pmvol 0.00020d0 pmvlpr 1.0d0 rmvol 0.1d0 tavol 0.5d0 pm2boxcbswap 0.015d0 pm2cbswmt 1.0d0 pm2cbswpr 1.0d0 pmcb 0.3d0 pmcbmt 1.0d0 pmall 1.0d0 pmtracm 0.65d0 pmtcmt 1.0d0 rmtrac 0.5d0 tatrac 0.5d0 pmrotate 1.0d0 pmromt 1.0d0 rmrot 1.0d0 tarot 0.5d0 cbmc_formulation 'Martin and Frischknecht 2006' cbmc_setting_style 'Martin and Frischknecht' #trappe-ua Thiophene input_style 'basic connectivity map' nunit 5 nmaxcbmc 5 lpdbnames .FALSE. forcefield 'TraPPE-UA' charge_assignment 'bond increment' unit ntype 1 'Scc(thio)' vibration 2 2 5 improper torsion 0 unit ntype 2 'CH' vibration 2 1 3 improper torsion 0 unit ntype 3 'CH' vibration 2 2 4 improper torsion 0 unit ntype 4 'CH' vibration 2 3 5 improper torsion 0 unit ntype 5 'CH' vibration 2 1 4 improper torsion 0 |
3樓2014-06-23 14:29:51
|
溫度為450°，臨界溫度579°。 我模擬的是純噻吩。 老是出現(xiàn)下面的問題？？ using the TraPPE-UA force field charge_assignment: bond increment Building the input file for molecule type: 1 unit: 1 name:Scc(thio) unit: 2 name:CH unit: 3 name:CH unit: 4 name:CH unit: 5 name:CH parameters not found for bond increment atoms 1 and 2 vib 1 names [Scc(thio) ][CH ] 希望你能幫我看看，萬分感謝�。� |
4樓2014-06-23 14:34:10
5樓2014-06-25 21:10:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 992648352 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
材料調(diào)劑
 5+3
|
想要一壺桃花水 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +8 | Mr. Z 2026-03-25 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +6 | 呼呼？~+123456 2026-03-24 | 6/300 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +3 | 小魚愛有機(jī) 2026-03-25 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +4 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+7
|
WWW西西弗斯 2026-03-24 | 7/350 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研]
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|