# Hybrid thermal Yagi-Uda nanoantennas for directional and narrow band long-wavelength IR radiation sources

Marco Centini, Maria Cristina Larciprete, Roberto Li Voti, Mario Bertolotti, Concita Sibilia, Mauro Antezza

► To cite this version:

Marco Centini, Maria Cristina Larciprete, Roberto Li Voti, Mario Bertolotti, Concita Sibilia, et al. Hybrid thermal Yagi-Uda nanoantennas for directional and narrow band long-wavelength IR radiation sources. Optics Express, 2020, 28 (13), pp.19334. 10.1364/OE.389837. hal-02871572

HAL Id: hal-02871572

https://hal.science/hal-02871572

Submitted on 14 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Hybrid thermal Yagi-Uda nanoantennas for directional and narrow band long-wavelength IR radiation sources

MARCO CENTINI,<fnref n="1" />,<fnref n="2" />,* MARIA CRISTINA LARCIPRETE,<fnref n="1" /> ID
ROBERTO LI VOTI,<fnref n="1" /> MARIO BERTOLOTTI,<fnref n="1" /> CONCITA SIBILIA,<fnref n="1" /> AND
MAURO ANTEZZA<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID

**Abstract:** We investigate the possibility of spatially and spectrally controlling the thermal infrared emission by exploitation of the Yagi-Uda antenna design. Hybrid antennas composed of both SiC and Au rods are considered and the contributions of emission from all the elements, at a given equilibrium temperature, are taken into account. We show that the detrimental effect due to thermal emission from the not ideal parasitic elements drastically affect the performances of conventional thermal Au antennas in the 12 µm wavelength range. Nevertheless, our results show that the hybrid approach allows the development of efficient narrow-band and high directivity sources. The possibility of exploiting the Yagi-Uda design both in transmission and in reception modes, may open the way to the realization of miniaturized, efficient, robust and cheap sensor devices for mass-market applications.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

The availability of low cost, integrated, radiation sources in the infrared range with narrow-band emission spectrum and good directional selectivity is of great importance in a variety of applications such as infrared sensing [1], thermophotovoltaics [2], radiation cooling [3], and thermal circuits [4]. With particular reference to gas sensing applications, narrow band emission sources with a good level of spatial coherence are required in order to optimize the signal to noise ratio and to increase selectivity and sensitivity. Quantum cascade lasers (QCL) [5] are probably the best choice for mid to far IR spectroscopy due to their high spectral density and single-mode operation. Using a set of QCL tuned at different wavelengths it could be possible to match individual absorption lines of different chemical species, for example. Nevertheless, the practical application of these sources for commercial sensing devices for home, automotive markets as well as smartphone/laptop accessories are strongly limited by high fabrication costs.

A common and inexpensive way to obtain infrared radiation is to take advantage of thermal emission from a heated body. However, the spectral and directional control of thermally emitted light is a challenging task due to its incoherent behavior (both spatially and temporally). In the last decade several approaches have been proposed to control the mechanisms of thermal radiation and heat exchange on the nanoscale [6–8] also in 2D and multilayered materials [9–14]. Enhancement and control of thermal radiation have been investigated in gratings [15], metamaterials [16–19], metasurfaces [20] and single nanoantennas [21]. It has been shown, both theoretically and experimentally, that sub-wavelength patterning of materials supporting surface polaritons results in narrow emission lobes. Enabling the exploitation of near field interactions

<footnote n="1">¹Sapienza University of Rome Department of Basic and Applied Sciences for Engineering, Via A. Scarpa 14, I-00161 Rome, Italy</footnote>

<footnote n="2">²Laboratoire Charles Coulomb (L2C), UMR 5221 CNRS-Université de Montpellier, F-34095 Montpellier, France</footnote>

<footnote n="3">³Institut Universitaire de France, 1 rue Descartes, F-75231 Paris, France</footnote>

<footnote>*marco.centini@uniroma1.it</footnote>
between different objects and the resonant coupling of evanescent modes into propagative ones, it is possible to obtain efficient and partially coherent emission. The mechanism at the base of these tailored emission is the excitation of resonant surface polaritons. Metals have been extensively studied in the visible and near IR range; enhancement and selectivity of thermal emission have been obtained thanks to resonant plasmon polariton excitation. Typically, these effects are stronger in a spectral range close to the plasma frequency of the metal and become weaker in the mid-wavelength IR (MWIR). For efficient sources in the long-wavelength IR (LWIR, 8-15 µm) it is possible to rely on polar materials supporting phonon polaritons. Thermal emission from silicon carbide (SiC) gratings or silica (SiO2) patterned structures has been intensively studied [22,23]. Moreover, the emission properties of a single SiC rod antenna have been studied by numerical calculations and experiments [24]. Recently it has also been shown that it is possible to use porous SiC to create narrow-band emitting/absorbing metamaterials [25] with a certain degree of tunability in order to match absorption lines of hazardous materials [26,27] in the 12-14 µm wavelength range.

