24小时热门版块排行榜    

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

padian

木虫 (著名写手)

[求助] 请教一个最低能量的问题已有1人参与

请问各位,我用QE计算不同a下Total Energy,由此来根据最低Total Energy来判断Lattice Parameter是多少,不知道这种思路是否正确?

QE计算的Total Energy包括很多部分,是看具体比较不同a下的某一个能量还是系统使用计算给出的Total Energy来进行比较呢?

-------------------------------
我基本上是用example01进行修改的,对Cu在FCC结构下,使用一个Cu原子,变化a从3.5~3.8A(7~7.6),来进行scf计算。现在的计算结果是单调的,并没有存在能量最低点。

请教大家!

谢谢
回复此楼

» 猜你喜欢

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

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

padian

木虫 (著名写手)

我又进一步扩充了a的范围,从3.0~4.0A(间隔为0.5A)
分别scf计算Total Energy
我只使用了1个Cu 原子
计算发现,能量最低点对应的晶格常数在3.2A附近,但是实际上Cu的晶格常数是在3.6A附近。
不知道为什么?我换了好几个赝势,都是差不多在3.2A能量最低。
我的文件如下:

-----------------------------------------------------------------------
    # self-consistent calculation  a = 3.0 A
    cat > Cu300.scf.david.in << EOF
&control
    calculation = 'scf'
    restart_mode='from_scratch',
    prefix='Copper300',
    tstress = .true.
    tprnfor = .true.
    pseudo_dir = '$PSEUDO_DIR/',
    outdir='$TMP_DIR/'
/
&system
    ibrav=  2, celldm(1) =6.00, nat=  1, ntyp= 1,
    ecutwfc =21.5,
    occupations='smearing', smearing='gaussian', degauss=0.02
/
&electrons
    diagonalization='david'
    mixing_mode = 'plain'
    mixing_beta = 0.7
    conv_thr =  1.0d-8
/
ATOMIC_SPECIES
Cu  63.546  Cu.pz-d-rrkjus.UPF
ATOMIC_POSITIONS
Cu 0.00 0.00 0.00
K_POINTS (automatic)
  12 12 12 0 0 0
7楼2014-04-23 22:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖


感谢参与,应助指数 +1
padian(fzx2008代发): 金币+1, 谢谢指导 2014-04-23 21:56:15
你需要找到能量最低点,若是没有,就需要扩展范围。最后得到一条抛物线,用状态方程拟合。
2楼2014-04-23 18:08:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

padian

木虫 (著名写手)

引用回帖:
2楼: Originally posted by mywai520 at 2014-04-23 18:08:22
你需要找到能量最低点,若是没有,就需要扩展范围。最后得到一条抛物线,用状态方程拟合。

所谓能量最低点,是比较每个例子计算出来的Total Energy么?

我是学习例子,在参考书上这个范围内是有最低点的,但是我算出来都是递增的。
3楼2014-04-23 18:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

优化构型应该用'relax'或者'vc-relax',不应该用scf。
4楼2014-04-23 21:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见