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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » LDA＋U
91/1返回列表
查看: 1121  |  回复: 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 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗？ +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 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
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
关闭
请填处理意见
关闭