应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 【求助】CuBTC MOF + CO2 MuSiC 问题
1025/11返回列表首页上一页345678下一页末页
查看: 6770  |  回复: 101
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 3 个 模拟EPI ，点击这里进行查看

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-11-05 21:18:08
引用回帖:
Originally posted by qphll at 2010-11-05 09:40:53:


这个不是问题, 我测试的压力都不算高. 事实上我注意到这个问题,
从文献读出的是压力值, 在该温度下, 压力和逸度的转换(RP方程) 如下, 单位均为KPa

压力         逸度
11.55       11.54
41.28       41 ...

我记得是不论你在哪个地方输入，都应该是逸度值，GCMC的输入肯定是逸度，如果程序能自动进行转换，那么势必需要物性参数，不然随便你输入什么吸附质，程序都能转换，从道理上说不通。Music这个程序不可能自带这么多种物质的物性参数，如果这样，想必它就不会免费了。
赞一下(2人) 回复此楼
41楼2010-11-05 09:47:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+2):谢谢 2010-11-05 22:34:53
还发现，atm-atm文件里关于CO2原子和CuBTC之间的库伦作用，你设置为了WFCOUL计算，和你使用的EMAP中的EWALD方法是冲突，不知道这样情况下，程序运行时会用哪个方法。
赞一下(2人) 回复此楼
42楼2010-11-05 19:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+10):谢谢 2010-11-07 11:38:16
引用回帖:
Originally posted by ghcacj at 2010-11-05 09:47:23:

我记得是不论你在哪个地方输入，都应该是逸度值，GCMC的输入肯定是逸度，如果程序能自动进行转换，那么势必需要物性参数，不然随便你输入什么吸附质，程序都能转换，从道理上说不通。Music这个程序不可能自带这 ...

逸度, 逸度, 逸度....

这是相关部分的源代码:

    !** array for fugacities
    Allocate(gcmcparams%fuglist(nsims, nspc), STAT=error)
    If (error/=0) Call allocErrDisplay(__FILE__,__LINE__,'fuglist')

    !** Get the fugacity
    Do i=1, nsims

      !** Set the pressure of each component
      Do j=1, gcmcparams%nspc
        spc = gcmcparams%gcmcspc(j)%spc

        !** pp is in kPa
        pp = gcmcparams%gcmcspc(j)%fuglist(i)%pressure
!        Write(*,*) spc, pp

        Call eos_setpressure(gcmcparams%eosparams, spc, pp)
      End Do

      !** Get the fugacity of each component
      Do j=1, gcmcparams%nspc
        spc = gcmcparams%gcmcspc(j)%spc
        fug = eos_getfugacity(gcmcparams%eosparams, spc) ! fug [=] kPa
        gcmcparams%gcmcspc(j)%fuglist(i)%fugacity = fug

        !** Get the excess chemical potential (B in Adams Notation)
        sivolume = volume*1.0e-30*Nav ! convert to m^3/mole
        B = Log(fug*1.0e3*sivolume/(Rgas*tk))

       !** This is Z/Omega, got from ideal parameters
      ! wrong:      Zig = eos_getConfInteg(gcmcparams%eosparams, j)
      ! we should be passing the sorbtype , not the index of gcmcmove
       Zig = eos_getConfInteg(gcmcparams%eosparams, spc)
       ratio = log(Zig)

       !** we need (PV)/(RTZ)
       B = B - ratio
!        Write(0,'(2a,i4,a,f16.2,6f16.10)') __FILE__,": ",__LINE__, &
!        " B ", ratio,B,Zig

        gcmcparams%gcmcspc(j)%fuglist(i)%B = B

        !** Get mu the chemical potential
        !** Calculate the DeBroglie wavelength
        mass = molecules_getmass(spc)
        mass = mass*1.0e-3              ! convert to kg
        Lambda = hplanck/Sqrt(twopi*mass/Nav*Rgas/Nav*tk)
        murti = (B - Log(sivolume/Nav/Lambda**3))
        gcmcparams%gcmcspc(j)%fuglist(i)%murti = murti
!        Write(*,*) fug, B, ratio, murti
      End Do
    End Do

由此可见, 不管是准备额外的压力文件, 还是直接在Control文件中写入数值, MuSiC都需要逸度.

我的错觉来源于这个源文件 virial.F90, 但是这个, 其实只是混合体系的再分配而已.

