Bandstructure calculation for Si

atom.config

        2
 LATTICE
      0.00000000     2.71535000     2.71535000
      2.71535000     0.00000000     2.71535000
      2.71535000     2.71535000     0.00000000
 POSITION
  14     0.00000000     0.00000000     0.00000000 0 0 0
  14     0.25000000     0.25000000     0.25000000 0 0 0

   1 4
   JOB = NONSCF
   IN.PSP1 = Si.SG15.PBE.UPF
   IN.ATOM = atom.config
   ECUT = 50
   IN.VR = T
   IN.KPT = T

      BAND                    # COMMENT line
      10                      # number of k-points between X and R
      0.500  0.000  0.000  X  # reciprocal coordinates; label 'X' for X point
      0.500  0.500  0.500  R
      10
      0.500  0.500  0.500  R
      0.500  0.500  0.000  M
      10
      0.500  0.500  0.000  M
      0.000  0.000  0.000  G
      10
      0.000  0.000  0.000  G
      0.500  0.500  0.500  R

      split_kp.x gen.kpt

Si.SG15.PBE.UPF

Si.SG15.PEB.UPF is the pseudopotential file.

1. You can submit PWmat tasks in different ways:

   mpirun -np 4 PWmat | tee output

   #!/bin/bash
   #PBS -N SCF
   #PBS -l nodes=1:ppn=4
   #PBS -q batch
   #PBS -l walltime=100:00:00
 
   ulimit -s unlimited
   cd $PBS_O_WORKDIR
 
   mpirun -np 4 PWmat | tee output

1. After NONSCF calculation, you can run “plot_band_structure.x” to obtain band structure in your current working directory. Then it will generate the following files: bandstructure.eps, bandstructure.png, bandstructure.pdf and bandstructure_1.txt (the data file of band structure), which can be used to plot band by origin or gnuplot. Note, for “spin=2”, another data file bandstructure_2.txt will be generated.

   plot_band_structure.x

Si_Band_Cal