1. Introduction

When light is incident from a high-refractive-index medium onto a low-index medium, it undergoes total internal reflection provided the angle of incidence is larger than a certain critical angle. Total internal reflection is a fundamental physical phenomenon with several famous applications; in particular modern telecommunications rely on optical fibers based on this phenomenon.

Since the tangential electric and magnetic fields must be continuous at the interface, there must be nonzero fields in the low-index medium, even though the incident wave is totally reflected. For lossless/gainless media, these evanescent fields decrease exponentially away from the interface. The presence of evanescent fields in the low-index medium suggests that the reflected wave will sense any perturbation induced there. In particular, if the low-index medium has gain, the reflection response will change compared to the lossless/gainless case. The problem of determining the correct electromagnetic response in the case of an active low-index medium is far from trivial, and has been discussed for 40 years without reaching consensus [1–11]. A key issue is whether the reflectivity may exceed unity (i.e., evanescent gain exists) when the active medium fills the entire half-space. Experiments have indicated that evanescent gain exists [12–15]. However, it has been argued that the amplified reflection may be due to backreflection from e.g., the boundaries of the active medium [11].

When the active medium has a finite thickness, it is well known that the overall reflection from the slab may exceed unity. This situation is fairly simple, as there is no need to determine the sign of the longitudinal wavenumber in solving Maxwell’s equations for this case; the two waves (with opposite signs) are present simultaneously.

Since there are no gain media with infinite thickness, why examine this case? The answer becomes clear if we formulate a similar question in terms of the refractive index: Why define the refractive index as a separate parameter, when the electromagnetic field in any realistic, bounded structure can be expressed in terms of the permeability and permittivity? While the refractive index or longitudinal wavenumber are not needed to obtain the formal solution to Maxwell’s equation in a finite slab, it is still useful since it immediately provides information about the involved physics. For example, it predicts whether the medium refracts positively or negatively [16]. Also, assuming darkness for time t < 0, the solution to Maxwell’s equations for a semi-infinite gain medium equals that of a finite slab for times t less than d/c, where d is the slab thickness and c is the vacuum velocity of light. Hence, understanding semi-infinite media helps explaining transient phenomena.

We will now summarize the existing controversy. Assuming well defined frequency-domain fields, Maxwell’s equations can be solved in the frequency-domain, using the sign convention exp(−iωt). With respect to Fig. 1 we define the transverse wavenumber (spatial frequency of the source) k_x. For simplicity we assume both media to be nonmagnetic. Let ɛ₁ and ɛ₂ be the relative permittivities of the high-index medium to the left and the low-index medium to the right, respectively. For plane waves, Maxwell’s equations require the longitudinal wavenumbers in the high-index and low-index media to be

 

Fig. 1 A wave is incident from a high-index medium to a low index medium with gain. The source produces a single, spatial frequency k_x. The electromagnetic boundary conditions require preservation of the wavenumber k_x parallel to the interface. The longitudinal wavenumbers are denoted k_1z and k_2z. Note that since the excitation is assumed to be causal, it contains a band of frequencies, and therefore also a band of k_1z’s and k_2z’s.

Download Full Size | PDF

 

Fig. 2 The two possible solutions for the wavenumber k_2z for monochromatic analysis and a gainy medium. The arrows indicate the two possible wavenumbers in the complex plane, as the gain tends to zero. For a lossy medium, we always have a solution that tends to the upper alternative $+i\sqrt{k_x^2-{\omega }_1^2/c^2}$ in the limit of zero loss. For simplicity we have taken Reɛ₂ = 1 here.

Download Full Size | PDF

In this work we will first go back to fundamental electromagnetics, to ensure that we use the principle of causality in its most primitive form: No signal can propagate faster than the vacuum velocity of light. After the general analysis in Section 2, we consider conventional, weak gain media in Section 3 and show that they provide evanescent gain. In Section 4 we present an example that demonstrates that not all gain media give evanescent gain; this depends on the medium’s global dispersion behavior.

2. Laplace transform frequency-domain analysis