For sensing applications, the growing requirement for device miniaturization and the demand for precisely controllable and tunable sources and detectors made the nanoantenna based approaches very appealing. Antenna-coupled infrared bolometers have been intensively studied to tailor and enhance the detector's wavelength, directional, and polarization dependence [28,29]. Concerning the realization of sources in the optical and near infrared ranges remarkable results have been obtained by using gold Yagi-Uda (YU) nanoantennas coupled with quantum dots emitters [30]. The YU design (sketched in Fig. 1) is widely used for radiofrequency waves. The efficiency of the device is related to the coherence of the signal emitted by the feeder and by the low losses induced by the parasitic elements. For example, the performance of the optical YU antenna presented in Ref. [30] relies on the narrow spectral line emission from the quantum dot that resonantly couples to the feeder and to the fact that all the other elements are passive.

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Typical scheme of a 5 elements Yagi-Uda antenna.</figcaption></figure>

Recently some possible scheme of metallic antennas and composite antennas for mid-IR range have been proposed [31] and experimentally investigated [32,33]. Concerning metallic antennas, it must be noted that the field enhancement and spectral selectivity of a resonant antenna drastically decreases at longer wavelengths. In Ref. [34] a gold YU design has been investigated for emission in the wavelength range around 7 µm with a quite large linewidth of Δλ=1 µm.
In this manuscript we propose a Hybrid thermal YU antenna composed of gold and SiC, operating in the LWIR range. In particular, we consider a SiC dipole antenna as the feeder for its narrow band high emissivity and gold rods acting as parasitic elements. With respect to previous work our numerical results take into account the emission contributions from the reflector and the directors that are, reasonably, at the same temperature of the feeder. We will show that a gold thermal YU antenna has a limited efficiency in LWIR range (around $\lambda=12\mu m$) and its behavior is mostly driven by the emission of the directors rather than the feeders. In order to overcame this limitation, we propose the hybrid configuration combining metallic and dielectric elements. We show that this is a reliable solution to improve the performances of YU antenna in the LWIR. We show that it is possible to obtain narrow band sources (with a $\Delta\lambda$ of 100 nm with respect to a peak of emission at 12 $\mu m$). We also show that the tolerances are extremely robust and the same antenna design can be used to tune the emission peak of the antenna in the reststrahlen band of the SiC by only changing the length of the feeder rod. These features make thermal YU hybrid antennas very appealing for use as multiple arrays of narrow band high directivity sources for different species of gas sensing.

## 2. Materials and methods

As a starting point we compare the thermal emission of gold and SiC dipole antennas to put into evidence the difference in the bandwidth of emission spectra.

The emitted power spectral density is evaluated by applying a previously developed model based on the fluctuational electrodynamics approach and on the discretization of the resulting volume integral equation to calculate emissivity and far field spatial emission pattern of nanoparticle ensembles [35]. More in details, the numerical model is based on the fluctuation dissipation theorem [36]: for a medium with relative permittivity $\epsilon_r$ at a given temperature T the cross correlation spectral density for the chaotic current density related to thermal fluctuations is given by the expression:

