Bloch theorem dictated wave chaos in microcavity crystals

Universality class of wave chaos emerges in many areas of science, such as molecular dynamics, optics, and network theory. In this work, we generalize the wave chaos theory to cavity lattice systems by discovering the intrinsic coupling of the crystal momentum to the internal cavity dynamics. The cavity-momentum locking substitutes the role of the deformed boundary shape in the ordinary single microcavity problem, providing a new platform for the in situ study of microcavity light dynamics. The transmutation of wave chaos in periodic lattices leads to a phase space reconfiguration that induces a dynamical localization transition. The degenerate scar-mode spinors hybridize and non-trivially localize around regular islands in phase space. In addition, we find that the momentum coupling becomes maximal at the Brillouin zone boundary, so the intercavity chaotic modes coupling and wave confinement are significantly altered. Our work pioneers the study of intertwining wave chaos in periodic systems and provide useful applications in light dynamics control.

ing substitutes the role of the deformed boundary shape in the ordinary single microcavity problem, providing a new platform for the in-situ study of microcavity light dynamics.The transmutation of wave chaos in periodic lattices leads to a phase space reconfiguration that induces a dynamical localization transition.
The degenerate scar-mode spinors hybridize and non-trivially localize around regular islands in phase space.In addition, we find that the momentum coupling becomes maximal at the Brillouin zone boundary, so the intercavity chaotic modes coupling and wave confinement are significantly altered.Our work pioneers the study of intertwining wave chaos in periodic systems and provide useful applications in light dynamics control.
The boundary deformation (BD) in microcavity has been considered the most common approach for phase space engineering to yield the desired optical properties.The main idea is based on the observation that BD sensitively reconfigures the underlying phase space transportation (30,31).Alternative approaches also have been proposed, such as tailoring phase space by defect (32) or including a circular hole inside the cavity (33,34).Despite these comprehensive efforts, previous approaches still require a holistic device change at the fabrication stage, which remains the main obstacle to further rapid progress.In this regard, developing a new platform on which the light dynamics can be studied by in-situ experimental control is genuinely appealing.
Here, we propose a two-dimensional periodic lattice structure consisting of multiple chaotic microcavities as a promising breakthrough.We reveal that the external crystal momentum coupling to the internal cavity dynamics can take over the role of the BD by breaking and restoring the inherent symmetry of the cavity.Both the direction and the amplitude of the external momentum can be readily controlled by steering the coherent excitation light sources.As a result, cavity lattice systems promise a hold for feasible control of wave chaos features such as dynamical localization and tunneling (35)(36)(37)(38)(39).
In this work, we explore chaotic signatures in periodic lattice dictated by the Bloch theorem for the first time.We study the adiabatic change of the internal cavity states due to the crystal momentum coupling as where Ψ tot (r) and k represent the total wave function over the lattice and the crystalline momentum, respectively.Here, ψ Int (r) describes the internal dynamics of the cavity states, which we focus on.The coupling between the internal phase space dynamics and the crystal momentum is observed explicitly.In particular, the momentum-induced coupling leads to the hybridizations of the wave function from the scar states to the regular orbit states; namely, we observe the crystal momentum-induced dynamical localization transition of the internal dynamics.On the lattice scale, we find that introducing the additional deformation leads to the finite Berry curvature and Hall effect triggered by skew scattering events.Our work firstly promotes applications to topological optical transport utilizing chaotic states.Finally, we discuss the possibility of extending our studies to various lattice systems.
We consider a square lattice consisting of a single dielectric cavity with a relative refractive index, n in /n out = 10, inside the cavity [see Fig. 1].With this high refractive index, we can find energetically well-isolated target modes by suppressing undesired overlaps with other extra modes.
However, our result is generally applicable to the broad range of the refractive index.The boundary of each cavity is determined by the four-fold rotational symmetric (C 4 ) deformation, which can be represented as, r(θ; ε) = r 0 [1 − ε cos(4θ)], where ε denotes the deformation strength.r 0 = R/ 1 + ε 2 /2, where R is the radius of the circle when ε = 0, is the normalization constant preserving a cavity area under the variation of ε.The Helemholtz equation, −∇ 2 ψ = n(x, y) ω c ψ, is solved to obtain resonant frequencies ω of the transverse-magnetic (TM)-polarized modes [ ψ = (0, 0, E z )], where c is a speed of light.For TM-modes, the waves and their normal derivative are set to be continuous across the cavity boundary.The two lattice vector, where we impose a periodic boundary condition accordingly, is given by |a To numerically obtain the modes both in a single cavity and in a lattice, we implement the boundary element method (40)(41)(42)(43).
iv Under a weak BD, invariant tori in the ray-dynamical phase space of the integrable system are destroyed according to the Kolmogorov-Arnold-Moser (KAM) and Poincaré-Birkhoff theorems (44,45).As the deformation increases, the regular and chaotic regions fill the phase space before it becomes fully chaotic and ergodic.Figure 2 shows the phase space in the Birkhoff- . Here, q/L is normalized boundary arclength where the ray bounces off (L is a cavity perimeter), and χ is the incident angle of the ray.The mixed phase space in the figure shows the island structures of regular motions surrounded by a chaotic sea.See Supplementary Materials for details.

