Next Article in Journal
Quantum Algorithm for Variant Maximum Satisfiability
Next Article in Special Issue
Universal Single-Mode Lasing in Fully Chaotic Billiard Lasers
Previous Article in Journal
DA-LSTM-VAE: Dual-Stage Attention-Based LSTM-VAE for KPI Anomaly Detection
Previous Article in Special Issue
Mather β-Function for Ellipses and Rigidity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ray–Wave Correspondence in Microstar Cavities

Institut für Physik, Otto-von-Guericke-Universität Magdeburg, Postfach 4120, D-39016 Magdeburg, Germany
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(11), 1614; https://doi.org/10.3390/e24111614
Submission received: 14 October 2022 / Revised: 1 November 2022 / Accepted: 3 November 2022 / Published: 5 November 2022

Abstract

:
In a previous work published by the authors in 2020, a novel concept of light confinement in a microcavity was introduced which is based on successive perfect transmissions at Brewster’s angle. Hence, a new class of open billiards was designed with star-shaped microcavities where rays propagate on orbits that leave and re-enter the cavity. In this article, we investigate the ray–wave correspondence in microstar cavities. An unintuitive difference between clockwise and counterclockwise propagation is revealed which is traced back to nonlinear resonance chains in phase space.

1. Introduction

Optical microcavities are fascinating systems to confine light on a very small scale [1,2]. Traditional whispering-gallery cavities rely on successive total internal reflections at the cavity’s dielectric interface which prevent a ray from leaving the cavity. Correspondingly, long-lived whispering-gallery modes exist as solutions of Maxwell’s equations. Such whispering-gallery cavities appear, e.g., as microdisks [3,4,5,6], microtoroid [7,8,9], microspheres [10], or microjet cavities [11]. They have attracted an immense amount of research interest, e.g., for quantum chaos [12,13,14,15] or exceptional points in non-Hermitian systems [16,17,18,19,20] and, additionally, offer a lot of applications, e.g., the generation of optical frequency combs [21], as microlasers with directional emission [22,23,24,25], optical sensors for nanoparticles [26,27,28] or rotating motion [29], or orbital angular momentum lasers [30].
In a former article [31], an alternative concept for light confinement was introduced which is completely different from the traditional whispering-gallery design. The perfect transmissions at Brewster’s angle are utilized such that a ray can leave (and reenter) the cavity without partial back reflection. Thus, via successive transmissions through the dielectric interface of the cavity, rays are guided along a periodic orbit without loss of intensity. This idea can be implemented in a star-shaped cavity as shown in Figure 1a or in Brewster-notched cavities [32].
In this article, we further investigate the ray–wave correspondence in the microstar cavity. In particular we look at (asymmetric) deformed microstar cavities where the formation of the modes is governed by the regular dynamics around an elliptical fixed point in phase space. In the semiclassical regime, we observe modes localizing along the regular structures in phase space, while this localization might be expected from quantum-chaos theory, we can also identify an unexpected difference between clockwise (CW) and counterclockwise (CCW) propagation in the star cavity. In the previous literature, a general imbalance of CW and CCW propagation due to asymmetric backscattering was discussed for microcavities [19,33,34,35]. However, in this article, we study the difference pattern between CW and CCW propagation that is due to a nonuniform loss across a resonance chain. Therefore, the revealed mechanism is a general feature of open systems and not restricted to asymmetric boundary deformations. In addition, a Frobenius–Perron operator (FPO) formalism is used as a connection between ray dynamics and optical modes.
The paper is organized as follows. In Section 2, the ray dynamics of the microstar cavity are issued with a special focus on the dynamics of CW and CCW propagation in phase space. The corresponding impacts on the wave dynamics are discussed in Section 3. A conclusion is given in Section 4.

2. Ray Dynamics and Phase-Space Description for Microstar Cavities

