24小时热门版块排行榜    

查看: 7249  |  回复: 101
本帖产生 3 个 模拟EPI ,点击这里进行查看

qphll

金虫 (正式写手)

★ ★ ★
ghcacj(金币+3):谢谢 2010-11-23 14:41:16
引用回帖:
Originally posted by ghcacj at 2010-11-23 13:56:39:
很好,LZ上次说模拟Cu-BTC吸附NH3的例子,有文献吗?

没. 光NH3本身的力场就是问题, 它和H2O的力场一样, 力场让人头疼. 值得做得东西还有就是将MOF的力场精细化, 一般采用的UFF太粗略了. 而且rigid处理总是不很恰当, 尤其是要研究分子在MOF里面的动态行为的时候, 比如扩散.

[ Last edited by qphll on 2010-11-23 at 14:25 ]
Life, Love, Laugh.
81楼2010-11-23 14:21:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+4):谢谢 2010-12-05 17:22:29
引用回帖:
Originally posted by qphll at 2010-11-23 14:21:33:


没. 光NH3本身的力场就是问题, 它和H2O的力场一样, 力场让人头疼. 值得做得东西还有就是将MOF的力场精细化, 一般采用的UFF太粗略了. 而且rigid处理总是不很恰当, 尤其是要研究分子在MOF里面的动态行为的时候, ...

NH3我不了解,但是H2O我知道,至今没有一个很确切的答案,反正力场就有很多种,各有各的优势,UFF是粗,但是很多人考虑到计算的结果还可以接受,也就都习惯用了,精细化的力场我看容易被人说移植性不好,反正是各有各的道理。Rigid我觉得也是无奈的选择,都知道Flexible好,尤其是大分子吸附时,可是对于不能并行的MC来说,Flexible的计算量很难想象。
82楼2010-11-23 14:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
问楼主个问题:当考虑库仑相互作用时,ALPHA@0.2  这个参数是什么意思啊?谢谢楼主
83楼2010-12-03 23:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★
zh1987hs(金币+3):谢谢 2010-12-05 17:22:54
update,

前1E7扔掉, 分析了后1E7的能量变化. 看起来, 在1.5E7以后, 体系差不多到达平衡态啦.

Life, Love, Laugh.
84楼2010-12-05 13:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★ ★ ★
zh1987hs(金币+4):谢谢 2010-12-05 17:23:05
引用回帖:
Originally posted by gjh123 at 2010-12-03 23:05:45:
问楼主个问题:当考虑库仑相互作用时,ALPHA@0.2  这个参数是什么意思啊?谢谢楼主

回复晚啦, 望见谅.

你的这个问题, 我前面有讲过:

KAPPA 是和 EWALD SUM里面的Gaussian function有关的一个参数, 它定义了 Gaussian Function的带宽.
如果你看Frenkel&Smit 的第十二章, 这里程序中的KAPPA参数值等于Frenkel&Smit Chap.12中的alpha^0.5.

请参见 本帖 69楼 我的某次update, 原则上讲, 不怎么需要担心这个参数. 但是弄明白些总是好的.
Life, Love, Laugh.
85楼2010-12-05 13:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)

谢谢楼主的热心讲解!
86楼2010-12-05 14:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
问LZ一个问题,目前文献中计算吸附等温线,一般分两种,一种是直接用music之类的MC软件跑等温线,也就是说在每个压力点下,吸附剂起始状态都是空的,然后根据uVT计算吸附量。另外一种方法是除了第一个压力点用空吸附剂作为起始状态,其余均用前一个压力点作为起始状态,这样的好处是可以节省计算时间,不过这样做的人往往宣称这样的方法更贴合吸附的实际情况,请问是这样的吗?
87楼2010-12-11 10:45:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★ ★
zh1987hs(金币+2):谢谢 2011-01-06 22:50:06
引用回帖:
Originally posted by ghcacj at 2010-12-11 10:45:22:
问LZ一个问题,目前文献中计算吸附等温线,一般分两种,一种是直接用music之类的MC软件跑等温线,也就是说在每个压力点下,吸附剂起始状态都是空的,然后根据uVT计算吸附量。另外一种方法是除了第一个压力点用空吸 ...

正解!

其实两种方法都是正确的. 只要每种方法的精度都达到要求, 那么应该是给出相同的结果. 其实这也是为什么MC可以续算的原因.

另外, 我推荐第二种方法. 确实能省不少时间.
Life, Love, Laugh.
88楼2010-12-11 11:28:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+2):谢谢 2011-01-06 22:50:14
引用回帖:
Originally posted by qphll at 2010-12-11 11:28:38:


正解!

其实两种方法都是正确的. 只要每种方法的精度都达到要求, 那么应该是给出相同的结果. 其实这也是为什么MC可以续算的原因.

另外, 我推荐第二种方法. 确实能省不少时间.

的确,我尝试过,尤其是在吸附量平衡后,时间大幅降低,查看计算细节,也发现多是平移和旋转被接受,效率的确提高很多。不过这里有一个问题,就是如果用Music来这么做,就必须一个点一个点的跑,要是能通过修改程序,让程序在计算等温线时自动选择前一个点的restart文件,则方便很多。
89楼2010-12-11 14:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
90楼2011-01-02 14:30:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qphll 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 16/800 2026-08-20 19:26 by wuchongjun
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
信息提示
请填处理意见