| 查看: 4686 | 回复: 38 | ||||
pdl9527专家顾问 (小有名气)
|
[求助]
MS的Dmol3模块COSMO计算的问题。 已有1人参与
|
|||
|
我在别的区提过了，没人解答。。。。 http://muchong.com/bbs/viewthread.php?tid=7397241 |
回复此楼» 猜你喜欢
东方理-中科大联合博士生招聘
已经有0人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
吉林大学材料物理本科生求问调剂信息
已经有23人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- DMOL3求吸附能遇到问题
已经有15人回复
- MS discover模块计算问题
已经有6人回复
- MS计算失败显示“找不到指定的模块”求指教！谢谢
已经有5人回复
- 请教，关于Dmol3，如何在优化好的结构中计算能量和频率
已经有4人回复
- MS中sorption模块使用自定义力场出现的问题
已经有26人回复
- MS中sorption模块得到吸附量的单位问题
已经有11人回复
- Dmol3计算周期性体系振动频率时k点的设置问题
已经有3人回复
- MS中的castep模块计算Fe-Mn体系，总是失败。
已经有6人回复
- 请教dmol3和ms关系？如何下载dmol3？
已经有12人回复
- MS-dmol3电荷密度分析
已经有14人回复
- 【求助】Dmol3 晶体计算
已经有15人回复
- 【求助】MS中dmol3模块——两种分子的共吸附优化
已经有9人回复
- 【求助】MS5.0中Dmol3计算频率问题
已经有21人回复
- 【讨论】MS中sorption模块的MC计算
已经有24人回复
- 【求助】请高手帮忙看看，派送大红包。dmol3 模块
已经有18人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
经常有人数据换算后出错，我将网站上的源码做了修改能够适应于MS默认的单位，这样MS7.0版本的直接将cosmo的数据部分保存后导入程序中即可使用。 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Sigma-average_v2.exe
2016-06-30 05:05:26, 1000 K
24楼2016-06-30 05:09:32
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★ ★ ★
gkf高: 金币+4 2014-05-13 16:14:05
gkf高: 金币+4 2014-05-13 16:14:05
|
虽然我对溶剂化了解的不多，昨天刚刚和人讨论过这个问题，赶上得很是时候。 ppt结尾引用了文献Klamt, Andreas, G. Shüürmann的，这个似乎对应的是Dmol3老的溶剂化计算方法（Dmol3的help能够查到），这种方法计算sigma的时候和新方法（默认的ibs）不一样，我们发觉空穴孔洞面积确实不同。如果要在Dmol3中实现Klamt方法，应该先保存input文件（Dmol3 Calculations-(其他参数设定完成)-Files-save files，然后产生input文件），然后把COSMO那一行的ibs修改成on，然后在Dmol3 Calculations-Files-run files运行。 即便这样，我不太确定每个原子的溶剂化半径在两个程序（你给出的ppt和Dmol3）中是否相同。不过仍旧建议你做个测试进行比较。 |
2楼2014-05-13 13:16:48
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2109.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
|
不好意思我那个附件传错了，这地方又传看一边：http://muchong.com/bbs/viewthread.php?tid=7402423 请在帮忙看看，他的COSMO是ON的状态的~ |
3楼2014-05-13 18:25:45
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2014-05-13 19:59:31
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
5楼2014-05-13 20:03:51
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2109.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
|
# Task parameters Calculate energy Symmetry on Max_memory 2048 File_usage smart Basis_Version V4.0.0 # Environment keywords Cosmo on Cosmo_Grid_Size 1082 Cosmo_Segments 92 Cosmo_Solvent_Radius 1.300000 Cosmo_A-Matrix_Cutoff 7.000000 Cosmo_Radius_Incr 0.000000 Cosmo_A-Constant 1.882190 Cosmo_B-Constant 0.010140 Cosmo_RadCorr_Incr 0.150000 Cosmo_Atomic_Radii 1 1.300000 6 2.000000 7 1.830000 8 1.720000 9 1.720000 15 2.1200000 16 2.1600000 17 2.0500000 35 2.1600000 53 2.3200000 # Electronic parameters Spin_polarization restricted Charge 0 Basis dnp Pseudopotential none Functional vwn-bp Aux_density octupole Integration_grid fine Occupation fermi Cutoff_Global 3.7000 angstrom Scf_density_convergence 1.0000e-006 Scf_charge_mixing 2.0000e-001 Scf_iterations 50 Scf_diis 6 pulay # Print options Print eigval_last_it # Calculated properties |
6楼2014-05-13 21:19:21
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★ ★
gmy1990: 金币+3 2014-05-14 03:46:53
gmy1990: 金币+3 2014-05-14 03:46:53
|
刚开始我也出现这个问题，我对照了所有参数发觉都没有问题，后来我猜测是计算P(sigma)的公式不一样导致的，ppt中采用的MS2.2Dmol3单位默认的可能是Angstrom，而新版本MS7.0的Dmol3采用的是a.u.，即长度单位是Bohr，两者换算关系是1Bohr=0.529Angstrom，所以先进行换算单位，然后通过ppt提供的程序分析结果即可。 过程如下： 先把cosmo中的$segment information部分拷贝到记事本propane.txt进行保存，然后导入到origin或是excel，其中有几列，Area那一列要将原来的换算成Area*0.529*0.529，charge/Area直接用chage那一列将原来的值除以新换算的Area即可，Potential那一列要将原来的除以0.529，然后再导出为propane.txt. 然后用ppt提供的程序进行处理，在此之前先在C:\创建一个C:\PROFIlES的文件夹，然后启动程序，该程序首先让你输入完整的文件路径（建议不要是中文的），然后输入生成文件的名称（这个随便），然后输入片段数目（total number of segments，即cosmo文件中的），如例子中的438.然后就可以到C:\PROFIlES去找生成的txt文件，然后将这个txt文件导入origin或excel作图即可。 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Sigma-average_v2.exe
2014-05-14 03:26:23, 344.12 K
» 本帖已获得的红花（最新10朵）
pdl95277楼2014-05-14 03:30:36
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2109.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
8楼2014-05-14 10:00:22
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2109.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
9楼2014-05-14 10:39:11
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3495.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★
gmy1990: 金币+2 2014-10-12 21:28:44
gmy1990: 金币+2 2014-10-12 21:28:44
|
你对照的是下载回来的example？如果是vt2005-00003，那么我发现有一点需要注意，严谨起见，我看了他的输出文件，它采用的Orbital Cutoff不是fine默认的3.7Angstroms，而是5.5Angstroms，即10.397a.u.另外1Bohhr=0.529.......Angstrom，这个也会造成一些偏差。 溶剂化模型了解的不多，说不上来什么，不过有个建议，你可以看看ChemiAndy和sobereva两位前辈的帖子，直接给传送门 http://hi.baidu.com/sobereva/item/5e0b9bd4d49d7ae6b2f777d7 http://muchong.com/bbs/viewthread.php?tid=3965375 以前读了他们的帖子深受启发。你也可以直接问他们。 COSMO/CPCM的模型在你的两份ppt上面所列出的文献你可以读一下，这比问任何人都要有效也更重要。 |
10楼2014-05-14 11:39:32
10