A microstar is a quasi-two-dimensional cavity with homogeneous refractive index n [31]. It is characterized by the number of spikes ν . In order to guide a ray via successive transmissions at Brewster’s angle χ B = arctan ( 1 / n ) along a closed orbit, both the refractive index and the opening angle α of the spikes need to be adjusted as
n = 1 + sin ( π / ν ) cos ( π / ν )
α = ν 2 2 ν π .
In the following, a microstar with ν = 9 spikes is considered with a refractive index of n 1.43 and spike opening angle α = 70 . The microstar with such a polygonal shape supports a family of periodic orbits propagating without loss of intensity. However, these periodic orbits are only marginally stable. A slight boundary deformation of the microstar can be used to stabilize the dynamics around one central orbit as shown in Figure 1a. In polar coordinates ( r , ϕ ) , the boundary of the most right deformed spike can be expressed as
r ( ϕ ) = 1 n sin ( ϕ ) + cos ( ϕ ) ε 1 ϕ ϕ π ν for ϕ [ 0 , π / ν ]
r ( ϕ ) = 1 n sin ( ϕ ) + cos ( ϕ ) + ε 2 ϕ ϕ + π ν for ϕ [ π / ν , 0 ]
with ϵ 1 ( ϵ 2 ) being the deformation parameter of the upper (lower) segment of the spike. Thus, the boundary of the deformed microstar is given by periodic continuation. Choosing ε 1 ε 2 allows for an asymmetric boundary deformation without a mirror reflection symmetry of the cavity.
The ray dynamics in the microstar are conveniently expressed in a Poincaré phase-space section, in the following phase space for short, which is introduced in Figure 1. Whenever a ray propagating along the spikes is passing the abscissa, its position x > 0 and the transversal momentum p = sin θ are recorded. For a microstar cavity without deformed boundary, i.e., ϵ 1 = ϵ 2 = 0 , a family of marginally stable orbits exists at p = 0 , whereas in the deformed case, an elliptical fixed point with a region of regular dynamics can be achieved. For a boundary deformation with a mirror-reflection symmetry ( ϵ 1 = ϵ 2 ), the fixed point is either at p = 0 or occurs pairwise at ± | p | . Contrary, for an asymmetric deformation ( ϵ 1 ϵ 2 ), a single elliptical fixed point can occur at a finite nonzero momentum, as in Figure 1c with ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) . For this deformation, a 5:1 resonance chain exists in the regular region, which is of further interest below. In comparison, the deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) lead to an island of regular dynamics with a similar shape but without a (macroscopic) resonance chain, see Figure 1d.
Note that the phase space represents CCW propagating rays. The corresponding rays propagating in CW direction are given by time reversal of the (geometrical) dynamics. For the phase-space variable p describing the momentum along the x-axis, this relates to a change in the sign, cf. Figure 1b. A mirror reflection, however, would not change the sign of p. Therefore, even in a cavity with a mirror-reflection symmetry, CW and CCW dynamics can have differences, namely in the direction along the symmetry line.
Figure 1. In (a), a microstar cavity with ν = 9 asymmetric spikes is shown as a black curve. The simple periodic orbit is shown by the green curve. In (b), a magnification of one spike is shown. Propagation in CCW (CW) direction is illustrated by red (blue) arrows. Figure (c) shows a phase-space portrait of the microstar with deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) , as in (a,b). The simple periodic orbit is shown as a green dot. A period-5 orbit is shown as magenta dots. (d) Phase-space portrait for deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) .
Figure 1. In (a), a microstar cavity with ν = 9 asymmetric spikes is shown as a black curve. The simple periodic orbit is shown by the green curve. In (b), a magnification of one spike is shown. Propagation in CCW (CW) direction is illustrated by red (blue) arrows. Figure (c) shows a phase-space portrait of the microstar with deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) , as in (a,b). The simple periodic orbit is shown as a green dot. A period-5 orbit is shown as magenta dots. (d) Phase-space portrait for deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) .
Entropy 24 01614 g001
In the following, we analyze the difference between CW and CCW propagation. Therefore, initial conditions ( x 0 , p 0 ) on a uniform grid in phase space with unit intensity I 0 = 1 are iterated N turns turns in CW and in CCW direction. The remaining fractions of the intensity I CCW in CCW direction are shown in Figure 2a–d for N turns = 1 , 15 , 50 , 100 , respectively. As can be seen, the loss introduced by partial back-reflections from Fresnel formulas is not uniform in the regular region, while the central fixed point suffers almost no loss due to Brewster-angle transmissions, the border of the regular region has more loss, which results in a loss gradient over the resonance chain. A consequence is a distortion of the intensity pattern, which can be seen in Figure 2b,c. This distortion is different for both propagation directions, resulting in a remarkable difference I CCW I CW around the resonance chain; see Figure 2e–h. Note that the effect becomes considerably weaker for long times where the only remaining intensity is centered around the central fixed point.
Intuitively, the difference between both propagation directions can be understood by analyzing the dynamics around the resonance chain, as shown in Figure 3. Focusing on one unstable orbit of the resonance chain (see red dot in Figure 3), the dynamics on the left side push the orbits closer to the center of the regular region, whereas on the right side, the orbits are pushed away from the regular region’s center. Hence, orbits on the left side of the unstable orbit will suffer less loss then orbits on the right side in the next iterations. This behavior is inverted when the dynamics are reversed from CCW to CW, which leads to the characteristic pattern of the intensity difference in Figure 2f–h.
For comparison, in Figure 4, the system without (macroscopic) resonance chain is investigated. Note that here the difference between intensities I CCW and I CW is much weaker and not structured in phase space.

Frobenius–Perron Operator

The Frobenius–Perron operator (FPO) describes the evolution of density (intensity) distributions in phase space [36]. It has been used for, e.g., maps [37,38,39,40,41], billiard systems [37], and the asymmetric backscattering in deformed microcavities [42]. Here, the aim is to construct classical intensity eigenstates for CCW and CW propagation in the deformed microstar via the FPO. For this purpose, Ulam’s method [43] is used to construct a matrix approximation F of the FPO as follows. The part of the phase space between x [ r ( ϕ / ν ) , r ( 0 ) ] and p [ 0.2 , 0.2 ] covering the regular island is divided into a lattice of N × N uniform rectangular cells. In the cell i, random initial conditions with total intensity I i are iterated once in (C)CW direction. The transferred intensity I f to the phase-space cell f is then recorded. Thus, the two N 2 × N 2 matrices with the elements F f i ( C ) C W = I f / I i are an approximation for the intensity transport in CW and CCW direction, respectively. As such, the eigenvalues are inside the unit circle, as shown in Figure 5a,b. Of special interest are the eigenstates of the FPO with the eigenvalue of largest magnitude, as they represent states with a long lifetime. For the system with a resonance chain, these long-lived eigenstates for (C)CW propagation are shown in Figure 5c–k, in addition to their differences. As it can be seen, the localization of the eigenstates around the central fixed point depends on the fineness of the discretization given by N. This is consistent with the previous literature on FPO eigenstates of maps with a mixed phase space [39,44]. For an infinite discretization, the intensity is localized on the stable periodic orbit, which is a single point in phase space. However, for a more rough discretization, i.e., N 200 , the eigenstates extend over the resonance chain leading to a local intensity distortion, which is different for CW and CCW eigenstates. Consequently, the difference between CW and CCW eigenstates is significant around the resonance chain and consistent with the ray dynamics shown in Figure 2f,g. In particular, the pattern with alternating preference of CW or CCW propagation along the stable and unstable points of the resonance chain is visible.
Note that the shown difference pattern between CW and CCW propagation is not due to an asymmetry of the cavity’s boundary. In fact, for deformation parameters ϵ 1 = ϵ 2 = 0.98 , the microstar has a mirror-reflection symmetry, a similar 5:1 resonance chain in the phase space, and a similar alternating pattern of the eigenstate difference (not shown).
In Figure 6, the eigenstates corresponding to the second, third, and fourth eigenvalues next to unity are shown. As such an eigenstate has positive and negative parts, it does not represent a proper intensity distribution by itself. However, for the dynamics of an arbitrary (positive) phase-space intensity on intermediate timescales, these eigenstates still provide useful information. Especially, for the second (Figure 6a,b) and fourth (Figure 6e,f) eigenstates, the chirality around the resonance chain as well as the difference between CW and CCW propagation can be seen clearly.
For a comparison, the eigenvalues and eigenstates of the FPO for the system without a resonance chain are shown in Figure 7. The eigenstates also localize at the center of the regular region, but they do not show a significant chiral behavior, i.e., the intensity distribution along the regular orbits is uniform.
Note that in contrast to ref. [42], here, the true-time mapping is not included in the FPO. However, this incorporation of different time scales for the iteration is not relevant here, as the rays propagating through the spikes of the microstar have almost the same optical path length regardless of whether they propagate with a large or small radius.

3. Wave Dynamics

In this section, the wave dynamics of the microstar are analyzed with a focus on the semiclassical regime and the difference between CW and CCW propagation. For quasi-two-dimensional microcavities, Maxwell’s equations reduce to the scalar mode Equation [45]
Δ + n 2 k 2 ψ = 0
where k = ω / c is the complex frequency which determines the Q-factor (lifetime) of a mode via Q = Re k / ( 2 Im k ) . To utilize Brewster’s angle, here, transverse electrical polarization is chosen such that the wave function ψ represents the magnetic field as H = ( 0 , 0 , Re [ ψ e i ω t ] ) . Hence, along the cavities interface ψ , and its scaled normal derivative n 2 v ψ are continuous. In addition, ψ is required to fulfill Sommerfeld’s outgoing wave condition. For the numerical calculation of the long-lived modes in the microstar, the boundary element method [46] and FEM software COMSOL [47] are used. In Figure 8, the long-lived modes in the complex frequency plane in the semiclassical regime with 350 Re k R 600 are shown. As in the wave regime discussed in ref. [31], the most long-lived mode localizes along the stable periodic orbit passing through the spikes of the microstar; see Figure 9d. However, additional long-lived modes exist in the semiclassical regime that also pass through the spikes but have more intensity maxima in radial direction; see Figure 9a–c. Thus, a classification of these modes with a radial mode number l is suitable.
To compare the wave and ray dynamics in phase space, the Husimi function H ( x , p ; ψ , ν ψ ) of a mode ψ is used. Here, we adapt the boundary Husimi function [48] to the phase space taken along the x-axis (see Figure 1). Thus, for x [ 0.5 R , 1.1 R ] , the Husimi function (for CCW propagation) is constructed via
H ( x , p ; ψ , ν ψ ) = 1 N k F h ψ ( x , p ) + i F h ν ψ ( x , p ) 2
where F = n 1 p 2 and σ = R / ( n k ) . The functions h ψ and h ν ψ are defined via an overlap with a Gaussian wave packet with
h f ( x , p ) = x m i n x m a x ζ ( t ; x , p ) f ( t ) d t
ζ ( t ; x , p ) = exp ( x t ) 2 2 σ 2 i p t .
The normalization constant N in Equation (6) is chosen such that max x , p H = 1 . Note that for the magnetic field, the normal derivative along the x-axis can be computed via ν ψ = y H z = i ω ε 0 ε r E x .
In Figure 9e–h, the Husimi functions of the corresponding modes are shown in the right panel. As an expected result, the Husimi function of the most long-lived mode localizes in the center of the regular island, see Figure 9h, with a maximum at the finite momentum p 0.01 , which nicely reflects the stable periodic orbit. Furthermore, as an expected result, the modes with higher radial mode numbers l have their Husimi functions localized on outer tori of the regular island. Hence, these states can be seen as higher-order quantized states on the regular region. However, in this regime of Re k R 404 , the modes cannot localize on the smaller sub-islands of the resonance chain itself. Consequently, the Husimi functions of the related modes interpolate through the resonance chain.
To verify the different intensity distortion for CCW and CW propagation, a few individual modes with a larger dimensionless frequency Re k R > 1200 are calculated. Then, for a mode, the Husimi function for CCW propagation H ( x , p ; ψ , ν ψ ) and the corresponding Husimi function for CW propagation H ( x , p ; ψ , ν ψ ) are calculated to obtain the difference.
I ( x , p ; ψ ) = H ( x , p ; ψ , ν ψ ) H ( x , p ; ψ , ν ψ ) .
Note that in Equation (9), each Husimi function is normalized as discussed above. Hence, a preferred sense of rotation in the mode ψ does not impact the difference I ( x , p ; ψ ) . In Figure 10a–i, the Husimi function and the difference I ( x , p ; ψ ) are shown for the long-lived modes in the microstar with a resonance chain. It is remarkable that over the range of k R 1200 , k R 1280 and k R 1360 , each mode and its partner mode with a quasi-degenerate frequency show almost the same pattern in I ( x , p ; ψ ) . This pattern is organized around the stable and unstable fixed point of the resonance chain and in a very good agreement with the pattern of the FPO eigenstates, e.g., in Figure 5e.
To highlight the importance of the resonance chain for this phenomenon, a mode pair in the microstar without macroscopic resonance chain is shown in Figure 10j–l. Already within the mode pair, the pattern I ( x , p ; ψ ) is different and not organized around an obvious phase-space structure.
A simplified model system that explains the effect in an open wave (quantum) system is discussed in the Appendix A.

