计算化学公社

 找回密码 Forget password
 注册 Register

关于cp2k输出的mulliken文件的分析问题

查看数: 11815 | 评论数: 26 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-4-5 10:30

正文摘要:

下图是我把用cp2k计算MD中的一个片段截取出来做mulliken分析的结果输出,就想请教一下各位大佬,该怎么分析这个输出数据?

回复 Reply

chen0201 发表于 Post on 2023-8-28 21:20:26
sobereva 发表于 2023-8-28 19:59
闭壳层波函数当然不输出自旋布居,全都为0输出有什么意思。Multiwfn不做这么无意义的事

不可能无法打 ...

谢谢老师,已经完成了

结果.png (51.16 KB, 下载次数 Times of downloads: 8)

结果.png
sobereva 发表于 Post on 2023-8-28 19:59:53
chen0201 发表于 2023-8-28 17:22
谢谢老师回复,我是指结果没有得到自旋布居的值。另外因为文件太大了,无法直接打开更改参数导致的,请问 ...

闭壳层波函数当然不输出自旋布居,全都为0输出有什么意思。Multiwfn不做这么无意义的事

不可能无法打开,ultraedit打开几GB的文件也就瞬间的事
chen0201 发表于 Post on 2023-8-28 17:22:20
本帖最后由 chen0201 于 2023-8-28 17:25 编辑
sobereva 发表于 2023-8-28 16:53
语义不明
如果你是指看不到alpha和beta电子布居数,那是因为当前体系是闭壳层的
另外,total net charg ...

谢谢老师回复,我是指结果没有得到自旋布居的值。另外因为文件太大了,无法直接打开更改参数导致的,请问老师,不更改直接看原子布居可以吗?我只是想得到加入体系的额外电子的分布,想看看原子布居数的。
sobereva 发表于 Post on 2023-8-28 16:53:30
chen0201 发表于 2023-8-28 15:49
请问老师,我是用cp2k得到的molden文件按照您的帖子操作的,但是就是得不到您帖子和楼主这样的结果,请问 ...

语义不明
如果你是指看不到alpha和beta电子布居数,那是因为当前体系是闭壳层的
另外,total net charge数值那么离谱,明显你没按照下文对molden文件做恰当的修改
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651http://bbs.keinsci.com/thread-31827-1-1.html
chen0201 发表于 Post on 2023-8-28 15:49:38
本帖最后由 chen0201 于 2023-8-28 15:52 编辑
sobereva 发表于 2020-4-5 15:04
仔细看《原子电荷计算方法的对比》(http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml)了解原 ...

请问老师,我是用cp2k得到的molden文件按照您的帖子操作的,但是就是得不到您帖子和楼主这样的结果,请问是什么原因呢?文件太大,无法上传,附上图片

结果.png (37.61 KB, 下载次数 Times of downloads: 10)

结果.png
lihao2690 发表于 Post on 2023-4-20 14:40:21
sobereva 发表于 2023-4-19 04:27
参考北京科音CP2K第一性原理培训班的ppt

收到,感谢
sobereva 发表于 Post on 2023-4-19 04:27:41
lihao2690 发表于 2023-4-18 15:28
sob老师,每隔一定步数产生一次,关键词怎么写呢?

参考北京科音CP2K第一性原理培训班的ppt


lihao2690 发表于 Post on 2023-4-18 15:28:39
sobereva 发表于 2023-4-2 04:44
每隔一定步数产生一次就行了,比如10步一次。步长通常都很小,在极短尺度内原子电荷变化不会很显著。

sob老师,每隔一定步数产生一次,关键词怎么写呢?
lihao2690 发表于 Post on 2023-4-3 11:05:40
sobereva 发表于 2023-4-2 04:44
每隔一定步数产生一次就行了,比如10步一次。步长通常都很小,在极短尺度内原子电荷变化不会很显著。

好的,谢谢sob老师
sobereva 发表于 Post on 2023-4-2 04:44:30
lihao2690 发表于 2023-4-1 10:30
您说的确实如此,我发现如果每步都产生molden文件,计算速度会降低一倍。我是想看电荷随时间的演化,是不 ...

每隔一定步数产生一次就行了,比如10步一次。步长通常都很小,在极短尺度内原子电荷变化不会很显著。
lihao2690 发表于 Post on 2023-4-1 10:30:04
sobereva 发表于 2023-4-1 05:18
显然不应当每一步都产生,又降低速度又浪费硬盘又毫无必要

详谈Multiwfn的命令行方式运行和批量运行的 ...

您说的确实如此,我发现如果每步都产生molden文件,计算速度会降低一倍。我是想看电荷随时间的演化,是不是必须每一步都生成molden呢?如果是这样,我如何提高运行效率呢?谢谢sob老师
sobereva 发表于 Post on 2023-4-1 05:18:45
lihao2690 发表于 2023-3-30 10:29
谢谢sob老师,发现跑MD的话产生molden文件会使计算速率降低一倍左右,另外每一步会产生一个molden文件, ...

显然不应当每一步都产生,又降低速度又浪费硬盘又毫无必要

详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612http://bbs.keinsci.com/thread-24929-1-1.html
lihao2690 发表于 Post on 2023-3-30 10:29:04
sobereva 发表于 2022-7-28 19:16
用Multiwfn基于CP2K的波函数算CM5电荷,远比Mulliken电荷强得多

按下文产生molden文件

谢谢sob老师,发现跑MD的话产生molden文件会使计算速率降低一倍左右,另外每一步会产生一个molden文件,这些如何在multiwfn中进行批处理呢?谢谢
sobereva 发表于 Post on 2023-3-24 05:27:32
lihao2690 发表于 2023-3-20 22:24
您好,请问输入什么关键词能够给出mulliken电荷呢

Multiwfn创建CP2K输入文件界面里有个选项是计算原子电荷,里面选上Mulliken。另外,PRINT_LEVEL MEDIUM的时候也默认会输出。顺带一提,用Multiwfn载入CP2K的molden文件做Mulliken布居分析,能得到远比Mulliken原子电荷更丰富的信息。

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

GMT+8, 2024-11-26 03:36 , Processed in 0.168365 second(s), 25 queries , Gzip On.

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