[交流] 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化

» 收錄本帖的淘帖專輯推薦

VASP	VASP軟件從入門到高深

» 猜你喜歡

268求調(diào)劑已經(jīng)有3人回復
302求調(diào)劑已經(jīng)有5人回復
311求調(diào)劑已經(jīng)有10人回復
被我言中：新模板不強調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復
303求調(diào)劑已經(jīng)有5人回復
考研求調(diào)劑已經(jīng)有4人回復
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復
能源材料化學課題組招收碩士研究生8-10名已經(jīng)有6人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

關(guān)于VASP一些瑣事的求助（K點，吸附優(yōu)化等）已經(jīng)有12人回復
晶體結(jié)構(gòu)優(yōu)化和畫能帶圖的方法是否保持一致？？？已經(jīng)有12人回復
單軸施加應(yīng)變結(jié)構(gòu)優(yōu)化的方法有哪些？已經(jīng)有7人回復
求助vasp_晶格參數(shù)掃描已經(jīng)有5人回復
vasp計算DOS,以及如何判斷參數(shù)的設(shè)置是合理的，達到要求的精度已經(jīng)有7人回復
vasp計算中，如何選擇正確的K點已經(jīng)有7人回復
【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化已經(jīng)有142人回復
【求助】晶格常數(shù)大了以后 VASP的計算值準確嗎已經(jīng)有5人回復
【求助】vasp優(yōu)化晶格參數(shù)a 已經(jīng)有13人回復
【求助】簡單問題----關(guān)于vasp晶格優(yōu)化已經(jīng)有7人回復
【求助】侯老師的vasp入門的晶格優(yōu)化部分沒看懂，問一下已經(jīng)有25人回復
【求助】VASP結(jié)構(gòu)優(yōu)化已經(jīng)有7人回復
【求助成功】VASP中，如何固定磁矩做結(jié)構(gòu)優(yōu)化已經(jīng)有14人回復
【原創(chuàng)】晶格優(yōu)化中ISIF取值的探究（六方MgO為例）已經(jīng)有9人回復
【討論】vasp結(jié)構(gòu)優(yōu)化之前的測試已經(jīng)有11人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-12-07 09:44:41

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★ ★
gzqdyouxia(金幣+2):謝謝指點 2010-12-10 09:20:31

引用回帖:

Originally posted by bluegirlll at 2010-12-08 16:21:11:
你好，我看了，并且按照他做的，可是怎么出不來結(jié)果呢，編輯腳本有沒有需要注意的，另外那個如果有兩個元素的材料就要進行復雜的優(yōu)化嗎？還有他的那個腳本為什么要刪除rm WAVECAR。請高人指點

我也是個菜鳥，對優(yōu)化不清楚，但是計算好像都得優(yōu)化啊。還有那個rm WAVECAR一般是要有的，如果你的計算目錄下沒有WAVECAR時，你也可以不要這個命令行。如果有wavecar時，istart=0將是無效的（好像是，你看看手冊，我有點既不清楚了）。

10樓2010-12-08 16:49:21

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

first_yg

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 463.5
帖子: 132
在線: 81.8小時
蟲號: 374955

★
sunyang1988(金幣+1):謝謝指點 2010-12-07 09:55:44
bluegirlll(金幣+1): 2010-12-08 16:18:50

侯博士的中文說明里就有，你可以查查看

2樓2010-12-07 09:51:38

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ustbwgs

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 1449.6
帖子: 720
在線: 165.5小時
蟲號: 1141993

★ ★
zzy870720z(金幣+2):謝謝提示 2010-12-07 14:18:47
bluegirlll(金幣+1): 2010-12-08 09:27:38

先根據(jù)文獻或?qū)嶒瀰?shù)寫個差不多的，然后再在這個取值附近的范圍內(nèi)，寫個腳本計算單點能，做能量收斂測試。

3樓2010-12-07 14:11:34

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluegirlll

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 105.9
帖子: 57
在線: 20.6小時
蟲號: 1046082

謝謝您啊，我用的是這個個腳本：

#! /bin/bash
BIN=/home/jhm/VASP/vasp
for i in 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 ; do
cat >POSCAR < fcc:
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 2
cartesian
0 0 0
0.25 0.25 0.25
0.75 0.75 0.75
!
echo "a= $i" ; $BIN
E=·tail -l OSZICAR·; echo $i $E >>SUMMARY.fcc
done
計算得到的SUMMARY.fcc文件的內(nèi)容為：

3.5 tail -l OSZICAR
3.6 tail -l OSZICAR
3.7 tail -l OSZICAR
3.8 tail -l OSZICAR
3.9 tail -l OSZICAR
4.0 tail -l OSZICAR
4.1 tail -l OSZICAR
4.2 tail -l OSZICAR
4.3 tail -l OSZICAR

是我腳本那個地方錯了嗎？請高人指點

引用回帖:

Originally posted by ustbwgs at 2010-12-07 14:11:34:
先根據(jù)文獻或?qū)嶒瀰?shù)寫個差不多的，然后再在這個取值附近的范圍內(nèi)，寫個腳本計算單點能，做能量收斂測試。

4樓2010-12-08 09:34:25

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ustbwgs

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 1449.6
帖子: 720
在線: 165.5小時
蟲號: 1141993

★ ★
ellsaking(金幣+2):謝謝指導 2010-12-08 10:03:15
bluegirlll(金幣+2): 2010-12-08 11:02:34

引用回帖:

Originally posted by bluegirlll at 2010-12-08 09:34:25:
謝謝您啊，我用的是這個個腳本：

#! /bin/bash
BIN=/home/jhm/VASP/vasp
for i in 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 ; do
cat >POSCAR < fcc:
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0 ...

E=·tail -l OSZICAR·; echo $i $E >>SUMMARY.fcc
改為E=“tail -l OSZICAR”
echo $i $E >>SUMMARY.fcc
也就是把倒引改為雙引，你看對不對，我一般也不用這個命令，可以查energy without entropy 或者TOTEN，效果一樣啊

[ Last edited by ustbwgs on 2010-12-8 at 10:03 ]

5樓2010-12-08 10:00:26

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★
zzy870720z(金幣+1):謝謝提示 2010-12-08 12:01:52

引用回帖:

把引號改成··，就是tab鍵上的那個鍵，試試吧

6樓2010-12-08 10:36:06

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluegirlll

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 105.9
帖子: 57
在線: 20.6小時
蟲號: 1046082

我已經(jīng)按您說的改了，可是出來這樣的結(jié)果，它去了OSIZCAR文件中的好多數(shù)據(jù)
SUMMARY.fcc內(nèi)容如下：

3.5 DAV: 9 0.616371635786E+02 -0.21638E+01 -0.36614E+00 887 0.211E+01 0.379E+00 DAV: 10 0.598839270162E+02 -0.17532E+01 -0.22497E+00 946 0.164E+01 0.135E+00 DAV: 11 0.599095645060E+02 0.25637E-01 -0.14704E-01 931 0.481E+00 0.625E-01 DAV: 12 0.599206379143E+02 0.11073E-01 -0.57937E-02 824 0.307E+00 0.358E-01 DAV: 13 0.599130268782E+02 -0.76110E-02 -0.27145E-02 953 0.152E+00 0.107E-01 DAV: 14 0.599127542034E+02 -0.27267E-03 -0.25634E-03 1209 0.508E-01 0.870E-02 DAV: 15 0.599126009216E+02 -0.15328E-03 -0.17119E-04 615 0.166E-01 0.464E-02 DAV: 16 0.599123843746E+02 -0.21655E-03 -0.15291E-04 685 0.131E-01 0.242E-02 DAV: 17 0.599123887861E+02 0.44115E-05 -0.17822E-05 398 0.428E-02 1 F= 0.59912389E+02 E0= 0.59912389E+02 d E =0.000000E+00 mag= -0.0009
3.6 DAV: 8 0.486136355673E+02 -0.35626E+01 -0.34687E+00 890 0.236E+01 0.522E+00 DAV: 9 0.469129254896E+02 -0.17007E+01 -0.18613E+00 940 0.131E+01 0.343E+00 DAV: 10 0.461778722637E+02 -0.73505E+00 -0.13600E+00 876 0.131E+01 0.144E+00 DAV: 11 0.462037274218E+02 0.25855E-01 -0.10867E-01 882 0.402E+00 0.832E-01 DAV: 12 0.461914507126E+02 -0.12277E-01 -0.69380E-02 829 0.287E+00 0.323E-01 DAV: 13 0.461915837635E+02 0.13305E-03 -0.55215E-03 1206 0.675E-01 0.121E-01 DAV: 14 0.461917513859E+02 0.16762E-03 -0.15565E-03 1084 0.526E-01 0.796E-02 DAV: 15 0.461913353655E+02 -0.41602E-03 -0.36541E-04 788 0.163E-01 0.584E-02 DAV: 16 0.461912802755E+02 -0.55090E-04 -0.28427E-04 663 0.170E-01 1 F= 0.46191280E+02 E0= 0.46191280E+02 d E =0.000000E+00 mag= 0.0122
3.7 DAV: 21 0.345758323297E+02 -0.93737E-03 -0.94467E-04 929 0.306E-01 0.402E-01 DAV: 22 0.345795760961E+02 0.37438E-02 -0.48287E-03 922 0.658E-01 0.264E-01 DAV: 23 0.345802496016E+02 0.67351E-03 -0.12725E-03 968 0.325E-01 0.202E-01 DAV: 24 0.345826284269E+02 0.23788E-02 -0.54723E-04 796 0.218E-01 0.161E-01 DAV: 25 0.345871012547E+02 0.44728E-02 -0.30116E-03 938 0.545E-01 0.130E-01 DAV: 26 0.345872092251E+02 0.10797E-03 -0.30063E-04 728 0.187E-01 0.972E-02 DAV: 27 0.345858613196E+02 -0.13479E-02 -0.55383E-04 816 0.274E-01 0.609E-02 DAV: 28 0.345862324715E+02 0.37115E-03 -0.52166E-04 748 0.350E-01 0.271E-02 DAV: 29 0.345861649687E+02 -0.67503E-04 -0.13802E-05 395 0.461E-02 1 F= 0.34586165E+02 E0= 0.34586165E+02 d E =0.000000E+00 mag= 0.42683.8 DAV: 7 0.279114265739E+02 -0.16376E+01 -0.56337E+00 1022 0.204E+01 0.665E+00 DAV: 8 0.254025058027E+02 -0.25089E+01 -0.34546E+00 870 0.240E+01 0.374E+00 DAV: 9 0.248851128671E+02 -0.51739E+00 -0.47666E-01 992 0.551E+00 0.273E+00 DAV: 10 0.244305976128E+02 -0.45452E+00 -0.56551E-01 893 0.739E+00 0.128E+00 DAV: 11 0.243774407325E+02 -0.53157E-01 -0.99253E-02 1009 0.289E+00 0.572E-01 DAV: 12 0.243763728670E+02 -0.10679E-02 -0.99008E-03 878 0.101E+00 0.177E-01 DAV: 13 0.243764253820E+02 0.52515E-04 -0.16325E-03 1012 0.268E-01 0.420E-02 DAV: 14 0.243762900673E+02 -0.13531E-03 -0.12378E-04 502 0.108E-01 0.291E-02 DAV: 15 0.243762466643E+02 -0.43403E-04 -0.59665E-05 425 0.931E-02 1 F= 0.24376247E+02 E0= 0.24376247E+02 d E =0.000000E+00 mag= 0.00023.9 DAV: 6 0.186939152777E+02 -0.20100E+00 -0.56562E+00 1138 0.238E+01 0.105E+01 DAV: 7 0.175924468941E+02 -0.11015E+01 -0.48486E+00 1076 0.181E+01 0.604E+00 DAV: 8 0.159033008957E+02 -0.16891E+01 -0.34831E+00 842 0.257E+01 0.342E+00 DAV: 9 0.157252883319E+02 -0.17801E+00 -0.39195E-01 1270 0.648E+00 0.242E+00 DAV: 10 0.153946901067E+02 -0.33060E+00 -0.28033E-01 986 0.450E+00 0.127E+00 DAV: 11 0.153219436031E+02 -0.72747E-01 -0.12657E-01 974 0.338E+00 0.493E-01 DAV: 12 0.153217217700E+02 -0.22183E-03 -0.12968E-02 1019 0.155E+00 0.230E-01 DAV: 13 0.153228266672E+02 0.11049E-02 -0.25340E-03 1071 0.660E-01 0.796E-02 DAV: 14 0.153229062451E+02 0.79578E-04 -0.32392E-04 758 0.153E-01 1 F= 0.15322906E+02 E0= 0.15322906E+02 d E =0.000000E+00 mag= 0.00034.0 DAV: 6 0.931139208363E+01 0.56404E+00 -0.56782E+00 1180 0.241E+01 0.957E+00 DAV: 7 0.878331383251E+01 -0.52808E+00 -0.41165E+00 1150 0.170E+01 0.519E+00 DAV: 8 0.775427396386E+01 -0.10290E+01 -0.20414E+00 854 0.190E+01 0.279E+00 DAV: 9 0.755153875548E+01 -0.20274E+00 -0.26699E-01 1247 0.457E+00 0.199E+00 DAV: 10 0.734397362079E+01 -0.20757E+00 -0.13448E-01 918 0.347E+00 0.112E+00 DAV: 11 0.728475733198E+01 -0.59216E-01 -0.81750E-02 1000 0.288E+00 0.398E-01 DAV: 12 0.728444227142E+01 -0.31506E-03 -0.36268E-03 925 0.636E-01 0.201E-01 DAV: 13 0.728455095831E+01 0.10869E-03 -0.14001E-03 1122 0.340E-01 0.352E-02 DAV: 14 0.728458422560E+01 0.33267E-04 -0.81814E-05 452 0.917E-02 1 F= 0.72845842E+01 E0= 0.72845842E+01 d E =0.000000E+00 mag= -0.00014.1 DAV: 5 0.167530005868E+00 -0.14240E-01 -0.14230E-01 988 0.256E+00 0.177E+01 DAV: 6 0.125671077861E+01 0.10892E+01 -0.57949E+00 1187 0.243E+01 0.874E+00 DAV: 7 0.11036......................................................
請問這是怎么回事��？非常感謝！

引用回帖:

Originally posted by zhang_jaj at 2010-12-08 10:36:06:

把引號改成··，就是tab鍵上的那個鍵，試試吧

7樓2010-12-08 11:05:12

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★ ★
sunyang1988(金幣+2):謝謝指點 2010-12-08 14:32:48
bluegirlll(金幣+1): 2010-12-10 09:07:08

引用回帖:

Originally posted by bluegirlll at 2010-12-08 11:05:12:
我已經(jīng)按您說的改了，可是出來這樣的結(jié)果，它去了OSIZCAR文件中的好多數(shù)據(jù)
SUMMARY.fcc內(nèi)容如下：

3.5 DAV: 9 0.616371635786E+02 -0.21638E+01 -0.36614E+00 887 0.211E+01 0.379E+00 DAV: 10 0.59883927016 ...

看你這個好像把OSZICAR文件里的所有數(shù)據(jù)都讀取了，tail -1 OSZICAR是讓取最后一行的呀，你為什么不試試用OUTCAR里的數(shù)據(jù)，侯老師的資料你有
E=·grep 'TOTEN' OUTCAR | tail -1 awk '{printf "%12.6f \n" , $5 }'·
注意最外層的引號是那個tab上的鍵，還有你開始的輸出錯誤，不會是因為把tail -1打成tail -l了吧。。

8樓2010-12-08 13:28:00

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluegirlll

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 105.9
帖子: 57
在線: 20.6小時
蟲號: 1046082

你好，我看了，并且按照他做的，可是怎么出不來結(jié)果呢，編輯腳本有沒有需要注意的，另外那個如果有兩個元素的材料就要進行復雜的優(yōu)化嗎？還有他的那個腳本為什么要刪除rm WAVECAR。請高人指點

引用回帖:

Originally posted by first_yg at 2010-12-07 09:51:38:
侯博士的中文說明里就有，你可以查查看

9樓2010-12-08 16:21:11

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluegirlll

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 105.9
帖子: 57
在線: 20.6小時
蟲號: 1046082

您好，真是非常感謝，結(jié)果已經(jīng)出來了。可是總能怎么出來正負不一直的結(jié)果呢，這個肯定是錯誤的吧？您幫我看一下吧，謝謝
13.7200 0 24.376247 0
14.8300 0 15.322906 0
16.0000 0 7.284584 0
17.2300 0 0.133717 0
18.5200 0 -6.083034 0
19.8800 0 -11.386774 0
21.3000 0 -15.823056 0
22.7800 0 -19.469415 0
24.3300 0 -22.411644 0
25.9600 0 -24.743574 0
27.6500 0 -26.549735 0
29.4100 0 -27.985790 0
前邊的值是原胞體積，后邊是能量

引用回帖:

Originally posted by zhang_jaj at 2010-12-08 16:49:21:

我也是個菜鳥，對優(yōu)化不清楚，但是計算好像都得優(yōu)化啊。還有那個rm WAVECAR一般是要有的，如果你的計算目錄下沒有WAVECAR時，你也可以不要這個命令行。如果有wavecar時，istart=0將是無效的（好像是，你看看手冊 ...

11樓2010-12-10 09:15:16

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

Originally posted by bluegirlll at 2010-12-10 09:15:16:
您好，真是非常感謝，結(jié)果已經(jīng)出來了�？墒强偰茉趺闯鰜碚摬灰恢钡慕Y(jié)果呢，這個肯定是錯誤的吧？您幫我看一下吧，謝謝
13.7200 0 24.376247 0
14.8300 0 15.322906 0
16.0000 0 7.284584 0
17.2300 0 0.133 ...

你好，我真的很想幫你，但我也不太清楚，我記得手冊還是論壇的帖子里，有提過關(guān)于能量鋸齒狀的原因及解決辦法，這是關(guān)于vasp計算的經(jīng)驗方面的，我接觸才不久，所以幫不了你，希望你早點解決也給我說下，讓我學習下

。

12樓2010-12-10 10:33:09

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
sunyang1988(金幣+1):謝謝交流 2010-12-10 11:06:17

引用回帖:

還有，你的能量值應(yīng)該到后面得都是負的吧，如果這樣的話我覺得后面得應(yīng)該是正確的吧，但我也想不通為什么前面幾個參數(shù)出來的是正值，難道急trail太差了，但是這又為什么能收斂呢？奇怪了。你的結(jié)果應(yīng)該不是鋸齒狀的吧。
期待高手來解決

13樓2010-12-10 10:40:47

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

leo肖

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 702.3
帖子: 97
在線: 28.8小時
蟲號: 898986

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
gzqdyouxia(金幣+1):謝謝交流 2010-12-10 12:39:08

我也用的是 E='tail -1 OSZICAR'

讀出的是最后一個數(shù)據(jù)，用的是tab上面的單引號，你再試試吧，祝你好運

14樓2010-12-10 12:08:09

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

小鉆女

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 114.1
帖子: 8
在線: 4小時
蟲號: 1154849

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
gzqdyouxia(金幣+1):謝謝交流 2010-12-10 16:31:01

你算的是什么物質(zhì)？看你的script好像有兩類原子？是不是你的情況比較復雜所以含有多個參數(shù)？

我也是剛開始學，多個參數(shù)的那個（侯柱峰6.2.2）我也還沒看懂。

15樓2010-12-10 15:37:35

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluegirlll

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 105.9
帖子: 57
在線: 20.6小時
蟲號: 1046082

是那個a的值太小了，就出現(xiàn)正值了，這是為什么呢？我后來往后取值就都是負的能量值了。我現(xiàn)在對k點優(yōu)化，為什么k點要取奇數(shù)呢，取偶數(shù)可以嗎？我取奇數(shù)畫出一個這樣的圖，k點數(shù)和能量值的
數(shù)據(jù)是這樣的，覺得不對
5  10  -31.092686
7  20  -30.859865
9  35  -31.087415
11 56  -31.088405
13 84  -31.087389
15 120 -31.087021
第一列式k點的值，第三列是能量值請高手指點

引用回帖:

Originally posted by zhang_jaj at 2010-12-10 10:40:47:

還有，你的能量值應(yīng)該到后面得都是負的吧，如果這樣的話我覺得后面得應(yīng)該是正確的吧，但我也想不通為什么前面幾個參數(shù)出來的是正值，難道急trail太差了，但是這又為什么能收斂呢？奇怪了。你的結(jié)果應(yīng)該不是鋸齒 ...

16樓2010-12-10 21:26:21

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

Originally posted by bluegirlll at 2010-12-10 21:26:21:
是那個a的值太小了，就出現(xiàn)正值了，這是為什么呢？我后來往后取值就都是負的能量值了。我現(xiàn)在對k點優(yōu)化，為什么k點要取奇數(shù)呢，取偶數(shù)可以嗎？我取奇數(shù)畫出一個這樣的圖，k點數(shù)和能量值的
數(shù)據(jù)是這樣的，覺得不對 ...

確實a太小的話會出現(xiàn)反常的結(jié)果，因為說DFT是考慮的局域作用，太小的話，DFT會做出錯誤的判斷，還有我看你的能量變化是個下拋物線，確實感覺不太對勁，不知道原因，k點去級還是偶看看DFT（density function theory）那本書里有點解釋。搜搜論壇應(yīng)該有，汗啊，幫不了你

17樓2010-12-10 21:45:22

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang_jaj

木蟲 (正式寫手)

應(yīng)助: 14 (小學生)
金幣: 2135.7
帖子: 447
在線: 448小時
蟲號: 1015849

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

引用回帖:

我剛查了下，這種能量的變化趨勢也是有可能得，幫你找到了那本書
http://www.gaoyang168.com/bbs/viewthread.php?tid=1610162
這里的3.1.3里有講。