24小时热门版块排行榜    

查看: 2426  |  回复: 8

阿龙6678

银虫 (初入文坛)

[求助] Material Studio的Sorption package模拟与实验不符 请问该怎么调参数

最近在用sorption来做模拟,我做的是Xylene在MOF上吸附的模拟,实验值得出xylene是可以进入孔道被吸附的,但是在模拟中xylene根本进不去。请问我该调什么参数呢?刚开始学习做模拟不熟悉,望各位大侠赐教!

我做的过程是:
1. Load解出来的MOF晶体cif,然后用Build-Symetry-Supercell 把晶体变成2*2*3的P1格子
2. Load 用MS画好并clean的p-xylene结构
3. Sorption 计算部分 所有的参数我都放在附图里面了。
4. 计算信息:
    Task                            : Fixed pressure
    Monte Carlo method              : Configurational bias
    Random number seed              :  1307782553
   ---- Fixed pressure calculation parameters ----
   Equilibration steps             :  1000000
   Production steps                :  10000000
   ---- Energy parameters ----
  Forcefield                      : Universal
  Electrostatic terms:
  Summation method              : Ewald & Group
   Accuracy                      : 1e-005 kcal/mol
   Charge group cutoff           : 18.5 A
   Buffer width                  : 0.5 A
  van der Waals terms:
   Summation method              : Atom based
   Truncation method             : Cubic spline
   Cutoff distance               : 18.5 A
   Spline width                  : 1 A
   Buffer width                  : 0.5 A
   Sample interval                 : 25 steps
   ---- Configurational Bias Monte Carlo method parameters ----
|                  | Ratio       | Amplitude   |
|                  |                |                   |
|__________________________________________
| Exchange   |       3.000 |                  |
| Regrowth   |       2.000 |                  |
| Twist         |       1.000 | 5 deg         |
| Rotate       |       1.000 | 5 deg         |
| Translate   |       1.000 | 1 A            |

| Sorbate     | Head bias   | Torsion     | Torsion     |
|                  | steps         | twisters    | bias steps |
|________________________________________________________
| p-Xylene    |          10     |           0    |                |

---- Fixed pressure calculation ----
Framework charge                :  0.000 e per cell
| Sorbate     | Charge      | Fugacity    |
|                  |                 | kPa           |
________________________________________
| p-Xylene    |       0.000 |     101.000 |
Temperature                     :  1.000 K
Total Fugacity                  :  101.000 kPa

---- Monte Carlo steps ----
|Sorbate  | Loading | Creation | Destruction| Regrowth | Rotation| Translation|C/D Ratio|
|             |              | steps     | steps          | steps       | steps     |Steps        |
|__________________________________________________________________________
| p-Xylene|       0    |              |                    |                |              |               |              |
|Accepted|             |           0 |           0       |           0    |           0 |           0  |          0  |
|Attempted|           | 2142880|  2142520     |    2856647 |1428268 |  142968  |           0 |
|  Ratio      |           |     0.000 |       0.000     |       0.000 |     0.000 |      0.000 |             |

---- Loading ----
| Sorbate     | Average     | Minimum     | Maximum     |
|                  | loading      | loading        | loading        |
|________________________________________________________
| p-Xylene    |       0.000   |           0       |           0       |

---- Isosteric heats ----
| Sorbate     | Average     | Minimum     | Maximum     |
|             | kcal/mol    | kcal/mol    | kcal/mol    |
|________________________________________________________

---- Overall system energies ----
Average total energy            :  0.000 kcal/mol

5. 计算结果中的Etotal和Loading图也在附件里。

希望能够得到牛人指导!在此感激不尽!

Set up.png



Energy.png



Constraints Properties Job-Controls.png



Loading.png



Etotal.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
chaizhm: 金币+1, 谢谢~ 2013-01-15 10:41:47
ben_ladeng: 专家考核 2013-04-01 18:45:29
你的xylene用的是全原子模型还是粗粒化模型?一个CH2用一个小球表示的是粗粒化模型,就像一个CH4可以用一个小球表示。
好好学习,天天向上。
2楼2013-01-11 13:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿龙6678

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by zyj8119 at 2013-01-10 17:55:14
你的xylene用的是全原子模型还是粗粒化模型?一个CH2用一个小球表示的是粗粒化模型,就像一个CH4可以用一个小球表示。

用的是全原子模型啊!我还试过同学帮我用DFT算出来优化过结构的Xylene,也是全原子模型,结果也是loading=0.
3楼2013-01-11 22:01:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
阿龙6678: 金币+5 2013-01-13 06:20:02
引用回帖:
3楼: Originally posted by 阿龙6678 at 2013-01-11 22:01:48
用的是全原子模型啊!我还试过同学帮我用DFT算出来优化过结构的Xylene,也是全原子模型,结果也是loading=0....

那么你不用全原子模型看看,用细致原子模型看看。。。
好好学习,天天向上。
4楼2013-01-12 15:59:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿龙6678

银虫 (初入文坛)

非常感谢你的回复!! 细致原子模型是怎么得出来的呢?
5楼2013-01-13 06:19:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 阿龙6678 at 2013-01-13 06:19:48
非常感谢你的回复!! 细致原子模型是怎么得出来的呢?

就是CH2的每个原子都画出来。。。。
好好学习,天天向上。
6楼2013-04-06 20:33:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一二三吗

新虫 (小有名气)

请教一下楼主  你的sorption loading 图中的那个per cell 单位是如何转换成别的单位的呢 求指导!!
7楼2013-09-26 14:59:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ligesi

银虫 (小有名气)

引用回帖:
5楼: Originally posted by 阿龙6678 at 2013-01-13 06:19:48
非常感谢你的回复!! 细致原子模型是怎么得出来的呢?

你的C-C键是参考文献是值得吗/
新手驾到,O(∩_∩)O哈哈~
8楼2013-11-21 21:50:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feona911

银虫 (小有名气)

楼主这些你是在哪里学的啊
9楼2017-06-29 15:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿龙6678 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见