24小时热门版块排行榜    

查看: 2394  |  回复: 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

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
阿龙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的回帖
查看全部 9 个回答

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的回帖

阿龙6678

银虫 (初入文坛)

非常感谢你的回复!! 细致原子模型是怎么得出来的呢?
5楼2013-01-13 06:19:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见