Control of electronic band profiles through depletion layer engineering in core–shell nanocrystals

Fermi level pinning in doped metal oxide (MO) nanocrystals (NCs) results in the formation of depletion layers, which affect their optical and electronic properties, and ultimately their application in smart optoelectronics, photocatalysis, or energy storage. For a precise control over functionality, it is important to understand and control their electronic bands at the nanoscale. Here, we show that depletion layer engineering allows designing the energetic band profiles and predicting the optoelectronic properties of MO NCs. This is achieved by shell thickness tuning of core–shell Sn:In2O3–In2O3 NCs, resulting in multiple band bending and multi-modal plasmonic response. We identify the modification of the band profiles after the light-induced accumulation of extra electrons as the main mechanism of photodoping and enhance the charge storage capability up to hundreds of electrons per NC through depletion layer engineering. Our experimental results are supported by theoretical models and are transferable to other core-multishell systems as well. Surface states, and the combination of suitable materials, induce spatial gradients in the carrier density of doped metal oxide nanocrystals, affecting their electronic structure and plasmonic behavior. Here the authors demonstrate depletion layer engineering and control in ITO/In2O3 core–shell nanocrystals by tuning the shell thickness or by photodoping.

D oped metal oxide (MO) nanocrystals (NCs) are gaining the attention of the scientific community thanks to their unique properties, such as high electron mobility 1 , the tuneability of their carrier density over several orders of magnitude 2 , chemical stability 3 , and low toxicity 3 , as well as suitable operating temperature 1 , which makes them appropriate for a large plethora of applications, ranging from nanoelectronics and plasmonics to the next-generation energy storage [3][4][5][6][7][8][9][10][11] . In doped MO NCs, surface states, such as surface trap states, defects, vacancies, as well as surface ligands and other bound molecules induce Fermi level pinning causing an upward bending of the energetic bands 2,4,[12][13][14][15][16][17][18] . The spatially varying conduction band translates into a gradient in the carrier density (n e ), sufficient to suppress entirely the metallic behavior of carriers close to the nanocrystal surface. This depletion region effectively acts as a dielectric 2,12,17,19 . Hence, the homogeneous flat-band model, which neglects Fermi level pinning, is not sufficient to accurately describe the behavior of free carriers in MO NCs, as introduced by other groups 2,17,19,20 . In fact, the depletion layer formation considerably affects the conductivity of NC films and their plasmonic behavior, with direct implications on the electric field enhancement, the localized surface plasmon resonance (LSPR) modulation and its sensitivity to the surroundings 2,17,19 . Furthermore, the presence of a surface depletion region induces an important alteration to the particle dielectric function 2 .
Given the strong impact of depletion layers on the optoelectronic properties of nanoscale oxide materials, in this work, we aimed at exploiting the depletion layer formation to control energetic band profiles as a means to understand and improve material characteristics. We explore depletion layer engineering beyond surface states by introducing additional electronic interfaces and by dynamically modulating the carrier density via post-synthetic approaches. We experimentally exemplify this scheme with Sn-doped Indium Oxide (ITO)-In 2 O 3 core-shell NCs and the fine-tuning of the shell thickness (t s ) as well as capacitive charge injection with light (i.e., photodoping). Numerical simulations on both cases serve as a framework to describe in detail the nanoscale evolution of their electronic structure supported by an empirical model that describes the experimental optical properties of all NCs before and after photodoping. The empirical fit model together with electron counting experiments support the band structure calculations well. Through this combined theoretical and experimental work, we unveil that double band bending is a key characteristic of ITO-In 2 O 3 core-shell NCs, describing well also the dynamic introduction of extra electrons via photodoping, a process not fully explained yet. We found that the photo-induced band bending results in an increase in n e predominantly in shell, contradicting the previously reported explanation of a uniform rise of the flatband Fermi level as main mechanism for photodoping 11,16,21,22 . Furthermore, the observed band bending supports charge separation towards the NC interface and avoids possible recombination. We finally exploit depletion layer engineering to improve the capacitive charging process in doped metal oxide nanocrystals upon photodoping, resulting in an accumulation of more than 600 stored charges per nanocrystal of the same size.

