Effects of Applied Magnetic Field on the Optical Properties and Binding Energies Spherical GaAs Quantum Dot with Donor Impurity

Abstract

The screened modified Kratzer potential (SMKP) model is utilized to scrutinize the impacts of an applied magnetic field (MF) on the binding energies and linear and nonlinear optical properties spherical GaAs quantum dot with donor impurity (DI). To accomplish this goal, we have used the diagonalization method to numerically solve the Schrödinger equation under the effective mass approximation for obtaining the electron energy levels and related electronic wave functions. The expressions used for evaluating linear, third-order nonlinear, and total optical absorption coefficients and relative refractive index changes were previously derived within the compact density matrix method. It has been shown here that the MF and DI impacts the characteristics of the absorption coefficients and the refractive index changes. This study’s results will find application in optoelectronics and related areas.

1. Introduction

Confinement potential (CP) models are analytical mathematical expresssions that represent the physical forces that acts on a particle within a specific region in space. It has been employed to model several physical phenomena. This is because it provides us with a relatively less-expensive approach to simulate physical systems of interest compared to experimetal and advanced computational approaches [1,2,3,4]. A number of physical phenomena have been modelled and simulated by several potential models [5,6,7,8,9], such as the quarkonium interactions have been modelled by Cornell potential [10,11,12], Yukawa potential [13,14,15], extended Cornell potential [16], etc. Diatomic and Polyatomic molecules have been modelled by the Morse potential [17,18,19,20,21], improved Tietz potential [22], Deng-Fan-Eckart potential [23,24,25] and many others. Plasma interactions have been modelled by Debye Huckel potential [26,27,28], Hulthen potential [29,30,31] and a host of others [32,33]. The applications of the potentials have also been witnessed in the emerging and ever advancing field of nanotechnology, as researchers have also modelled quantum dots (QDs) using Gaussian potential [34,35], modified Gaussian potential [35,36], Rosen-Morse potential [37] and many others. It has also been showed that the application of these models also aids in gathering some useful information like the thermodyanmic, transport, magnetic and optical properties of the above mentioned systems [3,22,38,39,40,41]. However, it will be necessary to point out that the application of the models mentioned above are not limited entirely to a specific area, based on appropriate manipulations and fittings of paramters, a model used to study diatomic molecules can be applied to study QD and vice versa. These possibilities do exist but there are exceptions particularly for cases concerning the diatomic molecules because the potential studying such a system must possess explcitly the dissociation energy, equilibrrium bond length, etc [42,43].

In this article, the SMKP [2,44] is adopted to model QDs. This SMKP has been applied by several authors to study diatomic molecules and other physical systems [2,44]. In order to broaden the scope of its applicability, this model is employed to simulate QDs since it poses the potential to study such systems (see Figure 1). Moreover, our choice is inspired by the fact that potential models with many fitting parameters tends to study physical systems and/or fit with experimental data better than those with fewer parameters and SMKP falls within this class (exponential-type potentials).

The study of the optical properties of QD has received immense attention from researchers in recent times. This is due to the exigent need for advanced nano-fabricated technologies. Optical properties of zero-dimensional structures (QD) has applications in the design process of electronic devices such as laser diodes, etc. [45,46]. It is on this basis that optical properties of semiconductor heterostructures have attracted unmatched consideration [46,47,48,49]. The study of optical properties of QD supplies us with such information like the optical absorption coefficients (OACs) [47], refractive index changes (RICs) [48,50], dipole transition [46,51], oscillator strength [52,53,54] and photoionization cross-sections [52]. Furthermore, this have enticed the interest of scholars in experimental and theoretical studies in recent years [46].

External fields such the MF, Aharanov-Bohm (AB) field, etc have been shown to modify and twist the optical properties of QDs in recent times. Notable amongst them is the removal degeneracy from the energy spectra by the magnetic and AB fields [2,4]. By modelling the QD with finite CP, Al et al [55] pointed out in their study that MF affects the impurity binding energies greatly when the dot radii is large also and the OACs. Again, Al et al. [46] have noted that the presence of MF removes degeneracy from the energies and shift the energy higher (blue shift). In addition, it has been showed that the cumulative effect of donor position, MF and dot size have been observed to make significant changes to the binding energies and OACs.

Having carefully highlighted the above, it is certain that the relevant motive for the present study is established. Being relevant to the field of device applications, the QDs are systems that exhibit remarkable optical and electronic properties. These can be properly tuned by suitably modifying the composition, structural parameters, and the application of external probes. Here, the optical properties and binding energies of donor impurity of a spherical GaAs QD structure with SMKP confinement potential have been investigated under the influence of externally applied MF, as well as the change of structure parameters. This is done with the purpose of analyzing the effects on the energy spectrum of externally applied fields as well as structural modifications. From such an outcome, the coefficients of linear and nonlinear optical absorption and refractive index change will be separately analyzed. The organization of the article is as follows: in Section 2, the formulation and theoretical solution of the problem are given, numerical results and their discussions are presented in Section 3 and our conclusions are given in Section 4.

2. Theoretical Framework

For an electron in motion in a spherical QD with the nonexistence of DI in the effective mass (EM) approximation, the Hamiltonian is given as [55];

(1)$H^{\left(0\right)}={\displaystyle \frac{{\overrightarrow{p}}^2}{2m^{*}}}+V\left(r\right)$

where $\overrightarrow{p}$ is the momentum operator $m^{*}$ is the electron EM, and $V\left(r\right)=D_e{(q-r_e{e}^{-\alpha r}/r)}^2$ is the spherically symmetric SMKP confinement potential with $D_e$ is the dissociation energy, $r_e$ is the internuclear distance, $\alpha$ is the screening parameter, q is the control parameter and r is the internuclear distance [2,44].

The Hamiltonian of the system comprising of the MF and donor impurity (DI) atom located at the center of the finite spherical QD is given by [55]:

(2)$H={\displaystyle \frac{1}{2m^{*}}}{\left(\overrightarrow{p}+{\displaystyle \frac{e}{c}}\overrightarrow{A}\right)}^2+V\left(r\right)-{\displaystyle \frac{e^{2}Z}{\epsilon |\overrightarrow{r}|}},$

