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

24小时热门版块排行榜    

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

赵奇一

新虫 (小有名气)

[求助] 关于QE计算结果的疑问 已有1人参与

最近尝试用QE算了算体系的DOS，印象中费米能附近的DOS值应该为0对吗？我的结果如下：
the Fermi energy is     4.7935 ev
但是在这个能量点周围
#  E (eV)   dos(E)     Int dos(E)
4.190  0.2737E+01  0.1668E+02
  4.200  0.2750E+01  0.1670E+02
  4.210  0.2773E+01  0.1673E+02
  4.220  0.2849E+01  0.1676E+02
  4.230  0.2941E+01  0.1679E+02
  4.240  0.2976E+01  0.1682E+02
  4.250  0.2924E+01  0.1685E+02
  4.260  0.2845E+01  0.1688E+02
  4.270  0.2773E+01  0.1690E+02
  4.280  0.2820E+01  0.1693E+02
  4.290  0.2889E+01  0.1696E+02
  4.300  0.2960E+01  0.1699E+02
  4.310  0.3026E+01  0.1702E+02
  4.320  0.3052E+01  0.1705E+02
  4.330  0.3044E+01  0.1708E+02
  4.340  0.3036E+01  0.1711E+02
  4.350  0.3033E+01  0.1714E+02
  4.360  0.3035E+01  0.1717E+02
  4.370  0.3030E+01  0.1720E+02
  4.380  0.2997E+01  0.1723E+02
  4.390  0.2953E+01  0.1726E+02
  4.400  0.2914E+01  0.1729E+02
  4.410  0.2894E+01  0.1732E+02
  4.420  0.2876E+01  0.1735E+02
  4.430  0.2863E+01  0.1738E+02
  4.440  0.2869E+01  0.1741E+02
  4.450  0.2892E+01  0.1744E+02
  4.460  0.2980E+01  0.1747E+02
  4.470  0.2993E+01  0.1750E+02
  4.480  0.2958E+01  0.1753E+02
  4.490  0.2925E+01  0.1755E+02
  4.500  0.2898E+01  0.1758E+02
  4.510  0.2873E+01  0.1761E+02
  4.520  0.2856E+01  0.1764E+02
  4.530  0.2843E+01  0.1767E+02
  4.540  0.2836E+01  0.1770E+02
  4.550  0.2832E+01  0.1773E+02
  4.560  0.2845E+01  0.1775E+02
  4.570  0.2893E+01  0.1778E+02
  4.580  0.2964E+01  0.1781E+02
  4.590  0.3053E+01  0.1784E+02
  4.600  0.3151E+01  0.1787E+02
  4.610  0.3248E+01  0.1791E+02
  4.620  0.3276E+01  0.1794E+02
  4.630  0.3354E+01  0.1797E+02
  4.640  0.3447E+01  0.1801E+02
  4.650  0.3540E+01  0.1804E+02
  4.660  0.3662E+01  0.1808E+02
  4.670  0.3840E+01  0.1812E+02
  4.680  0.3909E+01  0.1816E+02
  4.690  0.3933E+01  0.1820E+02
  4.700  0.3904E+01  0.1824E+02
  4.710  0.3856E+01  0.1827E+02
  4.720  0.3808E+01  0.1831E+02
  4.730  0.3811E+01  0.1835E+02
  4.740  0.3821E+01  0.1839E+02
  4.750  0.3801E+01  0.1843E+02
  4.760  0.3578E+01  0.1846E+02
  4.770  0.3457E+01  0.1850E+02
  4.780  0.3384E+01  0.1853E+02
  4.790  0.3387E+01  0.1856E+02
  4.800  0.3566E+01  0.1860E+02
  4.810  0.3758E+01  0.1864E+02
  4.820  0.3972E+01  0.1868E+02
  4.830  0.4310E+01  0.1872E+02
  4.840  0.4679E+01  0.1877E+02
  4.850  0.4772E+01  0.1882E+02
  4.860  0.4730E+01  0.1886E+02
  4.870  0.4906E+01  0.1891E+02
  4.880  0.4997E+01  0.1896E+02
  4.890  0.5062E+01  0.1901E+02
  4.900  0.5113E+01  0.1906E+02
  4.910  0.5163E+01  0.1912E+02
  4.920  0.5212E+01  0.1917E+02
  4.930  0.5242E+01  0.1922E+02
  4.940  0.5054E+01  0.1927E+02
  4.950  0.4686E+01  0.1932E+02
  4.960  0.4392E+01  0.1936E+02
  4.970  0.4231E+01  0.1940E+02
  4.980  0.4085E+01  0.1944E+02
  4.990  0.3984E+01  0.1948E+02
  5.000  0.4289E+01  0.1953E+02
  5.010  0.4243E+01  0.1957E+02
  5.020  0.4191E+01  0.1961E+02
  5.030  0.4144E+01  0.1965E+02
  5.040  0.4112E+01  0.1969E+02
  5.050  0.4245E+01  0.1974E+02
  5.060  0.4637E+01  0.1978E+02
  5.070  0.4490E+01  0.1983E+02
  5.080  0.4272E+01  0.1987E+02
  5.090  0.4363E+01  0.1991E+02
  5.100  0.4432E+01  0.1996E+02
  5.110  0.4571E+01  0.2000E+02
  5.120  0.5174E+01  0.2006E+02
  5.130  0.5474E+01  0.2011E+02
  5.140  0.5755E+01  0.2017E+02
  5.150  0.5945E+01  0.2023E+02
  5.160  0.5695E+01  0.2028E+02
  5.170  0.5393E+01  0.2034E+02
  5.180  0.5194E+01  0.2039E+02
  5.190  0.5273E+01  0.2044E+02
  5.200  0.5383E+01  0.2050E+02
  5.210  0.5338E+01  0.2055E+02
  5.220  0.5289E+01  0.2060E+02
  5.230  0.5268E+01  0.2066E+02
  5.240  0.5237E+01  0.2071E+02
  5.250  0.5129E+01  0.2076E+02
  5.260  0.5042E+01  0.2081E+02
  5.270  0.5077E+01  0.2086E+02
  5.280  0.5177E+01  0.2091E+02
  5.290  0.5229E+01  0.2096E+02
  5.300  0.5211E+01  0.2102E+02
  5.310  0.5176E+01  0.2107E+02
  5.320  0.5127E+01  0.2112E+02
  5.330  0.5073E+01  0.2117E+02
  5.340  0.5020E+01  0.2122E+02
  5.350  0.4982E+01  0.2127E+02
  5.360  0.4952E+01  0.2132E+02
  5.370  0.4887E+01  0.2137E+02
  5.380  0.4745E+01  0.2142E+02
  5.390  0.4608E+01  0.2146E+02
  5.400  0.4535E+01  0.2151E+02
  5.410  0.4475E+01  0.2155E+02
  5.420  0.4490E+01  0.2160E+02
  5.430  0.4582E+01  0.2164E+02
  5.440  0.4554E+01  0.2169E+02
  5.450  0.4533E+01  0.2173E+02
