| 查看: 6845 | 回复: 101 | ||||||||
| 本帖产生 3 个 模拟EPI ，点击这里进行查看 | ||||||||
qphll金虫 (正式写手)
|
[交流]
【求助】CuBTC MOF + CO2 MuSiC 问题 已有12人参与
|
|||||||
|
在用MuSiC测试体系: CO2+CuBTC的时候, 发现了两个问题: (1) 采用与文献相同模拟, 相同参数, 相同软件的情况下, 吸附量差别很大; (2) 改变压力做吸附等温线, 不同压力下的吸附量竟然差别不大. 应该是哪里的设置有些问题. 让我细细道来, 大家评论一下. 拿来测试的文献是: Molecular Simulations and Experimental Studies of CO2, CO and N2 adsorption in Metal-Organic Frameworks, J. Phys. Chem. C, 114, 15735, 2010 更具体一点, 试图重复的数据是该文献Figure-2中, 298K下CO2在CuBTC中的吸附. 在GetData的帮助下, Figure-2中CO2在CuBTC的吸附量大致如下: Pressure, KPa Loading, mol/kg 41.28 2.22 120.67 5.40 321.05 10.93 530.09 12.94 720.96 14.03 1149.34 14.95 以41.28KPa这个压力点为例, 贴出所有我计算中用到的文件. |
回复此楼» 收录本帖的淘帖专辑推荐
 有机场效应晶体管及有机太阳能电池 | MOF相关 | 第一性原理和电化学 | Software Operation |
