捐赠本论坛 (Donate this site)
开启辅助访问 切换到宽版

计算化学公社

 找回密码 Forget password
 注册 Register
在11月21-24日期间，本论坛仅限等级≥level 4的成员发贴，带来的不便请谅解！其他成员若有学术问题请等到能发贴的时候再发帖。此期间不要给管理员sobereva发私信问学术问题
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › 为什么非高对称性K点的轨道系数也是实数？
返回列表 Return 发新帖
Views: 39|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 为什么非高对称性K点的轨道系数也是实数？

[复制链接 Copy URL]
Freeman

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)
跳转到指定楼层 Go to specific reply
楼主
发表于 Post on 2024-11-8 15:58:12 | 只看该作者 Only view this author 回帖奖励 |倒序浏览 Reverse view |阅读模式 Reading model
大家好。

我知道高对称性的K点（例如Gamma点）的轨道系数可以是实数，但其他K点的则一般是复数。可是我用
  1. &PRINT
  2.   &MO
  3.     COEFFICIENTS
  4.     FILENAME MO
  5.   &END MO
  6. &END PRINT
复制代码

关键词打印出来的轨道系数文件，所有K点的轨道系数都是实数。这是为什么呢？

K点信息如下：(单位是2pi/bohr)
  1. K-point Symmetry Information
  2.        Number of Special K-points:                                            14
  3.                    Wavevector Basis                                       Weight
  4.           1  -0.33333  -0.33333  -0.33333                                      2
  5.           2  -0.33333  -0.33333   0.00000                                      2
  6.           3  -0.33333  -0.33333   0.33333                                      2
  7.           4  -0.33333   0.00000  -0.33333                                      2
  8.           5  -0.33333   0.00000   0.00000                                      2
  9.           6  -0.33333   0.00000   0.33333                                      2
  10.           7  -0.33333   0.33333  -0.33333                                      2
  11.           8  -0.33333   0.33333   0.00000                                      2
  12.           9  -0.33333   0.33333   0.33333                                      2
  13.          10   0.00000  -0.33333  -0.33333                                      2
  14.          11   0.00000  -0.33333   0.00000                                      2
  15.          12   0.00000  -0.33333   0.33333                                      2
  16.          13   0.00000   0.00000  -0.33333                                      2
  17.          14   0.00000   0.00000   0.00000                                      1

  19.        K-point Mesh:                                               3     3     3
  20.                    Wavevector Basis                   Special Points    Rotation
  21.           1  -0.33333  -0.33333  -0.33333              1           1           1
  22.           2  -0.33333  -0.33333   0.00000              2           2           1
  23.           3  -0.33333  -0.33333   0.33333              3           3           1
  24.           4  -0.33333   0.00000  -0.33333              4           4           1
  25.           5  -0.33333   0.00000   0.00000              5           5           1
  26.           6  -0.33333   0.00000   0.33333              6           6           1
  27.           7  -0.33333   0.33333  -0.33333              7           7           1
  28.           8  -0.33333   0.33333   0.00000              8           8           1
  29.           9  -0.33333   0.33333   0.33333              9           9           1
  30.          10   0.00000  -0.33333  -0.33333             10          10           1
  31.          11   0.00000  -0.33333   0.00000             11          11           1
  32.          12   0.00000  -0.33333   0.33333             12          12           1
  33.          13   0.00000   0.00000  -0.33333             13          13           1
  34.          14   0.00000   0.00000   0.00000             14          14           1
  35.          15   0.00000   0.00000   0.33333             13          13           2
  36.          16   0.00000   0.33333  -0.33333             12          12           2
  37.          17   0.00000   0.33333   0.00000             11          11           2
  38.          18   0.00000   0.33333   0.33333             10          10           2
  39.          19   0.33333  -0.33333  -0.33333              9           9           2
  40.          20   0.33333  -0.33333   0.00000              8           8           2
  41.          21   0.33333  -0.33333   0.33333              7           7           2
  42.          22   0.33333   0.00000  -0.33333              6           6           2
  43.          23   0.33333   0.00000   0.00000              5           5           2
  44.          24   0.33333   0.00000   0.33333              4           4           2
  45.          25   0.33333   0.33333  -0.33333              3           3           2
  46.          26   0.33333   0.33333   0.00000              2           2           2
  47.          27   0.33333   0.33333   0.33333              1           1           2
