Husimi functions at gradient index cavities designed by conformal transformation optics

Inbo Kim; Jinhang Cho; Yushin Kim; Bumki Min; Jung-Wan Ryu; Sunghwan Rim; Muhan Choi

doi:10.1364/OE.26.006851

1. Introduction

Maxwell's equations, the governing equations for electromagnetic fields in space-time including media, have fundamental symmetries under various transformations such as well-known Lorentz transformation and U(1) gauge transformation. Maxwell's equations are also form-invariant under general coordinate transformations [1–3]. Exploiting the form invariance of Maxwell's equations for general coordinate transformations, Pendry et al. theoretically presented transformation optics (TO) [4] which is a general methodology for designing electromagnetic materials including optical invisibility cloaks, and various photonic devices which manipulate the path of light waves [5, 6]. Particularly, very recently, applying the TO to the design of optical dielectric resonators, Kim et al. reported 2-dimensional (2D) gradient index dielectric cavities called transformation cavities by using conformal TO [7–9]. Surprisingly, in their transformation cavities, directional whispering gallery modes (WGMs) with very high quality factors (Q-factor) are possible, providing a great potential for cutting-edge optical device design.

Conventional 2D optical dielectric cavities with uniform refractive indices have been regarded as an open system version of closed quantum billiard systems and also as a potential candidate of low-threshold micro-lasers, so there have been vast and intensive works on these platform systems over the last twenty years [10, 11]. In order to analyze the characteristics of optical modes supported in dielectric cavities, it is often more useful to represent the optical modes in phase space than merely depict the mode intensity patterns in real physical space. Husimi function, one of the popular phase space representations of a wave function can be regarded as a quasi-probability distribution in phase space. It is defined by the overlap of the wave function with a coherent state of a minimum uncertainty wave packet [12]. The so-called Poincaré Husimi functions on the billiards boundary have been widely used when studying quantum-classical correspondence in quantum billiards [13,14]. For the case of optical dielectric cavities, Hentschel et al. have derived four different Husimi functions in a reduced phase space, i.e., on the Poincaré surface of section (PSOS) at the dielectric interfaces by using a stationary phase approximation [15].

However, since the Husimi functions of Hentschel et al. are valid exclusively for homogeneous dielectric interfaces, it is inappropriate to apply their Husimi functions directly to the transformation cavities to see the phase space structure of an optical mode formed there. In this paper, we show that the Husimi functions can still be useful for optical modes supported in transformation cavities with inhomogeneous refractive indices. To this end, we construct a virtual space which we call a reciprocal virtual (RV) space, through the inverse of the conformal mapping, then both polarized resonant modes supported in the transformation cavity of the physical space and the corresponding virtual modes in the unit disk cavity of the RV space turn out to be identical. Besides, the refractive index of the unit disk cavity in the RV space is uniform. Thus, based on these facts, one can obtain a phase space description of the internal waves in a transformation cavity through the conventional Poincaré Husimi functions calculated in the RV space.

This paper is organized as follows. In section 2, starting from a uniform index unit disk cavity in a fictitious space (called as an original virtual space) we describe the gradient index transformation cavity model. In section 3, RV space is introduced to apply conventional Husimi functions for the internal waves in a uniform index cavity. As an illustration of our method, Husimi functions for a directional high-Q WGM and a directional low-Q mode in a limaçon-shaped transformation cavity are calculated and analyzed in section 4. In the last section comes our conclusion.

2. Gradient index cavities designed by optical conformal mapping

We start with a homogeneous dielectric cavity of an infinite cylinder with a unit disk cross section in an original virtual (OV) space with Cartesian coordinates $(u, v, w)$ , taking $w$ -axis along the axis of the cylinder as shown in Fig. 1(a). In the OV space, the uniform refractive index of the cavity is $n$ and the refractive index of the exterior region is 1 as the refractive index of vacuum or practically that of air. In the physical space with Cartesian coordinates $(x, y, z)$ , a realizable 2D gradient index cavity model (hereafter called ‘transformation cavities’) can be built by applying a conformal coordinate transformation only to the inside region ( $Ω_{1}^{O V}$ ) of the unit disk cavity in the OV space and assigning the uniform refractive index of $n_{0} = 1$ for the outside region $(Ω_{0}$ ) of the transformed cylindrical cavity ( $z = w$ ) with a deformed cross section as shown schematically in Fig. 1(b). For these transformation cavities with deformed shapes, it was demonstrated that directional WGMs called ‘conformal WGMs’ (cWGMs) can be formed [7]. The conformal coordinate transformation from the OV space to the physical space is given by an analytic function $ζ = f (η) = x (u, v) + i y (u, v)$ of a complex variable $η = u + i v$ , that satisfies the Cauchy–Riemann equations,