$$ \langle \vec{J}_{\omega}^{0}(\vec{r}_1) \vec{J}_{\omega}^{0\dagger}(\vec{r}_2) \rangle = \frac{\omega}{\pi} \epsilon_0 \mathrm{Im}(\epsilon_{r,\omega}) \Theta_{\omega}(T) \delta(\vec{r}_1 - \vec{r}_2) \delta(\omega - \omega') \bar{I}; \quad (1) $$

Where the brackets represent a statistical ensemble average, $\delta$ is the Dirac function, $\omega$ is the angular frequency, $\bar{I}$ is the identity operator and:

$$ \Theta_{\omega}(T) = \frac{\hbar\omega}{\exp(\hbar\omega/k_B T) - 1}; \quad (2) $$

being $k_B$ the Boltzmann constant and $\hbar$ the reduced Planck constant. The symbol $\dagger$ stands for the operation of Hermitian conjugation. In a deterministic problem, the emission and scattering process form extended sources can be studied by considering a generalized current density:

$$ \vec{j}_{\omega}(\vec{r}) = -i\omega\epsilon_0(\epsilon_r(\vec{r}) - 1)\vec{E}_{\omega}(\vec{r}) + \vec{j}_{\omega}^{0}(\vec{r}); \quad (3) $$

Where the first term on the right side is related to the response of the medium (scattering) and the second term takes into account for impressed sources.

In Ref. [35] we have shown that the cross spectral density defined as:

$$ \overline{\overline{W}}_{\omega}(\vec{r}_1, \vec{r}_2)\delta(\omega - \omega') = \langle \vec{E}_{\omega}(\vec{r}_1) \vec{E}_{\omega'}^{\dagger}(\vec{r}_2) \rangle; \quad (4) $$

can be numerically evaluated on a discretized mesh as a product between an auxiliary deterministic matrix, $\overline{K}$ taking into account for the scattering problem, and the cross correlation function of the
fluctuating sources $\langle \bar{B}\bar{B}^\dagger \rangle$ given by the fluctuation dissipation theorem (see eq.1):

$$ \overline{\overline{W}}_{\omega} = \overline{\overline{K}}^{-1} \langle \bar{B}\bar{B}^{\dagger} \rangle \overline{\overline{K}}^{-1\dagger}; \qquad (5) $$

Where $\langle \bar{B}\bar{B}^\dagger \rangle$ is the numerical evaluation, on a given discrete mesh describing the system, of the following expression:

$$ \langle \bar{B}_j \bar{B}_k^\dagger \rangle = \mu_0^2 \omega^2 \int \overline{\overline{G}}_\omega(\vec{r}_j, \vec{r}') \overline{\overline{\overline{G}}}_\omega^\dagger(\vec{r}_k, \vec{r}'') \overline{\overline{\overline{J}}}_\omega^0(\vec{r}') \overline{\overline{\overline{J}}}_\omega^0(\vec{r}'') d\vec{r}' d\vec{r}''; \qquad (6) $$

where $j,k=1, \dots, N$; are the points of the mesh. All the details about the numerical integration are described in Ref. [35].

In the far-field regime, the diagonal elements of the cross spectral density (Eq. (5)) are proportional to the mean radiated intensity spectral density. Thus, the average power spectral density emitted by the source and collected at a given distance $d$ from the center, under a solid angle $\Omega_0$ can be obtained as:

$$ P_{\Omega_0}(\omega) = \frac{1}{2} c \epsilon_0 d^2 \int_{0}^{\Omega_0} \overline{\overline{W}}_{\omega}(\vec{r}_1, \vec{r}_1) d\Omega; \qquad (7) $$

where $\vec{r}_1$ are points on the sphere of radius $d$ and $c$ is the speed of light in vacuum.

We consider a SiC rod, 1.62 µm long, with a square cross section of 325X325 nm² and a Au rod, 3.44 µm long, with the same cross section. The refractive indices for both Au and SiC have been taken from Ref. [37]. The difference in length between the two rods is due to the different effective wavelength for the two materials, corresponding to a $\lambda/2$ dipole antenna in vacuum with the desired peak of emission around $\lambda=12$ µm. For simplicity all calculations are performed in vacuum. The adopted method allows to consider the system embedded in a homogenous medium provided it is transparent in the wavelength range under examination. According to the experimental results reported in Ref. [33] CaF₂ could be a good candidate to host our composite antenna for emission in the 12 µm range. Optimization of the antenna and prototype design is out of the scope of this work. Here we focus the attention on the peculiar features of the hybrid design with respect to pure metallic antennas. Once the material and the cross section are selected, longer rods will emit at longer wavelengths. We note that in this wavelength range the enhanced emission for the Au rod is only related to the fulfillment of half dipole condition because field's localization effects arising from excitation of localized surface plasmon polaritons are very weak (being far from Au plasma frequency). On the other hand, the enhanced signal and narrow spectral line from the SiC emitter is related to the resonant excitation of surface phonon polariton in its reststrahlen band. For both cases we calculated the far field power spectral density collected at $d=1$ mm distance from the antenna over a cone of maximal half-angle $\theta=15^\circ$ (numerical aperture $NA = \sin(\theta) = 0.26$) in forward (left) and backward (right) emission with respect to the antenna main axis (x-axis, as sketched in Fig. 1) when the antennas are at an equilibrium temperature of 400 K. The calculation has been performed by using cubic volume cells of size $\Delta V=65 \times 65 \times 65$ nm³. In order to check reliability and convergence of our results we performed a calculation on the single rod with a finer discretization. Increasing the number of cells by a factor 8 ($\Delta V=32.5 \times 32.5 \times 32.5$ nm³) the resonance peak shifts of 0.5% and the relative error on the emitted power is less than 5%, showing that further discretization is not required.

In Fig. 2 we report the normalized power spectral density for both cases in forward emission. Obviously, due to the symmetry of the geometry, the signal collected in forward and backward direction in the case of a single rod is the same. In order to provide a reference, all the quantities have been normalized with respect to the maximum value obtained for the SiC rod. The peak of emission for gold is about five times weaker than for SiC. We also note that there is a single
line at 12.02 µm and its full width at half maximum (FWHM) is approximately 100 nm for the SiC antenna. On the other hand, the emission for the gold rod at the same temperature has a much broader spectrum with a FWHM of approximately 2 µm. For both cases the radiative pattern evaluated at the peak of emission has the typical dipole pattern shape. For illustration purposes we report in the inset of Fig. 2 the radiation pattern calculated for the SiC antenna at the resonant wavelength of 12.02 µm. The z-axis is the axis of the dipole antenna. For the gold rod the radiation pattern is qualitatively the same. According to our calculations the estimated power emitted by a single SiC rod at 400 K in the forward direction within a cone corresponding to NA=0.26 is approximatively 10 pW.

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Normalized forward emitted spectral power density collected under a NA=0.26 calculated for the SiC dipole antenna (red) and gold antenna (blue) described in the text. The inset shows the emission pattern at the peak wavelength for the SiC antenna. The emission pattern for the gold antenna is qualitatively the same.</figcaption></figure>

As mentioned earlier, according to standard YU specifications all the director elements are half-wavelength antennas slightly detuned towards shorter wavelengths with respect to the feeder. Similarly, the reflector is detuned to longer wavelengths. We choose gold as material for the rods that will behave as directors and reflector. For simplicity and with the main objective to keep a further realization procedure as easy as possible, we kept the same section for all the elements (325 × 325 nm²) and selected the resonance wavelength by changing the length of the rods. Following the standard YU criteria [35,36], we started with the ratio between the length of the reflector (director) and the gold feeder equal to 1.15 (0.9). A fine optimization could be performed by changing either the length and the section of these elements but it is beyond the scope of this work. Moreover, our calculations show that small deviations from these values do not significantly affect the final result.

In order to quantify the amount of the background noise introduced by the emission of the not ideal parasitic elements, we firstly evaluated the thermal emissivity of just one reflector element at a temperature T=400 K. We depict in Fig. 3 the normalized power spectral density emitted by the reflector (dashed line) collected within a cone corresponding to the same NA=0.26
and normalized with respect to the spectral power emitted by the single SiC rod. We note that the reflector does not contribute much to the overall emission in the wavelength range under investigation. On the other hand, performing the same calculations on a set of one, two and three Au directors, separated by a typical 0.3λ distance (according to standard YU specifications [38,39]) we note a very different behavior. Indeed, the single director is by all means a gold antenna tuned at 10.5 µm and its emission is comparable to the gold feeder. Moreover, a system composed by two or more directors is a composite antenna, all the elements are strongly coupled, being separated by a subwavelength distance. Increasing the number of directors, resonant modes at longer wavelengths are created due to the increased size of the system and the resulting amount of “noise” is even higher than the emission level of the SiC feeder with only 3 directors. Since the emission from each element is partially coherent we note that there is an enhancement of emission directivity and a line narrowing. However the coupling between elements is weak due to poor field localization effects in this wavelength range and the line width is much larger compared to the spectral emission from the SiC element. Indeed, we will show that the SiC feeder narrow emission linewidth is crucial for the boost effect of the total signal from the YU antenna.

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Normalized emitted spectral power density for a single reflector (dashed line), one director (solid blue line), 2-(solid black) and 3-(solid red) directors at a given temperature of 400 K.</figcaption></figure>

**3. Results**

We considered the thermal emission from a 5 elements (reflector, SiC feeder and 3 directors) YU antenna, at an equilibrium temperature of 400 K. We performed a scan to find the optimal distances of reflector and directors, starting from the standard specifications. We considered only two scanning parameters: *d1* being the distance between the feeder and the reflector and *d2* the distance between the feeder and the director, as well as the distance between adjacent directors. In order to find optimized performances, we defined a function of merit taking into account two main objectives: 1) the ratio between the forward (*P*<sub>Ω₀</sub>*(forward*))and backward (*P*<sub>Ω₀</sub>*(backward*)) emitted power under the same numerical aperture NA=0.26 (function of merit
F1); 2) the previously investigated gain enhancement with respect to the single SiC rod at the peak wavelength of 12.02 µm (function of merit F2). Thus we want to optimize the function:

$$F = F1 * F2 = \frac{P_{\Omega_0}(\text{forward})}{P_{\Omega_0}(\text{backward})} * \frac{P_{\Omega_0}(\text{forward})}{P_{\text{SiC},\Omega_0}(\text{forward})}; \quad (8)$$

A map of F as a function of d1 and d2 is depicted in Fig. 4. We found optimized results by choosing 3.8 µm (0.31λ) as the distance between feeder and director and between each director, while the optimized distance between the feeder and the reflector is 2.4 µm (0.2 λ). With these values the function of merit reaches the maximum of about F=60. We note that performances of the antenna are not strongly affected by small deviations from optimal parameters. Indeed, we obtain already F>50 in a wide zone for the parameters d1 ranging from 1.9-3.1 µm and d2 from 3.5-4.1 µm.

<figure><img src="image_6.png" /><figcaption>Fig. 4. Map of the function of merit F defined in Eq. (7) as a function of the scanning parameters d1 (distance between the feeder and the reflector) and d2 (distance between the feeder and the director, distance between adjacent directors).</figcaption></figure>

This feature is extremely interesting in terms of fabrication tolerances of the antenna. In Fig. 5 we report the normalized power spectral density collected with a NA of 0.26 in forward emission for both cases of SiC and gold feeder. In the first case we note that despite of the noise introduced by parasitic elements the signal of the feeder is boosted by a factor of 5 and the narrow linewidth (approx. 100 nm) is preserved while the background emission from the directors is partially suppressed.

