【求助】关于QHA模块的使用 - 第一性原理 - QE(Pwscf) - 第2页小木虫论坛-学术科研互动平台

阿什顿联

11楼2010-06-01 11:24:04

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

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

★ ★
identation(金币+14): 2010-06-01 21:01:38
identation(金币+5): 2010-06-02 10:09:04
zzy870720z(金币+2):谢谢交流 2010-06-03 08:21:07

第一个问题我不知道，不过按照你的推测，很可能是带隙，你应该可以查到资料。
第二个问题：确实只需要ldau=.true.和Hubbard_U(i)就足够了，starting_ns_eigenvalue(m,ispin,I)一般不需要特别指定，除非你算出来的这个矩阵值不一样。钙钛矿而言，我估计前面二者足够了，关键是你要选好U值

12楼2010-06-01 20:46:06

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

identation

金虫 (著名写手)

1ST强帖: 16
应助: 46 (小学生)
贵宾: 2.516
金币: 253.4
散金: 2603
红花: 44
帖子: 1283
在线: 500.5小时
虫号: 842155
注册: 2009-09-06
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by goldenfisher at 2010-06-01 20:46:06:
第一个问题我不知道，不过按照你的推测，很可能是带隙，你应该可以查到资料。
第二个问题：确实只需要ldau=.true.和Hubbard_U(i)就足够了，starting_ns_eigenvalue(m,ispin,I)一般不需要特别指定，除非你算出来的 ...

U值的指定是不是在给定的U值计算出的带隙=实验测定的即可？
但是我把U指定成4.6eV，计算得到的带隙并没提高，这是怎么回事呢？
还有，磁性材料在费米能级以下有一段TDOS为零，这是怎么回事呢（下图）？

13楼2010-06-02 10:08:37

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

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

identation(金币+5): 2010-06-07 16:25:04

你这是磁性材料？
你设的U值是对Ti原子？至于为什么没提高，你得自己看了，具体的对象我不了解。
TDOS为零应该就是说这是一个半导体啊。怎么了，不是吗？

14楼2010-06-02 22:00:39

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

identation

金虫 (著名写手)

1ST强帖: 16
应助: 46 (小学生)
贵宾: 2.516
金币: 253.4
散金: 2603
红花: 44
帖子: 1283
在线: 500.5小时
虫号: 842155
注册: 2009-09-06
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by goldenfisher at 2010-06-02 22:00:39:
你这是磁性材料？
你设的U值是对Ti原子？至于为什么没提高，你得自己看了，具体的对象我不了解。
TDOS为零应该就是说这是一个半导体啊。怎么了，不是吗？

对了，还得请教一下，像下图这种，立方和四方相BaTiO3中，Ti原子偏离体心发生位移会使能量下降，这个曲线通过简单地移动Ti原子的位置(例如让Ti的晶体学坐标为0.5，0.5，0.501)比较和Ti在0.5，0.5，0.5位置时的能量就可以实现么？可是我计算的结果是在体心处能量最低呀，这是怎么回事呢？

[ Last edited by identation on 2010-6-7 at 16:30 ]

15楼2010-06-07 16:29:00

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

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

正确的code算出来肯定应该是Ti原子在体心处最低啊，这才是理想的晶格结构，当然应该最低。
我不知道这个图是怎么算出来的，我估计有点问题啊。

16楼2010-06-07 17:12:45

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

hanyanli0475

铜虫 (小有名气)

应助: 11 (小学生)
金币: 173.1
散金: 32
红花: 3
帖子: 296
在线: 130.5小时
虫号: 1160645
注册: 2010-12-02
性别: GG
专业: 凝聚态物性 II ：电子结构

★
identation(金币+1):谢谢参与

引用回帖:

1楼: Originally posted by identation at 2010-05-23 15:49:07:
以AlAs为例，想请教一下具体计算不同温度的热力学性质怎么实现。

以下是espresso-4.2中给出的用QHA计算alas例子中的Edit_Me文件

####################################################################### ...

您好，这几天也在琢磨QHA计算热学性质，我想请教一下，您给出的步骤：
1，自洽
2，声子
3，用QHA计算
我想问的，不用计算相关的声子太密度吗？我怎么查文献时，好像说是得计算声子总的态密度呢？请指教，呵呵！

什么时候能算完！！！！

17楼2011-11-08 19:32:00

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

fzuvivi

金虫 (小有名气)

应助: 4 (幼儿园)
金币: 1168
散金: 93
红花: 1
帖子: 224
在线: 330.1小时
虫号: 881678
注册: 2009-10-23
专业: 凝聚态物性I:结构、力学和

★
identation(金币+1):谢谢参与

引用回帖:

LZ，请问一个问题：用QHA得到的热力学性质，比如振动熵能与实验值相比较吗，也就是能吻合得上吗？？怎么比较？？据我了解，热力学上的熵可是总熵值，也就是它可能包括了振动熵值，电子熵等等；LZ，怎么考虑的？谢谢！！