Table of Contents
ISRN Optics
VolumeΒ 2012, Article IDΒ 847532, 8 pages
Research Article

Impurity Modulated Static Linear and First Nonlinear Polarizabilities of Doped Quantum Dots

1Department of Physics, Suri Vidyasagar College, Suri, Birbhum 731101, West Bengal, India
2Department of Chemistry, Physical Chemistry Section, Visva Bharati University, Santiniketan, Birbhum 731235, West Bengal, India

Received 7 January 2012; Accepted 26 February 2012

Academic Editors: G.Β Bellanca, D.Β Cojoc, E.Β Lidorikis, and G.Β Ma

Copyright Β© 2012 Nirmal Kr Datta and Manas Ghosh. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We explore the pattern of linear and first nonlinear optical (NLO) response of repulsive impurity doped quantum dots harmonically confined in two dimensions. The dopant impurity potential chosen assumes Gaussian form. The quantum dot is subject to a static electric field. For some fixed values of transverse magnetic field strength (πœ”π‘), and harmonic confinement potential (πœ”0), the influence of impurity strength (𝑉0), impurity stretch (πœ‰), and impurity location (π‘Ÿ0) on the diagonal components of static linear (𝛼π‘₯π‘₯ and 𝛼𝑦𝑦), and the first NLO (𝛽π‘₯π‘₯π‘₯ and 𝛽𝑦𝑦𝑦) responses of the dot are computed through linear variational route. The investigation reveals the crucial roles played by the various impurity parameters in modulating these optical responses. Interestingly, maximization in the first NLO responses has been observed for some particular dopant location and at some particular value of spatial stretch of the dopant.

1. Introduction

The study of impurity states in low-dimensional heterostructures has emerged as an important aspect to which many theoretical and experimental works have been dedicated. As a result, nowadays the researches on doped semiconductor devices come out to be ubiquitous [1]. The doped system has quantized properties making them ideal objects for scientific study and robust applications. The dopants modify the chemical potential of a material. The control of optoelectronic properties of a wide range of semiconductor devices [2–4] now turns out to be the most fascinating aspect of impurity doping in such materials.

Miniaturization of semiconductor devices reaches the bottom of the avenue with the advent of so-called low-dimensional structures such as quantum dots (QDs), with QD’s new perspectives and delicacies in the field of impurity doping emerge, due to the mingling of new confinement sources with impurity related potentials [5]. Such confinement, coupled with the dopant location, can dramatically alter the electronic and optical properties of the system [6, 7].

For quite sometimes, search for molecules and materials with high linear and nonlinear susceptibilities and ultrafast response time has been pursued all over the world [8–14]. The targeted molecules are important in communication technology, data storage, optical switching, and so forth [9, 11, 15]. As a result of numerous studies, certain broad designing clues have emerged. Based on these guidelines, most often one tries to design donor-acceptor molecules with large charge transfer from the donor to the acceptor moiety. NLO response is maximized for certain optimal combinations of charge transfer and hopping interactions and the length over which the charge is transferred [16–19]. Modulation in the molecular electronic parameters brought through subtle structural changes or change in substitution changes the electronic wave function and therefore the electron density distribution which ultimately shapes all the properties of atoms and molecules. It would be crucial to have at one’s disposal, systems, the electronic structures of which could be continuously tuned by adjusting suitable control parameters. In molecules, these variations are only discrete, having caused by either a change of substituent or by structural alteration or by both. Quantum dots are representatives of the systems where the electronic structural dispositions can be practically continuously varied as functions of a few of the system parameters. From the above discussion, it appears inevitable that variation in the electronic structure of the dot induced by the impurity could subtly shape all its properties. It is thus crucial to study the role of dopant impurity on the optical properties of the quantum dots. A thorough expedition through literature reveals some notable works on the optical properties of impurity doped quantum dots [20, 21]. Doped quantum dots, thus, beyond any doubt, possess the potential to exhibit high linear and nonlinear optical response properties and would thus be principally applied in the area of all-optical signal processing. A detailed investigation of various dopant parameters relevant to optical signal processing would continue to be a topic of active research. In what follows, of late, we have also studied the frequency dependent linear and nonlinear polarizabilities of doped quantum dots as functions of various dopant parameters [22, 23].

In the present manuscript, we inspect the diagonal components of linear (𝛼π‘₯π‘₯ and 𝛼𝑦𝑦) and first nonlinear polarizabilities (𝛽π‘₯π‘₯π‘₯ and 𝛽𝑦𝑦𝑦) of a repulsive impurity doped quantum dot subject to a static electric field. Following earlier works on the effects of a repulsive scatterer in multicarrier dots in the presence of magnetic field [24, 25], here we have considered that the QD is doped with a repulsive Gaussian impurity. When the impurity is doped in an on-center location, it does not destroy the inversion symmetry of the dot and consequently the emergence of 𝛽 value is ruled out. However, nonzero 𝛼 values are envisaged. At off-center dopant locations, however, emergence of both nonzero 𝛼 and 𝛽 values is observed owing to the destruction of inversion symmetry. We have found that, in conjunction with the dopant location, the strength and the spatial stretch of the dopant also affect the polarizability values quite prominently.

2. Method

