应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » LDA＋U
91/1返回列表
查看: 968  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ywyou

木虫 (小有名气)

[求助] LDA＋U

最近在学习PWSCF，想用一下里面的LDA＋U，去优化VO2，但是总是出现问题，下面是输出结果：
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 0
[cli_2]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 2
[cli_3]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 3
但当我把原子迟豫改成只进行电子自恰，或者保留原子迟豫而将＋U部分删除，程序都没有问题
求指点
回复此楼

» 猜你喜欢
1楼 2011-08-29 14:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
ywyou(金币+1): 2011-08-30 08:51:22
franch(金币+3): 谢谢回帖交流 2011-08-30 11:05:43
ywyou(金币+4): 2011-08-31 19:25:52
估计是内存或者什么机器的问题。加了U内存消耗更大。
赞一下(1人) 回复此楼
2楼2011-08-29 20:05:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ywyou

木虫 (小有名气)
引用回帖:
2楼: Originally posted by goldenfisher at 2011-08-29 20:05:31:
估计是内存或者什么机器的问题。加了U内存消耗更大。

我算的体系只有12个原子，机器是16G内存，不可能是内存不够的
赞一下 回复此楼
3楼2011-08-30 08:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
zzy870720z(金币+1): 鼓励交流 2011-08-31 09:34:49
那你的信息给出的问题只能是机器问题。要不你把input贴出来，可能哪里出错了。
赞一下(1人) 回复此楼
4楼2011-08-30 16:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ywyou

木虫 (小有名气)
引用回帖:
4楼: Originally posted by goldenfisher at 2011-08-30 16:40:05:
那你的信息给出的问题只能是机器问题。要不你把input贴出来，可能哪里出错了。

谢谢了。这是我的输入文件：
&control
calculation='relax'
restart_mode='from_scratch',
pseudo_dir = './',
outdir='./'
prefix='VO2-M'
tstress = .true.
tprnfor = .true.
/
&system
    ibrav = 12 ,
    celldm(1) = 10.1739 ,
    celldm(2) = 1.069 ,
    celldm(3) = 0.841 ,
    celldm(4) = -0.5394 ,
    nat=12 , ntyp=2 ,
    nbnd=60 ,
    ecutwfc = 60.0, ecutrho = 600.0
    occupations='smearing', smearing='gaussian', degauss=0.02
    nspin=1,
    lda_plus_u=.true., Hubbard_U(1)=3,
/
&electrons
diagonalization='david'
conv_thr = 1.0e-8
mixing_beta = 0.7
/
ATOMIC_SPECIES
V 50.9415   V.pbe-sp-van.UPF
O 15.9994   O.pbe-van_ak.UPF
ATOMIC_POSITIONS {crystal}
V  0.0264599993824959   0.2394700050354004    0.9789400100708007
V  0.4735400006175041   -0.2394700050354004    1.4789400100708008
V  -0.0264599993824959   -0.2394700050354004    -0.9789400100708007
V  0.5264599993824959   0.2394700050354004    -0.4789400100708007
O  0.2085900008678436   0.1061599999666214    0.2118500024080277
O  0.2988399863243103   0.4005100131034851    0.7025799751281738
O  0.2914099991321564   -0.1061599999666214    0.7118500024080278
O  0.2011600136756897   -0.4005100131034851    1.2025799751281741
O  -0.2085900008678436   -0.1061599999666214    -0.2118500024080277
O  -0.2988399863243103   -0.4005100131034851    -0.7025799751281738
O  0.7085900008678436   0.1061599999666214    0.2881499975919724
O  0.7988399863243103   0.4005100131034851    -0.2025799751281738
K_POINTS {automatic}
  4 4 4  0 0 0
赞一下 回复此楼
5楼2011-08-31 14:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
gzqdyouxia(金币+2): 谢谢指教 2011-09-02 15:16:53
ywyou(金币+5): 非常感谢 2011-09-03 08:53:23
尝试了一下，没有问题啊，输入文件：
CODE:
&control
    calculation='relax' ,
    restart_mode='from_scratch' ,
    pseudo_dir = '$PseudoDIR/' ,
    outdir='$TempDIR/' ,
    prefix='VO2-M' ,
    tstress = .true. ,
    tprnfor = .true. ,
    disk_io = 'low' ,