这是什么原因？能量间隔取得太大？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
天使可以飞起来是因为他们把自己看得很轻
1楼 2015-07-15 10:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-17 11:52:55
赵奇一: 金币+15, ★★★很有帮助 2015-07-18 11:03:56
赵奇一: 金币+5, ★★★很有帮助 2015-07-18 16:40:35
引用回帖:
9楼: Originally posted by 赵奇一 at 2015-07-16 19:44:26
&control
calculation='scf'
restart_mode='from_scratch',
pseudo_dir = './',
outdir='./'
prefix='C'
tstress = .true.
tprnfor = .true.
/
&system
ibrav = 14, celldm(1) =9.498 ...

你的坐标是从MS中抄过来的吗？
这个单位应该是晶体内坐标单位，默认是alat单位不是内坐标，所以结构错了。
ATOMIC_POSITIONS 这行应该加上单位 改成
ATOMIC_POSITIONS crystal
赞一下 回复此楼
10楼2015-07-16 21:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

KalaShayminS

铁杆木虫 (著名写手)
这是个金属体系，Fermi面附近肯定不是0，并没有什么问题。
赞一下 回复此楼
2楼2015-07-15 18:29:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)
引用回帖:
2楼: Originally posted by KalaShayminS at 2015-07-15 18:29:58
这是个金属体系，Fermi面附近肯定不是0，并没有什么问题。

我建立的模型是石墨烯...最后画出来的图在下边（图是我试着画的，很粗略，抱歉）。横坐标是能量，纵坐标为DOS值，感觉跟论坛里另一张图差别比较大，这样看来是不是我建的模型有问题？
关于QE计算结果的疑问
石墨烯.png
赞一下 回复此楼
天使可以飞起来是因为他们把自己看得很轻
3楼2015-07-15 21:24:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)
这是论坛里另一张图，是AesNaegea画的
关于QE计算结果的疑问-1
石墨烯vasp.jpg
赞一下 回复此楼
天使可以飞起来是因为他们把自己看得很轻
4楼2015-07-15 21:25:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +30 AI8RaGaPaSCR 2026-08-07 32/1600 2026-08-13 21:45 by Hd5o188FDaoB
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 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
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
关闭
请填处理意见
关闭