24小时热门版块排行榜    

查看: 5268  |  回复: 10

Prettyswan

木虫 (小有名气)


[交流] 【求助】如何使用LDA和GGA近似

http://muchong.com/bbs/viewthread.php?tid=1785190中对于如何使用LDA做计算的问题
引用回帖:
Originally posted by yuxiao_zihan at 2010-01-09 20:09:00:
不设定就可以了 本身GGA-tag只是在做GGA运算的时候采用的

引用回帖:
Originally posted by crystal_bj at 2010-02-17 17:40:07:
可以在INCAR中加入一行 LDA=CA, 或者什么都不写也可以,vasp会去读POTCAR的

1、是不是使用了LDA或者GGA赝势文件之后,除非在INCAR文件中指明要使用哪一种类型的GGA近似,否则系统会默认使用LDA。因为我使用了paw_pbe赝势文件,没有指定GGA-tag,结果屏显文字为
引用回帖:
LDA part: xc-table for Pade appr. of Perdew

2、在侯柱峰的入门指南中说,GGA的选择需要跟赝势类型匹配,采用LDA赝势时不能使用GGA,是不是说赝势是分LDA赝势和GGA赝势的?那么如何区分赝势文件的类型呢?我现在有pot、paw_pbe、paw、paw_gga、pot_GGA五种赝势。

[ Last edited by Prettyswan on 2010-11-16 at 10:52 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )
★ ★
ym23(金币+2):说的太对了 2010-11-16 15:31:36
引用回帖:
Originally posted by Prettyswan at 2010-11-16 09:09:24:

那么请问如何查看泛函类型,关键词是什么?
“最好还是泛函类型和赝势一致”是什么意思,一般来说赝势文件都用同一目录下的就可以了,那么“不同泛函的赝势文件是不一样的,即使能够计算,结果也不是很可靠”是 ...

希望你能看一些基础理论的知识
如果一点都不知道的话对你的计算会有很大的障碍
最多也就是一个只会计算的工具而已

下面是Cu的pw91超软赝势的开始部分
CODE:
  US Cu                                 
11.0000000000000000
parameters from PSCTR are:
   VRHFIN =Cu: d10 p1
   LEXCH  = 91
   EATOM  =  1426.8255 eV,  104.8687 Ry
   TITEL  = US Cu
   LULTRA =        T    use ultrasoft PP ?   
   IUNSCR =        1    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =     .980    partial core radius
   POMASS =   63.546; ZVAL   =   11.000    mass and valenz
   RCORE  =    2.480    outmost cutoff radius
   RWIGS  =    2.480; RWIGS  =    1.312    wigner-seitz radius (au A)
   ENMAX  =  233.729; ENMIN  =  175.297 eV
   EAUG   =  391.451

要学会看这些文件
   LEXCH  = 91泛函类型


楼上已经说明了,GGA是一种近似,在这种近似下面有很多泛函,比如pw91、pbe、wc等等
泛函指定需要和赝势文件中的一致
7楼2010-11-16 09:31:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


引用回帖:
Originally posted by zzy870720z at 2010-11-16 09:31:55:


希望你能看一些基础理论的知识
如果一点都不知道的话对你的计算会有很大的障碍
最多也就是一个只会计算的工具而已

下面是Cu的pw91超软赝势的开始部分
[code]
  US Cu                                 ...

关于“雅克比五个梯子”我倒是曾经看过,也了解到GGA可以分为PW91和PBE,不过印象不深。赝势文件的开头我也有时候看看,不过一些参数在vasp的说明手册也查不到所以暂时没有去细究。谢谢提醒和解答,
另外,LDA part: xc-table for Pade appr. of Perdew是什么意思?在输出的屏显文字中,哪里提示是用什么关联近似泛函呢?
再次感谢。

[ Last edited by Prettyswan on 2010-11-16 at 11:19 ]
8楼2010-11-16 10:51:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


GGA-tag的相关信息


zzy870720z(金币+1):谢谢提供 2010-11-17 21:25:09
今天,我在网上找到了GGA-tag的较为详细的信息,
http://blog.chinaunix.net/u2/65114/showart_575393.html中(省略了在VASP手册上出现的信息)
引用回帖:
…… ……
  如果header部分显示LEXCH = CA,那么就是LDA形式的赝势。
  如果header部分显示LEXCH=91 ,那么就是GGA形式的赝势。
