Ghostly Singularities in Light: How Boundaries Shape Wave Behavior

Author: Denis Avetisyan


New research reveals that the edges of photonic structures can create unusual ‘exceptional points’-singularities in light propagation-without needing energy loss or gain.

The study demonstrates that within a defect waveguide, variations in the odd index <span class="katex-eq" data-katex-display="false">n_0</span> induce a transition in light intensity decay, shifting from rapid dissipation at the exceptional point to emergent oscillatory dynamics as <span class="katex-eq" data-katex-display="false">n_0</span> surpasses a value of 5, a phenomenon corroborated by both numerical solutions of coupled-mode equations and an approximate analytical expression.
The study demonstrates that within a defect waveguide, variations in the odd index n_0 induce a transition in light intensity decay, shifting from rapid dissipation at the exceptional point to emergent oscillatory dynamics as n_0 surpasses a value of 5, a phenomenon corroborated by both numerical solutions of coupled-mode equations and an approximate analytical expression.

A semi-infinite photonic lattice with a side-coupled defect demonstrates boundary-induced delay as a mechanism for generating non-Hermitian singularities and non-Markovian dynamics.

While non-Hermitian physics typically relies on engineered gain, loss, or reservoirs, this work-’Boundary-Driven Exceptional Points in Photonic Waveguide Lattices’-demonstrates the surprising emergence of exceptional points within a conservative photonic system. Specifically, we predict and analyze how coherent reflections at the boundary of a semi-infinite photonic lattice with a side-coupled defect induce non-Markovian dynamics and give rise to memory-driven exceptional points. Our exact analytic treatment reveals that these singularities arise from boundary-induced delay, offering a tunable platform for exploring memory effects without external non-Hermitian elements. Could this boundary-driven approach unlock new avenues for realizing and controlling non-Hermitian phenomena in a wider range of conservative systems?


Beyond Static Control: Embracing the Imperfections of Light

Conventional photonic lattices, structures designed to guide and manipulate light, traditionally model light propagation using Bloch modes – solutions to the Schrödinger equation assuming a perfectly periodic potential. However, this approach breaks down when confronted with the intricacies of real-world systems exhibiting complex interactions such as absorption, amplification, or disorder. These imperfections introduce non-reciprocity and asymmetry, invalidating the assumptions behind Bloch modes and hindering accurate predictions of light behavior. The standard framework struggles to account for phenomena like localized states arising from loss, or the unusual directional sensitivity observed in certain lattice designs. Consequently, a more sophisticated mathematical description is needed to capture the full range of possibilities within these increasingly complex photonic environments, paving the way for exploration of non-Hermitian physics.

Photonic systems described by non-Hermitian Hamiltonians represent a significant departure from traditional wave mechanics, offering a means to model scenarios where energy is not conserved due to the presence of gain and loss. Unlike conventional Hermitian systems where eigenvalues remain real, these Hamiltonians allow for complex eigenvalues, directly influencing the propagation characteristics of light. This leads to unconventional wave dynamics, including phenomena such as non-reciprocal light transport, where light travels differently depending on the direction, and the emergence of exceptional points – singularities in the parameter space where waves undergo dramatic changes. The ability to engineer these gains and losses within photonic lattices opens avenues for novel devices with enhanced sensitivity, unidirectional invisibility, and tailored wave manipulation, exceeding the limitations of standard Bloch-based approaches. These systems are not merely theoretical curiosities; they provide a powerful platform for exploring fundamental physics and developing advanced photonic technologies.

The unusual behavior of light within non-Hermitian photonic lattices arises from a delicate interplay between fundamental symmetries – parity and time reversal. Typically, these symmetries ensure predictable wave propagation; however, introducing gain and loss deliberately breaks one or both, opening pathways for unprecedented control. When parity-time (PT) symmetry is preserved despite these imbalances, the system exhibits real-valued eigenvalues, allowing for lossless power oscillation and guiding of light along unconventional paths. Conversely, breaking PT symmetry leads to complex eigenvalues, manifesting as amplified or attenuated modes and the formation of exceptional points – singularities where the system’s response dramatically changes. Exploiting these symmetry-breaking scenarios offers the potential to create highly sensitive sensors, unidirectional invisibility cloaks, and novel optical switches, effectively transforming the way light is manipulated and utilized in photonic devices.

