应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】自洽计算的问题
51/1返回列表
查看: 1016  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

风轻云淡7

金虫 (小有名气)

[交流] 【求助】自洽计算的问题 已有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呢？不胜感激……呜呜……谢谢啊……
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-08-30 10:39:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风轻云淡7

金虫 (小有名气)
引用回帖:
Originally posted by valenhou001 at 2010-08-30 19:20:27:
需补一点基本的知识：
1）你计算的原胞里有2个C原子，每个C原子是4个价电子，所以总共8个价电子，电子从低填到高，所以电子填到第个能带，即每个k点的第个本征值。
“highest occupied, lowest unoccupied leve ...

恩恩……谢谢了啊……受益匪浅啊……嘿嘿就是本来我要设置为六角晶格的，但是结果是不对的，有十几个电子伏的能隙，所以请教了师兄，师兄要我设置为12试试，我就试了一下，结果还不错的……
赞一下 回复此楼
4楼2010-08-31 09:31:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

valenhou001

至尊木虫 (职业作家)
★ ★ ★ ★
小木虫(金币+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人) 回复此楼
2楼2010-08-30 19:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢补充 2010-08-31 08:05:46
xiaohunhun(金币+1):xiexie 2010-08-31 18:59:15
补充一点，石墨是半金属体系，你应该设置smearing等参数。
赞一下(3人) 回复此楼
3楼2010-08-30 20:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风轻云淡7

金虫 (小有名气)
引用回帖:
Originally posted by goldenfisher at 2010-08-30 20:59:43:
补充一点，石墨是半金属体系，你应该设置smearing等参数。

嗯嗯……好的，我改改的……谢谢了啊……
赞一下 回复此楼
5楼2010-08-31 09:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