0. Introduction
With the rapid advancement of miniaturized devices, antenna regulation of antennas has garnered significant attention. Dynamic modulation of the far-field to achieve versatile control is highly desirable. Specifically, negative refraction, tunable metasurfaces, optical switches, tunable cavities, and coherent perfect absorbers have been realized using indium tin oxide (ITO) as epsilon-near-zero (ENZ) materials and aluminum gallium arsenide (AlGaAs). [
1,
2,
3,
4]. However, there has been no report on temperature-controlled far-field radiation of antennas to date. Significantly, beyond purely optical considerations, the design of temperature-dependent nanoantennas must also account for their electromagnetic response. This can be addressed by solving dynamic regulation problems, such as solutions for exchange energy and anisotropies. The far-field radiation patterns are determined by the refractive index, which is temperature-dependent. By increasing the temperature of an antenna, multiople decomposition suggests that the temporal and spatial coherence of the far-field radiation can be engineered through careful material selection or dimensional adjustments. The engineering of far-field radiation is of great interest for applications in lighting, thermoregulation, energy harvesting [
5,
6], tagging and radiation therapy [
7].
Recent studies have demonstrated that ENZ materials and AlGaAs exhibit an exceptionally large temperature-dependent refractive index (see
Figure 1(a)) [
8,
9]. Consequently, ENZ and AlGaAs materials offer a new platform for optically tuning the material response within a subpicosecond timescale [
2]. One way to increase the temperature modulation speed is to forgo heating the material to obtain hot electrons that has been demonstrated using ultrafast pump pulses with modulation speeds in the 100 s of femiseconds in graphene [
10], and several picoseconds in metals such as gold and tungsten [
11].
III-V semiconductors have the advantage of a short electron lifetime compared to silicon, making them suitable for switching applications.
features a broad transparency windwo ranging from the near to mid-infrared, from 0.9 to 17
[
12]. Additionally,
exhibits a high relative permittivity. The resonance of ITO disk is weaker than that of
disk with the same diameter due to higher absorption losses of ITO. However, the absorption cross section of ITO disk remains evident. The electromagnetic response of the hybrid nanoantenna was further confirmed by measuring the optical Kerr effect as a function of temperature.
In this paper, we utilize ITO and to engineer the far-field control of radiation patterns. First, we present the basic formalism describing the refractive indices of ITO and , and study the scattering and far field properties to realize dynamically reconfigurable far field radiation pattern based on single nanophotonic structures. We then emphasize the hybridization of ITO and to exploit their properties for energy efficiency and sustainability, thereby engineering the far-field radiation effectively. The relevance of this study is primarily theoretical, as the relationship between far-field radiation distribution and the transition between super-absorption and super-scattering remains largely unexplored using the calibrated temperature-index models. We demonstrate that this aspect, although not explicitly accounted for, may be spatially transferable when applied to sufficiently large temperature range. Methods to achieve temperature control of far-field radiation have been implemented by heating, adjusting light intensity, or altering voltage as the refractive index changes with electron temperature. Notably, the far-field variation not only changes by half but also shifts the scattering direction from the original x-axis to the y-axis, resulting in a more significant control range. Clearly, nanoscale electromagnetics is one of the technologies with the potential to provide such functionality, including reception, emission, and more generally, the control of light at the nanoscale.
1. Materials and Methods
According to the optical properties of ITO as a function of electron temperature, the real and imaginary components of the refractive index at
nm are shown by the blue lines in
Figure 1. The change in the real part of the refractive index with intensity is approximately 0.22, and the linear refractive index is 0.45. The calculated refractive index of
at
is depicted in
Figure 1 by the red line. The change in the real part of the refractive index with intensity is approximately 0.22, while the linear refractive index is 3.32.
Figure 1.
Temperature-dependent materials. Refractive indices of and as a function of temperature at a wavelength of nm.
Figure 1.
Temperature-dependent materials. Refractive indices of and as a function of temperature at a wavelength of nm.
The temperature-dependent refractive index of ITO is expressed as,
where
is the high frequency permittivity,
N is the carrier density,
e is the electron charge,
and
are the electron effective mass and mobility, which vary with electron temperature. The plasma frequency
with the
. The parameters for the commercial ITO are:
,
,
,
,
. The mobility
. The real part of the dielectric constant
of ITO is zero at 1240 nm that is the ENZ wavelength
[
13,
14]. According to the optical properties of ITO as a function of electron temperature on the model, the real and imaginary components of the refractive index are plotted in
Figure 1 in blue at
[
13,
15]. The change of the real part of the refractive index with intensity is about 0.22, the linear refractive index is 0.45.
While, the Sellmeier equation of
with temperature dependence can be written as
The first brackets term is the refractive index at room temperature. The second brackets represents the temperature-dependent Sellmeier coefficient. The temperature
T, wavelength
enable us to calculate the index for
, where
. This Sellmeier can be used at 1240
and 2000
[
16]. The calculated index of
are plotted in
Figure 1 in red line at
[
16]. The change of the real part of the refractive index with intensity is about 0.22, the linear refractive index is 3.32.
2. Results
2.1. Tuning the Scattering and Far Field for Temperature-Dependent Materials
The proposed structure, consisting of three-dimensional (3D) ITO disk with height
nm and diameter
nm, is schematically shown in the inset of
Figure 2(a). The nonlinear material made of ITO is illuminated by a plane wave with electric field
, where
is the wavenumber,
is the angular frequency,
is the amplitude of incident field, and c.c. means complex conjugate.
is the free-space intensity of the incident plane wave.
is the electric field inside the ITO. The scattering properties of the ITO nanodisk, including its electric and magnetic dipole contributions, are depicted in
Figure 2(a). The nanodisk is discretized using a 10 nm hexagonal mesh following a convergence test. As shown in
Figure 2(a), the scattering cross-section exhibits minimal variation with temperature. The primary contribution to the multipole scattering arises from the electric dipoles (ED), while the contributions from magnetic dipole (MD), electric quadrupole (EQ), magnetic quadrupole (MQ) are relatively minor. Consequently, the far-field radiation pattern, influenced by the interference between the multipoles, remains largely unchanged. Due to the imaginary component of ITO, the antenna also exhibits an absorption cross-section.
Figure 2(b) illustrates that both the absorption and scattering cross-sections vary slightly with temperature, leading to a corresponding change in the extinction cross-section. It is evident from
Figure 2(b) that the absorption cross-section is dominant, approximately twice that of the scattering cross-section.
The optical response of radiation pattern can be explained by the multipole expansion of the induced displacement current
[
17,
18]. The excited multipole can be obtained by the induced nonlinear displacement current
. The exact multipole expansion introduced by [
19,
20]
where
,
,
,
, and
are the electric dipole (ED), magnetic dipole (MD), electric quadrupole (EQ), and magnetic quadrupole (MQ) multipole moments, respectively.
is the
nth spherical Bessel function. The total scattering cross section can be calculated by [
19,
20]
which is the contribution of each multipole moment (marked ’total’ in Figs. 1(c) and 2(c)). The interference between multipoles can satisfy the generalized Kerker condition. The far field corresponding to the radiation pattern (
) in space is
where
r,
, and
are the spherical coordinates.
k is the wavenumber,
c is the speed of light.
The temperature-induced change in the refractive index results in minor variations in the scattering. These variations are due to the contributions of multipoles, which remain relatively constant and exhibit an almost unchanged far-field radiation pattern, as illustrated in
Figure 2(c1)-(c2). In the study of the optical properties of a single ITO nanostructure, both scattering and absorption are influenced by temperature. However, these changes are minimal, primarily because the refractive index of ITO changes by less than 1, resulting in a very small difference between the refractive index of ITO and the background refractive index (set to 1).
Similarly, we use an
antenna of the same size to study its scattering and far-field characteristics.
allows for even greater control over the scattering properties of antennas. As a lossless dielectric, we focus on the scattering cross-section.
Figure 3(a) presents the scattering cross-section as a function of temperature for different height-to-diameter ratios (
). At a specific diameter and temperature, the scattering cross-section is not proportional to height variations, exhibiting irregular fluctuations. Furthermore, for a fixed size, the pattern of the scattering cross-section changing with temperature is not consistent. Taking
and
as an example, the scattering cross-section initially decreases and then increases. More diverse changes are shown in Supplementary Material Figure S2. When the temperature is below 866 K, the electric dipole contribution decreases while the magnetic dipole contribution gradually increases, leading to a reduction in far-field radiation in the
direction and an increase in the
direction. Approaching 866 K, the scattering is the smallest due to the smallest contribution of electric dipole. Correspondingly, the contribution of the magnetic quadrupole is the largest. When the temperature exceeds 866 K, the electric dipole contribution steadily increases, resulting in forward scattering in the far-field radiation pattern. The significant changes in electric dipoles before and after 866 K profoundly affect the far-field radiation patterns. The radiation pattern of
changes more than that of ITO due to its refractive index changing by 0.22 with temperature. Additionally, its refractive index is around 3.5, which is significantly different from the background refractive index.
2.2. Tuning the Scattering and Far Field of Hybrid Structures for Opposite Illumination
To realize the switching between tunable far-field radiation and superradiant absorption, we integrated both materials and investigated their optical characteristics under the opposite electric fields. A hybrid nonlinear antenna composed of ITO and lossless
is depicted in
Figure 4(a). Owing to the broken inversion symmetry, the hybrid antenna demonstrates a biaisotropic response.
In
Figure 4(b), it is observed that when light illuminates a material from two opposing directions, the extinction cross-section remains identical [
20]. Due to the loss characteristics of ITO, associated with optical bistability, both scattering and absorption cross-sections are contingent upon the direction of incident light. Initially, commencing from room temperature, the scattering cross-section of the hybrid structure exhibits a preference for illumination from below compared to illumination from above. However, nearing approximately 866 K, the scattering cross-sections converge, and with further temperature elevation, a reversal occurs, where the scattering cross-section for illumination from above surpasses that for illumination from below.
Figure 4c,d depict the multipole decomposition. Under bottom illumination, electric dipoles consistently dominate with variations in temperature. Under top illumination, electric dipoles mainly dominate at
K and
K, while both electric and magnetic quadrupoles contribute more at 866 K
K.
Figure 5.
Radiation patterns of hybrid antenna. (a1)-(a5) Evolution of representative radiation patterns for the bottom illumination direction, i.e., , increases 299 to 2000 . (b1)-(b5)The radiation patterns for the bottom illumination direction, i.e., , increases from 299 to 2000 .
Figure 5.
Radiation patterns of hybrid antenna. (a1)-(a5) Evolution of representative radiation patterns for the bottom illumination direction, i.e., , increases 299 to 2000 . (b1)-(b5)The radiation patterns for the bottom illumination direction, i.e., , increases from 299 to 2000 .
The induced electric, magnetic dipoles and quadrupoles interfere constructively (destructively) in the backward (forward) direction. Under bottom illumination, the radiation pattern initially consists of two lobes distributed symmetrically about the z-axis without backward illumination at room temperature. As temperature rises, backscattering gradually emerges on the radiation map. Further temperature increase induces bidirectional changes in the radiation pattern, with backward (forward) radiation decreasing (increasing) correspondingly. Ultimately, significant tunability of the induced multipole moments in the hybrid antenna by temperature enables control over radiation patterns from a bidirectional to unidirectional pattern. Specifically, under bottom illumination, the initially omnidirectional radiation at room temperature transforms into a three-lobe pattern distributed symmetrically about the y-axis and the upper side of the z-axis. With continued temperature elevation, the radiation pattern evolves from three lobes to four, predominantly radiating from the upper side of the z-axis. Subsequently, with further temperature increase, radiation in the y-direction of the four-lobe pattern diminishes, resulting in bidirectional radiation with forward scattering. Consequently, the nonreciprocal characteristics of hybrid structures arise from the distinct electromagnetic couplings generated.
3. Conclusion
In summary, this study investigated the influence of temperature variations on far-field radiation patterns through numerical simulations employing a spatially explicit temperature-index model. A hybrid structure comprising two materials, each exhibiting temperature-dependent refractive indices, was analyzed. The refractive index change from room temperature to 2000 K was determined to be 0.22. This structure capitalizes on the temperature-induced loss in ITO to achieve switching between absorption and scattering, while also leveraging the strong scattering characteristics resulting from the high refractive index of . Temperature adjustments induce diverse far-field radiation distributions. Furthermore, the operational wavelength range can be broadened significantly, and the law of refractive index changing with temperature has been given. Moreover, the methodology employed to determine the temperature-dependent refractive indices of ITO and at 1240 nm can be extended to other wavelengths, facilitating multi-parameter control of far-field radiation and enhancing regulation replication, thereby expanding the tunable wavelength range.
Author Contributions
The following statements should be used “Conceptualization, L. C.; methodology,Lin Cheng; software, F. W.; validation, L. C.; formal analysis, F. W.; investigation, L. C.; resources, F. W.; data curation, F. W.; writing—original draft preparation, L. C.; writing—review and editing, K. H.; visualization, K. H.; supervision, L. C.; project administration, all; funding acquisition, L. C., F. W., K. H. All authors have read and agreed to the published version of the manuscript.
Funding
We would like to acknowledges the support the National Natural Science Foundation of China (62305312), the Natural Science Foundation of Shanxi Province, China (202203021222021), Research Project Supported by Shanxi Scholarship Council of China(2312700048MZ) and the fellowship of China Postdoctoral Science Foundation (2022M722923). Shanxi Provincial Teaching Reform and Innovation Project (PX-624261). Shaanxi Fundamental Science Research Project for Mathematics and Physics (23JSQ020). Shanxi Provincial Teaching Reform and Innovation Project(J20221525).
Data Availability Statement
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
Conflicts of Interest
The authors declare no conflict of interest.
Appendix A
Appendix A.1. Far fields and Radiation Patterns
Figure A1.
The extinction, absorption and scattering cross section of the ITO antenna as a function of height H (a) with diameter nm. (b) with diameter nm. (c) When diameter nm, the contribution of different multipole moments ED, MD, EQ and MQ, respetively. (d) When diameter nm, the contribution of different multipole moments ED, MD, EQ and MQ, respetively.
Figure A1.
The extinction, absorption and scattering cross section of the ITO antenna as a function of height H (a) with diameter nm. (b) with diameter nm. (c) When diameter nm, the contribution of different multipole moments ED, MD, EQ and MQ, respetively. (d) When diameter nm, the contribution of different multipole moments ED, MD, EQ and MQ, respetively.
Using the induced multipole moments in Eq.(3), the far field can be found [
21,
22]
where
are the radial distance, polar angle, and azimuthal angle, respectively. Considering the contribution from all multipole moments (up to magnetic quadrupole), the electric field corresponding to the radiation pattern (
) in the
-plane, i.e.,
, can be written as
from which, the radiation pattern (
) can be obtained. Equations (
A1) and (
A2) are used to plot the radiation patterns in the main text.
Appendix A.2. Antennas Based on ENZ Materials
In this section, we discuss the temperature-dependant scattering cross sections of the ITO antenna.
We use finite difference time domain (FDTD) to analyze the two-dimensional (2D) far-field radiation patterns at
nm and
nm, respectively. Here, Periodic boundary conditions are employed along the
x and
y axes; Perfectly matched layer (PMLs) are applied along the
z axis. The structure is excited by a single plane wave source at 1240 nm for normal incidence. In order to obtain accuracy results, the mesh size is optimized after a converging test. The 2D radiation pattern are shown in
Figure A2.
Figure A2.
(a) The radiation patterns in (a) XZ plane (b) YZ plane (c) XY plane, with the diameter nm. (d)-(e) Same as (a)-(c) with the diameter nm.
Figure A2.
(a) The radiation patterns in (a) XZ plane (b) YZ plane (c) XY plane, with the diameter nm. (d)-(e) Same as (a)-(c) with the diameter nm.
Figure A1(a) and
Figure A1(b) shows an ITO antenna when illuminated by an
x-polarized plane wave propagating in the
z direction. The scattering, absorption and extinction cross section increases as the height of ITO increases.
Figure A1(c) and
Figure A1(d) show the contribution of different multipole moments. The proportion of the individual multipoles changes, and the coherent interference between the multipoles causes the far-field radiation to change from bidirectional to unidirectional (see
Figure A2(a)-
Figure A2(c)). When the
D increases to 1000 nm, the angle of radiation gradually decreases (see
Figure A2(d)-
Figure A2(f)). Consequently, the antenna exhibits a nearly unidirectional radiation pattern with very small backscattering [see
Figure A2(a)-
Figure A2(d)], a phenomenon known as the generalized Kerker effect [
22,
23].
Figure A3.
(a)-(l) When the diameter nm, the contribution of different multipole moments ED, MD, EQ, MQ, and the scattering cross section (normalized to ) as a function of the temperature for different height H.
Figure A3.
(a)-(l) When the diameter nm, the contribution of different multipole moments ED, MD, EQ, MQ, and the scattering cross section (normalized to ) as a function of the temperature for different height H.
Appendix A.3. Antennas Based on Al 0.18 Ga 0.82 As
In this section, we show the temperature-dependent multipole moments for different height. The scattering cross section increases when the height at nm. When nm, the scattering cross section decreases and then increases. The scattering cross-section at different height changes with temperature, which is due to the obvious change of the induced multipole with the change of height. These trends provide a reference for the design of hybrid structures.
Appendix A.4. Field Distribution of the Hybrid Nonreciprocal Nonlinear Antennas
When the hybrid antenna is illuminated from the top, compared to the bottom illumination, different field distribution in the
plane [see
Figure A4(a)-
Figure A4(d).
Figure A4.
(a) The field distribution in plane of the hybrid structure for opposite illumination for the bottom illumination direction when K.(b) The field distribution in plane of the hybrid structure for opposite illumination for the bottom illumination direction when K. (c)-(d) Same as (a)-(b) for the top illumination.
Figure A4.
(a) The field distribution in plane of the hybrid structure for opposite illumination for the bottom illumination direction when K.(b) The field distribution in plane of the hybrid structure for opposite illumination for the bottom illumination direction when K. (c)-(d) Same as (a)-(b) for the top illumination.
References
- N. Kinsey, C. Devault, J. Kim and A. Boltasseva, Optica2015.
-
M. Z. Alam, I. D. Leon and R. W. Boyd,Science, 2016, 352, 795.
- Argyropoulos, Christos, Chen, Pai-Yen, D’Aguanno, Giuseppe, Engheta, Nader, Alu and Andrea, Physical Review B Condensed Matter & Materials Physics,2012.
-
E. Mobini, R. Alaee, R. W. Boyd and K. Dolgaleva,ACS Photonics, 2021, 8, 3234–3240.
-
Jack, Ng, Huanyang, Chen, C., T. and Chan,Optics Letters, 2009.
- A. E. Miroshnichenko and M. I. Tribelsky, Physical Review Letters,2018,120, 033902.
-
S. Clement, J. M. Campbell, W. Deng, A. Guller, S. Nisar, G. Liu, B. C. Wilson and E. M. Goldys,Advanced Science, 2020, 7.
-
J. P. Kim and A. M. Sarangan,Optics Letters, 2007, 32, 536–538.
-
J. Paul, M. Miscuglio, Y. Gui, V. J. Sorger and J. K. Wahlstrand,Optics Letters, 2020.
- C. H. Lui, R. F. Mak, R. Shan and R. F. Heinz, Physical Review Letters,2010,105, 127404.
- M. A. A. O. S. Y. C. P. A. Fortov, JETP Letters,2015,101, 598–602.
- S. Gehrsitz and F. K. Reinhart, Journal of Applied Physics,2000,87, 7825–7837.
- K. J. Russell, T. L. Liu, S. Cui and E. L. Hu, Nature Photonics,2012,6, 459–462.
- R. Orad, G. Enno, A. M. Zahirul, D. L. Israel, U. Jeremy and R. W. Boyd, Optics Letters,2017,42, 3225.
- J. Paul, M. Miscuglio, Y. Gui, V. J. Sorger and J. K. Wahlstrand, Opt. Lett.,2021,46, 428–431.
- J. P. Kim and A. M. Sarangan, Opt. Lett.,2007,32, 536–538.
- R. Alaee, C. Rockstuhl and I. Fernandez-Corbaton, Advanced Optical Materials,2019,7, 1800783.1–1800783.14.
- R. Alaee, C. Rockstuhl and I. Fernandez-Corbaton, Optics Communications,2017,407.
- R. Alaee, C. Rockstuhl and I. Fernandez-Corbaton, Opt. Commun.,2018,407, 17 – 21.
-
R. Alaee, C. Rockstuhl and I. Fernandez-Corbaton,Adv. Opt. Mater., 2019, 7, 1800783.
-
J. D. Jackson,Classical Electrodynamics, Wiley, 1999.
- R. Alaee, R. Filter, D. Lehr, F. Lederer and C. Rockstuhl, Opt. Lett.,2015,40, 2645–2648.
- W. Liu and Y. S. Kivshar, Opt. Express,2018, 26, 13085–13105.
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).