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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】siesta输入文件中赝势文件的作用？
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 3736  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

watermall

金虫 (小有名气)

[交流] 【求助】siesta输入文件中赝势文件的作用？

最近要用siesta，但发现siesta输入文件关于描述原子的赝势有点不明白
例如在siesta home/Tests/graphite_c6_full 这个例子中，
输入文件为：
===============================
SystemLabel graphite_c6_full
NumberOfAtoms        4

NumberOfSpecies      1
%block ChemicalSpeciesLabel
    1   6 C.pbe
%endblock ChemicalSpeciesLabel

%block PS.lmax
   C.pbe 2
%endblock PS.lmax

LatticeConstant 1.0 Ang
%block LatticeParameters
    2.4610    2.4610    6.7080   90.0000   90.0000  120.0000
%endblock LatticeParameters

AtomicCoordinatesFormat Fractional
%block AtomicCoordinatesAndAtomicSpecies
    0.000000000    0.000000000    0.000000000     1
    0.000000000    0.000000000    0.500000000     1
    0.333333333    0.666666667    0.000000000     1
    0.666666667    0.333333333    0.500000000     1
%endblock AtomicCoordinatesAndAtomicSpecies

MM.Cutoff 20.0 Ang
%block MM.Potentials
1 1 C6 15.19 1.56
%endblock MM.Potentials

XC.functional GGA
XC.authors PBE

PAO.EnergyShift  0.0025 Ry
PAO.SplitNorm    0.30
PAO.SplitNormH   0.50

PAO.SoftDefault true
PAO.SoftPotential 50.0 Ry
PAO.SoftInnerRadius 0.80

kgrid_cutoff 16.0 Ang

PAO.BasisSize TZ2P

ElectronicTemperature 0 K

%block GridCellSampling
0.5    0.5    0.5
%endblock GridCellSampling

DM.MixingWeight      0.30
DM.NumberPulay       5

MD.TypeofRun CG
MD.VariableCell true
MD.NumCGsteps 100
MD.MaxForceTol 0.01 eV/Ang
MD.MaxStressTol 200 bar

LongOutput true
UseSaveData true
=================================
既然已经在block ChemicalSpeciesLabel 中定义了 C的赝势，为啥后面还要给出
XC.functional GGA
XC.authors PBE
PAO.EnergyShift  0.0025 Ry
PAO.SplitNorm    0.30
PAO.SplitNormH   0.50
PAO.SoftDefault true
PAO.SoftPotential 50.0 Ry
PAO.SoftInnerRadius 0.80
PAO.BasisSize TZ2P
这么多基组的信息？

我以前用的是pwscf，用的是基于平面波的赝势，计算的时候，选择好要用的赝势，再选择合适的截断能量，就可以计算了。有点搞不明白siesta用赝势和基组的原理，不知道我的问题有没有问清楚，但还是希望有人能够帮忙解决下。
回复此楼

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

第一性原理相关文档 SIESTA

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-16 21:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linso

银虫 (初入文坛)
★ ★ ★
watermall(金币+1):谢谢参与
zzy870720z(金币+2):谢谢讨论 2010-09-17 08:22:22
watermall(金币+10): 2010-09-17 13:31:12
以下纯粹是我的个人理解，如果有误还请不吝指正:-)
所有对问题的回答均以  ## --> 开头,  ##在fdf文件中是注释。


最近要用siesta，但发现siesta输入文件关于描述原子的赝势有点不明白
例如在siesta home/Tests/graphite_c6_full 这个例子中，
输入文件为：
===============================
SystemLabel graphite_c6_full
NumberOfAtoms        4

NumberOfSpecies      1
%block ChemicalSpeciesLabel
    1   6 C.pbe                              ##--> 程序将读取C.pbe.psf文件
%endblock ChemicalSpeciesLabel

%block PS.lmax
   C.pbe 2                                    ##--> for C.pbe, max KB projector = 2
%endblock PS.lmax

LatticeConstant 1.0 Ang
%block LatticeParameters
    2.4610    2.4610    6.7080   90.0000   90.0000  120.0000
%endblock LatticeParameters

AtomicCoordinatesFormat Fractional
%block AtomicCoordinatesAndAtomicSpecies
    0.000000000    0.000000000    0.000000000     1
    0.000000000    0.000000000    0.500000000     1
    0.333333333    0.666666667    0.000000000     1
    0.666666667    0.333333333    0.500000000     1
%endblock AtomicCoordinatesAndAtomicSpecies

MM.Cutoff 20.0 Ang
%block MM.Potentials
1 1 C6 15.19 1.56                           ##--> Grimme type Van der Vaals经验修正的参数  !! attention,  Grimme type Van der Vaals经验修正 现在到了Version 3(DFTD3), siesta这里是version 2  
%endblock MM.Potentials

