| 查看: 305 | 回复: 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时的能级值吗？ |
回复此楼» 猜你喜欢
安徽大学物质科学与信息技术研究院罗其全教授组2026年化学博士招生（计算催化组）
已经有1人回复
-
现招2026北京航空航天大学联培博士1名（北航学位）-材料AI方向
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有82人回复
-
普瓦捷环境与材料化学研究所 (IC2MP)博后招聘
已经有0人回复
-
香港科技大学(广州)方超课题组招聘全奖博士
已经有0人回复
-
宁夏大学化学工程与技术专业招收博士研究生
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
南京大学软物质计算模拟方向硕博连读/直博生招生
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
100