當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

wilsun7

銀蟲 (初入文壇)

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

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

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

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

第一步：自洽計(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í)會生成文件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)該就可以了。我的問題出現(xiàn)在最后一步，前面3步一切正常。最后一步 bands_FS.x > bands_FS.out 后就沒反應(yīng)了。。。希望大家?guī)臀铱纯磫栴}出在哪兒。謝謝

回復(fù)此樓

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

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

gemucai

木蟲 (正式寫手)

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

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

贊一下

回復(fù)此樓

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

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

查看全部 27 個(gè)回答

gemucai

木蟲 (正式寫手)

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

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

贊一下

回復(fù)此樓

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

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

wilsun7

銀蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線: 45.4小時(shí)
蟲號: 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)在這么多人做鐵硒，都看費(fèi)米面的nesting。但接下來如何做呢？

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

贊一下

回復(fù)此樓

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

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

goldenfisher

金蟲 (著名寫手)

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

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

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo