| 5 | 1/1 | 返回列表 |
| 查看: 2222 | 回復(fù): 4 | ||||
忽悠小皮銅蟲 (小有名氣)
|
[交流]
【求助成功】石墨烯能帶計(jì)算的問題 已有4人參與
|
|
大家好。我是siesta新手,最近學(xué)習(xí)計(jì)算石墨烯的能帶和態(tài)密度。我的計(jì)算結(jié)果很是怪異。請(qǐng)大家?guī)兔纯次业妮斎胛募惺裁磫栴}。謝謝! 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 |
至尊木蟲 (職業(yè)作家)
|
可能基矢弄錯(cuò)了。 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 你也需要檢查原子坐標(biāo)。 %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 ] |
銅蟲 (小有名氣)
金蟲 (小有名氣)
|
初學(xué)計(jì)算,計(jì)算的結(jié)果中沒有.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 |
鐵桿木蟲 (正式寫手)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +11 | chaojifeixia 2026-03-19 | 12/600 |
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|