Results
We performed numerical simulations based on the solution of Poisson's equation 2 within the parabolic band approximation to illustrate the band structure of NCs and their depletion layer formation (extended details on the calculations are reported in the Supplementary Information) 23,24 . Here, the depletion layer is defined as the region of the NC where n e drops below 10 26 m −3 (threshold value at which we can detect plasmonic features) 19 . In Fig. 1 we show the spatially dependent profile of the conduction band as a result of the upward band bending and its effect on the depletion layer width (W) for different parameters, such as surface potentials (E S ) (Fig. 1a), different materials (Fig. 1b) and the introduction of additional electronic interfaces ( Fig. 1c and Fig. 1d). In the first case, we consider the effect of surface states on the depletion layer formation. The effect of Fermi level pinning is modeled by a fixed surface potential (E S ¼ 0; 0:5; 1; 1:5; 2 eV), from which the band bending profile is derived. An ITO/surface electronic interface is formed. The value of E S can be found experimentally and it is a peculiar parameter for each material interface. It depends on several factors, such as specific densities of trap states, presence of defects and vacancies as well as surface ligands 2 . For increasing E S , we observe an increase in the depletion layer width, which affects a larger fraction of the NC volume (Fig. 1a). These results indicate the importance of surface control to engineer the band structure of NCs. Fig. 1b reports the effects of changing the composition of the NC while keeping E S at a fixed energy. The choice of material, the elemental composition, permittivity (ε), bandgap energy (E g ) and the control over doping levels are of fundamental importance. Different materials, as in this case ZnO, In 2 O 3 , ITO and CdO, have a specific impact on the band bending, showing that W is a unique feature of each system. Another powerful parameter to control the depletion layer and the energy level profile is the introduction of additional electronic interfaces beyond the surface of the nanoparticle. One example is the ITO-In 2 O 3 core-shell nanocrystal system (Fig. 1c).
By surrounding the core ITO NC with In 2 O 3 , two electronic interfaces are formed: ITO/In 2 O 3 and In 2 O 3 /surface. In this case, E S is approximately 0.2 eV below the conduction band minimum of In 2 O 3 , as reported in literature 2 . While in uniform NCs (ITO core only) the band profile is determined by the radial depletion region near the NC surface, the addition of shell layers with thickness t s strongly affects the band's profile ultimately resulting into a double bending of the conduction band. Hence, shell tuning is an effective tool to control W and the shape of the electronic bands inside the NC volume. This effect can be further extended when combining multiple materials sequentially together in coremultishell NC architectures. Figure 1d reports a heterostructure based on three different materials, introducing three electronic interfaces (other combinations of materials and structures are reported in the Supplementary Information, Supplementary Fig. 1). This leads to non-trivial bending and highlights that it is possible to design targeted band structures at the synthesis stage by combining two or more materials in core-shell or core-multishell heterostructures with varying width.
Effective control over NC geometry, size and doping level is crucial to make reliable quantitative assessment of W. To experimentally investigate depletion layer engineering predicted by numerical calculations, we synthetized ITO-In 2 O 3 core-shell NCs with varying shell thickness t s and induced a dynamic variation of their carrier density via photodoping (see Supplementary Information for further details on synthesis methods). Figure 2a shows the TEM images illustrating the progressive growth of the NC due to the formation of the In 2 O 3 shell around the ITO core. The crystalline integrity of the samples was confirmed by XRD measurements (Supplementary Fig. 2). We collected multiple aliquots during the synthesis at different stages of the growth resulting in a set of samples with the same physical core size (R core = 5.5 nm -C0) and various shell sizes (S1-S5, with t s = 1.15 nm, 1.9 nm, 2.9 nm, 4.25 nm). The successful achievement of core-shell structures was confirmed by a comparison of the Sn-dopant concentrations obtained by two different techniques: inductively coupled plasma mass spectrometry (ICP-OES) as volume sensitive technique and X-ray Photoelectron Spectroscopy (XPS) as a surface sensitive technique (Fig. 2b). These techniques probe the volumetric and surface content of Sn atoms, respectively, and have been shown to be effective methods in elucidating nanocrystal dopant distributions 25 . We observe a higher Sn-concentration from the volume-sensitive measurements (black curve in Fig. 2b) as compared to the surfacesensitive measurements (red curve in Fig. 2b) in all samples with shell. This indicates that the Sn dopants are localized in the core of the NCs without significant diffusion of Sn atoms into the shell (further analysis on diffusion effects can be found in Supplementary Fig. 3). The absorption spectra of the representative samples, normalized to the maximum are reported in Fig. 2c (dotted curves). The spectra are governed by intense resonances in the near-infrared (NIR) that are assigned to localized surface plasmon resonances (LSPRs) as a result of free electrons in the highly doped semiconductor (typically in the range of 10 27 m −3 ) 2,4,25-27 . The LSPR peak position ω LSPR and its peak profile are correlated to several factors, such as the NC geometrical features (e.g., R core ; t s ), n e , the depletion layer width (W), the dielectric constant of the surrounding medium (ε m ), as well as the structural defects and dopant concentration, providing a unique spectral signature of such parameters 21,28-30 . We will come back to this point later when describing the empirical fit model. From the modulation of the LSPR upon shell growth, we observe an initial blue shift of the LSPR (see Supplementary Information, Supplementary Fig. 4). This is ascribed to the activation of surface dopants with the growth of a thin In 2 O 3 layer, which results in an increased carrier density 25 . The following continuous red shift of the LSPR is due to the presence of an increasing shell thickness t s that modifies the dielectric surrounding of the NC 25 . Notably, in particles with a critical thickness t Ã s = 2.7 nm, a second shoulder appears in the spectrum. This indicates a more complex carrier density profile within the core-shell nanocrystals which induces an independent resonating mode, generated by a sufficiently high carrier density in the shell of the nanoparticle 31 .
To further study the electronic structure of core-shell NCs out of the equilibrium conditions, we post-synthetically alter the number of free carriers via photodoping 4,11,21,32 . Photodoping consists of introducing multiple free charge carriers via light absorption and suppressing carrier recombination by the quenching of the holes with hole scavengers 11,21,33 . The photodoping process in colloidal NCs has been recently investigated with optical and electrochemical (e.g., potentiometric titration) measurements 4,11,21 . Here, we induce the photodoping by exposing our colloidal NCs to light beyond the ITO band-gap in the ultraviolet (UV) region (300 nm-4.1 eV, FWHM = 20 nm) and an intensity of 36.8 mW cm −2 . Figure 2c shows the normalized absorbance of three representative examples before (dotted curves) and after (solid curves) the exposure to 20 min of UV light. After the introduction of extra photocarriers into the system the LSPR absorption increases in intensity (ΔI ¼ I photodoped À I as synthesized ) and its energy shifts ). The photoinduced effects progressively appear with the amount of light Increasing the E S results in the expansion of the depletion width, W (progressively from blue to red). b Impact of different materials on W at fixed E S . c Expansion of W and double bending of the depletion layer in a core-shell structure of ITO-In 2 O 3 with a core radius (R core ) of 5.5 nm and varying shell thickness (t s ¼ 0;1;2;3;4 nm, i.e., blue, light blue, orange, red and dark red, respectively). d Multiple shell system by combining an ITO core (R core = 5.5 nm) with a In 2 O 3 and ZnO shell, with total radius R = 9.5 nm. The band shows a complex profile with a triple bending (green curve). The gray curves illustrate the previously reported case of a uniform ITO NC (dark gray) and an ITO-In 2 O 3 core-shell NC (light gray) with total radius R = 9.5 nm for comparison.
absorbed ( Supplementary Fig. 5) in agreement with previous reports in the literature 4,16,21,34 . The introduced photoelectrons add to the initial free carrier density leading to a stronger interaction with the incoming radiation and hence an increased LSPR absorption. The impact of the photodoping on the LSPR modulation is extremely sensitive to t s , with ΔI almost doubling in the case of the biggest NCs. In Fig. 2d, the normalized absorption spectra for the sample having a t s = 4.25 nm are shown before (black dotted curve) and after (black continue curve) photodoping. In this case, it is possible to note a particularly strong splitting of the LSPR into two major contributions. These results display an enhanced sensitivity of the LSPR peak to photodoping by increasing t s and indirectly hint towards an increased number of stored photoelectrons for higher t s (since the LSPR absorption is proportional to n 2=3 e ) 19 . We now investigate the same system of ITO-In 2 O 3 NCs with varying thickness t s with numerical methods as introduced above. The values for t s were chosen equivalent to the size of synthesized NCs. To further investigate also the photodoping process, we numerically calculate the effects of additional free electrons in the system as a function of t s . To this aim, we introduced a generation function G R ð Þ ¼ I 0 αβe ÀαR , which extends the Poisson's equation by an additional term that represents the spatial distribution of the extra free carriers introduced into the system via photodoping. The intensity of incident photon flux is reported as I 0 , α denotes the photon absorption coefficient, and β denotes the quantum efficiency, respectively 35 . We target to identify how their presence modifies the energy bands and carrier density distribution of the system. In this way, we go beyond the results introduced in Fig. 1 and we assess the dynamic, post-synthetic variation of the electronic band profiles via light-induced charge injection, i.e., photodoping. A comparative study reporting electronic structures and carrier density profiles both before and after photodoping are shown in Fig. 3a and Fig. 3b. We first discuss the effects of shell formation on the electronic structure of the NCs (black curves in Fig. 3a). The Fermi level pinning anchors the depleted region to the surface of the nanocrystal at the same energy, irrespective of t s . However, with increasing t s it affects more strongly the In 2 O 3 shell region, which effectively shields the ITO core from depletion. Consequently, even if W increases, the depletion layer progressively shifts towards the outer region of the NC. An intermediate region between core and surface states is, thus, introduced, resulting into the expansion of the active core region (R active ), i.e., the region of the NC volume not affected by W, which is typically larger than R core . In fact, the spatial extent of these electronic features does not correspond to the as synthesized structural parameters (i.e.,R core , t s ). This expansion is not due to an introduction of extra donor atoms, nor to diffusion effects (more details in Supplementary Information, Supplementary Fig. 3). With increasing t s , a more pronounced bending of the bands occurs, and it extends for nanometers into the NC. The corresponding carrier density distribution (black curves in Fig. 3b) shows a non-trivial profile. The double bending can be explained by a leakage of carriers into the shell. The carrier density in extended regions of the undoped shell reaches values beyond~1 Á 10 26 m −3 . The presence of carrier density in this range in the undoped In 2 O 3 region indicates that for t s >t Ã s ¼ 2:7 nm it is not appropriate to approximate the ITO-In 2 O 3 system as a doped core-dielectric shell. Instead, it must be considered as a dual-plasmonic material with a specific carrier density in the core (n core = 1.1 Á 10 27 m −3 ) and enhanced carrier density in the shell with n shell < n core . This explains the experimentally observed double features in the LSPR (see Fig. 2c), which are reproduced by our simulated absorption spectra (Fig. 2d) when implementing the carrier density profile extracted from Fig. 3b. The observed double bending of the energetic bands becomes more pronounced upon photodoping in all samples (blue curves in Fig. 3a) with an immediate impact on the carrier density distribution (Fig. 3b). Indeed, after photodoping for samples beyond t Ã s the band profile approaches a step function with two distinct energy levels: the conduction band (CB) level in the core and an energetic level approx. 0.45 eV above the CB in the shell. This effect is observed in the carrier density profile as a (close to) two-step profile. In the samples with t s > t Ã s ¼ 2:7 nm the maximum n e reached in the NC shell differs from the one in the core, reaching values of around~4 Á 10 26 m −3 , while the core carrier density remains nearly constant. The light-induced modulation of the depletion layer width (ΔW ¼ W photodoped À W as synthesized ) increases with a ΔW $ t s 3 law (Supplementary Fig. 6). Since the photo-generated extra carriers tend to fill W, larger ΔW values of the assynthesized NCs justify the possibility to store more electrons in NCs with bigger shells. Hence, from our simulations we conclude that the filling of electronically depleted regions is the main mechanism behind the photodoping process of metal oxide NCs. These findings seem to be in contradiction with the literature reports on the experimentally observed uniform rise of the Fermi energy level as a result of the photo-induced accumulation of multiple photoelectrons, as shown by Schimpf et al. 11,16,21,22 . However, the plasmonic effects of double features observed after photodoping of core-shell nanocrystals are not explainable with a simple Fermi level rise. In a flat-band scenario a uniform rise of Fermi level would necessary imply a blue-shift of ω LSPR , while we experimentally observed photodoped NCs with no blue-shift, a red shift (see Fig. 4b, below) or even a splitting of the ω LSPR To further test our theory, we approach the photodoping process by applying an empirical model to fit the spectra of each sample. Representative fits are shown as orange curves in Fig. 2d. The plasmonic properties of doped MO NCs can be well described within the framework of the Mie scattering theory in the quasistatic approximation. For our samples, we found that the classical Drude model is sufficient to accurately describe the plasmon response, with quantum effects representing only a minor correction. However, for small NCs (e.g., R < 4 nm), and low free charge regime (e.g., N e <100), quantum mechanical effects cannot be neglected 36 (further discussion in the Supplementary Information). The optical response of metals and heavily-doped semiconductors is characterized by the polarizability of the free electrons, depicted in the Drude-Lorentz model with the complex dielectric permittivity ε ω ð Þ ¼ ε 1 À ω 2 P = ω 2 þ iωΓ À Á : Here, the bulk plasma frequency ω p ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi n e e 2 =ε 0 m Ã p is a function of the free carrier density (n e ) and the effective electron mass (m Ã ), Γ is a damping parameter accounting for electron-electron scattering and ε 1 is the high frequency dielectric constant. Within this picture, ω LSPR is directly linked to ω p of the material. The tuneability of the LSPR is provided by the proportionality to n e , which is related to the number of free charges over the active volume (n e $ N e =R 3 active ). Hence, we can link the absorption, which is our physical observable, to the electronic structure of the system. In previous works, the effect of depletion layers was addressed by Zandi et al., who introduced an effective dielectric function using a Maxwell-Garnett effective medium approximation (EMA) 2 . This approach shows that accumulating charges in the NC as a result of electrochemical doping has the effect of increasing the intensity and shifting the position of the LSPR peak as a result of the varying W 2 . We adapt this model by implementing a core-shell structure with frequency dependent core dielectric function ε core ω ð Þ (and constant carrier density, n core ) surrounded by a dielectric shell with ε DL ¼ 4 in the depletion layer. Outside the sphere a dielectric medium with fixed ε m ¼ 2:09 is present. Within this picture, we approximate the continuous carrier density profile n e ðRÞ with discrete regions of uniform density, while we define n e ¼ 0 inside the depletion region (Fig. 3c). We found that the two-layer model describes well the optical spectra when t s <t Ã s ¼ 2:7 nm (Fig. 3c, i and ii). Importantly, for t s >t Ã s and most photodoping cases, we found that it was necessary to extend this model in order to fit the spectra. To this aim, we developed a threelayer model based on the Maxwell-Garnett EMA with three concentric regions. In the first two regions the inner core and first shell region, which sum up to R active , have a frequency dependent dielectric constant of ε core ω ð Þ and ε shell ω ð Þ with constant carrier densities of n core R ð Þ and n shell ðRÞ, respectively. Surrounding the frequency dependent core and shell dielectric functions is an additional layer that accounts for the depletion of carriers in the shell, which was not previously considered in models found in literature 2,28,31 . Hence, these two concentric regions are surrounded by a third depleted layer of thickness W with fixed ε DL ¼ 4 and zero carrier density. The surrounding dielectric medium is ε m ¼ 2:09. By taking into account the formation of an additional depletion layer due to the electronic interface between shell and surface, our model goes beyond what has been implemented so far to describe capacitive charges in MO NCs 2,28,31 . In our study, for all values of t s , the most notable changes in n e after photodoping are observed to effectively increase R active and decrease W 28 . The core carrier density n core remains nearly constant with variations of less than 14%, while a significant variation occurs in the shell regions, with n shell around $ 5:4 Á 10 26 m −3 .
We give a quantitative comparison between the numerical and empirical approach by plotting the amount of stored carriers in the NC (ΔN e ), defined as the difference between the free carriers of the photodoped NC and the as-synthesized NC (ΔN e ¼ N photodoped e À N as synthesized e ). We observe a good agreement between both approaches, finding that ΔN e increases with t s with a ΔN e $ t s 3 trend, reaching values as high as 600 extra electrons (Fig. 4a). We advanced the studies of the NC stored carriers by using titration on photodoped NCs to count the number of stored electrons (further details in the Supplementary Information) 4,11,37 . By using molecular oxidants (F4TCNQ) to titrate the electrons, we directly measure the average number of electrons extracted per NC. F4TCNQ in this study acts as an electron acceptor. The optical features of F4TCNQ hold as a signature to quantify the extracted electrons. We observe an increase of the number of extracted Experimental comparison between the optical response of two samples with same size and doping concentrations but different electronic structure, before (dotted line) and after (solid line) photodoping (homogeneous ITO in blue, core-shell ITO-In 2 O 3 in red). The sensitivity of the LSPR modulation via photodoping is enhanced in the core-shell case. We highlight that the peak position of the LSPR after photodoping blueshifts in the homogeneous case, while it redshifts in the core-shell case indicating that depletion layer modulation is the main process of photodoping (see discussion above). photocarriers with increasing t s , in agreement with the trend reported for numerical simulations and empirical modeling (Fig. 4a). The discrepancy, up to a factor 2 in the case of large core-shell NCs, is most probably related to a reduced efficiency in the carrier extraction process. Nevertheless, the ΔN e $ t s 3 trend is reproduced displaying that the electron counting experiments together with the empirical fit model support the band structure calculations well.
Finally, we aim at isolating the impact of depletion layer engineering from the volume dependence of ΔN e as shown from Schimpf et al. 16 . From numerical simulations, we found that the number of stored carriers in NCs with a core-shell architecture are significantly larger than in the pure ITO case ( Supplementary  Fig. 7a and Supplementary Fig. 7b). To confirm this result experimentally, we perform a quantitative analysis of one specific t s and compare it to a similar NC without shell (only core) with all other parameters unchanged (i.e., total NC radius R, doping density N d , experimental conditions). The optical absorption spectra before and after photodoping are depicted in Fig. 4b. By applying our empirical model to this case, we obtain that core-shell NCs can accumulatẽ 40% more carriers than uniform ITO NCs of the same size. Our numerical simulations predict that this enhancement increases with increasing shell thickness. Hence, we demonstrate that depletion layer engineering can improve charge storage capacity and more generally that the band bending delivers an additional degree of freedom to artificially engineer the optoelectronic properties of MO NCs, both during synthesis and post-synthetically.