We consider the energy eigenstates of an electron subject to a harmonic confinement potential 𝑉(π‘₯,𝑦) and a perpendicular magnetic field 𝐡, where 𝑉(π‘₯,𝑦)=(1/2)π‘šβˆ—πœ”20(π‘₯2+𝑦2), πœ”π‘=𝑒𝐡/π‘šβˆ—π‘, and Landau gauge [𝐴=(𝐡𝑦,0,0)] have been used. πœ”0, πœ”π‘, and 𝐴 stand for harmonic confinement potential, cyclotron frequency (a measure of magnetic confinement offered by 𝐡), and vector potential, respectively. The Hamiltonian in our problem readsπ»ξ…ž0ℏ=βˆ’22π‘šβˆ—ξ‚΅πœ•2πœ•π‘₯2+πœ•2πœ•π‘¦2ξ‚Ά+12π‘šβˆ—πœ”20π‘₯2,+12π‘šβˆ—ξ€·πœ”20+πœ”2𝑐𝑦2βˆ’π‘–β„πœ”π‘π‘¦πœ•.πœ•π‘₯(1) Define Ξ©2=πœ”20+πœ”2𝑐 as the effective frequency in the 𝑦-direction. In real QDs, the electrons are confined in 3 dimensions that is, the carriers are dynamically confined to zero dimensions. The confinement length scales 𝑅1, 𝑅2, and 𝑅3 can be different in three spatial directions, but typically 𝑅3β‰ͺ𝑅1≃𝑅2≃100 nm. In models of such dots, 𝑅3 is often taken to be strictly zero and the confinement in the other two directions is described by a potential 𝑉 with 𝑉(π‘₯)β†’βˆž for |π‘₯|β†’βˆž, π‘₯=(π‘₯1,π‘₯2)βˆˆπ‘…2. A parabolic potential, 𝑉=(1/2)πœ”|π‘₯|2, is often used as a realistic and at the same time computationally convenient approximation. Assuming that the 𝑧-extension could be effectively considered zero, the electronic properties in these nanostructures have been successfully described within the model of the 1-electron motion in 2-d harmonic oscillator potential in the presence of a magnetic field [26–28]. Now, intrusion of impurity perturbation transforms the Hamiltonian to𝐻0ξ€·π‘₯,𝑦,πœ”π‘,πœ”0ξ€Έ=π»ξ…ž0ξ€·π‘₯,𝑦,πœ”π‘,πœ”0ξ€Έ+𝑉impξ€·π‘₯0,𝑦0ξ€Έ,(2) where 𝑉imp(π‘₯0,𝑦0)=𝑉imp(0)=𝑉0π‘’βˆ’πœ‰[(π‘₯βˆ’π‘₯0)2+(π‘¦βˆ’π‘¦0)2] with πœ‰>0 and 𝑉0>0 for repulsive impurity, and (π‘₯0,𝑦0) denotes the position of the impurity center. In QD’s, the electrons visit an almost equipotential surface wherein they are set free. These electrons are carrier electrons. In reality, regular dot occurs rarely as the deformation of the boundary of QD’s modifies the effective potential of the carrier electrons to bunch of uneven pockets. Sometimes, these electrons play the role of defects themselves (repulsive impurity) [29]. 𝑉0 is a measure of the strength of impurity potential whereas πœ‰ determines the extent of influence of the impurity potential. A large value of πœ‰ indicates that the spatial extension of impurity potential is highly restricted whereas a small πœ‰ accounts for spatially diffused one. Thus, a change in πœ‰ in turn causes a change in the extent of dot-impurity overlap that affects the excitation pattern noticeably [22, 23]. The presence of repulsive scatterer simulates dopant with excess electrons. The choice of Gaussian impurity potential is not arbitrary as it has been exploited by several investigators [30–35].

