Finite Element Method Linear Rectangular Element for Solving Finite Nanowire Superlattice Quantum Dot Structures

The International Journal of Engineering and Science (IJES)
|| Volume || 6 || Issue || 4 || Pages || PP 01-07 || 2017 ||
ISSN (e): 2319 – 1813 ISSN (p): 2319 – 1805
DOI: 10.9790/1813-0604010107 www.theijes.com Page 1
Finite Element Method Linear Rectangular Element for Solving
Finite Nanowire Superlattice Quantum Dot Structures
1
Eman Ali Hussain, 2
Jamil A. Al-Hawasy and 3
Lamyaa H. Ali
1
Asst.Prof., Dr., Department of Mathematics, college of Science, University of AL-Mustansiriyah, Iraq
2
Asst.Prof., Dr., Department of Mathematics, college of Science, University of AL-Mustansiriyah, Iraq
3
Lec., Department of Mathematics, college of Science, University of AL-Mustansiriyah, Iraq
--------------------------------------------------------ABSTRACT-----------------------------------------------------------
This paper concerned with the solution of finite nanowire superlattice quantum dot structures
with a cylindrical cross-section determine by electronic states in various type of layers in terms of wave
functions between structures containing the same number of barriers and wells (asymmetrical) or containing a
different number (symmetrical). The solution is considered with the Finite element method with different base
linear rectangular element to solve the one electron Ben Daniel-Duke equation. The results of numerical
examples are compared for accuracy and efficiency with the finite difference method of this method and finite
element method of linear triangular element. This comparison shows that good results of numerical examples.
---------------------------------------------------------------------------------------------------------------------------
Date of Submission: 17 February 2017 Date of Accepted: 07 April 2017
---------------------------------------------------------------------------------------------------------------------------
I. INTRODUCTION
The modeling of the electron states in semiconductor nanostructures remains a difficult computational task.
In 1966, Ben Daniel and Duke, [1] are used space-charge effects on electron tunneling of the one electron model
is formulated to describe Ben Daniel-Duke equation. There are various methods to solve this type of problems
such as in 1988 Ghatak et. al., [2] used transfer matrix method. In 1990, Paasch et. al., [3] used envelope
equation and wave function matching for narrow-gap semiconductors. In 1995, Mathine et. al., [4] applied
computational Fourier series solution of the BenDaniel-Duke Hamiltonian for arbitrary shaped quantum wells.
In 2003, Melnik et al, [5] applied finite element analysis of Nanowire superlattice structures. Whereas 2004,
Melnik et al, [6] using finite difference method in cylindrical polar coordinates Nanowire superlattice. In 2016,
Deyasi et. al., [7] applied numerically computed in presence of electric field using propagation matrix method.
In this paper we propose a Finite element method and using different base linear rectangular element for solving
Ben Daniel-Duke equation. This solution is considered with the solution of finite nanowire superlattice quantum
dot structures of a cylindrical cross-section and determine by electronic state in the two of
layers (the well) and (the barrier), three of layers with well layers at both ends and six of layers
with well layers at one end and barrier layers at the other end, in terms of wave functions, between structures
containing the same number of barriers and wells (asymmetrical) and structures where the number of barrier
layers is one above the number of well layers (symmetrical).
II. DESCRIPTION OF THE BEN DANIEL-DUKE EQUATION[5],[6]
The Nanowire superlattice (NWSL) is modeled as an ideal cylinder with sharp modulations in the longitudinal
(or -direction). The theory of the electron state used here is the one-band envelope function theory. In the
application, the effective mass of the electron is dependent on the position, but does not depend on energy. The
electronic structure can be found in general based on a multiband model for the envelope function:
where is the Hamiltonian, is called eigenfunction, and is called energy eigenvalue.
Quantum mechanical problems in semiconductor heterostructures are often treated using a single band effective
mass equivalent Hamiltonian. This means in a more general that only envelope functions as solutions of the
equivalent Ben Daniel-Duke equation:

