fd7 发表于 2024-9-28 17:58
老师，请问“C和O上的电子密度整体大幅减少，很大程度往Na的方向上转移了”这句话如果这样理解是否正确：电 ...

是

老师，请问“C和O上的电子密度整体大幅减少，很大程度往Na的方向上转移了”这句话如果这样理解是否正确：电子流动方向是从CO流向了Na，Na的正电荷会降低。

今天才发现 使用Multiwfn处理 cp2k-hartree-v_hartree-1_0 .cub文件获得平面平均曲线（plane-averaged curve)减去费米能就是功函数

yjb 发表于 2022-5-5 20:46
老师，Multiwfn可以画cp2k计算出来的cv的FES的图吗？ 谢谢老师

不能

老师，Multiwfn可以画cp2k计算出来的cv的FES的图吗？ 谢谢老师

sobereva 发表于 2022-4-15 02:58
预计第四季度有
www.keinsci.com首页预告栏都体现得很清楚

好的，到时一定参加，谢谢卢老师

静哥哥 发表于 2022-4-14 17:15
好的，谢谢卢老师。今天翻看了您关于杂化的帖子，真是太佩服了。不知道卢老师近期有没有计划举办cp2k的培 ...

预计第四季度有
www.keinsci.com首页预告栏都体现得很清楚

sobereva 发表于 2022-4-14 15:14
PBE0的HF杂化成份不是经验值，看原文，是理论推出来的。这点和很多其它泛函通过拟合、经验确定的HF成份截 ...

好的，谢谢卢老师。今天翻看了您关于杂化的帖子，真是太佩服了。不知道卢老师近期有没有计划举办cp2k的培训班？

静哥哥 发表于 2022-4-13 14:03
谢谢社长。
另外请教个问题，根据这个帖子https://www.cp2k.org/exercises:2015_ethz_mmm:tio2_gap
关 ...

PBE0的HF杂化成份不是经验值，看原文，是理论推出来的。这点和很多其它泛函通过拟合、经验确定的HF成份截然不同。
绝对不要用1，100% HF成份对于绝大多数问题都很不适合。
另外，那个网页里的做法属于拿可调参数凑数据，不建议那么搞。

HF成份是泛函直接定义的，看
不同DFT泛函的HF成份一览
http://sobereva.com/282（http://bbs.keinsci.com/thread-807-1-1.html）

适合用什么泛函看具体目的，比如对于带隙，HSE06就整体不错，而且是公认的，拿这个算带隙不会招致太大质疑。而自己凑一个HF成份，而且说不出来什么物理层面的意义，则明显更容易招质疑。

sobereva 发表于 2022-4-13 01:41
估计是表面区域电子结构复杂导致的
可以试试不同切的方式（界面处的原子怎么归属给两部分）
也可以尝试 ...

谢谢社长。
另外请教个问题，根据这个帖子https://www.cp2k.org/exercises:2015_ethz_mmm:tio2_gap
关于PBE0杂化函数中的HFX交换系数，是一个半经验值，自己改变并使结构符合自己的目标能带值。比如fe2o3我用hfx=1能很好地模拟气电子结构，比如dos和gap。那么，对于fe3o4或者fe2o3的缺陷模型也能用1吗？我大概搜了下，好像不能直接用1。对于不同体系这个系数可能有出入，那么我该怎么处理这个系数呢。抱歉，这块理论比较欠缺，请社长解惑下

静哥哥 发表于 2022-4-12 21:22
社长大大：
    打扰了。我按照您的帖子重复了下您的文件，都没问题，能很好的得到差分电荷密度。
    但 ...

估计是表面区域电子结构复杂导致的
可以试试不同切的方式（界面处的原子怎么归属给两部分）
也可以尝试CP2K版块里讨论的帮助收敛的方法，之前我也回复过很多次

上面的异质结切去一部分，我是随意切的，只是用于向您展示

社长大大：
    打扰了。我按照您的帖子重复了下您的文件，都没问题，能很好的得到差分电荷密度。
    但是涉及我自己的异质结构建结区差分电荷密度就出现了问题。我首先构建了异质结C=A+B，类似tio2孪晶，是3d模型，用multiwfn得到c.inp文件，cp2k结构优化后，再用energy进行计算得到cube文件，几十步就收敛了，这里一切正常。然后，我用multiwfn显示异质结C，记下A和B的原子位置，在c.inp里分别删除A或B 的原子。再次进行计算，想得到B和A的cube文件。但是，我发现在A和B进行energy 计算是，SCF无法收敛，跑了几百步都不收敛，我分别调节了 算法，MIXING下的各个参数，发现A或B都无法收敛。我的inp文件是默认的pbe.inp文件，按理说就算只剩A或者B应该很快收敛，可是波动很大。对于计算差分电荷密度，我又不能对A或者B再次结构优化。
     我和另外一个老师讨论了半天，觉得是不是因为我是直接切割的3D模型的cif，产生了很多悬挂键，导致不收敛？请教社长大大，我该怎么处理

今日对本文和Multiwfn程序进行了更新，将局部积分曲线和平面平均曲线做了区分，后者定义为前者除以截面面积。最新版Multiwfn的本文介绍的功能界面上已经多了计算plane-average curve的选项。