Going back to fundamental electromagnetics, we note that Maxwell’s equations, combined with appropriate causal constitutive relations, contain everything necessary to obtain a unique solution. To determine the correct solution, we must be certain that we consider the real, physical situation. The real physical fields are the ones in the time-domain. By requiring the fields to be zero for t < 0 (see Appendix B), we obtain the causal solution to Maxwell’s equations. The complex frequency-domain fields are usually found from the time-domain fields by a Fourier transform. However, when there is gain in the system, using the Fourier transform can be perilous, since the field may increase with time. At first sight, any instability seems to be convective in our case. This is however not true: A causal excitation involves an infinite band of frequencies. For a single spatial frequency k_x this means that modes with a wide range of incident angles are involved; in fact even the mode with k_2z = 0 may be excited. This “side wave” gets amplified and leads to infinite fields at the boundary. This instability is somewhat artificial, since its existence is dependent on infiniteness in the transverse direction; we will argue below how it can be ignored in certain situations. Nevertheless, within a linear medium framework, Fourier transforms do not necessarily exist. Therefore, as in electronics and control engineering, we generalize the analysis by using the Laplace transform,

In order to derive the Fresnel equations and determine the sign of k_2z, it is tempting to start with the response from a slab of finite thickness d, and then take the limit d → ∞. For finite d the solution to Maxwell’s equations is independent of the sign of k_2z in the slab [17, 18]. However, for an active slab the multiple reflections may diverge, especially for a large d. Thus, for real frequencies, the limit d → ∞ is not necessarily meaningful [11, 17]. A way around this, is to evaluate the fields for sufficiently large Imω, where the frequency-domain fields exist. There, an exponential increase is quenched by the exponential factor exp(−Imωt). As a result, we can take the limit d → ∞ [16]. For TE polarization, the Fresnel reflection coefficient ρ and the transmission coefficient τ (including the propagation factor exp(ik_2zz)) become [16, 17]

3. Weak gain media

To find the actual reflection and transmission response, we first consider conventional weak gain media, with the following assumptions or properties:

Properties 2 and 4 essentially mean that the gain is weak and the dispersion is small.

We now solve the equation

We therefore examine the region around k_xc in more detail. If there were two solutions ω_a and ω_b to the dispersion relation, then Eq. (8) would predict that

In the expression for k_2z and the Fresnel coefficients (4), these solutions appear as branch points. Hence, when evaluating the physical time-domain fields by the inverse Laplace transform, we must integrate above the associated branch cuts, from −∞ + iγ to +∞ + iγ, see Fig. 3. By path deformation this path is the same as the path from −∞ to ∞ plus the paths around the branch cuts in the upper half-plane (Fig. 3). Thus we may use the inverse Fourier transform to determine the time-domain fields, but only if we add the integrals around the branch cuts. Due to the exponential factor exp(−iωt), the integrals around the branch cuts diverge and dominate after some time.

 

Fig. 3 The complex ω-plane. For conventional, weak gain media, there are branch points right above ω = ±(k_xc +δ′) ≈ ±k_xc. The branch cuts can be chosen arbitrarily; however, the shown, vertical cuts minimize the integral around the part of the branch cuts in the upper half-plane.

Download Full Size | PDF

The divergence of the time-domain fields can be explained as follows. Any causal excitation involves an infinite frequency band. For example, the Laplace transform of a unit-step-function modulated cosine, u(t)cos(ω₁t), is $i\omega /({\omega }^2-{\omega }_1^2)$. Thus, it is nonzero for all finite ω ≠ 0. One of these frequencies is the branch-point frequency for which k_2z = 0, that is, ω ≈ k_xc + iδ. This frequency is complex; the imaginary part δ means that the associated eigenmode is a growing wave with envelope exp(δt). Physically, a wave with k_2z = 0 propagates along the boundary. Because the medium is gainy, this side wave picks up gain on its way. Consider a fixed observation point, e.g. the point z = 0⁺ and x = 0. Since the medium and the excitation are unbounded in the transverse x-direction, there are side waves that start arbitrarily far away from the observation point. Thus the field at the observation point diverges. As the field in medium 2 becomes infinite, the field in medium 1 is infinite as well. Since the field at a fixed point in space diverges and the instability is not a result of amplified, multiple reflections, the instability for the system in Fig. 1 can be classified as an absolute instability [17, 20, 21].