XC.functional GGA      
XC.authors PBE

PAO.EnergyShift  0.0025 Ry
PAO.SplitNorm    0.30
PAO.SplitNormH   0.50

PAO.SoftDefault true
PAO.SoftPotential 50.0 Ry
PAO.SoftInnerRadius 0.80

kgrid_cutoff 16.0 Ang

PAO.BasisSize TZ2P

ElectronicTemperature 0 K

%block GridCellSampling
0.5    0.5    0.5
%endblock GridCellSampling

DM.MixingWeight      0.30
DM.NumberPulay       5

MD.TypeofRun CG
MD.VariableCell true
MD.NumCGsteps 100
MD.MaxForceTol 0.01 eV/Ang
MD.MaxStressTol 200 bar

LongOutput true
UseSaveData true
=================================
既然已经在block ChemicalSpeciesLabel 中定义了 C的赝势，为啥后面还要给出
XC.functional GGA               ##--> 即使赝势是LDA的, 这里也可以指定为GGA,程序此时程序使用GGA进行计算，但使用LDA的赝势， 不对应对不对， 呵呵， siesta  从赝势文件中读取生成赝势的functional 类型后不一定以后用那种类型，缺省是LDA, 我的意思是， 即使是GGA赝势文件，如果这里不设的话，程序会用LDA function 计算，但使用的是GGA赝势(嗨，说得很拗口，不知你能理解这些话不，如果你看原代码得话，倒是很快能理解)
XC.authors PBE                 ##--> 原因同上， 缺省LDA的XC.authors 是PZ
                                         ##--> 如果赝势和程序中使用的functional 不对应， 输出中应该有类似 "WARNING: Pseudopotential generated with" 的
PAO.EnergyShift  0.0025 Ry     ##--> 这些参数都是用来
PAO.SplitNorm    0.30              ##--> 调整basis set的，可以控制basis的形状和 用于spilt 位置
PAO.SplitNormH   0.50             ##--> basis set 优化，
PAO.SoftDefault true               ##--> 详见 siesta 07年的workshop
PAO.SoftPotential 50.0 Ry        ##--> siesta\workshop\18-22 June 2007\19 June\Bases.pdf
PAO.SoftInnerRadius 0.80         ##--> 和
PAO.BasisSize TZ2P                  ##--> siesta\workshop\18-22 June 2007\19 June\BasisOptimization.pdf
这么多基组的信息？

我以前用的是pwscf，用的是基于平面波的赝势，计算的时候，选择好要用的赝势，再选择合适的截断能量，就可以计算了。
##--> 在siesta 中不一定， basis set和赝势都要optimise才会有很好的结果（一般定性结果的话，不优化是没问题的）
其中赝势影响小，basis set影响大。截断能量是kgrid_cutoff 或MeshCutoff，用来控制grid的密度， 这个cutoff 与pwscf 中ecutrho, ecutrho 或 ecutwfc相比，算是有类似关系， pwscf 的ecut{rho|wfc} 是控制kinetic energy 截尾的，也就是波函数\psi = \sum C_{k+G} | k+ G >,  (K+G)^2 <=  ecutwfc，siesta 的kgrid_cutoff 或MeshCutoff,用来控制R空间的格点密度，\delta X = 2 \pi / MeshCutoff

有点搞不明白siesta用赝势和基组的原理，
##--> siesta 的赝势是NC(normal conservation， 没有pwscf的uspp和gipaw), 一般用TM型(actually, Improved Troullier-Martins， siesta没有pwscf的RRKJ型)
##--> 基组basis set是NAO(numerical atomic orbital), 简单的说就是多个原子波叠加， 每一个L（角动量）可以有多个orbital, TZ2P 是每个占据L是3个orbital，外加Lmax+1 的2个orbital
以Ar[1s2 2S2 2P6 3S2 3P6 ] 原子为例， 价电子层是 3S2 3P6， 那么
TZDP 是 S(l=0, m=1), 3*1个； P(l=1, Lmax, m=-1,0,+1), 3*3个 ， D(l=2, Lmax+1, Polarise   orbital, m=-2,-1,0,1,2 ) 2*5个， 共 3+9+10 =22
对于你的C.pbe ， 以及 NumberOfAtoms        4
则是( 3*1+3*3+2*5 ) *4 共 88 orbital

不知道我的问题有没有问清楚，但还是希望有人能够帮忙解决下。
##--> 不知道我说清楚了没有，如果有错误还请不吝指正。

[ Last edited by linso on 2010-9-17 at 07:21 ]
赞一下(2人) 回复此楼
2楼2010-09-17 00:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
★
watermall(金币+1):谢谢参与
LS 高人，学习了
回复此楼
3楼2010-09-17 08:59:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)
★
watermall(金币+1):谢谢参与
引用回帖:
Originally posted by linso at 2010-09-17 00:43:04:
以下纯粹是我的个人理解，如果有误还请不吝指正:-)
所有对问题的回答均以  ## --> 开头,  ##在fdf文件中是注释。


最近要用siesta，但发现siesta输入文件关于描述原子的赝势有点不明白
例如在siesta hom ...

楼上的说的很细
回复此楼
4楼2010-09-17 21:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★
watermall(金币+1):谢谢参与
2楼的回答很详细，学习了。
赞一下(1人) 回复此楼
弘德明志博学笃行
5楼2010-09-18 09:25:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
★
watermall(金币+1): 谢谢参与
引用回帖:
337218楼: Originally posted by linso at 2010-09-17 00:43:04:
以下纯粹是我的个人理解，如果有误还请不吝指正:-)
所有对问题的回答均以  ## --> 开头,  ##在fdf文件中是注释。


最近要用siesta，但发现siesta输入文件关于描述原子的赝势有点不明白
例如在siesta hom ...

LDA计算使用GGA的赝势，这个该怎么理解？

难道是核电子用GGA，价电子用LDA？
赞一下 回复此楼
6楼2012-05-03 12:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljzhou86

金虫 (正式写手)
★
watermall(金币+1): 谢谢参与
studying ing
回复此楼
7楼2012-07-04 17:38:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
★
watermall(金币+1): 谢谢参与
引用回帖:
2396254楼: Originally posted by byin at 2012-05-03 12:01:39
LDA计算使用GGA的赝势，这个该怎么理解？

难道是核电子用GGA，价电子用LDA？...

个人理解，这两者并不矛盾啊，不管是LDA还是GGA都是描述原子核的势而已，两者只是在数值上有差异，并不影响价电子使用哪种近似方法。
赞一下 回复此楼
坠入无边的网~~
8楼2012-07-15 11:14:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

s09221109

金虫 (初入文坛)
★
watermall(金币+1): 谢谢参与
引用回帖:
2楼: Originally posted by linso at 2010-09-17 00:43:04
以下纯粹是我的个人理解，如果有误还请不吝指正:-)
所有对问题的回答均以  ## --> 开头,  ##在fdf文件中是注释。


最近要用siesta，但发现siesta输入文件关于描述原子的赝势有点不明白
例如在siesta home/Tes ...

TZDP 是 S(l=0, m=1)中m是不是应该=0呀？
上面不是说外加Lmax+1嘛，为什么你的计算中没有。我是初学者，什么都不懂，只是有些疑问。
赞一下 回复此楼
9楼2014-08-29 15:55:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 watermall 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工调剂303分，过四级 +14 栖梧待风 2026-04-02 14/700 2026-04-02 21:58 by 科研小专家
[考研] 366求调剂一志愿东北大学 +8 运气来得若有似� 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[考研] 0703化学 +5 goldtt 2026-04-02 7/350 2026-04-02 21:32 by dongzh2009
[考研] 一志愿武汉理工0856，初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 265求调剂 +12 林深温澜 2026-04-01 15/750 2026-04-02 19:14 by 一只好果子?
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 建环，能源，土木老师路过看一看！！！ +4 嘿嘿uu 2026-04-01 4/200 2026-04-01 20:42 by 无懈可击111
[考研] 290求调剂 +5 dfffsar 2026-03-29 5/250 2026-04-01 19:45 by 6781022
[考研] 0710生物学考研调剂 +3 李多米lee. 2026-03-27 4/200 2026-04-01 16:21 by zzchen2000
[考研] 309求调剂 +19 谁不是少年 2026-03-29 19/950 2026-04-01 15:47 by jp9609
[考研] 283求调剂 +9 A child 2026-03-28 9/450 2026-04-01 14:20 by Jaylen.
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! +9 真名有冰 2026-03-29 10/500 2026-03-31 18:01 by xhai2011
[考研] 070300化学354求调剂 +15 101次希望 2026-03-28 15/750 2026-03-31 17:58 by jp9609
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 343求调剂 +8 爱羁绊 2026-03-28 8/400 2026-03-31 16:12 by 不吃魚的貓
[考研] 0703化学 +20 妮妮ninicgb 2026-03-27 20/1000 2026-03-31 13:33 by 无际的草原
[考研] 调剂求院校招收 +7 鹤鲸鸽 2026-03-28 7/350 2026-03-31 11:21 by oooqiao
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
信息提示
关闭
请填处理意见
关闭