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)

回复此楼