4. Conclusions

In this article, we discussed the ray–wave correspondence in microstar cavities. From the ray dynamics, a (asymmetric) boundary deformation is used to manipulate the light around a periodic orbit with a regular island in phase space. Correspondingly, in the semiclassical regime, optical modes localize hierarchically on the regular island. Depending on the deformation of the microstar resonance, chains might appear in the regular island. Since the intensity loss is not uniform in the regular region, those resonance chains have an unusual effect on the dynamics. The propagation in CW and CCW direction shows an intensity distortion with a pattern that organizes along the resonance chain. This distortion is explained by the ray dynamics, and therefore manifests in the eigenstates of the FPO. Furthermore, the difference between CW and CCW propagation is revealed in the Husimi functions of the optical modes in the semiclassical regime.
Such an interplay between nonuniform loss and nonlinear dynamics is a rather general effect, as it can also be observed in a simple quantum model system. Hence, the effect is not restricted to microstar cavities. However, microstar cavities are ideal systems to study such dynamics, as a nontrivial connection between forward and backward dynamics in time manifests naturally in terms of CW and CCW propagation. Therefore, we believe our studies highlight microstar cavities as an interesting class of systems to unveil novel phenomena.

Author Contributions

Investigation, J.K.; Supervision, J.W.; Writing—original draft, J.K.; Writing—review and editing, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge support for the Book Processing Charge by the Open Access Publication Fund of Magdeburg University.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Difference in Forward and Backward Dynamics in a Generalized Pendulum Model

In this appendix, a simple model system is constructed to illustrate the interplay between a resonance chain and a loss gradient to cause a difference between two directions of dynamical propagation. In the microstar cavity, such two directions are given naturally by CW and CCW propagation. In the simplified one-dimensional model, the two directions are given by forward and backward propagation in time under Hamiltonian dynamics. The Hamiltonian is given by a generalized pendulum or normal form ansatz for a r:s resonance chain as [49,50,51]
H r : s ( θ , I ) = H 0 ( I ) + 2 V r : s I I r : s r 2 cos ( r θ ) .
Here, I 0 and θ [ 0 , 2 π ] are action-angle variables of the unperturbed part H 0 ( I ) = ( I I r : s ) 2 / ( 2 M r : s ) . The parameters are fixed to ( I r : s , M r : s , V r : s ) = ( 0.8 , 1 , 0.005 ) and r = 3 . Using a canonical transformation
M : ( θ , I ) ( Q , P ) = ( 2 I cos θ , 2 I sin θ ) ,
the Hamiltonian given in the coordinates ( θ , I ) can be related to a regular island in the phase space ( Q , P ) with an elliptical fixed point at the origin and a resonance chain around it; see Figure A1a,b.
Next, the Hamiltonian (A1) is quantized. Since the angle coordinate θ is 2 π -periodic, the action eigenstates | I n are a discrete basis of the Hilbert space with I n = ( n + 1 / 2 ) . Thus, the Hamiltonian (A1) has the matrix elements [52,53]
I m | H r : s | I n = ( I n I r : s ) 2 2 M r : s δ n m + V r : s I r : s r 2 n ! ( n r ) ! δ m , n r + ( n + r ) ! n ! δ m , n + r
and is therefore represented as an infinite matrix. In the numerical simulations, the reduced Planck’s constant is set to = 0.3 , and a cutoff for the action indices at N max = 90 is used. Hence, the (forward) time evolution operator U for over one unit time step is given as a finite matrix
U for = exp i H r : s .
Correspondingly, the operator for the dynamic one unit time step backward in time yields
U back = exp i H r : s .
Next, a nonuniform loss is introduced which has a gradient across the resonance chain and therefore mimics the situation in the microstar cavity. In the basis of action eigenstates, the loss is described by a reflection operator R with matrix elements
R ( I n ) = exp ( I n / I r : s )
in the diagonal, cf. Figure A1c. Hence, states with small I n suffer less loss than the states with large I n . The dynamics of the open system is therefore described by the sub-unitary time-evolution operator
U back , for ( open ) = R U back , for .
Note that U back ( open ) is not the time-reversed operator of U for ( open ) as loss is not turned into gain. This is consistent with the microstar cavity, where CCW is the time-reversed CW propagation, but both directions of the dynamics suffer loss from radiation and partial reflections. The eigenstates Ψ back , for ( open ) of U back , for ( open ) are calculated numerically. In particular, we are interested in the eigenstates with eigenvalues of the largest magnitude. Due to the loss profile, these eigenstates localize near I = 0 . Their Husimi function can be calculated as
H back , for ( θ , I ) = n = 0 N max α ( θ , I ) | I n I n | Ψ back , for ( open ) 2
with α ( θ , I ) being a minimal uncertainty wave packet at the phase-space point ( θ , I ) given by
α ( θ , I ) | I n = N ( I ) exp ( I n I ) 2 4 σ 2 + i ( I n I ) θ .
Here, σ = / 2 is a parameter and N ( I ) is a normalization factor ensuring α ( θ , I ) 2 = n = 0 N max | α ( θ , I ) | I n | 2 = 1 . The Husimi function for the long-lived eigenstates corresponding to the forward/backward dynamics are shown in Figure A1d,e. As it can be seen in the magnifications, the Husimi function’s contour lines are curved sinusoidally with a phase shift between forward and backward dynamics. Exactly this phase shift leads to the oscillatory pattern for the difference between both Husimi functions, which is organized along the stable and unstable fixed point of the resonance chain; see Figure A1f.
Figure A1. (a,b) Contour lines of the generalized pendulum Hamiltonian in the phase-space ( Q , P ) and ( θ , I ) are shown as black curves. The separatrix of the resonance chain is highlighted via a red curve. (c) The reflectivity function which defines the loss profile is shown via a thick black curve. Thin horizontal lines indicate the position of action eigenstates at I n . (d,e) The Husimi function for eigenstates of U for ( open ) ( U back ( open ) ) localizing predominantly at I 0 are shown via the red color gradient. Contour lines of the Husimi function are shown as red curves in the magnification plots. (f) Difference between the Husimi function for forward and backward propagation.
Figure A1. (a,b) Contour lines of the generalized pendulum Hamiltonian in the phase-space ( Q , P ) and ( θ , I ) are shown as black curves. The separatrix of the resonance chain is highlighted via a red curve. (c) The reflectivity function which defines the loss profile is shown via a thick black curve. Thin horizontal lines indicate the position of action eigenstates at I n . (d,e) The Husimi function for eigenstates of U for ( open ) ( U back ( open ) ) localizing predominantly at I 0 are shown via the red color gradient. Contour lines of the Husimi function are shown as red curves in the magnification plots. (f) Difference between the Husimi function for forward and backward propagation.
Entropy 24 01614 g0a1

