The energy levels of bound states of an electron in a quantum well with BenDaniel-Duke boundary condition are studied. Analytic, explicit, simple, and accurate formulae have been obtained for the ground state and the first excited state. In our approach, the exact, transcendental eigenvalues equations were replaced with approximate, tractable, algebraic equations, using algebraic approximations for some trigonometric functions. Our method can be applied to both type I and type II semiconductors and easily extended to quantum dots. The same approach was used for the semi-quantitative analyze of two toy models of Janus nanorods.
- type I and type II semiconductors
- BenDaniel-Duke boundary conditions
- Janus nanorods
- toy models
In the last three decades, nanophysics became a domain of increasing interest and intense research, due to the huge number of new effects produced at nanoscale level, in quantum wells (QWs), quantum dots (QDs), Janus nanoparticles, etc. These new effects are fascinating from the perspective of both applied and theoretical physics. The semiconductors provide the largest area of challenging subjects, due to their applications in nanoelectronic devices, multifunctional catalysis, (bio-)chemical sensors, data storage, solar energy conversion, etc.
An attractive aspect of nanophysics is the fact that a quite large number of interesting problems can be approached using quite simple theoretical tools, sometimes at the level of one-particle quantum mechanics. In some cases, the properties of nanostructures like quantum wells, quantum dots, or quantum rods can be explained by just solving the Schrodinger equation with simple potentials. For instance, the basic physical properties of a heterostructure consisting of a thin layer of a semiconductor sandwiched between two somewhat larger semiconductors of identical composition, can be obtained from the study of the movement of a particle with position-dependent mass (PDM) in a finite square well. This particle is, of course, a charge carrier in the semiconductor, and—in our case—will be an electron. As the effective mass of a charge carrier in a semiconductor depends on the charge carrier-lattice interaction, it changes if the lattice composition or the symmetry changes. So, excepting the case of a charge moving in a perfect crystal, the effective mass of an electron or hole is, rigorously speaking, position dependent.
But, simple as the theoretical tools needed for its investigation are, this problem of quantum mechanics involves two important issues: the position-dependent mass (PDM) quantum physics and semiconductor heterostructures. Let us shortly comment on these points.
The roots of the position-dependent effective mass concept are to be found in the pioneering works of Wannier (1937) and Slater (1949) (see Ref. 1 in ). Recent papers give explicit methods to obtain explicit solutions of the Schrödinger equation with PDM, for various forms of this dependence and for several classes of potentials [2, 3, 4].
However, in practical situations usually encountered in the physics of semiconductor junctions of two materials, and the simplest and more popular form of position dependence of the effective mass is a step function: the effective mass has a constant value in the material and another, constant value, in the material In such a case, the most convenient approach for obtaining the wave functions or the envelope functions in a heterostructure—for instance, a quantum well (QW) or quantum dot (QD)—is to solve the Schrödinger equation with BenDaniel-Duke boundary conditions for the wave function [5, 6].
The transition from the complex problem of a real semiconductor (for instance, Kane theory) to the simple problem of a particle moving in a square well with BenDaniel-Duke boundary conditions is indicated, for instance, in Chapter III of Bastard’s book . This simple problem provides, however, a realistic description of states near the high-symmetry points in the Brillouin zone of a large class of semiconductors. “It [i.e., ‘the simple problem’] often leads to analytical results and leaves the user with the feeling that he can trace back, in a relatively transparent way, the physical origin of the numerical results.” (, p. 63).
The boundary conditions for the wave functions or envelope functions at interfaces generate the eigenvalue equations for energy; of course, different boundary conditions generate different eigenvalue equations. They are transcendental equations, involving algebraic, trigonometric, hyperbolic, or even more complicated functions. With few exceptions (for instance, the Lambert equation ), their solutions, which cannot be expressed as a finite combinations of elementary functions, are not systematically studied.
However, in some situations, quite accurate analytical approximate solutions can be obtained. When a transcendental equation mixes algebraic and trigonometric functions, it might be possible to approximate the trigonometric functions with algebraic expression, and to transform, in this way, the exact transcendental equation into an approximate algebraic one. In its simplest form, for instance, in approximations like , for this “algebraization” is largely used. But what is really interesting is to use algebraic approximations of the trigonometric functions valid on their whole domain of definition, as de Alcantara Bonfin and Griffiths proposed in a recent paper ; such analytical approximations have been studied and extended by other authors .
In this chapter, we shall obtain approximate analytical results for the energy of electronic bound states in quantum wells and in simple models of Janus semiconductor nanorods. As the concept of Janus nanoparticle is less popular than the concept of QWs or QDs, we shall give here some short explanations.
Their name derives from the Roman god Janus: his head had two opposite faces. A Janus nanodot can be a sphere composed of two semispheres of different materials. A Janus nanorod can be a nanorod having the left half and the right half made of different materials. Due to their intrinsic duality, the opposite parts of Janus particles can be functionalized differently . Janus particles with an electron-donor and -acceptor side may be used in photovoltaics. As the Janus nanoparticles have lower symmetry than their homogenous counterparts, their theoretical description is more difficult. In this chapter, we shall propose toy models for semiconductor Janus nanorods.
The structure of this chapter is the following. We shall firstly formulate the basic theory for the quantum mechanical problem of a quantum well, composed of a thin semiconductor sandwiched between two massive ones. This heterojunction can de-modeled by a quantum well (QW), essentially a finite square well, with BenDaniel-Duke boundary conditions. Such a problem was recently discussed by several authors, like Singh et al. [11, 12], who replaced the trigonometric functions entering in the transcendental equations for the bound states energy by the first few terms of their series expansion; in this way, the equations become simple, tractable algebraic ones. Our approach is different, being based on a more sophisticate “algebraization” of trigonometric functions, as proposed by de Alcantara Bonfim and Griffiths . We shall obtain explicit formulas (series expansions) for the ground state energy and for the first excited state, very accurate if the well is not too shallow. Our results can be applied to both type I and type II semiconductors.
In the last part of our chapter, we shall study two toy models for semiconductor Janus nanorods; for the simplest one, we shall obtain analytical expressions for some energy eigenvalues of electronic bound states.
2. Basic theory
We shall solve the Schrodinger equation for an electron moving in a square well, described by the potential:
considering that its mass is position dependent. More exactly, the mass inside the well, , and the mass outside the well, are different:
So, the Schrodinger equation for bound states is:
Its physically acceptable solutions, that is, the wave functions, have to satisfy two conditions: (1) the continuity of the wave function and (2) the continuity of the probability currents density at the interface. The first one is encountered in all quantum mechanical problems, but the second one is specific to the case of the position-dependent mass , defined by the Eq. (2), and takes the form:
Eq. (4) is known as the BenDaniel-Duke boundary condition. The notations were used here to make more visible the physical content of this special boundary condition, and will not be maintained in the rest of the chapter.
The bound state has a unique energy, , but two wave vectors, one inside the well, , and another one outside, :
Due to the parity of the potential, the wave functions can be chosen to be symmetric or antisymmetric.
The symmetric wave functions, describing the even states, are:
The ground state wave function is, of course, . The antisymmetric wave functions, describing the odd states, are:
The continuity of these functions in gives:
So, the wave function outside the well is:
The wave functions are normalized if:
It is convenient to use the potential strength (introduced by Pitkanen , who actually used instead of )
and to define also , and as:
, and are dimensionless quantities; will be sometimes called dimensionless wave vector.
It is easy to see that:
If the mass is position dependent, according to (2), the eigenvalue equations obtained from the Schrodinger equations, using BenDaniel-Duke boundary conditions have the form:
We shall consider that both are positive; this corresponds to type I semiconductors. So, with , with replaced by for even states and by for odd states, we can put the Eqs. (25) and (26) in a more convenient form:
For they take the form of the well-known equations for the energy eigenvalues of the finite square well. Approximate analytical solutions of these equations were obtained for deep wells  and in the general case [8, 9, 16, 17].
In this chapter, we shall obtain precise analytical approximations for the energy of the first two states, that is, for the ground state and for the first excited state, considering the cases and separately. For moderate and deep wells, the formulae are both simple and accurate. In the limit we shall obtain the result of de Alcantara Bonfim and Griffiths, Eq. (17) of .
3. Approximate analytical solutions for eigenvalue equations
3.1. The first even state (the ground state)
According to Eq. (29), the dimensionless momentum of the first even state, which is also the ground state, is the smallest positive root of the equation:
We shall discuss separately the cases and
3.1.1. The case
It is useful to introduce the new parameters :
because the eigenvalue equation can be written in a simpler form:
In the most physically interesting cases, is quite large (the wells are quite deep), and according to (32), is even larger, so it is more convenient to use instead of as “large parameter”.
The precision of this approximation on various subintervals of depends on the exact value of with in our numerical evaluation, we shall use the value For a detailed discussion on this issue, see .
The algebraic approximation of the eigenvalue equation, we get with (34) is:
(35) can be written as:
and considering that the well is not too shallow:
we obtain for the physically interesting root the expression:
If the depth of the well increases indefinitely, and as requested. Indeed, in a finite well, the energy of a bound state is smaller than the corresponding energy in an infinite one, so the first term in in the previous formula is negative.
It is useful to write (39) in terms of more physical parameters, and In order to do this, let us notice that:
so Eq. (39) takes the form:
If the parameter cannot be considered “large,” the exact expression of the root can be obtained using the standard approach ; they are elementary, but cumbersome, and will be not given here; the interested reader can find them in  .
3.1.2. The case
If , the eigenvalue equation for the dimensionless wave vector is:
with the following definitions for the parameters:
a cubic equation:
Following the same steps as in the previous case, we find that the parameters enter into the various expressions needed for obtaining the cubic roots only through the monoms at various powers, and the roots of Eq. (46) can be obtained from the root (39) making the substitution:
3.2. The first odd state
3.2.1. The case
The exact eigenvalue equation for the first odd state, which is also the first excited state, can be written as:
As the shape of the function on the interval is quite similar with the shape of on the interval we can try an algebraization for similar to that proposed by de Alcantara Bonfim and Griffiths for
A detailed discussion of the precision of this approximation is given in  (see Fig. 3 and Eq. (88)). Following, exactly the same steps as in the case of the ground state, we find that
For both cases——in the limit of an infinitely deep root, as requested, and the first correction to this value is negative.
The relative errors of the formulas (39) and (50), with respect to the exact roots of the corresponding algebraic equations, are very small—of about for physically interesting values of the parameters In other words, the main contribution to the errors of our results is given by the approximation of trigonometric functions with algebraic ones, not by the approximation of the exact formulae of the roots of cubic equations with the low order terms of their series expansions.
As already mentioned, one of the physical motivations of the calculation of the energy of bound states in heterostructures is to explain their photoluminescence properties. In several cases (see for instance ), the authors use Barker’s formula for the energy levels in a square well . Much more precise analytical expressions for these energy are available in the literature [8, 9], for the case of constant mass; in this paper, we propose similar formulas, considering the case of position-dependent mass.
3.3. Higher-order states
In the previous subsections, we analyzed the ground state and the first excited state of a square well, with BenDaniel-Duke boundary conditions. For the de Alcantara Bonfim formula (34) can be extended to larger arguments:
but the eigenvalue equation, obtained in this way, is a sextic equation (which cannot be reduced to a cubic equation in ), so it cannot be solved. We meet similar difficulties if we try to use in the eigenvalue equation the algebraization of (see later on, Eq. (73) and (74) of the present paper). Even the “parabolic approximation” for or , in the sense used in , gives a quartic equation for the dimensionless wave vector. Its roots are given by complicated, but still elementary formulas.
3.4. Graphical illustration of our main results
where is the root of the cubic equations obtained after the algebraization of the transcendental eigenvalue equations for the ground state and for the first excited state According to the Eqs. (39) and (50), for a deep well, the root can be approximated with a quartic polynomial in , the inverse of the potential strength Let us mention that, if we replace in the definition of Eq. (16), with the free electron mass, we choose the length of the well and we express the potential in electron volts, we get:
We shall plot our main results, that is, the series expansions of the dimensionless wave vectors, and as functions of on the range when the condition of convergence is satisfactorily fulfilled (Figure 1). The energy is a monotonically increasing function of its values, for are obtained from Eqs. (42) and (50).
3.5. Applications to other nanostructures
Our calculations can be easily applied to type II semiconductors heterostructures, when one of the effective mass of the charge carrier is negative: (, chapter 3, Eqs. (35) and (36)); a detailed description of such heterointerfaces can be found for instance in , p. 66. So, instead of (25) and (26), the eigenvalue equations take the form:
and can be solved following exactly the same approach.
As already mentioned, the wave function in the Schrodinger Eq. (3) can be interpreted as an envelope function. This approximation works well when the materials constituting the heterostructures are perfectly lattice-matched and they crystallize in the same crystallographic structure (in the most cases, the zinc blend structure). Its application is restricted to the vicinity of the high-symmetry points in the host’s Brillouin zone Actually, most of the heterostructures’ energy levels relevant to actual devices are relatively closed to a symmetry point in the host’s Brillouin zone. A popular example is given by the lowest conduction states of GaAs-GaAlAs heterostructures with GaAs layer (typically, its thickness is about 100 Å or larger). A detailed description of the cases in which the envelope function model is successful is given in , p. 66 (See Figure 2).
As there are some similarities between QWs and QDs, our results are also relevant for these devices. The simplest remark is that the eigenvalues equations for the first odd state in a QW are identical to that corresponding to the state in a QD (see for instance , problem 63). Also, the eigenvalue equations for the wave vectors of the energy levels for a finite barrier rectangular shaped QD, Eq. (36) in , are quite similar to ours—(29) and (30), but somewhat more complicated. The ground state energy of electrons and holes in a core/shell QD is given by the Eq. (21) of , an equation similar to ours, just mentioned previously. Such results are important,
4. The infinite square well with two semiconductor slabs
4.1. The symmetric case
Let us consider an infinite 1D square well, delimited by two rigid walls situated in respectively containing two semiconductor slabs, of equal width, but of different materials. It is a toy model for a Janus nanorod, composed of two different semiconductors, with large work functions. We preferred to choose this particular case (equal width), in order to avoid too cumbersome mathematical calculations. The electron effective mass is position dependent, like in (2):
We want to investigate how the energies of the electronic bound states will be affected, compared to the situation when in the infinite well there is only one slab, with effective electron mass or As
we have, with (57):
The electronic wave function is obtained solving the Schrodinger equation, as in the case of a finite well, studied in Section 2:
The boundary conditions for the wave function give:
and the continuity in the origin:
The BenDaniel-Duke boundary condition means:
Together with the orthonormality condition for the wave function, Eqs. (61)–(63) and (65) form a system of five equations for five quantities, As the amplitudes are not of primary interest, we can combine (64) and (65) to obtain:
it can be written as:
If (67) gives:
so the solutions corresponding to the infinite well with an homogenous medium inside the walls.
Eq. (70) is a transcendental one, and its solutions cannot be expressed as a finite combination of elementary functions. A quite popular analytical approximation for the tangent function has been proposed by de Alcantara-Bonfim  and generalized by the present author :
In order to see how this approximation works, let us consider the first two roots of Eq. (68), if For we obtain (for instance, using the command FindRoot in Mathematica) and close to , respectively , that is, to the values corresponding to We shall discuss the case of the second root of Eq. (70). As and the two tangent functions appearing in Eq. (70) will be approximated by the two variants of Eqs. (73) and (74), the result being the following:
So, can be obtained as a root of a second order equation, namely:
We find that so quite close to the exact value, the error being:
However, due to the rapid variation of the tangent functions near its singularities, this approximation method must be used with utmost care, as it can easily give unacceptable results (this is the case of the first root, for ).
4.2. The asymmetric case
Let us consider now the case of a rectangular infinite asymmetric well, with the potential:
with containing, as in the previous example, two semiconductor slabs. It is also a toy model of a Janus nanorod, somewhat more realistic than that discussed in Section 4.1. We also chose a particular geometry (the same width for each slab) to avoid irrelevant mathematical complications. For an electronic bound state of energy the wave vectors (and the electronic effective masses) are different in different slabs, according to the relations:
Defining the wave vector by:
and following exactly the same steps as in the symmetric case, we obtain the following eigenvalue equation:
with the notations:
If (and, evidently, ) become imaginary, and The Eq. (82) and its hyperbolic counterpart are much complicated than (68); even if there are some methods of obtaining approximate analytical solutions, they will be not discussed here. The case of a finite asymmetric well, with two different semiconductor slices, can be studied following exactly the same approach, but now the complications are even more serious, as the wave function extends outside the wells.
In this chapter, we obtained approximate analytical solutions for the eigenvalue equation of the first two bound states in a semiconductor quantum well, in a particular case of position-dependent mass of the charge carrier—in fact, the simplest one, corresponding to BenDaniel-Duke boundary conditions. This position dependence can be characterized by the ratio of the mass inside, to the mass outside the well. Actually, we obtained quite simple expressions for the dimensionless wave vector, in terms of the potential strength and of Even if we solved this problem in terms of one-particle quantum mechanics, obtaining the wave function and the eigenvalues of the bound states, our results can be directly applied in the theory of envelope functions in the conduction band at heterointerfaces. Our approach is based on the “algebraization” of trigonometric functions present in the transcendental eigenvalue equations; in this way, they are transformed in tractable algebraic equations.
We also proposed two models for a semiconductor Janus nanorod—a system, which was not yet treated analytically.
Our results can be easily extended to more realistic (e.g., linear) position dependence of the mass carrier and to other nanosystems. For instance, the eigenvalue equations for the wave vectors of bound energy levels of a finite barrier rectangular-shaped quantum dot, Eq. (36) in , are quite similar to ours—(22), (23), but somewhat more complicated. The ground state energy of electrons and holes in a core/shell quantum dot is given by Eq. (21) of , an equation similar to ours, just mentioned previously. Such results are important, inter alia, for the interpretation of photoluminescence spectra of heterojunctions.
The author acknowledges the financial support of the IFIN-HH-ANCSI project PN 16 42 01 01/2016 and of the IFIN-HH-JINR grant 04-4-1121-2015/2017.