Advances in Mathematical Physics

Volume 2015, Article ID 281472, 10 pages

http://dx.doi.org/10.1155/2015/281472

## Leaky Modes of Waveguides as a Classical Optics Analogy of Quantum Resonances

^{1}UPIITA, Instituto Politecnico Nacional, Avenida IPN No. 2580, Colonia La Laguna Ticoman, C.P. 07340 Mexico D.F., Mexico^{2}Physics Department, Cinvestav, A.P. 14-740, 07000 Mexico D.F., Mexico

Received 14 April 2015; Accepted 25 June 2015

Academic Editor: Hagen Neidhardt

Copyright © 2015 Sara Cruz y Cruz and Oscar Rosas-Ortiz. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A classical optics waveguide structure is proposed to simulate resonances of short range one-dimensional potentials in quantum mechanics. The analogy is based on the well-known resemblance between the guided and radiation modes of a waveguide with the bound and scattering states of a quantum well. As resonances are scattering states that spend some time in the zone of influence of the scatterer, we associate them with the leaky modes of a waveguide, the latter characterized by suffering attenuation in the direction of propagation but increasing exponentially in the transverse directions. The resemblance is complete because resonances (leaky modes) can be interpreted as bound states (guided modes) with definite lifetime (longitudinal shift). As an immediate application we calculate the leaky modes (resonances) associated with a dielectric homogeneous slab (square well potential) and show that these modes are attenuated as they propagate.

#### 1. Introduction

In general, a* resonance* may be defined as the excitation of a system that results when one of its characteristic frequencies matches a particular value that is defined by either the boundary conditions, the external forces, or any other interaction or constraint applied to the system. The concept arises from the study of oscillating systems in classical mechanics and extends its applications to physical theories like electromagnetism, optics, acoustics, and quantum mechanics, among others [1]. In quantum and nuclear physics the resonances are associated with metastable states of a system which has sufficient energy to break up into two or more subsystems [2]. The whole range of scattering experiments provides a big number of examples in this matter [3, 4]. This last is because scattering includes the situation in which the impinging particle is “captured” for a while in the zone of influence of the scatterer, so that the projectile and the scatterer form a new unstable system (resonance state). The capturing occurs for specific energies of the projectile that are defined by the general properties of the scatterer. Eventually, the projectile is released and escapes from the interaction zone (decay process). In this model the “time of capture” corresponds to the lifetime of the decaying system that is formed of the scatterer and the projectile under the resonance condition. Collisions, in addition, are modeled as interactions localized in time and space. This implies that the involved potential vanishes rapidly enough in space, so that incoming and outgoing asymptotic states can be represented by wave packets in free motion. The problem is usually reduced to the analysis of one-dimensional effective potentials such that the main information is obtained from the transmission and reflection amplitudes of the scattering states; see, for example, [5]. In the case of one-dimensional square potential wells and barriers the resonances appear for specific values of the parameters that define the interaction. Indeed, for such parameters the resonances can be associated with the bell-shaped peaks of the corresponding transmission coefficient [6, 7]. The presence of resonances can be also stimulated either by shallowing a square well [8] or by adding to it a static electric field [9]. As resonances are special cases of scattering states, they are represented by irregular (not normalizable) vectors [10, 11] that satisfy the purely outgoing condition [12] (see also [1]). This last property of resonances is useful in constructing complex potentials that operate as optical devices with position-dependent complex refractive index [6, 13, 14].

Remarkably, quantum wells and barriers can be made with high precision and compositional control by modern epitaxial growth techniques [15]. This has opened a wide window of experimental possibilities to explore various quantum phenomena in real devices and actual laboratory settings. For example, it has been shown that the transmission time of a wave packet that is tunnelling through a one-dimensional barrier does not depend on the barrier thickness [16]; a similar property is found in one-dimensional potential wells [17, 18]. This condition implies anomalously large group velocities of the wave packets that impinge the barrier (well) since the velocity of the packet’s peak must increase with the width across the interaction zone. Such a superluminal phenomenon has been observed using electromagnetic analogues for evanescent modes [19] and microwave pulses [20]; at the same time this result has stimulated the designing of high-speed devices based on the tunneling properties of semiconductors (see, e.g., Chapters 11 and 12 of [21]). On the other hand, as the scattering by quantum wells attenuates the outgoing wave packets only because of the multiple reflections at the well boundaries, it is expected that the resonance condition implies a maximum in the time spent by the projectile in traversing the zone of influence of the scatterer [6]. Recent theoretical research includes the case of a rectangular well that is embedded in an environment formed of a zero potential energy (flat potential) at the right and a parabolic potential at the left of the well [22–24]. Thus, the connection between time delay and “time of capture” associated with the scattering processes would be verified in the laboratory by using semiconductor materials.

