 
nextnano^{3}  Tutorial
next generation 3D nano device simulator
1D Tutorial
GaN/AlN wurtzite structure: Strain, piezo and pyroelectric charges in
wurtzite
Author:
Stefan Birner
If you want to obtain the input files that are used within this tutorial, please contact stefan.birner@nextnano.de.
> 1Dwurtzite_nn3_noStrain_noPiezo_noPyro_noPoisson.in
 input file for the nextnano^{3} software
> 1Dwurtzite_nn3.in
/ *_nnp.in  input file for the nextnano^{3} and nextnano++ software
> 1Dwurtzite_nn3_Nface.in
/ *_nnp.in  input file for the nextnano^{3} and nextnano++ software
GaN/AlN structure
 This input file simulates a GaN/AlN/GaN wurtzite structure:
1Dwurtzite.in
The structure is grown pseudomorphically on GaN, i.e. the AlN is
tensilely strained, the GaN is unstrained. The growth direction [0001]
is along z, the interfaces are in the (x,y) plane. As we have [0001] growth
direction, we have Gapolar GaN (i.e. Gaface polarity).
 No strain.
$numericcontrol
simulationdimension
= 1
zeropotential
= yes
varshniparameterson = no !
Band gaps independent of temperature.
Absolute values from database are taken.
latticeconstantstempcoeffon = no !
Lattice constants independent of temperature.
Absolute values from database are taken.
piezoconstantszero
= yes ! Piezoelectric constants are set to
zero.
pyroconstantszero
= yes ! Pyroelectric constants are set to
zero.
$end_numericcontrol
 The pictures shows the conduction and valence band edges of the
heterostructure when no strain is applied (
straincalculation =
nostrain ). AlN is the barrier
for both electrons and holes.
Now we zoom into the valence band structure. In AlN the light hole lies above
the heavy hole, in GaN the situation is the opposite. Note that heavy and
light hole are not degenerate (as would be the case in zinc blende).
 Pseudomorphic strain.
straincalculation =
homogeneousstrain
pseudomorphicon = GaN !
lattice constant a=0.3189 nm (lattice constant c does not matter
for this growth direction)
The GaN is unstrained, the AlN (lattice constant a=0.3112 nm;
elastic constants c_{13}=108.0, c_{33}=373.0) is tensilely
strained as it has the smaller lattice constant.
The biaxial strain is tensile: eps_{ }= (a_{substrate} 
a)
/ a = 0.0247429
The uniaxial strain is compressive: eps_{__}=  2 c_{13}/c_{33}
eps_{ }=  0.0143283
The hydrostatic strain is positive which corresponds to an increase in volume
for AlN: eps_{hy }= Tr(eps_{ij}) = (2 eps_{ }+ eps_{__})
= 0.0351575
The strain leads to a shift in the conduction and valence band edges.
Conduction band at Gamma:
In wurtzite the crystal anisotropy leads to two distinct conduction band
deformation potentials for the Gamma point, one is parallel and the other one
is perpendicular to the c axis.
absolutedeformationpotentialscbs = a_{c,a axis}
a_{c,a axis } a_{c,c axis
} absolutedeformationpotentialscbs =
3.9d0 3.9d0 20.5d0 ! [eV]
Thus the conduction band edge energy including the hydrostatic
energy shift is given by:
E_{c}′ = E_{c} + a_{c,c axis }eps_{__ }
+ 2a_{c,a axis }eps_{} = 4.712 + (20.5
* (0.0143283) ) + 2 (3.9 ) *
0.0247429 = 4.712 + 0.10073553 = 4.81274 eV
So in this particular example, the barrier for electrons is increased.
uniaxcbdeformationpotentialscbs = 0d0 0d0 0d0
Data for uniaxial deformation potentials of other minima than
Gamma are not available yet. At the Gamma point, the uniaxial deformation
potential is zero.
Valence bands:
From a full treatment of the effect of strain on the sixband Hamiltonian
six valence band deformation potentials arise.
uniaxvbdeformationpotentials =
d1 d2
d3 d4 d5 d6 !
[eV]
uniaxvbdeformationpotentials =
17.1 7.9 8.8 3.9
3.4 3.4 ! [eV]
In contrast to zinc blende, an absolute deformation potential for the
valence band is not needed.
absolutedeformationpotentialvb =
0d0 ! (Not needed, should be removed from
wurtzite database in the future.)
The shifts of the valence bands are obtained by diagonalizing the BirPikus
strain Hamiltonian which is given in
Basics 2 (strain
effects). This is a general approach as it gives the correct shifts for
arbitrary orientations (However, it is only for valence bands).
In our example, the tensile strain in AlN shifts all holes upwards,
 the heavy hole by 0.31687 eV,
 the light hole by 0.63565 eV and
 the crystal field splitoff hole by 0.31717 eV,
