skip to main content


Title: Wave Excitation and Dynamics in Non-Hermitian Disordered Systems
Dynamic and steady-state aspects of wave propagation are deeply connected in lossless open systems ‎in which the scattering matrix is unitary. There is then an equivalence among the energy excited within ‎the medium through all channels, the Wigner time delay, which is the sum of dwell times in all ‎channels coupled to the medium, and the density of states. But these equivalences fall away in the ‎presence of material loss or gain. In this paper, we use microwave measurements, numerical ‎simulations, and theoretical analysis to discover the changing relationships among fundamental wave ‎properties with loss and gain, and their dependence upon dimensionality and spectral overlap. We ‎begin with the demonstrations that the transmission time in random 1D media is equal to the density ‎of states even in the presence of ultrastrong absorption and that its ensemble average is independent ‎of the strengths of scattering and absorption. In contrast, the Wigner time becomes imaginary in the ‎presence of loss, with real and imaginary parts that fall with absorption. In multichannel media, the ‎transmission time remains equal to the density of states and is independent of the scattering strength ‎in unitary systems but falls with absorption to a degree that increases with the strengths of absorption ‎and scattering, and the number of channels coupled to the medium. We show that the relationships ‎between key propagation variables in non-Hermitian systems can be understood in terms of the ‎singularities of the phase of the determinant of the transmission matrix. The poles of the transmission ‎matrix are the same as those of the scattering matrix, but the transmission zeros are fundamentally ‎different. Whereas the zeros of the scattering matrix are the complex conjugates of the poles, the ‎transmission zeros are topological: in unitary systems they occur only singly on the real axis or as ‎conjugate pairs. We follow the evolution and statistics of zeros in the complex plane as random ‎samples are deformed. The sensitivity of the spacing of zeros in the complex plane with deformation ‎of the sample has a square-root singularity at a zero point at which two single zeros and a complex ‎pair interconvert. The transmission time is a sum of Lorentzian functions associated with poles and ‎zeros. The sum over poles is the density of states with an average that is independent of scattering ‎and dissipation. But the sum over zeros changes with loss, gain, scattering strength and the number of ‎channels in ways that make it possible to control ultranarrow spectral features in transmission and ‎transmission time. We show that the field, including the contribution of the still coherent incident ‎wave, is a sum over modal partial fractions with amplitudes that are independent of loss and gain. The ‎energy excited may be expressed in terms of the resonances of the medium and is equal to the dwell ‎time even in the presence of loss or gain.‎  more » « less
Award ID(s):
2022629
NSF-PAR ID:
10319445
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physical review research
Volume:
4
Issue:
1
ISSN:
2643-1564
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Understanding vanishing transmission in Fano resonances in quantum systems and metamaterials and perfect and ultralow transmission in disordered media has advanced the knowledge and applications of wave interactions. Here, we use analytic theory and numerical simulations to understand and control the transmission and transmission time in complex systems by deforming a medium and adjusting the level of gain or loss. Unlike the zeros of the scattering matrix, the position and motion of the zeros of the determinant of the transmission matrix (TM) in the complex plane of frequency and field decay rate have robust topological properties. In systems without loss or gain, the transmission zeros appear either singly on the real axis or as conjugate pairs in the complex plane. As the structure is modified, two single zeros and a complex conjugate pair of zeros may interconvert when they meet at a square root singularity in the rate of change of the distance between the transmission zeros in the complex plane with sample deformation. The transmission time is the spectral derivative of the argument of the determinant of the TM. It is a sum over Lorentzian functions associated with the resonances of the medium, which is the density of states, and with the zeros of the TM. Transmission vanishes, and the transmission time diverges as zeros are brought near the real axis. Monitoring the transmission and transmission time when two zeros are close may open new possibilities for ultrasensitive detection. 
    more » « less
  2. Sirko, Leszek (Ed.)
    We demonstrate the ability to control the scattering properties of a two-dimensional wave-chaotic microwave billiard through the use of tunable metasurfaces located on the interior walls of the billiard. The complex reflection coefficient of the metasurfaces can be varied by applying a DC voltage bias to varactor diodes on the mushroom-shaped resonant patches, and this proves to be very effective at perturbing the eigenmodes of the cavity. Placing multiple metasurfaces inside the cavity allows us to engineer desired scattering conditions, such as coherent perfect absorption, by actively manipulating the poles and zeros of the scattering matrix through the application of multiple voltage biases. We demonstrate the ability to create on-demand coherent perfect absorption conditions at a specific frequency, and document the near-null of output power as a function of four independent parameters tuned through the coherent perfect absorption point. A remarkably low output-to-input power ratio P_{out}/P_{in} = 3:71 x 10^{-8} is achieved near the coherent perfect absorption point at 8.54 GHz. 
    more » « less
  3. Numerical difficulties associated with computing matrix elements of operators between Hartree–Fock–Bogoliubov (HFB) wavefunctions have plagued the development of HFB-based many-body theories for decades. The problem arises from divisions by zero in the standard formulation of the nonorthogonal Wick’s theorem in the limit of vanishing HFB overlap. In this Communication, we present a robust formulation of Wick’s theorem that stays well-behaved regardless of whether the HFB states are orthogonal or not. This new formulation ensures cancellation between the zeros of the overlap and the poles of the Pfaffian, which appears naturally in fermionic systems. Our formula explicitly eliminates self-interaction, which otherwise causes additional numerical challenges. A computationally efficient version of our formalism enables robust symmetry-projected HFB calculations with the same computational cost as mean-field theories. Moreover, we avoid potentially diverging normalization factors by introducing a robust normalization procedure. The resulting formalism treats even and odd number of particles on equal footing and reduces to Hartree–Fock as a natural limit. As proof of concept, we present a numerically stable and accurate solution to a Jordan–Wigner-transformed Hamiltonian, whose singularities motivated the present work. Our robust formulation of Wick’s theorem is a most promising development for methods using quasiparticle vacuum states.

     
    more » « less
  4. null (Ed.)
    Abstract We outline and interpret a recently developed theory of impedance matching or reflectionless excitation of arbitrary finite photonic structures in any dimension. The theory includes both the case of guided wave and free-space excitation. It describes the necessary and sufficient conditions for perfectly reflectionless excitation to be possible and specifies how many physical parameters must be tuned to achieve this. In the absence of geometric symmetries, such as parity and time-reversal, the product of parity and time-reversal, or rotational symmetry, the tuning of at least one structural parameter will be necessary to achieve reflectionless excitation. The theory employs a recently identified set of complex frequency solutions of the Maxwell equations as a starting point, which are defined by having zero reflection into a chosen set of input channels, and which are referred to as R-zeros. Tuning is generically necessary in order to move an R-zero to the real frequency axis, where it becomes a physical steady-state impedance-matched solution, which we refer to as a reflectionless scattering mode (RSM). In addition, except in single-channel systems, the RSM corresponds to a particular input wavefront, and any other wavefront will generally not be reflectionless. It is useful to consider the theory as representing a generalization of the concept of critical coupling of a resonator, but it holds in arbitrary dimension, for arbitrary number of channels, and even when resonances are not spectrally isolated. In a structure with parity and time-reversal symmetry (a real dielectric function) or with parity–time symmetry, generically a subset of the R-zeros has real frequencies, and reflectionless states exist at discrete frequencies without tuning. However, they do not exist within every spectral range, as they do in the special case of the Fabry–Pérot or two-mirror resonator, due to a spontaneous symmetry-breaking phenomenon when two RSMs meet. Such symmetry-breaking transitions correspond to a new kind of exceptional point, only recently identified, at which the shape of the reflection and transmission resonance lineshape is flattened. Numerical examples of RSMs are given for one-dimensional multimirror cavities, a two-dimensional multiwaveguide junction, and a multimode waveguide functioning as a perfect mode converter. Two solution methods to find R-zeros and RSMs are discussed. The first one is a straightforward generalization of the complex scaling or perfectly matched layer method and is applicable in a number of important cases; the second one involves a mode-specific boundary matching method that has only recently been demonstrated and can be applied to all geometries for which the theory is valid, including free space and multimode waveguide problems of the type solved here. 
    more » « less
  5. The non-Hermitian models, which are symmetric under parity (P) and time-reversal (T) operators, are the cornerstone for the fabrication of new ultra-sensitive optoelectronic devices. However, providing the gain in such systems usually demands precise control of nonlinear processes, limiting their application. In this paper, to bypass this obstacle, we introduce a class of time-dependent non-Hermitian Hamiltonians (not necessarily Floquet) that can describe a two-level system with temporally modulated on-site potential and couplings. We show that implementing an appropriate non-Unitary gauge transformation converts the original system to an effective one with a balanced gain and loss. This will allow us to derive the evolution of states analytically. Our proposed class of Hamiltonians can be employed in different platforms such as electronic circuits, acoustics, and photonics to design structures with hiddenPT-symmetry potentially without imaginary onsite amplification and absorption mechanism to obtain an exceptional point.

     
    more » « less