On the other hand, if the gold feeder is considered, we note that the signal emitted by the feeder is not boosted, the thermal emission is ruled by the directors and the overall effect is an enhancement of the emission of the directors. There is a well-defined emission line of approximately 1 µm bandwidth tuned at 11.45 µm over the background noise. These results suggest that the study of a gold thermal antenna must be performed by considering it as an array of coupled emitters/scatterers. The emission peak cannot be evaluated a priori by selecting the resonance of the feeder. In order to put into evidence the different behavior for the case of gold and hybrid antenna we calculated the signal emitted by a system composed only of reflector and
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Normalized forward emitted spectral power for a 5-element gold (black) and hybrid gold-SiC feeder (red) Yagi-Uda thermal antenna. The blue line shows the normalized forward emitted spectral power obtained from the parasitic elements only, without any feeder.</figcaption></figure>

directors, without feeder. The result is shown in Fig. 5. It is clear that the reflector is too far from the directors to act efficiently and it actually reduces the emission with respect to the system composed by 3 directors only (Fig. 3). We also note that the introduction of the SiC feeder does not modify the contribution from the other elements because it is acting on a very narrow band and it does not efficiently couple with the others. On the other hand, when we add the gold feeder, it behaves as a reflector for the directors, boosting the signal and shifting the wavelength peak due to coupling with the other elements. A finer optimization could be performed for the pure metal antenna nevertheless the large bandwidth related to gold antennas remains the main drawback. Finally, we plot the corresponding radiative pattern for both antennas at the peak of emission in Fig. 6(a-b). We note that in both cases we find the typical directional Yagi Uda emission pattern although the hybrid antenna is more efficient. The radiation patterns for both cases could be object of further optimization. However, this task goes beyond the scope of this work.

A deeper insight is provided by calculating the out of equilibrium case in which only the feeder or the directors or the reflector are at a given temperature of 400 K while the other elements are not emitting (corresponding to an ideal temperature of 0 K). In this way it is possible to see how each element contributes to the emission in the presence of the others acting only as scatterers by taking advantage of the superposition principle. The results are reported in Fig. 7 for both metallic and hybrid antennas. In Fig. 7(a) it is depicted the result for the gold YU antenna. As expected, the contribution from the directors plays the dominant role. However, the emission from the feeder is responsible for a shift of the wavelength peak. We also note that, with respect to the emission of the single elements calculated in Fig. 3 and Fig. 2 the YU scheme is responsible for the enhancement of emission of every part of the antenna (reflector, feeder and directors). Vice versa, with the hybrid YU (Fig. 7(b)), only the emission of the feeder is enhanced giving as a result a higher signal to noise ratio.
<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Normalized radiation pattern at the peak wavelength for the hybrid gold with SiC feeder Yagi-Uda antenna (a) and for the gold with gold feeder antenna (b). The quantities are normalized with respect to the maximum emission of a single SiC dipole antenna of the same size of the feeder.</figcaption></figure>
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 7.** Out of equilibrium emission from the previously studied Yagi-Uda antenna when only the feeder is emitting (red), only the reflector (blue) and only the directors (black) for gold (7a) and hybrid (7b) antenna. As a reference we also show the equilibrium case in which all the elements are at T=400 K (magenta)</figcaption></figure>
Finally we note that the contribution from the thermal emission of the parasitic elements is particularly evident if we want to maximize the contrast between forward and backward emission. In Fig. 8 it is shown the radiation pattern in the xy-plane at the peak wavelength for the case of the hybrid antenna. We compare the results obtained if only the emission from the feeder is considered (red) with respect to the equilibrium case (black). We note that the backward radiation in the acceptance angle ±15 ° is strongly suppressed if only the feeder is heated. The calculated forward to backward ratio is 25. In the equilibrium case the ratio is reduced to 11. The same argumentation could be performed from the metal YU antenna but in this case the cause of a poor forward to backward ratio is due to the fact that all the elements contribute to the emission and the reflector cannot reflect with the same efficiency the signals coming from elements placed at different distances. In Fig. 8 we show as a reference the radiation pattern at the peak for the all-gold YU antenna (green) exhibiting a pronounced backward emission.

<figure><img src="image_11.png" /><figcaption>**Fig. 8.** Emission pattern in the x-y plane at the peak wavelength for the hybrid YU antenna at the equilibrium temperature of 400 K (black) and for the out of equilibrium case in which only the feeder is heated (red). We note a reduction of backward emission if the contribution from the parasitic elements is not considered. As a reference we also report the emission pattern at peak wavelength for the gold antenna (green).</figcaption></figure>