thus reducing the barrier for the holes.
 Pyroelectric polarization (spontaneous polarization)
$numericcontrol
pyroconstantszero = no
The wurtzite materials GaN, AlN and InN are pyroelectric materials and
thus show a spontaneous polarization.
Thy pyroelectric field P^{py} is along the hexagonal z
direction (along c axis) and it is always negative. P^{py} is
directed from anion to cation. The positive z direction (i.e. [0001]) is from
cation to anion.
GaN: P^{py} =  0.034 C/m²
AlN: P^{py} =  0.090 C/m²
At the interfaces we have a discontinuity of P^{py}(z).
The pyroelectric polarizations at the interfaces are determined as follows:
1^{st} interface at 100 nm (GaN/AlN): P^{py}(left point)
 P^{py}(right point) =  0.034 + 0.090 = 0.056 C/m²
2^{nd} interface at 117 nm (AlN/GaN): P^{py}(left point)  P^{py}(right
point) =  0.090 + 0.034 =  0.056 C/m²
The interface charges can be found in the output file
interface_densitiesD.txt ($outputdensities ).
The interface charge of 0.056 C/m² corresponds to 34.952*10^{12}
electrons /cm².
Once having determined the pyroelectric polarization P^{py} one
is able to compute the pyroelectric charge density: rho_{py}(x) =
 div P^{py}(x).
If the hexagonal axis is oriented along the z axis as in our example, this
equation reduces to: rho_{py}(z) =  d/dz P^{py}(z)
The pyroelectric charge density per cm³ is given in density1Dpyro.dat .
In pyro_polarization.dat the pyroelectric constant as given in
the database or input file is plotted for each grid point.
 Piezoelectric polarization
$numericcontrol
piezoconstantszero = no
If AlN is strained, piezoelectric fields arise.
In GaN the piezoelectric polarization is zero as there is no strain.
AlN: e_{33 }= 1.79, e_{31 }=  0.50 (e_{15} is not
relevant for [0001] growth direction)
In AlN the piezoelectric polarization is:
P^{pz} = e_{33}eps_{__} + e_{31 }(eps_{
}+ eps_{}) = 1.79 * ( 0.0143283)  0.50 * 2 * 0.0247429 =

0.050390 C/m²
P^{pz} is directed in this example parallel to P^{py}
as the piezoelectric polarization is negative for tensilely strained AlN grown
on GaN.
GaN: P^{pz} = 0 C/m² (no strain)
AlN: P^{pz} =  0.050390 C/m²
At the interfaces we have a discontinuity of P^{pz}(z).
The piezoelectric polarizations at the interfaces are determined as follows:
1^{st} interface at 100 nm (GaN/AlN): P^{pz}(left point)
 P^{pz}(right point) = 0 + 0.050390 = 0.050390 C/m²
2^{nd} interface at 117 nm (AlN/GaN): P^{pz}(left point)  P^{pz}(right
point) =  0.050390  0 =  0.050390 C/m²
The interface charges can be found in the output file
interface_densitiesD.txt ($outputdensities ).
The interface charge of  0.050390 C/m² corresponds to 31.451*10^{12}
electrons /cm².
Once having determined the piezoelectric polarization P^{pz}
one is able to compute the piezoelectric charge density: rho_{pz}(x)
=  div P^{pz}(x).
The piezoelectric charge density per cm³ is given in density1Dpiezo.dat .
In our 1D example, this equation reduces to: rho_{pz}(z) =  d/dz P^{pz}(z)
In piezo_polarization.dat the piezoelectric constants (e_{33},
e_{31}, e_{15}) as given in the database or input file are
plotted for each grid point.
 Electrostatic potential including piezo and pyroelectric charges.