\frac{\partial x}{\partial u} = \frac{\partial y}{\partial v}, \frac{\partial y}{\partial u} = - \frac{\partial x}{\partial v} .

Conformal mappings preserve angles locally; both angles and the shapes of infinitesimally small figures are preserved, but not necessarily their size or curvature. Also the handedness of coordinate system is unchanged.

Fig. 1 Schematic illustration of a transformation cavity and its pertaining two virtual spaces: (a) the unit disk cavity in an original virtual (OV) space, (b) a limaçon-shaped transformation cavity in the physical space, (c) the corresponding unit disk cavity in a reciprocal virtual (RV) space. (refractive index profiles are expressed by common color scale.)

Download Full Size | PDF

Due to the translation symmetry along the $z$ -axis of the aforementioned cylindrical geometry, Maxwell’s equations simplify to a 2D scalar wave equation for the two polarization of light with respect to the $z$ -axis [10]. Optical modes in a transformation cavity with a deformed boundary are found by solving the following 2D scalar wave equation,

[Δ + n^{2} (r) k^{2}] ψ (r) = 0,

where

Δ \equiv \partial_{x}^{2} + \partial_{y}^{2}

, position vector,

r = (x, y) = (r \cos θ, r sin θ)

and the refractive index

n (r)

is given by