As mentioned before, one of the most interesting aspect of YU configuration is its robustness with respect to design parameters. This feature is extremely appealing for easy realization of tunable sources with small wavelength differences and good directionality. As an example we report the emission spectra for slightly different schemes. We considered the same antenna, and the same NA for collecting the emitted signal. We only changed the length of the SiC feeder rod. Beside the previously investigated case with L1 = 1.62 µm, we performed calculations for L2 = 1.36 µm and L3 = 1.88 µm. From the calculated emission spectra depicted in Fig. 9 we notice that neither the linewidths nor the efficiencies are affected, there is only a shift in the wavelength of the emission peak. The increased signal for longer rods it is just related to the increasing amount of SiC (see Fig. 9(b)) being the calculated power spectral density normalized with respect to the maximum value obtained for the L1 SiC rod, as usual. Arrays of different YU antennas might be used to finely match emission/absorption lines of various hazardous gases and
used both as sources and as receivers to create compact low cost, miniaturized sensing devices. Finally we note that the hybrid antenna appears differently with respect to the typical YU antenna (as sketched in Fig. 9(a)). Indeed the SiC feeder is smaller than the other elements. This is related to the different effective wavelength for gold and SiC rods.

<figure><img src="image_12.png" /><figcaption>**Fig. 9.** (a) sketch of the hybrid YU antenna with the length of the SiC feeder equal to L2. (b) Normalized emitted spectral power for three hybrid YU antennas obtained by varying the length of the SiC feeder from L1 = 1.62 µm (blue), L2 = 1.36 µm (red) and L3 = 1.88 µm (black)</figcaption></figure>

Despite of the background noise introduced by the not ideal parasitic elements, the coherent emission provided by the SiC antenna feeder makes it possible to boost the signal directivity even with an increased number of directors. However, our calculations show that the best signal to noise ratio is obtained with 3 directors. Concerning a possible experimental verification of the emission from a single YU antenna, we mentioned that the estimated emitted power from the single SiC rod within a cone of NA=0.26 is approximatively 10 pW and our results show that the power emitted by the YU antenna is enhanced with respect to the single SiC rod by a factor of 5. In the supporting info of Ref. [21] a detailed discussion is devoted to the sensitivity of their experimental apparatus and to the detection technique adopted to distinguish the signal emitted by the antenna from the background. The power emitted from the YU antenna is of the same order of the estimated power emitted from their patch antenna, thus it is reasonable to expect that a similar setup could be suitable for further experimental characterization. Depending on the requirements needed for a real application it could be worth developing multiple parameters optimization procedures to evaluate the best trade-off between fabrication costs and performances
and/or considering systems composed of periodic arrays of YU antennas to increase the emitted power. However, in this case the interaction between the emitters as well as their coupling through the substrate should be considered and this will be subject for future work.

## 4. Conclusions

We show that despite of the chaotic nature of thermal radiation it is possible to obtain efficient highly directional and narrow bandwidth thermal antennas in the mid to far IR by adapting the Yagi-Uda scheme with a combination of metallic and polar materials. In particular, we studied hybrid SiC/Au antennas. Our calculations show that it is possible to tune the emission wavelength by properly choosing the resonance of the feeder and to boost the signal and its directivity. The detrimental effect due to thermal emission of the parasitic elements has been taken into account and we have shown that a pure metallic Yagi-Uda antenna could not efficiently work in the 12 µm range, especially if a narrow spectral linewidth of IR radiation is required. This is due to the fact that enhancement and selectivity of thermal emission in metals is related to resonant plasmon polaritons excitation which is efficient only in a wavelength range close to the plasma frequency of the metal (typically in the visible and near IR). Nevertheless, combining high selectivity of SiC antennas in the 12 µm wavelength range (due to resonant phonon polariton excitation) and high scattering cross section and low losses of metals in that wavelength range efficient hybrid thermal Yagi-Uda antennas can be designed. Other polar materials such as SiO2 and MgF2, instead of SiC, could be considered in order to increase the spectral range palette of available sources. This could be possible by simply adapting the proposed scheme and rescaling the size of the elements to fit the desired operation wavelengths according to the reststrahlen band of the considered material. The enhanced performances and the robustness with respect to design parameters make hybrid thermal Yagi-Uda antennas very promising in both transmission and reception mode for integrated, low cost and mass-market gas sensing devices.

## Acknowledgments

M.C. kindly acknowledges Prof. Constantin Simovski for stimulating discussions and suggestions on spatial and temporal coherence of thermal radiation from nanostructured systems. M. C. thanks the CNRS and the group Theory of Light-Matter and Quantum Phenomena of the Laboratoire Charles Coulomb for hospitality during his stay in Montpellier.

M. A. thanks the SAPIENZA University of Rome and the Department of Basic and Applied Sciences for Engineering for hospitality during his stay in Rome under the visiting professor program. M. A. acknowledges support from the Institute Universitaire de France, Paris - France (UE).

