24小时热门版块排行榜    

查看: 2423  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见(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
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +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
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见