Enhanced terahertz magnetic dipole response by subwavelength fiber

Dielectric sub-wavelength particles have opened up a new platform for realization of magnetic light. Recently, we have demonstrated that a dipole emitter by a subwavelength ﬁber leads to an enhanced magnetic response. Here, we experimentally demonstrate an enhanced magnetic dipole source in the terahertz frequency range. By placing the ﬁber next to the hole in a metal screen, we ﬁnd that the radiation power can be enhanced more than one order of magnitude. The enhancement is due to the excitation of the Mie-type resonances in the ﬁber. We demonstrate that such a system is equivalent to a double-ﬁber system excited by a magnetic source. This coupled magnetic dipole and optical ﬁber system can be considered a unit cell of metasurfaces for manipulation of terahertz radiation and is a proof-of-concept

Dielectric sub-wavelength particles have opened up a new platform for realization of magnetic light. Recently, we have demonstrated that a dipole emitter by a subwavelength fiber leads to an enhanced magnetic response. Here, we experimentally demonstrate an enhanced magnetic dipole source in the terahertz frequency range. By placing the fiber next to the hole in a metal screen, we find that the radiation power can be enhanced more than one order of magnitude. The enhancement is due to the excitation of the Mie-type resonances in the fiber. We demonstrate that such a system is equivalent to a double-fiber system excited by a magnetic source. This coupled magnetic dipole and optical fiber system can be considered a unit cell of metasurfaces for manipulation of terahertz radiation and is a proof-of-concept of a possibility to achieve enhanced radiation of a dipole source in proximity of a sub-wavelength fiber. It can also be scaled down to optical frequencies opening up promising avenues for developing integrated nanophotonic devices such as nanoantennas or lasers on fibers. © 2018 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/1.5010348

I. INTRODUCTION
Subwavelength dielectric structures provide a platform for strong light matter interaction with breakthrough application in nanoantennas, nanolasers, metasurfaces, and metadevices. 1 At optical frequencies, it has been shown that dielectric nanoparticles lead to strong light localization when excited by a plane wave due to excitation of Mie resonances. 1 These resonances show strong electric and/or magnetic responses that can be leveraged in designing metadevices. As a new twist, we have recently shown that a dipole emitter can also excite the resonances of a nanofiber and lead to strong electric and/or magnetic responses. 2 These responses were represented by the enhanced resonant emission dominated by magnetic dipole excitation. This novel 2D fiber-based platform, a Electronic mail: s.atakaramians@unsw.edu.au invariant along the fiber and confined with subwavelength scales in the other dimension, not only has enhanced electric and magnetic responses similar to subwavelength particles due to the excitation of Mie resonances but also allows harnessing the guided waves of fibers opening up a new horizon for designing integrated metadevices. Modification of the radiation characteristics of a source by the electromagnetic environment is one of the most important discoveries of the twentieth century. In particular, change of the emitter's lifetime by a resonator is associated with the Purcell effect, 3 which was originally demonstrated for magnetic dipole (MD) type sources in the microwave frequency range. While the modification of an electric dipole (ED) source in the optical domain by the photonic systems has been studied extensively, the modification of an MD source has gained momentum only recently with the aim of developing ultra-efficient optical devices using the magnetic nature of light. 1,4 There have been rigorous theoretical studies 5-10 and experimental confirmation [11][12][13][14][15] of ED emission demonstrating that the strong Purcell effect can be achieved when the emitter is placed within microcavities, 11,16 in photonic crystals, 11,12 in hot spots of nanoantennas, [17][18][19] fibers, 10 and metamaterials. 7, 13,15 Recently there has been great interest in investigating the MD emission enhancement, 20 with the aim of exploiting the Purcell effect in the visible spectrum. Several effects were observed, including strong MD transitions, 21-23 possibility of selective excitation of these transitions 24,25 and strong magnetic response of nanoparticles. 1 Similarly, the Purcell enhancement of MD emitters located in or near various photonic systems 20 including bulk materials, planar structures, 20 maser cavity, 26 nanoplasmonic 27,28 and all dielectric 29,30 nanoantennas, and metamaterials [31][32][33] has been studied in either microwave or optics. To the best of our knowledge, there is no report on the Purcell factor enhancement of MD emitters in the vicinity of fibers, while the enhancement of ED emitters in the vicinity of fibers has been studied. 2,9,10 Furthermore, there is no report on the MD Purcell effect in the THz range as the realization of MD antennas at higher frequencies (including THz range and optics) are difficult, both because of the reduced sizes of the antennas as well as difficulty associated with feeding such antennas.
In this paper, we experimentally demonstrate the magnetic dipole radiation enhancement for the structure containing a hole in a metallic screen 34 and a dielectric subwavelength fiber, as shown in Fig. 1. This experiment is not only the first proof-of-concept of radiation enhancement of an MD dipole source in the vicinity of a subwavelength fiber but also the first observation of the Purcell effect for a magnetic dipole source in the THz frequency range. This hybrid platform of aperture and fiber can be further expanded into an array configuration opening up a new type of hybrid metasurfaces for manipulating THz and higher frequency radiation. Moreover, utilizing guided modes of the fiber will open up new avenues in developing integrated devices.