| mof 材料 |
» 猜你喜欢
材料专硕322
已经有7人回复
-
求调剂
已经有20人回复
-
一志愿郑州大学085600求调剂
已经有20人回复
-
化学调剂
已经有17人回复
-
3-溴吡啶-4-甲醛合成
已经有9人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有3人回复
-
285求调剂
已经有9人回复
-
材料求调剂
已经有11人回复
-
调剂
已经有8人回复
-
材料调剂
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 截断半径设置问题
已经有7人回复
- 关于工业上对CO2的捕集方法的疑问
已经有8人回复
- 请问为什么CO2-TPD结果是这样？
已经有12人回复
- 如何建立MOF模型？CO2、N2等气体模型如何建立？
已经有6人回复
- MuSic post处理结果出错
已经有10人回复
- Co的计算问题！新手求助
已经有10人回复
- HTML源代码实现不了想要的效果，求助
已经有7人回复
- Excel 迭代问题
已经有4人回复
- 分子中插入CO的方法
已经有9人回复
- m文件调用问题
已经有4人回复
- 求助Mapmaker使用过程自动关闭
已经有3人回复
- 不同测试温度下的晶胞参数差异问题
已经有5人回复
- 吸附数据单位换算！！
已经有5人回复
- 为何用同样的方法，软件之间差异这么大？
已经有16人回复
- 【求助】有关截断能测试问题
已经有6人回复
- 【请教】稀土的MOF的一点小问题
已经有9人回复
- 【求助】MC方法计算吸附量
已经有14人回复
- 【求助】关于MOF的问题
已经有7人回复
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
分别生成的两个pmap文件命名为: C_CO2.M4.LJ111.pmap O_CO2.M4.LJ111.pmap 至于emap, 针对每种吸附质, 只需要生成一个, 其原理就是生成该吸附质的库仑力格点. (a) Control File 内容 ------ General Information ------------------------------------------ CuBTC_M4, to generate eamp@CuBTC 1 1 1 1 2 191283 4 CuBTC_M4.emap.res CuBTC_M4.emap.con ------ Atomic Types -------------------------------------------------- 7 Patom Patom.atm Ca_CuBTC_M4 Ca_CuBTC_M4.atm Cb_CuBTC_M4 Cb_CuBTC_M4.atm Cc_CuBTC_M4 Cc_CuBTC_M4.atm Cu_CuBTC_M4 Cu_CuBTC_M4.atm O_CuBTC_M4 O_CuBTC_M4.atm H_CuBTC_M4 H_CuBTC_M4.atm ------ Molecule Types ------------------------------------------------- 2 probe probe.mol CuBTC_M4 CuBTC_M4.mol ------ Simulation Cell Information ------------------------------------ CuBTC_M4 1, 1, 1 1, 1, 1 ------ Forcefield Information ------------------------------------------- BASIC ATM atm_atm_file mol_mol_file intramolecular_file ------ Mapmaker Information ----------------------------------------------- 1 # Number of maps to make # PLEASE KEEP THIS 'BLANK LINE'!! CuBTC_M4 # Name of the sorbent probe # Name of prob atoms of the sorbate COUL EWALD # Typer of interaction to be mapped 0.1 # Grid of spacing, Angstrom 500.0 # High end potential cutoff, KJ/mol AUTO # Map filename or AUTO ------ Configuration Initialization ------------------------------------- probe GCMC NULL CuBTC_M4 FIXED NULL (b) probe.mol 内容: Molecule_Name: probe CHARGED Coord_Info: Listed Cartesian Rigid 1 1 0.000000 0.000000 0.000000 Patom 1.0 0 0 Molecule_DOF: 3 (c) Patom.atm 内容: ##### Basic Atom Information Atom_Name: Patom Atom_Symbol: P Atom_SS_Charge: 0.0 Atom_SZ_Charge: 0.0 Atom_Mass: 14.0067 Atom_Valency: 0 此处生成的CuBTC的emap命名为: CuBTC_M4.111.emap 至此, 前期的准备工作, 包括: atm文件, mol文件, pmap文件, emap文件都准备完毕. [ Last edited by qphll on 2010-11-3 at 14:50 ] |
4楼2010-11-03 14:15:15
ghcacj
荣誉版主 (著名写手)
阿超
- 模拟EPI: 16
- 应助: 26 (小学生)
- 贵宾: 2.764
- 金币: 2543.5
- 散金: 8899
- 红花: 35
- 沙发: 3
- 帖子: 1376
- 在线: 1199小时
- 虫号: 513281
- 注册: 2008-02-27
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3): 2010-11-03 15:28:55
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3): 2010-11-03 15:28:55
|
看了你的atm-atm文件，里面的截断半径不对，文献里是12.8，你的是13。Locut最好选择0.1。 此外，制作map的时候，计算步数不能只有1步，太少了，应该多迭代几次，最少100W步吧。CO2制作emap的时候，应该还是将C和O原子的emap分开制作。 吸附量算的这么大，我怀疑是emap做的不对，建议你将Emap分开制作，然后再试试。 [ Last edited by ghcacj on 2010-11-3 at 15:10 ] |
13楼2010-11-03 14:35:37
ghcacj
荣誉版主 (著名写手)
阿超
- 模拟EPI: 16
- 应助: 26 (小学生)
- 贵宾: 2.764
- 金币: 2543.5
- 散金: 8899
- 红花: 35
- 沙发: 3
- 帖子: 1376
- 在线: 1199小时
- 虫号: 513281
- 注册: 2008-02-27
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
15楼2010-11-03 15:15:59
ghcacj
荣誉版主 (著名写手)
阿超
- 模拟EPI: 16
- 应助: 26 (小学生)
- 贵宾: 2.764
- 金币: 2543.5
- 散金: 8899
- 红花: 35
- 沙发: 3
- 帖子: 1376
- 在线: 1199小时
- 虫号: 513281
- 注册: 2008-02-27
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
17楼2010-11-03 15:23:50
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zh1987hs(金币+4):鼓励 2010-11-04 10:49:21
ghcacj(金币+10, 模拟EPI+1):鼓励，支持 2010-11-04 10:53:28
zh1987hs(金币+4):鼓励 2010-11-04 10:49:21
ghcacj(金币+10, 模拟EPI+1):鼓励，支持 2010-11-04 10:53:28
|
update 2: 重新生成了CO2的两个pmap文件, 和CuBTC的EMAP文件. 改动如下: (1) Pmap的control file中 cutoff阀值由原先500KJ/mol 修改为200KJ/mol (2) EMAP的control file中 cutoff阀值由原先500KJ/mol 修改为100KJ/mol (3) atm_atm_file中, CuBTC原子的vdw HICUT为13.1A, LOCUT为0.01A; CO2原子和CO2-CBTC的vdw HICUT为12.8A, LOCUT为0.01A; (4) PMAP的迭代次数由原先的1修改为1000000; EMAP的生成方法和原先一样, 还是用一个带单位电荷的探针, 迭代步数依然为1. 在map文件生成的情况下, 单点计算前面提到的那留个压力, 同时在不适用map文件的情况下计算41.28Kpa下的吸附 试验结果如下: (1) 图一张  尽管我和文献一样, 每点迭代2E7, 取后50%平均, 但是文献是用了3*3*3, 我只是用1*1*1, 这样看起来结果吻合得还算可以. (2) 同样的压力(41.28KPa), 不用map, 计算结果为2.24kg/mol, 使用PMAP+EMAP后的结果是 2.34kg/mol. 文献报导的改点吸附为 2.22kg/mol. 貌似不用map文件, 计算更加可靠一些? (假设文献报导的值是standard reference) 言下之意, PMAP或者EMAP还有需要改进的地方??? HOW???? (3) MAP文件对计算速度的提升还是很明显的, 依然是(2)中的比较, 不用map文件, 耗时7.1759 hrs, 使用了PMAP+EMAP, 耗时 58.714 min. 注意的是, 不用map文件时, 迭代步数为 1E7; 用Map文件下的计算, 都是 2E7 这样折算下来, 速度差别, 那是相当相当的大........ |
» 本帖已获得的红花（最新10朵）
rdfce628楼2010-11-04 10:40:02
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
35楼2010-11-05 00:13:28
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+10):谢谢 2010-11-11 14:11:03
ghcacj(金币+10):谢谢 2010-11-11 14:11:03
|
我应该还有一个update没有做, 就是有关能量的波动曲线. 今天是没有时间做那个事情了. 但是现在如果回头看看做的这些的尝试, 有多少事实上是不需要的呢, 或者有哪些是在计算之前就需要弄清楚的呢? (1) 体系的大小 (2) VDW参数 (3) CUTOFF数值, VDW, COUL (WFCOUL or EWALD ?) (4) 计算步数, 能量曲线, 接受概率 (5) MC MOVE的步长 (6) 文献 (实验, 计算) 注意: (a) 如果要修改CuBTC 为 2*2*2, 很有可能会遇到报错, 说有too many spec-spec interaction. 需要修改ssbasic.F90 Line 1532, 将那里的体系限制从10000改大. (b) 在 CuBTC 为 2*2*2的情况下, map文件的生成是一个问题, 不管是emap, 还是pamp, 除了需要修改一些参数, 例如, ewald.F90中是hard coded to 50000. 你需要增大这个阀值. 但是, 增大了也没有用. 因为你的内存很有可能是不够的, 我大概估算了一下, 在 0.1 angstrom 为步长的情况下, map文件的生成, 似乎需要大于5G的内存分配. 或者我疏漏了什么, 但是有谁做过这个尝试嘛??? 不管怎样, 将这些问题放在心里, 那么你就可以在重现文献数据, 练手熟练的情况下, 放手尝试新的体系了, 而且会对你自己的计算结果有信心! (欠一个能量图的分析......) |
60楼2010-11-11 13:58:15
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+15, 模拟EPI+1):very very good 2010-11-14 08:43:48
ghcacj(金币+15, 模拟EPI+1):very very good 2010-11-14 08:43:48
|
针对这个问题, 做一些update. 先来看一下emap的mol_mol_file: Lattice Lattice NCOUL OFF Lattice Lattice COUL OFF Probe Lattice NCOUL OFF Probe Lattice COUL SUM FAST FIXED EWALD SFACTOR KMAX@15 KAPPA@6.7 LOCUT@1e-10 Probe Probe NCOUL OFF Probe Probe COUL OFF 首先需要指出的是, 这里除了 Probe 与 Lattice的COUL作用之外, 其他的相互作用全部关闭. 来仔细看看这句: Probe Lattice COUL SUM FAST FIXED EWALD SFACTOR KMAX@15 KAPPA@6.7 LOCUT@1e10 (1) “Probe” 和 “Lattice”, 这个就像“路人甲”一样, 你想怎么称呼都可以; (2) “COUL”是从 ff.F90过来的, 在那个源程序里面定义了体系的能量从COUL和NONCOUL计算得到; (3) “SUM”, “FAST”, “FIXED”,“EWALD”来自于ssdrive.F90, 我们来一一解读一下. (a) SUM Line 94, 定于了 “SUM”: “Type(SSSumParams), Pointer :: sum”. 很容易理解, 在任意位置上, Probe与Lattice的COUL作用, 需要累加起来. (b) FAST Line 495,是注释行, 关于“FAST”是这样说的, “fast -- True (False) => evaluate "Fast" (Slow) interactions” 字面上的理解就是, 用了“FAST”那么计算的时候, 关于相互作用的计算就能“加速”,至于如何加速? 不知道, 也不需要知道. (c) FIXED 没有从源程序找到很好的对应, 但是这个是和Control文件中对于Lattice的处理对应的, 我们将Lattice文件Fixed处理, 那么这里也是需要告诉mapmake文件这个信息. 除此之外, 没啥用处. (d) EWALD 这个指定了如果生成EMAP, 前面有提到过, GCMC的EMAP只能通过EWALD得到, 而on-the-fly COUL计算(i.e. 不用EMAP, 在GCMC中直接计算COUL)只能WFCOUL而不能使用EWALD. Line542再次证明了这一点认识: “on the fly ewald sums are implemented for md only” 也就是说, MD时, on the fly COUL计算, 是可以采用EWALD SUM的, 但是MC(包括GCMC)不行.... (4) SFACTOR, KMAX, KAPPA, LOCUT是EWALD方法中的参数. 具体的原理,参见Frenkel&Smit 的第十二章, 或者MuSiC的相关文档. 你也可以从ewald.F90 的subroutine (Line 152 to Line 269) ewald_init看到程序中是怎样指定这些参数的. 另外, sssum.F90也很值得一看, 如果你想了解更多EWALD是怎样在MuSiC中实现的话. 下面一一分析一下这四个参数. (a) SFACTOR Line 205 to Line 207: Case ('sfactor') !** We want to use the structure factor. eparams%useSF = .TRUE. 它是一个逻辑变量, 默认值是.FALSE., 所以在这里需要加入这个关键词来激活它, 变为.TRUE. 激活这个逻辑变量以后, 程序做了一些什么呢? 参加 Line 284 那边的注释. !** Call the structure factor routine. This calculates the k space !** vectors we will use later when computing the Ewald summation. If !** the coordinates of the ions are fixed, this needs to be calculated !** only once. If they change, it must be recalculated. (b) KMAX Line 184: “Number of k "boxes" in each direction (x,y,z)”. 你可能注意到在程序里面, KAMX is hard coded to be 5000, 你可能回想增大这个数值, 来处理更大的体系. 但是, 不幸的是, 你的内存很有可能比KMAX参数更快地达到阀值, 所以, 修改这个参数在实际使用上, 没有意义, 我没有测试过, 但是我觉得这里 KMAX@15, 其实15这个数值, 是没什么意义的, 只是需要一个input而已, 你用15, 10还是8, 没有区别. (c) KAPPA 是和 EWALD SUM里面的Gaussian function有关的一个参数, 它定义了 Gaussian Function的带宽. 如果你看Frenkel&Smit 的第十二章, 这里程序中的KAPPA参数值等于Frenkel&Smit Chap.12中的alpha^0.5. (d) LOCUT 这个很好理解, 请看ewald.F90文件中的这一行: Line 200, “ !** Low end cutoff for the reciprocal space k vectors.” 接下来的问题就是, 我们应该如何选择合适的KAPPA和KMAX参数? 这两个参数的选择是基于EWALD SUM 计算时对于实空间和倒空间的“利用”和平衡. 一般而言, KMAX必须让SUMMATION在倒空间的计算满足随着KAPPA参数收敛, 同时尽量让实空间的计算也能快速收敛. 比方说, 减小KAPPA参数可以使得倒空间的收敛加快, 但是此时实空间的收敛就会减慢. 在非常小的KAPPA参数下, 实空间的EWALD SUM 和WFCOUL方法一样. 在其他MAP MAKE 参数相同的情况下, 测试不同的KAPPA数值, 看生成MAP需要的不同时间, 来找到最优值, 但是一般来讲, 这个对于MC(包括GCMC, 和MC里面相应的EMAP的生成)没有实际意义. 但是对于MD来讲, 由于on the fly COUL计算可以采用EWALD SUM 方法, 这样的测试还是能有明显提速功效的. 好了, 对于生成EMAP时EWALD 方法的介绍差不多了解这些就足够了. 还有值得提醒的就是, MuSiC中支持的EWALD 方法只是对正交盒子有效 (orthorhombic unit cell). 这个对于GCMC问题不大, zeolite, MOF等材料, 大多数还是orthombic的, 但是如果是自己生成的初始结构, 想要算吸附, 然后又想要使用EMAP的时候, 就需要注意这一点! |
66楼2010-11-13 22:11:52
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ghcacj(金币+10):精彩 2010-11-23 09:19:17
ghcacj(金币+10):精彩 2010-11-23 09:19:17
|
小update一下. 常常看到文献上有这样的说法: "At least 20 million trials were used in the simulations, among which the first half was used for equilibrium, and the last half was used to calculate the ensemble averages." 那么这里的"20 million trials" 和 "50%" 这两个关键词是怎样来的? 同样用我们这里的例子, 同样用类似前面某次update提到的体系能量变化图, 我们再仔细看看: 上图!  大家可以看到, 我这里也是跑了 20 million, 也是将前50%扔掉, 后50%系综平均. 但是很明显, 从大图(黑线)可以看出, 体系的总能量在(10E6~14E6)这个阶段还是有很明显的降低. 从小图1(红线, 18E6~20E6)看出, 就算是将总迭代(2E7)的80%扔掉, 体系的能量在最后的20%那里还是在降低. 从小图2(蓝线, 19E6~20E6) 看起来, 似乎体系还是在波动, 但是考虑到波动的绝对值小于 20KJ (注意, 这里是 20KJ for the whole system, NOT 20KJ/mol), 再结合大图(黑线)最后 19E6~20E6处体现出来的'平坦', 我们有理由说服自己, 在19E6步以后, 体系已经达到平衡态了. 差不多可以安全地做系综平均了. 但是, 为了百分百放心, 最好在这个基础上, 再往后算10E6步, 然后看看这个区间 19E6~30E6 的能量波动, 如果能量的波动像小图2(蓝线)所示范围震荡, 那么OK, 放心大胆系综平均, 得出你的结论吧. 反之, 如果体系能量还是像小图1(红线)所示, 有一路下降的趋势(哪怕斜率平坦), 那么对不起, 体系还没有平衡, 继续跑吧...... 以上小小update, 可以看作如何选择迭代总步数, 如果选择系综平均(到底扔掉前面的百分之多少)的一些依据. [ Last edited by qphll on 2010-11-23 at 07:02 ] |
76楼2010-11-23 06:56:08
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
(1) 生成CO2的两个pmap文件, 一个是O原子在CuBTC的VDW格点; 另外一个是C原子在CuBTC的VDW格点. 以C原子为例. (a). Control File ------ General Information ------------------------------------------ CO2 in CuBTC_M4, to generate C_CO2 Pamp@CuBTC 1 1 1 1 2 191283 4 C_CO2.CuBTC_M4.pmap.res C_CO2.CuBTC_M4.pmap.con ------ Atomic Types -------------------------------------------------- 7 C_CO2 C_CO2.atm Ca_CuBTC_M4 Ca_CuBTC_M4.atm Cb_CuBTC_M4 Cb_CuBTC_M4.atm Cc_CuBTC_M4 Cc_CuBTC_M4.atm Cu_CuBTC_M4 Cu_CuBTC_M4.atm O_CuBTC_M4 O_CuBTC_M4.atm H_CuBTC_M4 H_CuBTC_M4.atm ------ Molecule Types ------------------------------------------------- 2 CO2_C_prob CO2_C_prob.mol CuBTC_M4 CuBTC_M4.mol ------ Simulation Cell Information ------------------------------------ CuBTC_M4 1, 1, 1 1, 1, 1 ------ Forcefield Information ------------------------------------------- BASIC ATM atm_atm_file mol_mol_file intramolecular_file ------ Mapmaker Information ----------------------------------------------- 1 # Number of maps to make # PLEASE KEEP THIS 'BLANK LINE'!! CuBTC_M4 # Name of the sorbent CO2_C_prob # Name of prob atoms of the sorbate NCOUL LJ # Typer of interaction to be mapped 0.1 # Grid of spacing, Angstrom 500.0 # High end potential cutoff, KJ/mol C_CO2.LJ.M4.pmap # Map filename or AUTO ------ Configuration Initialization ------------------------------------- CO2_C_prob GCMC NULL CuBTC_M4 FIXED NULL (b) mol_mol_file CuBTC_M4 CuBTC_M4 NCOUL OFF CuBTC_M4 CuBTC_M4 COUL OFF CO2_C_prob CO2_C_prob COUL OFF CO2_C_prob CO2_C_prob NCOUL OFF CO2_C_prob CuBTC_M4 COUL OFF CO2_C_prob CuBTC_M4 NCOUL BASIC LJ FAST (c) intramolecular_file Intra: CO2_C_prob Intra: CuBTC_M4 (d) atm_atm_file ### CuBTC_M4, JPCB, 2006, 110, 17776 Ca_CuBTC_M4 Ca_CuBTC_M4 LJ SIG@3.75 EPS@44.91 HICUT@13.0 LOCUT@0.5 Cb_CuBTC_M4 Cb_CuBTC_M4 LJ SIG@3.55 EPS@35.23 HICUT@13.0 LOCUT@0.5 Cc_CuBTC_M4 Cc_CuBTC_M4 LJ SIG@3.55 EPS@35.23 HICUT@13.0 LOCUT@0.5 Cu_CuBTC_M4 Cu_CuBTC_M4 LJ SIG@3.11 EPS@2.52 HICUT@13.0 LOCUT@0.5 O_CuBTC_M4 O_CuBTC_M4 LJ SIG@2.96 EPS@73.98 HICUT@13.0 LOCUT@0.5 H_CuBTC_M4 H_CuBTC_M4 LJ SIG@2.42 EPS@15.10 HICUT@13.0 LOCUT@0.5 # L-B Mixture for Intra-CuBTC atoms Ca_CuBTC_M4 Cb_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Ca_CuBTC_M4 Cc_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Ca_CuBTC_M4 Cu_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Ca_CuBTC_M4 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Ca_CuBTC_M4 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cb_CuBTC_M4 Cc_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cb_CuBTC_M4 Cu_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cb_CuBTC_M4 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cb_CuBTC_M4 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cc_CuBTC_M4 Cu_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cc_CuBTC_M4 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cc_CuBTC_M4 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cu_CuBTC_M4 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 Cu_CuBTC_M4 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CuBTC_M4 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 ########################################################################## ########################################################################## ### CO2 C_CO2 C_CO2 LJ SIG@2.8 EPS@27.0 HICUT@13.0 LOCUT@0.5 O_CO2 O_CO2 LJ SIG@3.05 EPS@79.0 HICUT@13.0 LOCUT@0.5 # L-B Mixture for Intra-NH3 atoms C_CO2 O_CO2 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 ## CO2-CuBTC_M4, L-B Mixture C_CO2 Ca_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 C_CO2 Cb_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 C_CO2 Cc_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 C_CO2 Cu_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 C_CO2 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 C_CO2 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 Ca_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 Cb_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 Cc_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 Cu_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 O_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 O_CO2 H_CuBTC_M4 LJ SIG@LBMIX EPS@LBMIX HICUT@13.0 LOCUT@0.5 ########################################################################## ############################################################################## ##### Coulombic ############################################################################ ### CO2-CO2 ##### C_CO2 C_CO2 WFCOUL HICUT@13.0 ALPHA@0.2 O_CO2 O_CO2 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 O_CO2 WFCOUL HICUT@13.0 ALPHA@0.1 ### CO2-CuBTC_M4 ##### C_CO2 Ca_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 Cb_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 Cc_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 Cu_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 O_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 C_CO2 H_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 Ca_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 Cb_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 Cc_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 Cu_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 O_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 O_CO2 H_CuBTC_M4 WFCOUL HICUT@13.0 ALPHA@0.1 ################################################################################ [ Last edited by qphll on 2010-11-3 at 14:31 ] |
2楼2010-11-03 14:14:15
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
一些molecule文件和atom文件 (a) CuBTC_M4.mol文件, 不管是生成pmap, 还是emap, 还是GCMC的时候, 都是同一个, 其文件内容如下: 字数限制, 该文件内容不完整. ## Basic Molecule Information ## CuBTC Model-4 ## CuBTC, Total number of atom in CuBTC UnitCell: 624 ## The Charge Information is from, JPCB, 2006, 110, 17776 ## Cu: 1.098 * 48 ## O: -0.665 * 192 ## Ca: 0.778 * 96 ## Cb: -0.092 * 96 ## Cc: -0.014 * 96 ## H: 0.109 * 96 ## The initial coordinates (dehydrated form) are from "Science, 1999, 283, 1148" Molecule_Name: CuBTC_M4 CHARGED Coord_Info: Listed Cartesian Rigid NOTRANSFORM 624 1 8.341 6.405 -1.376 O_CuBTC_M4 -0.665 0 0 2 18.002 19.938 -1.376 O_CuBTC_M4 -0.665 0 0 3 18.002 6.405 1.376 O_CuBTC_M4 -0.665 0 0 4 8.341 19.938 1.376 O_CuBTC_M4 -0.665 0 0 5 -1.376 8.341 6.405 O_CuBTC_M4 -0.665 0 0 6 -1.376 18.002 19.938 O_CuBTC_M4 -0.665 0 0 7 1.376 18.002 6.405 O_CuBTC_M4 -0.665 0 0 ..... ..... 621 10.195 21.654 4.689 Cb_CuBTC_M4 -0.092 0 0 622 11.362 20.990 5.353 Ca_CuBTC_M4 0.778 0 0 623 10.195 4.689 21.654 Cb_CuBTC_M4 -0.092 0 0 624 11.362 5.353 20.990 Ca_CuBTC_M4 0.778 0 0 Fundcell_Info: Listed 26.34300 26.34300 26.34300 # unit cell edge lengths 90.00000 90.00000 90.00000 # unit cell angles 0.00000 0.00000 0.00000 # origin of unit cell 26.34300 26.34300 26.34300 # effective bound box size (b) CO2_C_prob.mol 文件 Molecule_Name: CO2_C_prob Coord_Info: Listed Cartesian Rigid 1 # Number of atoms in molecule 1 0.000000 0.000000 0.000000 C_CO2 0.0 0 0 Atom文件比较简单, 这里给出分别来自CuBTC和CO2的原子文件为例 (c) Cu_CuBTC.atm ##### Basic atom information Atom_Name: Cu_CuBTC Atom_Symbol: Cu Atom_SS_Charge: 0.0 Atom_SZ_Charge: 0.0 Atom_Mass: 63.546 Atom_Valency: 2 (d) C_CO2.atm ##### Basic Atom Information Atom_Name: C_CO2 Atom_Symbol: C Atom_SS_Charge: 0.4 Atom_SZ_Charge: 0.4 Atom_Mass: 12.0 Atom_Valency: 4 [ Last edited by qphll on 2010-11-3 at 14:42 ] |
3楼2010-11-03 14:14:48
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
GCMC的一些文件 (a) Contron File ------ General Information ------------------------------------------ CO2 in CuBTC_M4 50000000 1000000 1000000 500000 1 91283 3 CO2.CuBTC_M4.res CO2.CuBTC_M4.con ------ Atomic Types -------------------------------------------------- 8 Ca_CuBTC_M4 Ca_CuBTC_M4.atm Cb_CuBTC_M4 Cb_CuBTC_M4.atm Cc_CuBTC_M4 Cc_CuBTC_M4.atm Cu_CuBTC_M4 Cu_CuBTC_M4.atm O_CuBTC_M4 O_CuBTC_M4.atm H_CuBTC_M4 H_CuBTC_M4.atm C_CO2 C_CO2.atm O_CO2 O_CO2.atm ------ Molecule Types ------------------------------------------------- 2 CO2 CO2.mol CuBTC_M4 CuBTC_M4.mol ------ Simulation Cell Information ------------------------------------ CuBTC_M4 1, 1, 1 1, 1, 1 ------ Forcefield Information ------------------------------------------- BASIC MOL atm_atm_file mol_mol_file intramolecular_file ------ Ideal Parameters ----------------------------------------------- Ideal 1 CO2 ------ GCMC Information ----------------------------------------------- 1 298.0 Ideal Parameters 1 5000 1 ------------------------- CO2 41.1949 Null 4 1.0, 1.0, 1.0, 1.0 RINSERT RDELETE RTRANSLATE 0.2, 0 RROTATE 0.2, 0 ------ Configuration Initialization ------------------------------------- CO2 GCMC NULL CuBTC_M4 FIXED NULL -------- Main Datafile Information ------------------------------------- Energy, position, Velocity, pair_energy (b) mol_mol_file CuBTC_M4 CuBTC_M4 NCOUL OFF CuBTC_M4 CuBTC_M4 COUL OFF CO2 CO2 NCOUL BASIC LJ FAST CO2 CO2 COUL BASIC FAST WFCOUL CO2 CuBTC_M4 NCOUL MAP@CuBTC_M4 FAST C_CO2@PMAP@C_CO2.M4.LJ111.pmap O_CO2@PMAP@O_CO2.M4.LJ111.pmap CO2 CuBTC_M4 COUL MAP FAST all@EMAP@CuBTC_M4.111.emap (c) atm_atm_file 与前面生成pmap时使用的一致. (d) intramolecular_file Intra: CO2 Intra: CuBTC_M4 计算过程倒也太平, 略过不表. [ Last edited by qphll on 2010-11-3 at 14:55 ] |
5楼2010-11-03 14:15:36
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
|
计算后处理, post (a) POST Control File # this section is required for working with any post code # # ------------------------------------------------------------ ### Required section ###### -- Post Processor Information ------------ GCMC # Type of simulation GCMC, NVTMC , MD .... CO2.CuBTC_M4.con # basename for config files what your .con files are called in your gcmc folder 1,1 # The first and last numbers of the .con files that were created in the gcmc folder new.met # name for new ctrlfile that will be regenerated CO2.post # Base name for output files 50, 0 # Percentages of data to skipped at start and end # This is used if you know roughly how long the simulation takes to reach equilibrium # The sections below are necessary only if you want the corresponding # analysis performed # ---------------- ALL OF THEM ARE OPTIONAL ------------------------ #### This section is reqd for energy averages in your post code output files #### as of now only total enrgies vs sim. step ------ Post : Energy Average Info ----------------------------------- 20 # Number of blocks into which data should be divided for stats #### This section is reqd for Loading averages in your post code outputfiles #### as of now only species loading vs sim. step (for all species) ------ Post : Loading Average Info ----------------------------------- 20 # Number of blocks into which data should be divided for stats (b) POST中生成的isotherm文件, 由于我这里是每个压力分开计算, 其实是单点吸附. #------------OVERALL ISOTHERM------------- # Molecule : CO2 # Pressure loading # kpa(1 st gcmc species) (molec/uc) 41.194900000000 373.0625000 [ Last edited by qphll on 2010-11-3 at 14:59 ] |
6楼2010-11-03 14:15:54
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
7楼2010-11-03 14:16:13
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
8楼2010-11-03 14:16:37
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
9楼2010-11-03 14:16:48
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
10楼2010-11-03 14:17:02