Another accessible laboratory tool to study quantum phenomena is provided by waveguide structures. The connection is based on the analogies between light waves and Schrödinger wave-state functions that have been studied since the early years of quantum physics [25]. In this approach the amplitude of an electromagnetic field is associated with the probability amplitude of a quantum state with two degrees of freedom while the propagation of such field through the waveguide corresponds to the time evolution of the quantum state. Such optical-quantum analogy is representative of a wider scheme that includes diverse connections between classical and quantum systems [26] and has been successfully applied in studying Anderson localization [27], coherent enhancement [28] and destruction of tunneling [29], coherent population transfer [30], and decay of metastable states coupled to semi-infinite, tight-binding lattices [31], among others. Of special interest for the purposes of this paper, nanoseconds delay-time measurements have been performed into the electromagnetic framework by using microwave setups of narrowed waveguides [32]. The relevance of this last result obeys the fact that the waveguide simulates a quantum potential barrier so that the measured electromagnetic delay-time is in correspondence with the quantum tunneling time. Thus, the waveguide structures represent an emergent and versatile arena to perform realistic experiments addressed to get deeper insights into the understanding of resonances in quantum physics.

The aim of this paper is to get a classical optics analogy of the resonances associated with one-dimensional short range potentials in quantum physics. The centerpiece of our approach is the connection between the mathematical structure of the propagation of electromagnetic signals in waveguides and the dynamics of a time-dependent quantum state. In this form our interest is twofold: we want to connect the theoretical calculations of resonances with measurable properties of an electromagnetic wave that propagates in a waveguide and we would like to extend the application of the resonance mathematical techniques to the solving of problems associated with electromagnetic waveguides. A collateral result will be the connection between the resonance properties of electrons that are injected into heterostructures and the properties of the electromagnetic signals that propagate in optical waveguides. The paper is organized as follows. In Section 2 we revisit the main properties of the resonance states associated with one-dimensional short range potentials. In particular, the relevance of the Fock-Breit-Wigner distribution in a scattering process is shown. When the transition amplitude from an initial scattering state to its time-evolved version is ruled by such a distribution, one immediately identifies energy resonances between the projectile and the scatterer. It is also shown that analytical continuation of the eigenvalue equation to a complex plane of energies is required to get an appropriate mathematical description of resonance states. In Section 3 we include some of the well-known generalities of the mode fields in waveguide structures and their connection with the Schrödinger equation of a system of two degrees of freedom. The case of electromagnetic propagation through a finite, dielectric homogeneous material and its association with a one-dimensional square well potential is considered in Section 3.1 as immediate application. Section 4 is devoted to the connection between resonance states and a singular form of modal field in waveguides; the latter is known as* leaky mode* and is characterized by suffering attenuation in the direction of propagation but increasing exponentially in the transverse directions. We apply the resonance approach of [6] in the calculation of the leaky modes associated with the dielectric homogeneous slab described in Section 3.1. As an example we provide the* propagation constant* values (energy resonances) corresponding to the leaky modes supported by a slab of definite width and refractive index, and show that these modes are attenuated as they propagate. Section 5 includes the calculation of the longitudinal shifts for the leaky modes of Section 4; we show that the maxima of these shifts are in correspondence with the maxima of the transmission coefficient, a result that verifies the connection between times of capture and resonances in the scattering processes. Finally, Section 6 includes some of the main conclusions of the present work.

#### 2. One-Dimensional Scattering States and Resonances