A semi-infinite waveguide lattice with a side-coupled defect, modeled by a single-level Fano-Anderson system, enables analysis of resonant states via complex contour integration with a Bromwich path deformed into Hankel contours and a branch cut along the imaginary axis.
A semi-infinite waveguide lattice with a side-coupled defect, modeled by a single-level Fano-Anderson system, enables analysis of resonant states via complex contour integration with a Bromwich path deformed into Hankel contours and a branch cut along the imaginary axis.

Fano Resonance: A Signature of Interaction

Fano resonance occurs when a discrete quantum state interacts with a continuum of states, resulting in a characteristic asymmetric lineshape in the system’s response. This interference phenomenon deviates from the standard Lorentzian lineshape observed in isolated resonances; instead, the lineshape exhibits a spectral asymmetry, often characterized by a divergent slope on one side. The asymmetry arises because the continuum allows for energy exchange even at energies away from the discrete state’s resonance frequency, effectively broadening and distorting the spectral response. The strength of this interaction, and therefore the degree of asymmetry, is determined by the coupling between the discrete state and the continuum. \Gamma_{rad} represents the decay rate into the continuum, while \Gamma_{iso} represents the isolated resonance width, and their ratio dictates the lineshape’s asymmetry.

The Fano-Anderson model, originally developed for scattering theory in atomic physics, provides a quantitative description of the interaction between discrete quantum states and a continuum of states, and is particularly well-suited to analyze side-coupled defect waveguides. In these waveguides, the discrete state is represented by the localized mode within the defect, while the continuum arises from the propagating modes in the adjacent waveguides. The model describes the total transmission or reflection as a function of the energy, incorporating a coupling parameter Γ that quantifies the strength of interaction between the discrete and continuum states. By fitting experimental transmission spectra to the Fano profile, researchers can extract parameters characterizing the defect’s coupling strength and lifetime, providing insights into waveguide design and functionality. The model’s applicability stems from its ability to account for the interference effects arising from these coupled states, resulting in the characteristic asymmetric lineshape observed in the transmission spectrum.

Within the Fano-Anderson model, the self-energy \Sigma(E) quantifies the influence of the continuum on the discrete state, and its analysis is central to characterizing the coupling strength. Calculation of \Sigma(E) often involves the Lambert W function, a multi-valued complex function defined as the inverse of f(w) = w e^w . This arises from solving the equations derived from the model, specifically when determining the energy dependence of the coupling between the discrete level and the continuum. The real part of the self-energy contributes to the broadening of the discrete state, while the imaginary part dictates the shift in its resonant energy. Precise determination of these components, enabled by the Lambert W function, allows for quantitative assessment of the coupling parameters and provides insight into the interaction mechanism.

The behavior of resonance poles <span class="katex-eq" data-katex-display="false">s_1</span> and <span class="katex-eq" data-katex-display="false">s_2</span> shifts with increasing odd values of <span class="katex-eq" data-katex-display="false">n_0</span> and normalized coupling <span class="katex-eq" data-katex-display="false">g/J</span>, demonstrating that approximate poles calculated via Lambert functions closely match exact numerical computations up to the boundary-driven non-Markovian exceptional point where eigenvalues cross.
The behavior of resonance poles s_1 and s_2 shifts with increasing odd values of n_0 and normalized coupling g/J, demonstrating that approximate poles calculated via Lambert functions closely match exact numerical computations up to the boundary-driven non-Markovian exceptional point where eigenvalues cross.

Memory and Delay: The Emergence of System History

