Dear Siesta Users and Developers,

I'm trying to calculate band structure of monoclinic (MCL) cell of BiBO3.
Is my choice of points is correct? Also I concern in my choice of
MeshCutoff.

SystemName BiBO3
SystemLabel BiBO3
NumberOfAtoms 20
SpinPolarized .true.
NumberOfSpecies 3
%block ChemicalSpeciesLabel
1  5 B
2  83 Bi
3  8 O
%endblock ChemicalSpeciesLabel

XC.functional GGA
XC.authors PBE

PAO.BasisSize DZP
PAO.EnergyShift 65 meV

LatticeConstant 1.0 Ang
%block LatticeParameters
5.02700 6.58500 8.34900 71.09000 90.000 90.000
%endblock LatticeParameters

AtomicCoordinatesFormat  Fractional
%block AtomicCoordinatesAndAtomicSpecies
0.129000 0.666600 0.029900 3
0.871000 0.333400 0.970100 3
0.629000 0.333400 0.470100 3
0.371000 0.666600 0.529900 3
0.252700 0.792000 0.109200 1
0.747300 0.208000 0.890800 1
0.752700 0.208000 0.390800 1
0.247300 0.792000 0.609200 1
0.110300 0.930700 0.169800 3
0.889700 0.069300 0.830200 3
0.610300 0.069300 0.330200 3
0.389700 0.930700 0.669800 3
0.234090 0.257130 0.122880 2
0.765910 0.742870 0.877120 2
0.734090 0.742870 0.377120 2
0.265910 0.257130 0.622880 2
0.521400 0.756100 0.143000 3
0.478600 0.243900 0.857000 3
0.021400 0.243900 0.357000 3
0.978600 0.756100 0.643000 3
%endblock AtomicCoordinatesAndAtomicSpecies

DM.MixingWeight 0.05
DM.NumberPulay 8
DM.Tolerance 1.d-4
MaxSCFIterations 300

MeshCutoff 300.0 Ry
WriteCoorXmol yes

%block kgrid_Monkhorst_Pack
   6   0   0  0.0
   0   4   0  0.0
   0   0   4  0.0
%endblock Kgrid_Monkhorst_Pack

%block ProjectedDensityOfStates
-30.00 30.00 0.500 500 eV
%endblock ProjectedDensityOfStates

%block BandLines
1 0.000000 0.000000 0.000000 \Gamma
30 0.000000 0.000000 1.000000 Y
10 0.000000 0.831749 1.341764 H
30 0.000000 1.000000 1.000000 C
10 1.000000 1.000000 1.000000 E
30 1.000000 1.168251 0.658236 M1
10 1.000000 1.000000 0.000000 A
30 0.000000 1.000000 0.000000 X
10 0.000000 1.168251 0.658236 H1
0 1.000000 0.831749 1.341764 M
30 1.000000 0.000000 1.000000 D
10 1.000000 0.000000 0.000000 Z
0 0.000000 0.000000 1.000000 Y
30 1.000000 0.000000 1.000000 D
%endblock BandLines

UseSaveData           .true.

SaveRho                 .true.       # Write valence pseudocharge at the
mesh


Best regards,
Maxim Arsentev

Responder a