| 查看: 1824 | 回复: 2 | ||||
[交流]
投影态密度有负值
|
| 本人用smeagol算两电极系统的投影态密度时，某些原子上的PDOS出现负值，请问有人碰到这种情况吗？又如何解释呢？谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 石墨烯的总态密度与分波态密度
已经有25人回复
- DOSCAR画态密度
已经有16人回复
- JPCC修回意见中提到关于能级的问题，不甚理解，求指点!!!
已经有13人回复
- abinit能否处理一个位置被两种元素占据的固溶体的能带和态密度
已经有10人回复
- 有关态密度图像
已经有5人回复
- 【求助成功】怎么画出5个d轨道分别对应的态密度
已经有23人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【求助】abinit如何能得到m分量的投影态密度？
已经有4人回复
- 【求助】求助：由态密度图如此推导的导电机制，如何理解？
已经有19人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- 【求助】MS５.０ Castep 计算不出部分态密度（PDOS）怎么办？（站内有人提问MS4.3，�
已经有19人回复
- 电子结构和态密度以及能带之间有什么关系?
已经有6人回复
- 【求助】在siesta中怎么画电势，轨道，投影态密度的空间分布【完结】
已经有5人回复
- 【求助】怎样由能带结构和态密度看计算的结构是不是可能存在超导现象
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
2026年材料结构与力学性能国际会议（MSMP 2026
+1/826
-
2027年航天航空、材料与人工智能国际会议（AMAI 2027）
+1/786
-
2027年电气工程、自动化与应用工程国际会议（EEAE 2027）
+1/775
-
专利代理人
+1/97
-
中国科大-合肥国家实验室冷原子量子网络团队招聘启事
+2/86
-
南方科技大学力航系急招收纳米材料、力学、机械等方向推免硕士研究生两名
+1/66
-
华南师范大学胡勇军教授课题组招收2027年推荐面试攻读研究生（推免生）
+1/36
-
中山大学柔性电子方向2027级博士招生（材料、材料加工、化学、化工、机械等背景）
+1/20
-
浙江海洋大学2027年硕士研究生招生专业目录及自命题科目考试大纲
+1/14
-
废油资源化教育部工程中心 诚聘催化方向特聘研究员
+1/10
-
2027 International Conference on Logistics and Supply Chain Management
+1/7
-
【2027硕士招生】招收材料、化工、环境、冶金类硕士研究生
+1/5
-
实验室折腾了三年的自制设备要不要申请专利？跟导师争了一晚上没吵明白
+1/5
-
2027 International Conference on Healthcare, Medicine and Medical Science
+1/4
-
哈尔滨工业大学王志江教授课题组招收博士生
+1/3
-
澳门大学光电方向招收博士研究生
+1/3
-
浙江大学Said Mikki教授课题组招生2027秋博士生-电子科学与技术
+1/3
-
2027 World Conference on Economics, Finance and Management Sciences
+1/3
-
NK细胞：固有免疫体系中的核心杀伤细胞
+1/3
-
2027 World Conference on Business, Marketing and Communication (WCBMC 2027)
+1/2
★
futiliu(金币+1):谢谢参与
WDD880227(金币+1): 感谢交流提示 2012-02-28 15:35:57
futiliu(金币+1):谢谢参与
WDD880227(金币+1): 感谢交流提示 2012-02-28 15:35:57
|
这个是由于基组的不正交性导致的。严格说来，如果投影的时候采用一个已经正交的基组来定义的话，就会避免这个问题，但是那个时候PDOS的物理意义没有现在明确。所以有点“鱼和熊掌不可兼得”。 更加详细的解释，我下面引用如下： http://www.mail-archive.com/sies ... am.es/msg02951.html ------------------------------------------ Dear Siesta users, I have a problem understanding the output in a .PDOS output file for my system. I am getting negative values in the PDOS for the projection onto a lot of the basis functions. Perhaps my understanding is wrong, but surely this should not be? According to my understanding the PDOS onto a specific basis function should have peaks at the energies corresponding to the eigenfunctions to which it contributes. I restarted the PDOS calculation from the output files of a converged calculation and used 1 SCF iteration with 50x50x50 Monkhorst Pack grid. Regards, Rainer ------------------------------------------- Dear Rainer: When you are working with a non-orthogonal basis set, as is the case in Siesta, neither the PDOS nor the Mulliken population analysis are positive definite magnitudes. The reason is due to the fact that the PDOS is defined as: g_mu(eps) = sum_nu rho_mu,nu S_nu,mu delta(eps-eps_i) rho_mu,nu = sum_i C_mu,i C*_nu,i (i stands for the bands, S stands for the overlap matrix, and C stands for the coefficients of the wave function) and the non-diagonal elements of the density matrix might be negative. However, the term in the diagonal IS positive definite, and it is usually larger than the rest of the terms. Therefore, when the PDOS is negative the value should be small in absolute value. Hope this helps, Javier --------------------------------------- Dear Javier Thanks very much for your reply. This certainly helps me understand why there can be negative values in the PDOS, at least mathematically. But these negative values of PDOS are not physically meaningful, right? So how should I interpret them? Can one just invert them? I am getting negative values that are roughly the same size as the positive values, so I dont feel comfortable just disregarding them. Regards, Rainer ----------------------------------------- Dear Rainer: I do not think you can give any physical meaning to a large negative Projected Density Of States. The problem is rather fundamental. The solution should come from a redefinition of the PDOS as the projection over a basis of orthogonalized orbitals. That would imply a first transformation from our non-orthogonal basis set to a orthogonal one (I guess that you will have to multiply by the inverse of the overlap matrix somewhere). Hope this helps, Javier ================================= |
2楼2012-02-28 14:55:07
3楼2016-10-26 21:59:41
