Self-assembled quantum dot (QD) solids are a highly tunable class of materials with a wide range of applications in solid-state electronics and optoelectronic devices. In this perspective, the authors highlight how the presence of microscopic disorder in these materials can influence their macroscopic optoelectronic properties. Specifically, they consider the dynamics of excitons in energetically disordered QD solids using a theoretical model framework for both localized and delocalized excitonic regimes. In both cases, they emphasize the tendency of energetic disorder to promote nonequilibrium relaxation dynamics and discuss how the signatures of these nonequilibrium effects manifest in time-dependent spectral measurements. Moreover, they describe the connection between the microscopic dynamics of excitons within the material and the measurement of material specific parameters, such as emission linewidth broadening and energetic dissipation rate.

## I. INTRODUCTION

The optoelectronic properties of colloidal quantum dots (QDs) depend sensitively on their size, shape, and chemical composition.^{1,2} This dependence has inspired the development of a class of solid materials made up of self-assembled QDs that exhibit highly tunable optoelectronic properties. This tunability has been leveraged to enable a wide range of solid-state applications, such as light-emitting diodes,^{3,4} solar cells,^{5} lasers,^{6} photodetectors,^{7,8} and luminescent solar concentrators.^{9} Notably, the optoelectronic properties of QD solids can also depend on the spatial arrangement of QDs within the material. However, understanding this dependence has been a challenge because it conveys through collective interactions that are especially sensitive to heterogeneity in the QD population, arising due to the process of QD synthesis, and in the spatial distribution of QDs, arising due to the process of self-assembly.

In a typical QD solid, excited electrons and holes tend to localize on individual QDs. These oppositely charged carriers can colocalize on each QD to form neutral quasiparticles known as excitons. The dynamics of a localized exciton in a QD solid involve a series of hops, whereby the exciton moves from one QD to another through the resonant energy transfer process. The theory to describe this type of exciton dynamics is that of Förster resonance energy transfer or FRET.^{10,11} Many studies, both experimental and theoretical, have used FRET to gain a better understanding of the role of exciton dynamics in the macroscopic properties of QD solids.^{12–22}

The limitations in energy transport that are implied by FRET-like exciton dynamics can be overcome by exploiting the quantum mechanical effect of exciton delocalization. When an exciton delocalizes over many individual QDs, it can more readily explore space and can undergo enhanced supertransfer through the constructive interference of individual QD transition dipole moments. Unfortunately, achieving exciton delocalization in QD solids has emerged as a significant challenge due to the persistence of weak inter-QD electronic coupling. This weak coupling is a consequence of insulating ligands that passivate QD surfaces as well as mismatched energetic resonances arising through variations in QD sizes. Theoretical simulations and model studies are essential to developing QD solids that support delocalized excitations.

Here, we present a general model framework, based on the principles of Förster theory, for simulating exciton dynamics in QD solids. We use this framework to study the effects of energetic disorder on the dynamics of localized and delocalized excitons. This model demonstrates how energetic disorder leads to nonequilibrium effects, and how those effects manifest in experiment. We then apply this framework to analyze spectrally resolved transient photoluminescence measurements of CdSe QD solids. We then identify how material properties of QDs, such as energetic disorder, energy dissipation upon optical excitation, and emission linewidth broadening, can influence exciton diffusivity as well as the transient shift in average emission energy. By extending this framework to the case where inter-QD coupling is larger, we are able to highlight the enhancement in exciton transport that can arise through exciton delocalization. Finally, we summarize and discuss the potential of delocalization to improve the performance of technologies based on QD solids.

## II. MODEL OF ENERGY TRANSFER IN QUANTUM DOT SOLIDS

Exciton dynamics in QD solids can be understood by considering how an exciton is transferred from one QD to another. Förster theory provides a basis for describing this process in the incoherent limit.^{10,11} In Förster theory, the coupling that drives energy transfer originates from the interaction of the transition dipole vectors of the donor and acceptor QDs. Thus, the rate of energy transfer between any two QDs scales as $1/ d 6$, where $d$ is the separation distance between the QDs.^{16,17} Due to the detailed balance condition, Förster theory also predicts faster rates for exciton transfer that is downhill in energy, leading to a transient redshift of the average emission energy in inhomogeneously broadened QD solids.^{23,24} Such variations in QD energies have been attributed to size, shape, and stoichiometric variations between individual QDs.^{25,26} This transient redshift has been observed in QD solids via spectrally resolved transient photoluminescence measurements.^{12–14,18–20}

The downhill energetic migration of excitons in QD solids has more subtle dynamic consequences due to the fact that thermalized excitons have, on average, fewer possible downhill transitions than nonthermalized excitons. This can lead to an exciton diffusivity that decreases over time, which has been observed in studies of time-resolved optical microscopy applied to inhomogeneously broadened QD solids,^{23} agreeing with theoretical predictions of incoherent transport over disordered energy landscape.^{24} These experiments reveal that Förster theory correctly predicts scaling parameters that affect the energy transfer rate in QD solids.^{16,17,20,23,27}

