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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 用siesta算Li在Si中结合能，是否开自旋？
51/1返回列表
查看: 1952  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wtianxing

银虫 (小有名气)

[求助] 用siesta算Li在Si中结合能，是否开自旋？ 已有1人参与

用siesta计算单个Li在64个完整Si晶体中的结合能，公式如下：
E（64Si+Li）-E（64Si）-E（Li）
其中E（64Si+Li）是将Li放入含64个Si的晶体中，用周期边界条件算的，E（64Si）是直接计算含64个Si晶体的能量，E（Li）是将一个锂加三边真空层算的。
问题出在单个锂用真空层算能量上，考虑自旋（SpinPolaried   .true.）后的能量要比不考虑的低0.32eV,而其他两个能量用不用自旋没有影响。
按说考虑自旋后能量更低，应该是对的，但是算出的结合能与其他文献采用vasp（不知道考没考虑自旋）的对不上，就差0.3左右，而不考虑自旋就对上了。请高手帮忙，非常感谢。

算单个Li的输入文件：
Siesta Version: siesta-2.0.1
Architecture  : Ifort(10.1.012)+MKL(10.0.5.025)+OpenMPI(1.2.6)@openSUSE-10.3_X86-64
Compiler flags: mpif90 -O2 -i-static

SystemLabel     work
NumberOfAtoms  1
NumberOfSpecies  1
%block ChemicalSpeciesLabel
    1     3  Li
%endblock ChemicalSpeciesLabel
PAO.BasisType    split
PAO.BasisSize   DZP
PAO.EnergyShift    0.010000 Ry
LatticeConstant 1.0 Ang
%block LatticeParameters
21.6804  23.578  15.30271078  90.000000  90.000000  90.000000
%endblock LatticeParameters
AtomicCoordinatesFormat NotScaledCartesianAng
%block AtomicCoordinatesAndAtomicSpecies
    12.605130 5.880100 8.131973                   1
%endblock AtomicCoordinatesAndAtomicSpecies

#----------------------------------------
# Energy control parameters for SCF cycles
MeshCutoff        200.0 Ry     # Energy cut-off for the planewave
%block kgrid_Monkhorst_Pack    # k mesh points for integration
1   0   0      0.500
0   1   0      0.500
0   0   1      0.500
%endblock kgrid_Monkhorst_Pack
MaxSCFIterations    300        
DM.MixingWeight      0.1     
DM.NumberPulay       10         
DM.Tolerance         1.d-4     
DM.UseSaveDM                     
SolutionMethod      diagon   
ElectronicTemperature   300 K  
XC.functional        GGA
XC.authors           PBE

SpinPolarized    .true.                   开之后能量低0.3左右
#---------------------------------------
# Control the files printed
WriteKpoints        .false.
WriteForces         .true.
WriteCoorXmol       .true.
WriteCoorCerius     .true.
WriteCoorStep       .true.
#---------------------------------------
MD.TypeOfRun         CG   
MD.NumCGsteps        300     
MD.MaxCGDispl         0.15 Ang
MD.MaxForceTol        0.010 eV/Ang
MD.VariableCell       .false.   
Diag.DivideAndConquer  .false.
************************** End of input data file *****************************

[ Last edited by wtianxing on 2013-12-16 at 15:51 ]
回复此楼

» 猜你喜欢

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

sungjen

银虫 (著名写手)
引用回帖:
2楼: Originally posted by wtianxing at 2013-12-17 08:29:53
试着修改K点，真空层尺寸等，两种情况的差值基本没有变化。

很好奇你算的是什么，为什么算这个，但是应该是不加的
赞一下 回复此楼
希望有一天成为大神
3楼2013-12-17 18:19:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

wtianxing

银虫 (小有名气)
试着修改K点，真空层尺寸等，两种情况的差值基本没有变化。
赞一下 回复此楼
2楼2013-12-17 08:29:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wtianxing

银虫 (小有名气)
引用回帖:
3楼: Originally posted by sungjen at 2013-12-17 18:19:00
很好奇你算的是什么，为什么算这个，但是应该是不加的...

算的是单个Li在完整Si晶体中稳定存在时的结合能。背景是锂离子电池，如果采用Si作为负极材料，充电的时候锂就会运动到Si中，所以希望得到Li在硅晶体中的稳定性。
现在结果显示，自旋只对采用真空层算单个Li能量时候有影响，而且比没开自旋的低0.3ev。我咨询过几个人，大家有争论。但是开自旋后能量变低，理论上是不是应该开自旋？能解释下为什么不加吗？
赞一下 回复此楼
4楼2013-12-18 08:36:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
引用回帖:
4楼: Originally posted by wtianxing at 2013-12-18 08:36:07
算的是单个Li在完整Si晶体中稳定存在时的结合能。背景是锂离子电池，如果采用Si作为负极材料，充电的时候锂就会运动到Si中，所以希望得到Li在硅晶体中的稳定性。
现在结果显示，自旋只对采用真空层算单个Li能量时 ...

嗯，理论上是应该加的。为什么不加我也说不清
赞一下 回复此楼
希望有一天成为大神
5楼2013-12-18 10:40:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 332求调剂 +12 蕉蕉123 2026-04-10 12/600 2026-04-10 19:01 by hmn_wj
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 调剂 +4 青灯不负 2026-04-09 4/200 2026-04-10 09:55 by LHGeng
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 308求调剂 +21 倘若起风了呢 2026-04-05 21/1050 2026-04-10 08:13 by Sammy2
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考研] 293分求调剂，外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 285求调剂 +7 恶法大二的气味� 2026-04-05 10/500 2026-04-08 14:34 by zhq0425
[考研] 081700，311，求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 复试调剂 +14 呼呼？~+123456 2026-04-05 14/700 2026-04-06 22:50 by chenzhimin
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