| 查看: 324 | 回复: 0 | |||
| 当前主题已经存档。 | |||
hbnuwangshuling铜虫 (小有名气)
|
[交流]
【求助】对PBC计算时的结果如何处理？
|
||
|
PBC计算时的结果是什么意思？G03开始，Gaussian就开始支持周期性边界条件了。具体的关键词是tv,也就是translation vector的首字母缩写。完整输入文件如下（根据Gaussian手册，在PBC情况下，blyp之类的算法比HF更快）： #p blyp/6-311G**(d) PBC(NKPoint=20) iop(5/33=3) C atom 0,1 C 0.00 0.00 0.00 TV 0.00 0.00 1.44 此输入文件说明我们模拟的是一维无穷长的碳链，碳原子之间的距离为1.44A。为了简单起见，我们不考虑由于派尔斯相变引起的构型变化。 用如上的输入文件完成计算后，在输出文件的SCF done和Population Analysis之间会有这样一段输出 3 alpha orbitals occupied out of 30 functions. alpha orb energies kxyz= 0 0 0 0.000 0.000 -10.024 -0.741 -0.344 | -0.344 0.261 0.314 0.314 0.445 eV alpha orb energies kxyz= 40 1 1 0.000 0.000 -10.023 -0.486 -0.056 | 0.087 0.087 0.600 0.600 0.736 eV alpha orb energies kxyz= 39 1 1 0.000 0.000 -10.023 -0.491 -0.047 | 0.074 0.074 0.609 0.609 0.691 eV alpha orb energies kxyz= 38 1 1 0.000 0.000 -10.023 -0.503 -0.031 | 0.051 0.051 0.620 0.620 0.645 eV alpha orb energies kxyz= 37 1 1 0.000 0.000 -10.023 -0.518 -0.007 | 0.021 0.021 0.602 0.621 0.621 eV alpha orb energies kxyz= 36 1 1 0.000 0.000 -10.024 -0.536 -0.013 | -0.013 0.021 0.563 0.606 0.606 eV alpha orb energies kxyz= 35 1 1 0.000 0.000 -10.024 -0.555 -0.048 | -0.048 0.053 0.529 0.579 0.579 eV alpha orb energies kxyz= 34 1 1 0.000 0.000 -10.024 -0.576 -0.083 | -0.083 0.087 0.501 0.546 0.546 eV alpha orb energies kxyz= 33 1 1 0.000 0.000 -10.024 -0.596 -0.117 | -0.117 0.122 0.484 0.512 0.512 eV alpha orb energies kxyz= 32 1 1 0.000 0.000 -10.024 -0.616 -0.150 | -0.150 0.155 0.477 0.478 0.478 eV alpha orb energies kxyz= 31 1 1 0.000 0.000 -10.024 -0.635 -0.181 | -0.181 0.185 0.447 0.447 0.484 eV alpha orb energies kxyz= 30 1 1 0.000 0.000 -10.024 -0.654 -0.210 | -0.210 0.211 0.418 0.418 0.501 eV alpha orb energies kxyz= 29 1 1 0.000 0.000 -10.024 -0.670 -0.236 | -0.236 0.232 0.392 0.392 0.526 eV alpha orb energies kxyz= 28 1 1 0.000 0.000 -10.024 -0.686 -0.260 | -0.260 0.248 0.369 0.369 0.542 eV alpha orb energies kxyz= 27 1 1 0.000 0.000 -10.024 -0.699 -0.281 | -0.281 0.258 0.349 0.349 0.518 eV alpha orb energies kxyz= 26 1 1 0.000 0.000 -10.024 -0.711 -0.298 | -0.298 0.263 0.333 0.333 0.497 eV alpha orb energies kxyz= 25 1 1 0.000 0.000 -10.024 -0.721 -0.313 | -0.313 0.264 0.321 0.321 0.480 eV alpha orb energies kxyz= 24 1 1 0.000 0.000 -10.024 -0.729 -0.326 | -0.326 0.264 0.313 0.313 0.466 eV alpha orb energies kxyz= 23 1 1 0.000 0.000 -10.024 -0.735 -0.335 | -0.335 0.264 0.310 0.310 0.456 eV alpha orb energies kxyz= 22 1 1 0.000 0.000 -10.024 -0.739 -0.341 | -0.341 0.262 0.312 0.312 0.449 eV alpha orb energies kxyz= 21 1 1 0.000 0.000 -10.024 -0.741 -0.344 | -0.344 0.262 0.314 0.314 0.445 eV HOCO at k-point 4 kxyz= 37 1 1 -0.200927 eV LUCO at k-point 0 kxyz= 0 0 0 -9.361851 eV Indirect gap -9.160924 eV 在这里，我们举一个一维的碳链的例子(如上图所示)，来具体手把手的教你如何用Gaussian03计算固体的能带结构。同样的方法可以用来计算碳纳米管的能带结构等。以上部分包含的信息是每个K点上, HOMO和LUMO左右的十条轨道的能量。将他们连接起来(按照kxyz=0, 21, 22, 23,…,40的顺序)就是费米面附近的能带图了。 以上是网址为http://new.quantumchemistry.net/ ... lc/200603/1088.html的粒子。但我不知道输出是什么意思？第一行13个数是K=0时的能级值吗？ |
回复此楼» 猜你喜欢
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有236人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有65人回复
» 本主题相关商家推荐: (我也要在这里推广)
10