This instability could be eliminated (or converted into a convective instability) by limiting the extent of the gain medium in the transverse direction with an absorbing boundary. Alternatively, the incident wave itself could be limited in the x-direction, leading to an infinite spectrum of k_x modes (see Appendix A and Ref. [8]). Rather than imposing such remedies, we will simply calculate the time-domain fields by an inverse Laplace transform above the branch cuts. If the excitation frequency ω₁ is sufficiently remote from the branch points, the side wave with k_2z = 0 is only excited very weakly, and can be neglected up to a certain time. The condition that the excitation frequency is remote from k_xc means that the incident angle is not close to the critical angle. This condition is imperative in order to distinguish between the reflected wave, with an angle of reflection equal to the angle of incidence, and the wave associated with the growing side wave, with “reflection” (or propagation) angle equal to the critical angle.

The reflected time-domain field for the excitation u(t)exp(ik_xx – iω₁t), with Laplace transform exp(ik_xx)/(iω₁ – iω), is given by

We can now answer the question about the existence of evanescent gain. To obtain Eq. (13), we have only considered two branch cuts in the upper half-plane; these are the necessary branch cuts due to the zeros of ${\varepsilon }_2(\omega ){\omega }^2/c^2-k_x^2$. We must ensure that the integrand in Eq. (12) is analytic everywhere else in the upper half-plane. That is, the sign of k_2z must be determined such that k_2z is analytic everywhere, except at the two branch cuts in the upper half-plane. Since k_2z → +ω/c as ω → +∞, we can determine the sign by decreasing ω from +∞ to ω₁, ensuring that k_2z is continuous everywhere except at ω = k_xc where it changes sign. From Fig. 4 we find that Imk_2z > 0 at the observation frequency ω₁. Hence, for weak conventional gain media, provided the “reflected” field from the side wave can be ignored, evanescent gain is possible. This result is consistent with [7,8], and the time-domain simulations in [10] where the dispersion of the medium is discarded.

 

Fig. 4 The function $k_{2z}^2(\omega )={\varepsilon }_2(\omega ){\omega }^2/c^2-k_x^2$ for a typical gain medium, plotted in the complex $k_{2z}^2$-plane. To identify k_2z, we require it to be +ω/c at ω = ∞, continuous as ω decreases towards zero, except at the branch cut at ω = k_xc where it changes sign.

Download Full Size | PDF

In Fig. 5 we plot the reflected and transmitted electric field for a weak Lorentzian medium, after the transients have died out, and before the side wave dominates. The reflected field was computed by Eq. (12), including the propagation factor exp(−ik_1zz) in the integral. The transmitted field was computed with the same equation, but with τ instead of ρexp(−ik_1zz) in the integral (see Eq. (4)). For z > 0 we clearly see an evanescent decaying field, while the reflected field for z < 0 is larger than unity.

 

Fig. 5 The reflected electric field (solid line, z < 0) and transmitted electric field (dashed line, z > 0) for a plane wave incident to a weak Lorentzian gain medium (20). The parameters used: F = 0.01, Γ = 0.1ω₀, k_xc = 2ω₀, ω₁ = ω₀, γ = 0.001, and ɛ₁ = 4.7. The field is plotted for x = 0 and ω₀t = 10³, and normalized to the incident field. The amplitude of the reflected field is 1.01. Note that the field is discontinuous at z = 0 because the incident wave is not included.

Download Full Size | PDF

It is interesting to examine the situation when we approach the critical angle associated with the frequency ω₁. If we insist on using only the first term of Eq. (13) in this case, a simple calculation shows that the power reflectance would have been bounded by ${(\sqrt{2}+1)}^2\approx 5.83$ at the critical angle. Also, the wavenumber k_2z and the reflected field would be discontinuous as we pass the critical angle. This is clearly a paradox, as the branch cuts were chosen arbitrarily. The dilemma is resolved by noting that the entire Eq. (13) must be used in this domain; both terms naturally coexist and cannot be separated. As we approach the critical angle, ℰ_bc(0,t) becomes comparable to or larger than the first term in Eq. (13), for all times. A different choice of branch cuts will alter each of these contributions, but the sum remains the same. For finite transverse dimension, the side wave’s contribution to the “reflected” field does not necessarily diverge any more; however, the intensity of the reflected field can be arbitrarily large as the dimension is increased, or if the reflections from the transverse end facets are large.

4. General gain media

More sophisticated gain media can be constructed, at least in principle, that behave differently compared to the conventional weak gain media. We will here show that we can obtain a near-imaginary k_2z with negative imaginary part at an observation frequency ω₁ < k_xc. Consider the permittivity

 

Fig. 6 The real and imaginary parts of $k_{2z}(\omega )=\sqrt{{\varepsilon }_2{\omega }^2/c^2-k_x^2}$, with ɛ₂ given by Eq. (15). We have set ω₂ = k_xc, N = k_xc(6/10 – i/1000) and P = k_xc(7/10 – i/1000).

Download Full Size | PDF

Any realistic incident wave contains a spectrum of wavenumbers k_x. While there are no zeros of $k_{2z}^2$ in the upper half-plane for the particular k_x considered above, this is not the case for all possible k_x. Thus, also for this medium there are growing waves. The fact that the medium has large gain, and the presence of instabilities, mean that it is very challenging to observe the “anti-evanescent” response in practice. In principle, however, up to a certain time the amplitude of the instabilities can be limited by ensuring a narrowbanded spectrum of incident k_x’s. Formally, if σ is the width of the incident wave, and ℰ_σ(x,z,t) is the resulting electric field, lim_σ→∞ℰ_σ(x,z,t) tends to the “anti-evanescent” response as t → ∞, while lim_t→∞ℰ_σ(x,z,t) = ∞ for any finite σ.

The permittivity (15) has a double pole at ω = 0. While the medium is causal in principle, the medium might be easier to realize if the pole is moved slightly away from the origin, into the lower half-plane. It turns out that this modification does not alter the permittivity function signifiantly, in the frequency range of interest. Also, if desired, the behavior at ω = ∞ can be adjusted along the lines described in Ref. [17].

In Fig. 7 we plot the reflected and transmitted field for the gain medium (15), calculated with the inverse Laplace transform for a sufficiently large t when the transients have died out. Only a single k_x has been excited. The reflection amplitude is 0.98, and the transmitted field is an exponentially increasing function of z. While realizable in principle, the example is highly unrealistic: To observe a behavior similar to that in Fig. 7, t must be at least of the order of 10²(k_xc)⁻¹; otherwise the transients would disturb the picture. Any realistic gain medium has finite thickness. However, to act as a semi-infinite medium, the thickness d of the gain medium must satisfy d > ct, or k_xd ≳ 10², such that the light has not reached the back end. With the “anti-evanescent” growth rate in Fig. 7, this would imply unphysically large fields (or in practice, nonlinear gain saturation). Hence, if the “anti-evanescent” behavior is to be observed experimentally, one would need to construct a medium where the transients die out rapidly, and/or a medium which leads to a sufficiently small |Imk_2z|. At the same time the medium must violate the conditions in Sec. 3; that is, it must have large gain and/or large dispersion for some frequencies.

 

Fig. 7 The reflected electric field (solid line, z < 0) and transmitted electric field (dashed line, z > 0) for a plane wave incident to the gain medium described by Eq. (15) and Fig. 6. The field is plotted for x = 0 and k_xct = 10⁵, and normalized to the incident field. The amplitude of the reflected field is 0.98. The parameters used: ω₁ = 0.65k_xc, γ = 10⁻⁶, and ɛ₁ = 4.

Download Full Size | PDF

5. Conclusion

We have considered the case where light is incident from a high-index medium to a low-index medium with gain, generalizing the situation with total internal reflection.