## Disclosures

The authors declare no conflicts of interest.

## References

1. J. Hodgkinson and R. P. Tatam, "Optical gas sensing: a review," *Meas. Sci. Technol.* **24**(1), 012004 (2013).

2. A. Lenert, D. M. Bierman, Y. Nam, W. R. Chan, I. Celanović, M. Soljačić, and E. N. Wang, "A nanophotonic solar thermophotovoltaic device," *Nat. Nanotechnol.* **9**(2), 126-130 (2014).

3. A. P. Raman, M. A. Anoma, L. Zhu, E. Rephaeli, and S. Fan, "Passive radiative cooling below ambient air temperature under direct sunlight," *Nature* **515**(7528), 540-544 (2014).

4. G. Sidebotham, Thermal Circuits, In: *Heat Transfer Modeling*. (Springer, Cham, 2015)

5. J. Faist, F. Capasso, D. L. Sivco, C. Sirtori, A. L. Hutchinson, and A. Y. Cho, "Quantum cascade laser," *Science* **264**(5158), 553-556 (1994).

6. W. Li and S. Fan, "Nanophotonic control of thermal radiation for energy application," *Opt. Express* **26**(12), 15995 (2018).

7. R. Messina and M. Antezza, "Scattering-matrix approach to Casimir-Lifshitz force and heat transfer out of thermal equilibrium between arbitrary bodies," *Phys. Rev. A* **84**(4), 042102 (2011).
8. R. Messina, A. Noto, B. Guizal, and M. Antezza, "Radiative heat transfer between metallic gratings using adaptative spatial resolution," *Phys. Rev. B* **95**(12), 125404 (2017).

9. B. Zhao, B. Guizal, Z. M. Zhang, S. Fan, and M. Antezza, "Near-field heat transfer between graphene/hBN multilayers," *Phys. Rev. B* **95**(24), 245437 (2017).

10. R. Messina, P. Ben-Abdallah, B. Guizal, and M. Antezza, "Graphene-based amplification and tuning of near-field radiative heat transfer between dissimilar polar materials," *Phys. Rev. B* **96**(4), 045402 (2017).

11. B. Müller, R. Incardone, M. Antezza, T. Emig, and M. Krüger, "Many body heat radiation and heat transfer in the presence of a non-absorbing background medium," *Phys. Rev. B* **95**(8), 085413 (2017).

12. M. C. Larciprete, A. Belardini, R. Li Voti, and C. Sibilia, "Pre-fractal multilayer structure for polarization- insensitive temporally and spatially coherent thermal emitter," *Opt. Express* **21**(S3), A576-A584 (2013).

13. G. Hamaoui, N. Horny, C. L. Gomez-Heredia, J. A. Ramirez-Rincon, J. Ordonez-Miranda, C. Champeaux, F. Dumas-Bouchiat, J. J. Alvarado-Gil, Y. Ezzahri, K. Joulain, and M. Chirtoc, "Thermophysical characterisation of VO2 thin films hysteresis and its application in thermal rectification," *Sci. Rep.* **9**(1), 8728 (2019).

14. J. Ordonez-Miranda, Y. Ezzahri, J. A. Tiburcio-Moreno, K. Joulain, and J. Drevillon, "Radiative Thermal Memristor," *Phys. Rev. Lett.* **123**(2), 025901 (2019).

15. J.-J. Greffet, R. Carminati, K. Joulain, J.-P. Mulet, S. Mainguy, and Y. Chen, "Coherent emission of light by thermal sources," *Nature* **416**(6876), 61-64 (2002).

16. S. Campione, F. Marquier, J.-P. Hugonin, A. R. Ellis, J. F. Klem, M. B. Sinclair, and T. S. Luk, "Directional and monochromatic thermal emitter from epsilon-near-zero conditions in semiconductor hyperbolic metamaterials," *Sci. Rep.* **6**(1), 34746 (2016).

17. Y. Zhang, M. Antezza, H.-L. Yi, and H.-P. Tan, "Metasurface-mediated anisotropic radiative heat transfer between nanoparticles," *Phys. Rev. B* **100**(8), 085426 (2019).

18. M. C. Larciprete, M. Centini, R. Li Voti, and C. Sibilia, "Asymmetric and tunable infrared emission in metamaterials composed by oriented air voids into a polar material," *J. Mod. Opt.* **66**(3), 299-303 (2019).

19. C. Simovski, S. Maslovski, I. Nefedov, and S. Tretyakov, "Optimization of radiative heat transfer in hyperbolic metamaterials for thermophotovoltaic applications," *Opt. Express* **21**(12), 14988-15013 (2013).

