24小时热门版块排行榜    

查看: 975  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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部分删除,程序都没有问题
求指点
回复此楼
已阅   回复此楼   关注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

结果太长,就直接把文件贴上来了~
思想重于技巧,内涵重于表象
6楼2011-09-02 10:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

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内存消耗更大。
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贴出来,可能哪里出错了。
4楼2011-08-30 16:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见