Results
Two degenerate scar modes (47), the non-trivial localization on the unstable fixed points, are observed in the single cavity case for ωR/2πc ≈ 0.375 at = 0.05.A series of the thick blue curves in Fig. 1c show the successive regular and scar modes that a similar sequence has been analyzed in (48).The scar modes form a pair under the underlying C 4 -symmetry.In the lattice, with the identical cavity geometry, the spectra at zero momentum exhibit very similar energy eigenstates (Fig. 1d).We observe the pair of scar modes equally for almost the same energy as the single cavity case.Figures 2a and b show the Husimi distributions of the two degenerate scar states, which are localized on top of the distinct unstable fixed points.These fixed points correspond to period-2 (i.e., two bounces for one cycle of a classical orbit) unstable periodic orbits shown in the leftmost panels in Fig. 2. The black curves in Fig. 2 show the stable and unstable manifolds of the period-2 unstable orbit, and the unstable fixed points corresponding to this orbit are the crossing points of those manifolds.Note that the Husimi distributions are obtained for the cavity boundary wave, and the inside-incident version was employed (49,50).
We now consider the inclusion of the finite crystal momentum k = (k x , k y ).In the presence of the non-zero crystal momentum, we observe the degeneracy lifting of the two scar modes.This lifting indicates the state hybridization due to the C 4 -rotational symmetry breaking of the internal wave function with the non-zero momentum.The hybridized states are found to localize prominently around the period-4 bow-tie orbit (see the rightmost panels in Fig. 2), the satellite orbit of the central period-2 island.Figures 2c and d exemplify the Husimi distributions for these hybridized modes localized on the bow-tie orbits.The mode evolution depending on the crystalline momentum variations in cavity lattices can be equivalently realized in single cavities by applying the additional deformation perturbation.See Supplementary Materials for detailed demonstrations.
In addition, Figs.3b-e show the hybridized states that have strong directional anisotropy in the direction of the crystal momentum.To be specific, the hybridization of the two scar modes is numerically obtained and described by the effective Hamiltonian, given as (51), where σ represent the Pauli matrices for the scar state degree of freedom.represents the hybridization strength between the two scar modes, while V 0 = 0.213 × 10 −3 represents the relative energy shift of the two states, respectively.Since the two scar states are energetically well-decoupled from other states, we can consider the two scar states as the effective spinor (|↑ and |↓ ).Under this spinor representation, the spinor wind twice around the loop encircling the zero momentum (Fig. 3 a).This winding of the spinor gives rise to 2π-Berry phase: , where |ψ(k) indicates one of the scar states and i = √ −1.
The degeneracy of the spinor scar states manifests as the topological quadratic band touchings (QBT) protected by C 4 -symmetry in the momentum space (see Fig. 4a) (52).Here, the mode hybridization due to the crystalline momentum is not restricted in the scar modes and is valid for the general modes that exhibit a degeneracy at time-reversal invariant momenta.
viii At the X point in BZ, there exist only two possibilities that the resonant energy of the regular state can adiabatically deform.First, even parity regular state, ψ E (x, y) (Fig. 4b) is not compatible with the anti-periodic boundary condition unless the wave function vanishes at the unit-cell boundary.

