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

wilsun7

應(yīng)助: 2 (幼兒園)
金幣: 233.6
帖子: 45
在線(xiàn): 45.4小時(shí)
蟲(chóng)號(hào): 1863760
注冊(cè): 2012-06-19
專(zhuān)業(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ù)此樓

» 猜你喜歡

一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工調(diào)劑一志愿大連海事085600，349 已經(jīng)有6人回復(fù)
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程調(diào)劑已經(jīng)有11人回復(fù)
311求調(diào)劑一志愿合肥工業(yè)大學(xué) 已經(jīng)有4人回復(fù)
一志愿食品科學(xué)與工程083200求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有10人回復(fù)
一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有15人回復(fù)

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

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

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線(xiàn): 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

錯(cuò)誤來(lái)自于讀入“前一步的pw.x進(jìn)行本征值計(jì)算的out文件 ”�？赡苁恰扒耙徊降膒w.x進(jìn)行本征值計(jì)算的out文件 ”有錯(cuò)誤。

贊一下

回復(fù)此樓

20樓2013-03-13 10:36:24

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

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

gemucai

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

應(yīng)助: 16 (小學(xué)生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線(xiàn): 201.7小時(shí)
蟲(chóng)號(hào): 816478
注冊(cè): 2009-07-27
專(zhuān)業(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
在線(xiàn): 45.4小時(shí)
蟲(chóng)號(hào): 1863760
注冊(cè): 2012-06-19
專(zhuān)業(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
在線(xiàn): 351.5小時(shí)
蟲(chóng)號(hào): 117764
注冊(cè): 2005-11-26
性別: GG
專(zhuān)業(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的回帖

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

[考研] 環(huán)境工程 085701，267求調(diào)劑 +8	minht 2026-03-29	8/400	2026-03-31 03:32 by 蒙奇奇521
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專(zhuān)碩378求調(diào)劑 +3	塔比烏斯 2026-03-30	3/150	2026-03-30 22:55 by 無(wú)際的草原
[考研] 317求調(diào)劑 +8	十閑wx 2026-03-24	8/400	2026-03-30 21:00 by peike
[考研] 化工專(zhuān)碩求調(diào)劑 +5	question挽風(fēng) 2026-03-24	5/250	2026-03-30 18:34 by tianyyysss
[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +11	nnnnnnn5 2026-03-25	11/550	2026-03-29 20:42 by 無(wú)際的草原
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無(wú)際的草原
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專(zhuān)業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 085600，專(zhuān)業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 086502化學(xué)工程342求調(diào)劑 +6	阿姨復(fù)古不過(guò) 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 求調(diào)劑 +4	零八# 2026-03-27	4/200	2026-03-27 18:07 by yu221
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky