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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » LDOS结果和PDOS不一致
111/2返回列表12下一页
查看: 2021  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lemonare

新虫 (小有名气)

[求助] LDOS结果和PDOS不一致

我接触Siesta时间不长，利用Siesta计算纳米材料的电子结构，因为不喜欢把想要算的东西全部放入fdf文件中一次全部算完，所以通常是算完自洽弛豫后，再利用弛豫后的结果算想要的性质。
最近遇到了一个问题，我在计算自洽弛豫时通常会同时计算PDOS，然后根据PDOS，我再算LDOS，可是最近出现了PDOS显示有电子密度的能量范围，LDOS显示电子密度为0，两者的结果不一致，即使重新计算也没有变化，我不明白是那里出了错误...
不知道有没有人遇到过这种错误，大家讨论一下。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
每一天都是一个成长的过程！
1楼2012-02-20 10:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lemonare

新虫 (小有名气)
SystemName          CoreShell
SystemLabel         CoreShell

NumberOfAtoms       300
NumberOfSpecies     3

%block ChemicalSpeciesLabel
1  48  Cd
2  34  Se
3  52  Te
%endblock ChemicalSpeciesLabel

PAO.BasisSize       DZP
PAO.EnergyShift     50 meV

LatticeConstant     1.0 Ang
%block LatticeVectors
      55.0000        0.00000       0.00000
      0.00000        55.0000       0.00000
      0.00000        0.00000       7.61600
%endblock LatticeVectors

%block kgrid_Monkhorst_Pack
   1   0    0    0.0
   0   1    0    0.0
   0   0    4    0.0
%endblock kgrid_Monkhorst_Pack

#%block ProjectedDensityOfStates
#   -7.00  -3.00  0.050  500  eV
#%endblock ProjectedDensityOfStates

%block LocalDensityOfStates
     -4.90  -4.85  eV
%endblock LocalDensityOfStates

XC.functional           GGA
XC.authors              PBE
SpinPolarized           .false.
MeshCutoff              300. Ry
MaxSCFIterations        200
DM.MixingWeight         0.15
DM.Tolerance            1.d-5

DM.NumberPulay          5

SolutionMethod         diagon
ElectronicTemperature  0 meV

MD.TypeOfRun          cg
MD.NumCGsteps         200
MD.MaxCGDispl         0.10 Ang
MD.MaxForceTol        0.04 eV/Ang

WriteEigenvalues         .true.
UseSaveData              .true.

AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
-2.551311         0.000000         0.000000         1
1.275655         2.209500         0.000000         1
1.275655         -2.209500         0.000000         1
-1.275655         2.209500         3.606000         1
-1.275655         -2.209500         3.606000         1
2.551311         0.000000         3.606000         1
......

PDOS是-7~-3eV范围内，显示-4.90~-4.85eV内有电子密度，可是我用LDOS算时，却没有电子密度...感觉上是LDOS有错误，不过不知道是什么地方错了
赞一下(1人) 回复此楼
每一天都是一个成长的过程！
3楼2012-02-22 08:53:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

vv_china

铜虫 (小有名气)
可以给一个实例吗？把PDOS和LDOS的参数给出来看看。
赞一下 回复此楼
2楼2012-02-22 01:06:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
fzx2008(金币+1): 谢谢交流！ 2012-02-22 20:35:01
lemonare(金币+20): ★★★很有帮助 5 2012-02-23 08:55:20
建议先只算一下gamma点试试，看能否吻合。
赞一下 回复此楼
弘德明志博学笃行
4楼2012-02-22 19:58:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
解决问题了？
回复此楼
弘德明志博学笃行
5楼2012-02-23 12:02:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemonare

新虫 (小有名气)
没有，计算还没有得到结果
赞一下 回复此楼
每一天都是一个成长的过程！
6楼2012-02-24 12:27:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemonare

新虫 (小有名气)
引用回帖:
楼: Originally posted by zhangguangping at 2012-02-23 12:02:44:
解决问题了？

尝试了，依旧没能解决，还是搞不清楚问题所在...
赞一下 回复此楼
每一天都是一个成长的过程！
7楼2012-02-25 10:26:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)
奇怪哦。我想问问用siesta如何算体系中某几个原子的LDOS？
赞一下 回复此楼
8楼2012-02-27 00:10:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemonare

新虫 (小有名气)
%block LocalDensityOfStates
     -4.90  -4.85  eV
%endblock LocalDensityOfStates

LDOS是能量E和位置r的函数，所以对能量积分的话，就能得到电子密度，Siesta就是对一个能量范围积分，得到的结果就是电子密度，不过计算的是全部原子，而不是几个原子，几个原子的DOS，可以用PDOS功能
赞一下 回复此楼
每一天都是一个成长的过程！
9楼2012-02-27 09:16:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)
fzx2008(金币+1): 谢谢回帖交流！ 2012-02-28 12:29:02
但是我看到你的PDOS参数也只是给出了能量范围。
赞一下 回复此楼
10楼2012-02-27 22:49:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lemonare 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