Maximal momentum coupling
The vanishing wave function manifests as the flat band with zero group velocity (see the upper band along the X-M interval in Fig. 4a).In addition, since ψ E (x, y) is confined more densely inside the cavity, where the refractive index is higher, the effective wavelength is reduced, i.e., even parity states form higher energy states (point "i" in Fig. 4a).
On the other hand, the odd parity regular state, ψ O (x, y) (Fig. 4c), shows the compensating intensity enhancement at the unit-cell boundary since the even and odd wave functions must form the complete basis of the two original scar modes.Hence, it results in lengthened effective wavelength (i.e., more waves are outside the cavity, the lower refractive index region) and manifests as the lower energy resonant states (point "ii" in Fig. 4a).Finally, when the momentum vector reaches the BZ corner [M = (π, π)], the global C 4 -rotational symmetry is restored.In this case, the hybridization disappears to form the degenerate QBT again, which can be interpreted as a revived dynamical localization of the spinor scar modes.
where j = x, y.
Since the optical microcavity array can be excited by the coherent beam source with a welldefined incident momentum, the Berry curvature dipole can be measured through the skew-symmetric beam transport.The semi-classical transport equation ( 55) describes the light dynamics of the Mie regime, and the non-zero Berry dipole manifests as the effective magnetic field as, where v k is the group velocity vector of the scar states having energy E k , and n(r) is the slowly varying local refractive index.The second term ( k × Ω xy ) in Eq. ( 4) gives rise to skew symmetric light transport that can be captured by control of the incident beam momentum k.The explicit demonstration of the skew symmetric light transportation can be found in Supplementary Materials.

Discussion
To conclude, we have studied the wave chaos of the deformed cavity coupled to the external crystalline momentum in a periodic cavity array.The external crystalline momentum now replaces the role of the boundary shape deformation.By controlling the momentum, we observe the dynamical localization transitions.Our work provides a new promising platform, enabling the in-situ study of various wave chaos phenomena.For example, if the additional higher energy and lower energy states are involved, it can induce crystal momentum-induced dynamical tunneling phenomena.
This will be the topic of future study.
In addition, we note that, contrary to the previous studies of the topological photonic crystals

Method
We solve the Helmholtz wave equation, which is deduced from Maxwell's equations without sources, for optical modes in single cavities and photonic crystals, where n(r) and ω = ck are, respectively, the refractive index of the piecewise homogeneous medium and the free-space temporal frequency with vacuum wavenumber k and speed of light  5).In addition to the dielectric boundary conditions, we impose a two-dimensional periodic condition at the unitcell boundary of the two-dimensional periodic lattice.A pure outgoing wave condition at infinity is applied for the modes in the single cavities.To solve Eq. ( 5) numerically, we employ the boundary element method (BEM) (40,41) and further implement the block Sakurai-Sugiura method (42,43) to compute the optical modes more efficiently. xiii

Figure 1 :
Figure 1: a Conceptual illustration of a photonic crystal consisting of a deformed dielectric microcavities.b Schematic diagram of a square lattice unit-cell with a lattice constant a. c Real part resonant frequencies in a single cavity as a function of the cavity deformation ε. d Energy eigenvalues in a cavity array lattice as a function of ε at zero crystalline momentum k = 0 [Γ in Fig. 4a].In c and d, successive Demkov-type couplings for the stable regular modes and the scar modes are highlighted by thick curves.The arrows at ε = 0.05 mark the degenerated scar modes (solid and dashed curves) we study.

Figure 2 :Figure 3 :
Figure 2: Husimi functions for modes in a cavity lattice superimposed on the ray-dynamical classical phase space (underlying gray dots) in Birkhoff-coordinate (q/L, p = sin χ), where q/L is a arclength normalized by a cavity perimeter L and χ a incident angle of ray.In a-d, the four large closed islands and the crossing points of stable and unstable manifolds (black curves) at p ∼ 0 correspond period-2 stable and unstable orbits, respectively.The orange dots in c and d represent the phase space points of the bow-tie orbit.a and b show the Husimi distributions for the scar states obtained when k = (k x , k y ) = (0, 0), while c and d are for the bow-tie orbit states induced by couplings of the two scar states when k = (π/2, 0).The leftmost and rightmost panels show corresponding classical orbits in real space.The four vertical lines in a-d indicate a quarter arclength interval L/4.When k = (0, 0) (a, b), the phase space (e.g., islands and manifolds) exhibits a L/4 translation symmetry, and the Husimi functions in a and b are distributed mutual exclusively.On the other hand, when k = (π/2, 0) (c, d), the translation symmetry of the phase space breaks, and the Husimi distributions in c and d have large overlaps around the vertical lines.Detailed correspondences between the Bloch momenta k and the phase space deformation can be found in Supplementary Materials.The continuous evolution of the Husimi function and the underlying phase space over the range [π/4 ≤ tan −1 (k y /k x ) ≤ 9π/4] for |k| = π/2 is shown in the supplementary animations: "Ani supple phc hus1 series.mov"and "Ani supple phc hus2 series.mov".

