| 查看: 1179 | 回复: 5 | ||||
[求助]
pwscf算費米面为什么计算很慢
|
| pwscf算費米面的原理和算能带不一样么?为什么我算费米面耗时这么长?还有我k点取多的话直接停止中间过程的nscf运算,也没报错.大家算費米面耗时长么?遇到过nscf停止没? |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算 |
» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教大师们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人回复
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1551
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
2楼2011-04-23 20:12:24
phys
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 3336.5
- 散金: 209
- 红花: 1
- 帖子: 372
- 在线: 1087.9小时
- 虫号: 1162047
- 注册: 2010-12-03
- 专业: Physics
3楼2011-04-23 21:14:13
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