| 查看: 1013 | 回复: 4 | |||
[交流]
【求助】自洽计算的问题 已有3人参与
|
|
&control calculation='scf' restart_mode='from_scratch', prefix='c' pseudo_dir = '/home/espresso-4.0.4/pseudo/', outdir='./', / &system ibrav= 12, celldm(1)=4.705417712,celldm(2)=1.0000000, celldm(3)=7.589260513,celldm(4)=0.500000000, nat= 2, ntyp= 1, ecutrho = 500.D0, ecutwfc =40.D0,nbnd=12, / &electrons diagonalization = 'cg' electron_maxstep = 200 conv_thr = 1.0d-5 mixing_beta = 0.70 / ATOMIC_SPECIES C 12.000000 C.pz-vbc.UPF ATOMIC_POSITIONS {crystal} C 0.000000000 0.000000000 0.000000000 C 0.333333333 0.333333333 0.000000000 K_POINTS {automatic} 24 24 1 0 0 0 大家好，这是我计算石墨的输入文件， k = 0.4167-0.7217 0.0000 ( 2908 PWs) bands (ev): -16.2608 -15.3256 -10.8233 -4.3372 -1.1748 6.9231 8.7710 9.0284 9.1780 10.3615 10.6452 10.6547 k = 0.4167-0.6736 0.0000 ( 2903 PWs) bands (ev): -16.1926 -14.9544 -11.5617 -4.0872 -1.2757 8.0781 8.3037 8.8358 9.9380 10.3714 10.4093 11.4086 1032,0-1 86% 8.2287 8.7241 9.5418 10.1591 k =-0.5000 0.5774 0.0000 ( 2921 PWs) bands (ev): -17.6296 -13.0474 -11.5440 -5.3178 0.2057 6.6358 7.3617 7.9731 8.6476 9.2495 9.5167 9.7175 highest occupied, lowest unoccupied level (ev): -2.5882 -2.5878 这是我输出文件的一部分，我想问一下，几个问题，谢谢……第一个，作为导体的输出文件中是不是不应该用最高的占据态和最低的为占据态？第二，能量的输出中应该负值的数值比较多才对，当时我的正数比较多，请问问题出在哪里？第三，我们知道石墨是单层的六角形了，我这里设置的为单斜精细，那我的晶格常数a应该是实际的c-c之间的距离呢，还是根号3在乘以a呢？不胜感激……呜呜……谢谢啊…… |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有68人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp磁性计算问题
已经有8人回复
- vasp算法的一些基础问题
已经有10人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- 静态自洽计算遇到的问题
已经有5人回复
- 静态自洽计算能量时EDIFFG的问题
已经有14人回复
- 【求助】我在用ATK做自洽计算的时候，出现如下的错误：
已经有3人回复
- 【求助】求助MS计算式任务进行顺序的问题
已经有7人回复
- 【求助】结构优化之后的静态自洽计算
已经有12人回复
- 【讨论】自洽计算与非自洽计算的含义、所做的工作
已经有10人回复
- 【求助】Hse06计算时KPOINTS出错
已经有30人回复
- 【求助】自洽计算与非自洽计算的区别
已经有4人回复
- 【求助成功】计算能带结构时出现错误，请教
已经有19人回复
- 【求助】静态自洽计算读WAVECAR有用吗？
已经有9人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【求助】态密度计算究竟是自洽计算，还是非自洽计算？
已经有24人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+3):谢谢侯老师指点 2010-08-31 08:05:36
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+3):谢谢侯老师指点 2010-08-31 08:05:36
|
需补一点基本的知识： 1）你计算的原胞里有2个C原子，每个C原子是4个价电子，所以总共8个价电子，电子从低填到高，所以电子填到第个能带，即每个k点的第个本征值。 “highest occupied, lowest unoccupied level (ev): -2.5882 -2.5878”给出了最高占据态的能量值，以及最低未占据态的能量值，即价带顶和导带底的值。 你总共计算了12个能带（由nbnd=12设置来确定的)。 本征值的绝对值是没有意义的。 你所提到的正数 多，是正常，没有问题。 那你为什么要设置“ ibrav= 12,”, 晶格常数a与C-C的键长关系来，你应按你所设置 ibrav= 12,所对应的晶系来确定。你可以用xcrysden --pwi 你的pwscf输入文件， 通过xcrysden可视化你的晶体结构，然后通过xcrysden上的工具测出原子的间距。 |
2楼2010-08-30 19:20:27
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
3楼2010-08-30 20:59:43
就是本来我要设置为六角晶格的，但是结果是不对的，有十几个电子伏的能隙，所以请教了师兄，师兄要我设置为12试试，我就试了一下，结果还不错的…… 4楼2010-08-31 09:31:55
5楼2010-08-31 09:32:21