References

  1. Vahala, K.J. Optical Microcavities; Advanced Series in Applied Physics; World Scientific: Singapore, 2004; Volume 5. [Google Scholar]
  2. Chang, R.K.; Campillo, A.J. Optical Processes in Microcavities; Advanced Series in Applied Physics; World Scientific: Singapore, 1996; Volume 3. [Google Scholar]
  3. Kippenberg, T.J.; Kalkman, J.; Polman, A.; Vahala, K.J. Demonstration of an erbium-doped microdisk laser on a silicon chip. Phys. Rev. A 2006, 74, 051802(R). [Google Scholar] [CrossRef] [Green Version]
  4. McCall, S.L.; Levi, A.F.J.; Slusher, R.E.; Pearton, S.J.; Logan, R.A. Whispering-gallery mode microdisk lasers. Appl. Phys. Lett. 1992, 60, 289–291. [Google Scholar] [CrossRef]
  5. Levi, A.F.J.; Slusher, R.E.; McCall, S.L.; Glass, J.L.; Pearton, S.J.; Logan, R.A. Directional light coupling from microdisk lasers. App. Phys. Lett. 1993, 62, 561–563. [Google Scholar] [CrossRef] [Green Version]
  6. Tamboli, A.C.; Haberer, E.D.; Sharma, R.; Lee, K.W.; Nakamura, S.; Hu, E.L. Room-temperature continuous-wave lasing in GaN/InGaN microdisks. Nat. Photonics 2007, 1, 61–64. [Google Scholar] [CrossRef]
  7. Ilchenko, V.S.; Gorodetsky, M.L.; Yao, X.S.; Maleki, L. Microtorus: A high-finesse microcavity with whispering-gallery modes. Opt. Lett. 2001, 26, 256–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Armani, D.K.; Kippenberg, T.J.; Spillane, S.M.; Vahala, K.J. Ultra-high-Q toroid microcavity on a chip. Nature 2003, 421, 925–928. [Google Scholar] [CrossRef]
  9. Zhang, X.; Choi, H.S.; Armani, A.M. Ultimate quality factor of silica microtoroid resonant cavities. Appl. Phys. Lett. 2010, 96, 153304. [Google Scholar] [CrossRef]
  10. Collot, L.; Lefevre-Seguin, V.; Brune, M.; Raimond, J.; Haroche, S. Very high-Q whispering-gallery modes observed on fused silica microspheres. Europhys. Lett. 1993, 23, 327–334. [Google Scholar] [CrossRef]
  11. Yang, J.; Moon, S.; Lee, S.B.; Lee, J.H.; An, K.; Shim, J.B.; Lee, H.W.; Kim, S.W. Development of a deformation-tunable quadrupolar microcavity. Rev. Sci. Instrum. 2006, 77, 083103. [Google Scholar] [CrossRef] [Green Version]
  12. Stöckmann, H.J. Quantum Chaos; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  13. Cao, H.; Wiersig, J. Dielectric microcavities: Model systems for wave chaos and non-Hermitian physics. Rev. Mod. Phys. 2015, 87, 61–111. [Google Scholar] [CrossRef]
  14. Nöckel, J.U.; Stone, A.D. Ray and wave chaos in asymmetric resonant optical cavities. Nature 1997, 385, 45–47. [Google Scholar] [CrossRef] [Green Version]
  15. Sunada, S.; Shinohara, S.; Fukushima, T.; Harayama, T. Signature of Wave Chaos in Spectral Characteristics of Microcavity Lasers. Phys. Rev. Lett. 2016, 116, 203903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Peng, B.; Özdemir, S.K.; Lei, F.; Monifi, F.; Gianfreda, M.; Long, G.L.; Fan, S.; Nori, F.; Bender, C.M.; Yang, L. Parity–time-symmetric whispering-gallery microcavities. Nat. Phys. 2014, 10, 394–398. [Google Scholar] [CrossRef] [Green Version]
  17. Lee, S.B.; Yang, J.; Moon, S.; Lee, S.Y.; Shim, J.B.; Kim, S.W.; Lee, J.H.; An, K. Observation of an Exceptional Point in a Chaotic Optical Microcavity. Phys. Rev. Lett. 2009, 103, 134101. [Google Scholar] [CrossRef] [Green Version]
  18. Zhu, J.; Özdemir, Ş.K.; He, L.; Yang, L. Controlled manipulation of mode splitting in an optical microcavity by two Rayleigh scatterers. Opt. Express 2010, 18, 23535–23543. [Google Scholar] [CrossRef] [Green Version]
  19. Peng, B.; Özdemir, Ş.K.; Liertzer, M.; Chen, W.; Kramer, J.; Yilmaz, H.; Wiersig, J.; Rotter, S.; Yang, L. Chiral modes and directional lasing at exceptional points. Proc. Natl. Acad. Sci. USA 2016, 113, 6845. [Google Scholar] [CrossRef] [Green Version]
  20. Richter, S.; Zirnstein, H.G.; Zúñiga Pérez, J.; Krüger, E.; Deparis, C.; Trefflich, L.; Sturm, C.; Rosenow, B.; Grundmann, M.; Schmidt-Grund, R. Voigt Exceptional Points in an Anisotropic ZnO-Based Planar Microcavity: Square-Root Topology, Polarization Vortices, and Circularity. Phys. Rev. Lett. 2019, 123, 227401. [Google Scholar] [CrossRef]
  21. Kippenberg, T.J.; Holzwarth, R.; Diddams, S.A. Microresonator-Based Optical Frequency Combs. Science 2011, 332, 555–559. [Google Scholar] [CrossRef]
  22. Nöckel, J.U.; Stone, A.D.; Chen, G.; Grossman, H.L.; Chang, R.K. Directional emission from asymmetric resonant cavities. Opt. Lett. 1996, 21, 1609–1611. [Google Scholar] [CrossRef] [Green Version]
  23. Wiersig, J.; Hentschel, M. Combining directional light output and ultralow loss in deformed microdisks. Phys. Rev. Lett. 2008, 100, 033901. [Google Scholar] [CrossRef]
  24. Wang, Q.J.; Yan, C.; Yu, N.; Unterhinninghofen, J.; Wiersig, J.; Pflügl, C.; Diehl, L.; Edamura, T.; Yamanishi, M.; Kan, H.; et al. Whispering-gallery mode resonators for highly unidirectional laser action. Proc. Natl. Acad. Sci. USA 2010, 107, 22407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Albert, F.; Hopfmann, C.; Eberspächer, A.; Arnold, F.; Emmerling, M.; Schneider, C.; Höfling, S.; Forchel, A.; Kamp, M.; Wiersig, J.; et al. Directional whispering gallery mode emission from Limaçon-shaped electrically pumped quantum dot micropillar lasers. Appl. Phys. Lett 2012, 101, 021116. [Google Scholar] [CrossRef]
  26. Armani, A.M.; Kulkarni, R.P.; Fraser, S.E.; Flagan, R.C.; Vahala, K.J. Label-Free, Single-Molecule Detection with Optical Microcavities. Science 2007, 317, 783–787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Vollmer, F.; Arnold, S. Whispering-gallery-mode biosensing: Label-free detection down to single molecules. Nat. Methods 2008, 5, 591. [Google Scholar] [CrossRef]
  28. Lu, T.; Lee, H.; Chen, T.; Herchak, S.; Kim, J.H.; Fraser, S.E.; Flagan, R.C.; Vahala, K. High sensitivity nanoparticle detection using optical microcavities. Proc. Nat. Acad. Sci. USA 2011, 108, 5976–5979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Sunada, S.; Harayama, T. Design of resonant microcavities: Application to optical gyroscopes. Opt. Express 2007, 15, 16245–16254. [Google Scholar] [CrossRef]
  30. Miao, P.; Zhang, Z.; Sun, J.; Walasik, W.; Longhi, S.; Litchinitser, N.M.; Feng, L. Orbital angular momentum microlaser. Science 2016, 353, 464–467. [Google Scholar] [CrossRef] [Green Version]
  31. Kullig, J.; Jiang, X.; Yang, L.; Wiersig, J. Microstar cavities: An alternative concept for the confinement of light. Phys. Rev. Res. 2020, 2, 012072. [Google Scholar] [CrossRef] [Green Version]
  32. Kullig, J.; Wiersig, J. Microdisk cavities with a Brewster notch. Phys. Rev. Res. 2021, 3, 023202. [Google Scholar] [CrossRef]
  33. Wiersig, J.; Kim, S.W.; Hentschel, M. Asymmetric scattering and nonorthogonal mode patterns in optical microspirals. Phys. Rev. A 2008, 78, 053809. [Google Scholar] [CrossRef]
  34. Wiersig, J.; Eberspächer, A.; Shim, J.B.; Ryu, J.W.; Shinohara, S.; Hentschel, M.; Schomerus, H. Nonorthogonal pairs of copropagating optical modes in deformed microdisk cavities. Phys. Rev. A 2011, 84, 023845. [Google Scholar] [CrossRef] [Green Version]
  35. Wiersig, J. Non-Hermitian Effects Due to Asymmetric Backscattering of Light in Whispering-Gallery Microcavities. In Parity-Time Symmetry and Its Applications; Christodoulides, D., Yang, J., Eds.; Springer: Singapore, 2018; pp. 155–184. [Google Scholar]
  36. Beck, C.; Schögl, F. Thermodynamics of Chaotic Systems: An Introduction; Cambridge Nonlinear Science Series; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  37. Altmann, E.G.; Portela, J.S.E.; Tél, T. Leaking chaotic systems. Rev. Mod. Phys. 2013, 85, 869–918. [Google Scholar] [CrossRef] [Green Version]
  38. Frahm, K.M.; Shepelyansky, D.L. Ulam method for the Chirikov standard map. Eur. Phys. J. B 2010, 76, 57–68. [Google Scholar] [CrossRef]
  39. Weber, J.; Haake, F.; Šeba, P. Frobenius-Perron Resonances for Maps with a Mixed Phase Space. Phys. Rev. Lett. 2000, 85, 3620–3623. [Google Scholar] [CrossRef] [Green Version]
  40. Ermann, L.; Shepelyansky, D.L. Ulam method and fractal Weyl law for Perron-Frobenius operators. Eur. Phys. J. B 2010, 75, 299–304. [Google Scholar] [CrossRef]
  41. Carlo, G.G.; Rivas, A.M.F.; Spina, M.E. Classical to quantum correspondence in dissipative directed transport. Phys. Rev. E 2015, 92, 052907. [Google Scholar] [CrossRef] [Green Version]
  42. Kullig, J.; Wiersig, J. Frobenius–Perron eigenstates in deformed microdisk cavities: Non-Hermitian physics and asymmetric backscattering in ray dynamics. New J. Phys. 2016, 18, 015005. [Google Scholar] [CrossRef]
  43. Ulam, S.M. A Collection of Mathematical Problems; Interscience Publishers: New York, NY, USA, 1960. [Google Scholar]
  44. Weber, J.; Haake, F.; Braun, P.A.; Manderfeld, C.; Seba, P. Resonances of the Frobenius-Perron operator for a Hamiltonian map with a mixed phase space. J. Phys. A Math. Gen. 2001, 34, 7195–7211. [Google Scholar] [CrossRef]
  45. Jackson, J.D. Classical Electrodynamics; John Wiley and Sons: New York, NY, USA, 1962. [Google Scholar]
  46. Wiersig, J. Boundary element method for resonances in dielectric microcavities. J. Opt. A Pure Appl. Opt. 2003, 5, 53. [Google Scholar] [CrossRef]
  47. COMSOL Multiphysics® v. 6.0; COMSOL AB: Stockholm, Sweden, 2022; Available online: www.comsol.com (accessed on 2 November 2022).
  48. Hentschel, M.; Schomerus, H.; Schubert, R. Husimi functions at dielectric interfaces: Inside-outside duality for optical systems and beyond. Europhys. Lett. 2003, 62, 636–642. [Google Scholar] [CrossRef]
  49. Lebœuf, P.; Mouchet, A. Normal Forms and Complex Periodic Orbits in Semiclassical Expansions of Hamiltonian Systems. Ann. Phys. 1999, 275, 54–112. [Google Scholar] [CrossRef] [Green Version]
  50. Lichtenberg, A.J.; Liebermann, M.A. Regular and Stochastic Motion; Springer: New York, NY, USA, 1983. [Google Scholar]
  51. Keshavamurthy, S.; Schlagheck, P.E. Dynamical Tunneling: Theory and Experiment, 1st ed.; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  52. Löck, S.; Bäcker, A.; Ketzmerick, R.; Schlagheck, P. Regular-to-Chaotic Tunneling Rates: From the Quantum to the Semiclassical Regime. Phys. Rev. Lett. 2010, 104, 114101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Mertig, N.; Kullig, J.; Löbner, C.; Bäcker, A.; Ketzmerick, R. Perturbation-free prediction of resonance-assisted tunneling in mixed regular-chaotic systems. Phys. Rev. E 2016, 94, 062220. [Google Scholar] [CrossRef] [PubMed]
