The CdTe QDs employed as reconfigurable elements, were obtained by a hydrothermal synthesis method. In a typical synthesis, 0.0533 g (0.2 mmol) of Cadmium acetate dihydrate (Cd(CH₃COO)₂ 2H₂O, 99.5%) were dissolved in 50 ml of deionized (DI) water in a 125-ml Erlenmeyer flask. Next, 18 µl of thioglycolic acid (TGA) were added and 1 M NaOH solution was incorporated dropwise to adjust the pH to a value between 10.5 - 11. In a separate flask, 0.0101 g (0.04 mmol) of potassium tellurite (K₂TeO₃, 95%) were dissolved in 50 ml of DI water and added to the first solution. Subsequently, 0.08 g of sodium borohydride (NaBH₄, 99.99%) was added to the mixture and the reaction was allowed to stir for 5 min. The solution was transferred to a three-neck round flask and attached to a Liebig condenser and refluxed under atmospheric conditions, at a temperature of 100˚C. QD size is determined by the reflux time. The solution was then washed with acetone and redispersed in deionized water to obtain a clean colloidal solution. A high-resolution TEM micrograph of the synthesized QDs is shown in the inset of Figure 1(b). Since the QDs are capped with TGA, their surfaces contain negatively-charged carboxylate groups [30] [31]. This makes CdTe QDs a suitable candidate for a layer-by-layer deposition employing positively and negatively charged polyelectrolytes, namely, Poly(diallyldimethylammonium chloride) and Poly(styrene sulfonate), respectively.

Quantum dot optical absorption characterization was performed employing an Ocean Optics Flame-S-UV-VIS spectrometer (Figure 1(a)). Bandgap energy and absorption are related by the following expression [32]

α h υ = A ( h ν − E g ) ( 1 / r )

where α is the experimentally collected absorption, h is Planck’s constant, υ is the photon’s frequency, E_g is the quantum dot bandgap energy and A is a proportionality constant. The nominal value of r reflects the nature of the optical transition. For a direct allowed transition, r is set to be 2. The optical bandgap can then be estimated graphically employing a method described by Tauc [33] (Figure 1(b)). The extracted bandgap value of CdTe QDs was estimated to be 3.11 eV. Alternatively, after a mathematical manipulation, the above equation can be rewritten as [34]

d ( ln ( α h ν ) ) d ( h ν ) = 1 r ( h ν − E g )

According to the above equation, a plot of d ( ln ( α h ν ) ) / d ( h ν ) vs. h ν should exhibit a discontinuity when h ν − E g → 0 . However, due to impurities and QD size inhomogeneity, the discontinuity could be replaced by a maximum. Figure 2 illustrates this method, from which a bandgap value of 3.18 eV was obtained, in reasonable agreement with the previously estimated value employing Tauc’s graphic method.

Additionally, QD size can be estimated using Brus’ quantum mechanical approximation [35]

E Q D = E B u l k + ℏ 2 π 2 2 r 2 ( 1 m e * + 1 m h * ) − 1.8 e 2 4 π ϵ ϵ 0 r

where r is the radius, ϵ = 10.2 is the relative permittivity of bulk CdTe, m e * = 0.11 m 0 and m h * = 0.4 m 0 are the respective electron and hole effective

masses, e is the electron charge, and E Q D and E B u l k are the QD and bulk bandgap values, respectively. Utilizing these parameters, the radius of the synthesized QDs was estimated to be 1.56 nm. Moreover, CdTe QDs size can also be approximated by the empirical equation [30]

D = ( 9.8127 × 10 − 7 ) λ 3 − ( 1.7147 × 10 − 3 ) λ 2 + ( 1.0064 ) λ − 194.84

where D is the diameter of the QDs and λ is the wavelength of the first excitonic absorption peak. As shown in Figure 1(a), the first excitonic absorption peak is located at a wavelength of 540 nm, indicating a diameter of 3.2 nm, in close agreement to the quantum mechanical approach. The confinement state of the generated electron-hole pairs can be determined by comparing the previous value to the Bohr radius corresponding to the first excitonic state which is given by the equation [36]

