$calculatekp

$calculatekp is the command to solve 6x6 k.p Hamiltonian for nitride or wurtzite based wide bandgap materials.

$calculatekp
all
$\theta$
$\phi$
kx_nmax  ky_nmax dk 
A1 A2 A3 A4 A5 A6 A7  
A1 A2 A3 A4 A5 A6 A7
A1 A2 A3 A4 A5 A6 A7
A1 A2 A3 A4 A5 A6 A7
..... (n layers)
Weighting $\Delta_1$  $\Delta_2$  $\Delta_3$ a c $c_{11}$ $c_{12}$ $c_{13}$ $c_{33}$ $c_{44}$ $c_{66}$
Weighting $\Delta_1$  $\Delta_2$  $\Delta_3$ a c $c_{11}$ $c_{12}$ $c_{13}$ $c_{33}$ $c_{44}$ $c_{66}$
Weighting $\Delta_1$  $\Delta_2$  $\Delta_3$ a c $c_{11}$ $c_{12}$ $c_{13}$ $c_{33}$ $c_{44}$ $c_{66}$
Weighting $\Delta_1$  $\Delta_2$  $\Delta_3$ a c $c_{11}$ $c_{12}$ $c_{13}$ $c_{33}$ $c_{44}$ $c_{66}$
..... (n layers)
$D_1$ $D_2$ $D_3$ $D_4$ $D_5$ $D_6$
$D_1$ $D_2$ $D_3$ $D_4$ $D_5$ $D_6$
$D_1$ $D_2$ $D_3$ $D_4$ $D_5$ $D_6$
$D_1$ $D_2$ $D_3$ $D_4$ $D_5$ $D_6$
.....(n layers)

Weighting is usually defined as 1 if each layer's parameters are calculated by user or GUI program. 

However, if it is not equal to 1, all parameters such as $\Delta_1$ $\Delta_2$ $\Delta_3$ a c $c_{11}$ $c_{12}$ $c_{13}$ $c_{33}$ $c_{44}$ $c_{66}$ $D_1$ $D_2$ $D_3$ $D_4$ $D_5$ $D_6$ will be weighted by

 n_layer_parameters = Weighting(n_layer_parameters) + (1-Weighting)*last_layer_parameters  

For example, if we have 5 layers, and weighting of the second layer is 0.5, then  

 $\Delta_{1} (2)$ =$\Delta_{1} (2) *  Weighting_2  + \Delta_{1} (5) *  (1-Weighting_2)$