where e, $\overrightarrow{A}$, $\epsilon$, $\overrightarrow{r}$ and c are charge of the electron, vector potential (VP) of the MF, the stationary dielectric constant of the QD material, electron position vector and the speed of light in vacuum, respectively with $Z=0(Z=1)$ accounts for absence (presence) of impurity. This VP $\overrightarrow{A}$ is selected such that the MF is $\overrightarrow{B}=\overrightarrow{▿}\times \overrightarrow{A}$. Invoking the Coulomb gauge ($\overrightarrow{▿}.\overrightarrow{A}=0$), the VP satisfies the relationship $\overrightarrow{A}=\overrightarrow{B}\times \overrightarrow{r}/2$. Considering that the MF is applied along the z-axis, the Hamiltonian (Equation (2)) in its dimensionless expression is given as [46]

(3)$H={\displaystyle \frac{1}{r^2}}{\displaystyle \frac{\partial }{\partial r}}\left(r^2{\displaystyle \frac{\partial }{\partial r}}\right)+{\displaystyle \frac{l(l+1)}{r^2}}+V\left(r\right)+m\gamma +{\displaystyle \frac{{\gamma }^2{r}^2{sin}^2\left(\theta \right)}{4}}-{\displaystyle \frac{2Z}{r}},$

where $\gamma =e{a}_B^2\overrightarrow{B}/\hslash c$ with $a_B={\hslash }^2\epsilon /m^{*}{e}^2$ is the effective Bohr radius. Furthermore, l and m are the angular and magnetic momentum quantum numbers, correspondingly. The energy levels in the presence ($Z=1$) and absence $(Z=0)$ of DI and their associated wave functions are derived by solving the eigenvalue corresponding to the Hamiltonian in Equation (3) is given as [53];

(4)$H{\psi }_{nlm}(r,\theta ,\varphi )=E_{nlm}{\psi }_{nlm}(r,\theta ,\varphi ),$

where n is the principal quantum number and $\varphi$ is the angle between the $\overrightarrow{r}$-projection on the xy-plane and the x-axis. Under the influence of MF and in the presence of DI atom in the QD system, Equation (4) is not solveable analytically. The solutions of Equation (4) can be obtained by employing the approach that is founded on the expansion of the electronic states over a complete orthogonal basis as presented as follows [55,56]

(5)${\psi }_{nlm}(r,\theta ,\varphi )=\sum _{j=1}^{\sigma }{\mathcal{C}}_{nlm}{\psi }_{n_j{l}_{j}m}^{\left(0\right)}(r,\theta ,\varphi ),$

where $\sigma$ is the size of the basis set with an accurateness of 0.001 meV of the energy value. ${\mathcal{C}}_{n_j{l}_j}$ is the j-th coefficient of the expansion (where the electronic states are categorized matching to a particular m value) and ${\psi }_{n_j{l}_{j}m}^{\left(0\right)}(r,\theta ,\varphi )$ are the exact solutions of the infinite QD (where the radius is three times that of the QD with finite confinement potential $V\left(r\right)$). We point out here that when the total Hamiltonian matrix elements are computed, The wave function expansion will incorporate the terms of the finite confinement potential $V_0\equiv D_e$ that exists in the interval $a<r<3a$. The MF is known to eliminate the spherical symmetry of the problem, the orbital quantum number is an invalid quantum number, nonetheless three quantum numbers $(n,l,m)$ are utilized to illustrate the new states suitably. In the absence of MF and DI [55]

(6)${\psi }_{n_j{l}_{j}m}^{\left(0\right)}(r,\theta ,\varphi )=R_{nl}^{\left(0\right)}\left(r\right)Y_{lm}(\theta ,\varphi ),$

where $R_{nl}^{\left(0\right)}\left(r\right)$ is the radial part of the non-correlated electron eigenfunction. $Y_{lm}(\theta ,\varphi )=e^{im\varphi }{P}_{lm}(cos\left(\theta \right))$ with $P_{lm}(cos\left(\theta \right))$ being the associated Legendre polynomial. Notice that m is an integer number which guarantee the azimuthal symmetry with respect to the axial z-axis. The impurity binding energy is given by $E_{nlm}^b=E_{nlm}(Z=0)-E_{nlm}(Z=1)$, where $E_{nlm}(Z=0)$ and $E_{nlm}(Z=1)$ are the energies of the electron in the absence and presence of DI, respectively. Intersubband transitions involves the ground and the first excited subband, once the incident photon energy $\hslash \omega$ is equal to the difference in energy between subbands, this is mathematically given as; $\hslash \omega =E_{21}=E_2-E_1$, where $E_1$ and $E_2$ are the ground and first excited states energy levels. In this study, we used z direction and x-y plane polarized EM radiation. A medium polarization is stimulated by the subband transition as a result of the ordering of the electric dipole moments in the direction of EM field [35,37,57]. The mathematical expressions of the linear ${\chi }^{\left(1\right)}$ and the third-order nonlinear ${\chi }^{\left(3\right)}$, using the density matrix approximation, susceptibility coefficients are presented as follows

(7)$\begin{array}{ccc}\hfill {\epsilon }_0{\chi }^{\left(1\right)}& =& {\displaystyle \frac{{\sigma }_{\nu }{\left|M_{21}\right|}^2}{E_{21}-\hslash \omega -i\hslash {\Gamma }_{12}}},\hfill \end{array}$

(8)$\begin{array}{ccc}\hfill {\epsilon }_0{\chi }^{\left(3\right)}& =& {\displaystyle \frac{{\sigma }_{\nu }{\left|M_{21}\overline{E}\right|}^2}{E_{21}-\hslash \omega -i\hslash {\Gamma }_{12}}}\left. [{\displaystyle \frac{4|M_{21}{|}^2}{{(E_{21}-\hslash \omega )}^2+{(\hslash {\Gamma }_{12})}^2}}-{\displaystyle \frac{4|M_{22}-M_{21}{|}^2}{{(E_{21}-i\hslash {\Gamma }_{12})}^2+{(E_{21}-\hslash \omega -i\hslash {\Gamma }_{12})}^2}}\right],\hfill \end{array}$