Let be the Hamiltonian defined by a short range one-dimensional potential . Assuming that the spectrum of is formed of discrete and continuous eigenvalues, one can consider the following sets of state vectors: where and are, respectively, the identity operator in and . Consistently, the spectrum and Hilbert space of are and , respectively. As usual, the discrete and continuous eigenvectors (1) are orthonormal in their respective subspace according to the rules Hereafter we pay detailed attention to the scattering states (i.e., vectors in ). The description of bound states (i.e., vectors in ) is straightforward by using the well-known relationship . The reason is that, as we have said in the introduction, resonance states arise quite naturally in the scattering processes.

Any scattering state can be expressed in the basis of continuous eigenvectors and is such that In -representation the above expressions read as with standing for the complex conjugate of . The above expressions allow writing the energy distribution as follows: Let (3), equivalently (5), be the initial state of a stationary system. The transition amplitude from to is given by the inner product It is clear that can be investigated in terms of spatial coordinates or as a function of the energy distribution (7). In this work we will assume that is given by the following expression: where is the energy of the initial state and is a parameter having units of energy. Note that the bell-shaped function (9) is centered at , for which , and has a half-width at half-maximum equal to . In addition, is a meromorphic function that has two poles . This function is known as either Cauchy (mathematics), Lorentz (statistical physics), or Fock-Breit-Wigner (FBW) (nuclear and particle physics) distribution. Using (7) we can write the Fourier coefficient in the expansion (3) as follows: After the analytic continuation of to the complex -plane one can introduce (9) into (8) to arrive, up to a global constant factor, at the expression This last result has the form of a transient oscillation [1] and means that the transition from to is an exponential decreasing function of time whenever is given by (10). That is, the initial state decays according to the time-dependent probability As we can see, the width of the distribution (9) represents an indirect measure of the lifetime of because at time the probability that this state has not yet decayed is reduced approximately in , that is, . It is clear that and are correlated; the smaller the value of the larger the lifetime . In general, one says that there is a resonance of width when either the energy distribution of the system at time is given by the FBW function (9) or the Fourier coefficients of the expansion (3) are given by (10).

We would like to emphasize that the dependence of , , and on the complex number , as this has been expressed in (9)–(12), is not merely aesthetic; this corresponds to the formulation developed by Gamow, Fock, Breit, and Wigner in the late 1920s and along the 1930 decade (see details in [1]). It is due to Siegert [12] that a resonance is introduced as a solution of the Schrödinger equation associated with the complex eigenvalue and satisfying the purely outgoing condition [6]: where and It can be shown that the above conditions lead automatically to the FBW distribution (9) with a number in the fourth quadrant of the complex plane ; see [1]. In turn, the distribution of the complex points characterizes a Riemman surface of two sheets, each one cutting along the positive real axis. These complex numbers define directly the behavior of the propagator (12). For instance, consider a very narrow distribution , that is, (the point is close to the positive real axis). As defines the vicinity , we see that “large lifetimes” means narrow widths , so that the energy of the system at time is very close to the initial energy . In the extremal case , the distribution is as narrow and singular as . That is, the related lifetime is large enough to consider the initial vector a nondecaying (stable) state. In this context, bound states are considered as decaying states with infinitely large lifetimes and vice versa, unstable (decaying) states correspond to bound states with a definite lifetime.

To get more insights into the complex eigenvalue let us consider the Schrödinger equation Assuming separation of variables , with and the separation constant, the probability density does not depend on time . Now, let us make the analytic continuation of the eigenvalue equation in order to include complex eigenvalues , that is, . In this case we have , and the probability density is an exponential decreasing function of time . Clearly, the spatial part of the wave-function is no longer a square-integrable function because this belongs to the complex eigenvalue and the Hamiltonian is still Hermitian. Indeed, it can be proved that diverges exponentially as . Thus, the density increases exponentially for either large or large negative values of . The usual form to avoid some of the complications connected with the limit is to consider the long lifetime limit ; see, for example, [33]. Therefore, we consider the condition This last means that the level width must be much smaller than the level spacing in such a way that closer resonances must imply narrower widths (longer lifetimes). In general, the main difficulty is precisely to find the adequate and . However, it has been shown [6] that for one-dimensional stationary short range potentials the superposition of a denumerable set of FBW distributions (each one centered at each resonance , ) entails an approximation of the transmission coefficient such that the larger the number of resonances involved, the higher the precision of the approximation:

#### 3. Mode Fields in Waveguide Structures

Consider the Helmholtz equation with the wave number in vacuum. Equation (18) can be achieved from the Maxwell equations for transversal electric fields in presence of inhomogeneous dielectric materials, neglecting the magnetic polarizability and the dispersive properties of the medium. Indeed, (18) corresponds to an electric wave propagating along the positive -direction and polarized in the -direction, assuming that the refractive index depends only on the -coordinate. In the paraxial regime we use the ansatz to transform (18) into the linear second-order differential equation known as the paraxial Helmholtz equation. In (19)-(20) the number is a reference refractive index that may be taken as the maximum of , and we have assumed that the weakly guiding condition () holds. In this approach the parameter , known as the propagation constant, is the value of the -component of the generalized linear momentum [34] and defines the slope of the light beam as follows:

In general, modal fields associated with open waveguides are classified into guided (trapped) modes and radiation modes, the latter propagating outside the core [35–38]. The guided modes are bounded fields that correspond to real values of the propagation constant fulfilling . They are perfectly guided through the waveguide as they suffer no attenuation along the optical axis and decay exponentially in the transverse direction. Thus, as guided modes are trapped into the core of the waveguide they add no energy flux into the infinite surrounding clad. In turn, propagating radiation modes correspond to real propagation constants fulfilling , so that they are oscillatory fields in the transverse as well as in the longitudinal directions [35]. They appear as a consequence of the infinite extension of the clad cross section. Guided and radiation modes constitute a complete set of modal fields in the sense that any physical solution to the paraxial Helmholtz equation can be expressed as a linear superposition of a finite number of guided modes and a continuum of radiation modes.

At this stage it is important to stress the resemblance of the electric field expression (19) with the solutions of the Schrödinger equation (15). This allows the identification In adittion, (20) would correspond to the eigenvalue equation , provided the following additional identification holds: The identification (22)-(23) is well known in the literature; see, for example, [34]. This allows to construct the classical optics analogy of diverse quantum phenomena and vice versa, the quantum analogy of classical optical phenomena [26, 34]. The analogy includes the sets of vectors introduced in (1) because bound states correspond to “trapped” waves that have no chance of escaping from the zone of influence of the scatterer (as we have seen in the previous section, these states can be interpreted as decaying states with infinitely large lifetime). Thus, quantum bound states correspond to guided electromagnetic modes. Particles in scattering states, on the other hand, can be temporally trapped by the scatterer in such a form that their presence in the zone of interaction decreases in time. In this form, scattering states correspond to radiation modes. The orthonormalization properties of guided and radiation modal fields can be stated in a form similar to that of the quantum mechanical bound and scattering states (see, e.g., [35, 37]). Hence, the total radiated power can be expressed as the sum of two terms: the first one associated with the power that propagates inside the core, given by a finite sum containing the contribution of guided modes, and a second term corresponding to the power radiated from the core, expressed as an integral over all radiation modes. In this form, denoting by and the guided and the radiation mode amplitudes, respectively, the total power of the radiated field in the waveguide axis direction is given by [36] On the other hand, as we have seen, resonances are scattering states of complex eigenenergy that spend definite intervals of time in the interaction zone. So resonances are not connected to radiation modes because these last do not belong to complex eigenvalues. Next we will take full advantage of the analogy indicated above to get optical models for the resonances discussed in Section 2.

##### 3.1. Wave Propagation through a Dielectric Homogeneous Slab

As an immediate application let us analyze the simplest case of electromagnetic propagation through a finite, dielectric homogeneous material. It is convenient to rewrite the eigenvalue equation (20) as follows: Considering a combination of media distributed along the -axis such that there is vacuum in and a medium of constant refractive index in , (25) is in correspondence with the eigenvalue equation defined by a linear square well of depth and width . The solutions to this equation are well known and can be found elsewhere, though the analysis of resonances and the resolution of the transmission coefficient in terms of the FBW distribution (see (17)) can be consulted in [6]. Hereafter we will label the regions of the media (vacuum, medium, vacuum) from left to right as , , and .

