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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助成功】石墨烯能带计算的问题
51/1返回列表
查看: 2152  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

忽悠小皮

铜虫 (小有名气)

[交流] 【求助成功】石墨烯能带计算的问题 已有4人参与

大家好。我是siesta新手，最近学习计算石墨烯的能带和态密度。我的计算结果很是怪异。请大家帮忙看看我的输入文件有什么问题。谢谢！


SystemName       graphene band structure
SystemLabel      g-b
NumberOfAtoms    2
NumberOfSpecies  1

%block ChemicalSpeciesLabel
    1    6  C
%endblock ChemicalSpeciesLabel


MeshCutoff           300. Ry

SolutionMethod        diagon

DM.UseSaveDM         .true.
UseSaveData          .true.

SpinPolarized         false
xc.functional         LDA
xc.authors            CA

SaveDeltaRho                  .true.
SaveRho                       .true.
SaveTotalPotential            .true.
WriteCoorStep                 .true.

%block LocalDensityOfStates
-8.5 1.50 eV
%endblock LocalDensityOfStates

%block ProjectedDensityOfStates
  -8.5  1.5  0.1  800   eV
%endblock ProjectedDensityOfStates

BandLinesScale  ReciprocalLatticeVectors
WriteBands    .true.
WriteKbands  .true.

%block BandLines
1     0.333333   0.333333  0.0000   K
100   0.000000   0.000000  0.0000   \Gamma
100   0.000000   0.500000    0.000    M
100   0.333333   0.333333  0.0000   K
%endblock BandLines

LatticeConstant     2.46    Ang
%block LatticeParameters
    1.0   1.0   10.0    90.  90.  60.
%endblock LatticeParameters

AtomicCoordinatesFormat     Fractional

%block AtomicCoordinatesAndAtomicSpecies
     0.333333     0.333333     0.0000     1
     0.666666     0.666666     0.0000     1
%endblock AtomicCoordinatesAndAtomicSpecies

%block kgrid_Monkhorst_Pack
     30    0    0    0.0
      0   30    0    0.0
      0    0    1    0.0
%endblock kgrid_Monkhorst_Pack

[ Last edited by 忽悠小皮 on 2010-5-20 at 17:29 ]
回复此楼

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

SIESTA

» 猜你喜欢

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

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

valenhou001

至尊木虫 (职业作家)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指导 2010-05-20 08:25:46
可能基矢弄错了。

LatticeConstant     2.46    Ang
%block LatticeParameters
    1.0   1.0   10.0    90.  90.  60.
%endblock LatticeParameters

需改为：


%block LatticeParameters
    2.46   2.46   10.0    90.  90.  120.0
%endblock LatticeParameters

或者：

LatticeConstant     2.46    Ang
%block LatticeVectors
                1.0    0.0    0.0
                -0.5    0.866025   0.0
                0.0    0.0       4.06504
%endblock  LatticeVectors

你也需要检查原子坐标。

%block AtomicCoordinatesAndAtomicSpecies
      0.0000000000000000  0.0000000000000000  0.0000000000000000     1
     0.3333333333333333  0.6666666666666667  0.0000000000000000    1
%endblock AtomicCoordinatesAndAtomicSpecies

[ Last edited by valenhou001 on 2010-5-19 at 21:53 ]
赞一下(2人) 回复此楼
2楼2010-05-19 21:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忽悠小皮

铜虫 (小有名气)
谢谢valenhou001 的提示。
现在搞清楚了，BandLines k点的坐标和我取的基矢不对应。
DOS能量范围也太小。
kgrid也不够精密。
赞一下 回复此楼
3楼2010-05-20 15:03:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jixiao0526

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
初学计算，计算的结果中没有.band文件，哪位高手帮我看看问题出在哪

NumberOfAtoms       9
NumberOfSpecies     2
%block ChemicalSpeciesLabel
1  6   C
2  3   Li
%endblock ChemicalSpeciesLabel
XC.functional  LDA
XC.authors     CA

# =====================================================================================
# Cell & K-MP-grid
LatticeConstant    4.9200 Ang   # Lattice constant alat
%block LatticeVectors               # in units of alat
1.0      0.0000     0.0
-0.5      0.866      0.00
0.0      0.0        3.048780488    # c=15.08000A
%endblock LatticeVectors

#Basis set generation
PAO.BasisSize        DZP
PAO.BasisType        split
PAO.EnergyShift 100 meV

#Molecular Dynamic Calculation
#MD
#MD.TypeOfRun         Broyden
#MD.NumCGsteps        100
#MD.MaxForceTol       0.05 eV/Ang

Functional and Solution Method(Order-N/diagoalization)
SpinPolarized        false
SolutionMethod        diagon
MeshCutoff        150.0 Ry


%block kgrid_Monkhorst_Pack
8 0 0   0.0
0 8 0   0.0
0 0 1   0.5
%endblock kgrid_Monkhorst_Pack
# =====================================================================================
SCF options
MaxSCFIterations      100           # Maximum number of SCF iter
DM.NumberPulay         3            # One Pulay every 5 iter
DM.MixingWeight       0.25          # New DM amount for next SCF cycle
DM.Tolerance          1.d-4         # Tolerance in maximum difference
DM.UseSaveDM          true            # to use continuae..ion files
SolutionMethod        diagon        # OrderN or Diagon

DM.UseSaveDM         true          # Restart
DM.AllowReuse          T
DM.AllowExtrapolation  F

MeshCutoff           150. Ry         # OK for carbon (see Emilio)
OccupationFunction   FD
ElectronicTemperature  250 meV       # Temp. for Fermi smearing


# =====================================================================================
%block BandLine
1   0.0000000000     0.0000000000     0.0000000000     \GAMMA
30   -0.333     0.667     0       F
20   0.0        0.5    0     Q
20   0     0       0.0    Z
%endblock BandLines
# ============================================================================
%block ProjectedDensityOfStates
-10.0  10.0    0.05   100    eV
%endblock ProjectedDensityOfStates
# =====================================================================================
# Output Options
LongOutput   true
BandLinesScale    ReciprocalLatticeVectors

# =====================================================================================
AtomicCoordinatesFormat    Ang  # Format for inates
%block AtomicCoordinatesAndAtomicSpecies < aa.xyz
赞一下(1人) 回复此楼
4楼2011-12-10 13:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
4楼: Originally posted by jixiao0526 at 2011-12-10 13:29:10:
初学计算，计算的结果中没有.band文件，哪位高手帮我看看问题出在哪

NumberOfAtoms       9
NumberOfSpecies     2
%block ChemicalSpeciesLabel
1  6   C
2  3   Li
%endblock ChemicalSpeciesLabel
...

加个WriteBands             T 试试，一般有BandLine 的block 就会有吧
赞一下(1人) 回复此楼
5楼2011-12-16 13:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 忽悠小皮 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