We write the trial wave function πœ“(π‘₯,𝑦) as a superposition of the product of harmonic oscillator eigenfunctions πœ™π‘›(𝛼π‘₯) and πœ™π‘š(𝛽𝑦), respectively, as follows:ξ“πœ“(π‘₯,𝑦)=𝑛,π‘šπΆπ‘›,π‘šπœ™π‘›(𝛼π‘₯)πœ™π‘š(𝛽𝑦),(3) where 𝐢𝑛,π‘š are the variational parameters and βˆšπ›Ό=π‘šβˆ—πœ”0/ℏ and βˆšπ›½=π‘šβˆ—Ξ©/ℏ. The matrix elements of 𝑉imp are given by𝑉imp𝑛,π‘š;𝑛′,π‘šβ€²=𝑉0ξ‚¬πœ™π‘›(𝛼π‘₯)πœ™π‘š(|||𝑒𝛽𝑦)βˆ’πœ‰[(π‘₯βˆ’π‘₯0)2+(π‘¦βˆ’π‘¦0)2]|||πœ™π‘›β€²(𝛼π‘₯)πœ™π‘šβ€²(𝛽𝑦)=𝑉0𝐼1𝐼2,(4) where𝐼1=ξ‚¬πœ™π‘›|||𝑒(𝛼π‘₯)βˆ’πœ‰(π‘₯βˆ’π‘₯0)2|||πœ™ξ…žπ‘›ξ‚­ξ€œ(𝛼π‘₯)=𝐴+βˆžβˆ’βˆžπ»π‘›(𝛼π‘₯)𝐻𝑛′(𝛼π‘₯)π‘’βˆ’π›Ό2π‘₯2π‘’βˆ’πœ‰(π‘₯βˆ’π‘₯0)2𝐼𝑑π‘₯,2=ξ‚¬πœ™π‘š(|||𝑒𝛽𝑦)βˆ’πœ‰(π‘¦βˆ’π‘¦0)2|||πœ™ξ…žπ‘šξ‚­ξ€œ(𝛽𝑦)=𝐡+βˆžβˆ’βˆžπ»π‘š(𝛽𝑦)π»π‘šβ€²(𝛽𝑦)π‘’βˆ’π›½2𝑦2π‘’βˆ’πœ‰(π‘¦βˆ’π‘¦0)2𝑑𝑦,(5) with 𝐴=𝛼/(2𝑛+𝑛′𝑛!π‘›ξ…ž!πœ‹)1/2 and 𝐡=𝛽/(2π‘š+π‘šβ€²π‘š!π‘šξ…ž!πœ‹)1/2. With the transformations 𝛿21=𝛼2+πœ‰, 𝛿22=𝛽2+πœ‰, πœ†1=exp[βˆ’(πœ‰π‘₯20(𝛿21βˆ’πœ‰))/𝛿21], and πœ†2=exp[βˆ’(πœ‰π‘¦20(𝛿22βˆ’πœ‰))/𝛿22], one can write𝐼1=π΄πœ†1ξ€œ+βˆžβˆ’βˆžπ»π‘›ξ€·π›Όβˆ—π‘’ξ€Έπ»π‘›β€²ξ€·π›Όβˆ—π‘’ξ€Έπ‘’βˆ’(π‘’βˆ’πœŒ1)2𝐼𝑑𝑒,2=π΅πœ†2ξ€œ+βˆžβˆ’βˆžπ»π‘šξ€·π›½βˆ—π‘£ξ€Έπ»π‘šβ€²ξ€·π›½βˆ—π‘£ξ€Έπ‘’βˆ’(π‘£βˆ’πœŒ2)2𝑑𝑣,(6) where, 𝜌1=πœ‰π‘₯0/𝛿1, 𝜌2=πœ‰π‘¦0/𝛿2, π›Όβˆ—=𝛼/𝛿1, π›½βˆ—=𝛽/𝛿2, 𝑒=𝛿1π‘₯, and 𝑣=𝛿2𝑦. With the help of the standard integral [36], it is now easy to write𝐼1=𝐷1min(𝑛,𝑛′)βˆ‘π‘˜=0π‘“ξ€·π‘˜,𝑛,π‘›ξ…žξ€Έ,𝐼2=𝐷2min(π‘š,π‘šβ€²)βˆ‘π‘™=0𝑔𝑙,π‘š,π‘šξ…žξ€Έ,(7) whereπ‘“ξ€·π‘˜,𝑛,π‘›ξ…žξ€Έ=2π‘˜β‹…π‘˜!β‹…π‘›πΆπ‘˜β‹…π‘›β€²πΆπ‘˜β‹…(1βˆ’π›Όβˆ—2)(𝑛+𝑛′)/2βˆ’π‘˜β‹…π»π‘›+π‘›β€²βˆ’2π‘˜ξ€·π›Ό1𝜌1ξ€Έ,𝑔𝑙,π‘š,π‘šξ…žξ€Έ=2𝑙⋅𝑙!β‹…π‘šπΆπ‘™β‹…π‘šβ€²πΆπ‘™β‹…(1βˆ’π›½βˆ—2)(π‘š+π‘šβ€²)/2βˆ’π‘™β‹…π»π‘š+π‘šβ€²βˆ’2𝑙𝛽1𝜌2ξ€Έ,(8) with 𝐷1=π΄πœ†1πœ‹1/2/𝛿1 and 𝐷2=π΅πœ†2πœ‹1/2/𝛿2. Thus, finally, we obtain(𝑉imp)𝑛,π‘š;𝑛′,π‘šβ€²=𝑉0⋅𝐷1⋅𝐷2β‹…min(𝑛,𝑛′)ξ“π‘˜=0min(π‘š,π‘šβ€²)𝑙=0π‘“ξ€·π‘˜,𝑛,π‘›ξ…žξ€Έξ€·β‹…π‘”π‘™,π‘š,π‘šξ…žξ€Έ.(9)𝐻𝑛(π‘₯) stands for the Hermite polynomials of 𝑛th order. The 𝑝th eigenstate of the system in this representation can be written asπœ“π‘(π‘₯,𝑦)=𝑖𝑗𝐢𝑖𝑗,π‘ξ€½πœ™π‘–(𝛼π‘₯)πœ™π‘—(ξ€Ύ,𝛽𝑦)(10) where 𝑖, 𝑗 are the appropriate quantum numbers, respectively, and (𝑖𝑗) are composite indices specifying the direct product basis.

In presence of a static electric field of strength πœ€, the effective Hamiltonian of the system (𝐻eff) becomes𝐻eff=𝐻0βˆ’ξ€½πœ€π‘₯|𝑒|π‘₯+πœ€π‘¦ξ€Ύ|𝑒|𝑦.(11) We have determined the values of ⟨𝐸(Β±πœ€π‘₯,πœ€π‘¦=0)⟩, ⟨𝐸(Β±2πœ€π‘₯,πœ€π‘¦=0)⟩, ⟨𝐸(πœ€π‘₯=0,Β±πœ€π‘¦)⟩, and ⟨𝐸(πœ€π‘₯=0,Β±2πœ€π‘¦)⟩ and we used the data to compute the direct components of polarizabilities of the dot by the following relations obtained by numerical differentiation [37–41]:𝛼π‘₯π‘₯πœ€2π‘₯=524⟨𝐸(0)βŸ©βˆ’3ξ€Ίξ€·πœ€βŸ¨πΈπ‘₯ξ€Έξ€·βŸ©+βŸ¨πΈβˆ’πœ€π‘₯ξ€ΈβŸ©ξ€»+1ξ€Ίξ€·12⟨𝐸2πœ€π‘₯ξ€Έξ€·βŸ©+βŸ¨πΈβˆ’2πœ€π‘₯ξ€ΈβŸ©ξ€»,(12) and a similar expression for π›Όπ‘¦π‘¦πœ€2𝑦:𝛽π‘₯π‘₯π‘₯πœ€3π‘₯=ξ€Ίξ€·πœ€βŸ¨πΈπ‘₯ξ€Έξ€·,0βŸ©βˆ’βŸ¨πΈβˆ’πœ€π‘₯ξ€ΈβŸ©ξ€»βˆ’1,02ξ€Ίξ€·βŸ¨πΈ2πœ€π‘₯ξ€Έξ€·,0βŸ©βˆ’βŸ¨πΈβˆ’2πœ€π‘₯ξ€ΈβŸ©ξ€»,,0(13) and a similar expression is used for computing 𝛽𝑦𝑦𝑦 component. In these expressions, ⟨𝐸(Β±πœ–π‘₯)⟩ represents average energy of the system when πœ€π‘¦=0 and πœ€π‘₯β‰ 0. Analogously, we have also written ⟨𝐸(Β±πœ–π‘¦)⟩. The above relations indicate different combinations of electric field intensities and orientations.