r B = ϵ ϵ 0 h 2 π e 2 ( 1 m e * + 1 m h * )

Since the calculated Bohr radius for the first excitonic state is 7.5 nm for CdTe, and, therefore, larger than the size of the synthesized QDs, it evinces the strong confinement state of electrons and holes [37].

In order to conduct THz measurements employing colloidal quantum dots, water must be exchanged as the host solvent, due to its extremely high absorption coefficient [38]. A thiol exchange scheme was employed to make the QDs solvable on toluene, which has a much lower THZ absorption coefficient [39], [40]. To this end, 1 ml of 1-dodecanethiol was added to 1 ml of a clean solution of CdTe QDs. Subsequently, 2 ml of acetone were added to the mixture followed by 2 minutes of vortexing. The solution was then heated to 60˚C until the organic phase was reduced to 1 ml. After the QDs had been transferred to the organic phase, they were diluted with toluene and precipitated with methanol. The product was decanted and diluted once more with toluene.

Quantum dot concentration was measured employing a quartz crystal microbalance (QCM) instrument. The technique is based on the resonance frequency shift of a piezoelectric crystal. The deposition of a liquid drop of colloidal QDs on top of an oscillating quartz crystal provokes a shift in its resonance frequency, which is related to the mass of the residual QDs film after the solvent was evaporated by the following equation [41]

Δ m = Δ f C f

where Δ f is the frequency shift and C f is the crystal sensitivity factor (56.6 Hz∙μg⁻¹∙cm²). The concentration of CdTe QDs dissolved in toluene was measured by dropping 15 μl of solution on the crystal’s active area. After the solvent was evaporated, the residual mass was measured. This process was repeated 9 times and the results were linearly fitted to obtain the concentration in units of μg/μl (see Figure 3).

Quantum Dot Terahertz characterization was carried out by employing a THz time domain spectroscopy (THz-TDS) technique, operating in transmission mode. The system makes use of a photoconductive antenna to generate an ultrafast, broadband THz pulse when is activated by a femtosecond laser with a wavelength of 800 nm. A second photoconductive antenna is utilized as a detector. Once proper alignment was achieved, the system was thermalized for 2 hours. An ultraviolet lamp (365 nm, 20 mW/cm²) was employed to excite the sample.

The time domain Terahertz signal is transformed into the frequency domain by a Fast Fourier Transform to obtain the amplitude and phase change. Air and the empty cuvette were used as references to calculate the complex permittivity

of toluene and colloidal quantum dots [42]. Subsequently, linear effective medium theory was employed to obtain the permittivity of the Quantum dots through the relation [43]

ϵ m i x = c ϵ Q D s + ( 1 − c ) ϵ h

where c is the volume fraction of colloidal QDs and ϵ h is the permittivity of the host solvent, which, for CdTe QDs was toluene. Figure 4 shows the real and imaginary part of the permittivity of the QDs, after applying the aforementioned effective medium theory. Ostensibly, the real and imaginary parts of ϵ increase when the nanoparticles are excited by an ultraviolet source. This change is thought to be related to the large exciton polarizability of the QDs, due to their atom-like response, which is present in the strong confinement state [44].

Finite Difference Time-Domain (FDTD) simulations were performed in order to elucidate the electromagnetic response of the metamaterial geometry. This technique solves Maxwell’s equations numerically in discretized time and space and the results are the electric and magnetic field distributions following the interaction of the source fields with the material distribution. The simulations were carried out employing a commercial solver (Lumerical Inc., ND). The dimensions of the simulated metamaterial unit cell are shown in Figure 5(a) employing a design reported elsewhere [45]. The outer black hexagon represents the boundary of the unit cell and the red lines represent a Perfect Electric Conductor (PEC).

