Christos S.
Garoufalis
*a,
David B.
Hayrapetyan
bc,
Hayk A.
Sarkisyan
b,
Paytsar A.
Mantashyan
bc,
Zaiping
Zeng
d,
Iosif
Galanakis
a,
Gabriel
Bester
e,
Torben
Steenbock
e and
Sotirios
Baskoutas
*ae
aMaterials Science Department, University of Patras, 26504 Patras, Greece. E-mail: garoufal@upatras.gr; bask@upatras.gr
bDepartment of General Physics and Quantum Nanostructures, Russian-Armenian University, 123 Hovsep Emin Str., Yerevan 0051, Armenia
cInstitute of Chemical Physics after A.B. Nalbandyan of NAS RA, 5/2 Paruyr Sevak St., Yerevan 0014, Armenia
dKey Laboratory for Special Functional Materials of Ministry of Education, Collaborative Innovation Center of Nano Functional Materials and Applications, and School of Materials Science and Engineering, Henan University, Kaifeng, Henan 475001, China
eInstitut für Physikalische Chemie, Universität Hamburg, Grindelallee 117, 20146 Hamburg, Germany
First published on 2nd April 2024
Quantum dots are widely recognized for their advantageous light-emitting properties. Their excitonic fine structure along with the high quantum yields offers a wide range of possibilities for technological applications. However, especially for the case of colloidal QDs, there are still characteristics and properties which are not adequately controlled and downgrade their performance for applications which go far beyond the simple light emission. Such a challenging task is the ability to manipulate the energetic ordering of exciton and biexciton emission and subsequently control phenomena such as Auger recombination, optical gain and photon entanglement. In the present work we attempt to engineer this ordering for the case of InP QDs embedded in polymer matrix, by means of their size, the dielectric confinement and external electric fields. We employ well tested, state of the art theoretical methods, in order to explore the conditions under which the exciton–biexciton configuration creates the desired conditions either for optical gain or photon entanglement. Indeed, this appears to be feasible for QDs with small diameters (1 nm, 1.5 nm) embedded in a host material with high dielectric constant and additional external electric fields. These findings offer a new design principle which might be complementary to the well-established type II core–shell QDs approach for achieving electron–hole separation.
Various techniques have been employed in order to decrease the Auger rates. The most promising efforts include either the utilization of elongated quantum dots (QDs)6 of core/shell hetero-QDs,7 and the application of external electric fields.8 These approaches favor the reduction of exciton–exciton coupling while in the same time they retain the advantages of strong confinement. The most innovative approach to address the issue of Auger decay involves the introduction of techniques and/or configurations which allow for optical-gain in the single-exciton regime, wherein Auger recombination is rendered inactive. A successful implementation of this idea has been achieved by Klimov et al.,9 by taking advantage of the band alignment in type-II core/shell hetero-quantum dots. Another successful approach to suppress Auger recombination is to use large QDs where the wave function overlap between electron and holes is reduced and non-local effects appear.10
In this particular case, the electrons and holes are confined in different parts of the nanostructure. The electrons remain at the core while the holes are confined in the shell. The resulting spatial separation between the two charge distributions leads to the emergence of a strong local field, which affects the lowest energy transitions by means of a Stark shift. If EX and EXX are the exciton (X) and biexciton (XX) total energies and ħωXX and ħωX denote the energy of the XX → X and X → 0 transitions, then the energy difference between these peaks, is represented by the symbol ΔXX = ħωXX − ħωX = Jee + Jhh − 2Jeh. The sign of this quantity plays an important role in the optical gain of nanocrystals since it indicates the relative positions between exciton and biexciton peaks. A negative value of ΔXX implies that the X–X interaction is attractive in nature forcing the XX → X transition to shift towards a lower energy. Such an effect can potentially increase the absorption probability of a photon emitted from transition (i.e. X → 0) that lies energetically just above the biexciton absorption peak (i.e. X → XX), having a negative impact on the possible existence of any lasing properties. On the other hand, a positive value of ΔXX would indicate a X–X repulsion which, if comparable to or greater than the transition line width, would enhance the lasing properties since the photon emitted by the X → 0 transition has lower energy than the one required for the X → XX transition.9 Positive values of ΔXX have been observed for small strained self-organized InAs/GaAs QDs.11
In parallel, it is widely recognized that the creation of a photon source that produces polarization entangled photons on demand is a vital prerequisite for applications in the fields of quantum information and quantum optics.12,13 To this end, QDs offer the possibility to take advantage of the biexciton–exciton transitions14,15 and create such entangled photon pairs. However, the favorable conditions for entanglement are destroyed by the fine structure splitting (FSS) which is present in the intermediate exciton state and significantly affects the radiative cascade. Various approaches have been successfully employed to mitigate the FSS, including the application of suitable electric,16 magnetic17 or strain18 as well as a combination of electric and strain fields,19 engineering of QD size and composition20,21 and controlling the growth process in order to produce highly symmetric QDs.22
An alternative way which could help to overcome the FSS and lead to the generation of entangled photon pairs, involves the fine tuning of the biexciton binding energy to zero.23,24 This has been experimentally achieved with techniques which include the introduction of lateral strain in the material25 or suitable electric fields.24,26
Motivated by the aforementioned open issues in the field and our recent findings regarding the ways that the dielectric confinement affects the exciton properties in wurtzite colloidal ZnO quantum dots,27 we chose to investigate possible routes to address the challenges mentioned above. In particular, it has been previously found27 that for the case of ZnO QDs embedded in a host material the exciton energy increases and the exciton binding energy decreases as the dielectric constant of the surrounding material becomes larger. This behavior is induced by the presence of a self-polarization potential due to the dielectric mismatch on the QD's surface. Specifically, this potential primarily affects the localization of the hole charge density, causing it to migrate towards the center of the QD, while in the same time the electron density appears to be almost unaffected. This different behavior leads to a subsequent reduction of the overlap between the electron and hole and inevitably reduces the exciton binding energy. The self canceling arrangement of positive and negative densities is destroyed by this charge separations and the emerging internal local field modifies the conditions for the biexciton creation. As a result, the exciton and biexciton properties of QDs might be manipulated by suitably engineering their host environment and their dielectric mismatch. This provides us with an extra degree of freedom which, when complemented with the additional application of external fields, can greatly enhance our capabilities to fine tune the exciton and biexciton ordering as desired.
Following this line of work, we turn our attention to zinc blende InP quantum dots and investigate the impact of the surrounding host material and the externally applied electric field on the exciton and biexciton fine structure. Actually InP QDs/polymer nanocomposites have been widely used as a means to improve the stability of InP QDs against heat and moisture28 and improve their quantum yield.29 To this purpose we employ large-scale atomistic pseudopotential calculations which incorporate the effects of multiband coupling, multivalley coupling, spin–orbit interaction, and excitonic correlation via the configuration interaction method.27,30–33 The primary objective is to explore the possible existence of favorable conditions (QD size, host material and external electric fields) which may actually allow us to control the exciton–biexciton ordering. The accuracy of the observed trends through the described computational approach is also tested against a small set of benchmark calculations performed with Density Fuctional Theory (DFT) methods coupled to our own CI package.
![]() | (1) |
In this formula, the QD potential is represented by a superposition of screened atomic pseudopotentials centered at atom sites
and vSOα represents the spin–orbit term. E is the applied electric field and |q| is the electron charge in atomic units (i.e. |q| = 1). In general the atomic pseudopotentials used in the EPM are thoroughly optimized in order to accurately reproduce the band structure of the bulk material, including the values of the hole and electron effective masses. In the particular case the adopted potentials for InP have been taken from literature.34 Furthermore, in order to consider the dielectric confinement and the external electric field, the model Hamiltonian has been complemented by the terms Vs (self-polarization) and
(electric field) respectively. The analytic form of the self-energy term can be found in our previous work.27
Once the eigenfunctions of the system have been calculated, a screened configuration interaction (CI) scheme is employed in order to obtain the excitations of the QD. The excitonic wavefunctions are expanded in terms of singly excited Slater determinants constructed from the single-particle wavefunctions produced by the solution of the corresponding Schrödinger equation. Since the CI expansion is limited only to single excitations and the active space includes only a small number of electron and hole states, it is necessary to screen all the Coulomb interactions.30 As for the exchange interactions, it has been argued35 that in the limiting case where transitions are considered over a wide energy range, they should be unscreened. However, in our case where only a few lowest energy transitions are considered, the interaction should be screened by the full dielectric function. Thus, in the present work we will employ the previously used approximation,27 in which the inverse dielectric constant consists of an electronic high-frequency contribution36 and an ionic low-frequency contribution.37,38
In order to find the optimal size and composition of the active space, an initial set of benchmark calculations was performed. The results indicated that the required accuracy as well as the desired computational efficiency can be achieved by including in our CI treatment the 10 highest occupied states and 8 lowest unoccupied states.
For the case of biexciton calculations the CI expansion includes only doubly excited determinants (i.e. two electrons and two holes). The optical transitions are calculated within the dipole approximation.
Although the adopted computational scheme has been proved extremely accurate for a wide range of materials and properties throughout the years, it has been considered necessary to further validate its trends by means of selective first principles calculations. For this reason we have performed some additional benchmark DFT calculations for the smallest QDs in order to verify that the ordering change in the position of exciton and biexciton peaks is validated by a completely different methodology. In particular, all structures considered in this step were optimized at a PBE/def2-SVP39,40 level of theory. In the relaxed geometries, the solvation effects (with the desired dielectric constants) were taken into consideration using the COnductor-like Screening MOdel (COSMO) model.41 On top of this, we carried out two-component, single-point calculations42,43 using the x2c-SVPall-2c basis set44 and the PBE0 hybrid functional.45–49 The excited states were then calculated using the same screened CI methodology50 used for the EPM calculations.27
The exciton energy obtained through the CI method for various quantum dot diameters across four different dielectric matrices is displayed in Fig. 1a. The anticipated trend of the exciton energy decreasing as a function of the quantum dot (QD) diameter is consistently observed in all scenarios, irrespective of the surrounding environment. Furthermore, it is evident that the dielectric constant of the host material causes a shift in the exciton energy towards higher values across all size ranges. The observed trend can be safely attributed to the self-polarization potential, as has been found both experimentally52–54 and theoretically55,56 in the case of InP quantum dots. The variation of the exciton binding energy across the various quantum dot (QD) diameters within the four distinct dielectric matrices is illustrated in Fig. 1b. The exciton binding energy exhibits a decreasing behavior as the quantum dot diameter increases across all host matrices as expected from quantum confinement. However, the most interesting feature so far, is related to the dependence of the binding energy of the excitons on the dielectric constant of the surrounding material. As the ε of the host increases the binding energy decreases significantly (initially) but then it seems as if it saturates. For example as we move from small values (e.g. M4) to larger ones (e.g. M15) the change is quite prominent but as we keep increasing ε (e.g. from M15 to M40 or M80) the shift of the curve seem to converge to a specific limit.
![]() | ||
Fig. 1 (a) Exciton energy as a function of QD diameter for all the host materials considered. (b) Dependence of the exciton binding energy on size and surrounding material. |
In order to gain some insight into this behavior, we constructed density difference isosurfaces for the holes and electrons participating in the formation of excitons and focus on the charge redistribution induced by the presence of the matrix material, as described in ref. 26. The density differences defined by the formula ρdiff = ρ(εin≠εout) − ρ(εin=εout) can give us a visual representation of any rearrangement of the hole or electron charge distributions. In particular, areas of the ρdiff isosurface with positive values correspond to charge accumulation areas while the negative ones to charge depletion (see Fig. 2). Interestingly it is found that only the hole density is significantly affected by the presence of the host material (with increased εout) resulting in a charge migration towards the center of the QD. On the other hand, the electron density exhibits migration towards the surface but with a very small percentage. After all, it would seem that the dielectric mismatch at the dot's interface with the host material induces some mild charge separation within the dot, which in turn reduces the Coulomb interaction and accordingly the exciton binding energy. Moreover, this effect appears to quickly reach saturation, and from a certain point onward, an increase in the value of ε does not result in any additional charge separation or decrease in the exciton's binding energy.
The biexciton energy behaves similarly to the exciton energy and decreased with increasing QD size and with decreasing εout (see Fig. 3a). In the same time, ΔXX appears to increase as we move to hosts with larger ε and interestingly, when the value of ε is adequately large, it might even change sign (see Fig. 3b). This is actually the case for the smallest diameters (i.e. D = 1.0 nm and D = 1.5 nm). Actually the same behavior has also been observed for small strained self-organized InAs/GaAs QDs.11 This is rather intriguing for two reasons: the first one is that as it has already been stated, positive values of ΔXX imply optical gain while zero values of ΔXX fulfill the condition for photon entanglement.24 It is noteworthy that in Fig. 3b the M4 plot exhibits a non-monotonic behavior. This can be attributed to the relative values of εin and εout. In particular, the value of εin varies from approximately 2.9 to almost 5 (for the specific size ranges). As a result, the difference εin–εout for the case of M4 changes sign as the size of the QDs increases. Since this difference is important for the self-polarization effects27 induced by the matrix environment, we deduce that it is tightly connected with the observed behavior. In all other matrices the εout is significantly larger than εin and as a result this trend does not appear.
![]() | ||
Fig. 3 (a) Variation of ħωXX as a function of size. (b) ΔXX = ħωXX − ħωX as a function of the QD diameter in different host materials. |
This behavior can be understood if one considers explicitly the terms which contribute to the ΔXX: ΔXX = Jee + Jhh − 2Jeh. The values of all right hand terms depend on the localization patterns of the hole and electron charge densities. Fig. 4 makes this point very clear: when there is a dielectric mismatch, the hole density moves towards the center of the QD, but the electron density is essentially insensitive (moves very slightly towards the surface). As a result, the effect of the host material on the charge distribution of the carriers indicates that although Jee is only marginally affected, the term Jhh increases as expected (since the holes tend to localize to the center) and the term Jeh decreases since the average e–h distance increases. As is also clearly seen from Fig. 3b, above a specific threshold diameter (e.g. D = 3 nm for M4 and M15, D = 4 nm for M40, M80) the biexciton binding energy decreases with decreasing confinement. Actually the radius which defines the onset of this threshold diameter moves to higher values as the dielectric constant of the host materials increases.
![]() | ||
Fig. 4 Biexciton charge polarization due to the presence of dielectric mismatch. The color coding and the description is similar to the one of Fig. 2. In the specific case the percentage of hole charge movement towards the center is ∼13% while the corresponding value for the electron (towards the surface) is 1%. |
This means that, as it could be expected, for smaller dots the repulsive term Vhh is significantly enhanced due to the combined effects of strong spatial confinement and hole migration towards the center of the QD, hence affecting the binding characteristics of the biexciton. As the size of the QD increases this combined effect becomes looser and the attractive term Veh takes the leading role.
Due to the fact that the observed changes in the ordering of exciton and biexciton peaks (i.e. ΔXX > 0 or ΔXX < 0) is a rather important finding, we chose to verify the trend by means of a completely different computational approach. In particular, we performed some benchmark DFT calculations, as described in the previous section, for the smallest nanoparticles with different dielectric constants using the COSMO solvation model. Indeed it was verified that the ordering change of exciton and biexciton peaks takes place for small QDs and large values of ε (e.g. for QD diameter equal to 1 nm, ΔXX = 18.8 meV for M40 and ΔXX = 19.7 meV for M80).
As a means to further increase our ability to manipulate the relative position of the X and XX peaks, the application of external electric field was also considered. This was achieved by incorporating into the Hamiltonian of eqn (1), the last term. The electric field was applied along the [001] direction, and set to 100, 200, 300 and 400 kV cm−1. Its impact on the exciton can be expressed by means of the induced Stark shift (see Fig. 5a which depicts its dependence on the QD diameter and the strength of the field). As is clearly shown the absolute value of the shift is an increasing function of both the QD diameter and the field strength. Beyond that, our results revealed that, indeed, the application of an external electric field offers some extra ability to manipulate the relative peak positions. This can be clearly seen in Fig. 5b, where the field seems to shift the whole ΔXX curve towards larger values. In this way even the initially bound biexciton of the 2.0 nm dot appears to move close to the zero binding energy threshold. Moreover, for the smaller dot, the field induced enhancement of ΔXX (for the specific host) appears to become as large as 14 meV, which is quite close to the experimental line widths of fluorescence line narrowing (FLN) spectra (15–30 meV) of InP quantum dots.57
This journal is © The Royal Society of Chemistry 2024 |