Finite Element Method Linear Rectangular Element for Solving Finite Nanowire Superlattice ..
where is the effective mass in each superlattice layer, is the wave function, is the potential
experienced by electron it can taken to be infinite outside the wire (vacuum) and piecewise inside, and is the
energy.
The NWSL is modeled as a cylinder of infinite length, with alternating layers of (the well) and
(the barrier) see Fig.(1). The NWSL's currently being grown are free standing. Eq. (2) is supplemented by
Dirichlet boundary conditions (on the surface of the nanowire structures) which correspond to the case where
the nanostructure environment is vacuum.
(a) (b)
Fig. (1) (a)Nanowire superlattice model for a three layers, with well layers at both ends and (b) The domain for
a unit cell of nanowire superlattice.
In this work we study a finite sequence of wells and barriers, with infinite potential barriers at the ends, three
patterns are possible:
1. The ends are both wells,
2. At one end there is a well and at the other there is a barrier,
3. The ends are both barriers.
The structure in the first and third case is symmetric with respect to a center plane, and the wavefunctions are
symmetric or antisymmetric. While the structure is asymmetric in the in the second case with respect to a center
plane.
III. FINITE ELEMENT SOLUTION
To solve the one-electron Ben Daniel–Duke equation we write the gradient and the divergence operators in
cylindrical coordinates, i.e.
Since the cylindrical is symmetry, and where is an integer. The Ben
Daniel–Duke equation, Eq. (3), becomes
with the boundary condition
where is a bounded open domain in the plane with boundary
The Ben Daniel–Duke equation governing physical process in a cylindrical Nanowire superlattice model are
described in terms of cylindrical coordinates. When the geometry, loading, and boundary conditions are
independent of the circumferential direction ( -coordinate in Fig. (1)), the problem is said to be axisymmetric
and the Ben Daniel–Duke equation becomes in two-dimensional representing by and
3.1 Weak Form
Assume that is a typical element, whether rectangular or quadrilateral, of the finite element mesh, and we
develop the finite element model of (4) over . Various two-dimensional elements will be discussed in the
sequel.
There are three steps in the development of the weak form of (4) over the typical element .

The first step is formulated (4) by the following weak form
Where
and
In the second step we distribute the differentiation among and equally. To achieve this, we integrate the first
terms in (4a) by parts. We note the identities
Next, we use the component form of the gradient (or divergence) theorem on the boundary and using (6),
and (7) in (5), we obtain
where is an outer unite vector normal on and is the length of an infinitesimal line element
along the boundary.
From an inspection of the boundary integral in (8), we note that the specification of constitutes the essential
boundary condition, and hence is the primary variable. The specification of the coefficient of the weight
function in the boundary expression, i.e. we can let
The third and last step of the formulation is to substitute the definition (9) in (8) and write the weak form of (4)
as
3.2 Finite Element Model of Ben Daniel-Duke equation
Suppose that is approximated over a typical finite element by the expression
where is the value of at the node of the element and are the Lagrange interpolation
functions.
Substituting the finite element approximation (11) for into the weak form (10), we get
This equation must hold for every admissible choice of weight function Since we need independent
algebraic equations to solve for the unknowns, we choose linearly independent functions for
For each choice of we obtain an algebraic relation among . We label the
algebraic equation resulting from substitution of into (12):

Problem (13) is supplemented by Neumann's boundary conditions if . In the following subsections, we
will solve the above meant form of the Ben Daniel-Duke equation by FEM using different base linear
rectangular element and also in the next subsection we will solve it using FEM with linear rectangular element.
3.3 Linear Rectangular Elements
A solution of the weak form of Ben Daniel-Duke equation (10) and finite element matrices in (12) shows that
should be at least four linearly independent term and is linear in and , with a bilinear term in and in
is of the form
This polynomial requires an element with four nodes. The linear rectangular element is a compatible element
because on any side varies only linearly and there are two nodes to uniquely define it.
Here consider an approximation of the form (14) and use a rectangular element with sides and . For the sake
of convenience, choose a local cylindrical polar coordinate system to derive the interpolation functions.
Assume that
And require
where
or, in concise form,
where are the cylindrical polar coordinate of node The interpolation function have the
properties
Now, to solve Ben Daniel–Duke equation, using linear rectangular element, the problem will reduce to the
following generalized eigenvalue problem.
Where matrices and are matrices, is the vector of unknowns of dimensionality , and is the
number of nodes.
where
and

To determine the element coefficient matrices and in (19) and (20), consider the rectangular element show
the Fig (2).
Fig.( 2) The rectangular element.
As a first choice the uniform mesh are used by four linear rectangular element, as shown in Fig.( 3) to represent
the domain (mesh Rect1), to determine the element coefficient matrices and of Eq. (18).
The elements 1, 2, 3, and 4 as shown in Fig.( 3) are identical in orientation as well as geometry. Then all four
elements have the same element matrices, and it is necessary to compute them only for element 1.
To evaluation of the integral in Eq. (19) and Eq. (20), it is possible to obtain the closed form for the and
matrices by carrying out the integrals in Eq. (19) and Eq. (20) exactly respectively. We consider element 1 as
the typical element. Hence, the element coefficient and matrices are
Where and
The element matrices in Eq. (21) and Eq. (22) are valid for the Ben Daniel–Duke equation (4) on any right-angle
triangle with sides and .
Fig.( 3) The uniform mesh Rect1 of four linear rectangular element.
Then the mesh shown in Fig. (3), will be as
where