3. Results and Discussion

3.1. System Parameters

The model Hamiltonian (cf. (1)) can be made to represent a 2-d quantum dot with a single carrier electron [28, 42]. The form of the confinement potential indicates lateral electrostatic confinement of the electrons in the π‘₯-𝑦 plane. Thus, it is disklike which is normally fabricated by molecular beam epitaxy (MBE) technique. Because of unequal confinements in π‘₯ and 𝑦-directions, the symmetry is nonhomogeneous. π‘šβˆ— is the effective electronic mass appropriate for describing the motion of the electrons within the lattice of the material to be used. We have used π‘šβˆ—=0.5π‘š0 and set ℏ=𝑒=π‘š0=π‘Ž0=1. The radial position of dopant (π‘Ÿ0) has been varied from 0.0 a.u. (on-center) to 70.71 a.u. (off-center) positions. In the linear variational calculation, we have used basis functions (cf. (3)) with 𝑛,π‘š=0βˆ’20 for each of the directions (π‘₯,𝑦). The direct product basis spans a space of (21Γ—21) dimension. We have checked that the basis functions span the 2-d space effectively completely, at least with respect to representing the observables under investigation. We have made the convergence test with still greater number of basis functions.

We have made some attempt to reasonably connect our theoretical parameters to the real life-doped QD. The parameter πœ‰ in the impurity potential can be correlated to 1/𝑑2, where 𝑑 is proportional to the width of the impurity potential [24, 25]. Table 1 gives the 𝑑 values in nm corresponding to different πœ‰ values to have a feelings of the actual extension of the impurity domain. The π‘šβˆ— value that we have used (π‘šβˆ—=0.5 a.u.) closely resembles Ge quantum dots (π‘šβˆ—=0.55 a.u.). The maximum value of the dopant strength (𝑉0) was limited to ~10βˆ’3 a.u. or 27.2meV and the applied magnetic field is of the order of milliTesla (mT). Here, we want to mention that our method is not strictly a perturbation theory in that sense where some perturbation is added to the original Hamiltonian through the agency of a perturbation parameter. In our work, we have exploited linear variation principle to determine the approximate eigenstates of the system as it was not possible to solve the time-independent Schrodinger equation containing impurity in the Hamiltonian as well as in the presence of static electric field. Thus, firstly, we have evaluated the approximate eigenstates of the system taking harmonic oscillator eigenfunctions as the basis function. This involves modulation of the linear combination coefficients. Obviously, the modulation was not arbitrary but was governed by the diagonalization of the Hamiltonian matrix containing the contributions from impurity potential and static electric field. Thus, the method is actually a basis function expansion technique with properly adjusted linear combination coefficients. The validity of the method lies in the fact that how far it can endure the increase in 𝑉0 values so that normalization is not hampered. We did not make any forceful normalization of the wave function and checked that inherent normalization is well maintained even with a 𝑉0 value as large as 1.0Γ—10βˆ’3 a.u. for a very large basis size (21Γ—21) indicating that our method is quite robust so far as stability of the method is concerned. The static electric field has an intensity of 2.57Γ—104 Vmβˆ’1. We believe that these values might give some realization of real systems containing actual impurity.

Table 1: The width of the impurity domain in nm for different πœ‰ values.
3.2. Aspects of Polarizability Components
3.2.1. Influence of Impurity Location

Figure 1 depicts the variation of components of linear polarizability (𝛼π‘₯π‘₯ and 𝛼𝑦𝑦) in π‘₯ and 𝑦 directions as functions of radial position (π‘Ÿ0) of the dopant. Both the components are found to decrease monotonically with increase in π‘Ÿ0. An increase in π‘Ÿ0 takes the dopant away from the dot confinement and consequently there occurs a reduction in dot-impurity overlap. As a result, the extent of repulsive interaction also diminishes which is reflected in the observed behaviors of linear polarizability components. From the plot, it has also been observed that the 𝛼𝑦𝑦 value is lower than 𝛼π‘₯π‘₯ component. This could be due to varied extent of confinement in π‘₯ and 𝑦 directions. A close look at (1) reveals that the effective harmonic frequency along 𝑦-axis is Ω=πœ”20+πœ”2𝑐 while that along π‘₯-axis is πœ”0. Thus, the 𝑦-direction is under stronger confinement in comparison to the π‘₯-direction and naturally enjoys less dispersive character. Consequently, the linear polarizability components which are intimately connected with the dispersive nature of the system assume different values in π‘₯ and 𝑦 directions. It needs to be mentioned that macroscopic polarizability should be isotropic for a symmetric QD in a perpendicular magnetic field. In our investigation, we have found anisotropy in the values of 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 on a microscopic scale because of a nonzero πœ”π‘. We feel that on a macroscopic scale, the anisotropies arising out of various microscopic components may balance each other leading to no net anisotropy of the system. However, a more extensive and exhaustive computation is required for an in depth investigation.

Figure 1: Plots of 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 versus π‘Ÿ0 with 𝑉0=1.0Γ—10βˆ’6 a.u. and πœ‰=0.001 a.u. (i) for 𝛼π‘₯π‘₯ and (ii) for 𝛼𝑦𝑦.