3 viiFigure 4 :
Figure 4: a Energy bands as a function of the crystalline momentum (k x , k y ) for a fixed cavity deformation ε = 0.05.The zoomed bands correspond to the ones indicated by the arrow in Fig. 1d.b Spatial distribution of the even parity mode (ψ E ) in a unit-cell for (k x , k y )=(π, 0) (i in a).The upper panel depicts the wave projection onto the x-axis showing the vanished wave at the unitcell boundary and dense confinement of wave inside the cavity region (shaded rectangle).c The same as b but for the odd parity mode (ψ O ) (ii in a).The upper panel depicts the intense wave distribution outside the cavity, including the unit-cell boundary.
When the momentum becomes maximal near the Brillouin zone (BZ) boundary [X = (π, 0) and Y = (0, π)], the wavefunction characteristics in the whole unit-cell domain deviate significantly due to the change in the boundary condition.For instance, at X point, the internal wave function satisfies the anti-periodic boundary condition on the unit-cell boundary: ψ Int (x + a, y) = −ψ Int (x, y) and ψ Int (x, y + a) = ψ Int (x, y).The linear combination of the two scar states forms two regular states depending on its parity, [ψ(x, y) = ±ψ(−x, y)], due to C 4 -rotational symmetry.

Figure 5 :
Figure 5: Schematic illustration of BDs a original C 4 -symmetric boundary with QBT and b C 4breaking boundary.C 4 -breaking boundary perturbation splits the single QBT into a pair of Dirac cones.Subsequently, C 2 -breaking boundary perturbation gaps out the Dirac cones and induces the Berry curvature dipole.The Berry curvature dipole induces the skew scattering depending on the incident momentum, k inc .c Distribution of non-zero Berry curvature in the momentum space induced by scar states.

(
which mainly focus on the rayleigh regime), the governing dynamics of the chaotic state in the xi Mie regime is described by the semi-classical transport.In this work, we propose, for the first time, the possibility of realizing Berry curvature-induced transport phenomena that utilize the intrinsic wave property of the chaotic states.The crossover between Rayleigh and Mie regime of Berry curvature induced transport would pioneer a new aspect of wave-particle correspondence in the field of wave chaos.Our work can be generalized in a few different aspects.Though we have considered the simple square lattice, many other two-dimensional lattices are expected to show qualitatively different behaviors.For instance, away from a simple Bravais lattice, the deformed cavity in Lieb and Kagome lattice is expected to produce a flat band.The intrinsic localization properties would give rise to a stronger coupling within a unit-cell.Further studies on the various lattice systems and wave chaos would also be an intriguing topic for future study.

c.
Given the cavity boundary shape, R(θ; ε) = R 0 [1 − ε cos(4θ)], where R 0 = ρ/ 1 + ε 2 /2, we set n(r) = 10 for |r| < R(θ; ε) and n(r) = 1 otherwise.ρ is the radius of a circle when ε = 0.The lattice constant is given as |a x | = |a y | = a = 2.2ρ.We focus on the transverse magnetic [TM; ψ = (0, 0, E z )] polarization of modes.The TM-polarized modes fulfill the dielectric boundary condition that E z and ν • ∇E z are continuous across the boundary interface of two different refractive index domains.Here, ν is an outward normal vector of the boundary.Note xii the transverse-electric [TE; ψ = (0, 0, H z )] modes satisfy a different boundary condition: H z and 1/n 2 ν • ∇H z are continuous, yet, still the solution of the same wave equation, Eq. (