_{1}

^{*}

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.3}Ga
_{0.7}As 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 [_{0.3}Ga_{0.7}As 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 [

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_{0}: Dielectric constant, z_{0}: 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 [

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 [

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.3}Ga_{0.7}As QW of width d = 125A. Such structure has been calculated by R. Chen et al. [

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^{2}) 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 [

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