计算化学公社

 找回密码 Forget password
 注册 Register
Views: 45|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
大家好。

我知道高对称性的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

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


5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

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!

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 01:04 , Processed in 0.163366 second(s), 24 queries , Gzip On.

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