计算化学公社

标题: cp2k可以计算出的自旋向上和向下的态密度一样，怎么办？ [打印本页]

作者
Author: TJU_RYC 时间: 2024-4-15 19:41
标题: cp2k可以计算出的自旋向上和向下的态密度一样，怎么办？
本帖最后由 TJU_RYC 于 2024-4-17 19:20 编辑

想计算出类似这样的态密度但查了好几天就只会这个，但计算出来没有区分自旋
现在会计算自旋了但是一样的，不知道怎么破，希望会的老师或师兄师姐指点一二，谢谢啦

下面是输入文件，我希望能算出和别人文献中的一样，自旋上下态密度不一样

作者
Author: sobereva 时间: 2024-4-16 03:17
Multiwfn基于CP2K的Molden文件，对开壳层体系绘制TDOS/PDOS/OPDOS时可以直接选对什么自旋绘制，参考下文
使用Multiwfn绘制态密度(DOS)图考察电子结构
http://sobereva.com/482（http://bbs.keinsci.com/thread-13242-1-1.html）
详谈使用CP2K产生给Multiwfn用的molden格式的波函数文件
http://sobereva.com/651（http://bbs.keinsci.com/thread-31827-1-1.html）
北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里也都给了具体绘制态密度的例子

最新版的Multiwfn直接载入CP2K输出的含有各个k点能级信息的输出文件，输入cp、5，界面里也可以选择对什么自旋绘制DOS，但目前仅限TDOS

CP2K自己的&PDOS基本没有用处。

作者
Author: TJU_RYC 时间: 2024-4-16 08:46

sobereva 发表于 2024-4-16 03:17
Multiwfn基于CP2K的Molden文件，对开壳层体系绘制TDOS/PDOS/OPDOS时可以直接选对什么自旋绘制，参考下文
...

老师，最新版指的是2024.1吗？我现在是2023.2需要更新不？

作者
Author: TJU_RYC 时间: 2024-4-16 09:42

sobereva 发表于 2024-4-16 03:17
Multiwfn基于CP2K的Molden文件，对开壳层体系绘制TDOS/PDOS/OPDOS时可以直接选对什么自旋绘制，参考下文
...

看歪了，我先试试，谢谢老师

作者
Author: highlight 时间: 2024-4-16 14:25
你用的 ROKS 吧，加了 LSD 应该自己就分开了呀

作者
Author: TJU_RYC 时间: 2024-4-16 14:44

highlight 发表于 2024-4-16 14:25
你用的 ROKS 吧，加了 LSD 应该自己就分开了呀

老师，我看了一下没有发现有ROKS字段，我上传了输入文件（网上找的案例），可以麻烦您帮我看看吗？

作者
Author: highlight 时间: 2024-4-16 16:38
&FORCE_EVAL
METHOD Quickstep
&DFT
   LSD
   BASIS_SET_FILE_NAME  BASIS_MOLOPT
   ......

作者
Author: TJU_RYC 时间: 2024-4-16 20:22
自旋向上和向下的态密度是一样的

作者
Author: TJU_RYC 时间: 2024-4-16 20:22

highlight 发表于 2024-4-16 16:38
&FORCE_EVAL
METHOD Quickstep
&DFT

感谢老师，成功了，只是又有个问题现在自旋向上和向下的态密度是一样的，是哪里错了呀，麻烦老师在帮我看看？

作者
Author: highlight 时间: 2024-4-17 15:05
先确定你算的物质，是不是应该有不同自旋的态密度的分离，你的 d/e 两幅图不也没有吗？
要是确定有，可以试试改变初猜，我记得论坛里有几个帖子讨论过

作者
Author: TJU_RYC 时间: 2024-4-17 19:23

highlight 发表于 2024-4-17 15:05
先确定你算的物质，是不是应该有不同自旋的态密度的分离，你的 d/e 两幅图不也没有吗？
要是确定有，可以 ...

谢谢老师，那我先看看论坛的帖子，我也更新了我的输入文件和问题描述，老师有空的话或有也帮我看看呗，谢谢啦

欢迎光临计算化学公社 (http://bbs.keinsci.com/)