A. Subwavelength aperture as a magnetic dipole source
The diffraction of a plane wave by an aperture in a perfectly conducting screen can be represented as a radiation of a series of multipole sources. 35,36 When the dimension of the aperture is small compared to the wavelength, the aperture acts as the ED and MD source of radiation with the contributions of higher multipoles being negligible. 34,35,37 The effective ED and MD moments can be calculated in terms of integrals of the tangential electric field over the opening of the aperture as follows: 34 where E tan is the tangential electric field in the aperture, x is the position vector in the xy-plane,ẑ is the unit vector normal to the aperture plane, and the integration is over the area of aperture. It has been shown that the approximation of a subwavelength aperture with effective dipole sources is only valid when the diameter is less than 0.3λ. 38 Figure 2 shows the normalized ED and MD moments of a 300 µm diameter aperture in the z = 0 plane as a function of frequency. We use the commercial software CST Microwave Studio to calculate the tangential component of the electric field in the aperture, when it is excited with an x-polarized plane wave propagating in the z direction. The relative aperture diameter to the wavelength will be 0.3 at f = 0.3 THz, beyond which the approximation breaks down [gray shaded region in Fig. 2(b)]. The MD components are normalized to the maximum of y magnetic moment MD max , while the electric component is normalized to MD max /c to be unit less, where c is the speed of light in vacuum. As one clearly observes, the MD y component is the dominant component (MD y /MD x > 100 and MD y /ED z > 1000). This is consistent with what has been reported in the literature, 34,35,37 i.e., for a normal incident plane wave, where the electric and magnetic fields of excitation are parallel to the film, the scattering properties of a sub-wavelength aperture in a metallic film are mainly determined by a magnetic dipole, which is directed appositely to the magnetic field component of the incident plane wave [ H in Fig. 2(a)].
Although one could argue that the total radiation pattern of the system is angle dependent, i.e., perpendicular to the metallic plane MD y moment is maximum and ED z is zero, while at grazing angle MD y is zero and ED z is in its maximum, it is safe to consider the system as a MD y source for this work as, first, MD y is more than three orders of magnitude stronger and, second, as discussed below, we are investigating the power flowing through a flat surface parallel to the metallic film, which is mainly along the direction of maximum of MD y and zero of ED z radiations.