To understand the dynamics of excitons in QD solids, let us consider a model QD solid system that is illustrated in Fig. 1(a). This model system includes QDs assembled in a two-dimensional hexagonally closed packed lattice. QDs are inhomogeneously broadened such that exciton energy at a given QD site is drawn randomly from a Gaussian distribution with mean $ \epsilon \xaf$ and standard deviation $ \sigma ih$ (the inhomogeneous linewidth). Moreover, each QD is assigned a fixed transition dipole vector, $ \mu ^$, oriented randomly on the surface of a unit sphere, assuming that orientations of the transition dipole moments are isotropic.

^{28}we can simplify Eq. (2) as

The two-dimensional model QD system, depicted in Fig. 1(a), was specifically parameterized to reflect the characteristic of the experimental system investigated by Akselrod *et al.* in Ref. 23, in particular, a thin film of colloidal CdSe/ZnCdS core–shell quantum dots having CdSe cores with an average diameter of 4.2 nm and Zn $ 0.5$Cd $ 0.5$S shells with an average thickness of less than 1 nm. The surfaces of these colloidal QDs were capped with short aromatic capping ligands (benzylphosphonic acid). We find that a model system of two-dimensional QD solids is sufficient to understand the general qualitative effects of material specific parameters such as energetic disorder on exciton transport that would occur in three-dimensional QD solids used in most experimental settings.^{23}

By applying the FRET rate equation [Eq. (1)] to this model, it is possible to generate an entire transition rate matrix for excitons within the QD solid. This rate matrix can then be used to simulate the energy transport properties of model materials or to aid in the interpretation of experiments. To accomplish this, we use a chemical master equation as described in Sec. III.

## III. LOCALIZED EXCITON HOPPING PICTURE

^{28}we present an analytical expression for average energy at equilibrium in the case of FRET for the first time. We recall that at equilibrium, exciton population satisfies the detailed balance condition

^{23}

We apply our model to analyze spectrally resolved transient photoluminescence measurements of CdSe/ZnCdS core–shell colloidal QD assembly, as shown in Fig. 2. For details about the sample and the measurement technique, we refer readers to Akselrod *et al*.^{23} Briefly, QDs are excited at 405 nm (3.06 eV) using a laser diode producing pulses $\u223c500$ ps in duration with a repetition rate of 10 MHz and operating at a low laser fluence to probe dynamics of noninteracting excitons. As seen in Fig. 2(b), emission spectra of this QD sample have an asymmetric lineshape, with an elongated tail toward low energy ( $<2.06$ eV). Moreover, the ratio of the photoluminescence intensity at the low energy tail to that at high energy (between 2.1 and 2.2 eV) increases over time. This asymmetric emission lineshape has been also observed in other colloidal QD systems, in which the low energy tail is attributed to subband-edge states whose origin is under debate.^{29}

In our analysis, we only consider the band-edge exciton state whose emission peak is fit to a normal distribution as illustrated in Fig. 2. This emission peak has a total linewidth of about 28 meV that stays relatively constant throughout the measurement. The average exciton energy saturates to a value that is 12 meV lower than the initial value within the first 20 ns. Based on these observations and using Eq. (10), we estimate inhomogeneous and homogeneous linewidths of CdSe QD from the ensemble measurements to be 14 and 25 meV, respectively, provided that Stokes shift of this sample has been measured to be 38 meV.^{23}

The dynamic redshift of emission energy due to energetic disorder reveals information about the spatiotemporal dynamics of excitons in colloidal QD solids. Solving Eq. (5) using a kinetic Monte Carlo algorithm,^{24} we relate the transient energetic to exciton diffusivity by varying inhomogeneous linewidth, homogeneous linewidth, and Stokes shift as shown in Fig. 3. Since energetic disorder leads to a time-dependent diffusivity,^{24} we report mean exciton diffusivity to be the value determined when exciton population equilibrates to a thermalized distribution. Based on results plotted in Fig. 3, mean exciton diffusivity can be enhanced by decreasing the net loss in initial excitation energy, which can be achieved by reducing energetic disorder (inhomogeneous linewidth), increasing available thermal energy for exciton hopping (homogeneous linewidth), and minimizing Stokes shift. Figure 3(b) highlights that increasing homogeneous emission linewidth can mitigate the net negative effect of energy dissipation on mean exciton diffusivity until $ \sigma h\u2248 \sigma ih \Delta ss$, beyond which diffusivity remains constant. For transient redshifts [Figs. 3(d)–3(f)], the final average energy follows the prediction by Eq. (10) except in cases where $ \sigma h\u226a \sigma ih$ [Fig. 3(e)] and $ \Delta ss\u226b \sigma h 2 + \sigma ih 2$ [Fig. 3(f)]. In these situations, excitons never reach the thermal equilibrium because the probability of exciton hopping to neighboring QDs is lower compared to that of exciton decaying back to the electronic ground state.

