版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

wilsun7

銀蟲(chóng) (初入文壇)

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

[求助] 計(jì)算費(fèi)米面時(shí)遇到的問(wèn)題

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

以下附上本人的計(jì)算過(guò)程，希望和大家一起研究

第一步：自洽計(jì)算（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)生要計(jì)算費(fèi)米面時(shí)的密集網(wǎng)格k點(diǎn)（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
這時(shí)會(huì)生成文件kvecs_fese 包含了所有k點(diǎn)

第三步：進(jìn)行一次非自洽計(jì)算，計(jì)算這些k點(diǎn)的本征值(cat scf.in kvecs_FeSe > nscf.in)
把 scf 的輸入文件和生成的k點(diǎn)附加到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將計(jì)算的k點(diǎn)以及相應(yīng)的本征值，轉(zhuǎn)換成xcrysden軟件的bxsf格式。
首先把 nscf-FeSe.out 復(fù)制到 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
然后應(yīng)該就可以了。我的問(wèn)題出現(xiàn)在最后一步，前面3步一切正常。最后一步 bands_FS.x > bands_FS.out 后就沒(méi)反應(yīng)了。。。希望大家?guī)臀铱纯磫?wèn)題出在哪兒。謝謝

回復(fù)此樓

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

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

gemucai

木蟲(chóng) (正式寫(xiě)手)

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

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

贊一下

回復(fù)此樓

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

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

wilsun7

銀蟲(chóng) (初入文壇)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

goldenfisher

金蟲(chóng) (著名寫(xiě)手)

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

【答案】應(yīng)助回帖

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

貌似需要運(yùn)行的命令不對(duì)？bands_FS.x > bands_FS.out
不需要輸入文件？

贊一下

回復(fù)此樓

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

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

wilsun7

銀蟲(chóng) (初入文壇)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

gemucai

木蟲(chóng) (正式寫(xiě)手)

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

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(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，但不知道為什么沒(méi)反應(yīng)，輸入命令我是依照前人的，請(qǐng)參考：http://valenhou.blog.edu.cn/2006/133295.html...

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

我用著方法做過(guò)鐵硒費(fèi)米面，沒(méi)問(wèn)題，祝好運(yùn)。

贊一下

回復(fù)此樓

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

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

wilsun7

銀蟲(chóng) (初入文壇)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

gemucai

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時(shí)
蟲(chóng)號(hào): 816478
注冊(cè): 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
謝謝，輸入文件應(yīng)該沒(méi)問(wèn)題，但就是不知道最后一步為什么不行，能否把你做的最后兩步的過(guò)程發(fā)一下，我看看我問(wèn)題出在哪兒，非常感謝！...

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

贊一下

回復(fù)此樓

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

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

wilsun7

銀蟲(chóng) (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時(shí)
蟲(chóng)號(hào): 1863760
注冊(cè): 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ě)漏了。...

非常感謝。但我運(yùn)行后又出現(xiàn)了新問(wè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

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

贊一下

回復(fù)此樓

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

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

goldenfisher

金蟲(chóng) (著名寫(xiě)手)

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

【答案】應(yīng)助回帖

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

應(yīng)該是out這種格式。所以括號(hào)里應(yīng)跟input文件

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

贊一下

回復(fù)此樓

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

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

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

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 281求調(diào)劑（0805） +3	煙汐憶海 2026-03-16	8/400	2026-03-18 00:47 by 煙汐憶海
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 0703化學(xué)336分求調(diào)劑 +4	zbzihdhd 2026-03-15	5/250	2026-03-17 17:33 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 333求調(diào)劑 +3	球球古力 2026-03-11	3/150	2026-03-13 21:27 by JourneyLucky
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu

24小時(shí)熱門(mén)版塊排行榜

[求助] 計(jì)算費(fèi)米面時(shí)遇到的問(wèn)題

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

» 猜你喜歡

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

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