20. A. Lochbaum, Y. Fedoryshyn, A. Dorodnyy, U. Koch, C. Hafner, and J. Leuthold, "On-Chip Narrowband Thermal Emitter for Mid-IR Optical Gas Sensing," *ACS Photonics* **4**(6), 1371-1380 (2017).

21. C. Li, V. Krachmalnicoff, P. Bouchon, J. Jaeck, N. Bardou, R. Haïdar, and Y. De Wilde, "Near-Field and Far-Field Thermal Emission of an Individual Patch Nanoantenna," *Phys. Rev. Lett.* **121**(24), 243901 (2018).

22. C. Arnold, F. Marquier, M. Garin, F. Pardo, S. Collin, N. Bardou, J.-L. Pelouard, and J.-J. Greffet, "Coherent thermal infrared emission by two-dimensional silicon carbide gratings," *Phys. Rev. B* **86**(3), 035316 (2012).

23. K. Joulain, Y. Ezzahri, J. Drevillon, B. Rousseau, and D. de Sousa Meneses, "Radiative Thermal Rectification between SiC and SiO2," *Opt. Express* **23**(24), A1388-A1397 (2015).

24. J. A. Schuller, T. Taubner, and M. L. Brongersma, "Optical antenna thermal emitters," *Nat. Photonics* **3**(11), 658-661 (2009).

25. M. C. Larciprete, M. Centini, R. Li Voti, and C. Sibilia, "Selective and tunable thermal emission in metamaterials composed of oriented polar inclusions," *J. Opt. Soc. Am. B* **34**(7), 1459-1464 (2017).

26. F. Pristera, M. Halik, A. Castelli, and W. Fredericks, "Analysis of explosives using infrared spectroscopy," *Anal. Chem.* **32**(4), 495-508 (1960).

27. R. Furstenberg, C. A. Kendziora, J. Stepnowski, S. V. Stepnowski, M. Rake, M. R. Papantonakis, V. Nguyen, G. K. Hubler, and R. A. McGill, "Stand-off detection of trace explosives via resonant infrared photothermal imaging," *Appl. Phys. Lett.* **93**(22), 224103 (2008).

28. E. Tucker, J. D' Archangel, M. Raschke, E. Briones, F. J. González, and G. Boreman, "Near-field mapping of dipole nano-antennacoupled bolometers," *J. Appl. Phys.* **114**(3), 033109 (2013).

29. J. A. Bean, A. Weeks, and G. D. Boreman, "Performance Optimization of Antenna-Coupled Al/AlOx/Pt Tunnel Diode Infrared Detectors," *IEEE J. Quantum Electron.* **47**(1), 126-135 (2011).

30. A. G. Curto, G. Volpe, T. H. Taminiau, M. P. Kreuzer, R. Quidant, and N. F. van Hulst, "Unidirectional Emission of a Quantum Dot Coupled to a Nanoantenna," *Science* **329**(5994), 930-933 (2010).

31. T. Wang, P. Li, D. N. Chigrin, A. J. Giles, F. J. Bezares, O. J. Glembocki, J. D. Caldwell, and T. Taubner, "Phonon-Polaritonic Bowtie Nanoantennas: Controlling Infrared Thermal Radiation at the Nanoscale," *ACS Photonics* **4**(7), 1753-1760 (2017).

32. H. Li, H. Peng, C. Ji, L. Lu, Z. Li, J. Wang, Z. Wu, Y. Jiang, J. Xu, and Z. Liu, "Electrically tunable mid-infrared antennas based on VO2," *J. Mod. Opt.* **65**(15), 1809-1816 (2018).

33. Y. Xiang, S. Amarie, W. Cai, W. Luo, W. Wu, M. Ren, X. Zhang, and J. Xu, "Real-space mapping of mid-infrared near-field of Yagi-Uda antenna in the emission mode," *Opt. Express* **27**(4), 5884-5892 (2019).

34. B. Yu, J. Li, and S. Shen, "Directional control of narrow-band thermal emission from nanoantennas," *J. Photonics Energy* **9**(3), 1 (2019).

35. M. Centini, A. Benedetti, M. C. Larciprete, A. Belardini, R. Li Voti, M. Bertolotti, and C. Sibilia, "Midinfrared thermal emission properties of finite arrays of gold dipole nanoantennas," *Phys. Rev. B: Condens. Matter Mater. Phys.* **92**(20), 205411 (2015).

36. S. M. Rytov, Y. A. Kravtsov, and V. I. Tatarskii, *Principles of Statistical Radio-Physics* 3 (Springer-Verlag, Berlin Heidelberg, 1989).
37. E. D. Palik, *Handbook of Optical Constants of Solids* (Academic Press, Inc., Orlando, Florida, 1985).

38. Constantine A Balanis, *Modern Antenna Handbook*. (John Wiley and Sons, New York, 2008)

39. D. Cheng and C. Chen, "Optimum element spacings for Yagi-Uda arrays," *IEEE Trans. Antennas Propag.* **21**(5), 615-623 (1973).