所以, 逸度, 逸度, 逸度.
赞一下(1人) 回复此楼
Life, Love, Laugh.
43楼2010-11-07 05:02:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+8):谢谢 2010-11-07 11:38:06
引用回帖:
Originally posted by ghcacj at 2010-11-05 19:56:27:
还发现，atm-atm文件里关于CO2原子和CuBTC之间的库伦作用，你设置为了WFCOUL计算，和你使用的EMAP中的EWALD方法是冲突，不知道这样情况下，程序运行时会用哪个方法。

这个没有问题. 原因有二:

(1) 目前MuSiC只是支持生成 LJ Pamp和 EWALD Emap, 参加官方文档.

虽然这个文档是零八年三月更新的, 但是在最新的版本源代码中, 我没有找到针对这个的更新. 所以应该还是只是支持这两类? (欢迎讨论!)

http://zeolites.cqe.northwestern ... ntation/node39.html

Current code makes pmaps and emaps for only two types of interactions. pmap for LJ, emap for EWALD. There are no combined maps yet. Non-orthogonal cases are still being tested (as of 21 march 2008)

(2) 在生成emap文件以后, 在实际计算时, 程序只是考虑emap文件中的数值. atm_atm_file中关于库仑计算的WFCOUL参数, 已经影响不到实际计算了.

总结一下,
在我这里, 如果有用到EMAP, 那么EMAP的生成是EWALD方法(只能选择这个方法), GCMC计算时对库仑力的处理, 从EMAP来, 也是EWALD方法; 如果不用EMAP, 我对体系库仑力的处理采用WFCOUL.
赞一下(1人) 回复此楼
Life, Love, Laugh.
44楼2010-11-07 05:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by qphll at 2010-11-07 05:14:16:



这个没有问题. 原因有二:

(1) 目前MuSiC只是支持生成 LJ Pamp和 EWALD Emap, 参加官方文档.

虽然这个文档是零八年三月更新的, 但是在最新的版本源代码中, 我没有找到针对这个的更新. 所以应该还是只是 ...

目前对于gcmc来说，想使用ewald方法，就必须用emap，而wfcoul方法只支持on the fly的计算，所以目前我都采用吸附质-吸附质之间的静电力用wfcoul，吸附剂-吸附质的静电力用ewald方法。（也是被逼的）
你的学习能力和态度都很强，佩服佩服。
赞一下(1人) 回复此楼
45楼2010-11-07 11:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
能将你生成EMAP时用到的Sorb-Sorb（分子-分子文件）贴出来吗？
赞一下(1人) 回复此楼
46楼2010-11-08 13:36:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
引用回帖:
Originally posted by ghcacj at 2010-11-08 13:36:12:
能将你生成EMAP时用到的Sorb-Sorb（分子-分子文件）贴出来吗？

CuBTC   CuBTC   NCOUL   OFF
CuBTC   CuBTC    COUL   OFF

probe     probe    NCOUL OFF
probe     probe    COUL  OFF


probe     CuBTC   NCOUL OFF  
probe     CuBTC    COUL  SUM FAST FIXED EWALD SFACTOR  KMAX@15  KAPPA@6.7  LOCUT@1e-10
赞一下 回复此楼
Life, Love, Laugh.
47楼2010-11-11 11:38:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
这两天在开会, 几天没有更新, 待会写一些update.
赞一下 回复此楼
Life, Love, Laugh.
48楼2010-11-11 11:43:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
那很奇怪的是我也照着文献做的，计算文件几乎和你一样，为何结果偏小？
赞一下(1人) 回复此楼
49楼2010-11-11 12:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)
引用回帖:
Originally posted by qphll at 2010-11-04 11:02:01:


说得有道理. VDW cutoff的调整, 也是我其中的一个尝试. 回头再来update, 先做会正事....

积累了好多东西没有做, 开会就像休假, 早出晚归...

先update一下这个有关VDW的问题.

上图:



由此可以看到, VDW CUTOFF是蛮能影响结果的.

正像前楼说的,  如果要和文献比较, 那么就要完全一致. 如果是自己来摸索, 一定也要做一些测试, 直到CUTOFF几乎不影响最终的结果为止.

[ Last edited by qphll on 2010-11-11 at 12:37 ]
赞一下 回复此楼
Life, Love, Laugh.
50楼2010-11-11 12:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qphll 的主题更新
1025/11返回列表首页上一页345678下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 化学工程321分求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-15 18:49 by a不易
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 泣血叩求调剂恩，愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
关闭
请填处理意见
关闭