B. Aperture-fiber system: Experiment and simulation
We use THz time domain spectroscopy (THz-TDS) to conduct our experiment, enhanced MD radiation in the vicinity of sub-wavelength fibers. The TDS system comprises two photoconductive antennas (a common H-shaped antenna with a Si lens on top for generation and a near-field microprobe tip 39 for detection) pumped by a mode-locked Ti:sapphire laser. The near-field microprobe introduces a minimum distortion to the field structure as was shown in previous studies. [40][41][42] The generated THz electric field is linearly polarized [in this case, x-polarized as shown in Fig. 3(a)] and it is focused on the center of a 300 µm aperture in a very thin (<100 µm) copper plate using two plastic lenses. The detector is moved in the xy-plane to generate a near-field raster image of the transmitted electric field. The scan is conducted over a square area of 850 × 850 µm 2 with 0.03 mm step size in both directions. The same step size was used for time delay leading into 14 GHz frequency resolution. A 300 µm diameter soft glass microfiber, known as F2, is used for this experiment. The chosen soft glass has a relatively medium refractive index (n = 2.7) and material loss <10 cm 1 at the interested frequency range. 43 It is worth highlighting that the numerical modeling has shown that the maximum enhancement in the transmission is achieved when the diameter of the fiber is chosen to be identical to that of the aperture. We attribute this to the effective coupling of the aperture emission into resonance modes of the system.
The experiment is conducted in three steps. First, the aperture is positioned at the focal point of the beam and the center of the raster scan is determined, where the maximum signal is recorded at the detector. Second, the fiber is positioned on the aperture and the raster scan is taken at 10 µm distance from the fiber (310 µm from the brass plate). Third, the fiber is removed and a reference raster scan is taken from the aperture only at 310 µm distance from the brass plate. All three steps are repeated for two orientation of the fiber [ Fig. 3(a)], i.e., fiber parallel to the excitation (fiber along x) and fiber perpendicular to the excitation (fiber along y). Figure 3(b) shows the measured normalized power integrated over the square area of the raster scan (850 × 850 µm 2 ) for aperture only (black solid line), and when the fiber is oriented along the x-axis (red solid line) and along the y-axis (green solid line). The curves are normalized to the maximum integrated power through the aperture. Figure 3(b) also shows the simulated normalized power integrated on the same area (dotted lines), where similarly the curves are normalized to the maximum integrated power through the aperture. As expected for a bare aperture, the power gradually increases and reaches its maximum when the diameter is equal to half of the operating wavelength (a = λ/2), i.e., f = c/(2a) = 0.5 THz. With the fiber, the transmitted signal is enhanced for frequencies below 0.36 THz for both orientations. As discussed in our previous studies, 2,10 the enhancement in transmission is due to excitation of Mie resonances in the fiber cross section, which is also the position of the first transverse electric (TE) and magnetic (TM) whispering gallery modes (WGMs) formed at the cross section of the fiber.
The vertical arrows in Fig. 3(b) indicate the points where the signal level falls within the noise level of the system for each measurement and thus experimental results beyond the arrows cannot be trusted. The shift towards lower frequencies of the arrows is consistent with the enhancement of the transmission through the fibers. The experimental and simulation results do not match well above 0.4 THz, the region which is off interest to this work. We attribute this discrepancy between the simulation and experiment to the fact that at higher frequencies the beam width (FWHM = 500 µm at 0.4 THz) becomes comparable to the aperture size and therefore the field uniformity assumption in the aperture is not valid anymore. Figure 3(c) shows the measured (first column) and numerically calculated (second column) normalized electric field intensity at 0.3 THz for aperture only (first row), microfiber along x (second row), and microfiber along y (third row). The enhancement of the transmitted signal through the aperture-microfiber system can also be visually observed in Fig. 3(c). Moreover, there is a focusing effect perpendicular to the direction of the microfiber due to the focusing effect in the near field. 44 See the supplementary material for the time evolution of THz radiation. Figure 3(d) shows the normalized integrated power for both orientations of microfibers, where curves are normalized to aperture only [black curves in Fig. 3(b)]. We have extrapolated the normalized power of the bare aperture (using the simulation results) for the frequencies in the noise level of the system so that we can calculate the enhancement ratio up to the aperture-fiber system noise level for each case. A relatively good agreement between measured (solid lines) and numerical (dashed lines) results is observed. Two peaks at 0.15 and 0.32 THz correspond to the microfiber along x. At the position of the first peak (at 0.15 THz), more than 12 times enhancement is achieved experimentally, while the numerical simulation predicts more than 20 times enhancement. The enhancement (7 times) and position (at 0.32 THz) of the second peak are in good agreement with the numerical results. When the fiber is along y, there exists a peak at 0.244 THz where the signal is enhanced 5 times, while the numerical simulation predicts 7 times enhancement at 0.23 THz. Note that the near-field probe measures the x-component of the electric field; therefore all the simulation results presented so far are also considering the x-component of the electric field only.
There are a few factors that can contribute to the discrepancy between the measured (solid curves) and simulated (dashed curves) results: the relative alignment of the microfiber and aperture, the extrapolation of aperture only curve ( f < 0.25 THz), and the signal being close to the noise level of the system. The simulation results confirm that a 30-50 µm misalignment of the microfiber with respect to the aperture can lead into 6%-10% drop in the enhancement ratio. The slight horizontal shift of the peaks is less than the frequency resolution of the scan. Although increasing the time delay might give better resolution, it will add the uncertainty of power changes of the laser source into the measurements. As stated earlier, observed enhancement is due to excitation of Mie resonances. However, the positions of the peaks in the enhanced ratio (also in the scattering cross section of the system) are shifted compared to the positions of the first resonances of the Mie scattering coefficient of a single isolated cylinder. 45 As we explained in Sec. II C, this is due to the existence of the mirror image of the fiber, which results into near-field coupling of electric field of the system resonances and leads to the shift of the resonance.