#### 4. Leaky Modes and Resonances

Besides the guided and radiation modes discussed in the previous section, it is well known the existence of solutions to the Helmholtz equation with complex propagation constants that have an “unusual behavior”: they increase exponentially in the transverse direction and suffer attenuation in the longitudinal one [35]. These solutions are called leaky modes and have attracted much attention over the time. The leaky modes do not belong to the set of guided and radiation modes of infinite transverse waveguides, though they constitute a discrete set of modes that can be used to approximate the continuous mode contribution to the radiation field, in regions close to the core and far from the source, in order to avoid integration overall radiation modes [36–40]. Yet, leaky modes have properties that are very similar to the ones of guided modes, but the former modes loss energy as they propagate along the optical axis [37]. In slab waveguides this leakage is due to refraction at the waveguide core-clad boundaries (refracting leaky modes [36]). In waveguides with curved cross section the leakage can be also produced by the partial reflection at the boundaries due only to the curvature of the cross section (tunneling leaky modes) [35, 36]. It can be shown that the leakage is slower in the tunneling modes compared to the refractive modes [35]. In any case, these modes behave much like guided modes as they are able to propagate large distances inside the guide before their power loss started to be significant. Different methods have been applied to determine the leaky modes of planar [37, 41] as well as cylindric [35, 39] waveguides involving isotropic and anisotropic dielectrics [42]. These techniques include numerical methods as mode-matching [43], finite differences and finite element [38, 42, 44], approximate methods as WKB [45], transfer matrix [41] and differential transfer matrix methods [40, 46]. In contrast to the diversity of mathematical techniques to determine the propagation constants of leaky modes there is still a conceptual discussion concerning them, the exponential growing in the transversal direction of these nearly guided modes make them appear as nonphysical solutions to the Helmholtz equation. Yet, this behavior is required if the energy flux must be preserved: as the flux is damped in the longitudinal direction, there must be an increasing flux in the transverse direction [37].

The problem of finding leaky modes is mathematically equivalent to determine the resonant states in quantum mechanics for a particle in a short range interaction potential. This is convenient because we may apply different methods that have been addressed to obtain the resonance energies and the corresponding wave functions in order to determine the complex propagation constants and the modal fields. To this end, let us regard the square well potential analogy indicated in Section 3.1 with a complex propagation constant, . The corresponding complex wave number is related to throughThe energy density flux in the transverse direction is not null because (see [6]). A wave coming from region enters into resonance with the optical medium and is thus trapped inside the core (region ); then it propagates and travels a finite distance along the optical axis. After a while, the wave is finally emitted into the clad. The result is a wave propagating to the left in region and another one propagating to the right in region . The correct sign for is thus positive and the whole process obeys the pure outgoing condition. The expression for the electric energy density fulfillsmeaning that the power decreases exponentially as grows whenever and . The optical-quantum analogy requires also a purely outgoing wave (see Section 2) so the amplitudes of the external electric fields increase exponentially as . Leaky modes then correspond to wave numbers in the fourth quadrant of the complex plane. The factor can be interpreted as the mean distance traveled by the beam, along the optical axis, before its intensity decays by of its intensity at a focus. In this context represents a measure of the transient confinement of the field into the core. In the limit we will have a quasistationary mode propagating inside the slab in a very similar way as a guided mode. In this limit, it is possible to establish an analytic approach to calculate the propagation constant. Bearing this in mind, following (16), we may writewhere is the separation between two consecutive resonances. To get the position and the width of the resonances we follow the approximations reported in [6]. After some calculations one arrives at the expressionswith an integer number fulfillingThus, for each combination of the parameters , there is a finite number of resonances given by (29). Table 1 shows the propagation constants corresponding to the leaky modes supported by a slab waveguide of refractive index and width . The behavior of the real part of the field amplitude of these modes, as they propagate in the longitudinal direction through the waveguide, is depicted in Figure 1 for (left), (center), and (right). Observe that the modes with smaller values of attenuate slower compared to those with larger values of .