In principle, it is apparent that both solutions (k_2z in the second and fourth quadrant of the complex plane) can be attained with a suitably engineered medium. In other words, evanescent gain may or may not be the case, dependent on the detailed permittivity function. This demonstrates the fact that the sign of k_2z cannot be determined from the electromagnetic parameters at a single frequency, but must be identified from the entire frequency domain dependence, after a check of possible non-analytic points (instabilities) in the upper half-plane of complex frequency.

For conventional, weak gain media, we have seen that there is an absolute instability associated with infinite transverse dimensions. In some cases this instability can be eliminated or ignored; then evanescent gain prevails.

A. Finite incident beam and finite size medium

To origin from a realistic source, an incident beam should not only be causal, but also be of finite width. We will here describe how to model an incident beam using standard Fourier optics, and argue that even for active media, we are allowed to interchange the order of integration with respect to transverse wavenumber k_x and frequency ω. Thus we can treat a causal excitation of each k_x separately.

Let ℰ (x,t) be the incident TE field at the interface between the high-index medium and the active low index medium. Performing a Laplace transform t → ω followed by a Fourier transform x → k_x, we obtain the transformed field E(k_x,ω). The inverse transform is given by

(18)$$\begin{array}{lll}ℰ(x,t)\hfill & =\hfill & \frac{1}{{(2\pi )}^2}{\int }_{-\infty +i\gamma }^{\infty +i\gamma }\text{d}\omega e^{-i\omega t}{\int }_{-\infty }^{\infty }\text{d}{k}_{x}E(k_x,\omega )e^{i{k}_{x}x}\hfill \\ \hfill & =\hfill & \frac{1}{{(2\pi )}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }E(k_x,{\omega }^{\prime }+i\gamma )e^{-i{\omega }^{\prime }t+\gamma t}{e}^{i{k}_{x}x}\text{d}{k}_x\text{d}{\omega }^{\prime }.\hfill \end{array}$$

By Fubini’s theorem we may interchange the order of integration in Eq. (18), provided E(k_x,ω′ + iγ) is absolute integrable with respect to k_x and ω′. This is the case assuming that the incident field is sufficiently smooth with respect to t and x. For example, taking the incident wave to be a(x)e^iK_xxb(t)e^−iω₁t, the transformed field becomes E(k_x,ω′ + iγ) = A(k_x – K_x)B(ω′ – ω₁ + iγ), where A is the Fourier transform of a, and B is the Laplace transform of b. Here we assume that b(t) = 0 for t < 0. If a and b are continuous, A and B are absolute integrable.

We can repeat the above argument for the total field (incident + reflected, and transmitted). Assuming no superexponential instabilities, the total field is uniformly bounded:

(19)$$\left|ℰ(x,z,t)\right|\le C\text{exp}(\gamma t),$$

for positive constants C and γ. Then the transforms t → ω followed by x → k_x exist, and we can express the total field in the form (18). The total field is determined using the wave equation. In order to consider each mode k_x separately, we interchange the order of integration for each term in the wave equation. To do so, we require the second order derivatives with respect to t and x to be continuous.

It remains to prove that our solution is consistent with this requirement. From the theory in Sec. 2, we find the solution for each k_x, given a sufficiently smooth incident field. For this solution, the Fresnel equations show that the reflection and transmission coefficients tend to zero and unity, respectively, as |ω′| → ∞ or |k_x| → ∞. Therefore the reflected and transmitted field in the (ω,k_x)-domain adopt any absolute integrability property from the incident field.

In our analysis the incident field u(t)exp(ik_xx–iω₁t) is not continuous. Hence, strictly speaking, the above described method cannot be used. However, by smoothing the discontinuity around t = 0, we can make the field and its second order derivative continuous. This modification will not affect the discussion in general, since a slower transient will reduce the bandwidth. Thus the side waves are excited weaker, such that inequality (14) is satisfied with an even larger margin.

In a real experiment, not only the beam width, but also the size of the active medium itself, must be finite. Provided the fields never reach the end of the structure in the time window of interest, the fields will be identical to those in a semi-infinite active medium. Thus we can consider a setup as in Fig. 8, where the least distance from the incident beam to the boundary is d. For t < d/c the fields will be the same as if the finite-size medium were replaced by a semi-infinite medium.

 