Optical waveguide systems, due to the finite time required for light to traverse their length, demonstrate non-instantaneous responses resulting in memory effects. This means the current state of signal propagation is not solely determined by the present input, but is also dependent on the history of excitation within the waveguide. Specifically, photons traveling within the system experience a delay related to the waveguide’s characteristics; previously excited photons continue to propagate and interact with subsequent inputs. This contrasts with Markovian systems where future behavior is independent of past states. The manifestation of these memory effects is observable in phenomena like pulse broadening and spectral distortion, and necessitates modeling approaches that account for the temporal evolution of the electromagnetic field within the waveguide.

Non-Markovian dynamics within optical waveguide systems are mathematically represented by differential-delayed equations. These equations extend standard differential equations by incorporating terms that depend on the field at previous times. Specifically, the current rate of change is not solely determined by the instantaneous field value, but also by the field’s history, represented by a time-delayed component. This delayed component accounts for the finite time it takes for information to propagate through the system and influence the current state. The general form includes both instantaneous feedback terms-standard in Markovian descriptions-and integral terms representing the influence of past excitation, effectively modeling memory effects within the propagation process. \frac{d}{dt}A(t) = F[A(t), A(t-\tau)] , where A(t) is the field at time t, τ is the delay, and F represents the system’s functional dependence.

The characteristic delay time τ within an optical waveguide system is determined by the ratio of the initial number of excitations n_0 to the coupling strength J, expressed as \tau = n_0/J. This parameter governs the temporal separation between an excitation and its delayed response, directly influencing the system’s memory effects. Abrupt changes in the waveguide environment, such as the introduction of a new element or alteration of system parameters, can be mathematically represented using the Heaviside step function. This function models instantaneous transitions, allowing for the simulation of step-change inputs and the analysis of the system’s response to discontinuous environmental modifications.

Relaxation dynamics, focusing on the final stage, demonstrate how increasing the coupling strength <span class="katex-eq" data-katex-display="false">g_1 = \alpha g</span> modulates the system's behavior across a range of α values from 0 to 0.4.
Relaxation dynamics, focusing on the final stage, demonstrate how increasing the coupling strength g_1 = \alpha g modulates the system’s behavior across a range of α values from 0 to 0.4.

Mathematical Tools for Untangling Temporal Behavior

The Laplace transform converts differential equations into algebraic equations, simplifying the solution process, particularly for systems exhibiting time delays. This is achieved by transforming time-domain functions, f(t), into the frequency domain, F(s), where ‘s’ is a complex frequency variable. Delay terms, represented as f(t-\tau), where τ is the delay, are transformed into e^{-s\tau} in the Laplace domain. This algebraic manipulation allows for solutions to be found using techniques like partial fraction decomposition and inverse Laplace transforms, effectively handling the complexities introduced by time delays in dynamic systems. The method is particularly valuable in control systems, signal processing, and areas where initial conditions and transient responses are critical.

The inverse Laplace transform, used to determine the time-domain function from its Laplace-domain representation, frequently necessitates complex integration along a contour in the complex plane. The Bromwich integral, a common technique for evaluating this inverse transform, involves integrating the function F(s) multiplied by e^{st} along a vertical line in the complex plane, often denoted by \sigma = \text{Re}(s). The contour is typically chosen to ensure the integral converges and encompasses any poles of F(s) that lie to the right of the contour. Closing the contour, often with a semicircle of infinite radius, allows application of the residue theorem to compute the integral and thus determine the inverse Laplace transform. The location and residues of the poles significantly influence the resulting time-domain function.

Solutions to certain time-delay differential equations, when obtained via Laplace transforms, can be expressed using special functions such as the Bessel function of the first kind, denoted as J_n(x), where n is the order of the function and x is the argument. The appearance of these functions indicates an inherent symmetry or oscillatory behavior within the system being modeled. Specifically, Bessel functions arise frequently when dealing with problems involving cylindrical coordinates or radial symmetry, and their presence in the solution suggests that the time-delay effects are coupled with such spatial or rotational characteristics. The order of the Bessel function is often directly related to the magnitude of the delay, providing a quantitative link between the delay parameter and the solution’s symmetry.

