24小时热门版块排行榜    

查看: 1774  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见