Power circulation via negative energy-flux wormholes in optical nanowaveguides

: We show that energy circulation within a pulse is possible when it propagates in a high-contrast dielectric nanowire. This process is accomplished through electromagnetic “wormhole” regions, in which the Poynting vector associated with the guided mode is negative with respect to the direction of propagation. For demonstration purposes this mechanism is elucidated in AlGaAs and silicon nanowaveguides. The effect of dispersion on the power circulation is also considered.


Introduction
Controlling the flow of light on a nanometer scale may have important implications in both science and technology [1]. Progress in the area of nano-photonics is expected to impact communication and computing technologies as well as bio-photonics. In the last few years we have witnessed rapid advances in the fields of near-field microscopy and spectroscopy and sub-wavelength imaging [2]. The fabrication of high contrast dielectric nanowaveguides is yet another important development in this area. Because of their strong index contrast, these subwavelength waveguides or nanowires can exhibit considerably altered dispersion characteristics and enhanced nonlinearity [3][4][5]. When implemented with semiconductor materials such as GaAs or AlGaAs, these nano-waveguides show promise for efficient second harmonic generation and self-phase modulation at very low optical powers [6][7][8].
In the last few decades or so, the electromagnetic problem associated with wave propagation in cylindrical dielectric waveguides has been intensely investigated in a number of studies [9][10][11]. Given the history of this topic, at a first glance it may seem that all is known regarding the guiding behavior of these structures. An unusual effect associated with the electromagnetic properties of the fundamental 11 HE mode in high-contrast dielectric rods was theoretically found by Gillespie in 1960 [12]. This effect, which has so far received little if any attention, arises from the fact that in such high-contrast structures there can be two regions (close to the core-cladding boundary) where the Poynting vector is negative with respect to the direction of propagation. We emphasize that in all cases (for lossless fibers), the phase velocity and group velocity vectors of the 11 HE mode are parallel, i.e. the negative power never exceeds the positive power and thus no backward waves are associated with the presence of these regions [12]. An important facet of this intriguing effect that still remains unresolved arises when one examines the energy balance equation of the 11 HE mode during pulse propagation, e.g., by considering [13] (1) over a virtual volume V (surrounded by a surface A ) that is totally embedded in one of these regions. In Eq. (1), S is the energy flux or Poynting vector, and W is the total electromagnetic energy in the volume V . More specifically, if a pulse first excites the front face of the virtual surface enclosing this volume (where the energy flux is negative), then how is the escaping energy balanced in this region? To address this issue we will investigate the behavior of this system in the time domain.
In this paper we show, that in high-contrast dielectric optical nanowires, power circulation is possible during pulse propagation via negative energy-flux wormholes. This phenomenon is in fact universal since it occurs in a variety of waveguide geometries as long as the waveguide index-contrast exceeds a critical value. For demonstration purposes this mechanism is examined in circular AlGaAs and silicon nanowaveguides. We theoretically demonstrate that under pulse conditions these wormholes are constantly exchanging energy with their surroundings and at "infinity" (with respect to the pulse center) they become strongly warpedin a way similar to that expected in negative-index waveguides. For completeness dispersion effects are also accounted in our analysis and pertinent examples are provided. Methods to detect the presence of these negative energy-flux regions or wormholes are also suggested.