We now investigate the diagonal components of first hyperpolarizability (𝛽π‘₯π‘₯π‘₯ and 𝛽𝑦𝑦𝑦) of the doped dot. For the emergence of nonzero (𝛽), doping of impurity at off-center position is essential. An impurity doping at on-center position does not destroy the inversion symmetry of the Hamiltonian and is therefore unable to generate non-zero 𝛽 value. Both the 𝛽 components exhibit maximization at a particular dopant location as the impurity is shifted away from the dot center (Figure 2). A shift of the dopant away from the dot confinement center is endowed with opposite consequences. On one hand, the intensity of dot confinement certainly decreases thereby favoring the high optical response. On the other hand, the said shift also decreases the extent of dot-impurity repulsive interaction and discourages emergence of high optical response. So far as generation of 𝛽 values are concerned; it appears that the interplay of these two opposing factors are responsible for the aforesaid maximization. For high NLO response, therefore, the designer quantum dots appear to be not hard to find. The requirement being the controlled incorporation of the dopant at definite site to break the symmetry of the confinement potential. In this regard, we should mention that of late there are some excellent experiments which show the mechanism of dopant incorporation [43] and how such incorporations can be controlled [44].

Figure 2: Plots of |𝛽π‘₯π‘₯π‘₯| and |𝛽𝑦𝑦𝑦| versus π‘Ÿ0 with 𝑉0=1.0Γ—10βˆ’6 a.u. and πœ‰=0.001 a.u. (i) for |𝛽π‘₯π‘₯π‘₯| and (ii) for |𝛽𝑦𝑦𝑦|.
3.2.2. Influence of Impurity Strength

We have now varied 𝑉0 over a wide range to observe its influence on the polarizability components. Figure 3 exhibits the 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 profiles as function of dopant strength for on and two off-center (π‘Ÿ0=42.43 a.u. and 70.71 a.u.) dopant locations. The plots reveal how a variation in the dopant strength can influence the polarizability components. We notice that at all dopant locations both the linear polarizability components exhibit a saturation at high dopant strength. Within extremely small dopant strength regime (𝑉0β†’0), there occurs an initial uprise in the polarizability values with increase in 𝑉0. It seems that in this regime the dot-impurity interaction is very feeble. Small increase in dopant strength in this domain primarily causes some sort of development in the dot-impurity repulsive interaction and we observe an initial surge in the polarizability values. At large 𝑉0 the plots show some kind of saturation indicating a steady extent of dot-impurity interaction. As observed earlier now also we find that 𝛼𝑦𝑦 values are somewhat lower in magnitude than 𝛼π‘₯π‘₯ values. Moreover, from the plots, we find that while the values of a particular component for on and near off-center dopant locations (π‘Ÿ0=42.43 a.u.) are in close proximity, the value drops substantially at far off-center dopant location (π‘Ÿ0=70.71 a.u.) owing to a diminished dot-impurity interaction.

Figure 3: Plots of 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 versus 𝑉0 with πœ‰=0.001 a.u. at on and two off-center dopant locations: (i) for 𝛼π‘₯π‘₯ at π‘Ÿ0=0.0 a.u., (ii) for 𝛼π‘₯π‘₯ at π‘Ÿ0=42.43 a.u., (iii) for 𝛼π‘₯π‘₯ at π‘Ÿ0=70.71 a.u., (iv) for 𝛼𝑦𝑦 at π‘Ÿ0=0.0 a.u., (v) for 𝛼𝑦𝑦 at π‘Ÿ0=42.43 a.u., and (vi) for 𝛼𝑦𝑦 at π‘Ÿ0=70.71 a.u.

A similar saturation behavior is also envisaged in |𝛽π‘₯π‘₯π‘₯| and |𝛽𝑦𝑦𝑦| profiles as a function of 𝑉0 (Figure 4). Since the |𝛽| components assume zero value at on-center dopant location the figure depicts their profiles at near off-center (π‘Ÿ0=42.43 a.u.) and far off-center (π‘Ÿ0=70.71 a.u.) dopant locations. Here, also we notice a lower magnitude of |𝛽𝑦𝑦𝑦| component in comparison to its counterpart along π‘₯-direction. Interestingly, each of the |𝛽| component exhibits greater magnitude in the far off-center dopant location in comparison to that of near off-center position. The behavior happens to be quite obvious as already we have come across position dependent maximization in |𝛽| components.

Figure 4: Plots of |𝛽π‘₯π‘₯π‘₯| and |𝛽𝑦𝑦𝑦| versus 𝑉0 with πœ‰=0.001 a.u. at two off-center dopant locations: (i) for |𝛽π‘₯π‘₯π‘₯| at π‘Ÿ0=42.43 a.u., (ii) for |𝛽π‘₯π‘₯π‘₯| at π‘Ÿ0=70.71 a.u., (iii) for |𝛽𝑦𝑦𝑦| at π‘Ÿ0=42.43 a.u., and (iv) for |𝛽𝑦𝑦𝑦| at π‘Ÿ0=70.71 a.u.
3.2.3. Influence of Impurity Spread

We now turn our attention towards inspecting the influence of spatial stretch of impurity (πœ‰) on the polarizability components for on and two off-center (π‘Ÿ0=42.43 a.u. and 70.71 a.u.) dopant locations (Figure 5). The 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 components initially delineate sort of decrease with increase in πœ‰ (in small πœ‰ domain) but ultimately culminate in some steady value with further increase in πœ‰ value. A small πœ‰ value implies that the impurity potential is diffused over a long spatial region, so also its influence. In the very low πœ‰ regime, this extreme diffusive nature of the impurity domain results in substantial dot-impurity interaction and consequently a large dispersive nature of the system. The polarizability components thus register somewhat large value. The scenario changes completely in the high πœ‰ regime. Now, a spatially quenched impurity potential undergoes very weak overlap with the dot confinement center resulting in small dispersive nature of the system. This causes a fall in the linear polarizability values. At very high value of πœ‰, the spatial stretch of impurity becomes highly condensed and the decreasing dispersive nature of the system reaches its limit. The outcome being the steady values of the linear polarizability components. Also, as before, the 𝛼𝑦𝑦 components are smaller in magnitude in comparison to their π‘₯-directional counterpart. Furthermore, for each components, the extent of fall in its magnitude with increase in πœ‰ becomes progressively more pronounced with shift of the dopant from on- to off-center location. This is simply because of the fact that since a high πœ‰ value makes dot-impurity overlap already weak, a shift of dopant to more off-center location further aggravates the situation and hastens the fall of polarizability. Peculiarly, for the intermediate off-center impurity positions, the curve is not monotonous. It seems difficult to understand why in this case we do not encounter an otherwise straightforward curve. It appears that at intermediate dopant locations, along with dot-impurity overlap, the extent of dot confinement also plays a role in shaping the polarizability components. At these locations, a reduced dot-impurity overlap decreases the polarizability values. On the other hand, because of considerable dot-impurity separation the dot confinement also becomes less stringent making the system more flexible which enhances the polarizability values. The competitive behavior between these two opposite influences could be the cause behind the departure of the aforesaid profile from an otherwise monotonous plot.

