Effect of electric and magnetic fields on impurity binding energy in zinc-blend symmetric InGaN/GaN multiple quantum dots

The binding energy of ground state for hydrogenic impurity in multiple quantum dots is calculated in the framework of effective-mass approximation and using a variational method. It is shown that the binding energy is a function of the size of dots, impurity position and external fields strength. The binding energy has a maximum value when the impurity is located on the center of dots and decreases for other impurity positions. The external electric and magnetic fields change the magnitude and the position of peaks. PACS Codes 73.20.D; 71.21.La; 71.55.Eq


Background
The study of confined quantum systems has been the interesting subject of investigation since the beginning of quantum theory. The interest in the study of the physical properties of confined quantum systems such as quantum wells, wires, and dots, has increased, with the recent progress in semiconductor nanotechnology [1][2][3][4][5].
In recent years, theoretical and experimental investigations have been performed into the issue of the hydrogenic binding of an electron to a donor impurity which is confined within low-dimensional heterostructures [6,7]. The understanding of the electronic and optical properties of impurities in such systems is important because the optical and transport properties of devices made from these materials are strongly affected by the presence of shallow impurities.
The wide-band gap GaN material systems have attracted much attention for their applications in optoelectronic devices [8]. Built-in electric field is absent in zinc-blend (ZB) GaN structures because the spontaneous polarization does not exist in the ZB GaN due to the higher crystal symmetry. There has been a lot of work devoted to understanding of hydrogenic impurity states in ZB GaN quantum dots and quantum wire [9][10][11][12][13][14]. In all of those calculations, they are all base on single lowdimensional quantum structures. The impurity effects in ZB GaN-based multiple QDs have also been investigated theoretically [15,16]. In theses systems, wave function would penetrate more to the adjacent quantum dots if the barrier height or the barrier thickness is reduced.
On the other hand, the application of an external electric field can provide much valuable information about the confined impurities. Recent theoretical investigations predicated both the field induced level shifts and the field dependence of the carrier lifetime. Therefore, the impurity and the applied electric field effects on the optical properties of QDs are of great interest for fundamental physics and device applications [17,18].
To our knowledge, there have not been theoretical investigations on impurity states in ZB symmetric multiple GaN QDs under external electric and magnetic fields. In this paper the variational method is used for calculating the impurity binding energy in symmetric InGaN/GaN multiple quantum dots. In this regard, a trial wave function based on the carrier wave function in cylindrical quantum dot is introduced and the energy is calculated. In Section 2, the Hamiltonian and the calculation method http://www.nanoconvergencejournal.com/content/1/1/25 are given. The numerical calculations and discussion on typical InGaN/GaN material are presented in Section 3.

Methods
We consider a cylindrical ZB symmetric In x Ga 1−x N multiple QD of radius R and length L d surrounded by two large energy gap materials In y Ga 1−y N in the radial direction and GaN in the z-direction, Figure 1. Within the framework of effective mass approximation, the Hamiltonian description of an electron in the presence of an external field and hydrogenic impurity is given aŝ where m * , F and V (ρ, z) are the effective mass of charge carriers, electric field strength and the confinement potential. The confining potential is Using the separation of variable, the eigenfunction ofĤ 0 is written as The radial wave function of the ground state can be expressed as where barrier region (7) where Matching the wave functions and their derivatives at the boundaries, the normalization constants and energy eigenvalues, The wave function forĤ is obtained by variational method, where α and β are the variational parameters and ρ i and z i are the position of the impurity. The Hamiltonian for the system in the presence of magnetic field iŝ where A is the vector potential of magnetic field. For a homogeneous magnetic field B(0, 0, B) directed along z-axis, the vector potential is chosen as A = B × r /2. The wave function for H 0 can be written as Eq. (5), where the ground state radial wave functions f (ρ) are Whittaker functions and h(z) is the linear combination of sin(γ z) and cos(γ z) functions, where 2 2 ) and μ = eB 0 /(2m * cR * ). The impurity binding energy, E b , is obtained as

Results and discussion
In this study, the numerical calculations are carried out on a typical GaN/In  For F = 0 the binding energy has a maximum value when the impurity is located at the center of the inner dot. Also, the binding energy has two weaker symmetric local maxima for the impurity in the outer dots. The reason is that the electron wave function is mainly distributed inside the well region of the dots, and the Coulomb interaction is considerably large when the impurity is at the center of dots. The results are similar to works done by Wei, Chang and Zheng [16,19]. For F = 0, the electron is under electric force and its wave function is pushed to left side of the system. In this case the variation of binding energy is not symmetric, it has only a maximum value when the impurity is located at the center of QD1.
The variation of binding energy in terms of impurity position under magnetic field along z-direction is presented in Figure 3. As it is seen, the binding energy increases as the magnetic field is applied. Because of, the magnetic field causes the wave function of electron localizes inside the dot regions and, therefore, the interaction of electron and impurity located in the center of QDs increases. Figure 4 displays the binding energy as a function of the height of dot, L d , for different impurity positions. When the impurity is at the center of dots (z i = 0, z i = L d + L b ), the binding energy increases , it reaches a maximum for a given l 0 and then decreases as the dot height increases. It is due to the electron wave function is localized inside the dot for L d ≤ l 0 2.3 nm, and then the Coulomb interaction is increased but for L d > l 0 the distance between electron and impurity is increased and the binding energy decreases. For impurity position at the center of barrier (z i = (L d + L b )/2), the binding energy decreases as L d increases, this is because the interaction between electron and impurity decreases. The binding energy for impurity position at the end of QDs (z i = L d /2, 3L d /2 + L b ) decreases as L d increases. This is also because the wave function is distributed at the dot and the distance between electron and impurity increases as L d increases.
The effect of external electric field on the binding energy on symmetric multiple QDs is considered in Figure 5. As it is seen, the maximum value of binding energy for impurity at the center of dot decreases and is shift to left (L d < 2 nm). The presence of electric field causes the wave function is pushed to the left direction and the binding energy for impurity positions at the left side of system, such as z i = −(L d + L b ) increases, while it decreases for impurity positions at the right side.
In Figure 6, the variation of the binding energy as a function of L d and in the presence of magnetic field is shown. The binding energy increases for all impurity positions. For z i = 0 and z i = L d + L b , the binding energies have the

Conclusion
The binding energy in a multiple cylindrical quantum dots using the variational method and appropriate wave function is calculated for ZB GaN structures in the presence of electric and magnetic fields. The results clearly showed that the binding energy has three peaks, that are around the center of dots, and decreases as the dot size increases. The electric field varies the value and the position of binding energy peaks according to impurity