where $M_{ij}$ is the matrix elements of the electric dipole moment, $E_{ij}=E_i-E_j$ is the energy difference between the two electronic states, $\hslash \omega$ is the incident photon energy, ${\sigma }_v$ is the carrier density and ${\Gamma }_{12}$ is the relaxation rate. The matrix elements of electric dipole moment can be evaluated using $M_{ij}=⟨{\psi }_j|ersin\left(\theta e^{\pm i\varphi }\right)|{\psi }_j⟩$ for $i,j=1,2$. Therefore, the dipole transition matrix element in general is written as

(9)$M_{ij}^{\left(xyz\right)}=\left\{\begin{array}{cc}{M}_{ij}^{\left(z\right)}\hfill & if\mathsf{\Delta }m=0\hfill \\ M_{ij}^{\left(xy\right)}\hfill & if\mathsf{\Delta }m=\pm 1.\hfill \end{array}\right.$

The intersubband optical transition probability amplitude between the two electronic states, such as the lowest (initial) and the first excited (final), is proportional to the matrix element of electric dipole moment. Because the electric dipole moment is unaffected by spin, optical transitions can only occur between spin states with the same spin. When the selection rules $\mathsf{\Delta }l=\pm 1$ and $\mathsf{\Delta }m=0,\pm 1$ satisfied, the matrix elements of electric dipole moment are non-vanishing [46,48]. The OAC-$\alpha \left(\omega \right)$ is obtained from the susceptibility $\chi \left(\omega \right)$ using the relation given as $\alpha \left(\omega \right)=\omega \sqrt{{\mu }_0/\epsilon }\Im \left[{\epsilon }_0\chi \left(\omega \right)\right]$, where ${\mu }_0$ vacuum permeability and $\epsilon =n_r^2{\epsilon }_0$ permittivity of the material with $n_r$ is the refractive index. The linear and third-order nonlinear OACs can be expressed as

(10)$\begin{array}{ccc}\hfill {\alpha }^{\left(1\right)}\left(\omega \right)& =& \omega \sqrt{{\mu }_0/ϵ}{\displaystyle \frac{\hslash {\sigma }_v{\Gamma }_{12}{\left|M_{21}\right|}^2}{{(E_{21}-\hslash \omega )}^2+{\hslash }^2{\Gamma }_{12}^2}},\hfill \\ \hfill {\alpha }^{\left(3\right)}(\omega ,I)& =& {\displaystyle \frac{-\hslash I\omega \sqrt{{\mu }_0/ϵ}{\sigma }_v{\Gamma }_{12}{\left|M_{21}\right|}^2}{2{\epsilon }_0{n}_{r}c{[{(E_{21}-\hslash \omega )}^2-{\hslash }^2{\Gamma }_{12}^2]}^2}}\left\{4|M_{21}{|}^2\right.\hfill \end{array}$

(11)$\begin{array}{ccc}& -& \left.{\displaystyle \frac{|M_{22}-M_{11}{|}^2[3E_{21}^2-4\hslash \omega E_{21}+{\hslash }^2({\omega }^2-{\Gamma }_{12})]}{E_{21}^2+{\hslash }^2{\Gamma }_{12}^2}}\right\},\hfill \end{array}$

where $I=2{\epsilon }_0{n}_{r}c{\left|\overline{E}\right|}^2$ represents the incident optical intensity. Because of the spherical symmetry of the system ($M_{ii}=0$), the third-order nonlinear OAC-${\alpha }^{\left(3\right)}(\omega ,I)$ is negative and also it is proportional to the incident optical intensity-I. Then, the total OAC $\alpha (\omega ,I)$ is given by $\alpha (\omega ,I)={\alpha }^{\left(1\right)}(\omega ,I)+{\alpha }^{\left(3\right)}(\omega ,I)$. Since $|M_{jj}-M_{ii}|$ equals to zero for the on-center donor impurity, the third order nonlinear OAC becomes ${\alpha }^{\left(3\right)}(\omega ,I)=-\sqrt{{\mu }_0/{\epsilon }_r}2\hslash \omega I{\sigma }_v{\Gamma }_{12}{\left|M_{21}\right|}^4/{\epsilon }_0{n}_{r}c[{(E_{21}-\hslash \omega )}^2+{\hslash }^2{\Gamma }_{12}^2]$. Since the value of the $\hslash {\Gamma }_{ij}$ can not be ignored and the corresponding resonance condition is given by $\hslash \omega \sqrt{E_{ij}^2+{\hslash }^2{\Gamma }_{ij}^2}$. The RIC $\mathsf{\Delta }n\left(\omega \right)$ can be obtained from susceptibility $\chi \left(\omega \right)$ by $\mathsf{\Delta }n\left(\omega \right)/n_r=\Re [\chi \left(\omega \right)/2n_r^2]$ [38,46]. Then the linear $\mathsf{\Delta }{n}^{\left(1\right)}\left(\omega \right)$ and third-order nonlinear $\mathsf{\Delta }{n}^{\left(3\right)}\left(\omega \right)$ RICs are given by

(12)$\begin{array}{ccc}\hfill \frac{\mathsf{\Delta }{n}^{\left(1\right)}\left(\omega \right)}{n_r}& =& {\displaystyle \frac{{\sigma }_v{\left|M_{21}\right|}^2}{2n_r^2{\epsilon }_0}}{\displaystyle \frac{E_{21}-\hslash \omega }{{(E_{21}-\hslash \omega )}^2+{\hslash }^2{\Gamma }_{12}}},\hfill \\ \hfill \frac{\mathsf{\Delta }{n}^{\left(3\right)}\left(\omega \right)}{n_r}& =& {\displaystyle \frac{{\sigma }_v{\left|M_{21}\right|}^2{\mu }_{0}cI}{4n_r^2{\epsilon }_0}}\left\{{\displaystyle \frac{4\left(E_{21-\hslash \omega }\right){\left|M_{21}\right|}^2}{{(E_{21}-\hslash \omega )}^2-{\hslash }^2{\Gamma }_{12}^2}}\right.\hfill \end{array}$

(13)$\begin{array}{ccc}\hfill & -& \left.{\displaystyle \frac{|M_{22}-M_{11}{|}^2{E}_{21}-\hslash \omega \left(E_{21}{(E_{21}-\hslash \omega )}^2-{\hslash }^2{\Gamma }_{12}^2\right)-{\hslash }^2{\Gamma }_{12}^2\left(2E_{21-\hslash \omega }\right)}{[{(E_{21}-\hslash \omega )}^2-{\hslash }^2{\Gamma }_{12}^2][{(E_{21}-\hslash \omega )}^2+{\hslash }^2{\Gamma }_{12}^2]}}\right\}.\hfill \end{array}$

It is clear from Equations (12) and (13), the third-order nonlinear RIC-$\mathsf{\Delta }{n}^{\left(3\right)}/n_r$ has contrasting sign with the sign of $\mathsf{\Delta }{n}^{\left(1\right)}/n_r$ and it is proportionate to the incident optical intensity-I. Therefore, the OACs (Eqautions (10) and (11)) and RICs (Eqautions (12) and (13)), depends only on the matrix element $M_{21}$ of the dipole moment. Once a photon with energy-$\hslash \omega$ is released into the system, the matrix element $M_{21}$ provides the probability amplitude of one electron optical transition between the state $|1⟩$ and $|2⟩$ described by the wave functions.

3. Results and Discussion

This study presents the optical properties of the GaAs QD structure using SMKP. Numerical computations are carried out for this QD structure. Fitting parameters used to achieve our set objectives are presented as follows: $\mu =0.067$ me, ${\sigma }_v=5\times {10}^{22}$ m${}^{-3}$ and $n_r=3.2$ [35]. Consequently, ${\epsilon }_0=8.854\times {10}^{-12}$ F/m is the vacuum permittivity, ${\mu }_0=4{\pi }^{-7}$ H/m is the vacuum permeability, and ${\Gamma }_{12}=1/T_{12}$ where $T_{12}=0.14$ ps, $D_e=224.46$ meV, $r_e=1.5$, $q=1$, $\alpha =0.05$, $m_e=0.067m_0$, $\epsilon =13.18$, $\sigma =5\times {10}^{22}$ m${}^{-3}$, $I=200$ MW/m${}^2$ and $a_B=104.24$ A.

The experimental value for the depth of the potential well-$D_e$ and radius of the dot- $R_0$ are selected in line with what obtains in literature [57]. It is important to point out here that ${\sigma }^{+}$ corresponds to transition between $1p_{+1}$ and $1s$ state ($1p_{+1}\to 1s$), ${\sigma }^{-}$ corresponds to transition between $1p_{-1}$ and $1s$ ($1p_{-1}\to 1s$) and $\pi$ transition between $1p_0$ and $1s$ ($1p_0\to 1s$). In addition, p corresponds to $l=1$, while s corresponds to $l=0$. l is the angular momentum quantum number, $-1$, 0 and $+1$ are m values and m is the magnetic momentum quantum number.

The SMKP is the adopted confining potential for this study, it has finite depth and range. This potential is a continuous function that has four parameters that can be adjustable at will to fit a physical system of interest. In this study, our goal is to study the optical properties of a QD with the SMKP with DI in the presence of an applied MF. Historically, it has been established that for a potential model to be used to study intersubband transition occurrences and related properties of QDs, the potential must be of a quadratic-like form just like the SMKP (see Figure 1). Figure 1 shows the confining SMKP as a function of r. The shape of this model is in agreement with the shape of other confinement potential adopted previously for such considerations [36,37,58]. However, our model is an improvement to the aforementioned in the sense that our potential has more fitting parameters and it is a general case of the Kratzer potential [59]. Figure 1 also shows the probability density, ${\left|\varphi \left(r\right)\right|}^2$ for $1s$ and $1p$ states and in the presence and absence of DI. It is seen that probability density increases for higher states but there is no discernable difference (overlap) in the peak positions and this implies that the particle has an intense localization. Probability density offers information on the concentration of the density distribution.

Figure 2a,b show the electron energies versus the MF in the absence and presence of DI for $1s$, $1p_{-1}$, $1p_0$ and $1p_{+1}$. However, for $1p_{-1}$, $1p_0$ and $1p_{+1}$ states, there is a convergence at $B=0T$ and with MF there is a splitting which definitely implies elimination of degeneracy due to Zeeman effect. Figure 2c shows the binding energies versus the MF for $1s$, $1p_{\pm 1}$ and $1p_0$. The energy increases in the order presented as follows; $1p_0$, $1p_{\pm 1}$ and $1s$. The binding energy is also seen to rise with rising MF.

Figure 3a shows the OACs as a function of the photon energy for ${\sigma }^{\pm }$$(Z=0)$, ${\sigma }^{\pm }$$(Z=1)$, $\pi$$(Z=0)$ and $\pi$$(Z=1)$ transitions. To study the OACs of this system, the DI position is fixed at $r_d=14$ nm. It is observed that the linear and total OACs first increase and then decrease as the photon energy increases while negative nonlinear OAC shows the contrary variation pattern, this observed in both transitions. A quasi-resonance peak is observed about photon energy of $\approx 8$ meV for total OAC. The OAC is higher in the presence of DI than in the absence for ${\sigma }^{\pm }$-transition, while it is lower for $\pi$-transition. The DI shifts the peak position of linear and nonlinear OACs slightly towards blue (high energy) due to the increase in transition energy, while the absence of the DI causes these peaks to shift to red (low energy) in the ${\sigma }^{\pm }$-transition and the converse is observed in $\pi$-transition. We may conclude here that the influence of these DI in the case of the ${\sigma }^{\pm }$ and $\pi$-transitions on the OACs is great and can be instrumental in the manipulation of the OACs. In Figure 3b we plot RICs as a function of the photon energy for ${\sigma }^{\pm }$$(Z=0)$, ${\sigma }^{\pm }$$(Z=1)$, $\pi$$(Z=0)$ and $\pi$$(Z=1)$ transitions. The RICs firstly decrease and then slightly increases with rising of photon energy as shown in Figure 3b. When the DI is present, the RICs increase (decrease) for ${\sigma }^{\pm }$-$(\pi$-)-transition and their positions move to the greater (smaller) energies. This is attributed to the fact that the energy difference between the two electronic states increases for ${\sigma }^{\pm }$-transition, while it decreases for $\pi$-transition.

Figure 4a,d show the OACs as a function of the photon energy for ${\sigma }^{-}$-transition in the present presence and absence of DI respectively with varying MF. The linear and total OACs rise first and decrease with rising photon energy, but their resonant peaks decrease with rising MF. However, there is no discernable difference in the OACs for all values of B-field between the presence and absence of DI for the considered transition. In Figure 4b,e, the linear and total OACs again rises first and decrease with rising photon energy, but they do not show a discernable variation pattern as the MF rises for the transition in the presence and absence of DI In addition, the presence of DI reduces the peak amplitude of the $\pi$-transition while red-shifting the peak position. Figure 4c,f shows the variation of the OACs as a function of the photon energy with varying MF in the absence and presence of DI for ${\sigma }^{+}$-transition. The linear and total OACs first increase and then decrease as the photon energy increases while negative nonlinear OAC shows the opposite change. The increase in B-field parameters shifts the peak position of OACs slightly towards blue (high energy) due to the increase in transition energy. In addition, the peak amplitudes of linear and nonlinear OACs increase with the increase of B-field, due to more overlapping of the wavefunctions of the ground and first excited states and hence the increase in the dipole matrix element. Moreover, the peak amplitude of the total OAC also increase with the increase of the MF, since the positive increase in the peak amplitude of the linear OAC is dominant on the negative increase in the peak amplitude of the nonlinear OAC. However, there is no discernable difference in the OAC for the absence and presence of DI. We may conclude here that the influence of the effects of MF in ${\sigma }^{\pm }$-transitions and the effects of DI in $\pi$-transitions are abundant and these effects play an influential role in the control of the behavior of OACs.

Figure 5 shows the RICs as a function of the photon energy for ${\sigma }^{\pm }$$(Z=0)$, ${\sigma }^{\pm }$$(Z=1)$, $\pi$$(Z=0)$ and $\pi$$(Z=1)$ transitions. In Figure 5a,c,d,f, the RICs firstly decrease and then increase with increasing photon energy but slightly shifts upward as a function of the photon energy for ${\sigma }^{+}$-transitions when the MF is rising and but the contrary is observed in the ${\sigma }^{-}$-transitions. This is attributed to the fact that the energy difference between the two electronic states changes as shown in Table 1. In Figure 3b,e shows the same behavior with rising of photon energy for total RIC in the $\pi$-transition in the presence and absence of DI. There is again no discernable difference in the trend for $\pi$-transition as the MF rises. Again, we note here that the impact of the effects of MF in ${\sigma }^{\pm }$-transitions and the effects of DI in $\pi$-transitions are ample and these effects play a significant function in the regulation of the behavior of OACs.

4. Conclusions

This study presents the effects of applied MF on the optical properties and binding energies (for ${\sigma }^{\pm }$ and $\pi$) of donor DI in spherical GaAs quantum dot. In view of this, the SMKP potential has been utilized. To derive the eigenvalues equations and wave functions, the Schrodinger equation in spherical coordinate was solved by diagonalization approach. On variation of the B-field in the presence and absence of DI, the absorption coefficients and refractive index changes switch from lower to higher energies (blue shift) and also from higher to lower energies (red shift) for several transitions considered. It can be inferred that the MF and DI are greatly instrumental for manipulation of the optical properties of the GaAs quantum dot. Our results are vital in the design of nano-scale electronic devices such as quantum dot solar cells, quantum dot lasers, quantum dot photodetectors, quantum dot infrared photodetectors, phototransistors, and light-emitting diodes.

Author Contributions

C.O.E., E.B.A. and F.U.: Conceived and designed the study, acquired, analyzed and interpreted the data, and handled the review, R.E., N.R. and N.A. handled the computational analysis. N.R., M.A. and S.A.A.: revised the manuscript for important intellectual content, and contributed in general recommendations and data interpretation. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

C.O.E. acknowledges eJDS (ICTP) and M. Asjad acknowledges Khalifa University for financial support in the project No. 8474000358 (FSU-2021-018).

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. The confining SMKP (V(r)) and probability density [Forumla omitted. See PDF.] as a function of r.

Figure 2. (a) Electron energies versus the MF in the absence of impurity; (b) Electron energies versus the magnetic field in the presence of DI. (c) Binding energies versus the MF.

Figure 3. (a) The plot of optical absorption coefficient as a function of the photon energy. (b) The Refractive index changes as a function of the photon energy. Both plots are for [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.][Forumla omitted. See PDF.] transitions. In (a,b), the dashed, doted and solid curves are for linear, third-order nonlinear and total optical absorption coefficients (Refractive Index Changes) respectively.

Figure 4. Plot of optical absorption coefficients as a function of the photon energy for [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.][Forumla omitted. See PDF.] transitions and varying MF (a–f). The dashed, doted and solid curves represent linear, third-order nonlinear and total optical absorption coefficients respectively.

Figure 5. Plot of refractive index changes as a function of the photon energy for fixed values of [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.], [Forumla omitted. See PDF.][Forumla omitted. See PDF.] and [Forumla omitted. See PDF.][Forumla omitted. See PDF.] transitions and varying MF (a–f). The dashed, doted and solid curves account for linear, third-order nonlinear and total refractive index changes respectively.

Table 1

The subband energy differences and dipole moment matrix elements for (a), (b), (c), (d), (e) and (f) respectively, in Figure 4 and Figure 5.

B	0	2T	4T
$E_{ij}$	4.49174	2.83055	1.29578
$M_{ij}$	13.9468	13.9282	13.8737
$E_{ij}$	4.49174	4.39326	4.10013
$M_{ij}$	9.86191	9.85168	9.82146
$E_{ij}$	4.49174	6.28157	8.19782
$M_{ij}$	13.9468	13.9282	13.8737
$E_{ij}$	5.4793	3.81847	2.28483
$M_{ij}$	13.824	13.8066	13.7554
$E_{ij}$	2.75256	2.65404	2.36068
$M_{ij}$	9.74208	9.73287	9.70561
$E_{ij}$	5.4793	7.26949	9.18687
$M_{ij}$	13.824	13.8066	13.7554

References

1. Jia, C.-S.; Jia, Y. Relativistic rotation-vibrational energies for the Cs2 molecule. Eur. Phys. J. D; 2017; 71, 3. [DOI: https://dx.doi.org/10.1140/epjd/e2016-70415-y]

2. Edet, C.O.; Ikot, A.N. Effects of Topological Defect on the Energy Spectra and Thermo-magnetic Properties of CO Diatomic Molecule. J. Low Temp. Phys.; 2021; 203, 84. [DOI: https://dx.doi.org/10.1007/s10909-021-02577-9]

3. Jia, C.-S.; Zhang, L.-H.; Wang, C.-W. Thermodynamic properties for the lithium dimer. Chem. Phys. Lett.; 2017; 667, 211. [DOI: https://dx.doi.org/10.1016/j.cplett.2016.11.059]

4. Edet, C.O.; Ikot, A.N.; Onyeaju, M.C.; Okorie, U.S.; Rampho, G.J.; Lekala, M.L.; Kaya, S. Thermo-magnetic properties of the screened Kratzer potential with spatially varying mass under the influence of Aharanov-Bohm (AB) and position-dependent magnetic fields. Phys. E Low-Dimens. Syst. Nanostruct.; 2021; 131, 114710. [DOI: https://dx.doi.org/10.1016/j.physe.2021.114710]

5. Edet, C.O.; Mahmoud, S.; Inyang, E.P.; Ali, N.; Aljunid, S.A.; Endut, R.; Ikot, A.N.; Asjad, M. Non-Relativistic Treatment of the 2D Electron System Interacting via Varshni–Shukla Potential Using the Asymptotic Iteration Method. Mathematics; 2022; 10, 2824. [DOI: https://dx.doi.org/10.3390/math10152824]

6. Edet, C.O.; Osang, J.E.; Ali, N.; Agbo, E.P.; Aljunid, S.A.; Endut, R.; Ettah, E.B.; Khordad, R.; Ikot, A.N.; Asjad, M. Non-Relativistic Energy Spectra of the Modified Hylleraas Potential and Its Thermodynamic Properties in Arbitrary Dimensions. Quantum Rep.; 2022; 4, pp. 238-250. [DOI: https://dx.doi.org/10.3390/quantum4030016]

7. Edet, C.O.; Ettah, E.B.; Aljunid, S.A.; Endut, R.; Ali, N.; Ikot, A.N.; Asjad, M. Global Quantum Information-Theoretic Measures in the Presence of Magnetic and Aharanov-Bohm (AB) Fields. Symmetry; 2022; 14, 976. [DOI: https://dx.doi.org/10.3390/sym14050976]

8. Edet, C.O.; Nwbabuzor, P.O.; Ettah, E.B.; Duque, C.A.; Ali, N.; Ikot, A.N.; Mahmoud, S.; Asjad, M. Magneto-transport and thermal properties of the Yukawa potential in cosmic string space-time. Results Phys.; 2022; 39, 105749. [DOI: https://dx.doi.org/10.1016/j.rinp.2022.105749]

9. Edet, C.O.; Amadi, P.O.; Ettah, E.B.; Ali, N.; Asjad, M.; Ikot, A.N. The magnetocaloric effect, thermo-magnetic and transport properties of LiH diatomic molecule. Mol. Phys.; 2022; 120, e2059025. [DOI: https://dx.doi.org/10.1080/00268976.2022.2059025]

10. Abu-Shady, M.; Edet, C.O.; Ikot, A.N. Non-relativistic quark model under external magnetic and Aharanov–Bohm (AB) fields in the presence of temperature-dependent confined Cornell potential. Can. J. Phys.; 2021; 99, 1024. [DOI: https://dx.doi.org/10.1139/cjp-2020-0101]

11. Hall, R.L.; Saad, N. Schrödinger spectrum generated by the Cornell potential. Open Phys.; 2015; 13, pp. 83-89. [DOI: https://dx.doi.org/10.1515/phys-2015-0012]

12. Fulcher, L.P.; Chen, Z.; Yeong, K.C. Energies of quark-antiquark systems, the Cornell potential, and the spinless Salpeter equation. Phys. Rev. D; 1993; 47, 4122. [DOI: https://dx.doi.org/10.1103/PhysRevD.47.4122]

13. Yukawa, H. On the interaction of elementary particles. I. Proc. Phys.-Math. Soc. Jpn.; 1935; 17, 48.

14. Yukawa, T. Exchange potential between “nearly” identical nuclei. Phys. Rev. C; 1974; 9, 864. [DOI: https://dx.doi.org/10.1103/PhysRevC.9.864]

15. Edet, C.O.; Lima, F.C.E.; Almeida, C.A.S.; Ali, N.; Asjad, M. Quantum Information of the Aharanov–Bohm Ring with Yukawa Interaction in the Presence of Disclination. Entropy; 2022; 24, 1059. [DOI: https://dx.doi.org/10.3390/e24081059]

16. Abu-Shady, M.; Khokha, E.M. Bound state solutions of the Dirac equation for the generalized Cornell potential model. Int. J. Mod. Phys. A; 2021; 36, 2150195. [DOI: https://dx.doi.org/10.1142/S0217751X21501955]

17. Morse, P.M. Diatomic molecules according to the wave mechanics. II. Vibrational levels. Phys. Rev.; 1929; 34, 57. [DOI: https://dx.doi.org/10.1103/PhysRev.34.57]

18. Morse, P.M.; Stueckelberg, E.C.G. Diatomic molecules according to the wave mechanics I: Electronic levels of the hydrogen molecular ion. Phys. Rev.; 1929; 33, 932. [DOI: https://dx.doi.org/10.1103/PhysRev.33.932]

19. Girifalco, L.A.; Weizer, V.G. Application of the Morse potential function to cubic metals. Phys. Rev.; 1959; 114, 687. [DOI: https://dx.doi.org/10.1103/PhysRev.114.687]

20. Gomez, I.S.; Santos, E.S.; Abla, O. Morse potential in relativistic contexts from generalized momentum operator: Schottky anomalies, Pekeris approximation and mapping. Mod. Phys. Lett. A; 2021; 36, 2150140. [DOI: https://dx.doi.org/10.1142/S0217732321501406]

21. Okoi, P.O.; Edet, C.O.; Magu, T.O. Relativistic treatment of the Hellmann-generalized Morse potential. Rev. Mex. Fis.; 2019; 66, pp. 1-13. [DOI: https://dx.doi.org/10.31349/RevMexFis.66.1]

22. Jia, C.-S.; Wang, C.-W.; Zhang, L.-H.; Peng, X.-L.; Zeng, R.; You, X.-T. Partition function of improved Tietz oscillators. Chem. Phys. Lett.; 2017; 676, 150. [DOI: https://dx.doi.org/10.1016/j.cplett.2017.03.068]

23. Deng, Z.H.; Fan, Y.P. A Potential Function of Diatomic Molecules. J. Shandong Univ. (Nat. Sci.); 1957; 1, 11.

24. Dong, S.-H. Relativistic treatment of spinless particles subject to a rotating Deng—Fan oscillator. Commun. Theor. Phys.; 2011; 55, 969. [DOI: https://dx.doi.org/10.1088/0253-6102/55/6/05]

25. Edet, C.O.; Okorie, U.S.; Osobonye, G.; Ikot, A.N.; Rampho, G.J.; Sever, R. Thermal properties of Deng–Fan–Eckart potential model using Poisson summation approach. J. Math. Chem.; 2020; 58, 989. [DOI: https://dx.doi.org/10.1007/s10910-020-01107-4]

26. McChesney, M. Debye-Huckel Plasma Corrections. AIAA J.; 1963; 1, 1666. [DOI: https://dx.doi.org/10.2514/3.1881]

27. Kar, S. Dynamic dipole polarizability of the helium atom with Debye-Hückel potentials. Phys. Rev. A; 2012; 86, 062516. [DOI: https://dx.doi.org/10.1103/PhysRevA.86.062516]

28. Zammit, M.C.; Fursa, D.V.; Bray, I. Electron scattering in a helium Debye plasma. Chem. Phys.; 2012; 398, 214. [DOI: https://dx.doi.org/10.1016/j.chemphys.2011.04.034]

29. Mahdavi, M.; Akbarian, H. Calculation of effective potential of Hulthen interaction for a nonideal dense plasma using dielectric response method. Chin. J. Phys.; 2021; 77, pp. 1381-1386. [DOI: https://dx.doi.org/10.1016/j.cjph.2021.12.007]

30. Yoon, J.H.; Yun, Y.S. Hulthen potential for the plasma-screening effect on the continuum electron capture rate in the solar core. J. Korean Phys. Soc.; 2000; 37, 73.

31. Bahar, M.K.; Soylu, A.; Poszwa, A. The hulthén potential model for hydrogen atoms in debye plasma. IEEE Trans. Plasma Sci.; 2016; 44, 2297. [DOI: https://dx.doi.org/10.1109/TPS.2016.2604421]

32. Sedehi, H.R.; Khordad, R. Magnetocaloric effect, magnetic susceptibility and specific heat of tuned quantum dot/ring systems. Phys. E Low-Dimens. Syst. Nanostruct.; 2021; 134, 114886. [DOI: https://dx.doi.org/10.1016/j.physe.2021.114886]

33. Ikot, A.N.; Okorie, U.; Ngiangia, A.T.; Onate, C.A.; Edet, C.O.; Akpan, I.O.; Amadi, P.O. Bound state solutions of the Schrödinger equation with energy-dependent molecular Kratzer potential via asymptotic iteration method. Eclet. Quim. J.; 2020; 45, 65. [DOI: https://dx.doi.org/10.26850/1678-4618eqj.v45.1.2020.p65-77]

34. Boyacioglu, B.; Chatterjee, A. Dia-and paramagnetism and total susceptibility of GaAs quantum dots with Gaussian confinement. Phys. E Low-Dimens. Syst. Nanostruct.; 2012; 44, 1826. [DOI: https://dx.doi.org/10.1016/j.physe.2012.05.001]

35. Khordad, R. Pressure effect on optical properties of modified Gaussian quantum dots. Phys. B Condens. Matter; 2012; 407, 1128. [DOI: https://dx.doi.org/10.1016/j.physb.2012.01.095]

36. Gharaati, A.; Khordad, R. A new confinement potential in spherical quantum dots: Modified Gaussian potential. Superlattices Microstruct.; 2010; 48, 276. [DOI: https://dx.doi.org/10.1016/j.spmi.2010.06.014]

37. Khordad, R.; Mirhosseini, B. Linear and nonlinear optical properties in spherical quantum dots: Rosen-Morse potential. Opt. Spectrosc.; 2014; 117, 434. [DOI: https://dx.doi.org/10.1134/S0030400X14090100]

38. Ebeling, W.; Sokolov, I. Statistical Thermodynamics and Stochastic Theory of Nonequilibrium Systems; World Scientific Publishing Company: Singapore, 2005; Volume 8.

39. Oyewumi, K.J.; Falaye, B.J.; Onate, C.A.; Oluwadare, O.J.; Yahya, W.A. Thermodynamic properties and the approximate solutions of the Schrödinger equation with the shifted Deng–Fan potential model. Mol. Phys.; 2014; 112, 127. [DOI: https://dx.doi.org/10.1080/00268976.2013.804960]

40. Jia, C.-S.; Zhang, L.-H.; Peng, X.-L.; Luo, J.-X.; Zhao, Y.-L.; Liu, J.-Y.; Guo, J.-J.; Tang, L.-D. Prediction of entropy and Gibbs free energy for nitrogen. Chem. Eng. Sci.; 2019; 202, 70. [DOI: https://dx.doi.org/10.1016/j.ces.2019.03.033]

41. Jia, C.-S.; Wang, C.-W.; Zhang, L.-H.; Peng, X.-L.; Tang, H.-M.; Zeng, R. Enthalpy of gaseous phosphorus dimer. Chem. Eng. Sci.; 2018; 183, 26. [DOI: https://dx.doi.org/10.1016/j.ces.2018.03.009]

42. Jia, C.-S.; Diao, Y.-F.; Liu, X.-J.; Wang, P.-Q.; Liu, J.-Y.; Zhang, G.-D. Equivalence of the Wei potential model and Tietz potential model for diatomic molecules. J. Chem. Phys.; 2012; 137, 014101. [DOI: https://dx.doi.org/10.1063/1.4731340]

43. Hamzavi, M.; Ikhdair, S.M.; Thylwe, K.-E. Equivalence of the empirical shifted Deng–Fan oscillator potential for diatomic molecules. J. Math. Chem.; 2013; 51, 227. [DOI: https://dx.doi.org/10.1007/s10910-012-0075-x]

44. Edet, C.O.; Ikot, A.N.; Okorie, U.S.; Ramantswana, M.; Rampho, G.J.; Horchani, R.; Abdullah, H.Y.; Zahran, H.Y.; Obagboye, L.F.; Abdel-Aty, A.-H. Eigenfunctions, uncertainties and thermal properties of diatomic molecules under screened modified Kratzer potential. Indian J. Phys.; 2022; 1, pp. 1-20. [DOI: https://dx.doi.org/10.1007/s12648-022-02292-5]

45. Donfack, B.; Mbognou, F.C.; Tedondje, G.T.; Cedric, T.M.; Fotue, A.J. Cumulative Effects of Laser and Spin–Orbit Interaction (SOI) on the Thermal Properties of Quantum Pseudo-dot. J. Low Temp. Phys.; 2021; 206, pp. 63-79. [DOI: https://dx.doi.org/10.1007/s10909-021-02623-6]

46. Al, E.B.; Kasapoglu, E.; Sari, H.; Sokmen, I. Zeeman splitting, Zeeman transitions and optical absorption of an electron confined in spherical quantum dots under the magnetic field. Philos. Mag.; 2021; 101, 117. [DOI: https://dx.doi.org/10.1080/14786435.2020.1821112]

47. Acosta, R.E.; Morales, A.L.; Duque, C.M.; Mora-Ramos, M.E.; Duque, C.A. Optical absorption and refractive index changes in a semiconductor quantum ring: Electric field and donor impurity effects. Phys. Status Solidi (B); 2016; 253, 744. [DOI: https://dx.doi.org/10.1002/pssb.201552514]

48. Al, E.B. Effect of size modulation and donor position on intersubbands refractive index changes of a donor within a spherical core/shell/shell semiconductor quantum dot. Cumhur. Sci. J.; 2021; 42, 694. [DOI: https://dx.doi.org/10.17776/csj.927289]

49. M’zerd, S.; Haouari, M.E.; Talbi, A.; Feddi, E.; Mora-Ramos, M.E. Impact of electron-LO-phonon correction and donor impurity localization on the linear and nonlinear optical properties in spherical core/shell semiconductor quantum dots. J. Alloys Compd.; 2018; 753, 68. [DOI: https://dx.doi.org/10.1016/j.jallcom.2018.04.185]

50. Duque, C.M.; Mora-Ramos, M.E.; Duque, C.A. On-center donor impurity-related nonlinear corrections to optical absorption and refractive index in a two-dimensional quantum ring. Opt. Commun.; 2012; 285, 5456. [DOI: https://dx.doi.org/10.1016/j.optcom.2012.07.119]

51. Rastegar-Sedehi, H.R. Magnetocaloric effect in Rashba spin-orbit coupling and Zeeman splitting of a narrow nanowire quantum dot. Eur. Phys. J. Plus; 2021; 136, 514. [DOI: https://dx.doi.org/10.1140/epjp/s13360-021-01532-7]

52. Baghramyan, H.M.; Barseghyan, M.G.; Kirakosyan, A.A.; Laroze, D.; Duque, C.A. Donor-impurity related photoionization cross section in GaAs/Ga1^- xAlxAs concentric double quantum rings: Effects of geometry and hydrostatic pressure. Phys. B Condens. Matter; 2014; 449, 193. [DOI: https://dx.doi.org/10.1016/j.physb.2014.05.034]

53. Al, E.B.; Kasapoglu, E.; Sakiroglu, S.; Sari, H.; Sokmen, I.; Duque, C.A. Binding energies and optical absorption of donor impurities in spherical quantum dot under applied magnetic field. Phys. E Low-Dimens. Syst. Nanostruct.; 2020; 119, 114011. [DOI: https://dx.doi.org/10.1016/j.physe.2020.114011]

54. Ozmen, A.; Yakar, Y.; Çakır, B.; Atav, U. Computation of the oscillator strength and absorption coefficients for the intersubband transitions of the spherical quantum dot. Opt. Commun.; 2009; 282, 3999. [DOI: https://dx.doi.org/10.1016/j.optcom.2009.06.043]

55. Al, E.B.; Kasapoglu, E.; Sari, H.; Sokmen, I. Optical properties of spherical quantum dot in the presence of donor impurity under the magnetic field. Phys. B Condens. Matter; 2021; 613, 412874. [DOI: https://dx.doi.org/10.1016/j.physb.2021.412874]

56. Wu, S.; Wan, L. Electronic structures in a CdSe spherical quantum dot in a magnetic field: Diagonalization method and variational method. J. Appl. Phys.; 2012; 111, 063711. [DOI: https://dx.doi.org/10.1063/1.3695454]

57. Mathe, L.; Onyenegecha, C.P.; Farcas, A.-A.; Pioras-Ţimbolmas, L.-M.; Solaimani, M.; Hassanabadi, H. Linear and nonlinear optical properties in spherical quantum dots: Inversely quadratic Hellmann potential. Phys. Lett. A; 2021; 397, 127262. [DOI: https://dx.doi.org/10.1016/j.physleta.2021.127262]

58. Khordad, R.B. Application of Tietz potential to study optical properties of spherical quantum dots. Mirhosseini Pramana; 2015; 85, 723. [DOI: https://dx.doi.org/10.1007/s12043-014-0906-3]

59. Berkdemir, C. Application of the Nikiforov-Uvarov method in quantum mechanics. Theoretical Concepts of Quantum Mechanics; BoD—Books on Demand: Norderstedt, Germany, 2012; 225.