卡开发发 发表于 2024-3-22 14:45
看上去还算合理，因为DOS处理/绘制的时候牵扯展宽。

谢谢老师！

不会扣篮的后卫 发表于 2024-3-22 11:39
老师，您好。关于之前画片段DOS图的问题，开发者回复是QE老版本分析时侯要和计算的时候使用相同数量的CPU ...

看上去还算合理，因为DOS处理/绘制的时候牵扯展宽。

卡开发发 发表于 2024-3-1 14:32
正常来说，用一个整体的态密度投影到原子轨道上，当我们想知道一些片段的贡献时，用这组投影的结果直接求 ...

感谢老师的解答！我再去和程序开发者讨论一下。

不会扣篮的后卫 发表于 2024-3-1 14:17
所以老师您的意思是该程序只能画出总的DOS图，不能画出pDOS图？

正常来说，用一个整体的态密度投影到原子轨道上，当我们想知道一些片段的贡献时，用这组投影的结果直接求和来估计片段的贡献（虽然因为一些划分方式存在问题可能不严格具备加合性），这样做是没问题的。但是如果你要把整个体系拆分成各个片段分别再进行自洽，存在这样的问题：
1、本身各个片段之间的能级相对位置就不定（因为电子数不同于整体的情况），此时大家都平移到Fermi能级，这样实际的能量标准都是不同的。
2、重新自洽的后的能级和各个能级的占据数也不同于整体的，或者反过来说，片段之间有可能在组合后形成成键，亦即能级发生分裂生成成键和反键，而此时密度也发生变化（也就是对应着整体和片段和的电子密度差）。

当然你计算的结果却是存在着一些很怪的问题，例如你的各个片段中fdos.dat.pdos_tot的数据DOS和PDOS做出来差异似乎很大，这个本身似乎也不太正常。

卡开发发 发表于 2024-3-1 14:04
那样可能不太行，一个是各个片段重新自洽之后Fermi能级其实会发生移动，另外如果牵扯到成键，其实能级的 ...

所以老师您的意思是该程序只能画出总的DOS图，不能画出pDOS图？

不会扣篮的后卫 发表于 2024-3-1 12:59
不好意思，老师。我之前的信息没有介绍清楚。

1，8个片段组成了这个体系。通过基于QE5.1.0改写的程序 ...

那样可能不太行，一个是各个片段重新自洽之后Fermi能级其实会发生移动，另外如果牵扯到成键，其实能级的相对位置也会发生偏移，这样原理上就不具备加合性。

卡开发发 发表于 2024-3-1 12:32
我不太确定，大概这几个问题：
1、你的8个片段和整体是什么关系？是这8个片段是整体划分出来的？
2、这 ...

不好意思，老师。我之前的信息没有介绍清楚。

1，8个片段组成了这个体系。通过基于QE5.1.0改写的程序来切割周期性的体系，变成8个周期性的小片段。
2，片段的PDOS获取的过程是各自的片段重新进行了自洽计算再各自计算PDOS。
3，8个片段的fdos.dat.pdos_atm#xxxxxx不能各自能加和出来对应的fdos.dat.pdos_tot的结果。

不会扣篮的后卫 发表于 2024-3-1 12:13
老师，您好。我按照您之前的建议，用了这个输入文件，

&PROJWFC

我不太确定，大概这几个问题：
1、你的8个片段和整体是什么关系？是这8个片段是整体划分出来的？
2、这些片段的PDOS获取的过程是怎样？是各自的片段重新进行了自洽计算再各自计算PDOS，还是这些片段的PDOS直接从整体通过某种算法投影出来的？
3、8个片段的fdos.dat.pdos_atm#xxxxxx是否各自能加和出来对应的fdos.dat.pdos_tot的结果？

卡开发发 发表于 2024-2-21 08:03
1、我不确定DOS能量范围更低是否有更大的值。
2、如果能量范围足够宽，你可以积分到Fermi能级看看是否和赝 ...

谢谢老师的建议，我再去研究一下。

1、我不确定DOS能量范围更低是否有更大的值。
2、如果能量范围足够宽，你可以积分到Fermi能级看看是否和赝势的价电子数一致，因为这里的PDOS是按照Lowdin划分过的，本来原子的价电子数目也不多所以都会比较小。
3、你可以看看总DOS的情况，然后用sumpdos.x的工具把PDOS加加看看看差距多大。