The above assembled coefficient matrix for finite element mesh is of order . In this case the sub matrices
(the integrals) of compute by an analytic method. The MATLAB program is used to solve the sub matrices
(the integrals) of .
It is important to mention that the elements of the above sub matrices (the integrals) of are calculated by an
analytic method
Notes (1):
In the following examples, the representative set of parameters were used for the calculations were as:
Radius These
parameters correspond to a structure.
Example (1):
Consider the Nanowire superlattice model, with the two layers, at one end there is a well and at the other there is
a barrier (asymmetrical).
In this example the linear rectangular element was used to solve Ben Daniel-Duke equation. Equation (18) is
solved using MATLAB program at . The results are shown in Fig (4).
Fig. (4) Ground state for the NWSL calculated with the FEM of linear rectangular element, FEM of linear
triangular element, and the FDM, discretization of the two layers, with at one end there is a well and at the other
there is a barrier, asymmetrical.
Example (2):
Consider the Nanowire superlattice model for a three layers, with well layers at both ends (symmetrical).
solved using MATLAB program at . The results are shown in Fig (5).

triangular element, and the FDM, discretization of the three layers, with well layers at both ends, symmetrical.
Example (3): Consider the Nanowire superlattice model for a six layers, at one end there is a well and at the
other there is a barrier (asymmetrical).
solved using MATLAB program at . The results are shown in Fig. (6).
triangular element, and the FDM, discretization of the six layers asymmetrical.
IV. CONCLUSION
The finite element method with using different base linear rectangular element are used for solving the one
electron Ben Daniel-Duke equation with the solution of finite nanowire superlattice quantum dot structures
with a cylindrical cross-section and determine by electronic states in various type of layers in
terms of wave functions between structures containing the same number of barriers and wells (asymmetrical) or
containing a different number (symmetrical). The results of numerical examples are gives a good accuracy and
efficiency of this method compare the results between finite difference methods and FEM of linear triangular
element.
REFERENCES

[1]. D. J. BenDaniel and C. B. Duke ‘‘space-charge effects on electron tunneling’’J. of Physical Review, Vol. 152, No. 2 (1966), pp.
683-692.
[2]. A. K. Ghatak, K. Thyagarajan and M. R. Shenoy ‘‘A novel numerical technique for solving the one dimensional schroedinger
equation using matrix approach ‑ application to quantum well structures’’IEEE Journal of Quantum Electronics, Vol. 24, No.
8 (1988), pp. 1524‑1531.
[3]. G. Paasch, P. H. Nguyen and G. Gobsch ‘‘Envelope Equation and Wave Function Matching for Narrow-Gap
Semiconductors’’Physica Status Solidi (b), Vol. 162, No. 1 (1990), pp. 155–163.
[4]. D. L. Mathine, S. K. Myjak and G. N. Maracas ‘‘A Computational Fourier Series Solutin of the BenDaniel-Duke Hamiltonian
for Arbitrary Shaped Quantum Wells’’IEEE Journal of Quantum Electronics, Vol. 31, (1995), pp. 1216‑1222.
[5]. M. Willatzence, R. Melnik, C. Galeriu and L. C. Lew Yan Voon ‘‘Finite Element Analysis of Nanowire Superlattice
Structures’’Springer, Vol. 31, (2003), pp. 755‑763.
[6]. C. Galeriu, L. C. Lew Yan Voon, R. Melnik and M. Willatzence ‘‘Modeling a Nanowire Superlattice Using the Finite
Difference Method in Cylindrical Polar Coordinates’’Computer Physics Communications, Vol. 157, (2004), pp. 147‑159.
[7]. D. Sarkara, A. Deyasib ‘‘Field Induced Tuning of DOS and Eigenstates in Double Quantum Well Structure having Gaussian
Geometry’’American Scientific Publishers , Vol. 5, No. 1, (2016), pp. 138‑143.

Finite Element Method Linear Rectangular Element for Solving Finite Nanowire Superlattice Quantum Dot Structures

More Related Content

What's hot (20)

Similar to Finite Element Method Linear Rectangular Element for Solving Finite Nanowire Superlattice Quantum Dot Structures (20)

Recently uploaded (20)

Finite Element Method Linear Rectangular Element for Solving Finite Nanowire Superlattice Quantum Dot Structures