Analysis
To demonstrate our results we consider wave propagation in a cylindrical optical nanowire of core radius a . This dielectric rod is surrounded by air/vacuum and its refractive index is taken here to be n . The vacuum wavelength of the optical wave exciting this nano-waveguide is 0 λ . In this case, the electromagnetic field associated with fundamental 11 HE mode of this cylindrical structure is given by [14] ( where its propagation constant β can be obtained from the first root of the corresponding eigenvalue equation. In this study, without any loss of generality, the fields of Eq. (2) were judiciously chosen so as the transverse electric field of the of the fundamental mode is primarily x-polarized ( x HE 11 ) [14]. As usual the transverse components ( in the core and cladding regions respectively. Similarly, the magnetic fields are given by: where A,B,C,D, u, and w, are all frequency dependent quantities [14]. The time-averaged Poynting vector can then be obtained from continuous-wave (CW) conditions has only a ẑ -component, i.e., is parallel to the direction of propagation.
As an example, Fig. 1 . Figure 1(b) clearly shows that in this case the intensity distribution in this nanowaveguide is elliptical-like (with an aspect ratio of 80 %). More importantly, at the core-air interface (around the x-axis or 0 = θ ) there are two regions or wormholes where the energy flux density or Poynting vector becomes negative. What is quite interesting, is that, this effect results and coexists with the process of total internal reflection (TIR) which is responsible for the existence of the guided mode. These two areas have a crescent-like shape as shown in Fig. 1(c). In general, these negative z S features arise from the severe disruption of the electromagnetic field lines at the core-cladding interface and are only possible in high-index-contrast structures. It is worth noting that this effect does not occur for TE and TM modes since in these cases one can directly show that z S is everywhere positive. Figure 2 Figure   2 also indicates that these effects can become considerably more pronounced at higher index contrasts. This may be more relevant in the microwave region where materials with very high refractive indices of 10 8 − ≈ n are known to exist [15]. We emphasize that these negative energy flux regions occur in a universal fashion, e.g. they are not uniquely associated with cylindrical waveguides. Figure 3 depicts the electric field distributions as well as the negative Poynting vector regions as obtained using finite element methods for three different waveguide geometries when  nm × Note that these wormhole domains are always enhanced around edges and regions of high geometrical curvature. Thus these effects are a byproduct of high-index contrast and geometrical effects, around which the electromagnetic field lines are severely distorted.

Results and discussion
Even though the CW analysis indicates that such negative energy flux regions are possible in high-index optical material systems the question remains how these wormholes get established in the first place and how is their presence consistent with the energy considerations of Eq. (1). To address these issues we will examine this problem in the time domain. To do so we express both the electric and magnetic field vectors through a Fourier superposition, i.e., where 0 ω is the carrier angular frequency of the wavepacket and 0 β is the corresponding propagation constant of the 11 HE mode.
In Eq. (4), ) , ( T z φ represents the envelope of the pulse which is given by where in Eq. (5) the envelope was expressed in polar form, e.g.
Let us now analyze this latter result. The first term in the bracket corresponds to a modulated ẑ power flow (CW-like part). The second term in Eq. (5) describes an energy redistribution among spectral frequency components (different frequencies see a different 11 HE distribution). From the fact that  Even more importantly, at a certain time-distance (ahead or behind the pulse center) where the power tends to diminish, these wormholes become strongly warped. This warping effect of the power-flow density is schematically illustrated in Fig. 5, in the T x − plane. Our results shed light on how such an energy balance is accomplished around these negative z S zones. More specifically, at the front of the pulse, these wormholes are fed from the positive z S regions (from core and cladding) whereas at the back of the pulse this energy is returned. Even under CW conditions, this same reasoning is still applicable if one considers the time history (pulse-fronts) of this field at ±∞ → T . It is important to note that the mechanism responsible for setting up the power circulation in these nanowires is fundamentally different from the one involved in similar effects in negative-index (metamaterial) waveguides [16]. As we pointed out before these effects arise because of highindex contrast and geometrical effects.
One possible avenue to detect the presence of these wormholes is to observe the effect of negative z S on nanoparticles given that these regions primarily exist in the air/vacuum cladding. Since the longitudinal radiation force component is parallel to the z S vector [17][18][19], then particles located within the wormholes are expected to move towards the source. This is in contrast to what happens in similar configurations using low contrast waveguide structures (particles are repelled away from the source). From the previous discussion the direction of the force is expected to be polarization dependent. Other methods to detect this "negative wind" from dielectric waveguides are also currently considered.

Conclusion
In conclusion we have shown that that power circulation is possible during pulse propagation in high-contrast optical nanowires. This process is accomplished via negative energy-flux wormholes that are constantly exchanging energy with their surroundings. Our results explain how these zones are formed and provide a detailed account as to their energy balance. The implications of our results on other aspects concerning wave propagation in optical nanowires are currently investigated.