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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】.psf 的格式
31/1返回列表
查看: 1137  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

leila3409

捐助贵宾 (正式写手)

[交流] 【求助】.psf 的格式 已有2人参与

这是H的SIESTA赝势PSF文件
H ca nrl nc  
ATM 3.2.2  9-APR-07 Troullier-Martins
1s 1.00  r= 1.33/2p 0.00  r= 1.33/3d 0.00  r= 0.37/4f 0.00  r= 1.33/
   4  0  864  2.478752176666E-03  1.250000000000E-02    1.000000000000
Radial grid follows
  3.117886413544E-05  6.274991010250E-05  9.471807094134E-05  ...
--------------------------------------------------------------------------------------------------------
请问第一行的第三个字符串nrl是什么意思？是no-relativistic吗？nc呢？是不是no-core correction?谢谢！

另外，用GGA的赝势了，在INPUT里还可以设
xc.functional     LDA
xc.authors     PZ   
吗？这样算起来是没问题，但是合理吗？

[ Last edited by leila3409 on 2011-3-1 at 11:08 ]
回复此楼

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

siesta问题

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-03-01 10:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+2): 你还是那么热心 2011-03-01 14:56:59
引用回帖:
Originally posted by leila3409 at 2011-03-01 03:42:33:
这是H的SIESTA赝势PSF文件
H ca nrl nc  
ATM 3.2.2  9-APR-07 Troullier-Martins
1s 1.00  r= 1.33/2p 0.00  r= 1.33/3d 0.00  r= 0.37/4f 0.00  r= 1.33/
   4  0  864  2.478752176666E-03  1.250000 ...

我觉得你的猜测是对的。另外如果赝势是GGA的，那么计算的时候最好也是GGA的。既然你觉得GGA的好，为什么还要用LDA的泛函呢？如果不统一的话，你的输出文件里面应该有warning的。当然从计算的可行性上来讲，赝势用GGA的价电子的泛函用LDA的，一样算，互不相干。但是从物理上什么意思呢？到底是用的什么呢？
另外在siesta里面也有一种杂化泛函，就是可以将GGA和LDA的混合使用。你的这个情况就属于这种混用的一个极端。
只要你觉得合理，有道理，算的结果好，这个未必不可吧。
赞一下(2人) 回复此楼
弘德明志博学笃行
2楼2011-03-01 14:22:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
Originally posted by zhangguangping at 2011-03-01 14:22:57:
我觉得你的猜测是对的。另外如果赝势是GGA的，那么计算的时候最好也是GGA的。既然你觉得GGA的好，为什么还要用LDA的泛函呢？如果不统一的话，你的输出文件里面应该有warning的。当然从计算的可行性上来讲，赝势 ...

张大牛来了~~~谢谢！
因为我遇到这样一个问题，我在算H2的一些性质以验证我的势和基函数的精确性。我却发现用LDA和GGA的势函数的计算结果几乎一致，但SIESTA的输入文件里面xc.functional 如果选GGA就会出现与实验结果差别很大，如果xc.functional=LDA，计算结果就会与实验结果很稳合。我测试了很多其他的参数，没有找到别的影响因素，只有这个xc.functional 。
可是我将来要算H在fe中的性质，fe的赝势是GGA的，所以我想保持一致，都用GGA的势。
赞一下 回复此楼
3楼2011-03-01 14:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leila3409 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