24小时热门版块排行榜

返回列表

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

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)

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

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

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

★ ★
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-04-07 10:13:40

引用回帖:

13楼: Originally posted by goldenfisher at 2012-04-07 02:19:32:
卡在“Dense grid xxxx G-vectors”正式由于计算资源不足的表现，那是因为使用的k点太多所致。你单机算一个例子，如果你使用很密的k点，也很容易导致这个问题。何况楼主的体系，对单机而言，体系不算小，电子数目 ...

用台机（一直都用台机，老板到现在不批钱至少买个高端台机）算过60来个原子体系的，用7x7x7的K点，USPP赝势，45Ry+450Ry，0.01Ry的MP Smearing，耗费内存10GB左右，使用六个核心，卡在开始生成波函的地方十几分钟以上才开始进入SCF步~不过限制G向量网格时真没碰到问题~
原子中有Ti和W，都是考虑了半芯态生成的~整个电子数目还是非常多的~
这个过程有点类似VASP中，如果是在FFT planing那里挂的，可能就是编译问题多点；如果是reading WAVCAR挂的，可能是编译问题，也可能是资源问题了，要看具体情况了~
个人理解~