Green’s function technique is used to obtain the solution of Shredinger equation for impurity states in a quantum well (QW) under the magnetic field. Binding energy of impurity states is defined as poles of the wave function. We studied effects of the magnetic field magnitude and impurity position on the binding energy. The calculations were performed for both ground and excited states. The dependences of binding energies versus impurity position and magnetic field are presented for GaAs/Al 0.3Ga 0.7As QW.
Since Bastard calculated [
This paper is devoted to the investigations of the magnetic field effects on shallow impurity states (IS) in semiconductor quantum well. Variational techniques [1-5] are most frequently used for this aim. But it is necessary to note that the accuracy of the variational method cannot be estimated. Besides, most of variational calculations are restricted by weak-field and strong-field cases. In [
To calculate wave functions and binding energy of impurity states in a QW under magnetic field we need to solve the appropriate Shredinger equation. Then for energy gap (below the first QW level) we will receive the solution describing localized impurity states. Impurity states belonging to the second and higher QW levels are in the continuum and, therefore, resonant. In this case the solutions of Shredinger equation will correspond to a superposition of impurity and free electron states. And as a result we have band structures and wave functions modified by impurity potential [
Within the framework of an effective-mass approximation, the Hamiltonian in the presence of a magnetic field B that is perpendicular to the plane of a QW can be written as
where m^(m||) are the transverse (longitudinal) electron masses,
Coulomb interaction between an electron and the impurity ion, V(z): Square-well confinement potential, e: Unit charge, k: Permittivity, k0: Dielectric constant, z0: Position of impurity atom.
The origin of coordinates in the plane of QW coincides with the position of the impurity (the growth axis of QW is taken as the z axis). Cylindrical coordinates are used: (R, q, z) are the distances and angle describing the position of an electron.
Because of the axial symmetry of the system, angular momentum projection onto the z axis is conserved Lz = h m (m = 0, ±1, ±2,··· (Magnetic quantum number)), and their eigen functions exp(imz) determine the dependence of the unknown electron wave function versus an angle q:. Solving the one-dimensional Schrödinger equation
we get the basis for an expansion of the WF:
Substituting (2) in (1), we obtain a system of differential equations.
where.
In (3) we use dimensionless variables, introduced as units of distance (Bohr radius ab), energy (Rydberg R*) and measure of the magnetic field respectively
, ,
where—frequencies.
To construct the solutions of (3) we introduce the Green’s function GN(R, R') via the equation
Thus the expression for fN(R) can be written as
The Green’s function GN(R, R′) is given by
where I, K are the standard Bessel function of order m/2.
The Green’s function can be written directly for the left part of Equation (3). In this case it would dependent on the energy and would have different formulas for the gap and continuum (see for example [
The integrand in (4) is a decreasing function of distance R, so that for sufficiently large R the integral can be cut at some value Rmax and replaced by the finite sum
where j is the number of equidistant pieces of length DR on the interval [0, Rmax].
Expression (5) can be rewritten as
where the operator A presents the matrix, whose terms are defined as
Unlike the work [
The determinant of matrix A should be equal to zero for the existence of the solution of the system (5):
The solution of the Equation (6) allows finding the impurity binding energy for various impurity location and values of the magnetic fields.
The impurity binding energy was calculated for GaAs/ Al0.3Ga0.7As QW of width d = 125A. Such structure has been calculated by R. Chen et al. [
We used the values dR = 0.1 and Rmax = 10 in the calculations. The increase of Rmax and (or) change of dR did not change the result. In general a change in the result with a change of calculation parameters is determined by the oscillations of the Bessel functions. The argument of Bessel’s function (gR2) predicts that as a magnetic field increases n times the value of Rmax should increase times. Taking into account the fact that the values of magnetic fields generally used in the experiments are more than 1 Тesla, for which the result is steady at Rmax = 10, such situation is not seen as critical. Let us notice, that the similar situation is observed in a plane wave basis expansion method [
The comparison between the results of the present work and [
The position of the impurity states for the case of edge doping is also shown.
A well known Coulomb expression tells us about the inverse relationship between charge energy interactions and distance. Since the maximum of the electron wave function at the first QW level is in the center, then the offset of the impurity atom to the edge will increase a distance between the ion and the electron, and thus the binding energy will decrease (as it is seen in
The increase of a magnetic field increases the localization of impurity wave function (
In conclusion, Green’s function technique is used to obtain the solution of Schrödinger equation for impurity states in a square QW under magnetic field. Results of our calculations are in agreement with the data of other authors. We studied the effects of the impurity position
and magnitude of magnetic field on the binding energy. A displacement of the impurity atom from the center to the edge has a different effect on the ground and excited states. The binding energy of the ground state reduces significantly, while for the excited states it changes insignificantly. The increase of the magnetic field increases the space localization of the wave functions.
Our method does not depend on either additional parameters or physical assumptions; it allows taking into account other QW subbands, demonstrating stable and steady results at the increase or change of calculation parameters Rmax and dR.
This work was supported in part by the Russian-Ukrainian project grant 06-02-12 (U) and by grants from the project of Introducing Innovation Research Team by Guangdong Province No. CXB 2011092 10101A, CXB 2011092 10103A.