| 查看: 271 | 回复: 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时的能级值吗？ |
回复此楼» 猜你喜欢
材料科学与工程 或 材料与化工（专硕）调剂 要求：初试考408
已经有0人回复
-
298分 070300求调剂
已经有15人回复
-
物理化学论文润色/翻译怎么收费?
已经有62人回复
-
湖南大学张世国团队招26年第二批博士—离子液体、高分子与凝胶材料
已经有0人回复
-
湖南农业大学能源材料创新团队招收2025化学专业研究生（学硕）调剂
已经有0人回复
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
已经有9人回复
-
调剂信息
已经有0人回复
-
河北农大生物与医药专业招收调剂，名额充足
已经有0人回复
-
齐鲁工业大学 国家重点实验室-新能源、锂离子电池、材料、化学方向
已经有0人回复
-
杰青/长江团队招收调剂
已经有0人回复
-
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