The metallic incursions are placed on top of a silicon slab with a refractive index of n = 3.42 + 0.01i [46]. Periodical boundary conditions were set in the directions orthogonal to the electromagnetic source, and a Perfectly Matched Layer (PML) was set in the direction parallel to it. The electromagnetic source

was a plane wave with the electric field aligned orthogonal to the capacitor plates in the center of the unit cell. An optical image of the microfabricated metamaterial structure is shown in Figure 5(b).

Figure 6 shows the simulated and experimental transmission graphs featuring a resonance at 0.62 THz. The electric field distribution at the resonant frequency is illustrated in Figure 6(b). Ostensibly, the electric field is strongly confined within the capacitor area of the unit cell.

For Metamaterial fabrication a standard image reversal lift-off procedure was performed on a <100> high resistivity silicon substrate, with a nominal thickness of 400 μm. Silver was chosen as the material for the metallic array due to its relatively low-loss properties [47]. To this end, thermal evaporation was employed to deposit a silver layer with a thickness of 100 nm on top of a 20 nm layer of

chromium that was included for adhesion purposes. Terahertz Time-Domain Spectroscopy (THz TDS) in transmission mode was employed to experimentally explore the resonant properties of the metamaterial. The measurements were performed at room temperature with a Menlo Systems TeraSmart K15 equipment. Five thousand waveforms were averaged for each sample, and the time domain waveforms were passed through a Gaussian window to remove the Fabry-Perot oscillations in order to obtain a cleaner signal [48]. Subsequently, a Fast Fourier Transform (FFT) was applied to obtain the amplitude and phase of the transmitted beam as a function of frequency. In order to extract the effect of the metamaterial, we define the transmission as the ratio of the amplitude of the electric field through the metamaterial structure to the one through a blank, 400 um thick silicon wafer. Prior to deploying the synthesized QDs, the substrate with the metamaterial geometry was exposed to an O₂ plasma (30 W for 10 minutes) to activate the surface, and two bilayers of alternating positively charged Poly(diallyldimethylammonium chloride) (PDDA) (1.5 wt%) and negatively charged Poly(styrene sulfonate) (PSS) (0.3 wt%) were spin coated at 2500 rpm for 2 min. Subsequently, CdTe QDs were then deployed on the metamaterial by employing a bilayer scheme (Figure 7) [49] in which a layer of PDDA was spin cast on the surface, followed by a layer of the aqueous solution of negatively charged CdTe QDs. This last alternating bilayer step was repeated 6, 10, 13 and 15 times to evaluate the performance as a function of the number of bilayers employed. The sample surface was rinsed with DI water and dried with a gentle stream of N₂ after each individual spin coated layer. Profilometry measurements indicated that, on average, for every pair of PDDA + CdTe QDs layers, the thickness increased by approximately 10 nm.

Transmission modulation (TM) was measured employing an ultraviolet source with a wavelength of 365 nm, which corresponds to an energy of 3.4 eV, i.e., larger than the bandgap energy of CdTe QDs. The UV source was set to illuminate all samples with a different number of bilayers while varying the incident power density from 6 to 20 mW/cm². At the resonant frequency of 0.62 THz, the TM was observed to increase with all incident power density increases (Figure 8(a)). This is thought to be due to the susceptibility of the photoexcited QDs dependence on the excitation density. That is, the number of incident photons per unit area has an impact on the polarizability of the excitons confined in the QDs [44].

In terms of the number of bilayers employed, the TM value increased as the number of bilayers increased from 6 to 10 and decreased for subsequent increases in the number of bilayers, i.e., 13 and 15 bilayers. Specifically, for 6 bilayers, the modulation increased from 0.2% to 1.3%, the modulation observed

with 10 bilayers increased from 1.84% to 2.68%, and with 13 bilayers, the modulation increased from 1.2% to 1.6%. Lastly, for 15 bilayers, the modulation increased from 0.13% to 0.48%, being the smallest modulation attained. Quenching and shadowing effects are thought to be responsible for the observed decrease in transmission modulation depth when employing 13 or more bilayers which is presented in Figure 8(b).