C. MD-double fiber: Equivalent system and discussion
We have demonstrated that a subwavelength fiber excited by an ED can be modeled in the vicinity of the first radiation peak by an effective ED or MD depending on the relative orientation of the dipole source and the fiber. 2 This happens due to excitation of the electric and magnetic Mie resonances overlapping with the position of the first TE and TM WGMs formed at the cross section of the fiber. 10 It is expected that an MD source also excites Mie resonances (TE and TM WGMs) of the fiber leading into effective induced ED and/or MD modes (depending on the orientation of the source relative to the fiber) similar to what has been observed for transverse oriented ED and MD sources in the vicinity of the nanospheres. 29 In general, the coupling efficiency of an electric and magnetic dipole with TE and TM WGMs is described through p · E and m · H, respectively ( p and m are electric and magnetic moments, respectively, and E and H are the electric and magnetic fields of the mode). 10 Hence depending on the direction of dipole, different coupling and enhancement are observed for different source and fiber orientations. In this section, we investigate and compare the numerical results of an equivalent system, consisting of a magnetic dipole as an excitation source, for the aperture-fiber system. Replacing the aperture by a magnetic dipole and using image theory to eliminate the perfect conducting plane lead into a double-fiber equivalent system shown in Fig. 4(a).
Based on the image theory, the reflection of ED and MD parallel to a perfect conductor is antiparallel and parallel to the source, respectively. To take care of image sources after removing the perfectly conducting plane, we introduce a second fiber, corresponding to a mirror image of the fiber with respect to the conductor plane. For an MD source located in the middle and parallel to the fibers axis (along z), the induced MD inside the fibers will be along the same directions, while for an MD source transverse to the fibers (along y), the induced ED in the both fibers will be antiparallel.
We have designed the THz MD source using a subwavelength loop (7 µm diameter) excited by a current source, Fig. 4(a). The fibers are positioned 9 µm apart resulting into a 1 µm gap between MD and each fiber (negligible compared to the wavelength). The normalized integrated power on a similar square area as the raster scan (850 × 850 µm 2 and 10 µm away from the fiber) for the MD source only and MD-double fiber system (when fibers are along the x and y axes) are shown in Fig. 4(b). The curves for all three cases are normalized to the integrated power of the MD source at 0.5 THz. The results exhibit a similar trend to that of the simulation results for the aperture-fiber system up to 0.3 THz (a/λ = 0.3). As discussed above, beyond 0.3 THz [gray-shaded region in Figs. 4(b) and 4(c)], the approximation of the aperture acting as an MD source does not hold anymore. 38 This is apparent from the comparison of black dashed and marked solid lines. The enhancement ratio for both orientations is shown in Fig. 4(c), where the results are in good agreement for the first peaks. This indicates that the system is equivalent to an MD and double-fiber system. The vertical dashed black lines show the position of Mie resonance of an infinitely long fiber, 2,45 which as mentioned do not overlap with the enhancement points of the system. The peaks observed in the system are shifted due to existence of the mirror image of the fiber leading into coupling of the Mie resonances of the two fibers in the near field. Note that we have observed an 8%-10% decrease in the first peak magnitudes when the power is integrated over a large sphere or infinite plane instead of the square plane of integration. This indicates that the enhancement ratio in Fig. 3(d) can be a reasonable first order approximation of the MD Purcell enhancement of the aperture-fiber or equivalently dipole-double fiber system.

III. CONCLUSION
We have experimentally demonstrated that a subwavelength fiber in front of a sub-wavelength aperture unidirectionally enhances the forward emission more than one order of magnitude. We attribute this enhanced emission to the excitation of Mie-type/WGM resonances formed in the cross section of the microfiber. We have demonstrated that the aperture-fiber system is equivalent of an MD-double fiber system. This confirms that the experimentally measured enhancement ratio is a reasonable approximation of the THz Purcell enhancement of an MD source in the vicinity of subwavelength fibers. Presented results confirm that the approximation of field uniformity in the aperture and representation of the aperture as an MD source is a valid description of the system.
The enhancement of forward emission depends on two factors: loss of the fiber material and refractive index. Lower material loss and higher refractive index will enhance the emission further. An ideal material for THz would be silicon (Si) with lower losses and high refractive index compared to the soft glass in this work. However, the required diameter to observe enhancement due to the first Mie resonances in Si fibers is less than 100 µm for THz, which is extremely difficult to achieve with fiber drawing. It is worth noting that if the fiber is replaced by a sphere of the same diameter, similar enhancement of emission is expected due to excitation of Mie resonances of the sphere. 46 The advantage of the fiber system is that it enables the selection of the resonances.
Such a hybrid aperture-fiber system can be used as a unit cell of metasurfaces for tailoring electromagnetic radiation specifically in the THz part of the spectrum which is less developed compared to other frequency ranges. Moreover, this THz experiment is a proof-of-concept of a possibility to achieve enhanced radiation of an MD source coupled to a nanofiber in optics, which has a potential in developing integrated nanophotonic devices such as integrated nanoantennas on fibers and lasers when the fiber is doped with an active material.

SUPPLEMENTARY MATERIAL
See supplementary material for the time evolution of THz radiation passing through aperture only, aperture and fiber along x, and aperture and fiber along y.