Figure 2. The intensity distribution ( x 0 , p 0 , I CCW ) of rays with initial conditions ( x 0 , p 0 , I 0 = 1 ) after 1 , 15 , 50 , 100 turns in CCW direction is shown in (ad). Correspondingly, (eh) show the distribution ( x 0 , p 0 , I CCW I CW ) of the difference between propagation in CCW and CW direction. Note the adjusted scale in the colormap for (eh).
Figure 2. The intensity distribution ( x 0 , p 0 , I CCW ) of rays with initial conditions ( x 0 , p 0 , I 0 = 1 ) after 1 , 15 , 50 , 100 turns in CCW direction is shown in (ad). Correspondingly, (eh) show the distribution ( x 0 , p 0 , I CCW I CW ) of the difference between propagation in CCW and CW direction. Note the adjusted scale in the colormap for (eh).
Entropy 24 01614 g002
Figure 3. Magnification of the phase-space portrait from Figure 1c around the 5:1 resonance chain. Iterations for selected (nonperiodic) orbits are shown as dots colored from black to orange. The unstable period-5 orbit is represented by dots with a cyclic colormap. The arrows illustrate the dynamics around one of the unstable fixed points. The deformation parameters of the microstar are ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) .
Figure 3. Magnification of the phase-space portrait from Figure 1c around the 5:1 resonance chain. Iterations for selected (nonperiodic) orbits are shown as dots colored from black to orange. The unstable period-5 orbit is represented by dots with a cyclic colormap. The arrows illustrate the dynamics around one of the unstable fixed points. The deformation parameters of the microstar are ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) .
Entropy 24 01614 g003
Figure 4. Similar as in Figure 2, (ad) show the intensity distributions ( x 0 , p 0 , I CCW ) and (eh) show the distributions ( x 0 , p 0 , I CCW I CW ) , but for deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) , where no macroscopic resonance chain occurs in the phase space. Note the adjusted scale in the colormap for (eh).
Figure 4. Similar as in Figure 2, (ad) show the intensity distributions ( x 0 , p 0 , I CCW ) and (eh) show the distributions ( x 0 , p 0 , I CCW I CW ) , but for deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) , where no macroscopic resonance chain occurs in the phase space. Note the adjusted scale in the colormap for (eh).
Entropy 24 01614 g004
Figure 5. (a) The 10,000 largest-magnitude eigenvalues λ of the FPO with N = 200 are shown in the complex plane. Orange squares (blue crosses) represent the eigenvalues of the CCW (CW) FPO. Note that they are often on top of each other. (b) Shows a magnification of (a) around unity. (ce) show the eigenstates corresponding to the eigenvalue of the largest modulus for CCW direction, CW direction, and the difference between both eigenstates for N = 100 . The same sequence is shown in (fh) for N = 200 and in (ik) for N = 400 .
Figure 5. (a) The 10,000 largest-magnitude eigenvalues λ of the FPO with N = 200 are shown in the complex plane. Orange squares (blue crosses) represent the eigenvalues of the CCW (CW) FPO. Note that they are often on top of each other. (b) Shows a magnification of (a) around unity. (ce) show the eigenstates corresponding to the eigenvalue of the largest modulus for CCW direction, CW direction, and the difference between both eigenstates for N = 100 . The same sequence is shown in (fh) for N = 200 and in (ik) for N = 400 .
Entropy 24 01614 g005
Figure 6. Eigenstates of the FPO ( N = 200 ) with (a,b) second, (c,d) third, and (e,f) fourth eigenvalues next to unity. The left (right) panel shows the eigenstates for CCW (CW) propagation.
Figure 6. Eigenstates of the FPO ( N = 200 ) with (a,b) second, (c,d) third, and (e,f) fourth eigenvalues next to unity. The left (right) panel shows the eigenstates for CCW (CW) propagation.
Entropy 24 01614 g006
Figure 7. (a) Eigenvalues close to unity of the FPO ( N = 200 ) for the system without a resonance chain characterized by ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) (cf. Figure 4). (b) Eigenstate of the FPO for the CCW dynamics with largest-modulus eigenvalue. (c) Difference between CCW and CW eigenstate. (df) Next eigenstates of the FPO in CCW direction with eigenvalues close to unity.
Figure 7. (a) Eigenvalues close to unity of the FPO ( N = 200 ) for the system without a resonance chain characterized by ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) (cf. Figure 4). (b) Eigenstate of the FPO for the CCW dynamics with largest-modulus eigenvalue. (c) Difference between CCW and CW eigenstate. (df) Next eigenstates of the FPO in CCW direction with eigenvalues close to unity.
Entropy 24 01614 g007
Figure 8. Dimensionless frequencies k R of the long-lived modes in the microstar cavity that can be classified by their mode number in radial direction by (blue) l = 1 , (red) l = 2 , (green) l = 3 , (black) l = 4 , (yellow) l = 5 , and (cyan) l = 6 . The gray line indicates the frequency of the modes shown in Figure 9.
Figure 8. Dimensionless frequencies k R of the long-lived modes in the microstar cavity that can be classified by their mode number in radial direction by (blue) l = 1 , (red) l = 2 , (green) l = 3 , (black) l = 4 , (yellow) l = 5 , and (cyan) l = 6 . The gray line indicates the frequency of the modes shown in Figure 9.
Entropy 24 01614 g008
Figure 9. (ad) shows the intensity pattern of the modes with radial mode number l = 4 to l = 1 and Re k R 404 (see gray line in Figure 8). (eh) represents the Husimi projection of the corresponding mode.
Figure 9. (ad) shows the intensity pattern of the modes with radial mode number l = 4 to l = 1 and Re k R 404 (see gray line in Figure 8). (eh) represents the Husimi projection of the corresponding mode.
Entropy 24 01614 g009
Figure 10. (a) Husimi function for CCW propagation, (b) difference I ( x , p ; ψ ) between Husimi functions for CCW and CW propagation, see Equation (9). (c) Same as (b) but for a partner mode with almost the same frequency k R as in (b). The panels (di) follow the same sequence as (a,b) but for modes with different frequency k R in the same system given by deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) . Modes for the microstar given by ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) are shown in (jl).
Figure 10. (a) Husimi function for CCW propagation, (b) difference I ( x , p ; ψ ) between Husimi functions for CCW and CW propagation, see Equation (9). (c) Same as (b) but for a partner mode with almost the same frequency k R as in (b). The panels (di) follow the same sequence as (a,b) but for modes with different frequency k R in the same system given by deformation parameters ( ϵ 1 , ϵ 2 ) = ( 0.76 , 1.2 ) . Modes for the microstar given by ( ϵ 1 , ϵ 2 ) = ( 0.78 , 1.212 ) are shown in (jl).
Entropy 24 01614 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kullig, J.; Wiersig, J. Ray–Wave Correspondence in Microstar Cavities. Entropy 2022, 24, 1614. https://doi.org/10.3390/e24111614

AMA Style

Kullig J, Wiersig J. Ray–Wave Correspondence in Microstar Cavities. Entropy. 2022; 24(11):1614. https://doi.org/10.3390/e24111614

Chicago/Turabian Style

Kullig, Julius, and Jan Wiersig. 2022. "Ray–Wave Correspondence in Microstar Cavities" Entropy 24, no. 11: 1614. https://doi.org/10.3390/e24111614

APA Style

Kullig, J., & Wiersig, J. (2022). Ray–Wave Correspondence in Microstar Cavities. Entropy, 24(11), 1614. https://doi.org/10.3390/e24111614

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop