版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

wilsun7

銀蟲 (初入文壇)

應助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時
蟲號: 1863760
注冊: 2012-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 計算費米面時遇到的問題

大家好，本人新手，在計算費米面時遇到一個問題：在運行 bands_FS.x > bands_FS.out 時沒有任何反應，查看系統(tǒng)監(jiān)視器，bands_FS.x 根本沒有運行。希望大神幫我看看問題處在那。計算是參考的這個網(wǎng)站： http://valenhou.blog.edu.cn/2006/133295.html

以下附上本人的計算過程，希望和大家一起研究

第一步：自洽計算（pw.x < scf-FeSe.in > scf-FeSe.out）
輸入文件：
&CONTROL
               calculation = 'scf' ,
            restart_mode = 'from_scratch' ,
               wf_collect = .true. ,
                  outdir = './tmp' ,
               pseudo_dir = '/home/wilson/Downloads/espresso-4.3.1/Psudoptential_files',
                  prefix = 'FeSe' ,
               verbosity = 'high'
/
&SYSTEM
                     ibrav = 0,
               celldm(1) = 7.1148,
                     nat = 4,
                     ntyp = 2,
                     nspin = 1,
                  ecutwfc = 30.0 ,
                  ecutrho = 300.0 ,
               occupations = 'smearing' ,
                  degauss = 0.01 ,
                  smearing = 'methfessel-paxton' ,
/
&ELECTRONS
electron_maxstep = 500,
                  conv_thr = 1.0e-10 ,
               mixing_beta = 0.1 ,
         diagonalization = 'david' ,
/
CELL_PARAMETERS
   1.000000000 0.000000000 0.000000000
   0.000000000 1.000000000 0.000000000
   0.000000000 0.000000000 1.465600000
ATOMIC_SPECIES
Fe  55.85 Fe.pbe-sp-van.UPF
Se  78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS crystal
Fe 0.250000 0.75000 0.00000
Fe 0.750000 0.25000 0.00000
Se 0.250000 0.25000 0.26000
Se 0.75000 0.75000 0.74000
K_POINTS automatic
8  8  6 1 1 1

第二步：用kvecs_FS.x產(chǎn)生要計算費米面時的密集網(wǎng)格k點（kvecs_FS.x < kvecs_FS.in > kvecs_FS.out）
kvecs_FS.in 文件如下：
kvecs_FS.in
1.000000 0.000000 0.000000
0.000000 1.000000 0.000000
0.000000 0.000000 0.682314
16 16 16
FeSe

然后 kvecs_FS.x < kvecs_FS.in > kvecs_FS.out
這時會生成文件kvecs_fese 包含了所有k點

第三步：進行一次非自洽計算，計算這些k點的本征值(cat scf.in kvecs_FeSe > nscf.in)
把 scf 的輸入文件和生成的k點附加到nscf中。

輸入文件：
&CONTROL
               calculation = 'nscf' ,
            restart_mode = 'from_scratch' ,
               wf_collect = .true. ,
                  outdir = './tmp' ,
               pseudo_dir = '/home/wilson/Downloads/espresso-4.3.1/Psudoptential_files',
                  prefix = 'FeSe' ,
               verbosity = 'high'
/
&SYSTEM
                     ibrav = 0,
               celldm(1) = 7.1148,
                     nat = 4,
                     ntyp = 2,
                     nspin = 1,
                  ecutwfc = 30.0 ,
                  ecutrho = 300.0 ,
               occupations = 'smearing' ,
                  degauss = 0.01 ,
                  smearing = 'methfessel-paxton' ,
/
&ELECTRONS
electron_maxstep = 500,
                  conv_thr = 1.0e-10 ,
               mixing_beta = 0.1 ,
         diagonalization = 'david' ,
/
CELL_PARAMETERS
   1.000000000 0.000000000 0.000000000
   0.000000000 1.000000000 0.000000000
   0.000000000 0.000000000 1.465600000
ATOMIC_SPECIES
Fe  55.85 Fe.pbe-sp-van.UPF
Se  78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS crystal
Fe 0.250000 0.75000 0.00000
Fe 0.750000 0.25000 0.00000
Se 0.250000 0.25000 0.26000
Se 0.75000 0.75000 0.74000
K_POINTS
  4913
0.000000 0.000000 0.000000  1.00
0.000000 0.000000 0.042645  1.00
0.000000 0.000000 0.085289  1.00
   .....             .....          ......          ...
再run pw.x < nscf-FeSe.in > nscf-FeSe.out

第四步：用bands_FS.x將計算的k點以及相應的本征值，轉(zhuǎn)換成xcrysden軟件的bxsf格式。
首先把 nscf-FeSe.out 復制到 Bands.out
然后新建文件input_FS ：
20 24
8.8449
FeSe
16 16 16
1.000000 0.000000 0.000000
0.000000 1.000000 0.000000
0.000000 0.000000 0.682314

再run bands_FS.x > bands_FS.out
然后應該就可以了。我的問題出現(xiàn)在最后一步，前面3步一切正常。最后一步 bands_FS.x > bands_FS.out 后就沒反應了。。。希望大家?guī)臀铱纯磫栴}出在哪兒。謝謝

回復此樓

1樓 2012-06-27 06:44:38

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

還是等牛人吧，我很驚訝現(xiàn)在這么多人做鐵硒，都看費米面的nesting。但接下來如何做呢？

贊一下

回復此樓

2樓2012-06-27 08:46:07

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

wilsun7

銀蟲 (初入文壇)

應助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時
蟲號: 1863760
注冊: 2012-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by gemucai at 2012-06-27 08:46:07
還是等牛人吧，我很驚訝現(xiàn)在這么多人做鐵硒，都看費米面的nesting。但接下來如何做呢？

我也只是看paper時看到有人做，然后練習一下而已。。

贊一下

回復此樓

3樓2012-06-27 22:13:52

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★ ★
感謝參與，應助指數(shù) +1
wilsun7: 金幣+1 2012-06-27 23:24:31
xueht987: 金幣+1, 謝謝回帖，鼓勵交流 2012-06-29 22:04:01

貌似需要運行的命令不對？bands_FS.x > bands_FS.out
不需要輸入文件？

贊一下

回復此樓

4樓2012-06-27 23:15:32

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

wilsun7

銀蟲 (初入文壇)

應助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時
蟲號: 1863760
注冊: 2012-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by goldenfisher at 2012-06-27 23:15:32
貌似需要運行的命令不對？bands_FS.x > bands_FS.out
不需要輸入文件？

輸入文件是input_FS，但不知道為什么沒反應，輸入命令我是依照前人的，請參考：http://valenhou.blog.edu.cn/2006/133295.html

贊一下

回復此樓

5樓2012-06-27 23:24:08

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應助回帖

★ ★ ★
感謝參與，應助指數(shù) +1
WDD880227: 金幣+1, 多謝提示 2012-06-28 09:36:32
wilsun7: 金幣+2 2012-06-28 10:43:15

引用回帖:

5樓: Originally posted by wilsun7 at 2012-06-27 23:24:08
輸入文件是input_FS，但不知道為什么沒反應，輸入命令我是依照前人的，請參考：http://valenhou.blog.edu.cn/2006/133295.html...

對了，不要完全按照他的來做，費米面都是我自己弄的。
他的shell腳本可要可不要，你自己復制粘貼是一樣的效果。其中的a1,a2,a3其實就是你*.scf.out里的倒空間基矢b1,b2,b3，復制粘貼就行了。nabc之類的他都給你了，費米能級不用我說應該就知道上哪兒找吧？

我用著方法做過鐵硒費米面，沒問題，祝好運。

贊一下

回復此樓

6樓2012-06-28 08:26:37

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

wilsun7

銀蟲 (初入文壇)

應助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時
蟲號: 1863760
注冊: 2012-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

6樓: Originally posted by gemucai at 2012-06-28 08:26:37
對了，不要完全按照他的來做，費米面都是我自己弄的。
他的shell腳本可要可不要，你自己復制粘貼是一樣的效果。其中的a1,a2,a3其實就是你*.scf.out里的倒空間基矢b1,b2,b3，復制粘貼就行了。nabc之類的他都給你了 ...

謝謝，輸入文件應該沒問題，但就是不知道最后一步為什么不行，能否把你做的最后兩步的過程發(fā)一下，我看看我問題出在哪兒，非常感謝！

贊一下

回復此樓

7樓2012-06-28 10:43:03

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

gemucai

木蟲 (正式寫手)

應助: 16 (小學生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
fzx2008: 金幣+1, 謝謝回帖！ 2012-06-28 20:50:09

引用回帖:

7樓: Originally posted by wilsun7 at 2012-06-28 10:43:03
謝謝，輸入文件應該沒問題，但就是不知道最后一步為什么不行，能否把你做的最后兩步的過程發(fā)一下，我看看我問題出在哪兒，非常感謝！...

最后一步是
bands_FS.x bands_FS.out
侯老師寫漏了。

贊一下

回復此樓

8樓2012-06-28 11:03:07

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

wilsun7

銀蟲 (初入文壇)

應助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時
蟲號: 1863760
注冊: 2012-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

8樓: Originally posted by gemucai at 2012-06-28 11:03:07
最后一步是
bands_FS.x <bands.out> bands_FS.out
侯老師寫漏了。...

非常感謝。但我運行后又出現(xiàn)了新問題。。
輸出結(jié)果是：
forrtl: severe (59): list-directed I/O syntax error, unit 5, file stdin
Image             PC             Routine          Line       Source
bands_FS.x       000000000048005A  Unknown             Unknown  Unknown
bands_FS.x       000000000047EB56  Unknown             Unknown  Unknown
bands_FS.x       0000000000441070  Unknown             Unknown  Unknown
bands_FS.x       000000000040C35E  Unknown             Unknown  Unknown
bands_FS.x       000000000040B89F  Unknown             Unknown  Unknown
bands_FS.x       00000000004289B3  Unknown             Unknown  Unknown
bands_FS.x       00000000004271B7  Unknown             Unknown  Unknown
bands_FS.x       000000000040566C  MAIN__                   353  bands_FS.f90
bands_FS.x       0000000000402BCC  Unknown             Unknown  Unknown
libc.so.6       00007F3D8EED1C4D  Unknown             Unknown  Unknown
bands_FS.x       0000000000402AC9  Unknown             Unknown  Unknown

感覺像是系統(tǒng)編譯的問題，但不知道具體怎么回事。。。

贊一下

回復此樓

9樓2012-06-28 22:49:58

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★ ★
wilsun7: 金幣+1, ★有幫助 2012-06-29 11:24:42
xueht987: 金幣+1, 謝謝回帖！ 2012-06-29 22:04:28

應該是out這種格式。所以括號里應跟input文件

[ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

贊一下

回復此樓

10樓2012-06-28 23:17:17

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wilsun7 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 280求調(diào)劑 +4	咕嚕曉曉 2026-03-18	4/200	2026-03-18 10:27 by zeng1010
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	10/500	2026-03-18 09:21 by 咪嗚喵嗚
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[考研] 268求調(diào)劑 +7	好運連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0703 物理化學調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 297求調(diào)劑 +4	學海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 學碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 0703化學求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu

24小時熱門版塊排行榜

[求助] 計算費米面時遇到的問題

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

» 猜你喜歡

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

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