24小时热门版块排行榜    

查看: 3530  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

» 猜你喜欢

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

已阅   回复此楼   关注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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见