计算化学公社

 找回密码 Forget password
 注册 Register

【原创】分享一个基于DMol3计算电子密度差的Matlab程序(未经授权,禁止转载)

查看数: 8336 | 评论数: 12 | 收藏 Add to favorites 20
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
wsz
发布时间: 2022-4-8 02:37

正文摘要:

本帖最后由 wsz 于 2022-4-8 02:52 编辑 程序(见附件)及示例文件:「ElectronDensityDifference4DMol3」https://www.aliyundrive.com/s/KbmTSqbBiWm 程序使用的电子密度差的计算公式:

回复 Reply

wsz 发表于 Post on 2024-3-28 17:52:04
polymer_ck 发表于 2024-3-28 16:23
老师,您的代码帮了大忙了,请问如何致谢或引用您的工作?

不需要引用,可以直接用哈
polymer_ck 发表于 Post on 2024-3-28 16:23:23
老师,您的代码帮了大忙了,请问如何致谢或引用您的工作?
wsz 发表于 Post on 2024-3-21 16:33:50
丁真xxx 发表于 2024-3-21 15:52
老师,Dmol3我看是不能计算部分的差分点密度,请问您是怎么得到三个文件,

分别自洽得到
丁真xxx 发表于 Post on 2024-3-21 15:52:49
老师,Dmol3我看是不能计算部分的差分点密度,请问您是怎么得到三个文件,
sdc5070163 发表于 Post on 2023-12-6 13:58:24
wsz 发表于 2023-11-30 12:56
把这几个脚本和相关文件放到同一个目录,修改Main文件import_grd对应的文件名(4~6行),直接run就行。

感谢老师
wsz 发表于 Post on 2023-11-30 12:56:14
sdc5070163 发表于 2023-11-27 19:39
老师您好 非常感谢您的分享 不过是否可以提供一下详细的操作流程?对Matlab实在是操作不起来

把这几个脚本和相关文件放到同一个目录,修改Main文件import_grd对应的文件名(4~6行),直接run就行。
sdc5070163 发表于 Post on 2023-11-27 19:39:52
老师您好 非常感谢您的分享 不过是否可以提供一下详细的操作流程?对Matlab实在是操作不起来
Pengplus 发表于 Post on 2023-9-10 14:10:08
对于CASTEP的差分电荷,我的做法是纯手工: 下载三个计算结果文件(.check文件 在remote review中有)→c2x转化check文件为.xsf文件(http://www.c2x.org.uk/)→VESTA打开处理  即可得到VESTA画风的差分电荷图   每次操作都略显繁琐但是又无可奈何
kaikai 发表于 Post on 2023-9-7 13:39:02
wsz 发表于 2023-9-7 10:33
没有专门发表,可以直接用哈。

多谢老师
wsz 发表于 Post on 2023-9-7 10:33:59
kaikai 发表于 2023-9-7 08:35
老师,你的这篇文章中介绍的方法帮了我的大忙,请问这个程序发表到文章中了吗?我想引用一下。

没有专门发表,可以直接用哈。
kaikai 发表于 Post on 2023-9-7 08:35:38
老师,你的这篇文章中介绍的方法帮了我的大忙,请问这个程序发表到文章中了吗?我想引用一下。
sobereva 发表于 Post on 2022-4-8 15:45:47
顺带一提,Multiwfn也可以直接载入grd文件,之后可以用主功能13的子功能11对当前载入的和另一个grd文件做各种数学运算(20多种,求差包括在内),算完的格点数据可以直接在主功能13里用选项-2观看等值面,也可以导出成cub文件用其它程序绘制。

另外,Multiwfn的settings.ini里的iuserfunc设为-1(线性插值)或-3(B-spline插值)后,在启动并载入一个grd文件后,可以用主功能3和主功能4分别在特定路径上和特定平面上绘制曲线图和各种形式的平面图,数值对应载入的格点数据插值出的数值。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 23:53 , Processed in 0.212732 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list