Figure 5: Plots of 𝛼π‘₯π‘₯ and 𝛼𝑦𝑦 versus πœ‰ with 𝑉0=1.0Γ—10βˆ’6 a.u. at on and two off-center dopant locations: (i) for 𝛼π‘₯π‘₯ at π‘Ÿ0=0.0 a.u., (ii) for 𝛼𝑦𝑦 at π‘Ÿ0=0.0 a.u., (iii) for 𝛼π‘₯π‘₯ at π‘Ÿ0=42.43 a.u., (iv) for 𝛼𝑦𝑦 at π‘Ÿ0=42.43 a.u., (v) for 𝛼π‘₯π‘₯ at π‘Ÿ0=70.71 a.u., and (vi) for 𝛼𝑦𝑦 at π‘Ÿ0=70.71 a.u.

|𝛽π‘₯π‘₯π‘₯| and |𝛽𝑦𝑦𝑦| components, on the other hand, exhibit prominent maximization as a function of πœ‰ (Figure 6). Here, we present the corresponding profiles at near off-center (π‘Ÿ0=42.43 a.u.) and far off-center (π‘Ÿ0=70.71 a.u.) dopant locations and as usual the |𝛽| components assume zero value at on-center dopant location. The continual decrease in dopant’s spatial stretch with increase in πœ‰ value might have some contrasting consequences on |𝛽| components. Although it reduces the confines of dot potential on the dopant but simultaneously the dot-impurity repulsive interaction is also diminished. We feel that a changeover in the mutual dominance of these two contrasting factors could give rise to maximization of |𝛽𝑦𝑦𝑦| components at some particular πœ‰ value. As expected, here also we notice a lower magnitude of |𝛽𝑦𝑦𝑦| component in comparison to its counterpart along π‘₯-direction. Interestingly, each of the |𝛽| component exhibits greater magnitude in the far off-center dopant location in comparison to that of near off-center position. The behavior happens to be quite obvious as already we have visualized position dependent maximization in |𝛽| components.

Figure 6: Plots of |𝛽π‘₯π‘₯π‘₯| and |𝛽𝑦𝑦𝑦| versus πœ‰ with 𝑉0=1.0Γ—10βˆ’6 a.u. at two off-center dopant locations: (i) for |𝛽π‘₯π‘₯π‘₯| at π‘Ÿ0=42.43 a.u., (ii) for |𝛽π‘₯π‘₯π‘₯| at π‘Ÿ0=70.71 a.u., (iii) for |𝛽𝑦𝑦𝑦| at π‘Ÿ0=42.43 a.u., and (iv) for |𝛽𝑦𝑦𝑦| at π‘Ÿ0=70.71 a.u.

4. Conclusions

The diagonal components of linear and first nonlinear polarizabilities of repulsive impurity doped quantum dots subject to a static electric field reveal intriguing features. For an in-depth analysis, we examine the roles played by impurity spread, impurity strength and most importantly impurity location meticulously on these components. We have found that the linear polarizability components decrease with increase in separation of dopant location from that of dot confinement center and also with an increase in spatial shrinkage of dopant potential. On the other hand, we envisage an increase in the said components with impurity strength. However, in all cases finally some steady behavior has been observed and 𝑦 direction has been found to be under more stringent confinement than the π‘₯ direction. The first nonlinear polarizability components evince maximization at some particular dopant location and for some particular value of spatial spread of dopant in absolute sense. The maximization appears to be due to the conflict between two opposing factors that foment and hinder the effective confinement of the dot on the dopant. However, as a function of dopant strength, we observe a persistent increase in |𝛽| values in absolute sense which culminates in a steady magnitude. We expect that the results obtained could have important implications in optical applications of quantum dot nanodevices.

At the fag end of the discussion, it appears to be quite significant to highlight the new findings in the present investigation in the light of the results of [23]. First of all, we want to mention that in the present study we have explored the role of impurity location rigorously which was absent from our earlier study on frequency-dependent polarizability (FDP) [23]. Secondly, in [23], we observed that in the limitingly small dopant strength domain an increase in dopant strength causes an initial lull in the FDP values. However, as the dopant strength is increased as bit FDP values begin to increase smoothly with dopant strength culminating in a steady value. In the present study, although we have found some kind of similar trend at high 𝑉0 limit, but behavior is absolutely different in small 𝑉0 domain. The role of dopant’s spatial spread comes out to be quite distinct in the two investigations. In the earlier study, we found that squeezing the spatial expansion of the dopant domain results in persistent increase in FDP values whereas in the present enquiry a reverse behavior has been envisaged. We hope to explore a rigorous investigation on the frequency-dependent and static polarizability components on a comparative basis in near future.

In the present investigation, we did not consider the influence of size on the optical properties. Although in principle the dot wave function can stretch up to infinity but in practice, it actually terminates at some finite values. Thus, the size effect would be important at length scales within the actual termination of wave function.

It is quite expected that donor and acceptor impurity would exhibit distinct impacts on the NLO properties. Recently, Hazra et al. have investigated the role of donor and acceptor impurities in a slightly different context. It needs further study to precisely understand the distinct roles of acceptor and donor impurity pertinent to the present investigation.


