24小时热门版块排行榜    

查看: 1666  |  回复: 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的回帖

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的回帖

mywai520

铁杆木虫 (著名写手)


引用回帖:
3楼: Originally posted by padian at 2014-04-23 18:35:35
所谓能量最低点,是比较每个例子计算出来的Total Energy么?

我是学习例子,在参考书上这个范围内是有最低点的,但是我算出来都是递增的。...

是的
5楼2014-04-23 21:40:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+1, 谢谢指导 2014-04-23 21:56:30
padian: 金币+19, ★★★★★最佳答案, 还有问题请教您 2014-04-23 21:57:02
你现在是确定晶格常数做优化,且只有两个原子,用 scf是没有错的。不需要做原子的弛豫,那样没有意义。
6楼2014-04-23 21:41:38
已阅   回复此楼   关注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的回帖

padian

木虫 (著名写手)

引用回帖:
6楼: Originally posted by mywai520 at 2014-04-23 21:41:38
你现在是确定晶格常数做优化,且只有两个原子,用 scf是没有错的。不需要做原子的弛豫,那样没有意义。

求大侠帮我看一下

模拟和实验差距很大,这是因为QE的Cu赝势,还是我计算输入存在问题?
8楼2014-04-23 22:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


引用回帖:
8楼: Originally posted by padian at 2014-04-23 22:02:13
求大侠帮我看一下

模拟和实验差距很大,这是因为QE的Cu赝势,还是我计算输入存在问题?...

你的0.5值太大了,还要小。0.1试一试
9楼2014-04-23 22:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

padian

木虫 (著名写手)

引用回帖:
9楼: Originally posted by mywai520 at 2014-04-23 22:20:27
你的0.5值太大了,还要小。0.1试一试...

那么计算量就比较多,我先试一试

[ 发自手机版 http://muchong.com/3g ]
10楼2014-04-23 22:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 padian 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见