复制代码

轨道系数如下

  1. MO| EIGENVALUES, OCCUPATION NUMBERS, AND SPHERICAL EIGENVECTORS FOR K POINT 1
  2. MO|
  3. MO|                             1            2            3            4
  4. MO|                         -0.726292    -0.726278    -0.726278    -0.726278
  5. MO|
  6. MO|                          2.000000     2.000000     2.000000     2.000000
  7. MO|
  8. MO|    1     1 Pb  2s        0.000098     0.000118     0.000167     0.000119
  9. MO|    2     1 Pb  3s       -0.000042    -0.000063    -0.000089    -0.000063
  10. MO|    3     1 Pb  3py       0.000022     0.000003     0.000000    -0.000003
  11. MO|    4     1 Pb  3pz       0.000022    -0.000083     0.000000     0.000082
  12. MO|    5     1 Pb  3px       0.000022     0.000003     0.000000    -0.000003
  13. 。。。。。。
  14. MO| EIGENVALUES, OCCUPATION NUMBERS, AND SPHERICAL EIGENVECTORS FOR K POINT 2
  15. MO|
  16. MO|                             1            2            3            4
  17. MO|                         -0.726298    -0.726283    -0.726283    -0.726269
  18. MO|
  19. MO|                          2.000000     2.000000     2.000000     2.000000
  20. MO|
  21. MO|    1     1 Pb  2s       -0.000146     0.000145    -0.000145     0.000024
  22. MO|    2     1 Pb  3s        0.000045    -0.000050     0.000050    -0.000011
  23. MO|    3     1 Pb  3py      -0.000022    -0.000002    -0.000041     0.000083
  24. MO|    4     1 Pb  3pz      -0.000000     0.000000     0.000000    -0.000000
  25. MO|    5     1 Pb  3px      -0.000022     0.000041     0.000002     0.000083
  26. MO|    6     1 Pb  4py      -0.000004     0.000006    -0.000027     0.000048
  27. 。。。。。。
  28. MO| EIGENVALUES, OCCUPATION NUMBERS, AND SPHERICAL EIGENVECTORS FOR K POINT 3
  29. MO|
  30. MO|                             1            2            3            4
  31. MO|                         -0.726292    -0.726278    -0.726278    -0.726278
  32. MO|
  33. MO|                          2.000000     2.000000     2.000000     2.000000
  34. MO|
  35. MO|    1     1 Pb  2s        0.000183     0.000122     0.000123    -0.000122
  36. MO|    2     1 Pb  3s       -0.000078    -0.000065    -0.000065     0.000065
  37. MO|    3     1 Pb  3py       0.000009     0.000006    -0.000002    -0.000076
  38. MO|    4     1 Pb  3pz      -0.000009     0.000084     0.000000     0.000000
  39. MO|    5     1 Pb  3px       0.000009     0.000005     0.000076     0.000002
  40. MO|    6     1 Pb  4py      -0.000001    -0.000012     0.000014    -0.000051
复制代码

pristine.inp (7.73 KB, 下载次数 Times of downloads: 0)


回复 Reply

举报 Report

sobereva

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长
2#
发表于 Post on 2024-11-9 06:07:32 | 只看该作者 Only view this author
展开系数只有一套，只不过实部和虚部的前因子不一样而已。
北京科音自然科学研究中心（http://www.keinsci.com）致力于计算化学的发展和传播，长期开办高质量的各种计算化学类培训：初级量子化学培训班、中级量子化学培训班、高级量子化学培训班、量子化学波函数分析与Multiwfn程序培训班、分子动力学与GROMACS培训班、CP2K第一性原理计算培训班，内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路！培训各种常见问题见《北京科音办的培训班FAQ》。
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章！
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”：1号：18616395，2号：466017436，3号：764390338，搜索群号能搜到哪个说明目前哪个能加，合计9000人。北京科音培训班的学员在群中可申请VIP头衔，提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog：http://sobereva.com（发布大量原创计算化学相关博文）
Multiwfn主页：http://sobereva.com/multiwfn（十分强大的量子化学波函数分析程序）
ResearchGate：https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!
回复 Reply

举报 Report

返回列表 Return 发新帖

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 06:41 , Processed in 0.169261 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list