We present 4D topological textures in (quasi)monochromatic nonparaxial optical lattices that contain all possible polarization ellipses with every combination of ellipticity and orientation in 3D space. These fields span the nonparaxial polarization space (a complex projective plane) and a 4-sphere within specific spatiotemporal regions, forming 4D skyrmionic structures. Constructed from five plane waves with adiabatically varying relative amplitudes, they are experimentally realizable in free space by focusing a temporally variant beam with a high numerical aperture lens.

We demonstrate pixelation-free real-time widefield endoscopic imaging through an aperiodic multicore fiber (MCF) without any distal opto-mechanical elements or proximal scanners. Exploiting the memory effect in MCFs the images in our system are directly obtained without any post-processing using a static wavefront correction obtained from a single calibration procedure. Our approach allows for video-rate 3D widefield imaging of incoherently illuminated objects with imaging speed not limited by the wavefront shaping device refresh rate.

Polarization scramblers are essential for many imaging applications involving polarization sensitive instruments and partially polarized fluxes. In such cases, the light must be depolarized to allow properly calibrated measurements. Several types of depolarizers are already in use, but none is optimal due to the inevitable image degradation associated with the scrambling process. Here, we present a device based on an all-dielectric metasurface using anisotropic scatterers capable of generating multiple polarization states by varying their orientation angle. Our new scrambling solution allows a massive reduction in the integrated degree of polarization and thus the spatial depolarization of any incident linear polarization, while allowing easier integration into the instrument design and reducing the impact on its image quality.

Phase retrieval is an inverse problem that, on one hand, is crucial in many applications across imaging and physics, and, on the other hand, leads to deep research questions in theoretical signal processing and applied harmonic analysis. This survey paper is an outcome of the recent workshop Phase Retrieval in Mathematics and Applications (PRiMA) (held on August 5--9 2024 at the Lorentz Center in Leiden, The Netherlands) that brought together experts working on theoretical and practical aspects of the phase retrieval problem with the purpose to formulate and explore essential open problems in the field.

We describe a phase-retrieval-based imaging method to directly spatially resolve the vector lattice distortions in an extended crystalline sample by explicit coupling of independent Bragg ptychography data sets into the reconstruction process. Our method addresses this multi-peak Bragg ptychography (MPBP) inverse problem by explicit gradient descent optimization of an objective function based on modeling of the probe-lattice interaction, along with corrective steps to address spurious reconstruction artifacts. Robust convergence of the optimization process is ensured by computing exact gradients with the automatic differentiation capabilities of high-performance computing software packages. We demonstrate MPBP reconstruction with simulated ptychography data mimicking diffraction from a single crystal membrane containing heterogeneities that manifest as phase discontinuities in the diffracted wave. We show the superior ability of such an optimization-based approach in removing reconstruction artifacts compared to existing phase retrieval and lattice distortion reconstruction approaches.

Resonances are common in wave physics and their full and rigorous characterization is crucial to correctly tailor the response of a system in both time and frequency domains. However, they have been conventionally described by the quality factor, a real-valued number quantifying the sharpness of a single peak in the amplitude spectrum, and associated with a singularity in the complex frequency plane. But the amplitude of a physical signal does not hold all the information on the resonance and it has not been established that even the knowledge of the full distribution of singularities carries this information. Here we derive a dimensionless quality function that fully characterizes resonances from the knowledge of the phase spectrum of the signal. This function is driven by the spectral derivative of the phase. It is equivalent to the imaginary part of the Wigner-Smith time delay but it has the advantage of being valid for arbitrary response functions, including the components of the S-matrix. The spectral derivative of the phase can be calculated numerically from simulations or experimental acquisitions of the phase spectrum. Alternatively, it can be retrieved from the distribution of poles and zeros in the complex frequency plane through an analytic expression, which demonstrates that singularities do not suffice to fully characterize the resonances and that both singularities and zeros must be taken into account to retrieve the quality function. This approach permits to extract all the characteristics of resonances from arbitrary spectral response functions without a priori knowledge on the physical system.

The lensless endoscope is a promising device designed to image tissues in vivo at the cellular scale. The traditional acquisition setup consists in raster scanning during which the focused light beam from the optical fiber illuminates sequentially each pixel of the field of view (FOV). The calibration step to focus the beam and the sampling scheme both take time. In this preliminary work, we propose a scanning method based on compressive sampling theory. The method does not rely on a focused beam but rather on the random illumination patterns generated by the single-mode fibers. Experiments are performed on synthetic data for different compression rates (from 10 to 100% of the FOV).

