Optical Theory Group

Abstract : Resonances of open non-Hermitian systems are associated with the poles of the system scattering matrix. Perturbations of the system cause these poles to shift in the complex frequency plane. In this work, we introduce a novel method for calculating shifts in scattering matrix poles using generalized Wigner-Smith operators. We link our method to traditional cavity perturbation theory and validate its effectiveness through application to complex photonic networks. Our findings underscore the versatility of generalized Wigner-Smith operators for analyzing a broad spectrum of resonant systems and provides new insight into resonant properties of non-Hermitian systems.

Abstract : Random matrix theory is a useful tool in the study of the physics of multiple scattering systems, often striking a balance between computation speed and physical rigour. Propagation of waves through thick disordered media, as arises in for example optical scattering or electron transport, typically necessitates cascading of multiple random matrices drawn from an underlying ensemble for thin media, greatly increasing computational burden. Here we propose a dual pool based bootstrapping approach to speed up statistical studies of scattering in thick random media. We examine how potential matrix reuse in this approach can impact statistical estimates of population averages. Specifically, we discuss how both bias and additional variance in the sample mean estimator are introduced through bootstrapping. In the diffusive scattering regime, the extra estimator variance is shown to originate from samples in which cascaded transfer matrices are permuted matrix products. Through analysis of the combinatorics and cycle structure of permutations we quantify the resulting correlations. Proofs of several analytic formulae enumerating the frequency with which correlations of different strengths occur are derived. Extension to the ballistic regime is briefly considered.

Abstract : In this work we present a method for generating random matrices describing electromagnetic scattering from disordered media containing dielectric particles with prescribed single particle scattering characteristics. Resulting scattering matrices automatically satisfy the physical constraints of unitarity, reciprocity and time reversal, whilst also incorporating the polarization properties of electromagnetic waves and scattering anisotropy. Our technique therefore enables statistical study of a variety of polarization phenomena, including depolarization rates and polarization-dependent scattering by chiral particles. In this vein, we perform numerical simulations for media containing isotropic and chiral spherical particles of different sizes for thicknesses ranging from the single to multiple scattering regime and discuss our results, drawing comparisons to established theory.