| 5 | 1/1 | 返回列表 |
| 查看: 2225 | 回復(fù): 4 | ||||
忽悠小皮銅蟲 (小有名氣)
|
[交流]
【求助成功】石墨烯能帶計(jì)算的問題 已有4人參與
|
|
大家好。我是siesta新手,最近學(xué)習(xí)計(jì)算石墨烯的能帶和態(tài)密度。我的計(jì)算結(jié)果很是怪異。請大家?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ā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 298求調(diào)劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +3 | NIFFFfff 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|