The classical solution to the Helmholtz wave equation in spherical coordinates is well known and has found many important applications in wave propagation, scattering, and imaging in optics and acoustics. The separable solution is comprised of spherical Bessel functions in the radial direction and spherical harmonics in the angular directions. The nature of the spherical Bessel functions includes a long asymptotic oscillatory tail at large radii, not conducive to applications where a tight concentration of wave amplitude around a ring is desired, for example in toroidal configurations. However, we have found that certain practical bandpass spectral shapes, centered around a peak frequency, can create a superposition of spherical Bessel functions that effectively concentrate the wave amplitude around a defined ring at the time instant of coherent addition, avoiding the long tail asymptotic oscillations of the single frequency solution. Theoretical solutions are shown for different bandpass spectra applied to the spherical Bessel functions, along with numerical solutions of transient wave propagation using practical hemispherical source shapes. These findings introduce a framework by which ring or toroidal concentrated waves can be produced with a simple bandpass superposition applied to hemispherical source shapes and with reference to the classical solutions in spherical coordinates.

Fiber Fabry--Perot (FFP) resonators of a few centimeters are optimized as a function of the reflectivity of the mirrors and the dimensions of the intra-cavity waveguide. Loaded quality factor in excess of 10^9, with an optimum of 4___x___10^9, together with an intrinsic quality factor larger than 10^10 and intrinsic finesse in the range of 10^5 have been measured. An application to the stabilization of laser frequency fluctuations is presented.

Speckle based imaging consists of forming a super-resolved reconstruction of an unknown sample from low-resolution images obtained under random inhomogeneous illuminations (speckles). In a blind context where the illuminations are unknown, we study the intrinsic capacity of speckle-based imagers to recover spatial frequencies outside the frequency support of the data, with minimal assumptions about the sample. We demonstrate that, under physically realistic conditions, the covariance of the data has a super-resolution power corresponding to the squared magnitude of the imager point spread function. This theoretical result is important for many practical imaging systems such as acoustic and electromagnetic tomographs, fluorescence and photoacoustic microscopes, or synthetic aperture radar imaging. A numerical validation is presented in the case of fluorescence microscopy.

The Singularity Expansion Method Parameter Optimizer - SEMPO - is a toolbox to extract the complex poles, zeros and residues of an arbitrary response function acquired along the real frequency axis. SEMPO allows to determine this full set of complex parameters of linear physical systems from their spectral responses only, without prior information about the system. The method leverages on the Singularity Expansion Method of the physical signal. This analytical expansion of the meromorphic function in the complex frequency plane motivates the use of the Cauchy method and auto-differentiation-based optimization approach to retrieve the complex poles, zeros and residues from the knowledge of the spectrum over a finite and real spectral range. Both approaches can be sequentially associated to provide highly accurate reconstructions of physical signals in large spectral windows. The performances of SEMPO are assessed and analysed in several configurations that include the dielectric permittivity of materials and the optical response spectra of various optical metasurfaces.

The Lorenz-Mie scattering of a wide class of focused electromagnetic fields off spherical particles is studied. The focused fields in question are constructed through complex focal displacements, leading to closed-form expressions that can exhibit several interesting physical properties, such as orbital and/or spin angular momentum, spatially-varying polarization, and a controllable degree of focusing. These fields constitute complete bases that can be considered as nonparaxial extensions of the standard Laguerre-Gauss beams and the recently proposed polynomials-of-Gaussians beams. Their analytic form turns out to lead also to closed-form expressions for their multipolar expansion. Such expansion can be used to compute the field scattered by a spherical particle and the resulting forces and torques exerted on it, for any relative position between the field's focus and the particle.

In this second chapter, we analyse transmission problems between a dielectric and a dispersive negative material. In the first part, we consider a transmission problem between two half-spaces, filled respectively by the vacuum and a Drude material, and separated by a planar interface. In this setting, we answer to the following question: does this medium satisfy a limiting amplitude principle? This principle defines the stationary regime as the large time asymptotic behavior of a system subject to a periodic excitation. In the second part, we consider the transmission problem of an infinite strip of Drude material embedded in the vacuum and analyse the existence and dispersive properties of guided waves. In both problems, our spectral analysis enlighten new and unusual physical phenomena for the considered transmission problems due to the presence of the dispersive negative material. In particular, we prove the existence of an interface resonance in the first part and the existence of slow light phenomena for guiding waves in the second part.

