应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何使用LDA和GGA近似
111/1返回列表
查看: 4704  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP soft study

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-11-15 22:58:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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 ]
赞一下(1人) 回复此楼
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
赞一下(2人) 回复此楼
9楼2010-11-17 18:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zzy870720z

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

另外，在POTCAR文件中有泛函的类型，有关键词的
赞一下(1人) 回复此楼
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 才需要
赞一下(1人) 回复此楼
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的回帖

zzy870720z

荣誉版主 (文坛精英)
★ ★
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等等
泛函指定需要和赝势文件中的一致
赞一下(2人) 回复此楼
7楼2010-11-16 09:31:55
已阅   回复此楼   关注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里面的讲稿。
赞一下(2人) 回复此楼
11楼2013-10-20 10:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