n (r) = {\begin{matrix} n {| \frac{d ζ}{d η} |}^{- 1}, r Ω_{1} (i n t e r i o r) \\ 1, r Ω_{0} (e x t e r i o r) \end{matrix},

with the radiation boundary condition,

ψ (r) ~ h (θ, k) \frac{e^{i k r}}{\sqrt{r}} for r \to \infty,

where

h (θ, k)

is the far-field angular distribution of the radiation emission. The radiation boundary condition leads to solutions exponentially decaying in time with discrete complex eigenvalues

k

with

Im (k) < 0

. Complex angular frequency is

ω = c k

, where

c

is the speed of light in a vacuum given by

1 / \sqrt{ε_{0} μ_{0}}

. The lifetime

τ

of these so-called ‘resonant modes’ or ‘resonances’ is given by the imaginary part of the wave number as

τ = - 1 / 2 Im (ω)

; the quality factor of each resonance is defined by

Q = 2 π τ / T = - Re (ω) / 2 Im (ω),

where the oscillation period of light wave

T = 2 π / Re (ω)

. The complex wave function

ψ (r)

represents

E_{z}

and

H_{z}

with regard to the transverse magnetic (TM) and the transverse electric (TE) polarizations, respectively, i.e., the TM polarized time-harmonic electric field is given by

E (r, t) = E_{z} (r, t) \hat{z} = Re [ψ (r) e^{- i ω t}] \hat{z}

and the TE polarized

H

-field is given by

H (r, t) = H_{z} (r, t) \hat{z} = Re [ψ (r) e^{- i ω t}] \hat{z}

. The boundary conditions at the cavity–air interface are neither of Dirichlet nor of Neumann type but a mixed type given by [10, 16]

ψ_{1} = ψ_{0}, \partial_{ν} ψ_{1} = \partial_{ν} ψ_{0} for TM modes,

ψ_{1} = ψ_{0}, \partial_{ν} ψ_{1} = \partial_{ν} ψ_{0} for TM modes,

where

ψ_{j}

are wave functions on the cavity boundary from the interior region (

j = 1

) and the exterior region (

j = 0

) and

\partial_{ν} ψ_{j}

are their normal derivatives, respectively;

\partial_{ν} \equiv v (r) \cdot \nabla |_{r}

and

v (r)

is the outward normal unit vector on the cavity boundary curve

Γ

at point

r

and

n_{1}

is the gradually-varying refractive index along the boundary, i.e.,

n_{1} = n {| d ζ / d η |}^{- 1}

evaluated at the boundary.

3. Poincaré Husimi function at transformation cavities

The Husimi function, a quasiprobability distribution in phase space is originally defined as the overlap of the wave function with a coherent state that represents a minimal-uncertainty wave packet. Crespi et al. derived so-called a Poincaré Husimi function at the boundary of 2D closed billiard systems by projecting the conventional Husimi function from full 4-dimensional phase space onto a 2D reduced phase space, namely on the PSOS at the system boundary [13]. It has been widely used in 2D closed billiard systems in the context of quantum chaos, especially in order to study the semiclassical regime [14, 17]. For the case of dielectric systems with piecewise constant refractive indices, Hentschel, et al. have derived four different Poincaré Husimi functions, corresponding to incident and emerging waves on inside and outside of the interface by using a semiclassical (saddle point) approximation with a Green’s function technique [15]. Among them, the Husimi functions for the internal waves have been useful when exploring ray-wave correspondence in optical dielectric cavities.

Their expression for circular cavity can be readily extended to cavities with deformed boundaries as below. Four different Husimi functions which represent the intensities of the incident and emerging waves from arc length $s$ along the cavity boundary $Γ$ in the directions $χ_{j}$ ( $j = 1$ ; interior region, $j = 0$ ; exterior region), respectively are (see Fig. 2)

H_{j}^{i n c (e m)} (s, \sin χ_{j}) = \frac{k_{j}}{2 π} {| {(- 1)}^{j} F_{j} h_{j} (s, \sin χ_{j}) + (-) \frac{i}{k_{0} F_{j}} h_{j}^{'} (s, \sin χ_{j}) |}^{2}

with weighting factors

F_{j} = \sqrt{n_{j} \cos χ_{j}}

and the functions

h_{j}

,

h_{j}^{'}

given by

h_{j} (s, \sin χ_{j}) = \oint_{Γ} d s^{'} ψ_{j} (s^{'}) ξ (s^{'}; s, \sin χ_{j}),

{h^{'}}_{j} (s, s i n χ_{j}) = \oint_{Γ} d s^{'} \partial_{ν} ψ_{j} (s^{'}) ξ (s^{'}; s, \sin χ_{j}),

which are the overlaps of the wave function

ψ

and its normal derivative

\partial_{ν} ψ

, taken on the respective side

j

of the cavity boundary, with the minimum-uncertainty wave packet

ξ (s^{'}; s, \sin χ_{j}) = \frac{1}{\sqrt[4]{σ π}} \sum_{l \in ℤ} exp [\frac{- {(s^{'} - s + 2 π l)}^{2}}{2 σ} - i k_{j} \sin χ_{j} (s^{'} + 2 π l)],

periodic in

s^{'}

and centered around

(s, \sin χ_{j})

. The uniform refractive index of interior region is

n_{1} = n

and the refractive index of exterior region is

n_{0} = 1

; the

k_{j}

are the wave number in each region. The width of the wave packet in the

s

-direction along with the uncertainty in

\sin χ_{j}

can be controlled with the parameter

σ = \sqrt{2} / k_{1}

.

Fig. 2 Incident and emerging rays at a dielectric interface with uniform refractive indices.

Download Full Size | PDF

The boundary Husimi function of each region is valid when the refractive index is constant in each region. Therefore, in the case of the transformation cavities the Eq. (6) can be used for the exterior region where the refractive index is 1, but the formula cannot be directly applied to the interior region. However, fortunately, if we construct a virtual space through inverse conformal transformation, a transformation cavity is mapped back to the original unit disk cavity with a constant refractive index, so the boundary Husimi function for the resonant modes in the interior region can be calculated there by the above formula. In other words, it becomes possible to look at the phase space structure that shows the emission direction of the resonant modes and the location on the boundary where the major emission leaks out.

So, we introduce the reciprocal virtual (RV) space with Cartesian coordinates $(u, v, w)$ , which is obtained from the physical space by the inverse conformal mapping, $η = f^{- 1} (ζ)$ , as shown in Fig. 1(c). It is noted that the inverse conformal mapping is not a one-to-one mapping and therefore the RV space corresponding to the OV space should be selected. Functions, differential operators, and relevant symbols in the RV space will be expressed with tildes. Under the inverse conformal mapping, the scalar wave Eq. (2) is transformed to

[\tilde{Δ} + {\tilde{n}}^{2} (\tilde{r}) k^{2}] \tilde{ψ} (\tilde{r}) = 0,

where

\tilde{Δ} \equiv \partial_{u}^{2} + \partial_{v}^{2}

, position vector

\tilde{r} = (u, v)

, and refractive index

\tilde{n} (\tilde{r})

is given by

\tilde{n} (\tilde{r}) = {\begin{matrix} n, \tilde{r} \in {\tilde{Ω}}_{1} (i n t e r i o r) \\ {| \frac{d η}{d ζ} |}^{- 1} \tilde{r} \in {\tilde{Ω}}_{0} (e x t e r i o r) \end{matrix} .

Notice that the intracavity refractive index in the RV space is same as that of the OV space but the exterior index varies spatially as shown in Fig. 1(c). Incidentally, if we take the RV space as an OV space, the transformation cavity in the physical space can be obtained from the RV space without artificial setting of

n_{0} = 1

for the outside region.

In the RV space, the Husimi functions for the incident and the emerging waves in the unit disk cavity with uniform index $n$ can be written as below.

H_{1}^{i n c (e m)} (\tilde{s}, \sin {\tilde{χ}}_{1}) = \frac{k_{1}}{2 π} {| F_{1} h_{1} (\tilde{s}, \sin {\tilde{χ}}_{1}) - (+) \frac{i n}{k_{1} F_{1}} {h^{'}}_{1} (\tilde{s}, \sin {\tilde{χ}}_{1}) |}^{2}

with weighting factor

F_{1} = \sqrt{n \cos {\tilde{χ}}_{1}}

and the functions

h_{1}

,

{h^{'}}_{1}

given by

h_{1} (\tilde{s}, \sin {\tilde{χ}}_{1}) = \oint_{\tilde{Γ}} d {\tilde{s}}^{'} {\tilde{ψ}}_{1} ({\tilde{s}}^{'}) ξ ({\tilde{s}}^{'}; \tilde{s}, \sin {\tilde{χ}}_{1}),

{h^{'}}_{1} (\tilde{s}, \sin {\tilde{χ}}_{1}) = \oint_{\tilde{Γ}} d {\tilde{s}}^{'} {\tilde{\partial}}_{ν} {\tilde{ψ}}_{1} ({\tilde{s}}^{'}) ξ ({\tilde{s}}^{'}; \tilde{s}, \sin {\tilde{χ}}_{1}),

which are the overlaps of the wave function

\tilde{ψ}

and its normal derivative

{\tilde{\partial}}_{ν} \tilde{ψ}

, taken on the inside (

j = 1

) of the unit disk cavity at arc length

{\tilde{s}}^{'}

along the boundary

\tilde{Γ}

in the RV space, with the minimum-uncertainty wave packet

ξ ({\tilde{s}}^{'}; \tilde{s}, \sin {\tilde{χ}}_{1}) = \frac{1}{\sqrt[4]{σ π}} \sum_{l \in ℤ} exp [\frac{- {({\tilde{s}}^{'} - \tilde{s} + 2 π l)}^{2}}{2 σ} - i k_{1} \sin {\tilde{χ}}_{1} ({\tilde{s}}^{'} + 2 π l)],

which is periodic in

{\tilde{s}}^{'}

and centered around

(\tilde{s}, \sin {\tilde{χ}}_{1})

. The parameter

σ = \sqrt{2 /} k_{1}

controls the uncertainty in

\tilde{s}

and

\sin {\tilde{χ}}_{1}

. For a resonant mode,

k_{1} = n Re [k]

where

k

is a complex wave number of the resonant mode.

In order to calculate the Husimi function in the RV space, it is necessary to obtain $\tilde{ψ}$ , ${\tilde{\partial}}_{ν} \tilde{ψ}$ on the boundary of the unit disk cavity. To obtain these from the $ψ$ , ${\tilde{\partial}}_{ν} ψ$ in physical space, we need to know the transformation rules between them. According to TO [4,18], 3-dimensional transformation rules for fields and gradient operators between the physical space and the RV space can be expressed as follows,

E = {(Λ^{- 1})}^{T} \tilde{E}, H = {(Λ^{- 1})}^{T} \tilde{H}, \tilde{\nabla} = {(Λ^{- 1})}^{T} \tilde{\nabla} .

where the Jacobian matrix

Λ

and its inverse

Λ^{- 1}

are given by

Λ = [\begin{matrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} & 0 \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} & 0 \\ 0 & 0 & 1 \end{matrix}], Λ^{- 1} = [\begin{matrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} & 0 \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} & 0 \\ 0 & 0 & 1 \end{matrix}]

because

z = w

,

\partial x / \partial w = \partial y / \partial w = 0

, and

\partial z / \partial u = \partial z / \partial v = 0

. Due to the Cauchy–Riemann equations, the scale factor of conformal mapping

ζ = f (η)

is given by

\det (Λ^{- 1})

which equals

{| d ζ / d η |}^{- 2}

.

First of all, for TM/TE polarizations, wave functions in the physical space and their counterpart functions in the RV space are identical as can be shown by using the above field transformation rules;

(\begin{matrix} 0 \\ 0 \\ ψ \end{matrix}) = [\begin{matrix} \frac{\partial u}{\partial x} & \frac{\partial v}{\partial x} & 0 \\ \frac{\partial u}{\partial y} & \frac{\partial v}{\partial y} & 0 \\ 0 & 0 & 1 \end{matrix}] (\begin{matrix} 0 \\ 0 \\ \tilde{ψ} \end{matrix}) .

So,

ψ (x, y) = \tilde{ψ} (u, v)

with

x + i y = f (u + i v)

between the physical space and the RV space. Secondly, from the complex analysis [19], the transformation rule between unit normal vectors

ν_{x} + i ν_{y}

and

ν_{u} + i ν_{v}

at the cavity boundaries in the respective spaces can be written as

ν_{x} + i ν_{y} = (ν_{u} + i ν_{v}) \frac{d ζ}{d η} {| \frac{d ζ}{d η} |}^{- 1} .

Using the Cauchy–Riemann equations, Eq. (17) can be rewritten in 3-D Cartesian vector notation with

ν = {(ν_{x} ν_{y} 0)}^{T}

and

\tilde{v} = {(ν_{u} ν_{v} 0)}^{T}

. as

ν = Λ \tilde{ν} / \sqrt{\det Λ}

.

Thus, normal derivative transforms as

\partial_{ν} \equiv ν \cdot \nabla = \tilde{ν} \cdot \frac{Λ^{T} {(Λ^{T})}^{- 1}}{\sqrt{\det Λ}} \tilde{\nabla} = {| \frac{d ζ}{d η} |}^{- 1} {\tilde{\partial}}_{ν} .

Therefore, through the relations

ψ (s) = \tilde{ψ} (\tilde{s})

and

\partial_{ν} ψ (s) = {| d ζ / d η |}^{- 1} {\tilde{\partial}}_{ν} \tilde{ψ} (\tilde{s})

, one can obtain the

\tilde{ψ}

,

{\tilde{\partial}}_{ν} \tilde{ψ}

on the boundary of the unit disk cavity in the RV space from the

ψ

,

{\tilde{\partial}}_{ν} \tilde{ψ}

on the boundary of a transformation cavity and calculate the Husimi functions (Eq. (11)) for the internal waves in the RV space.

4. Examples: a limaçon-shaped transformation cavity

As examples, we calculate the Husimi function for incident waves in the unit disk cavity in the RV space for a high-Q resonance (cWGM) and a low-Q resonance in a limaçon-shaped transformation cavity. The conformal mapping which transforms the unit circle to the limaçon [7,20], is given by

ζ = β (η + ϵ η^{2}),

where

η = u + i v

and

ζ = x + i y

are complex variables that denote positions in the respective complex planes,

ϵ

is a deformation parameter, and

β

is a positive scaling factor. We use

ϵ = 0.23

,

β = 0.685

and

n = 1.8

, the refractive index of unit disk cavity in the OV space. Wave functions of resonant modes formed in the transformation cavity can be calculated by a boundary element method (BEM) exploiting the RV space [7,21] or finite element method (FEM) based electromagnetic field solvers e.g., COMSOL Multiphysics.

Using the BEM, we obtained a high-Q TM mode which is identified as a cWGM with a complex wave number $k = 12.2333 - i 0.0022$ . The Q-factor of the resonance is 2776.23. The mode intensity distribution, the far field pattern and the intracavity Husimi function for incident waves are depicted in Fig. 3 ( $\tilde{s}$ is normalized with $2 π$ , the total arc length of the unit disk cavity). As can be seen in Fig. 3(a), the cWGM is well confined by total internal reflection of light around the rim of the limaçon-shaped transformation cavity and this feature is reflected in the Husimi plot as the two band-type intensities located in the regions $| \sin {\tilde{χ}}_{1} | > {| d η / d ζ |}_{\tilde{s} (\tilde{r})}^{- 1} / n$ where the r.h.s. of the inequality represents the critical lines of total internal reflection depicted by two yellow solid curves in Fig. 3(c) (the $\pm$ signs in $\sin {\tilde{χ}}_{1}$ denote counterclockwise (CCW) and the clockwise (CW) circulations of light waves, respectively). The bi-directional far field emission shown in Fig. 3(b) can be explained by the Husimi plot in the RV space in Fig. 3(c); the closest position between the intensity bands and critical lines in $\tilde{s}$ is $\tilde{s} = 0$ or 1 which corresponds to the position $ζ = ($ 0.84255, 0) on the boundary of the limaçon-shaped transformation cavity where maximum evanescent leakage or tunneling emissions come out in tangential directions i.e., parallel to the $y$ -axis [7].

Fig. 3 A high-Q TM resonance (cWGM) in a limaçon-shaped transformation cavity ( $k = 12.2333 - i 0.0022$ ) (a) the mode intensity distribution; two yellow arrows denote the directions of tunneling emissions at $\tilde{s} = 0$ in Fig. 3(c), (b) the far field pattern, (c) the intracavity Husimi plot for incident waves in the RV space; two yellow solid curves are critical lines of total internal reflection and the $\pm$ signs in $\sin {\tilde{χ}}_{1}$ denote the counterclockwise (CCW) and clockwise (CW) circulations of light waves, respectively. ( $\tilde{s}$ is normalized with $2 π$ , the total arc length of the unit disk cavity.)

Download Full Size | PDF

Also, we obtained a low-Q TM mode with $k = 12.5960 - i 0.1302$ whose intensity pattern is depicted in Fig. 4(a). The Q-factor of the resonance is 48.37. The far field emission shown in Fig. 4(b) is refractive emission as can be seen clearly from the Husimi plot in Fig. 4(c) where the intensity peaks of CW and CCW waves are located below the critical lines, i.e., in the regions $| \sin {\tilde{χ}}_{1} | < {| d η / d ζ |}_{\tilde{s} (\tilde{r})}^{- 1} / n$ . To analyze this refractive emission approximately, we resort to Snell’s law, $n_{1} \sin χ_{1} = \sin χ_{0}$ in the physical space. Because of conformality of the transformation, one can take $\sin χ_{1} = \sin {\tilde{χ}}_{1} ~ - 0.385$ , the peak position of lower intensity in the Husimi plot. The peak position in $\tilde{s}$ coordinate is $\tilde{s} = 0.1256$ which can be converted into $s = 0.1732$ , normalized value with the whole boundary length of the limaçon. The intracavity refractive index at this position is $n_{1}$ = 1.9314. So, the refraction angle, $χ_{0} = - {48.04}^{°}$ is obtained from $\sin χ_{0} ~ - 0.7436$ . In polar coordinate, angular coordinate of the main emission point on the boundary of the limaçon-shaped transformation cavity is 0.929 (rad) or ${53.226}^{°}$ and the outward normal direction at the point is ${58.93}^{°}$ . Thus, we can conclude that the main refractive emission comes out at the boundary point of angular coordinate ${53.226}^{°}$ in the direction ${10.89}^{°}$ . By the mirror symmetry with respect to the $x$ -axis, the peak of upper intensity in the Husimi plot corresponds to the refractive emission at the boundary point of angular coordinate $- {53.226}^{°}$ in the direction $- {10.89}^{°}$ . Thus the far field pattern in Fig. 4(b) can be understood as the interference of the two dominant (out of phase) refractive emission.

Fig. 4 A low-Q TM resonance in a limaçon-shaped transformation cavity ( $k = 12.5960 - i 0.1302$ ) (a) the mode intensity distribution; two yellow arrows designate dominant refractive emissions at the corresponding peak positions $\tilde{s} = 0.1256$ and $\tilde{s} = 0.8744$ in the Husimi plot of Fig. 4(c), (b) the far field pattern, (c) the intracavity Husimi plot for incident waves in the RV space; two yellow solid curves are critical lines of total internal reflection same as in Fig. 3(c) and two yellow arrows point the intensity peaks of CW and CCW waves ( $\tilde{s}$ is normalized with $2 π$ , the total arc length of the unit disk cavity.)

Download Full Size | PDF

5. Conclusion

A method to find a widely used phase space representation i.e., the Husimi representation for the internal waves on the boundary of a transformation cavity was presented by constructing a virtual space with a uniform dielectric cavity via inverse conformal transformation. As a verification of the method, Husimi plots for a high-Q and a low-Q resonance were obtained in the virtual space and their emission characteristics could be revealed from the plots. The phase space description of the corresponding modes in the virtual space agrees well with the far-field directionality of the resonant modes of a transformation cavity in the physical space. We expect that our method will be useful in studying waves in transformation cavities.

Funding

National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2017R1A2B4012045 and No. 2017R1A4A1015565); Center for Advanced Meta-Materials (CAMM) funded by Korea Government (MSIP) as Global Frontier Project (CAMM 2014M3A6B3063709); Institute for Basic Science in Korea (IBS-R024-D1).

Acknowledgments

We appreciate our previous collaborator, late Professor Dr. S.-Y. Lee.

References and links

1. E. J. Post, Formal Structure of Electromagnetics (Wiley, 1962).

2. J. Plebanski, “Electromagnetic waves in gravitational fields,” Phys. Rev. 118(5), 1396–1408 (1960). [CrossRef]

3. I. V. Lindell, Methods for Electromagnetic Field Analysis (Oxford University, 1992).

4. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]

5. O. Ozgun and M. Kuzuoglu, “Form Invariance of Maxwell’s Equations: The Pathway to Novel Metamaterial Specifications for Electromagnetic Reshaping,” IEEE Antennas Propag. Mag. 52(3), 51–65 (2010). [CrossRef]

6. Y. Liu and X. Zhang, “Recent advances in transformation optics,” Nanoscale 4(17), 5277–5292 (2012). [CrossRef] [PubMed]

7. Y. Kim, S. H. Lee, J. W. Ryu, I. Kim, J. H. Han, H. S. Tae, M. Choi, and B. Min, “Designing whispering gallery modes via transformation optics,” Nat. Photonics 10(10), 647–652 (2016). [CrossRef]

8. U. Leonhardt, “Optical Conformal Mapping,” Science 312(5781), 1777–1780 (2006). [CrossRef] [PubMed]

9. L. Xu and H. Chen, “Conformal transformation optics,” Nat. Photonics 9(1), 15–23 (2015). [CrossRef]

10. H. Cao and J. Wiersig, “Dielectric microcavities: Model systems for wave chaos and non-Hermitian physics,” Rev. Mod. Phys. 87(1), 61–111 (2015). [CrossRef]

11. T. Harayama and S. Shinohara, “Two-dimensional microcavity lasers,” Laser Photonics Rev. 5(2), 247–271 (2011). [CrossRef]

12. K. Husimi, “Some formal properties of the density matrix,” Proc. Phys.-Math. Soc. Jpn. 22, 264–314 (1940).

13. B. Crespi, G. Perez, and S.-J. Chang, “Quantum Poincaré sections for two-dimensional billiards,” Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics 47(2), 986–991 (1993). [CrossRef] [PubMed]

14. H.-J. Stöckmann, Quantum Chaos (Cambridge University, 1999).

15. M. Hentschel, H. Schomerus, and R. Schubert, “Husimi functions at dielectric interfaces: Inside-outside duality for optical systems and beyond,” Europhys. Lett. 62(5), 636–642 (2003). [CrossRef]

16. W. C. Chew, Waves and Fields in Inhomogeneous Media (Wiley-IEEE, 1999).

17. A. Bäcker, S. Fürstberger, and R. Schubert, “Poincare Husimi representation of eigenstates in quantum billiards,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(3), 036204 (2004). [CrossRef] [PubMed]

18. D. Liu, L. H. Gabrielli, M. Lipson, and S. G. Johnson, “Transformation inverse design,” Opt. Express 21(12), 14223–14243 (2013). [CrossRef] [PubMed]

19. T. Needham, Visual Complex Analysis (Clarendon, 2002).

20. M. Robnik, “Classical dynamics of a family of billiards with analytic boundaries,” J. Phys. A 16(17), 3971–3986 (1983). [CrossRef]

21. In preparation (to be published elsewhere).

Husimi functions at gradient index cavities designed by conformal transformation optics

Abstract

1. Introduction

2. Gradient index cavities designed by optical conformal mapping

3. Poincaré Husimi function at transformation cavities

4. Examples: a limaçon-shaped transformation cavity

5. Conclusion

Funding

Acknowledgments

References and links

Cited By

Figures (4)

Equations (22)

Optics Express