Exciton-photon interactions in semiconductor nanocrystals: {radiative transitions, non-radiative processes,} and environment effects

In this review we discuss several fundamental processes taking place in semiconductor nanocrystals (quantum dots, QDs) when their electron subsystem interacts with electromagnetic (EM) radiation. The physical phenomena of light emission and EM energy transfer from a QD exciton to other electronic systems such as neighbouring nanocrystals and polarisable 3D (semi-infinite dielectric or metal) and 2D (graphene) materials are considered. The cases of direct (II-VI) and indirect (silicon) band gap semiconductors are compared. We also cover the relevant non-radiative mechanisms such as the Auger process, electron capture on dangling bonds and interaction with phonons. The emphasis is on explaining the underlying physics and illustrating it with calculated and experimental results in a comprehensive, tutorial manner.

In this review we discuss several fundamental processes taking place in semiconductor nanocrystals (quantum dots, QDs) when their electron subsystem interacts with electromagnetic (EM) radiation. The physical phenomena of light emission and EM energy transfer from a QD exciton to other electronic systems such as neighbouring nanocrystals and polarisable 3D (semi-infinite dielectric or metal) and 2D (graphene) materials are considered. The cases of direct (II-VI) and indirect (silicon) band gap semiconductors are compared. We also cover the relevant non-radiative mechanisms such as the Auger process, electron capture on dangling bonds and interaction with phonons. The emphasis is on explaining the underlying physics and illustrating it with calculated and experimental results in a comprehensive, tutorial manner.

I. INTRODUCTION
Semiconductor nanostructures form a basis for modern electronic technologies. Nowadays, they are employed in a wide range of applications in the fields of optoelectronics, photonics, photovoltaics, biosensing, photocatalysis, etc. By virtue of the quantum confinement effect, their electronic spectra and, consequently, the related optical properties depend on the nanostructure size. This feature is most pronounced for zero-dimensional (0D) objects, nanocrystals (NCs), where the conduction band electron motion is fully localized in all directions. Consequently, at least a part of the NC energy spectrum is completely discrete. In the limiting case of strong quantum confinement effect, when the NC size is much less than the effective exciton Bohr radius, the former strongly influences the electron and hole energies. It leads to the sizedependent energies of photons emitted or absorbed by the nanocrystals and allows one to control their optical spectra. For this reason, they sometimes are called "artificial atoms" 1 . A more scientific term to distinguish semiconductor NCs with strongly size-dependent electronic and optical properties is Quantum Dot (QD), which will be used in this article. 238 After the pioneering works by Efros and Efros 2 and Brus 3 explaining the origin of size-dependent optical spectra of nanocrystals, the interest to these objects grew exponentially through the 1980s and 1990s. A variety of methods has been used for preparation of crystallites with sizes not greater than several nanometers, such as ion implantation 4-6 , chemical vapour deposition from a gas phase, 7 magnetron sputtering 8 , colloidal synthesis 9 , electron beam epitaxy 10 , etc. One of these methods, chemical growth in colloidal solutions, allows for obtaining good quality NCs of a broad range of semiconductors with an almost spherical shape and a rather narrow size distribution, characterized by rather narrow emis-sion bands, which is, as a rule, extremely desirable for applications. This method was first suggested for II-VI NCs by Murray et al. 11 and became widely employed later (see, e.g., chapters by Kudera et al., Reiss, Gaponik and Rogach in the book in 12 for review). The colloidal chemistry methods work fairly well for II-VI and III-V materials such as CdSe, CdS, CdTe, InAs, InP, etc. The authors of 12,13 and several other groups have been able to produce good QDs of some IV-VI materials 14,15 and silicon QDs (see, e.g., reviews in 16,17 ).
Colloidal QDs of II-VI materials (especially CdSe and CdTe) are probably the most studied. During the last decade, a considerable research activity has been focused on the synthesis and investigation of the optical properties of ZnSe QDs. These nanocrystals exhibit large blue shift of the photoluminescence 18 and high enough (up to 50 percent) quantum yield 19 . Similarly to CdSe NCs, the ZnSe QDs can be synthesized with a sufficiently narrow size distribution, which is possible to control by temperature 20 . Due to ZnSe's lower toxicity, compared to many other II-VI semiconductors, these QDs turn out to be attractive objects also for biosensing 21 .
Further improvement of the light emission properties of such QDs was achieved by fabrication of core-shell structures, successful for several pairs of II-VI 22 and IV-VI 23 materials, where the shell made of a wider band gap material provides a better protection of the quantized electronic states in the QD core 12 . For ZnSe/ZnS core-shell QDs it was revealed, in particular, that thermodynamical (slow) growth of the ZnS-shell on the colloidal ZnSe-core leads to the quantum yield increase because of decreasing amount of the traps at the core-shell interface 24 .
The maturity of this technology of synthesis of semiconductor nanocrystals 239 is witnessed by the incorporation of colloidal QDs into real-world products such as colour displays 29 and light-emitting diodes for lighting 30 , both already commercialized by large companies such as Samsung 31 and OSRAM 32 , respectively (see Ref. 33 for a recent review covering both colloidal and epitaxial QDs). More on the scientific research side, a broad variety of nanostructures can be prepared using colloidal NCs as building blocks, in particular, multilayer structures of QDs of different average size, deposited on different substrates 12 . Combining these structures with other materials, such as organic dielectrics 34 , epitaxial quantum well (QW) heterostructures 35 , metallic nanoparticles 36,37 , patterned metallic surface 38 or graphene 39,40 can result in new interesting effects and applications, which physics is related to the coupling of QD excitons with elementary excitations (such as surface plasmons or QW excitons) in the surrounding materials. One interesting possibility is the excitation of colloidal QDs by pumping energy through a nearby epitaxial QW, as demonstrated in 35 . Another research topic with QDs, popular in the past two decades is the modification of their emission properties caused by the strong light-matter interaction, which has been achieved for a variety of composite structures that include point emitters embedded in Fabry-Perot microcavities, micropillars and photonic crystals as recently reviewed in 41 .
In this article, we overview several aspects of the exciton-photon interactions in nanocrystal QDs with focus on the role of non-radiative processes and environment, which affect the light emission. The QD photoluminescnce is a result of competition between the radiative and various non-radiative processes, such as capture on dangling bonds, multi-phonon intra-band relaxation, Auger recombination, etc. In the past decade, a growing interest has been paid to the multi-exciton dynamics in nanocrystals because of their high potential for photovoltaic applications. Such processes as carrier multiplication (or multi-exciton generation) in nanocrystals, as well as the Auger recombination (which is just a fast reverse process), are widely discussed in the literature 23,43-45 and we will consider these effects below. Special attention is traditionally paid to nanostructured silicon 25,46-49 because of its widest use in microelectronics, high purity, natural abundance, low cost, and non-toxicity; we shall also dedicate some space to the exciton-photon interactions and main non-radiative processes mentioned above, in Si NCs and related nanostructures.
As a rule, one deals with nanocrystal ensembles (rather than with isolated QDs). Therefore, non-radiative energy exchange between the nanocrystals is possible via the Förster-type 50 and Dexter-type 51 exciton migration. Both mechanisms were originally proposed for fluorophores. The second one, based on the exchange interaction between electrons located on different sites can be neglected in the presence of allowed dipole transitions 51 . Recently, a third, so-called exciton tandem tunneling mechanism was proposed 52 , which is specific of connected nanocrystals. The most universal process named Förster Resonant Energy Transfer (FRET) was first observed for QDs by Kagan et al. in specially designed films containing two different sizes of nanocrystals acting as donors and acceptors, respectively 53 . Later it was shown in a number of works 54-59 performed on systems composed of two different QD species that the luminescence of the smaller dots (donors) is quenched by the large dots (acceptors), whose emission in turn is enhanced. These studies demonstrated a dependence of the FRET effect on the NC density and spatial arrangement. Results of other experiments performed on multilayer SiO x /SiO 2 structures 60-62 , porous Si 63 , three-dimensional (3D) ensembles of Si 64,65 crystallites were interpreted as a manifestation of exciton migration via FRET-type mechanisms. The effect has a potential interest for photonics (e.g., photoluminescence upconversion in QD ensembles 66 ), sensing 67 , lighting and energy harvesting (e.g., by unidirectional energy transfer in size-gradient layered QD assemblies 54,55 or fractal aggregates 68,69 ). We shall discuss it below, in particular, addressing the question whether the FRET rates and length scales can be tuned by the photonic environment 70 .
The influence of the environment on the excitonrelated optical properties of NCs is the third topic that will be discussed in this overview. The most common effect is the energy transfer between a photoexcited QD and surrounding materials, which can be reversible or not. It has the same physical nature as FRET. In most cases, it is responsible for the photoluminescence (PL) quenching 71,72 . However, the environment can be used for engineering the photonic density of states (DOS) in the vicinity of the QDs, often referred to as the Purcell effect 73 . Observed for the first time in QDs about twenty years ago 74 , it has recently been discussed with respect to radiative decay rates of embedded point emitters 75 . Moreover, strong near-field effects associated with localized surface plasmon resonance (LSPR) in metallic (nano-)structures can enhance PL emission 38,76 . Moreover, one can think of energy transfer from a recombining QD exciton to propagating surface plasmons that would carry the energy over a large distance and then eventually transfer it to another QD (by creating an exciton), thus allowing for a long-range exciton transport between two dots, much more efficient than if it occurred directly. Below we consider the influence of a flat interface between two media on the PL emission and FRET rates for a QD emitter located in the vicinity of the interface 240 .
The article is organized as follows. In Section II, we introduce basic notions of radiative transitions and discuss their rates (i.e., probabilities per unit time) for NCs of materials possessing direct or indirect band structure. In particular, the rates are calculated for intrinsic and doped silicon NCs, which are either hydrogen-coated or halogenated. Section III is devoted to non-radiative transitions and considers the Auger recombination, danglingbond traps and phonon-assisted relaxation of hot carriers. In Section IV, multiple exciton generation initiated by a highly excited electron or hole is considered. Section V is dedicated to the exciton transfer processes between two QDs via the FRET mechanism. In Section VI, emission decay and FRET rates near a plane inter-face between two dielectrics or a dielectric and a metal are discussed. The last section is left for summary and conclusions.

A. Spontaneous Emission Rate
The emission rate of a point dipole emitter located in an infinite dielectric medium with the permittivity ε 1 is given by 79 where k 1 = √ ε 1 ω/c is the modulus of the wavevector, ω is the oscillation frequency, c is the velocity of light in vacuum,h denotes the Planck constant and d is the transition dipole moment matrix element. This expression, with some adaptation, describes the radiative decay rate of an exciton in a QD made of a direct band gap semiconductor. We shall write it as the inverse of the radiative lifetime: Here, m 0 and −e are the free electron mass and charge, respectively, and E if and p if denote the transition energy between the initial and final states and the momentum matrix element, respectively. The latter is more natural to use in the electronic band theory of crystals; it can be related to d appearing in Equation (1) as p = −im 0 ωd/e. The parameter κ in Equation (2) arises due to the difference in the dielectric constants of the crystallite (ε) and its surrounding (ε 1 ) and takes account of the so called depolarization effects 80 . Hereafter, we consider the crystallite as a sphere with the radius R and dielectric constant ε embedded in a homogeneous medium with the dielectric constant ε 1 . In this case, κ is given as 81 This parameter varies within a wide range of values depending on the dielectric constants of the two materials.

B. Nanocrystals of Direct-Band-Gap Materials
The majority of III-V and II-VI semiconductors are direct band gap materials. In this case, in the framework of the envelope function approximation, where p cv is the standard bulk-like momentum matrix element of inter-band transitions, and F c,i and F v,f , in case of the photon emission, stand for the envelope functions of the initial state in the conduction band and the final state in the valence band, respectively. The scalar product of these envelope functions determines selection rules for such transitions.
In the simplest case, the QD is treated as an infinitely deep spherical potential well with radius R and the conduction and valence bands are assumed simple and parabolic. Consequently, the electron and hole ground states are described by the same envelope function and the scalar product in (4) equals to unity. It means that the radiative recombination rate turns out to be almost independent of the QD radius. This dependence manifests only in the transition energy, E if . In the limiting case of a strong quantum confinement regime, where R a ex (a ex denotes the effective exciton Bohr radius), the exciton binding energy can be neglected. Then, the model of infinitely deep spherical potential well yields E if (R) = E ∞ g +h 2 π 2 /2µR 2 , where E ∞ g is the band gap energy of the bulk semiconductor and µ is the reduced effective mass of the electron-hole pair.
The above assumption concerning the extreme strong confinement regime is not quite realistic; usually the QD radius can be just slightly smaller than a ex (for instance, a ex = 5.6 nm for CdSe 12 ) and the exciton binding energy can be estimated as 3 E b ≈ 1.79e 2 /(εR). Furthermore, the most common semiconductors have a degenerate valence band with light and heavy hole sub-bands, so that the confined holes are described by two envelope functions, which are completely different from the electron's one 82 . Thus, the matrix element (4) does depend on R (see, e.g., chapter by Vasilevskiy in the book 12 ). Still, the dependence of τ −1 R on the QD size is rather weak for direct-gap materials as compared to, e.g., silicon NCs considered below. Using typical parameters for CdSe QDs, ε ≈ 10, ε 1 ≈ 2 and p 2 cv /m 0 ≈ 10 eV 82 , Thus, τ R is of the order of a nanosecond, in accordance with experimental observations.
It is well known that introduction of shallow impurities is capable of modifying electronic properties of bulk silicon. Similarly, doping with shallow impurities influences the electronic structure of silicon NCs [110][111][112][113][114][115][116][117][118] , which, in turn, affects the electron-hole radiative recombination. It was revealed that doping of Si nanocrystals with P or Li is (under certain conditions) capable of improving their emittance 85,88,90,119 . In particular, the performed calculations 88,100,107,108 show that doping with P or Li can essentially increase the radiative transition rates.
For NCs with R > ∼ 1 nm, semi-empirical methods are usually employed for computing their electronic structure, such as envelope function approximation 120-122 or tight binding model 123 . The performed calculations of the phonon-assisted radiative recombination rates in Si crystallites in the SiO 2 matrix yielded the values varying from ∼ 10 2 to ∼ 10 5 s −1 as the crystallite radius decreases from 3 to 1 nm. Within a simple model of an infinitely deep spherical potential well, the dependence of the radiative decay rate on the crystallite radius is [122]. For the no-phonon radiative transitions, the rates turn out to be much lower, by one to three orders of magnitude with increasing the size within the same range. In this case, τ −1 R ∝ R −8 [122], i.e., the rate sharply drops as R increases. At the same time, doping with phosphorus provides a several orders of magnitude increase of the radiative recombination rates for the no-phonon transitions. The rates become even 1-3 orders greater than those of the phonon-assisted transitions. They slowly decrease with increasing the NC size, especially for high concentration of phosphorus 100 .
As calculations show, the radiative transitions including both no-phonon and phonon-assisted ones become faster in small Si crystallites. However, for crystallites with R < ∼ 1 nm, all semi-empirical methods of calculation turn out to be essentially less accurate. In this case, first-principle calculations have to be carried out. In Figure 1, the radiative recombination rates computed at room temperature T within the Casida's version 124  Here, NCs are considered in vacuum, therefore, ε 1 = 1, while for Si ε = 12. In this case, the rate values become ∼16 times lower compared to the case of Si NCs embedded in silicon dioxide due to the factor κ √ ε 1 .
As seen in Figure 1, generally, doped crystallites demonstrate 1-2 orders of magnitude faster transitions than intrinsic ones. For Si crystallites doped with phosphorus, increasing rates of the radiative transitions have been explained by the efficient mixing of the electronic states in the Γ and X points of the Brillouin zone (Γ-X mixing) caused by the short-range field of the phosphorus ion 88,100 . Meanwhile, the radiative decay in the Lidoped crystallites becomes faster due to the high density of states above the inter-band energy gap (in the rangē Surface chemistry is another method capable of modifying electronic structure of nanocrystals 125,126 . This method is especially efficient for small crystallites, where the role of surface rises. It has been demonstrated experimentally 105,127,128 that chemical synthesis allows to fabricate Si nanocrystals with various ligand coatings, which emit light in visible range with typical nanosecond radiative lifetimes 105,129,130 . Calculations based on the tight-binding model 106 performed for CH 3 -capped Si nanocrystals confirmed the conclusion on the increase of the radiative recombination rate due to the formation of direct-like electronic structure of the crystallite. Let us now consider some other kind of the surface reconstruction of the crystallites, the halogenation. Table  1 presents calculated rates of radiative transitions between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in completely halogenated with Cl-, Br-and H-coated silicon crystallites. Table 1. Calculated rates (s −1 ) of the radiative HOMO-LUMO transitions in H-, Cl-, and Br-passivated silicon nanocrystals.
It is possible to see from Table 1 that the surface halogenation significantly slows down the radiative recombination in silicon nanocrystals, especially in the case of the bromine coating. Such a slowing down is due to the spatial separation of the electron densities of the HOMO and LUMO states, which takes place in halogenated Si crystallites [131][132][133] . The separation turns out to be more pronounced in the Br-passivated crystallites, in which the HOMO states are strongly concentrated near the Br atoms. Presumably, the bromine atoms produce strong lattice distortions, and the stress fields have an additional localizing effect on the electron density. An effect of the rate reduction was also found for completely fluorine-passivated Si nanocrystals 103 .

III. NONRADIATIVE PROCESSES
As is well known, the emissivity of any light-emitting entity depends not only on the radiative transitions' rate, but also on the speed of nonradiative decay processes. In the case of NCs, these include Auger recombination, capture of photo-electrons on dangling bonds and cooling of hot carriers. Below, we briefly discuss all these processes.

A. Auger Recombination
In order for the Auger process to occur, at least one negative or positive trion (exciton + one electron or exciton + one hole) has to be initially excited in a nanocrystal, as shown in Figure 2. The trion turns into a highenergy single-particle excitation as indicated by the arrows in the figure.
The rate of an Auger-type transition can be written as where U if is the matrix element of the two-particle Coulomb interaction operator, the two-electron wave functions are approximated by products of the singleelectron ones and the δ-function expressing the energy conservation in the Fermi's Golden Rule has been broadened with the half-width Γ = 10 meV. All the possible final states of the Auger electron with energy E f (shown in Figure 2) are summed up. The initial states were chosen so that electrons and holes are in the LUMO and HOMO states, respectively. Usually, the Auger recombination is an efficient process that can "shunt" other photon-assisted processes such as light emission or multiple exciton generation. For instance, in Si crystallites, the Auger rates can be sufficiently high, as the measurements 134,135 and calculations 107,[136][137][138][139] show. Therefore, suppression of the Auger process is an important practical problem.
It turns out that coating with halogens can slow down the Auger processes in Si NCs. Below we present the calculated Auger rates and demonstrate their reduction due to the surface halogenation of Si crystallites. The calculations take into account both the long-range (U 0 (r 1 , r 2 )) and the short-range (U 1 (r 1 , r 2 )) parts of the carriercarrier interaction. The first one is a macroscopic electrostatic field modified by the nanocrystal boundary, while the second one describes at a microscopic level the point charges interaction at short distances: Here, r = |r 1 − r 2 |, A = 1.142, α = 0.82/a B and β = 5/a B 100 , where a B stands for the Bohr radius, and we set 1 = 1 as before. The calculated values of τ −1 A are listed in Table 2 for both eeh and ehh trion annihilation in the Si n X m crystallites considered here. The surface halogenation of Si nanocrystals also slows down the Auger recombination according to the data presented in Table 2. In order to quantitatively analyse the effect of halogenation on the Auger rates, we have calculated the average decimal logarithms of τ −1 A for the Br-, Cl-and H-passivated nanocrystals. In the case of the eeh Auger recombination one has obtained: lg(τ −1 A(H) ) = 11.66, lg(τ −1 A(Cl) ) = 10.55 and lg(τ −1 A(Br) ) = 10.60. In the case of the ehh Auger recombination the rates also decrease but not so strongly: lg(τ −1 A(H) ) = 11.86, lg(τ −1 A(Cl) ) = 11.44 and lg(τ −1 A(Br) ) = 11.20. It should be noted also that strong reduction of the Auger recombination rates was theoretically shown by Califano 140 for GaSb NCs whose surfaces were passivated with atoms of electronegative elements.

B. Capture on Dangling Bonds
There is one more relatively fast process inhibiting efficient light emission from NCs, which is photo-carrier trapping on surface defects called P b -centres or dangling bonds. The P b -centres produce rather deep energy levels within the nanocrystal band gap, which act as electron traps. The capture on neutral dangling bonds is a multi-phonon process, and its rate strongly depends on temperature. The earlier performed calculations 42 for Si crystallites yielded τ −1 C shown in Figure 3 as function of the nanocrystal radius in comparison with the rates of the Auger recombination and radiative transitions.
It can be seen from Figure 3 that for the NC radius greater than ∼1 nm, the capture rate becomes larger than the radiative recombination rate. Moreover, starting from R close to 1.5 nm, the trapping on dangling bonds turns out to be even faster than the Auger recombination. Obviously, for nanocrystals whose radii are greater than 1.5 nm, the nonradiative trapping on the surface defects suppresses all the other competing relaxation processes inside the NC if the surface dangling bonds are not passivated.

C. Phonon-Assisted Relaxation of Hot Carriers
It is well known that highly excited ("hot") electrons, holes or excitons may relax to the lower states in order to minimize the system energy. In bulk semiconductors and QWS, where the electronic energy spectra are continuous, such a relaxation is accompanied by the emission of phonons caused by the interaction between electrons and lattice vibrations. This interaction also plays a major role in QDs and it is generally accepted that the mechanisms are essentially the same in nanocrystals and bulk materials.
The most universal mechanism of coupling between the electrons (or holes) and long-wavelength acoustic phonons is through the volume deformation potential. The bottom of a non-degenerate band (e.g., conduction band in materials with zinc-blend structure) is shifted proportionally to the (local) relative variation of the volume 141 where a c is the bulk deformation potential constant, and u ν denotes the atomic displacement operator corresponding to an acoustic phonon mode ν (it is expressed in terms of phonon creation and annihilation operators in usual way). For optical phonons, there are two mechanisms: (i) a universal optical deformation potential (ODP) coupling and (ii) the Fröhlich one characteristic of polar materials. In bulk semiconductors with cubic structure, for long-wavelength optical phonons, the ODP coupling vanishes by symmetry for any non-degenerate band but it is non-zero for holes near the top of the valence band 141 . Its expression can be found in the book 12 (Chapter 8) or in 142 . The Fröhlich interaction Hamiltonian for electrons is simply given bŷ whereφ ν is the electrostatic potential operator corresponding to an optical phonon mode ν. With the known electron and hole wavefunctions and phonon displacements obtained from lattice dynamics equations, the calculation of the coupling matrix elements is straightforward. Even though the Fröhlich interaction usually dominates in QDs made of polar semiconductors, the ODP mechanism is important for some of them (e.g., InP) as revealed by modelling of QD resonant Raman spectra 142 .
The so-called rigid-ion model is more suitable when describing the electron-phonon interaction in multi-valley semiconductors, such as Si and Ge 143 , where the processes of inter-valley scattering are often important. Within this model, the electron-phonon interaction is described by the following operator (which describes both acoustic and optical phonons), Here, V at (r − R ns ) stands for the electron potential energy in the field of the s-th atom situated in the nth primitive cell of the lattice and R ns is the positionvector of this atom in equilibrium. In practice, atomic pseudopotentials are used.
The discreteness of the energy spectrum in quantum dots implies that the rate of the single-phonon emission should strongly depend on the level spacing 144 ; in particular, electronic transitions causing the phonon emission should be impossible if the spacing is larger than the phonon energy. This simple idea gave rise to the "phonon bottleneck" concept, which predicts the inefficiency of hot carrier relaxation by emission of phonons in QDs 145 . However, this prediction relies on the assumption that the phonon emission is irreversible, with a probability described by the perturbation theory (Fermi's Golden Rule), which may not be reliable in QDs. The electronphonon interaction can be enhanced in nanocrystals because of the spatial confinement of both electrons and phonons and, therefore, multiple scattering processes are important. It means that the electron-phonon interaction in QDs, in principle, must be treated in a nonperturbative and non-adiabatic way, leading to the energy spectra described by polaronic quasi-particle excitations [146][147][148][149][150][151][152] . In other words, virtual transitions between different electronic levels, assisted by phonons and not requiring energy conservation may be important enough to guarantee a significant modification of the exciton energy spectrum and dynamics 12 . Thus, does the phonon bottleneck in QDs really exist? Experimental results still are not completely conclusive. On the one hand, an efficient relaxation of optically created excitons was reported in a number of works studying self-assembled QDs (SAQDs) [153][154][155] , with both PL emission rise time 153,154 and photo-induced intraband absorption decay time 155 below or of the order of 10 ps. Studies performed on nanocrystal QDs 156,157 , where exciton energy level spacings are larger than in SAQDs, also revealed ultrafast intraband relaxation with a characteristic time in the picoseconds' domain, which is characteristic of CdSe and CdTe QDs in general 12 . On the other hand, there are published experimental results indicating that the relaxation of optically created excitons can be slow (in the nanoseconds' range), both in self-assembled 158,159 and nanocrystal 160,161 QDs.
It is important to realize that the polaron concept by itself cannot explain intraband relaxation of carriers in QDs as polaron is a stationary state of an electron (or exciton) coupled to optical phonons. Some additional interactions should therefore be responsible for the polaron relaxation 152 . Several possible mechanisms of hot carrier relaxation in QDs have been proposed: (i) The polaron has a rather short lifetime 162 because of the anharmonic effects that lead to a fast decay of nanocrystal's optical phonons forming the polaron.
(ii) Acoustic phonons can provide the possibility of transitions between different (exciton-) polaron states formed mostly by the interaction with optical phonons; the polaron spectrum is discrete but relatively dense owing to the non-adiabaticity of this interaction 148 . If the acoustic phonon spectrum is continuous, this additional interaction would drive the polaron dynamics towards equilibrium.
(iii) In the strong confinement regime where the electron's kinetic energy is larger than the electronhole interaction, the electron (eventually dressed by phonons and forming the polaron) can relax by an Auger-type mechanism. The excess energy is first transferred from the electron to the hole through their Coulomb interaction and the subsequent hole's cooling occurs via emission of acoustic phonons 163,164 . It can be feasible because the hole level spacings are relatively small in QDs and match the continuum of acoustic phonon energies.
There is no final conclusion concerning the relevance of each of the three mechanisms to the breaking of the phonon bottleneck in QDs; however, the first one seems to be the most popular in the literature, even though calculations yielded too low relaxation rates for this mechanism, of the order of 1 ns −1 , for a typical quantum dot considering anharmonic (Grüneisen) parameters characteristic of bulk materials and the Fröhlich interaction with spherically confined optical phonons and electrons 152 . Moreover, the electron-phonon interaction in nanocrystals of non-polar materials (see below) may be too weak to make the polaron effect important. If this is the case, the relaxation proceeds via multi-phonon processes (like temperature-induced interband excitation of charge carriers in bulk semiconductors). Usually, multi-phonon processes are much slower compared to the single-phonon transitions. Therefore, the bottleneck effect does not disappear completely.
Estimations of the single-phonon relaxation rates in Si nanocrystals 165,166 yield typical values of the order of 10 10 -10 12 s −1 . At the same time, the multiphonon relaxation rates calculated within the Huang-Rhys model 42,167,168 for various transitions vary from 10 7 s −1 to 10 11 s −1 for nanocrystals whose diameters do not exceed 5 nm. Note that the rates sharply reduce as the nanocrystal diameter decreases. Similar values of the multi-phonon relaxation rates were obtained experimentally for InGaAs/GaAs SAQDs 169 . Therefore, the precise reason for very fast exciton dynamics characteristic of CdSe and CdTe colloidal QDs remains an open question, in our view.

IV. MULTIPLE EXCITON GENERATION
As demonstrated above, the surface halogenation suppresses both radiative and Auger recombination in nanocrystals. These processes are reverse ones with respect to the carrier multiplication (or multi-exciton generation): they tend to decrease the number of excitons in a system, while the process of multi-exciton generation, shown schematically in Figure 4, has an opposite trend. Initially, a high-energy photon creates a highly excited electron-hole pair which then reduces its energy creating one more electron-hole pair with lower energy. As a result, two (or even more) excitons can arise in the system after absorption of a single photon.
The multi-exciton generation is a fundamental process for photovoltaics, where light energy transforms into electric current. Its occurrence in NCs has been experimentally confirmed 43,134,[170][171][172][173] . In order to be more efficient, this process should be faster than other competitive processes taking place along with exciton generation, such as inter-band radiative recombination or Auger recombination. From this point of view, slowing the latter down caused by surface halogenation is an extremely positive factor. It is important to understand how halogen passivation influences the multi-exciton generation itself.
The rates of the multi-exciton generation in the Si 317 X 172 crystallite (X = H, Cl, and Br) can be calculated using the following relation, whereÛ =Û 0 +Û 1 , as before, while Ψ i1i2 and Ψ f1f2 are the products of the single-particle Kohn-Sham wave functions of the initial or final electron states participating in the transition, as shown in Figure 4. Here, we present calculated the rates for the Si 317 X 172 crystallite within the range of excess energy 0 < ∆E < 0.6 eV, where ∆E = E i1 − E LU M O − E g for the process initiated by highly excited electron, and ∆E = E HOM O −E f1 −E g for the hole-initiated process. The calculated rates are shown in Figure 5. It is evident that the rates rise globally as the excess energy increases, because of the considerable increase in the number of possible states participating in the transitions with increasing ∆E, which opens up many new channels for the realization of exciton generation.
It is important to emphasize that the bromination of a Si crystallite increases the exciton generation rates compared to an H-passivated crystallite, especially if the process is initiated by a highly excited hole. The rates of exciton generation in chlorinated Si crystallite turn out to be lower than those in the hydrogenated one at small excess energies. Meanwhile, upon approaching ∆E ∼ 0.5 eV, τ −1 G in the Si 317 Cl 172 crystallite increases and tends to the typical values observed in Si 317 H 172 crystallite. Accordingly, it is possible to conclude that halogen coating of Si crystallites, at least, does not reduce their ability to generate excitons, particularly when the excess energies are not too small. This is in contrast with the radiative and Auger recombination processes, where the rates became substantially lower due to the halogenation.  It means that the halogenation of Si NCs can increase the efficiency of the photon-to-exciton conversion, which is defined by an excess of the number of created excitons (n) over the number of absorbed photons (N ) 43,45 : η = n/N > 1 (internal quantum efficiency). There are many experimental works in which external quantum efficiency was measured in Si crystallites 134,[170][171][172][173] , as well as in crystallites of IV-VI, II-VI or III-V semiconductors 43,[174][175][176][177][178] . The authors reported on the observation of multi-exciton generation in the investigated systems.
On the theoretical side, the consideration of the exciton kinetics in the halogen-coated Si crystallites 132 has revealed strong dependence of η on the quantitative relationship between the rates τ −1 G and τ −1 A . According to the obtained results the decrease in the Auger rate (caused by the halogenation), and its absence in the multi-exciton generation rate, is accompanied by a gradual increase in the quantum efficiency of the order of a few tens of percentage points.

V. FRET IN ENSEMBLES OF NANOCRYSTALS
All the processes considered above may occur in isolated nanocrystals. Meanwhile, usually in experiments, as was already pointed out in the Introduction, one deals with ensembles of NCs, where a non-radiative energy exchange between them takes place and strongly influences the ensemble photoluminescence [54][55][56][57][58][59][179][180][181] . Such an energy transfer occurs by means of the electron tunnelling if the NCs are connected 52 or through the Förster-type migration of excitons 50 if the NCs are separated in space. Below we consider the latter mechanism in some detail.
The Förster resonant exciton transfer (FRET) takes place mainly through the dipole-dipole interaction 241 of two QDs (a donor and an acceptor). The Quantum Electrodynamics theory of FRET developed in 183 reproduces the results obtained in a simpler way by Förster 50 and Dexter 51 who considered the electrostatic interaction of two dipoles in two adjacent crystallites 184 , where b is the inter-crystallite centre-to-centre vector. In one crystallite, the electron-hole pair annihilates and transfers its energy into a neighbouring crystallite where a new electron-hole pair is excited. Thus, virtual transfer of excitons between two crystallites can be realized without charge transfer. In order to calculate the rate of the Förster exciton transition from a QD with a radius R 1 into a neighbouring one with a radius R 2 , we use, as before, the Fermi's Golden Rule: Here,V is the dipole interaction operator identical to (10), Ψ i = ψ c (r 1 )ψ v (r 2 ) is the wave function of the initial two-particle state with the energy coinciding with the energy gap of the first crystallite, E g (R 1 ), and Γ is a phenomenological damping parameter. Initially, there is one electron in the conduction band of the first crystallite, with the wave function ψ c (r 1 ), while a hole exists in the valence band. In the second QD, the electron occupies a valence band state described by the wave function ψ v (r 2 ). In the final state, the system has the wave function Ψ f = ψ v (r 1 )ψ c (r 2 ), corresponding to the electronhole pair transferred to the second QD.
The result is essentially the same as if one considered the interaction between two transient point dipoles, d D and d A (located at points r 1 and r 2 ) assuming that they are stationary. The square of the matrix element in (11) (the donor-acceptor coupling parameter) can be written as 185 where ν is a factor of the order of unity that takes into account the relative orientation of the dipoles 186 , ν = 2 3 if these orientations are completely random.
The broadening of the exciton transfer resonance, Γ can be related to the spectral overlap of the emission and absorption spectra of the donor and acceptor, respectively 51 , and the transfer rate can be expressed as 187 where γ (D) 0 is the spontaneous emission rate of the donor, I A (ω) is the absorption lineshape function of the acceptor, L D (ω) is the emission lineshape function of the donor (both are normalized to unity) and Q A is the frequencyintegrated absorption cross section of the acceptor QD, with γ (A) 0 being the spontaneous emission rate of the acceptor. Notice that the depolarization factor κ 2 , Equation (3), which distinguishes QD donors and acceptors from, e.g., molecules, has been included explicitly in Equation (13).
The FRET rate decreases rather quickly with the donor-acceptor distance, k F = γ (D) 0 ( b0 b ) 6 , where b 0 is a parameter of the dimension of length named "Förster radius" whose definition is clear from Equation (13).
In ensembles formed of direct-band-gap II-VI or III-V semiconductor crystallites (QDs), the exciton transfer has the rates ∼10 8 -10 9 ·s −1 , as measurements 53,54,56,57 and calculations 184,186,188 show. Notice that k F depends on b and the above values probably correspond to somewhat different donor-acceptor distances. Experimental observations of the exciton transfer were carried out with ensembles of closely packed monodisperse CdSe nanocrystals as well as of two-size three-dimensional mixtures 56 and bilayered CdSe nanocrystal systems 57 . It was shown that the inter-layer transfer in bilayered ensembles with controlled donor-acceptor separation turns out to be more efficient than in the 3D ensembles of the monodisperse and the two-size-mixed crystallites 57 . Note that the measured and computed FRET rates are of the same order of magnitude as the radiative recombination rates in high-density ensembles of colloidal II-VI crystallites. Moreover, a few measurements of the Förster radius have been performed by controlling the distance between two different groups of QDs of different size, acting as donors and acceptors 189 , or by preparing a homogeneously blended solid-state films 190 composed of two groups of dots. In both works, CdSe/ZnS core-shell QDs were used and the reported results for the Förster radius are 14-22 nm and 8-9 nm, respectively; notice that the former work 189 used larger nanocrystals. Despite the uncertainty in experimental conditions and difficulty to evaluate the spectral overlap in Equation (13) for individual QDs, there is a consensus that R F typically is of the order of 10 nm for highly luminescent colloidal nanocrystals.
In ensembles of silicon crystallites the exciton transfer turns out to be much slower than the radiative recombination, its rates are two to three orders of magnitude lower 184,191,192 (smaller than ∼10 3 ·s −1 ). Doping of silicon nanocrystals with phosphorus allows to increase the rates up to the values comparable with those of radiative recombination 193 ( < ∼ 10 7 s −1 for nanocrystal radius R > ∼ 1 nm). Nevertheless, these values remain still smaller than the FRET rates for colloidal QDs made of direct-gap semiconductors.
For the sake of completeness, Förster-type resonant energy transfer from a QD exciton can occur also when the final elementary excitations have a different nature, for instance, they are phonons or molecular vibrations. The latter can involve ligand molecules present in the vicinity of a colloidal QD or solvent molecules and is particularly important for quantum dots emitting in the infrared 194,195 . Such a process was termed electronic-to-vibrational energy transfer (EVET) 194 . The importance of EVET was demonstrated by comparing the luminescence properties of HgTe QDs dissolved in two chemically identical solvents: H 2 O and D 2 O 195 .

VI. QD EMITTERS NEAR A FLAT INTERFACE
Now, we shall discuss the influence of a flat interface between two media on the emission and FRET rates for a QD emitter located at a certain distance from it. If the characteristic distances (such as that between the emitter and the interface) are small compared to the EM wavelength, one can treat the problem in the electrostatic approximation where one neglects both the retardation effects and the magnetic field associated with the electric field present in the media; then, the image dipole method can be used to take into account the effect of surface polarization induced by the dipole 196 . A more general approach consists in using the dyadic Green's function formalism 197 . If both dielectrics are dispersionless, the spontaneous decay and FRET are affected through the photonic density of states (DOS) renormalization due to the reflection of electromagnetic (EM) waves at the interface. However, even in this geometrically simplest situation, there is a less trivial effect produced by EM waves created by the polarization in the second medium, totally reflected at the interface so that only an exponentially decreasing field amplitude of theirs reaches the dipole. Yet, this effect can be dominating at small distances, not only for radiative decay 198 but also for FRET. When the second medium is a metal, several contributions to both decay and transfer processes arise that can be associated with (i) propagating EM waves (so-called radiative losses), affected by the presence of the second medium; (ii) coupling to propagating surface plasmons (SPs); and (iii) Ohmic losses (when the exciton energy is irreversibly transferred to heat via electron scattering in the metal).
Even though the subject has been studied during several decades 71,198,199 and is described in topical reviews 200 and books 79 , we find some significant effects that were overlooked in previously published discussions and, in our opinion, are important for potential applications of FRET in photonics and energy harvesting. In particular, we will show that (1) if a QD embedded in a dielectric matrix is located near an interface with another dielectric (a substrate) having a lower dielectric constant, its emission is polarized parallel to the interface; (2) in the opposite case of substrate's dielectric constant higher than that of the matrix, the emission and transfer rates are strongly enhanced (without dissipation) near the interface; and (3) the resonant coupling between SPs propagating along a metal/dielectric interface and excitons confined in QDs located at a distance of the order of the light wavelength from the interface, can be used for long-range FRET. We present these rates calculated for typical colloidal nanocrystal QDs.

A. Radiative Lifetime Near Interface
We shall consider the QD emitter as a point dipole. 242 The radiative decay rate in the presence of other bodies is determined by the total field that acts on the emitter (created by itself and scattered by the bodies) and given by 79 where E 0 is the dipole field in infinite space. Here, d is the complex-conjugate of the classical dipole moment, the necessary adaptation for the quantum mechanical transition dipole matrix element consists in multiplying it by 2 in the final result, d → 2d 79 . The electric component of the scattered field created by an emitting dipole located in the origin of the reference frame can be written in terms of the Green's dyadic 79 , which is specific of system's geometry: The 3 × 3 dyadic matrix G(r , r) is determined by the Fresnel coefficients, r (p) and r (s) , and its explicit expression can be found in 197 .
From Equation (15), one obtains for a dipole perpendicular to the interface (i.e., "vertical", see inset in Figure  6): and for a dipole parallel to the interface (i.e., "horizontal"): Here, γ 0 is the decay rate in an infinite medium with dielectric constant ε 1 , given by Equation (1). The integrals in Equations (17) and (18) with respect to the normalized in-plane wavevector, s = q/k 1 , can be divided into two parts, one from 0 to 1 corresponding to propagating waves in the upper half-space and the other from 1 to ∞ representing evanescent waves with imaginary wavevector component along z axis perpendicular to the interface, k 1z . The latter type of waves exist in two cases: (i) if ε 1 < ε 2 and both dielectric constants are positive and (ii) if ε 2 ≡ Re(ε 2 ) < 0, i.e., when the second medium is a metal.
In the case of two non-dispersive dielectrics, γ given by Equations (17) and (18) represents only radiative losses of dipole's energy, renormalized by the back action of the scattered waves. Yet, for ε 1 < ε 2 there are EM waves (excited by the dipole) propagating in the lower half-space (with real z-component of the wavevector inside medium 2, k 2z ) and experiencing total internal reflection at the interface with the upper medium, possessing imaginary  Figure 6: Dependence of the radiative decay rate, γ, on the dielectric constant ratio ε1/ε2 for two two different distances to the interface between them, k1h = 0.1 and 0.2. The dotted shows the dependence when the contribution of evanescent waves (that exists for 0 < ε1/ε2 < 1) is neglected; it is independent of h. The inset shows two orthogonal orientations of the dipole moment. The upper plot is for the vertical dipole and the lower for the horizontal orientation. k 1z . The Fresnel coefficient r (p) has an imaginary part associated with these waves. If the interaction of the dipole with these evanescent waves is neglected (dotted lines in Figure 6), i.e., we set the integral from 1 to ∞ in Equations (17) and (18) equal to zero, then γ becomes independent of dipole's distance from the interface. As can be seen from this figure, the emission decay rate is strongly enhanced, γ >> γ 0 , at small distances (k 1 h < 1) for both dipole orientations because of the strong interaction with the evanescent waves. It has been pointed out for plane interfaces by Lukosz and Kunz 198 , and it is similar to the coupling to the whispering gallery modes (a kind of surface EM waves) observed for emitters placed inside a micrometer-size dielectric sphere 202,203 . As we will show in the next section, this coupling can be used for FRET enhancement.
For ε 1 > ε 2 only propagating waves in the upper medium (with k 1z real) are present and γ does not depend on h. As can be seen from Figure 6, the emission rate of a dipole oriented perpendicular to the interface is suppressed (γ → 0), while that parallel to the interface has γ close to its value in an infinite dielectric, γ 0 . As known, the emission of a spherical QD made of a semiconductor material with cubic crystal lattice does not have any preferential polarization in empty space (if excited by non-polarized light). If the dot is embedded in a matrix with sufficiently high ε 1 and located near the surface, one can expect its luminescence to be strongly polarized. Note that this is not a quenching effect because there is no dissipation in the system (both ε 1 and ε 2 are real).

B. Non-Radiative Losses to a Metal Substrate
If the substrate is lossy, e.g., a metal, several additional decay channels arise. Mathematically, they are related to the imaginary part of the Fresnel coefficients. Direct evaluation of the integrals in Equations (17) and (18) yields the overall effect of all these channels. However, their relative contributions depend on the distance and the real and imaginary parts of ε 2 . At very small distances (k 1 h << 1), the dominating contribution is due to Ohmic losses 71,200 , The coupling to propagating surface plasmons mathematically is described by the pole of r (p) , which yields an additional imaginary part of the integral from 1 to ∞ in Equations (17) and (18). Indeed, the equation (r (p) ) −1 = 0 determines the dispersion relation of ppolarized SPs. The other Fresnel coefficient, r (s) , has no poles, and accordingly there are no s-polarized SPs at a metal-dielectric interface. If we neglect the imaginary part of ε 2 , the SP wavevector is real. The resonant coupling occurs for q sp = ω c ε 2 ε 1 ε 2 + ε 1 (where ε 2 ε 1 < 0 and ε 2 + ε 1 < 0). The pole contribution to the integral from 1 to ∞ in Equations (17) and (18) is calculated explicitly: where u = ε 2 /ε 1 and s sp = q sp /k 1 . Taking ε 2 (ω) within the simple Drude model for bulk gold, we calculated γ Ol , γ sp , and the radiative decay rate (neglecting ε 2 ). They are presented in Figure 7. As expected, at small distances the Ohmic losses dominate (see Figure 7, lower panel). However, as the distance to the interface increases, for k 1 h ∼ 1, the coupling to SPs becomes the main mechanism that determines the lifetime (see Figure 7, upper panel). Although this effect decreases with the distance exponentially, this decrease can be rather slow, as |ε 2 |/ε 1 >> 1 and s sp ∼ 1, so it is nearly constant over tens of nanometres. Of course, at still large distances γ tot → γ 0 .
Thus, for intermediate distances of the order of a hundred nanometres the radiative decay is mostly due to the coupling to SPs within a broad spectral range (approximately from 2.5 to 4.5 eV). In this range, relevant for colloidal QDs, one can expect also plasmonenhanced FRET. Even though the enhancement is moderate, the coupling to SPs is not dissipative unless we approach the electrostatic SP resonance frequency, ω spr = ω p / ε ∞ 2 + ε 1 (here, ω p is the bulk plasma frquency and ε ∞ 2 is the background dielectric constant of the metal). In other words, excitation can be transferred from the QD to SPs and back many times without dissipation. It can lead to strong coupling regime characteristic of plasmonic microcavities 41,204 .

C. Energy Transfer to a 2D Material
Graphene, an atomic-thick monolayer of carbon, is a semimetal (or gapless semiconductor) with unusual electronic properties first demonstrated by Geim and Novoselov 205 . Low-energy excitations in graphene are massless, chiral, Dirac fermions. In neutral graphene, the chemical potential (hereafter called Fermi level) crosses exactly the Dirac point 206 . Pristine graphene is transparent in a broad spectral range from the infrared (IR) to the ultraviolet (UV), with the residual absorption of ≈2.3% originating from interband transitions 207 . The Fermi level, E F , can be shifted by up to ≈1 eV with respect to the Dirac point by applying an appropriate gate voltage. The usual doping via adsorbed impurities is also possible 206 . One can say that graphene is a transparent conductor with a tunable conductivity, both static and frequency-dependent (for the frequencies up to the mid-IR) 208 . It supports propagating surface plasmons whose dispersion relation can be controlled via gate voltage 209 . Interestingly, a specific type of quantum dots can form in graphene, namely, mass profile QDs, a system within current experimental reach 210 , for which FRET can also be important 211 .
Numerous experiments have shown that the nanoparticle or dye molecule emission is quenched in a broad spectral range by a single graphene sheet 39,[212][213][214][215] . It is understandable for the lower-energy part of the spectrum, as graphene is a conductor, although with a relatively low free carrier concentrations (which can reach 10 13 cm −2 ), with the plasma frequency in the THz-to-IR range; it also is supported by calculations. However, the situation is more complex than in the case of a point emitter in the vicinity of a usual metal.
In addition to the Ohmic losses and energy transfer to surface plasmons, there is an additional mechanism that may be named exciton transfer (ET) to graphene. It is similar to FRET and consists in energy transfer to electron-hole pairs generated by inter-band transitions, symmetric with respect to the Dirac point (see Figure 8). Doped graphene absorbs EM radiation due to inter-band transitions only for energies above 2E F , therefore, it also applies to the ET. The main contribution of this mechanism scales with distance as h −4 . It can be demonstrated in a rather simple way that elucidates its connection to FRET.
Let us assume that acceptor dipoles are distributed within a plane (representing the graphene sheet), then the ET rate can be obtained by integrating (13) over this plane and written as follows, where ρ A is the surface density of acceptor dipoles. The emitter lineshape may be taken as a δ-function. As ρ A Q A I A = β(ω) is the absorbance of the graphene, we obtain This formula can be derived in a more rigorous way 216 , with just slightly different numerical coefficients.
The h −4 scaling law was verified in several works and it seems to be obeyed for molecular emitters (down to h ≈ 5 nm) 212 but not so much for QDs where deviations were found for small distances 213 . The above theory does not take into account the finite size of the QDs. As a first approximation, one can assume that the point dipole is located at the centre of the QD. In this case, the distance of the dipole from the graphene sheet can be approximated as h * = (h + R QD + s) where R QD is the radius of the nanocrystal core of the QD and s is thickness of the organic shell for a colloidal QD. The latter can be estimated as the chain length of the capping agent, s ∼ 1.5 nm 12 . Using this correction, h → h * , indeed improves the agreement between the theory and experiment 213 . However, one needs to bear in mind also the possibility of charge transfer from an excited QD to graphene if the distance is sufficiently small for tunnelling. Indeed, the demonstration of hybrid graphene-QD phototransistors 217 and solar cells 218 where the light is absorbed in the QDs and changes the electric current in graphene implies that such processes can take place in specially designed structures.
Calculations including both intra-band (i.e., Drude plasmons) and inter-band transitions in graphene show that there can be a no-quenching spectral window for an emitter over a strongly doped graphene sheet 72 . Indeed, the graphene absorbance, β(ω), due to Drude plasmons decreases with the frequency, while the inter-band transitions contribute to β (and, therefore, to the quenching) only for ω > 2E F /h 209 . As the Fermi level in graphene can be tuned electrically, it opens an interesting possibility of electrical switching of the QD emission by controlling the quenching rate; it has been demonstrated experimentally in 219,220 .
Energy transfer (ET) from QDs to other 2D materials, namely, few-atomic-layer-thick semiconductors (transition metal dichalcogenides (TMDs)) has also been investigated in the recent years 215,[221][222][223][224] . These materials support robust excitons, which determine their optical properties in the visible range (along with the inter-band transitions) and they depend strongly on the number of monolayers in the material 225 . The ET from a QD exciton to the continuum of states above the band gap in the 2D material should work similar to the case of graphene (Equation (23)). However, it was found that the donor QD emission is quenched most strongly near monolayer MoS 2 and then the effect decreases with the number of monolayers (presumably keeping the QD distance to the semiconductor surface fixed), in contrast with the case of few-layer graphene 215,222 . Meanwhile, the tendency similar to graphene (quenching increases with the number of monolayers) was found for the less studied TMD SnS 2 223 . On the theoretical side, it has been suggested that the total spontaneous emission rate of a quantum emitter can be enhanced several orders of magnitude due to the excitation of surface exciton-polariton modes supported by the 2D semiconductor 224 . In the vicinity of an excitonic transition, the real part of the dielectric function of the semiconductor material is negative, like in a metal, so one can indeed expect an enhancement of the local photonic DOS and, consequently, the Purcell effect 226 . However, it would require matching of the energies of the QD and 2D excitons.

D. FRET between QDs Near Interface
In the vicinity of a polarisable material, such as graphene, the dipole-dipole interaction responsible for the FRET between two QDs becomes renormalized because of the additional interaction between the induced charges. In the electrostatic limit, these additional interactions can be modelled with image dipoles (see Figure  9). This physical mechanism, of course, is not unique to graphene and can take place in the vicinity of plasmonic nanoparticles 36,37,227 . The effect of plasmonic enhancement of FRET has been detected experimentally by observing a decrease of the emission decay rate of donor QDs and a corresponding increase of the luminescence intensity for acceptor dots when they are placed in the vicinity of gold nanoparticles 37 . Of course, if the concentration of the nanoparticles becomes too high, the emission is quenched for both donor and acceptor QDs because of the energy transfer to lossy modes in the metal (the Ohmic losses, Section VI B). Graphene is a conductor with electrically controllable conductivity and, in principle, it promises a possibility of controlling the inter-dot FRET rate. The exciton energy transfer between two QDs located near a graphenecovered interface occurs not just through the direct dipole-dipole interaction but also through polarization charges that they induce on the interface. One can think that the image dipoles of Figure 9 have magnitudes that depend on the optical conductivity of graphene and, therefore, on its Fermi energy 196 . This interaction can be written in the form withT =T (0) +T (1) , wherê is the usual dipole-dipole interaction tensor (n = b/b) and T (1) (b) = d 2 q (2π) 2 (−ie q +e z )⊗(ie q +e z )qA(q)e −2qh+iq·b (24) is the part due to the image dipoles. In Equation (24), e q = q/q, e z is the unit vector along z-axis, and σ g (ω) is the optical conductivity of graphene. The components ofT (1) have an oscillatory dependence upon the inter-dot distance 228 , which arises from the Bessel functions appearing as the result of the angular integration in Equation (24). It means that instead of the usual monotonic b −6 dependence of the FRET rate one may have an oscillatory behaviour, controllable by the gate voltage through the graphene Fermi energy (which determines σ g (ω)). Beyond the non-monotonic dependence upon the inter-dot distance, the emission/transfer frequency also should affect it in a rather complex way. If this interaction is strong enough, the emitters (QDs in our case) can be coupled to create a collective radiative mode, a phenomenon called superradiance 229 , which has been observed for epitaxial quantum dots 230 . As the coupling between the emitters can be strongly enhanced or suppressed in the vicinity of graphene, depending on the distance between them, the frequency and the graphene Fermi energy, both superradiance and subradiance regimes can be expected 228 . So far, we are not aware of an experimental demonstration of such effects. Losses associated with the real part of graphene's optical conductivity can be one possible reason for this. The influence of a polarisable surface (such as graphene) on the irreversible FRET should be an easierto-observe effect than the superradiance. It can be seen as a surface plasmons' effect. Huge FRET enhancement factors of the order of 10 6 (compared to vacuum) have been predicted for two dipole emitters on a graphene sheet, based on numerical electrodynamics calculations 40 , reaching a maximum when the distance between the dipoles equals twice the graphene plasmon propagation length L p . This prediction is at variance with the results presented in 228 according to which the effect must be considerably smaller and varying on scale of the surface plasmon wavelength, λ p << L p .

VII. CONCLUDING REMARKS
To summarize, we would like to emphasize, once again, that the radiative properties and efficiency of light emission in semiconductor nanostructures are determined not only by the radiative transitions, but also by various nonradiative processes which can be even more intense than the radiative transitions themselves. As a result, in these cases, the radiative recombination can be significantly suppressed. On the other hand, even the non-radiative processes can be used in a constructive way. For instance, the exciton migration is accompanied by a non-radiative energy transfer (FRET) that can be channelled in a desired direction via creating a certain "architecture" of the nanocrystals in the ensemble, which allows to concentrate and illuminate the energy inside a given area 59,191 or direct it to a certain layer of a funnel-type heterostructure of QD monolayers 231 . Exciton energy transfer and emission with upconversion because of simultaneous absorption of optical phonons has been suggested for QDassisted cooling 232 . The multi-exciton effects make it possible to efficiently transform the absorbed photons into the rising number of electron-hole pairs capable of participating in the electric current. Such processes underlie the operating principle of photosynthetic light harvesting systems 233 and can be mimicked in solar cells 234 .
The exciton transport can be strongly enhanced if me-diated by polaritons arising under strong coupling of the excitons to light in microcavities 235,236 . Combining entities supporting localized excitons, such as QDs, dye molecules or J-aggregates, with microcavities in order to achieve the strong coupling regime is an important area of research envisaging various applications 41 , in particular, in the context of controllable quantum emitters. A plasmonic surface considered in Section VI can be seen as a near-field microcavity if the lossy channels are not the dominating ones. A detailed discussion of the strong coupling between surface plasmon-polaritons and quantum emitters can be found elsewhere 237 .

Acknowledgments
Funding from the Ministry of Science and Higher Education of the Russian Federation (State Assignment No 0729-2020-0058), the European Commission within the project "Graphene-Driven Revolutions in ICT and Beyond" (Ref. No. 696656), from the Portuguese Foundation for Science and Technology (FCT) in the framework of the PTDC/NAN-OPT/29265/2017 "Towards high speed optical devices by exploiting the unique electronic properties of engineered 2D materials" project and the Strategic Funding UID/FIS/04650/2019 is gratefully acknowledged.