Since Bastard calculated [1] the impurity binding energy in QW for the very first time many works on this subject have been published. Their main purpose is to determine the binding energy of impurity states. By changing this important parameter, for instance by external fields, it is possible to control the free carrier’s concentration, and thus affect many physical properties of the structure: low-temperature behavior of the kinetic characteristics, the type of conductivity, electronic mobility, optical properties.

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 [6] the numerical method is used for finding the solution for wave function (WF) as expansion on plane wave basis. Then the problem is reduced to a task of determination of eigen vectors (envelope functions) and eigen values (binding energy) of rather large matrices. An alternative nonvariational method to calculate the binding energy of impurity states in a QW was developed in [7]. The electron wave function was constructed through the Green’s function of corresponding Schredinger equation. Similar scheme is used in the present paper. Our method allowed to receive the binding energy and wave functions of IS in a QW under magnetic field of arbitrary strength. Also, the presented method led to recalculation of results for QW levels in a magnetic field. We used our method for various impurity location and values of magnetic fields in GaAs/Al_0.3Ga_0.7As QW.

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 [8]. Positions of impurity states, which are on the background of the continuum, can be defined as the poles of the WF. This approach is implemented in the present work.

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, k₀: Dielectric constant, z₀: 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 L_z = 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 a_b), energy (Rydberg R^*) and measure of the magnetic field respectively

, ,

where—frequencies.

To construct the solutions of (3) we introduce the Green’s function G_N(R, R') via the equation

Thus the expression for f_N(R) can be written as

The Green’s function G_N(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 [7]). From this point of view the above presented expressions have their own advantages: The Green function has an identical kind for all energy range; by equating the impurity potential to zero we can receive the solution for magnetic subbands.

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 R_max and replaced by the finite sum

where j is the number of equidistant pieces of length DR on the interval [0, R_max].

Expression (5) can be rewritten as

where the operator A presents the matrix, whose terms are defined as

Unlike the work [6] energy here is not eigenvalue of the matrix, but is a parameter on which the matrix elements depend.

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/ Al_0.3Ga_0.7As QW of width d = 125A. Such structure has been calculated by R. Chen et al. [5]. This means that we can compare our results with those obtained in that work.

We used the values dR = 0.1 and R_max = 10 in the calculations. The increase of R_max 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 (gR²) predicts that as a magnetic field increases n times the value of R_max 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 R_max = 10, such situation is not seen as critical. Let us notice, that the similar situation is observed in a plane wave basis expansion method [6], however this method operates with a larger matrix, whose size is 20 times more than in the present works. Therefore, the presented method will be more effective for the cases of small magnetic fields.

The comparison between the results of the present work and [5] is shown in Figure 1. As it is seen they are in qualitative and quantitative agreement.

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 Figure 1 for the states m = 0). However, the binding energy of the excited states is reduced slightly when impurity atom move to the edge. The thing is that the maximum of the wave functions for excited states is shifted from the center both along the z-axis and the xy-plane. Therefore the displacement of the impurity atom only along the z axis, within the width of QW, has no essential effect on the binding energy.

The increase of a magnetic field increases the localization of impurity wave function (Figure 2). The effect is negligible for the ground state, and clearly exhibits itself for the excited state. This is not surprising since the electron WF in a QW under magnetic field is described by the factor, which will be essential especially for excited states, as they have a maximum for large r.

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 R_max 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.