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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 关于QE计算结果的疑问
181/2返回列表12下一页
查看: 2327  |  回复: 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

铁杆木虫 (著名写手)
这是个金属体系，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的回帖

KalaShayminS

铁杆木虫 (著名写手)
引用回帖:
3楼: Originally posted by 赵奇一 at 2015-07-15 21:24:14
我建立的模型是石墨烯...最后画出来的图在下边（图是我试着画的，很粗略，抱歉）。横坐标是能量，纵坐标为DOS值，感觉跟论坛里另一张图差别比较大，这样看来是不是我建的模型有问题？

石墨烯.png
...

这个确实不对，先检查检查模型吧
赞一下 回复此楼
5楼2015-07-16 09:41:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朱十三

铁杆木虫 (正式写手)
嗯，知道了，谢谢你的解答
赞一下 回复此楼
话说，以前穷人是不取名的，以出生日期命名，如此算来本来该叫廿三。无奈竟忘记还有廿字，遂减去十，故称十三
6楼2015-07-16 11:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)
引用回帖:
5楼: Originally posted by KalaShayminS at 2015-07-16 09:41:00
这个确实不对，先检查检查模型吧...

在您的指点下，这两周一直尝试通过QE进行求解模型DOS值。我主要通过MS进行建模，然后将原子坐标，a,b,c以及α,β,γ的值代换到QE输入文件中去。计算过程没什么问题，但是结果差别很大。除了这个模型外还有单层MoS2，我直接通过MS进行计算，得到的DOS和一篇文献用VASP结果对照。在相同的能量区间里，走势没有什么区别，无非就是相同能量点一个DOS值大，一个DOS值小。听别人说这是选择K点不同造成的，两个都对（这种说法正确吗）。
但是我用QE进行计算的结果与文献和MS得到的结果差别很大。很郁闷。是不是MS建立的模型在向QE输入文件转换的过程中有什么细节方面要注意而我不知道？还处于初学阶段，希望您能不吝赐教，感激不尽...
赞一下 回复此楼
天使可以飞起来是因为他们把自己看得很轻
7楼2015-07-16 18:16:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
引用回帖:
7楼: Originally posted by 赵奇一 at 2015-07-16 18:16:57
在您的指点下，这两周一直尝试通过QE进行求解模型DOS值。我主要通过MS进行建模，然后将原子坐标，a,b,c以及α,β,γ的值代换到QE输入文件中去。计算过程没什么问题，但是结果差别很大。除了这个模型外还有单层MoS2 ...

贴一个石墨烯的输入文件吧。
赞一下 回复此楼
8楼2015-07-16 19:40:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)
引用回帖:
8楼: Originally posted by KalaShayminS at 2015-07-16 19:40:52
贴一个石墨烯的输入文件吧。...

&control
calculation='scf'
restart_mode='from_scratch',
pseudo_dir = './',
outdir='./'
prefix='C'
tstress = .true.
tprnfor = .true.
/
&system
ibrav = 14, celldm(1) =9.498009, celldm(2) =1.0, celldm(3) =1.9383, celldm(4) =0.0, celldm(5) =0.0, celldm(6) =-0.5, nat= 18, ntyp= 1,
ecutwfc = 25.0, ecutrho = 300.0
occupations='smearing', smearing='gaussian', degauss=0.02
/
&electrons
diagonalization='david'
conv_thr = 1.0e-8
mixing_beta = 0.7
/
ATOMIC_SPECIES
C  12.0  C.pz-rrkjus.UPF
ATOMIC_POSITIONS
  C   0.0000000000000000   0.0000000000000000   0.2499999999999999
  C   0.1111111111111111   0.2222222222222220   0.2500000000000000
  C   0.3333333333333333   0.0000000000000000   0.2500000000000000
  C   0.4444444444444444   0.2222222222222220   0.2500000000000000
  C   0.6666666666666666   0.0000000000000000   0.2500000000000000
  C   0.7777777777777778   0.2222222222222221   0.2500000000000000
  C   0.0000000000000000   0.3333333333333333   0.2500000000000000
  C   0.1111111111111111   0.5555555555555554   0.2500000000000000
  C   0.3333333333333333   0.3333333333333333   0.2500000000000000
  C   0.4444444444444444   0.5555555555555555   0.2500000000000000
  C   0.6666666666666666   0.3333333333333333   0.2500000000000000
  C   0.7777777777777778   0.5555555555555555   0.2499999999999999
  C   0.0000000000000000   0.6666666666666666   0.2500000000000000
  C   0.1111111111111111   0.8888888888888887   0.2500000000000000
  C   0.3333333333333333   0.6666666666666666   0.2500000000000000
  C   0.4444444444444444   0.8888888888888887   0.2499999999999999
  C   0.6666666666666666   0.6666666666666666   0.2499999999999999
  C   0.7777777777777778   0.8888888888888887   0.2499999999999999
K_POINTS (automatic)
8 8 8 0 0 0
谢谢你...
赞一下 回复此楼
天使可以飞起来是因为他们把自己看得很轻
9楼2015-07-16 19:44:26
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 赵奇一 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