| 查看: 1190 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
pwscf算費米面为什么计算很慢
|
|||
| pwscf算費米面的原理和算能带不一样么?为什么我算费米面耗时这么长?还有我k点取多的话直接停止中间过程的nscf运算,也没报错.大家算費米面耗时长么?遇到过nscf停止没? |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算 |
» 猜你喜欢
读博
已经有5人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教大师们pwscf和vasp的费米能级分别是如何定义的，一样吗？
已经有6人回复- 护照为什么这么慢？
已经有6人回复
- 关于pwscf的ibrav=0和不为0时的问题
已经有12人回复
- 新手求助pwscf的进行磁性运算和LDA+U运算
已经有35人回复
- 为什么pwscf单节点8核运行比1个核还要慢?
已经有7人回复
- pwscf是如何计算debye温度的
已经有7人回复
- pwscf并行计算时总是隔一段时间重新开始
已经有5人回复
- 计算费米面时遇到的问题
已经有26人回复
- 请问用pwscf可以计算光学性质吗？
已经有8人回复
- pwscf为什么要用ev.x进行物态拟合？
已经有6人回复
- 求PWSCF使用说明书
已经有3人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- 【交流】空穴有费米面吗？
已经有6人回复
- 【原创】PWSCF中的nbnd和degauss
已经有11人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
- 【其他】pwscf中的原胞体积计算！
已经有4人回复
4楼2011-04-23 21:39:36
5楼2011-04-23 21:41:31
★
franch(金币+1): 鼓励交流 2011-04-23 22:51:31
franch(金币+1): 鼓励交流 2011-04-23 22:51:31
|
大概的输入文件: # 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楼2011-04-23 21:58:17