Discussion
In this work, we demonstrate that depletion layer engineering is an important tool to design and control energetic band profiles in metal oxide NCs. Our results are based on a combination of theory and experiment: we implement a numerical model that is able to account for additional free carriers in the MO NCs, we developed an empirical three-layer model that describes the optical response of the (photodoped) core-shell ITO-In 2 O 3 NCs, and confirmed our results with experimental electron counting experiments through reaction with F4TCNQ. From this combined theoretical and experimental approach, we found that, first, double bending of the bands dominates the electronic structure of (photodoped) core-shell ITO-In 2 O 3 NCs and that the depletion layer predominantly affects the In 2 O 3 shell. Second, the electronic rearrangement of energy bands and the filling of electronically depleted regions resulting in the evolution of different levels of carrier density in core and shell are the main mechanism behind the photodoping process of metal oxides NCs. Third, depletion layer engineering allows enhancing the charge storage capability of ITO NCs of the same size. We can extend this model to other systems as well demonstrating the validity of our approach. Our results show that the modulation of the depletion layer represents an interesting avenue to design and improve the properties of MO NCs and their core-shell or core-multishell structures. We foresee multiple practical applications ranging from energy storage to sensing for ITO-and similar metal oxide nanocrystalsbased devices that will benefit from the control of electronic band profiles through depletion layer engineering.