## IV. SUPERTRANSFER FOR DELOCALIZED EXCITONS

Due to inhomogeneous emission, low oscillator strength, and fast exciton dephasing,^{30–32} excitons in II–VI QDs are thought to be localized on individual colloidal QDs with diffusion length measured between 5 and 35 nm.^{23,33,34} Recent developments in colloidal QD synthesis, however, have achieved ensemble emission linewidth as narrow as the homogeneous linewidth, leading to highly ordered superlattices of colloidal QDs.^{35,36} By inducing favorable alignment of transition dipole moments of neighboring QDs, the superlattice structure could offer unique optoelectronic properties of delocalized exciton. Recent studies have also focused on enabling exciton delocalization by using electronically conductive surface ligands.^{37–41} If an exciton is delocalized over several QD sites, colloidal QD solid can achieve superradiance or superfluorescence as observed in molecular aggregates^{42} and epitaxially grown QDs.^{43} Most recently, superfluorescence has been reported in colloidal QD solids made from cesium lead halide perovskite (CsPbX $ 3$, X = Cl, Br),^{44} leading to speculations of enhanced exciton diffusion lengths in these systems through supertransfer.^{45–47}

In Secs. II and III, we have considered localized excitons, whose dynamics evolves via incoherent, hopping-type transport. In the presence of strong inter-QD interactions, however, electronic excitations can be delocalized across multiple QDs, leading to excitonic states that are superpositions of individual QD wave functions. Here, we discuss the potential implication of exciton delocalization on the overall exciton transport.

^{48}

^{49}In our model, exciton hopping from the donor (D) to the acceptor (A) eigenstate is captured by the generalized Förster theory described by the Fermi’s golden rule

^{50}

There are several ways to experimentally increase electronic couplings in QD solids, including reducing inter-QD distances, using longer-wavelength emitters, and decreasing temperature. Moreover, several reported experiments have aimed to facilitate exciton delocalization through ligand exchange by replacing standard insulating ligands with exciton-delocalizing ligands such as methylthiophenolate^{38} and phenyldithiolcarbamate.^{39,41} Quantifying the electronic coupling in these strongly coupled systems requires care because the physics that drive these couplings are generally not well approximated using simple Förster theory.^{51} Furthermore, our model does not include some of the changes that can affect exciton transport followed by ligand exchange, including the reduction in inter-QD spacing distances, higher probability of exciton dissociation in strongly coupled QDs, and increase in surface traps and exciton quenching sites. Therefore, associating the parameter $ J c$ with experimentally measurable QD systems and observing enhanced exciton transport upon electronic coupling would be experimentally challenging.

In this section, we have used the same material parameters as CdSe QD solids to model delocalized exciton transport as done for localized exciton transport in earlier sections to systematically predict the enhancement in exciton diffusivity upon exciton delocalization. In doing so, we have assumed that all material parameters have been held constant (inter-QD spacing, homogeneous/inhomogeneous linewidth, and Stokes shift) except the coupling constant $ J c$. We have found that an order-of-magnitude increase in exciton diffusivity is expected even when the average exciton delocalization size is only $\u223c4$ QDs. Therefore, as long as there is sufficiently strong electronic coupling between a few neighboring QDs, there is a possibility of significantly enhanced exciton transport—even without perfect alignment of transition dipole vectors.

## V. CONCLUSIONS AND OUTLOOK

Disorder is an intrinsic property of QD solids and can manifest as spatiotemporal variations in excitation energies, otherwise known as inhomogeneous broadening. In general, an excitation of inhomogeneously broadened ensemble of QDs leads to an energetic relaxation of exciton population, which results in a mean exciton diffusivity that is lower than expectations for a perfectly ordered material. Our work highlights two directions in which mean exciton diffusivity in disordered QD solids can be increased. First, exciton diffusivity can be maximized by specifically balancing the effects of homogeneous broadening, inhomogeneous broadening, and Stokes shift. Second, under ambient conditions, an order-of-magnitude increase in exciton diffusivity can be achieved by harnessing the effects of exciton delocalization. Notably, we have observed that the remarkable energy transport enhancements arising due to exction delocalization are robust against disorder in the energetic and orientational arrangements of QDs in the material. While there exist several experimental challenges associated with observing exciton delocalization, our theoretical study motivates continued efforts to achieve exciton delocalization by controlling QD–QD coupling strength either at the individual QD level through surface chemistry or at the collective level by fabricating highly ordered QD arrays.

## ACKNOWLEDGMENTS

E.M.Y.L. acknowledges helpful discussions with Hendrik Utzat. The submission of this work has been supported by the 2017 AIChE Annual Meeting’s Electronic and Photonic Material Graduate Student Award, sponsored by the *Journal of Vacuum Science and Technology*. This work has been funded by the Center for Excitonics, an Energy Frontier Research Center funded by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award No. DE-SC0001088 (MIT).

## REFERENCES

*Chemical Dynamics in Condensed Phases*