Harnessing Memory for Next-Generation Photonics

Photonic lattices, structures that guide and manipulate light, are gaining enhanced capabilities through the incorporation of memory effects. This approach moves beyond static light control, enabling devices to ‘remember’ and respond to past optical signals, thereby creating more complex and dynamic functionalities. By strategically introducing these memory elements-often achieved through nonlinear interactions or feedback loops within the lattice-researchers are developing systems capable of advanced signal processing, optical switching, and even computation. These advancements promise a shift from passive photonic components to active, intelligent devices, potentially revolutionizing fields like optical communications and sensing, as the lattice’s behavior is no longer solely determined by the present input but also by its recent history.

Delay-coupled semiconductor lasers present a versatile architecture for sculpting light into intricate patterns and generating signals beyond the capabilities of single lasers. These systems exploit the feedback loop created by delaying and re-injecting a portion of the laser’s output, fostering interactions that give rise to complex behaviors like self-organization and chaos. This approach allows researchers to move beyond static optical fields, creating dynamic photonic structures and waveforms with potential applications in optical computing, sensing, and secure communication. By carefully controlling the delay time and coupling strength, it becomes possible to engineer the laser’s output, generating novel optical signals characterized by unique temporal and spatial properties-a significant advancement in manipulating light for advanced technologies.

Recent research highlights the creation of boundary-driven exceptional points (EPs) within a photonic lattice, representing a critical juncture where conventional optical behavior breaks down. These EPs, characterized by a specific condition of (g/J)^2 * n_0, signify a point of enhanced sensitivity and accelerated decay of light. The study demonstrates that at a coupling strength of g/J \approx 0.35, light experiences its fastest possible decay rate, offering potential for novel optical devices and enhanced light-matter interactions. This precise control over light decay, enabled by the engineered photonic lattice and the emergence of EPs, opens pathways for applications ranging from advanced sensing technologies to the development of more efficient optical switches and modulators.

The research highlights how exceptional points arise not from designed non-Hermiticity, but from the system’s inherent boundary conditions and resulting delays. This echoes a deeper principle: robustness emerges, it cannot be designed. The study demonstrates that complex behavior – in this case, non-Hermitian singularities – doesn’t necessitate external control or engineered dissipation. As Isaac Newton observed, “We build too many walls and not enough bridges.” The photonic lattice, constrained by its semi-infinite nature, generates these singularities internally, suggesting that system structure is stronger than individual control. Memory effects, arising from the boundary, become integral to the system’s behavior, a consequence of its architecture rather than imposed regulation.

The Horizon of Influence

The demonstration that boundary conditions alone can sculpt non-Hermitian behavior within a conservative system subtly shifts the discourse. The search for engineered reservoirs of gain and loss, for actively driven dissipation, may have been a misplaced emphasis. Order doesn’t need architects; it emerges from local rules – in this case, the simple physics of reflection and delay at an edge. The observed exceptional points aren’t imposed, but revealed by the lattice’s own limits.

However, the semi-infinite geometry, while elegant, presents a clear constraint. The forest evolves without a forester, yet follows rules of light and water – but this forest has a definite, and somewhat artificial, edge. Future work will undoubtedly explore the influence of more complex boundaries, finite-size effects, and disorder. Can these seemingly benign perturbations amplify the memory effects, or will they wash them out, restoring the predictable flow of energy?

The study highlights the often-overlooked role of the observer – or, more precisely, the system’s self-awareness of its own limits. The lattice ‘remembers’ through its boundaries. The challenge now lies in understanding how to harness this intrinsic memory, not to control it, but to influence its emergence and shape the resulting landscape of exceptional points. Control is an illusion; influence is real.


Original article: https://arxiv.org/pdf/2603.07212.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

See also:

2026-03-11 03:17