Methods
Core-shell nanocrystal synthesis. ITO/In 2 O 3 core/shell nanocrystals were synthesized in a continuous growth approach with the following step-by-step procedure 19,25,38 . A precursor solution was prepared mixing in a flask tin(IV) acetate and indium(III) acetate in a 1:9 Sn:In ratio. Subsequently, oleic acid was added in a 1:6 metal to acid ratio to yield a 10% Sn doped ITO precursor solution. The flask was stirred and left at 150°C under N 2 for 3 h for degassing The ITO nanocrystals (core) were first prepared by adding the ITO precursor solution via a syringe pump (drop-by-drop at a rate of 0.35 mL/min) to 13.0 mL of oleyl alcohol at 290°C. During the slow-injection procedure a flow of 130 mL/min of N 2 gas was kept in the reaction flask to quickly remove any water vapor formed during the reaction. The ITO cores, stabilized with oleic acid ligands, were continuously grown to a size of 5.5 nm (radius) and isolated by precipitating with 12 mL ethanol. The solid part was collected by centrifugation at 7300 rpm (5540×g) for 10 min, washed twice more with ethanol and dispersed in hexane.
Then, part of the cores was kept for analysis and the rest of the solution reintroduced in fresh oleyl alcohol. For shelling, a second precursor solution was prepared by following the same procedure. In order to yield an undoped indium oleate precursor solution, indium(III) acetate was mixed with oleic acid in a 1:6 molar ratio. Undoped indium oleate was added with the same slow-injection procedure described above. Core-shell samples were then washed with ethanol and the process repeated several times until a final size of~10 nm (radius) was reached. All experiments were performed on samples collected at different stages of the shell growth, and hence sharing the very same ITO core.
Structural characterization of core-shell NCs. The structural characterization of the samples with different shell thickness were analyzed by transmission electron microscopy (TEM) to determine the size and confirm the successful formation of nanocrystals. TEM measurements were performed with a JEOL JEM-1400Plus operating at 120 kV and using lacey carbon grids supported by a copper mesh. The size distribution of the NCs was extracted using ImageJ tools on the images collected 39 .
X-ray Diffraction (XRD) analyses were carried out on a PANanalytical Empyrean X-ray diffractometer equipped with a 1.8 kW Cu Kα ceramic X-ray tube and a PIXcel3D 2×2 area detector, operating at 45 kV and 40 mA. Specimens for the XRD measurements were prepared by dropping a concentrated NCs solution onto a zero-diffraction silicon substrate. The diffraction patterns were collected under ambient conditions using a parallel beam geometry and the symmetric reflection mode. XRD data analysis was carried out using the HighScore 4.1 software from PANalytical.
X-ray Photoemission Spectroscopy (XPS) measurements were performed on a Kratos Axis Ultra DLD spectrometer, using a monochromatic Al Kα source (15 kV, 20 mA). Specimens were prepared by dropping a concentrated NCs solution onto a highly ordered pyrolytic graphite (HOPG, ZYA grade) substrate. High resolution spectra of the Sn 3d and In 3d regions were acquired at pass energy of 10 eV, and energy step of 0.1 eV, over a 300 ×700 microns area. The photoelectrons were detected at a take-off angle of ϕ = 0°with respect to the surface normal. The pressure in the analysis chamber was maintained below 7 × 10 −9 Torr for data acquisition. The data was converted to the VAMAS format and processed using the CasaXPS software, version 2.3.24 40 . The binding energy (BE) scale was internally referenced to C 1 s peak (BE for C-C = 284.8 eV). For the quantitative analysis, the areas of In 3d and Sn 3d peaks were calculated after applying the appropriate background correction across the binding energy range of the peaks of interest. The relative atomic concentrations were then estimated, using the so-called relative sensitivity factors (RSF) provided by Kratos (RSF In 3d = 7.265, RSF Sn 3d = 7.875).
Inductively coupled plasma mass spectrometry (ICP-OES) was performed on all samples to estimate the doping levels and concentrations of the ITO NCs. The elemental analysis was carried out via inductively coupled plasma optical emission spectroscopy (ICP-OES) on an iCAP 6000 Series ICP-OES spectrometer (Thermo Scientific). In a volumetric flask, each sample was dissolved in aqua regia [HCl/ HNO 3 3:1 (v/v)] and left overnight at RT to completely digest the NCs. Afterward, Milli-Q grade water (18.3 MΩ cm) was added to the sample. The solution was then filtered using a polytetrafluorethylene membrane filter with 0.45 μm pore size. All chemical analyses performed by ICP-OES were affected by a systematic error of about 5%.
Optical measurements. Optical measurements were carried out on a Cary5000 UV−vis−NIR Spectrometer. Spectra were collected in anhydrous toluene in the range 280-3200 nm with a scan resolution of 1 nm. After drying the solvent, ITO NCs were transferred in anhydrous toluene (Sigma-Aldrich) in a nitrogen filled glove box. Rectangular anaerobic cuvettes with a sealed screw cap (Starna Scientific) were used for photodoping and titration experiments.
Photodoping process. Before photodoping, the ITO-In 2 O 3 NCs were dissolved in anhydrous toluene, as described above. Subsequently, the photodoping process on the NCs is achieved by illuminating the quartz cuvette containing the solution with the NCs with a UV LED (central wavelength: 300 nm, bandwidth: 20 nm). The cuvette was placed at a distance of 12 mm from the cuvette window (Thorlabs M300L4). UV power density at the front window of the cuvette was 36.8 mW cm −2 .
Redox titration. The titrant was prepared by dissolving 0.34 mg of F4TCNQ (2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane) in 30 mL of anhydrous toluene. The NCs solution was prepared in anhydrous toluene and photodoped in the same manner as described above (typical concentration~0.1• 10 −9 mol/L). The titrant addition steps were carried out in the inert environment of a nitrogen-filled glove box to avoid any contact with ambient oxygen. Electron counting was performed after photodoping by spectroscopic analysis of the neutral, anion and dianion forms of the F4TCNQ molecules 37 . In details, the amount of F4TCNQ (n F4TCNQ , in moles) added at each step of the experiment was calculated from volume of titrant introduced in the cuvette (V), the titrant concentration (C ¼ 0.085 mg/mL), and the titrant molecular mass (276.15 g/mol): n F4TCNQ ðVÞ ¼ CV=276:15. The number of NCs present in solution (n NC ) was calculated from the mean NC size (from TEM images) and the average weight of a NC (from ICP-OES measurements). Thus, the number of F4TCNQ molecules reacted per NC was calculated as: n reacted ðVÞ ¼ n F4TCNQ ðVÞ=n NC . To calculate the number of extracted electrons, n reacted was then multiplied to a factor one or two according to the kind of reaction involved, corresponding to the formation of dianion or anion species, respectively. Two volumes were identified: V 1 , corresponding to the saturation of the dianion reactions (with exclusively two-electron transfers occurring), and V 2 , corresponding to the appearance of neutral peaks signatures (signaling the presence of non-reacted titrants). The growth of anion peaks (i.e., transfer of one electron) between V 1 and V 2 indicate that electron transfer reactions keep occurring after V 1 . In this study, two-electron transfers were considered up to the midpoint between V 1 and V 2 (V mid ¼ ðV 1 þ V 2 Þ=2 and the transfer of one extra electron was considered between V mid and V 2 . The total number of extracted charges was estimated as 2e À Á n reacted V mid À Á þ 1e À Á n reacted V 2 À Á À n reacted V mid À Á À Á . Error bars are representative to the distance between V 1 and V 2 . The effects of titrants on assynthesized ITO-In 2 O 3 samples were tested, showing no sign of interaction in the spectrum. Multi-layers fitting model for LSPRs. The distinct dielectric response of the core-shell NCs is implemented as an effective dielectric function ε eff ω ð Þ based on the Maxwell-Garnett effective medium approximation (EMA). This model is further extended to consider multiple shell regions and corresponding dielectric environments. We fit the experimental data with a particles warm optimization algorithm in MATLAB (R2020a. Natick, Massachusetts: The MathWorks Inc.) and we extract the carrier densities n e;core and n e;shell and spatial extensions (R core , R shell ) of the core and shell regions, respectively, for each NC of increasing shell thickness before and after photodoping.
COMSOL simulations. Simulations for the energy band and carrier density profiles were solved numerically for spherical NCs using a finite-element method. Poisson's equation was solved with the software COMSOL Multiphysics v5.6 (Comsol Inc., Burlinghton MA USA) using a finite-element scheme (see Supplementary Information for details).