The authors N. K. Datta and M. Ghosh thank D. S. T.-F. I. S. T (Government of India) and U. G. C.-S. A. P (Government of India) for financial support.


  1. P. M. Koenraad and M. E. FlattΓ©, β€œSingle dopants in semiconductors,” Nature Materials, vol. 10, pp. 91–100, 2011. View at Publisher Β· View at Google Scholar
  2. H. J. Queisser and E. E. Haller, β€œDefects in semiconductors: some fatal, some vital,” Science, vol. 281, no. 5379, pp. 945–950, 1998. View at Google Scholar Β· View at Scopus
  3. B. Γ‡akir, Y. Yakar, A. Γ–zmen, M. Γ–. Sezer, and M. Şahin, β€œLinear and nonlinear optical absorption coefficients and binding energy of a spherical quantum dot,” Superlattices and Microstructures, vol. 47, no. 4, pp. 556–566, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  4. Y. Yakar, B. Γ‡akir, and A. Γ–zmen, β€œCalculation of linear and nonlinear optical absorption coefficients of a spherical quantum dot with parabolic potential,” Optics Communications, vol. 283, no. 9, pp. 1795–1800, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  5. J. L. Movilla and J. Planelles, β€œOff-centering of hydrogenic impurities in quantum dots,” Physical Review B, vol. 71, no. 7, Article ID 075319, pp. 1–7, 2005. View at Publisher Β· View at Google Scholar Β· View at Scopus
  6. M. J. Kelly, Low-Dimensional Semiconductors, Oxford University Press, Oxford, UK, 1995.
  7. C. P. Poole Jr and F. J. Owens, Introduction to Nanotechnology, Wiley, New York, NY, USA, 2003.
  8. J. L. BrΓ©das, C. Adant, P. Tackx, A. Persoons, and B. M. Pierce, β€œThird-order nonlinear optical response in organic materials: theoretical and experimental aspects,” Chemical Reviews, vol. 94, no. 1, pp. 243–278, 1994. View at Google Scholar Β· View at Scopus
  9. P. N. Prasad and D. J. Williams, Introduction to Non-linear Optical Effects in Molecules and Polymers, Wiley, New York, NY, USA, 1991.
  10. G. I. Stegemen, A. Miller, and J. E. Midwinter, Eds., Photonics in Switching volume 1 Background and Components, Academic Press, Boston, Mass, USA, 1993.
  11. D. S. Chemla and J. Zyss, Eds., Non-linear Optical Properties of Molecules and Crystals, vol. 1, Academic Press, New York, NY, USA, 1983.
  12. A. Leclercq, E. Zojer, S. H. Jang et al., β€œQuantum-chemical investigation of second-order nonlinear optical chromophores: comparison of strong nitrile-based acceptor end groups and role of auxiliary donors and acceptors,” Journal of Chemical Physics, vol. 124, no. 4, Article ID 044510, pp. 1–7, 2006. View at Publisher Β· View at Google Scholar Β· View at Scopus
  13. Z.-Y. Hu, A. Fort, M. Barzoukas, A. K. Y. Jen, S. Barlow, and S. R. Marder, β€œTrends in optical nonlinearity and thermal stability in electrooptic chromophores based upon the 3-(dicyanomethylene)-2,3-dihydrobenzothiophene-1, 1-dioxide acceptor,” Journal of Physical Chemistry B, vol. 108, no. 25, pp. 8626–8630, 2004. View at Publisher Β· View at Google Scholar Β· View at Scopus
  14. G. Ramos-Ortiz, M. Cha, S. Thayumanavan, J. C. Mendez, S. R. Marder, and B. Kippelen, β€œThird-order optical autocorrelator for time-domain operation at telecommunication wavelengths,” Applied Physics Letters, vol. 85, no. 2, pp. 179–181, 2004. View at Publisher Β· View at Google Scholar Β· View at Scopus
  15. G. J. Ashwell, R. C. Hargreaves, C. E. Baldwin, G. S. Bahra, and C. R. Brown, β€œImproved second-harmonic generation from Langmuir-Blodgett films of hemicyanine dyes,” Nature, vol. 357, no. 6377, pp. 393–395, 1992. View at Google Scholar Β· View at Scopus
  16. S. R. Marder, J. E. Sohn, and G. D. Stucky, Eds., Materials for Non-Linear Optics: Chemical Perspectives, ACS Sympos, vol. 455, American Chemical Society, Washington, DC, USA, 1991.
  17. R. Sen, D. Majumdar, S. P. Bhattacharyya, and S. N. Bhattacharyya, β€œModeling hyperpolarizabilities of some TICT molecules and their analogues,” Journal of Physical Chemistry, vol. 97, no. 29, pp. 7491–7498, 1993. View at Google Scholar
  18. A. D. Buckingham, E. P. Concannon, and I. D. Hands, β€œHyperpolarizability of interacting atoms,” Journal of Physical Chemistry, vol. 98, no. 41, pp. 10455–10459, 1994. View at Google Scholar Β· View at Scopus
  19. K. L. C. Hunt and A. D. Buckingham, β€œThe polarizability of H2 in the triplet state,” The Journal of Chemical Physics, vol. 72, no. 4, pp. 2832–2840, 1980. View at Google Scholar Β· View at Scopus
  20. F. Furtmayr, M. Vielemeyer, M. Stutzmann, A. Laufer, B. K. Meyer, and M. Eickhoff, β€œOptical properties of Si- and Mg-doped gallium nitride nanowires grown by plasma-assisted molecular beam epitaxy,” Journal of Applied Physics, vol. 104, no. 7, Article ID 074309, 2008. View at Publisher Β· View at Google Scholar Β· View at Scopus
  21. O. Weidemann, P. K. Kandaswamy, E. Monroy, G. Jegert, M. Stutzmann, and M. Eickhoff, β€œGaN quantum dots as optical transducers for chemical sensors,” Applied Physics Letters, vol. 94, no. 11, Article ID 113108, 2009. View at Publisher Β· View at Google Scholar Β· View at Scopus
  22. K. Sarkar, N. Kumar Datta, and M. Ghosh, β€œFrequency dependent linear and non-linear response properties of electron impurity doped quantum dots: influence of impurity location,” Physica E, vol. 42, no. 5, pp. 1659–1666, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  23. N. K. Datta and M. Ghosh, β€œImpurity strength and impurity domain modulated frequency-dependent linear and second non-linear response properties of doped quantum dots,” Physica Status Solidi B, vol. 248, no. 8, pp. 1941–1948, 2011. View at Publisher Β· View at Google Scholar
  24. V. Halonen, P. HyvΓΆnen, P. PietilΓ€inen, and T. Chakraborty, β€œEffects of scattering centers on the energy spectrum of a quantum dot,” Physical Review B, vol. 53, no. 11, pp. 6971–6974, 1996. View at Google Scholar Β· View at Scopus
  25. V. Halonen, P. Pietilinen, and T. Chakraborty, β€œOptical-absorption spectra of quantum dots and rings with a repulsive scattering centre,” Europhysics Letters, vol. 33, no. 5, pp. 377–382, 1996. View at Publisher Β· View at Google Scholar
  26. R. Turton, The Quantum Dot. A Journey into Future Microelectronics, Oxford University Press, New York, NY, USA, 1995.
  27. L. Jacak, P. Hawrylak, and A. Wojos, Quantum Dots, Springer-Verlag, Berlin, Germany, 1998.
  28. M. A. Reed, Quantum Dots, Scientific American, Ocala, Fla, USA, 1993.
  29. Y. Alhassid, β€œThe statistical theory of quantum dots,” Reviews of Modern Physics, vol. 72, no. 4, pp. 895–968, 2000. View at Google Scholar Β· View at Scopus
  30. J. Adamowski, A. KwaΕ›niowski, and B. Szafran, β€œLO-phonon-induced screening of electron-electron interaction in D-centres and quantum dots,” Journal of Physics Condensed Matter, vol. 17, no. 28, pp. 4489–4500, 2005. View at Publisher Β· View at Google Scholar Β· View at Scopus
  31. S. Bednarek, B. Szafran, K. Lis, and J. Adamowski, β€œModeling of electronic properties of electrostatic quantum dots,” Physical Review B, vol. 68, no. 15, Article ID 155333, 9 pages, 2003. View at Google Scholar
  32. P. D. Siverns, S. Malik, G. McPherson et al., β€œScanning transmission-electron microscopy study of InAs/GaAs quantum dots,” Physical Review B, vol. 58, no. 16, pp. R10127–R10130, 1998. View at Google Scholar Β· View at Scopus
  33. S. Bednarek, B. Szafran, and J. Adamowski, β€œTheoretical description of electronic properties of vertical gated quantum dots,” Physical Review B, vol. 64, no. 19, Article ID 195303, pp. 1953031–19530313, 2001. View at Google Scholar Β· View at Scopus
  34. A. P. Alivisatos, β€œPerspectives on the physical chemistry of semiconductor nanocrystals,” Journal of Physical Chemistry, vol. 100, no. 31, pp. 13226–13239, 1996. View at Google Scholar Β· View at Scopus
  35. A. A. Guzelian, U. Banin, A. V. Kadavanich, X. Peng, and A. P. Alivisatos, β€œColloidal chemical synthesis and characterization of InAs nanocrystal quantum dots,” Applied Physics Letters, vol. 69, no. 10, pp. 1432–1434, 1996. View at Google Scholar Β· View at Scopus
  36. I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, Corrected and Enlarged Edition, Academic Press, City, State, USA, 1980.
  37. J. J. P. Stewart and K. M. Dieter, β€œCalculation of the nonlinear optical properties of molecules,” Journal of Computational Chemistry, vol. 11, no. 1, pp. 82–87, 1990. View at Google Scholar
  38. P. Dutta and S. P. Bhattacharyya, β€œOn exact calculation of response properties of oscillators in static electric field: a Fourier grid Hamiltonian approach. I. One-dimensional systems,” International Journal of Quantum Chemistry, vol. 51, no. 5, pp. 293–305, 1994. View at Publisher Β· View at Google Scholar
  39. G. Maroulis and A. J. Thakkar, β€œMultipole moments, polarizabilities, and hyperpolarizabilities for N2 from fourth-order many-body perturbation theory calculations,” The Journal of Chemical Physics, vol. 88, no. 12, pp. 7623–7632, 1988. View at Google Scholar Β· View at Scopus
  40. G. Maroulis, β€œHyperpolarizability of H2O,” The Journal of Chemical Physics, vol. 94, no. 2, pp. 1182–1190, 1991. View at Google Scholar Β· View at Scopus
  41. G. Maroulis, β€œHyperpolarizability of H2O,” Journal of Chemical Physics, vol. 94, pp. 1182–1190, 1991. View at Google Scholar
  42. T. Chakraborty, Quantum Dots-A Survey of the Properties of Artificial Atoms, Elsevier, Amsterdam, The Netherlands, 1999.
  43. S. V. Nistor, M. Stefan, L. C. Nistor, E. Goovaerts, and G. Van Tendeloo, β€œIncorporation and localization of substitutional Mn2+ ions in cubic ZnS quantum dots,” Physical Review B, vol. 81, no. 3, Article ID 035336, 2010. View at Publisher Β· View at Google Scholar Β· View at Scopus
  44. P. A. Sundqvist, V. Narayan, S. StafstrΓΆm, and M. Willander, β€œSelf-consistent drift-diffusion model of nanoscale impurity profiles in semiconductor layers, quantum wires, and quantum dots,” Physical Review B, vol. 67, no. 16, Article ID 165330, pp. 1653301–1653308, 2003. View at Google Scholar Β· View at Scopus