SEMICONDUCTOR QUANTUM WELL PHYSICS
16/02/1999
Overview
We are mainly concerned with (in this course) solids. Inparticular what happens when we reduce the dimensions of the structures. The physics changes dramatically as we change the number of dimensions of a solid. The dimensions can either be macroscopic (1 cm (108 Å) for example) or microscopic (10 Å). A 3d macroscopic solid means that it is 1cm*1cm*1cm. This is refered to as a bulk material. A 2d is classed as 1cm*1cm*100Å. Hence a 1d solid is classed as 1cm*100Å*100Å, and finally 0d is classed as 100Å*100Å*100Å. Remember though, dimensions here are refering to macroscopic dimensions.
A 2d structure is called a quantum well. A 1d structure is called a quantum wire. A 0d structure is called a quantum dot.
Quantum dot lasers have been built with dimensions in the order (10Å)³ and an output of 1 micro-watt. If they are all stacked together into a structure of (1cm)³. The power of this structure is 1021-6=1015 W, which is a large amount of power.
"It is win all the way with lowered dimensions!"
3D Solids
Say we have two Hydrogen atoms separated by infinity.
fig 1
Both atoms have their own (equal) energy level. If they are brought together a Hydrogen molecule is formed, this molecule has two energy levels, one is higher than the atomic energy level and the other is lower than the atomic energy. The lower is called the bonding level and the higher is called the antibonding level. This is typical energy level splitting. If a structure of 1 mole of atoms is brought together there are roughly 1023 energy levels together, the difference between these levels are tiny, the levels are called an energyband. When the gap between levels is so small we refer to it as a quasi-continuous energy band. An electron with energy kT at room temperature can easily move around quasi-continuous energy states/bands.
We need to know what the relationship is between the energy and the momentum of charge carriers such as electrons. These are two quantities which are always conserved. Recall that a particle with momentum P has a De Broglie wavelength lambda associated with it.
P=h/lambda
P=(h/lambda).(2pi/2pi)
P=2pi.hbar/lambda
=hbar.k
(k=2pi/lambda)
So if we want to connect the energy and the momentum of the partical we can equally say we want to connect the energy and the wave vector, k.
It has taken theoreticians a long time to solve this from first principals, experimentalists solved it earlier with a slightly flawed method.
Suppose an electron is moving around in free space, it has no potential energy, so all of it's energy is kinetic (it is also moving around at non-relativistic speeds). We can express the energy as,
E=½mv²
E=P²/2m
E=(hbar²/2m)k²
In a vector relation we can say,
k²=kx²+ky²+kz²
Therefore for the electron moving around in free space,
E=(hbar²/2m).(kx²+ky²+kz²)
Just consider,
ðE/ðkx=(hbar²/m).kx
ð²E/ðkx²=hbar²/m=ð²E/ðky²=ð²E/ðkz²
ð²E/ðkxðky=0=ð²E/ðkxðkz=ð²E/ðkyðkz
(We will see the significance of this later.)
We can generally say,
E=a(kx²+ky²+kz²)
Experimentalists made a guess at the solution that the above can be expanded as a power series in the wave vector. However, the series could go on forever and it is unknown how far to go in the series.
E=ao+a1kx+a2ky+a3kz+a4kx²+a5ky²+a6kz²+a7kxkz+a8kxky+a9kykz
It is known in principal why this seems okay. The momentum operator commutes with the Hamiltonian which means the momentum components are good quantum numbers (they are conserved and constant). This means that one good qantum number can be expanded with other good quantum numbers. The only problem is where to truncate the series.
The argument which led to the power series solution can be simplified by symmetry arguments.
Symmetry Arguments
Consider a cubic crystal of Zinc Sulphide, with a zinc atom in the centre and eight sulphur atoms around it forming a cube. This is the unit cell of the macroscopic crystal, this is called point group symmetry. The point group is a set of operations that sends the unit cell into itself. If the cell is observed before and after such an operation no difference is seen, the energy remains the same.
Now let us consider a rotation by pi/2.
fig 2
If we rotate the orthornomal set of x,y and z, we see that for a rotation of pi/2,
Therefore applying this argument to the wave vector of the atom and electrons in the crystal,
kz goes into kz
kx goes into ky
ky goes into -kx
Hence we write the power series as,
E=ao+a1+ky+a2(-kx)+a3k2+a4(ky)²+a5(-kx)²+a6 (kz)²+a7kykz+a8ky(-kx)+a9(-kx)kz+...
Before After
ky a2 a1
kx a1 -a2
kz a3 a3
ky² a5 a4
kx² a4 a5
kz² a6 a6
kykz a9 a7
kykx a8 -a8
kxkz a7 -a9
From this we can say that,
a4=a5=a (say)
a6=a6=b (say)
a1=a2=a3=a7=a8=a9=0
Hence we write,
E=ao+a3kz+b(kx²+ky²)+a6kz²
Then sending x into z gives us a6=a5=a4, hence we get,
E=ao+b(kx²+ky²+kz²)
This looks like a free electron! This is called the "nearly free electron model".
The majority of solids that occur in nature have symmetry operations that send x into y and x into z. This means to the mathematics that there are no linear terms left in the expression which cause the expression to behave like the free electron model.
We choose the b parameter to resemble hbar²/2m, we actually get,
b=hbar²/2m*
Where m* is the effective mass. The electron's mass does not really change, but it is used so that the b parameter looks like a known physical quantity. This confusion is the main flaw of the experimentalist's solution !
If for a crystal x into y and x into z are both symmetry operations, we get the energy as,
E=ao+b(kx²+ky²+kz²)
This is called a direct band gap semiconductor.
If x into y and x into z are not both symmetry operations we get the energy as,
E=ao+a1kx+a2ky+a3kz+a4kx²+a5ky²+a6kz²+a7kxky+a8kykz+a9kzkx
This is called an indirect band gap semicoductor.
The confusion about effective mass also comes in when they rewrite the coefficients as,
hbar²/2mxx=a4kx²
hbar²/2myy=a5ky²
hbar²/2mzz=a6kz²
hbar²/2mxy=a7kxky
hbar²/2myz=a8kykz
hbar²/2mzx=a9kzkx
These are called the "effective mass tensors". For an electron,
1/mxx=1/myy=1/mzz
Experimentally it is found that for many solids (not all) this expansion can be truncated at the second order to a good approximation. There are some solids however where this is not a good approximation.
23/02/1999
Direct and Indirect bandgaps in bulk (ie 3D) materials
For simplicity we shal neglect all three dimensions in mathematical analysis. In general here,
E=Ao+Bkx²
We look for point group symmetry, ie there is no linear terms (no terms in kx).
When we form Zinc sulphide (for example), bringing together zinc and sulphur atoms. When atoms are brought together in a solid the individual atomic energy levels are broadend out into bands and generally the bands do not overlap in energy value.
In this example, zinc has filled inner shells and two electrons in an outer shell. When the energy levels are broadened, the filled ones must still be filled. Sulphur has filled inner shells and it's outer shell is short of two electrons from being full. The outer shell (band) of the sulphur atom is lower than the outer shell of the zinc atom, so the two electrons on the outer zinc shell go to fill up the outer shell of sulphur atom, they goto a lower energy state. The recently vacated outer zinc shell is now the next empty state of the compound. In general we do not consider the filled inner shells (apart from X-ray spectroscopy). All the physics in semiconductor physics is concerned with the highest filled electron energy level and the next unfilled energy level. The highest filled level is termed the valence band and the next unoccupied level is called the conduction band. The energy difference between the valence and the conduction band is called the band gap. In atomic terms an electron traversing the band gap actually means that it has travelled from one atom in the compound to a different atom. For our example an electron moves from a sulphur atom to a zinc atom when it is excitred to the conduction band. An electron in the condction band can move from one zinc atom to the next without any (internal) potential energy penalty, this forms an electrical current in the compound.
Electrons can move in the valence band by electrons from one sulphur atom filling the gap in another sulphur atom. Exerimentalists use a model to describe this as holes moving around, entities that have a positive charge and aniliate electrons, negative charge carriers, on contact. When electrons move in the valence band, many move at once, it is difficult to consider this many electrons moving at once, if we consider the movement of gaps in the electron population, it is easier to conceptualise. The movement of electrons and gaps (holes) are different and they even move in opposite directions, it is said they have different mobilities. The symbol n is used to denote the electron density and the symbol p is used to denote the hole cencentration. Both electrons and holes contribute to the electric current (remember, both are simply electron movements !).
Going back to the valence/conduction band model. If we have an electron high in the conduction band it will sink to the bottom of the conduction band. If we then were to excite an electron to the conduction and from the valence band then we create a hole in the valence band. All of the electrons in the valence band will fall and occupy lower empty states, this creates the effect of a hole rising up the valence band.
Recalling our energy expression
E=Ao+Bkx²
This is of the same form in either band, but the coefficient is different in each band bacause each band is acriss different atoms. If we were to plot the electron energy of both atoms we would have two parabolas. In k-space, at k=0, we have the direct energy gap between the parabolas
fig 1
However in the valence band, we now consider holes, so the energy is of opposite sign so the lower parabola is turned upside down!.
fig 2
So it can be clearly seen that band gap is the separation of the minima of the energy parabolas in k-space.
Indirect band gap
Here if we have linear terms due to asymmetry, we write,
E=A+Bk+Ck²
As k goes to zero E goes to A and k goes to ±infinity so does E. To find the value of A we find the turning point of the k parabola.
dE/dk=B+2ck=0
Therefore,
k=-B/2C
If we consider shifting the parabola along the k axis, ie we change coordinate systems.
k=k'-a
E=A+Bk+Ck²
E=A+B(k'-a)+C(k'-a)²
=A+Bk-Ba+C(k'²+2ak+a²)
=A+Bk-Ba+Ck²+2Cak+Ca²
Now let,
B=2Ca
a=B/2C
Doing this now gives us,
E=A'+Ck'²=A'+Ck²
Where,
A'=A+a²C-Ba
Now, the electron's energy in both bands are now displaced from each other, the minimums do not lie along the same k value. The band gap is now said to be indirect. It is indirect because the band gap now lies along a diagonal line in the E-k space. The use of this leads us to the next section.
Photoluminescence
Many devices in modern technology are concerned with emitting or absorbing photons.
Consider a photon hitting, or colliding, with a zinc atom. An outer electron is excited and moves off to another zinc atom, it moves through the conducton band, and another photon is emitted. This is a collision, therefore energy and momentum are always conserved.
fig 3
Absorbtion is simply a special case of a collision. In the case of absorbtion, the frequency of the emitted photon happens to be zero. This is also simply the reverse of emission, the two are the same process. All processes involving photons are scattering processes in which momentum is conserved.
Recall the De Brogliè relation,
p=h/lambda
=(h/lambda)/(v/v)
Also from Plank,
E=hv
Hence the two can be combined,
p=E/c
Because of the high magnitude of c, p is approximately zero, which implies virtually zero mass for a photon.
Consider the E-k diagram for a direct band gap in figure 4. If p=hbar.k then for pure photon absorbtion p is zero so for the transition from the valence band to the conduction band, there is no change in the momentum (because non is gained from photon) so the k vector remains constant and the transition is verticle in the diagram.
fig 4
For indirect band gap semiconductors emission and asorbtion is weak. This is because there is a difference in the k vector (and hence momentum) between the minimum energies and photons cannot provide enough momentum to give the transition.
fig 5
This means that indriect band gap semiconductors are very poor for photoluminescence applications. Direct band gap semicondutors as far better.
Einstein's theoy of stimulated emission and absorbtion
If an electron is desired to be sent from a lower energy to a higher energy. If there is a certain probability of an electron absorbing a photon to goto the higher energy. The more photons that are sent to simulate the electron the higher the chance of an absorption event. After the excitation event there are N-1 photons. If we let the process run in reverse we would see the electron come down and be stimulated by the N-1 photons to emit another photon. This is true even if N=1 we would see the electron being stimulated by the vacuum. Einstein called this stimulated emission. This is fundamental concepts of lasers.
02/03/1999
Scattering rates and Final density of states
(Transition rate)
Suppose we have some initial state, i, and we goto some final state, j, by means of some scattering process. We want to know how quickly we can get from state i to state j, i.e. the scattering rate. A general example would be how quickly an electron emits a photon after a collision. Quantum mechanics tells us that the scattering is proportional to the density of final states at j. The density of final states is conventionally written as þ.
We now need to take some tecniques from Statistical mechanics. Suppose we have a box with dimensions dx,dy and dz. It has some arbitrary position in space. We want to find the probability that in x to x+dx and y to y+dy and z to z+dz you would find a 'particle' whose momentum lies between px to px+dpx, py to py+dpy and pz to pz+dpz. The answer will be some number dN.
dN=(dx.dpx/h)*(dy.dpy/h)*(dz.dpz/h)
Where h is Plank's constant.
Heisenburg says that if we confine a particle to a region dx, the following is true,
dx.dpx=>h
So the minimum value of the product dx.dpx is h. The "dx.dpx", "dy.dpy", "dz.dpz" is called a phase space. In terms of the quantities included it is called a six dimensional 3-D phase space (for bulk material). So rewriting dN (for bulk material) we have,
dN=(dx.dy.dz)*(dpx.dpy.dpz)/h³
dN=V.dpx.dpy.dpz/h³
Where V is the total volume of the system and dN' is the number per unit volume
dN'=dN/V=dpx.dpy.dpz/h³
If we were to integrate over space we could use cartesian coordinates, dx,dy,dz or polar coordinates,
dV=dxdydz
dV=4pir²dr
In the second, polar coordinates, 4pir² is the surface of a sphere and dr is the thickness of a spherical shell. We can therefore write,
dN'=
In a solid we can write,
E=p²/2m*
(Where m* is the effective mass).
We can go onto to then say,
dE=pdp/m*
pdp=m*dE
but,
p=(2m*)½E½
dN'=(4pi/h³).(2m*)½.E½.m*.dE
dN'=(4pi.2½.pi.(m*)3/2.E½.dE/h³
dN'/dE=density of states per unit volume per unit energy interval.
þ(E)=dN'/dE=4*2½pi(m*)2/3.E½/h³
þ(E) a E½
If we set off with N atoms per unit volume and bring them together from infinity we have N energy levels before we and we should have N energy levels after, but they are now broadened.
Integ (E2,E1) {þ(E)dE} = N
The density of states at the bottom of the band is low (this is bad news for photoluminescence, eg lasers), but is very large for large values of E, so the density of states at the top of the energy band is very high (this is bad news for non-radiative decay processes, ie low quantum efficiency).
Suppose we have a laser where we exicte an electron into the conduction band, creating a electron-hole pair which then recombines with the hole in the valence band thus creating a photon. The density of states in the lower levels is low though, so the radiative transition rate is low. This means that we have fewer photons per unit time. To restate, this is because the density of final states is low.
Now look at another system. Let us suppose that we have a semiconductor which has a high energy electron which impacts the valence band. This collision hits a valence electron which is then exicted into the conduction band (there is no photon created), the imapacting electron rebounds away with a lower energy. If this process is ran backwards, this is known as a Auger process. In 3D Auger processes goto final electrons with high energy which means large density of states which means it is highly likely to happen. These processes are bad news. When the high energy electrons are released, the radiative rate goes down and energy is the dissipated in the lattice.
If we now go to two dimensional structures, we have,
dN=(dxdpx/h)(dydpy/h)
dN=A.dpxdpy/h²
We also now consider a circle instead of a sphere for polar coordingates.
dN'=dpxdpy/h²=dN/A
A=2pirdr
dN'=(2pi/h²)m*.dE
þ=dN'/dE=2pim*/h²=constant
Integ (E2,E1) {þ(E)dE} = N
Since the density of states is constant, there are more states at high levels and fewer at higher energies. So the radiative rate goes up by large factors and the nonradiative (Auger) goes down by large factors. All this is gained by having a two dimensional structure rather than three dimensional. Now let us take this further, to one dimension.
In 1D,
dN'=dxdpx/h
=Ldpx/h
=dN/L=dpx/h
E=p²/2m* > px=(2m*)½E½
dE=pxdpx/m*
dpx=m*dE/px
=(m*½/2½m*)E-½.dE
dN'/dE a E-½
So we can see the pattern,
Bulk Well Wire Dot
3D 2D 1D 0D
dN'/dE a E½ E0 E-½ E-1
In Galium Arsenide semiconductors, someone called Gunn found something unusual. In general when the voltage is increased the current increases. Gunn found that this was true upto a critical voltage, but beyond this he found that the current pulsed at high voltages. In most semiconductors the charge carriers have low energy. When carriers are sat near the bottom of the conductiton band (or the 'holes' are at the top of the valence band). These carriers are said to be COOL. Gunn was dealing with HOT carrieres.
Most experiments dealt with low k (low momentum) carriers, in general recall E=Ak². However the energy could be of the form,
E=ak²+bk³+ck4
But for small values k³ and k4 are tiny. As;
k->0 E->0
k->infin E->k4->infin
To determine when these occur we find the turning points,
dE/dk=2ak-3bk²+4ck³
=k(2a-2bk+4ck²)
=0
fig 1
Cool carriers only really exist around the minimum around k=0. Gunn pushed the carriers to explore the rest of the polynomial structure. Their momentum increased but thier energy dereased. We now wish to explore the consequences of this.
From classical mechanics,
v=p/m*=hbar.k/m*
E=hbar²k²/2m*=Ak²-bk³+ck4
1/m*=(2/hbar²k²)(ak²-bk³+ck4)
=(2/hbar²)(a-bk+ck²)
This is called non-parbolicity. m* is not a constant, it varies with k and hence varies with E. So writing our velocity expression again,
v=hbar.k/m*
v=(hbar.k/m*).(k/k)(2/2)
v=(hbar²k²/2m*).(2/k)=(2/k).E
v=(2/k)(ak²-bk³+ck4)
v=2(ak-bk²+ck³)
But if k<<1, v a k
k-> v->0
k->infin v->k³->infin
dv/dk=2(a-2bk+2ck²)=0
(k=k1 k=k2)
fig 2
From Newton's law, increasing force increases momentum. Also current is number of carriers times the average current. If the velocity decreaes then so does the current. At high voltages the momentum increases and the effective mass decreases so much that the velocity decreases which in turn reduces the current, even as the votage increases. Until the Gunn effect was discovered all conductances were positive, given dI/dV it was said that differential conductances were always positive. But in such semiconductors the differential conductance was negative.
Consider the equation of continuity from quantum mechanics,
(dp/dt)+(dj/dx)=0
and Poisson's equation,
ð²Ø/dx²=-þ/E
and Ohm's law,
j=-ò.dØ/dx
Putting these together gets us,
dþ/dt=-dj/dx
If ò is a constant,
dþ/dt=-d/dx(-ò.dØ/dx)
=ò.ð²Ø/dx²
dþ/dt=-òþ/E
þ=þo.exp(-òt/E)
If ò is negative the charge fluctuation grows exponentially quickly with the time. This is what happens in the Gunn effect. Very large internal dipoles
fig 3
09/02/1999
Donors - Acceptors and Impurities
Recall that in semiconductors we have a fillied valence band and an empty conduction band. In semiconductor devices we are concerned with electrical properties. Pure semiconductors has poor electrical conductivity. An electron can only participate in a physical process if if there are unoccupied energy levels available to it nearby. For example, Sodium is a metal with filled inner shells and one electron in it's conduction level (which we shall say has has a height of 1eV), whereas two would fill it. So the outer level is only half full. At room temperature the energy is of order kT which is about 1/40 eV. So the only electrons that could ever move are the ones that within 1/40 eV of the next unoccupied energy state. Statistically that is only about 1 to 2 percent. This for an anomaly obsetrved in the classicaql era. When measuring the conductivity of Sodium, since it had a single valence electron they expected to see a current of magnitude expected when there are as many electrons as atoms, the observed current was much lower. That was because there was of about 1/40 of the electrons could actually move.
Impurities
When a donor is added we imagine simply that we have a lattice of an origonal material (like silicon) and replace one atom with a donor atom. This donor atom has an extra electron and positive charge compared to the rest of everything. Say the extra electron would normally (meaning if it was free and not seeing the donor) goto inot the conduction band, but it can then see the extra positive charge on the donor atom, the charges interact and the energy of the donor electron is lowered, this creates a lower energy band which sits just below the conduction band. The opposite is true for acceptors, they create an energy level that is just above the valence band. The energy difference between the donor level and conduction band and acceptor level and valence band is the same. It is called the binding energy is denoted by nabla, and,
kT=>nabla
Chemical potential, µ
The chemical pootential, µ, is a measure of the of the energy required to create a 'particle' in the system. At zero Kelvin the chemical potential lies half way between the uppermost occupied energy level of any significance and the next unoccupied energy level of any significance. i.e. This energy is half way between the conduction and valence level, it is known to experiementalists is known as the fermi level.
The chemical potential always lies half between the levels, but for metals this energy lies at almost the same level at the uppermost occupied levels is because the energy levels are so close that the chemical potential looks like it is at the top of the uppmost occupied level.
What happens to µ if donors are added to the system ? A donors are added the chemical potential rises because it is harder to create/add electrons when there are donor due to the electrons orbiting the donor. Likewise if acceptors are added the chemical potenital falls because there are more 'holes' which make it easier for electrons to be crveated/added.
Increasing the dopant concentration unitl the orbits start to overlap leads to a state where we say the semiconductror is degenerately doped. This means that there are so many energy levels present we have a continuum or a band. When we have a high concentration of donors (for example) we start to get an energy band below the conduction band, adding more donors we get an overlapping of donor and conduction band. Since the chemical potential lies between occupied and unoccupied levels, it will begin to lie above the donor level. In a degenerately doped n-type semiconductor the chemical potential lies in the conduction band.
Esaki Tunnel Diode
Let us first revise our p-n junction knowledge.
Consider a p type semiconductor and an n type semiconductor, of which neither are degenerately doped.
fig 1
When the two are joined together, there will be a transfer of particles unitl the chemical potential is the same.
fig 2
That process gives rise to a contact potential. This is due to the transfer of particles. As electrons are forced to move across, it gets harder for each electron to move across because they are being repulsed by the electrons that went across before them. We get a depletion zone due to this particle transfer. The depletion zone width decreases as the density of either acceptor or donor concentration increases.
If we apply a reverse bias to the junction, we force electrons from right to left the energy difference between regions increases, this can continue until it is impossible for electrons to climb the potential hill. This is covered well in most text books.
Back to Esaki!
Esaki doped a pn degenerately, so much that the depletion region was so small that there was an appreciable tunneling current
fig 3
Becuase of the degenerate doping there are holes in the valence band and electrons in the conduction band.
fig 4
The tunnel current determined the behaviour in fig 3. We actually have current decreasing with increasing voltage. This meant we have differential resistance which lead to microwave oscillators. This works for an semiconductor you can both n-type and p-type degenerately dope.
16/03/1999
Quantum Wells
So far we have discussed 3 dimensional structures, ie three large dimensions (~ 1µm), now we discuss what happens when we reduce one of the dimensions (to around 100Å).
Consider Gallium Arsenide, then we replace some of the Gaillium with some aluminium
GaAs -> Ga1-xAlxAs
It is found that the band gap is then increased with the aluminium substitution. When this material is made in practice, it is grown in sandwiches, below is how it would look in real space.
GaAs|GaAlAS|GaAs|GaAlAS|GaAs|GaAlAS
The GaAS region is called the narrow(er) band gap and the GaAlAs is called the large(r) band gap. Each region is only 100Å thick, ie it is very small and very difficult to produce. Let us consider how these fit together on an energy basis. The smaller band gap can either fit inside the larger on or not. If it does it called a type 1 quantum well and if it doesn't it is called a type 2 quantum well.
fig 1
fig 2
Figure 2 shows the energy profile for an electron traveling through a type 1 semiconductor/quantum wells. Figure two also shows us that in a type 1 structure the electrons and the holes sit in the smae semiconductor, in the same region of space.
fig 3
Figure 3 shows the energy profile for an electron traveling through a type 2 semiconductor/quantum wells. Figure three shows us that in a type 2 semiconductor the electrons and the holes are spatially separated (in real space that is).
By drawing an analogy to bulk semconductors, in energy - k spave we have direct and indirect band gaps. In the former electrons and holes can recombine (radiative decay) very easily, an in the latter it is difficult, they need phonon assistence to perform radiative decay. So we can see that type 2 semiconductors are poor radiative emittors and type 1 semiconductors are goog radiative emittors.
Multiquantum wells and Super lattices
For this topic we will consider a type 1 semiconductor. While looking at the electrons in the conduction band of the type 1, we can see that we have a series of quantum wells, wells because of the energy profile and quantum because of the small size, the spaces between are called barriers. The wells are always of quantum order of size, but the barriers are not necessarily of quantum size.
If we perform a thought experiment on the barriers size now. Let us expand the barrier to a great distance. Quantum mechanics tells us that the wave function of a particle can tunnel into the barrier (while exponentionally dedcaying). If the barrier is of such a size that the wave function of particles in each well cannot penetrate the next well, we have multiquantum wells because the wave functions in each well are isolated. If however now, the barrier is reduced (to quantum dimensions) so that wave functions from eacah begin to overlap, particles in each well interact we have what is called a super-lattice, the band gap in a super-lattice is temed a mini-band.
An n well multiquantum well is quantum wells with a great separation and an n well super lattice is wells with a small separation. You can control the properties of the miniband in a semiconductor by varying the well and the barrier width.
We wish to have a quantitive understanding we now employ quantum mechanics. Say we have a type 1 well of depth apporximately 300 mev, we will assume the dpeth is infinite, this will simplify our analysis but still demonstrate the essential physics. The well also has a width d. For the infinite barrier model, Psi=0 at the edges of the well. For solutions here we must have harmonic/stnading waves that fit within the width of the well. So the condition for standing waves is,
d=n(lambda/2)
We have our zero energy at the bottom of the well and we also have no potential energy, so we just have kinetic energy relative to the bottom of the well. So,
E=KE+PE
E=KE
E=p²/2m
p=lambda/h De Broglie
Energy=h²/(2m*lambda²)
1/lamdba = n/2d
E=(h²/2m)*(n²/4d²)
E=(h²/8md²)n²
E=nabla*n²
nabla=h²/8md²
Envelope function Approximation
This is based on the following argument. Suppose we have a bulk material, atoms separated by say one angstrom. For an electron travelling along sucha lattice, there is only appoximate solutions, based on a plane wave,
A*exp(ikx)
It is supposed that there are no atoms present (this is not a real construct), we write,
Psisolid=Aexp(ikx)*f(x)
but,
f(x->infin)=1
and since we are in a periodic lattice we can say,
f(x+a)=f(x)
This is Bloch's theory. Now consider the given hamiltonian,
hamiltonian=p²/2m + V(x)
(where m is the free electron mass). So we can then write,
(p²/2m + V(x))Psi=E Psi
(p²/2m + V(x)) exp(ikx) f(x) = E exp(ikx) f(x)
through rigorous mathematical manipulation we can eliminate,
(p²/2m* + V(x)) exp(ikx) = E exp(ikx)
Where m* is the effective electron mass. the exp(ikx) is a slowly varying fuction because k is small.
We introduce a rapidly varying function, u(x), which varies rapidly describing the potential along each atom, we write,
exp(ikx)*u(x)
Plotting this would gives the periodic function (u) damped by the slowly varying function (exp), as an envelope, in harmonic terms.
The envelope approximation assumes that an electron in a quantum well has slowly varying function which can be solved for by using the effective mass approximation.
In 3D we solve the equation,
((-hbar²/2m)*ð²/ðx² + u(x))Psi = E Psi
-((hbar²/2m*)ð²/ðx² + V(x))Ø= E Ø
In 2D we have,
((-hbar²/2m)*ð²/ðx² + u(x))Psi = E Psi
((-hbar²/2m*)*ð²/ðx² + V(x))Ø = E Ø
Where Ø is the envelope function. Now for the infinite well approximation,
((-hbar²/2m*)*ð²/ðx² + V=0)Psi = E Psi
((-hbar²/2m*)*ð²/ðx² + V1)Psi = E Psi
((-hbar²/2m*)*ð²/ðx²)Psi = (E-V1) Psi= E' Psi
ð²Psi/ðx² = -a²Psi
a²=2mE'/-hbar²
Let Psi=y
ð²y/ðx²=-a²y
The solution is,
y=Asin(ax)+Bcos(ax)
y=0 when x=0, ie B=0.
Hence,
y=Asin(ax)
y=o when x=d,
0=Asin(ad) -> ad=npi (n is integer)
1)
(2mE'/hbar²)=n²pi²
E'=(n²m²hbar²/2pid²)
=(n²pi²/2md²)(hbar²/4pi²)
=(hbar²/8md²)n²
Psi = A sin (ax) where a varies with E.
THe total wave function is PSI=u(x).sin(ax), where u(x) varies rapidly.
For electron recombining, where A is a vector potential,
p -> p-(e/c)A
p²/2m -> p²/2m - (e/mc)Ap-(e/mc).divA
p²/2m -> p²/2m - e/mc A p
-> kinetic + operator
operator a ð/ðx
Transition rate a (matrix element)
Psiinitial=Psifinal
Matrix element = Integ () {Psi*f.ð/ðx Psii}
(if any of this was required in an exam, it would be given!!!).
Say we have an electron recombining with a hole, Psii describes it in the conduction band and Psif describes it i the valence band.
Psii=ue(x)*Øe
Psif=uh(x)*Øh
u varies rapidly and Ø varies slowly. Now we define the matrix element,
matrix element = Integ () {(uh(x)*Øh)ð/ðx(ue(x)*Øe)dx}
= Integ () {(uh*(x)*Øh)*(ð/ðx(ue(x)Øe+ue(x)ð/ðx(Øe)))dx}
= Integ () {(uh*(x)*ðue/ðx*Øh(x)*Øe(x)+uh*(x)*ue*(x)*Øh(x)*ðØe/ðx)dx}
= (Integ () {uh*e/ðx*dx} * Integ (ØhØehdx) {} ) + (Integ () {uh*(x)ue(x)dx} * Integ () {Øh(x)(ðØe/ðx)dx} )
V(x)=P * Integ () {ØhØedx}
This is the expression for the internad transition, which can verified by experiment.
In the Integ () {uh*(ðØe/ðx)} term is zero and ue is a funciotn of a free electron at u=0. The various energies are orthogonal to one another and therefore is zero.
We also consider the inter
Recap
Envelope functions approximation
Bulk - Bloch's theorem,
Psi(r)=exp(ik.r)*uk(r)
uk(r+a)=uk(r)
a=lattice vector
We wish to eliminate uk, we use the Schrödinger equation,
[(-hbar²/2mo)del² + U(r)]Psi = E Psi
[(-hbar²/2mo)del² + V(r)]Ø = E Ø
Ø=exp(ikr) - Bulk
V(r)=slowrly varying potential at the bottom of a given band
=const.
Doing this is known as effective mass theory
When we bring different semiconductors together, we get a periodic potential, the difference in the potentials being VB,
((-hbar²/2m*w)del² + V=0)Øw=EØw
((-hbar²/2m*B)del² + VB)ØB=EØB.
We then have,
Psi=ØU
Where Ø is the envelope function.
13/04/1999
We are interested in the optical properties of the semiconductor combinations. We wish to determine the "allowed" transitions, ie the selection rules.
fig 1
We use something called the Transition Operator which looks something like this,
A.þ
Where A is the vector potential of the E.M. field (ie photons). þ is the linear momentum operator,
Transition rate a |Integ (Psi*f . A.þ . Psii . dT)|²
a |Integ (Psi*f.ð/ðz Psii . dT)|²
Psif=Øfuf
Psii=Øiui
ð/ðz Psii = (ðØi/ðz)ui + Øi(ðui/ðz)
Integ (Psi*f.ð/ðzPsii) = Integ (Ø*fu*f)(ðui/ðz)uidz)
=Integ ((Ø*fu*f)Øi(ðui/ðz)dz)
Integ (f(z)g(z)dz)=Integ(f(z)dz)*Integ(g(z)dz)
If g(Z) is rapidsly varying f(z) is slowly varying
= Integ (Ø*f(ðØi/ðz)dz) * Integ (u*f.ui . dz) + Integ(Ø*fØidz) * Integ(u*f(ðui/ðz)dz)
Infra Red detectors or Infra red lasers.
Because the structure of such materials is periodic initial and final states of functions become equal (for inter sub band transitions), so we have the general quantum mechanics results that certina integral factors are equal to 1. So we come to the conclusion,
Transition Rate (TR) a Integ(Ø*f(ð/ðz)Øi)
Also recall,
Øeven(-Z)=Øeven(Z)
Øodd(-Z)=-Øodd(z)
fig 2
From the argument in figure two, the rule of thumb is that the allowed transitions are even to odd (odd and even energy levels that is).
Interband Transitions
Integ(u*fuidZ)=0 because uf and ui are eigen functions belonging to different energy levels
Integ (u*f(ð/ðz)ui) = some number.
Integ (Ø*fØidz)
Here odd to even transitions become zero. Even to Even transitions give some value and hence exist. Odd to Odd transions give a value also. This means that for interband transitions only odd to even transitions are phorbidden.
Quantum Well lasers - Solid State
Photons stimulate absorption and emission (Einstien). If we send in photons that exactly match the energy gap between levels. The stimulated absorption is proportional to the number of photons, N.
fig 3
fig 4
Einstein thought that emission was simply the reverse of absorption. However going through the N and N-1 photon levels, when N=1, we have spontaneous emission.
fig 5
Radiation - þ(v) this is energy density per unit frequency.
Einstien suppose he had a two level atom.
fig 6
The stimulated absorption rate is W12=B12þ(v)
The stimulated emission rate is W21=B21þ(v)
The spontaneous emission (down) rate is A21
down rate = N2(B21þ(v)+A21)
up rate = N1B12þ(v)
At equilibrium up rate=down rate,
N2(B21þ(v)+A21) = N1B12þ(v)
Einstein guessed B21, A21 and B12 are properties of the system, they do not depend on the radiation. If they do not depend on the radiation field choose a particular radiation field, in particular, black body radiation. This is because Plank had just determined the þ(v).
þ(v)=(8pihv²/c³)(1/(exp(hv/kT)-1))
This expression is for a system in thermal equilibrium. In which case Boltzmann statistics applies, which state the density of energy states is,
þ(E) a exp(-ßE)
ß=1/kT
Suppose we have density of states g2 and g1.
N2 a g2exp(-ßE2)
N1 a g1exp(-ßE1)
þ(v)=(8pihv²/c³)(1/(exp(hv/kT)-1))
N2/N1 = (g2/g1)exp(-ß(E2-E1))
So we can write,
B21þ(v)+A21=(N1/N2)ß12þ(v)
A21=((N1/N2)B12-B21)þ(v)
A21/((N1/N2)B12-B21)=þ(v)
A21/((g1/g2)B12exp(ßhv)-B21)=þ(v)
(8pihv²/c³)(1/(exp(hv/kT)-1))=A21/((g1/g2)B12exp(ßhv)-B21)
Therefore,
(g1/g2)B12=B21
þ(v)=A21/(B21(exp(ßhv)-1)
Therefore,
A21/B21=8pihv³/c³
Therefore,
A21=(8pihv³/c³)B21
Remember, c=lambda*v, hence,
A21=(8pih/lambda³)B21
For microwaves the spontaneous emission is negligible. For visible spectra, B21þ(v)>>A21 for the lazing state.
20/04/1999
Now we rewrite,
W21=B21þ(v)=(A21c³/8pihv³)þ(v)
A21 is known as the spontaneous emission rate, we write,
1/Tspont=A21
Experimentalists would say Tspont is the spontaneous radiative decay time.
W21=c³þ(v)/8pihv³Tspont
This is written with an additional term, with the line .shape function,
W21=(c³þ(v)/8pihv³Tspont)a(v)
Where,
Integ (+infin,-infin) {a(v)dv} = 1
Threshold conduction in a LASER
Consider two energy levels with N2 electrons and N1 electrons in them, (N2 in the upper course). The following relation can be used,
c.þ(v)=I(v)
Where I is the intensity of the beam. þ is the density of the beam, so in unit time c.þ is the energy which has flowed through our system.
Hence,
W21=c²I(v)a(v)/8pihv³Tspont
Now consider,
dI/dz=(W21N2-W12N1)hv
dI/dz=(N2-N1(W12/W21))W21hv
W12/W21=þ(v)ß12/þ(v)ß21
W12/W21=ß12/ß21
=g2/g1
dI/dz=(N2-(g2/g1)N1)W•21hv
dI/dz=(N2(g2/g1)N1)(c²I(v)a(v)/8pihv³Tspont)I(v)
dI/dT=y(v)I(v)
y(v)='gain'=(N2(g2/g1)N1)(c²I(v)a(v)/8pihv³Tspont)
dI/dz=yI
dI/I=ydz
Ln(I)=yz+c
I=exp(yz)exp(c)=Ioexp(yz)
We need to prove that g2=g1. If g2=g1,
y=N2(g2/g1)N1)(c²I(v)a(v)/8pihv³Tspont
Therefore,
y=>0 if N2=>N1
This is population inversion. Recall our 3D relation,
dN=(dx/h)(dy/h)(dz/h)dpx.dpy.dpz
a dpx.dpy.dpz
a 4pi.p².dp
p²=px²+py²+pz²
p=h/lambda=hbar.k
dN a k².dk
Solid State LASERs
Population inversion is eachived by passing an electric current through the system. First let us look at a bulk 3D solid state laser.
Take a semiconductor, one side n-type the other side p-type, and degenerately dope it. This means that energy levels will spread out,
fig 1
From this we would expect at best to get a photon emitted when the electron transits the band gap, this will not always happen, so it is not 100% efficient. Such lasers are roughly 70% efficient.
fig 2
Solid state lasers are important because they operate at voltages compatible with integrated circuits and can be monolithically integrated with them using photo-lithographic processes used in integrated ciruits.
Next we return to an earlier step, we need to determine N1 and N2. Consider the depletion zone of a semiconductor. It has a thickness t and cross-sectional area A. Suppose we put in a total number of N dopants. So the number of dopants per unit volume is,
N/At
If the lifetime of (minority) charge carriers is T, then it can be shown that,
t=(DT)½
Where D is a diffusion constant. The charge carriers obey Fermi-Dirac statistics, so we can write,
N2=(N/At)fc(E2)
N1=(N/At)fv(E1)
Where f is a fermi-dirac function for the conduiction and valence bands, indicated by the subscript.
So recal that,
y(v)=(N2(g2/g1)N1)(c²I(v)a(v)/8pihv³Tspont)
y=(N/At)(fc(E2)-fv(E1))(c²I(v)a(v)/8pihv³Tspont)
We now need to change this into a form which has experimental measurements. Let us consider the electric current, I. I=number of electrons per unit time multipllied by charge. Hence, I/e=number of electrons per unit time flowing through the system. For semiconductors this is the number of electrons that cross the bandgap per unit time. If N is the total number of electrons, then I/e is also N time the probability that an electron crosses the band gap. So we can write,
I/e=N(PR+PNR)
Where PR is the radiative rate and PNR is the non-raiavte rate (Auger, multi-phonon emission). We now write,
I/e=N(PR+PNR)(PR/PR)
=N((PR+PNR)/PR)(1/Tspont)
N/Tspont=(I/e).(PR/(PR+PNR))=(I/e)np
Where np is the radiative eifficeny of the system.
np=PR/(PR+PNR)
y=(1/At)(fc(v)+fv(v))(c²a(v)/8piv²e)(I/e)
y=(1/At)(fc(v)+fv(v))(c²a(v)/8piv²)I
(fc(v)+fv(v)) is callled the fill factor (it's maximum value is one). So ythreshold is required to before lasing begins.
Ithreshold a At/(fc(v)+fv(v)).((1/n)(1/a))
In 3D structures, A is in mm² and t is 10000 Å. In 2D structures A is in mm² and t is 100 Å.
Now let us consider the threshold conditions. If we have a depletion zone of width 1 µm. The wavelength of light is often larger than this, so we would have a source of loss by means of the n and p regions, we also have gain in the deletion zone. Now consider that if we have two mirros Seperated by a space L, and photons relfect between the two and amplify each other (imagine the depletion zone barriers are like mirrors!)
LGamma.yth=aa.LGamma+ac(1-LGamma)+(1/L)ln(1/e)
LGamma is the fraction of lasing mode confined to the active region. The Ln term is the distributed mirror loss term. aa and ac are the optical losses,
R=(R1R2)½
R1 is the reflectivity at mirror one and R2 is the reflectivity at mirror 2.
Now consider a wave of unit power, which travels a distance 2L (return trip from one mirror the other). When it returns it's power will be,
1*R1*R2*exp(g-a)2L
g is the gain and a is the loss, both per unit length. Threshold condition is when,
1=R1*R2*exp((g-a)2L)
Now take logs of both sides,
0=ln(R1R2)+2L(g-a)
now write,
g=a-(1/2L).Ln(R1R2)
also note,
Ln(1/R1R2)=ln(1)-ln(R1R2)=-Ln(R1R2)
g=a+(1/L).Ln((1/R1R2)½)
=a+(1/L).Ln(1/(R1R2)½)
So recall from above,
LGamma.yth=LGamma.aa+(1-LGamma)ac+(1/L).Ln(1/(R1R2)½)
ac is the loss in the cladding layer and aa is the loss in the active region.
Now we wish to find out what the power output is.
P=((I-Ith)/e)hv.n
P=((I-Ith)/e)hv.n.((1/L).Ln(R)/aaLGamma+acLGamma+(1/L).Ln(1/R))
27/04/1999
Ultra low pressure
This involves sending atoms onto a substrate. The Mean free path between collisions is large (this is reffered to as molecular flow). The general process this falls under is called Molecular Beam Epitaxy (MBE). THe processes in MBE are well understood. MBE allows us to grow monolayer by monolayer, which allows us to make atomicallu abbrupt inferfaces, both of these are advantagous. The monolayer transition region is small (such regions are termed interface roughness). The disadvantage with this is that is very expensive. Another disadvantage is that the growth rates are very slow, typically one Angstrom every 10 seconds. See fig 1 for how different pressures are used for different materials. Because of the nature of this process, there is no boundry layer formed.
fig 1
Vapour Phase Epitaxy
A mixture of gases flow across a substrate and interact with one another. The flow is viscous, and a boundry layer is formed. The chemical process involved are poorly understood. The interfaces are not (in general) as atomically abbrupt as in MBE. The main advantage is that it is very cheap in relation to MBE.
Now consider we were to grow a super lattice, with a potential as shown in fig 2
fig 2
The period of the structure is d (~= 100Å). The period of bulk-interatomic spacing is (1Å). The energy bands are split into what are known as Brillooin zones. The Brillooin zone in the bulk is written as,
-pi/a<=k<=pi/a
So for our superlattice,
-pi/d<=k<=pi/d
fig 3
Isaaki argued the following, a charge carrier can get trapped in a Brillooin zone, on the E-K curve of figure 3. It's momentum is then periodically reveresed/reflected. This will give us an AC current from applying a DC voltage. This has not been obersvered yet, but the search continues ! Such a trapped oscillating charge, is known as a Bloch Oscillator.
Consider the energy relation,
E=hbar².k²/2m
dE/dk=hbar²k/m=hbar(hbar.k/m)=hbar.(p/m)=hbar.v
v=(1/hbar)(dE/dk)
F.e=dp/dt
(field strength)
=hbar(dk/dt)
eF/hbar=dk/dt
Acc'n. = dv/dt =dv/dk * dk/dt
v= (1/hbar) (dE/dk)
dk/dt=eF/hbar
a=dv/dt=(dv/dk)(dk/dt)
=(1/hbar²)eF(d²E/dk²)
At the point of inflection on our E-k plot, as the momentum increases the acceleration reduces, that indicates negative differential resistence! This was another prediction from Isaaki that has not been realised yet.
Douisle Barrier Tunneling Diodes
Consider a conduction band function, with two barriers, the barriers are 10 Å wide and the spacing between them is 50 - 100 Å. To either side of the pair of barriers, we degenerately dope the substance.
fig 4
The electrons still cannot get past the barriers. Suppose we have an electron trapped between the barriers and we make the barriers infintely high. If the electron makes a standing wave, we have,
d = lambda / 2
lambda = 2d
p=h/lambda
E=p²/2m=h²/2mlambda²
=h²/8md²
This is called a resonace state, since the electron makes a resonant curve in the well.
For electrons in a bulk metal structure, we have the fermi energy,
Ef=hbar²kf²/2m
If we have the same barriers sturcture in our metal as described above. There is a resonant energy level in between the barriers. If we apply a voltage we rise the fermi level, once the fermi level is at the same height as the resonant level, they can tunnel through the barrier. On the outside of the barriers we have a 3D system and we have a fermi-sphere of electrons which have the same energy
fig 5
n-i-p-i Structures
fig 6
This is simply structure doped so thaat we have n type, an insulator, p type and another insulator. These structure give rise to Type II superlatticies. There use is good for non-linear optics. Suppose that the concentration of n types is the same as that for p types, this is called self-compensated. We have ionized donors and ionizes acceptors. Consider the potential profile.
fig 7
If we have a square well and choose our origin to be in the centre a charge carrier would see the same potential either side. We could also say for at that point,
dV/dx=0 (at x=0)
If we then solve for Poisson's equation,
ð²V/ðx²=4pie²n/E
n is the dopant concentration.
ðV/ðx=(4pie²n/E)x + C
ðV/ðx=0 at x=0 therefore C=0
ðV/ðx=(4pie²n/E)x
V=(4pie²nx²/2E)+C2
V=0 at x=0 hence C2=0
We now work out what V is at x=d/2 (d is the width of the well).
V=(4pie²n/2E)(d/2)²
The depth of the well depends on n and d.
the end (I think)