HSE bandstructure calculation for AsGa

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

etot.input

   4 1
   JOB = NONSCF
   IN.PSP1 = Ga.SG15.PBE.UPF
   IN.PSP2 = As.SG15.PBE.UPF
   IN.ATOM = atom.config
   XCFUNCTIONAL = HSE
   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

As.SG15.PBE.UPF, Ga.SG15.PBE.UPF

As.SG15.PEB.UPF and Ga.SG15.PBE.UPF are the pseudopotential files.

1. You can submit PWmat tasks in different ways:

   mpirun -np 4 PWmat | tee output

Run the command directly

   #!/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

Submit the task with a pbs script

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