…… ……
另外一个常用的GGA是RPBE,在VASP4.6的版本中,并没有专门针对RPBE的赝势。如果需要用RPBE的交换关联势,只需直接用PBE的POTCAR,然后在INCAR文件里指定GGA=RP
…… ……
在vasp的INCAR文件中,不写GGA关键词时,默认的是CA形式的LDA,在OUTCAR中并会出现警告的提示.CA形式指的是基于 Ceperley and Alder用quantum Monte Carlo计算出来的correlation energy的局域密度近似。参见Ceperley D.M. and Alder B.J. 1980 Phys.Rev.Lett.45 566

9楼2010-11-17 18:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

bingmou(金币+1):感谢交流! 2010-11-15 23:13:21
Prettyswan(金币+7):感谢解答 2010-11-17 17:52:58
pot、paw都是LDA赝势,另外三个是gga赝势
不同泛函的赝势类型不一样
最好还是泛函类型和赝势一致
不同泛函的赝势文件是不一样的
即使能够计算,结果也不是很可靠

另外,在POTCAR文件中有泛函的类型,有关键词的
2楼2010-11-15 23:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skyscorpioa

新虫 (初入文坛)



ym23(金币+1):谢谢提示 2010-11-16 09:08:05
Prettyswan(金币+3): 2010-11-17 17:53:07
PBE 不需要 在 INCAR 中指定 GGA, PW91 才需要
3楼2010-11-16 01:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


引用回帖:
Originally posted by zzy870720z at 2010-11-15 23:07:29:
pot、paw都是LDA赝势,另外三个是gga赝势
不同泛函的赝势类型不一样
最好还是泛函类型和赝势一致
不同泛函的赝势文件是不一样的
即使能够计算,结果也不是很可靠

另外,在POTCAR文件中有泛函的类型,有关键 ...

那么请问如何查看泛函类型,关键词是什么?
“最好还是泛函类型和赝势一致”是什么意思,一般来说赝势文件都用同一目录下的就可以了,那么“不同泛函的赝势文件是不一样的,即使能够计算,结果也不是很可靠”是什么意思呢?
4楼2010-11-16 09:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


引用回帖:
Originally posted by skyscorpioa at 2010-11-16 01:20:51:
PBE 不需要 在 INCAR 中指定 GGA, PW91 才需要

是说PBE不能使用GGA么?
5楼2010-11-16 09:10:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skyscorpioa

新虫 (初入文坛)


引用回帖:
Originally posted by Prettyswan at 2010-11-16 09:10:23:

是说PBE不能使用GGA么?

PBE 是 GGA 的一种, PW91 也是GGA的一种,不过这两种非常接近
6楼2010-11-16 09:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)


引用回帖:
8楼: Originally posted by Prettyswan at 2010-11-16 10:51:24
关于“雅克比五个梯子”我倒是曾经看过,也了解到GGA可以分为PW91和PBE,不过印象不深。赝势文件的开头我也有时候看看,不过一些参数在vasp的说明手册也查不到所以暂时没有去细究。谢谢提醒和解答,
另外,LDA p ...

我用的GGA赝势,而且也设置了GGA=91,但是屏幕上还是有LDA part: xc-table for。。。,这是为什么呢
10楼2013-10-18 16:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

米色空白

金虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
GGA可以认为是LDA的“改良版”,是处理电子有比较明显地局域化的时候所采用的一种近似,或者说是一套方法吧。对于这套方法,不同的人又对其进行了稍有区别的具体描述,这些就是泛函,笼统地讲就是那些POTCAR文件。其中每个POTCAR开头都会有“LEXCH”一行,我才应该就是Exchange Potential的意思吧,其实就是Kohn-Sham方程里面那个交换关联项的具体表述。这一行的内容可能写作“91”,代表的是GGA近似下的PW91泛函,如果写的是“PE”则代表GGA-PBE,当然还可能写作“CA”,那就是LDA的意思了。赝势种类不一而足,具体细节可以参照VASP官方举办的Workshop里面的讲稿。
11楼2013-10-20 10:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 22:18 by xxniao123
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见