24小时热门版块排行榜    

查看: 2431  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见