Fig. 8 The semi-infinite gain medium can be replaced by a finite size gain medium, provided we only consider times t < d/c.

Download Full Size | PDF

B. Instabilities in infinite media

It is convenient to divide instabilities into two categories, convective and absolute instabilities (see e.g. [20,21]). Media with absolute instabilities are often regarded as impractical for small-signal, linear applications, since for an unbounded medium the fields diverge even at a fixed point in space. In contrast, media with convective instabilities are useful in the linear regime. Here the fields do not diverge at a fixed point in space; the growing wave is rather convected away.

However, even in the case with only convective instabilities, there may be fundamental problems in the case where the medium occupies an infinite region or half-space: Any small perturbation may propagate an infinite distance, thus picking up an infinite amount of gain. In our analysis we assume that the active medium is dark for t < 0. It is not clear whether this is possible, not even in principle, since perturbations in the remote past would not die out but rather increase exponentially.

The remedy is motivated by practical considerations. In an experiment, the active medium must have finite size in all directions. For a medium without absolute instabilities and with a given maximum size d, there will be no instabilities provided the gain is sufficiently weak. Examples of such configurations include optical amplifiers, and laser resonators with pumping below threshold. When there are no instabilities, we can turn on the pump in remote past such that the perturbations have died out before t = 0. For 0 < t < d/c we can still regard the medium as semi-infinite, since, as seen from Fig. 8, it makes no difference.

C. Determining the reflected time-domain field

Here we will calculate the reflected field in the time-domain, when the gain medium is described by a weak, inverted Lorentzian function:

(20)$${\varepsilon }_2(\omega )=1-\frac{F{\omega }_0^2}{{\omega }_0^2-{\omega }^2-i\omega \Gamma }.$$

In Eq. (20) F, ω₀, and Γ are positive parameters, describing the resonance strength, frequency, and bandwidth, respectively. The physical, time-domain reflected field at z = x = 0 is given by the inverse Laplace transform (12), repeated for convenience here:

(21)$${ℰ}_{\rho }(0,t)=\frac{1}{2\pi }{\int }_{-\infty +i\gamma }^{+\infty +i\gamma }\frac{k_{1z}-k_{2z}}{k_{1z}+k_{2z}}\frac{\text{exp}(-i\omega t)}{i{\omega }_1-i\omega }\text{d}\omega .$$

The field can be interpreted by evaluating integral (21) by a generalized version of the residue theorem. We here recognize that integrating along path −∞ + iγ to +∞ + iγ, is the same as integrating around all branch cuts and poles. The denominator k_1z + k_2z does not have any zeros, provided the permittivity ɛ₁ can be considered constant and larger than unity in the frequency range of interest. Thus we only need to consider the branch cuts extending from branch points of k_1z and k_2z, and the pole at ω = ω₁. Note that the branch cuts are arbitrary, as long as they extend from the branch points. We let all branch cuts lie parallel to the imaginary axis, towards Imω = −∞. See illustration in Fig. 9. The branch points of k_1z are located far away from (and below) the real frequency axis, provided the medium’s bandwidth is sufficiently large. The wavenumber k_2z has two branch points in the upper half-plane, located immediately above ω = ±k_xc. In addition there are four branch points located below the real frequency axis, with imaginary parts −Γ/2; two simple zeros and two simple poles. The integrals around the latter four branch cuts decay with time constant at most 2/Γ. Thus, for t ≳ 2/Γ, the only contributing terms are the residue of the pole at ω₁, and the contribution ℰ_bc(0,t) from the two remaining branch cuts of k_2z:

(22)$${ℰ}_{\rho }(0,t)=\frac{k_{1z}-k_{2z}}{k_{1z}+k_{2z}}\text{exp}(-i{\omega }_{1}t)+{ℰ}_{\text{bc}}(0,t).$$

Here k_1z and k_2z have been evaluated at the frequency ω₁. We write ℰ_bc(0,t) = ℰ_bc−(0,t) + ℰ_bc+(0,t), where ℰ_bc−(0,t) and ℰ_bc+(0,t) are the contributions from the branch cuts in the left and right half-planes, respectively.

 

Fig. 9 Circles and stars mark zeros and poles of k_2z respectively. The cross marks the pole at ω = ω₁. Branch cuts are arbitrarily chosen to lie parallel to the imaginary axis, extending into the lower half-plane Im(ω) < 0. The integration path C is shown with the dashed line. Contributing branch cuts and poles are enclosed by paths C₁, bc− and bc+. For t ≳ 2/Γ, we have ∮_C f(ω)dω ≈ ∮_C1 f(ω)dω + ∫_bc− f(ω)dω + ∫_bc+ f(ω)dω.

Download Full Size | PDF

Assuming F ≪ 1, Γ ≪ ω₀ and $\sqrt{2}{\omega }_0$ < k_xc, the branch cut in the right half-plane extends from approximately ω = k_xc + iδ to ω = k_xc – i∞, where δ ≤ FΓ. Then, for t ≳ 2/Γ

(23)$$\begin{array}{l}{ℰ}_{\text{bc}+}(0,t)\approx \frac{1}{2\pi }{\int }_{k_{x}c-i\frac{\Gamma }{2}}^{k_{x}c+i\delta }{\rho }_1(\omega )\frac{\text{exp}(-i\omega t)}{i{\omega }_1-i\omega }\text{d}\omega \\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}-\frac{1}{2\pi }{\int }_{k_{x}c-i\frac{\Gamma }{2}}^{k_{x}c+i\delta }{\rho }_{\text{r}}(\omega )\frac{\text{exp}(-i\omega t)}{i{\omega }_1-i\omega }\text{d}\omega .\end{array}$$

Here subscripts l and r indicate that ρ(ω) is discontinuous when crossing the branch cut, denoting the left and right side of the branch cut respectively. We further define f_l,r(ω) = k_2z/k_1z. Since k_2z is small in the vicinity of k_xc, by first order approximation ρ_l,r(ω) = 1 – 2 f_l,r(ω), where f_r(ω) = − f_l(ω). The integral (23) can now be simplified:

(24)$${ℰ}_{\text{bc}+}(0,t)\approx \frac{-2}{i\pi }{\int }_{k_{x}c-i\frac{\Gamma }{2}}^{k_{x}c+i\delta }{f}_1(\omega )\frac{\text{exp}(-i\omega t)}{{\omega }_1-\omega }\text{d}\omega ,$$

In order to obtain a manageable expression for f_l(ω), it is useful to express $k_{2z}^2$ as a function of its zeros and poles. With poles denoted by subscript p, and zeros denoted by subscripts k_xc and ω₀ (indicating the location along the real frequency axis), $k_{2z}^2$ appears as

(25)$$k_{2z}^2=\frac{(\omega -{\omega }_{{\omega }_0})(\omega +{\omega }_{{\omega }_0}^{*})(\omega -{\omega }_{k_{x}c})(\omega +{\omega }_{k_{x}c}^{*})}{c^2(\omega -{\omega }_p)(\omega +{\omega }_p^{*})}.$$

Identifying δ = Im(ω_kxc) and ω_i = Im(ω), and recognizing that (ω – ω_ω0)/(ω – ω_p) ≈ 1 at ω = ω_kxc, Eq. (25) can be simplified: $k_{2z}^2\approx -i2k_x(\delta -{\omega }_i)/c$. This gives

(26)$$\left|f_1(k_{x}c+i{\omega }_i)\right|\approx \sqrt{2(\delta -{\omega }_i)/(k_{x}c({\varepsilon }_1-1))}.$$

For k_xc – ω₁ ≫ Γ, we can now find an upper bound of integral (24) by noting that $\sqrt{\delta -{\omega }_i}\le \sqrt{F\Gamma +\Gamma /2}$ for all ω_i considered. We can estimate ℰ_bc−(0,t) similarly, yielding the bound

(27)$$\left|{ℰ}_{\text{bc}}(0,t)\right|\le \frac{{\Gamma }^{3/2}}{(k_{x}c-{\omega }_1)\sqrt{k_{x}c({\varepsilon }_1-1)}}\left(e^{F\Gamma t}-e^{-\Gamma t/2}\right)$$

