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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » LDA＋U
91/1返回列表
查看: 1116  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
关闭
请填处理意见
关闭