Light interaction with optical cavities is of fundamental interest to enhance the light-matter interaction and to shape the spectral features of the electromagnetic fields. Important efforts have been carried out to develop modal theories of open optical cavities relying on an expansion of the fields on the eigen-fields of the cavity. Here, we show how such an expansion predicts the temporal dynamics of optical resonators. We consider a Fabry-Perot cavity to derive the full analytical expressions of the internal and scattered field on the quasi-normal modes basis together with the complex eigen-frequencies. We evince the convergence and accuracy of this expansion before deriving the impulse response function (IRF) of the open cavity. We benefit from this modal expansion and IRF to demonstrate that the eigen-modes of the open cavity impact the signals only during the transient regimes and not in the permanent regime.

We propose that wave propagation through a class of mechanical metamaterials opens unprecedented avenues in seismic wave protection based on spectral properties of auxetic-like metamaterials. The elastic parameters of these metamaterials like the bulk and shear moduli, the mass density, and even the Poisson ratio, can exhibit negative values in elastic stop bands. We show here that the propagation of seismic waves with frequencies ranging from 1Hz to 40Hz can be influenced by a decameter scale version of auxetic-like metamaterials buried in the soil, with the combined effects of impedance mismatch, local resonances and Bragg stop bands. More precisely, we numerically examine and illustrate the markedly different behaviors between the propagation of seismic waves through a homogeneous isotropic elastic medium (concrete) and an auxetic-like metamaterial plate consisting of 64 cells (40mx40mx40m), utilized here as a foundation of a building one would like to protect from seismic site effects. This novel class of seismic metamaterials opens band gaps at frequencies compatible with seismic waves when they are designed appropriately, what makes them interesting candidates for seismic isolation structures.

Resonances, also known as quasi normal modes (QNM) in the non-Hermitian case, play an ubiquitous role in all domains of physics ruled by wave phenomena, notably in continuum mechanics, acoustics, electrodynamics, and quantum theory. In this paper, we present a QNM expansion for dispersive systems, recently applied to photonics but based on sixty year old techniques in mechanics. The resulting numerical algorithm appears to be physically agnostic, that is independent of the considered physical problem and can therefore be implemented as a mere toolbox in a nonlinear eigenvalue computation library.

A general process is proposed to experimentally design anisotropic inhomogeneous metamaterials obtained through a change of coordinate in the Helmholtz equation. The method is applied to the case of a cylindrical transformation that allows to perform cloaking. To approximate such complex metamaterials we apply results of the theory of homogenization and combine them with a genetic algorithm. To illustrate the power of our approach, we design three types of cloaks composed of isotropic concentric layers structured with three types of perforations: curved rectangles, split rings and crosses. These cloaks have parameters compatible with existing technology and they mimic the behavior of the transformed material. Numerical simulations have been performed to qualitatively and quantitatively study the cloaking efficiency of these metamaterials.

Lensless illumination single-pixel imaging with a multicore fiber (MCF) is a computational imaging technique that enables potential endoscopic observations of biological samples at cellular scale. In this work, we show that this technique is tantamount to collecting multiple symmetric rank-one projections (SROP) of an interferometric matrix--a matrix encoding the spectral content of the sample image. In this model, each SROP is induced by the complex sketching vector shaping the incident light wavefront with a spatial light modulator (SLM), while the projected interferometric matrix collects up to $O(Q^2)$ image frequencies for a $Q$-core MCF. While this scheme subsumes previous sensing modalities, such as raster scanning (RS) imaging with beamformed illumination, we demonstrate that collecting the measurements of $M$ random SLM configurations--and thus acquiring $M$ SROPs--allows us to estimate an image of interest if $M$ and $Q$ scale log-linearly with the image sparsity level This demonstration is achieved both theoretically, with a specific restricted isometry analysis of the sensing scheme, and with extensive Monte Carlo experiments. On a practical side, we perform a single calibration of the sensing system robust to certain deviations to the theoretical model and independent of the sketching vectors used during the imaging phase. Experimental results made on an actual MCF system demonstrate the effectiveness of this imaging procedure on a benchmark image.

We inspect the propagation of shear polarized surface waves akin to Love waves through a forest of trees of same height atop a guiding layer on a soil substrate. We discover that the foliage of trees { brings a radical change in} the nature of the dispersion relation of these surface waves, which behave like spoof plasmons in the limit of a vanishing guiding layer, and like Love waves in the limit of trees with a vanishing height. When we consider a forest with trees of increasing or decreasing height, this hybrid "Spoof Love" wave is either reflected backwards or converted into a downward propagating bulk wave. An asymptotic analysis shows the forest behaves like an anisotropic wedge with effective boundary conditions.