24小时热门版块排行榜

>论坛更新日志 (4611)
>休闲灌水 (479)
>虫友互识 (388)
>导师招生 (343)
>文献求助 (331)
>论文投稿 (188)
>博后之家 (174)
>招聘信息布告栏 (134)
>硕博家园 (67)
>学术会议 (65)
>考博 (59)
>找工作 (45)
>基金申请 (40)
>教师之家 (34)
>育儿交流 (33)
>有机交流 (30)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

pariah

铜虫 (小有名气)

应助: 6 (幼儿园)
金币: 41.4
散金: 5
红花: 4
帖子: 135
在线: 171.3小时
虫号: 1320067
注册: 2011-06-10
专业: 光学

[求助] relax计算停住了。大家看看是怎么回事

输入文件如下，如果把mp网格取在3 3 3 1 1 1 那么程序能给出计算结果，如果是6 6 6 1 1 1 就会一直停在一个地方。

CODE:

 &control

    calculation='relax'

    restart_mode='from_scratch',

    pseudo_dir = '/home/Xiaole/Desktop/Pwscf/espresso-4.3.2/pseudo/',

    outdir='/home/Xiaole/Desktop/Pwscf/espresso-4.3.2/work/work_bamgf4/tmp/',

    prefix='bamgf4'

    tprnfor = .true.

    tstress = .true.

 /

 &system

    ibrav=  0, celldm(1) =1.889725989, nat=  12, ntyp= 3, ecutwfc =40.0,

    occupations='smearing',smearing='gaussian',degauss=0.001

 /

 &electrons

    conv_thr    = 1.D-6

    mixing_beta = 0.5D0

 /

 &IONS

    ion_dynamics='bfgs'

 /

ATOMIC_SPECIES

 Ba  137.33 Ba.pbe-nsp-van.UPF

 Mg   24.32 Mg.pbe-nsp-bpaw.UPF

  F   19.00  F.pbe-n-van.UPF

CELL_PARAMETERS {cubic}

    2.0594950000    7.2313000000    0.0000000000

   -2.0594950000    7.2313000000    0.0000000000

    0.0000000000    0.0000000000    5.8115100000

ATOMIC_POSITIONS {angstrom}

 Ba      0.0000000000    -2.1481299780    -2.6947971870

 Mg     -2.0594950000    -1.2510149000     0.0000000000

  F     -2.0594950000    -2.7999593600    -1.1448674700

  F     -2.0594950000    -2.3950065600     1.5144795060

  F     -2.0594950000    -0.3948289800    -1.7911073820

  F      0.0000000000    -1.1295290600    -0.0586962510

 Ba      0.0000000000     2.1481299780     0.2109578130

 Mg      0.0000000000    -5.9802851000    -2.9057550000

  F      0.0000000000    -4.4313406400     1.7608875300

  F      0.0000000000    -4.8362934400    -1.3912754940

  F      0.0000000000    -6.8364710200     1.1146476180

  F      0.0000000000     1.1295290600     2.8470587490

K_POINTS {automatic}

   6 6 6  0 0 0

输出的最后几行如下
      k( 61) = (  -0.1618520  -0.0230480  -0.0286787), wk = 0.0370370
      k( 62) = (  -0.1618520  -0.0230480  -0.0573574), wk = 0.0370370
      k( 63) = (  -0.1213890  -0.0115240  -0.0286787), wk = 0.0370370
      k( 64) = (  -0.1213890  -0.0115240  -0.0573574), wk = 0.0370370
      k( 65) = ( 0.0000000 0.0460959  -0.0286787), wk = 0.0185185
      k( 66) = ( 0.0000000 0.0460959  -0.0573574), wk = 0.0185185
      k( 67) = (  -0.2023150  -0.0115240  -0.0286787), wk = 0.0370370
      k( 68) = (  -0.2023150  -0.0115240  -0.0573574), wk = 0.0370370

   Dense  grid: 39989 G-vectors    FFT dimensions: (  60,  60,  45)

回复此楼

» 猜你喜欢

1楼 2012-04-06 09:40:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★ ★
fzx2008: 金币+2, 专家考核, 谢谢回帖交流！ 2012-04-07 10:12:38

引用回帖:

12楼: Originally posted by souledge at 2012-04-07 01:20:36:
不过一般计算资源不足不是容易卡在
“Starting wfc are xxxx randomized atomic wfcs”
这里么~
LZ的体系并不大，卡在Dense grid xxxx G-vectors这里，怎么看也不像是计算资源不足的样子，毕竟都使用了6个核心 ...

卡在“Dense grid xxxx G-vectors”正式由于计算资源不足的表现，那是因为使用的k点太多所致。你单机算一个例子，如果你使用很密的k点，也很容易导致这个问题。何况楼主的体系，对单机而言，体系不算小，电子数目也不少。
还有一个问题，lz使用的是USPP的赝势，却没有定义ecutrho. ecutrho默认是ecutwfc的四倍。而且体系中有Ba这个元素，ecutwfc=40貌似有些小。最好做一个测试。