24小时热门版块排行榜    

查看: 1417  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 321求调剂 +3 大米饭! 2026-03-15 3/150 2026-03-15 17:48 by 哈哈哈哈嘿嘿嘿
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
信息提示
请填处理意见