| 查看: 205 | 回复: 0 | |||
| 当前主题已经存档。 | |||
WATERpml木虫 (小有名气)
|
[交流]
【讨论】basis set不同，计算得到的单分子的unoccupied states会不同吗？
|
||
|
When comparing the DOS of isolated molecules as obtained by VASP with the results of other codes which use local basis sets (Siesta, Turbomole), it seems that the states match only up to a certain energy below the vacuum energy. For higher energies (but still below the vacuum energy), neither the number nor the position of the orbitals match necessarily. I focused on the molecule 4-Pyridin-4-yl-phenylamine. Using parameters as comparable as possible, I see 4 unoccupied states below the vacuum energy with Siesta and Turbomole (at approximately the same energies). Using VASP, the DOS matches those results very nicely up to the LUMO+2. 3D-representations of the orbitals are consistent as well. Between the states LUMO+2 and LUMO+3 (according to Siesta and Turbomole), VASP calculates an additional state. By looking at its 3d-representation I can see that this state looks like what is the LUMO+4 (the first state above the vacuum energy) in siesta and turbomole. A very quick test with another plane-wave code, abinit, finds a DOS very similar to what VASP finds. To sum up, in this example I see HOMO LUMO L+1 L+2 L+3 VacEn. L+4 in siesta and turbomole. Vasp finds: HOMO LUMO L+1 L+2 ex L+3 VacEn. where the peak "ex" corresponds to L+4 from siesta and turbomole, and VacEn is the vacuum energy. From 3d representations of the L+3 state I don't think they represent the same orbital. Is it possible that plane-wave and local-basis set codes behave somehow differently concerning unoccupied orbitals? If so, why? There is a related thread, http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?4.6181 There, however, it seems that the calculation of states above the vacuum energy was the issue (which, I guess, could in principle be solved by setting NBANDS accordingly). I want to emphasize that I see (for at least three molecules i looked at so far) a different DOS also in a certain range of unoccupied states below the vacuum energy. Thanks in advance PS：此问题我们也在VASP Forum发贴：http://cms.mpi.univie.ac.at/vasp ... pic.php?4.6486#6486，请各位高手不吝赐教！ |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