Consequently for 2/Γ ≲ t ≲ 1/FΓ, the field is well described by the first term in Eq. (22).

References and links

1. G. N. Romanov and S. S. Shakhidzhanov, “Amplification of electromagnetic field in total internal reflection from a region of inverted population,” JETP Lett. 16, 209–211 (1972).

2. A. A. Kolokolov, “Reflection of plane-waves from an amplifying medium,” JETP Lett. 21, 312–313 (1975).

3. P. R. Callary and C. K. Carniglia, “Internal-reflection from an amplifying layer,” J. Opt. Soc. Am. 66, 775–779 (1976). [CrossRef]  

4. W. Lukosz and P. P. Herrmann, “Amplification by reflection from an active medium,” Opt. Commun. 17, 192–195 (1976). [CrossRef]  

5. R. F. Cybulski and C. K. Carniglia, “Internal-reflection from an exponential amplifying region,” J. Opt. Soc. Am. 67, 1620–1627 (1977). [CrossRef]  

6. R. F. Cybulski and M. P. Silverman, “Investigation of light amplification by enhanced internal-reflection .1. theoretical reflectance and transmittance of an exponentially nonuniform gain region,” J. Opt. Soc. Am. 73, 1732–1738 (1983). [CrossRef]  

7. A. A. Kolokolov, “Determination of the reflection coefficient of a plane monochromatic wave,” J. Commun. Technol. Electron. 43, 837–845 (1998).

8. A. A. Kolokolov, “Fresnel formulas and the principle of causality,” Phys. Usp. 42, 931–940 (1999). [CrossRef]  

9. J. Fan, A. Dogariu, and L. J. Wang, “Amplified total internal reflection,” Opt. Express 11, 299–308 (2003). [CrossRef]   [PubMed]  

10. K. J. Willis, J. B. Schneider, and S. C. Hagness, “Amplified total internal reflection: theory, analysis, and demonstration of existence via FDTD,” Opt. Express 16, 1903–1913 (2008). [CrossRef]   [PubMed]  

11. A. Siegman, “Fresnel reflection, Lenserf reflection, and evanescent gain,” Opt. Photon. News 21, 38–45 (2010). [CrossRef]  

12. C. J. Koester, “Laser action by enhanced total internal reflection,” IEEE J. Quantum Electron. 2, 580–584 (1966). [CrossRef]  

13. B. Y. Kogan, V. M. Volkov, and S. A. Lebedev, “Superluminescence and generation of stimulated radiation under internal-reflection conditions,” JETP Lett. 16, 100–101 (1972).

14. S. A. Lebedev, V. M. Volkov, and B. Y. Kogan, “Value of the gain for light internally reflected from a medium with inverted population,” Opt. Spectrosc. 35, 565–566 (1973).

15. M. P. Silverman and J. R. F. Cybulski, “Investigation of light amplification by enhanced internal reflection. Part II. Experimental determination of the single-pass reflectance of an optically pumped gain region,” J. Opt. Soc. Am. 73, 1739–1743 (1983). [CrossRef]  

16. B. Nistad and J. Skaar, “Causality and electromagnetic properties of active media,” Phys. Rev. E 78, 036603 (2008). [CrossRef]  

17. J. Skaar, “Fresnel equations and the refractive index of active media,” Phys. Rev. E 73, 026605 (2006). [CrossRef]  

18. A. Lakhtakia, J. B. Geddes III, and T. G. Mackay, “When does the choice of the refractive index of a linear, homogeneous, isotropic, active, dielectric medium matter?” Opt. Express 15, 17709–17714 (2007). [CrossRef]   [PubMed]  

19. L. V. Ahlfors, Complex Analysis (McGraw-Hill International Editions, 1979).

20. P. A. Sturrock, “Kinematics of growing waves,” Phys. Rev. 112, 1488–1503 (1958). [CrossRef]  

21. R. J. Briggs, Electron-Stream Interactions with Plasmas (MIT Press, 1964).