The electrostatic potential is the solution of the nonlinear Poisson
equation.
 div (ε(r) grad phi(r)) = rho(r,phi)
The charge density rho contains the (static) piezo and pyroelectric
charge densities as well as the electron and hole charge
densities and ionized donors and acceptors. The latter depend on
the electrostatic potential phi but the piezo and pyroelectric charge
densities don't.
(zeropotential = no )
Here we plot the electrostatic potential in units of eV for three different
cases:
 including pyroelectric charges
(piezoconstantszero = no ,
pyroconstantszero = yes )
 including piezoelectric charges
(piezoconstantszero = yes ,
pyroconstantszero = no )
 including pyro and piezoelectric charges
(piezoconstantszero = no ,
pyroconstantszero = no )
In our example the piezo and pyroelectric contributions are of the same size.
The Fermi level is set to zero.
 The electrostatic potential in units of eV has to be substracted from the
conduction and valence band edges. The resulting band structure is plotted in
this figure (The piezo and pyroelectric charges were considered in this
calculation).
Note that the conduction band is pulled below the Fermi level (around 98 nm)
and the valence band above the Fermi level (around 118 nm).
Nface polarity
This input file simulates a GaN/AlN/GaN wurtzite structure (Nface polarity):
1Dwurtzite_Nface.in
The structure is the same as above. The only difference is the polarity:
Nface instead of Gaface.
The above examples correspond to Gaface polarity. In order to
investigate Nface polarity, we have to switch from [0001] growth
direction to [0001] growth direction (z > z ). In order to keep a
righthanded coordinate axes system we also invert the x direction (x >
x ). The third direction (here: the y direction) is calculated
automatically (internally).
$domaincoordinates
domaintype
= 0 0 1
zcoordinates =
0.0d0 119d0
growthcoordinateaxis = 0 0 1
! along z direction
pseudomorphicon =
GaN
!********************************
! This is along [0001] direction: Gaface polarity
! hkilxdirection = 1 0 1 0
! MillerBravais indices of x coordinate axis [ 1 0 1 0]
!!hkilydirection = 1 2 1 0
! MillerBravais indices of y coordinate axis [1 2 1 0]
! hkilzdirection = 0 0 0 1
! MillerBravais indices of z coordinate axis [ 0 0 0 1]
! Equivalent Bravais indices
! hklxdirectionzb = 1 0 0
! Miller indices of x coordinate axis
[1 0 0]
! hklydirectionzb = 0 1 0
! Miller indices of y coordinate axis [0 1 0]
! hklzdirectionzb = 0 0 1
! Miller indices of z coordinate axis
[0 0 1]
!********************************
!********************************
! This is along [0001] direction: Nface polarity
hkilxdirection = 1 0 1 0
! => x ! MillerBravais indices of x
coordinate axis [1 0
1 0]
! hkilydirection = 1 2 1 0
! MillerBravais indices of y coordinate axis [1 2 1 0]
hkilzdirection = 0 0 0 1
! => z ! MillerBravais indices of z
coordinate axis [ 0 0 0 1]
! Equivalent Bravais indices
! hklxdirectionzb = 1 0 0
! Miller indices of x coordinate axis
[1 0 0]
! hklydirectionzb = 0 1 0
! Miller indices of y coordinate axis [0 1 0]
! hklzdirectionzb = 0 0 1
! Miller indices of z coordinate axis [0 0 1]
!********************************
The figure shows the Fermi level (at 0 eV), the electrostatic potential and the
conduction and valence band edges (heavy and light hole) of the Nface
heterostructure.
(The piezo and pyroelectric charges were considered in this
calculation).
Note the difference in comparison to the figure above for the Gaface case.
The location of the 2DEGs and 2DHGs (twodimension hole gases) is reversed.
An alternative way of switching from Gaface to Nface would have been to
change the sign of the pyroelectric constants.
The boundary conditions for the Poisson equation in this simulation were
Neumann boundary conditions, i.e. the derivative of the electrostatic potential
is zero at the boundaries (flat band boundary conditions).
No doping was present in this simulation.
Note: The figures shown in this tutorial have been generated a long time
ago. The actual results may differ from these figures because the input files
now contain a modified grid space resolution, and modified layer widths, as well
as a different total structure size.
 Please help us to improve our tutorial! Send comments to
support
[at] nextnano.com .
