24小时热门版块排行榜    

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

didi5158

金虫 (小有名气)

[交流] 【求助】结构优化计算了几天,一直不停止,请大家帮忙! 已有7人参与

对晶体进行的结构优化,计算了将近4天了,一直不停止。不知道这是不是正常?如果不正常如何调整参数?有经验的朋友麻烦说道说道,谢谢了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

didi5158

金虫 (小有名气)

引用回帖:
Originally posted by 子虚乌有5388 at 2010-07-29 10:10:55:
把你的input贴出来看看才知道吧

&CONTROL
                 calculation = 'relax' ,
                      outdir = '/root/espresso-4.2/work/Dean/LSAFB/temp/' ,
                      wfcdir = '/root/espresso-4.2/work/Dean/LSAFB/wftemp/' ,
                  pseudo_dir = '/root/espresso-4.2/work/Dean/LSAFB/pp/' ,
                      prefix = 'LSAF' ,
                 lkpoint_dir = .true. ,
               etot_conv_thr = 1.0D-4 ,
               forc_conv_thr = 1.0D-3 ,
                       nstep = 100 ,
                     tstress = .true. ,
                     tprnfor = .true. ,
                          dt = 20 ,
                    dipfield = .true. ,
/
&SYSTEM
                       ibrav = 4,
                   celldm(1) = 5.0716,
                   celldm(3) = 2.0093,
                         nat = 18,
                        ntyp = 4,
                     ecutwfc = 80 ,
                 occupations = 'smearing' ,
                     degauss = 0.02 ,
                    smearing = 'gaussian' ,
/
&ELECTRONS
            electron_maxstep = 100,
                    conv_thr = 1.0D-6 ,
                 mixing_mode = 'plain' ,
                 mixing_beta = 0.7D ,
/
&IONS
                ion_dynamics = 'bfgs' ,
           pot_extrapolation = 'second_order' ,
           wfc_extrapolation = 'second_order' ,
                     upscale = 10D ,
                   bfgs_ndim = 1 ,
            trust_radius_max = 0.8D0 ,
            trust_radius_min = 1.0D-3 ,
            trust_radius_ini = 0.5D ,
/
ATOMIC_SPECIES
   Li    6.94100  Li.pbe-n-van.UPF
   Al   26.98154  Al.pbe-sp-van.UPF
   Sr   87.62000  Sr.pbe-nsp-van.UPF
    F   18.99840  F.pbe-n-van.UPF
ATOMIC_POSITIONS alat
   Sr      0.000000000    0.000000000    0.000000000   
   Sr      0.000000000    0.000000000    1.004630000   
   Li      0.500000000    0.288650000    1.500000000   
   Li      0.000000000    0.577370000    0.502310000   
   Al      0.500500000    0.288650000    0.502320000   
   Al      0.000000000    0.577380000    1.506900000   
    F      0.371500000    0.028400000    0.298380000   
    F      0.789650000    0.307530000    0.298380000   
    F      0.338850000    0.530090000    0.298380000   
    F      0.661150000    0.530090000    0.706260000   
    F      0.628500000    0.028400000    0.706260000   
    F      0.210350000    0.307530000    0.706260000   
    F      0.128500000    0.837620000    1.710890000   
    F     -0.289650000    0.558500000    1.710890000   
    F      0.161150000    0.335930000    1.710890100   
    F     -0.161150000    0.335930000    1.303010000   
    F     -0.128500000    0.837620000    1.303010000   
    F      0.289650000    0.558500000    1.303010000   
K_POINTS automatic
  8 8 8   1 1 1

大家麻烦看看

[ Last edited by didi5158 on 2010-7-29 at 10:51 ]
7楼2010-07-29 10:50:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

didi5158

金虫 (小有名气)

引用回帖:
Originally posted by zzy870720z at 2010-07-29 09:21:08:
可能是你体系太大, 精度太高
自己查一下原因
比如看看能量是不是已经收敛了

能量收敛的挺好的。结构变化也大。
8楼2010-07-29 10:53:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

didi5158

金虫 (小有名气)

引用回帖:
Originally posted by zd777 at 2010-07-29 10:36:30:
pwscf优化就是慢,你看看每一步的优化的结构是什么样的,是不是完全变样了

结构每次与上一次变化不大。
9楼2010-07-29 10:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

didi5158

金虫 (小有名气)

引用回帖:
Originally posted by goldenfisher at 2010-07-29 13:05:02:
bfgs_ndim = 1
            trust_radius_max = 0.8D0 ,
            trust_radius_min = 1.0D-3 ,
            trust_radius_ini = 0.5D ,
这些参数通通可以不要,用默认的即可。这可能是导致你不收敛的原因 ...

非常感谢你的建议,我刚刚接触计算不久,也有人告诉过我具体计算之前将参数测试一下,我也不是太明白如何测试。请教一下具体怎么操作,有无相关的资料提供?先谢谢了!
12楼2010-07-29 19:31:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

didi5158

金虫 (小有名气)

引用回帖:
Originally posted by 子虚乌有5388 at 2010-07-29 11:57:38:
ecutwfc = 80

这个条件太苛刻了吧,相当于1088eV,你确信用得着???

非常感谢你的建议
13楼2010-07-29 19:32:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 didi5158 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见