| 6 | 1/1 | 返回列表 |
| 查看: 1306 | 回復: 5 | |||
[求助]
pwscf算費米面為什么計算很慢
|
| pwscf算費米面的原理和算能帶不一樣么?為什么我算費米面耗時這么長?還有我k點取多的話直接停止中間過程的nscf運算,也沒報錯.大家算費米面耗時長么?遇到過nscf停止沒? |
第一性原理計算 |
金蟲 (著名寫手)
木蟲 (正式寫手)

|
大概的輸入文件: # self-consistent calculation cat > LaFeAsO.scf.in << EOF &control calculation='scf' restart_mode='from_scratch', prefix='LaFeAsO' pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 6, celldm(1) =7.539817722, celldm(3) =2.190380711, nat= 8, ntyp= 4, ecutwfc = 65, occupations='smearing', smearing='gauss', degauss=0.02 / &ELECTRONS conv_thr = 1.D-8, mixing_beta = 0.7 / ATOMIC_SPECIES La 138.9055 La.pbe-nsp-van.UPF Fe 55.845 Fe.pbe-sp-van.UPF As 74.92160 As.pbe-n-van.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS O 0.000000000 0.000000000 0.000000000 O 0.500000000 0.500000000 0.000000000 Fe 0.500000000 0.500000000 1.073662967 Fe 0.000000000 0.000000000 1.073662967 La 0.000000000 0.500000000 0.308554309 La 0.500000000 0.000000000 1.838905585 As 0.000000000 0.500000000 1.379693713 As 0.500000000 0.000000000 0.767184718 K_POINTS automatic 10 10 10 0 0 0 EOF $ECHO " running the scf calculation for Ni...\c" $PW_COMMAND < LaFeAsO.scf.in > LaFeAsO.scf.out check_failure $? $ECHO " done" # # prepare input file ni.fs.in # Sysname='LaFeAsO' nabc=' 16 16 16 ' n_start=33 n_last=36 # E_Fermi=`grep Fermi LaFeAsO.scf.out | cut -c 26-36` a1=`grep 'b(1)' LaFeAsO.scf.out | cut -c 24-54` a2=`grep 'b(2)' LaFeAsO.scf.out | cut -c 24-54` a3=`grep 'b(3)' LaFeAsO.scf.out | cut -c 24-54` cat > kvecs_FS.in < $a2 $a3 $nabc $Sysname EOF $KVECS_COMMAND < kvecs_FS.in > kvecs_FS.out cat > LaFeAsO.fs_NSP.in << EOF &control calculation='nscf' prefix='LaFeAsO', pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 6, celldm(1) =7.539817722, celldm(3) =2.190380711, nat= 8, ntyp= 4, ecutwfc = 65, nbnd= 46 / &electrons diagonalization='cg' conv_thr = 1.0e-10 mixing_beta = 0.7 / ATOMIC_SPECIES La 138.9055 La.pbe-nsp-van.UPF Fe 55.845 Fe.pbe-sp-van.UPF As 74.92160 As.pbe-n-van.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS O 0.000000000 0.000000000 0.000000000 O 0.500000000 0.500000000 0.000000000 Fe 0.500000000 0.500000000 1.073662967 Fe 0.000000000 0.000000000 1.073662967 La 0.000000000 0.500000000 0.308554309 La 0.500000000 0.000000000 1.838905585 As 0.000000000 0.500000000 1.379693713 As 0.500000000 0.000000000 0.767184718 K_POINTS EOF cat kvecs_$Sysname >> LaFeAsO.fs_NSP.in # input file ni.fs.in ready $ECHO " running the Fermi Surface calculation ... \c" $PW_COMMAND < LaFeAsO.fs_NSP.in > LaFeAsO.fs_NSP.out check_failure $? $ECHO " done" # prepare input data (input_FS, Bands.out) for bands_FS mv LaFeAsO.fs_NSP.out Bands_NSP.out cat > input_FS < $E_Fermi $Sysname $nabc $a1 $a2 $a3 EOF $BANDS_COMMAND < Bands_NSP.out > bands_fs.out check_failure $? mv Bands_FS.bxsf LaFeAsO.fs_NSP.bxsf $ECHO " Fermi surface plot: use 'xcrysden --bxsf ni.fs_NSP.bxsf' to plot ...\c" $ECHO " done" $ECHO $ECHO "$EXAMPLE_DIR: done" |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26申博 +4 | 加油沖! 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085602 307分 求調劑 +10 | 不知道叫什么! 2026-03-26 | 10/500 |
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調劑 +10 | 美味可樂雞翅 2026-03-26 | 11/550 |
|
|
[碩博家園] 求調劑 有機化學考研356分 +11 | Nadiums 2026-03-25 | 12/600 |
|
|
[考研] 一志愿華東師范大學有機化學專業(yè),初試351分,復試被刷求調劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
[考研] 346求調劑 一志愿070303有機化學 +11 | 蘿卜燉青菜 2026-03-28 | 12/600 |
|
|
[考研] 086000生物與醫(yī)藥 初試274求調劑 +4 | 小叮當來了 2026-03-30 | 4/200 |
|
|
[考研] 322求調劑 +3 | 熹僖XX 2026-03-31 | 3/150 |
|
|
[考研] 083000環(huán)境科學與工程調劑,總分281 +4 | 橙子(勝意) 2026-03-30 | 4/200 |
|
|
[考研] 一志愿南昌大學324求調劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 324求調劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
[考研] 329求調劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
[考研] 085602 化學工程專碩 340分求調劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 348求調劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 085600,材料與化工321分求調劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 331環(huán)境科學與工程求調劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 07化學303求調劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 網絡空間安全0839招調劑 +4 | w320357296 2026-03-25 | 6/300 |
|