/
&system
    ibrav = 12 ,
    celldm(1) = 10.1739 ,
    celldm(2) = 1.069 ,
    celldm(3) = 0.841 ,
    celldm(4) = -0.5394 ,
    nat = 12 ,
    ntyp=2 ,
    nbnd = 60 ,
    ecutwfc = 60.0 ,
    ecutrho = 600.0 ,
    occupations = 'smearing' ,
    smearing = 'gaussian' ,
    degauss = 0.02 ,
    nspin = 1,
    lda_plus_u = .true.,
    Hubbard_U(1) = 3,
/
&electrons
     diagonalization = 'david' ,
     conv_thr = 1.0e-8 ,
     mixing_beta = 0.7 ,
/
&ions
     ion_dynamics = 'bfgs' ,
/
ATOMIC_SPECIES
  V 50.9415   V.pbe-sp-van.UPF
  O 15.9994   O.pbe-van_ak.UPF
ATOMIC_POSITIONS {crystal}
  V  0.0264599993824959   0.2394700050354004    0.9789400100708007
  V  0.4735400006175041   -0.2394700050354004    1.4789400100708008
  V  -0.0264599993824959   -0.2394700050354004    -0.9789400100708007
  V  0.5264599993824959   0.2394700050354004    -0.4789400100708007
  O  0.2085900008678436   0.1061599999666214    0.2118500024080277
  O  0.2988399863243103   0.4005100131034851    0.7025799751281738
  O  0.2914099991321564   -0.1061599999666214    0.7118500024080278
  O  0.2011600136756897   -0.4005100131034851    1.2025799751281741
  O  -0.2085900008678436   -0.1061599999666214    -0.2118500024080277
  O  -0.2988399863243103   -0.4005100131034851    -0.7025799751281738
  O  0.7085900008678436   0.1061599999666214    0.2881499975919724
  O  0.7988399863243103   0.4005100131034851    -0.2025799751281738
K_POINTS {automatic}
  3 3 3  0 0 0

结果太长，就直接把文件贴上来了～
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
6楼2011-09-02 10:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ywyou

木虫 (小有名气)
引用回帖:
6楼: Originally posted by souledge at 2011-09-02 10:56:40:
尝试了一下，没有问题啊，输入文件：[code]
&control
    calculation='relax' ,
    restart_mode='from_scratch' ,
    pseudo_dir = '$PseudoDIR/' ,
    outdir='$TempDIR/' ,
    prefix='VO2-M'  ...

非常感谢，我刚刚开始学习这个软件，但是总找不到说明书，觉得这个软件过于的繁琐，请问你有没有好一点的说明书？能传一份吗？
赞一下 回复此楼
7楼2011-09-03 08:55:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
zzy870720z(金币+2): 谢谢指点 2011-09-20 12:38:39
我也试了一下，没有问题。我想肯定是你的并行环境有问题。
说明书就在源程序的Doc目录下可以找到，还有example, tutorial可以学习。很多可以慢慢学的资料
赞一下(1人) 回复此楼
8楼2011-09-20 11:07:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
zzy870720z(金币+3): 谢谢侯老师详细指点 2011-09-21 12:52:55
LDA+U基本上计算量的增加是可忽略的。在使用pwscf里面的LDA+U时，先得确定一下源代码里面LDA+U的部分（set_hubbard_l.f90和tabd.f90）是否支持了V这个元素。如果没有的话，得自己增加一下，再编译代码。

另外既然是VO2，应该采用考虑自旋极化的计算，即nspin=2。

另外做'relax'的计算，有关离子的优化方法也得进行设置。

最后，你提交作业的命令可能不大正确，最后也能贴出来进行一下确认。
赞一下(1人) 回复此楼
9楼2011-09-20 21:23:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ywyou 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