24小时热门版块排行榜    

查看: 5289  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见