| 查看: 2449 | 回復(fù): 17 | ||
趙奇一新蟲 (小有名氣)
|
[求助]
關(guān)于QE計(jì)算結(jié)果的疑問 已有1人參與
|
|
最近嘗試用QE算了算體系的DOS,印象中費(fèi)米能附近的DOS值應(yīng)該為0對嗎?我的結(jié)果如下: the Fermi energy is 4.7935 ev 但是在這個(gè)能量點(diǎn)周圍 # 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 這是什么原因?能量間隔取得太大? |

鐵桿木蟲 (著名寫手)
新蟲 (小有名氣)

新蟲 (小有名氣)

鐵桿木蟲 (著名寫手)
鐵桿木蟲 (正式寫手)

新蟲 (小有名氣)
|
在您的指點(diǎn)下,這兩周一直嘗試通過QE進(jìn)行求解模型DOS值。我主要通過MS進(jìn)行建模,然后將原子坐標(biāo),a,b,c以及α,β,γ的值代換到QE輸入文件中去。計(jì)算過程沒什么問題,但是結(jié)果差別很大。除了這個(gè)模型外還有單層MoS2,我直接通過MS進(jìn)行計(jì)算,得到的DOS和一篇文獻(xiàn)用VASP結(jié)果對照。在相同的能量區(qū)間里,走勢沒有什么區(qū)別,無非就是相同能量點(diǎn)一個(gè)DOS值大,一個(gè)DOS值小。聽別人說這是選擇K點(diǎn)不同造成的,兩個(gè)都對(這種說法正確嗎)。 但是我用QE進(jìn)行計(jì)算的結(jié)果與文獻(xiàn)和MS得到的結(jié)果差別很大。很郁悶。是不是MS建立的模型在向QE輸入文件轉(zhuǎn)換的過程中有什么細(xì)節(jié)方面要注意而我不知道?還處于初學(xué)階段,希望您能不吝賜教,感激不盡... |

鐵桿木蟲 (著名寫手)
新蟲 (小有名氣)
|
&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 謝謝你... |

鐵桿木蟲 (著名寫手)
|
你的坐標(biāo)是從MS中抄過來的嗎? 這個(gè)單位應(yīng)該是晶體內(nèi)坐標(biāo)單位,默認(rèn)是alat單位不是內(nèi)坐標(biāo),所以結(jié)構(gòu)錯(cuò)了。 ATOMIC_POSITIONS 這行應(yīng)該加上單位 改成 ATOMIC_POSITIONS crystal |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|