Wavefront-modified vector beams for THz cornea spectroscopy.

Terahertz spectroscopy is a promising method to diagnose ocular diseases, where the cornea is typically imaged by Gaussian beams. However, the beam's mismatch with the cornea's spherical surface produces a 5-10 % error in analysis. We investigate cornea spectroscopy with wavefront-modified vector beams, reducing the original analysis error to less than 0.5 %. Vector beams are synthesized by our developed 3D Angular Spectrum Method expanded to vector spherical harmonic presentation, allowing wavefront modification and scattering analysis from 100-layer cornea models. We show that wavefront-modified spherical vector beams possess increased accuracy and non-sensitive focusing on cornea spectroscopy compared to the Gaussian beams. Additionally, we investigate wavefront-modified cylindrical vector beams, which show frequency-dependent scattering power arising from s- and p-polarizations. As a result, these beams are unsuitable for cornea spectroscopy, although they have potential for optical force applications. Wavefront-modified vector beams can be applied to spherical target spectroscopy and optical force applications, such as medicine, medical imaging, and optical tweezers.

Mie scattering with 3D angular spectrum method.

Mie theory is a powerful method to model electromagnetic scattering from a multilayered sphere. Usually, the incident beam is expanded to its vector spherical harmonic representation defined by beam shape coefficients, and the multilayer sphere scattering is obtained by the T-matrix method. However, obtaining the beam shape coefficients for arbitrarily shaped incident beams has limitations on source locations and requires different methods when the incident beam is defined inside or outside the computational domain or at the scatterer surface. We propose a 3D angular spectrum method for defining beam shape coefficients from arbitrary source field distributions. This method enables the placement of the sources freely within the computational domain without singularities, allowing flexibility in beam design. We demonstrate incident field synthesis and spherical scattering by comparing morphology-dependent resonances to known values, achieving excellent matching and high accuracy. The proposed method has significant benefits for optical systems and inverse beam design. It allows for the analysis of electromagnetic forward/backward propagation between optical elements and spherical targets using a single method. It is also valuable for optical force beam design and analysis.

Curved boundary integral method for electromagnetic fields.

The angular spectrum method is a rigorous method to synthesize near and far-field electromagnetic beams from planar field distributions. However, this limitation of planar surfaces has restricted its applicability to beams with simple focal planes. We propose a curved boundary integral method (CBIM) to synthesize electromagnetic beams from arbitrary surfaces to address this limitation and expand the method's scope to synthesize beams from and between shaped objects. This study presents a detailed theoretical framework behind the CBIM and validates its effectiveness and accuracy with a comprehensive set of simulations. Additionally, we present mathematical proof to support our proposal. The proposed method satisfies Maxwell's equations and significantly benefits optical systems and inverse beam design. It allows for analyzing electromagnetic forward/backward propagation between optical elements using a single method. It is also valuable for optical force beam design and analysis.

650 GHz imaging as alignment verification for millimeter wave corneal reflectometry.

A system concept for online alignment verification of millimeter-wave, corneal reflectometry is presented. The system utilizes beam scanning to generate magnitude-only reflectivity maps of the cornea at 650 GHz and compares these images to a precomputed/measured template map to confirm/reject sufficient alignment. A system utilizing 5 off-axis parabolic mirrors, a thin film beam splitter, and 2-axis galvanometric mirror was designed, simulated, and evaluated with geometric and physical optics. Simulation results informed the construction of a demonstrator system which was tested with a reference reflector. Similarity metrics computed with the aligned template and 26 misaligned positions, distributed on a 0.5 mm x 0.5 mm x 0.5 mm mesh, demonstrated sufficient misalignment detection sensitivity in 23 out of 26 positions. The results show that positional accuracy on the order of 0.5 mm is possible using 0.462 mm wavelength radiation due to the perturbation of coupling efficiency via beam distortion and beam walk-off.

Vector spherical harmonic analysis and experimental validation of spherical shells illuminated with broadband, millimeter wave Gaussian beams: applications to corneal sensing.

Coupling to longitudinal modes of thin spherical shells, under Gaussian-beam illumination, was explored with a theoretical method based on Fourier-optics analysis and vector spherical harmonics and was scrutinized with an experimental setup. For the theory part, the illumination frequency band was fixed between 100-600 GHz and the outer spherical shell radius of curvature and thickness are 7.5 mm and 0.5 mm, respectively. The shell material was either the lossless cornea or an aqueous effective media representing the cornea. Six different beam-target strategies were introduced being potential candidates for maximum coupling. Two dispersion-tuned beam ensembles with strongly frequency-dependent phase center location have been created with a fixed incident beam 1/e radius and radius of curvature called forward strategies. These computations of different alignments were continued with four beam ensembles of frequency-invariant phase center, constructed from fits to experimental data, oriented at four different axial locations with respect to the spherical shell center of curvature, they are called reverse strategies. Coupling efficiency for all strategies was calculated for different targets including perfect electrical conductor (PEC) sphere, PEC core covered by a cornea loss-free layer and cornea. All scattering strategies contrasted to scattering from equivalent planar targets as a reference with maximum coupling. The results show that, under an ideal calibration, forward strategies are a closer approximation to the plane-wave condition for the cornea. An experimental setup was assembled to explore the simulation approach in a frequency range between 220 GHz to 330 GHz. Two different quartz samples with permittivity of 4.1 were mounted on a water core, acting for a cornea. The first and second quartz radius and thickness were 7.5 mm and 0.5 mm and 8 mm and 1 mm, respectively. An adequate agreement between theory and experiment was confirmed. A particle optimisation swarm algorithm was applied to extract the thickness and permittivity of quartz from the measured back-scattered field for reverse strategies.

Calibration Alignment Sensitivity in Corneal Terahertz Imaging.

Improving the longitudinal modes coupling in layered spherical structure contributes significantly to corneal terahertz sensing, which plays a crucial role in the early diagnosis of cornea dystrophies. Using a steel sphere to calibrate reflection from the cornea sample assists in enhancing the resolution of longitudinal modes. The requirement and challenges toward applying the calibration sphere are introduced and addressed. Six corneas with different properties are spotted to study the effect of perturbations in the calibration sphere in a frequency range from 100 GHz to 600 GHz. A particle-swarm optimization algorithm is employed to quantify corneal characteristics considering cases of accurately calibrated and perturbed calibrated scenarios. For the first case, the study is carried out with signal-to-noise values of 40 dB, 50 dB and 60 dB at waveguide bands WR-5.1, WR-3.4, and WR-2.2. As expected, better estimation is achieved in high-SNR cases. Furthermore, the lower waveguide band is revealed as the most proper band for the assessment of corneal features. For perturbed cases, the analysis is continued for the noise level of 60 dB in the three waveguide bands. Consequently, the error in the estimation of corneal properties rises significantly (around 30%).

Extreme Asymmetry in Metasurfaces via Evanescent Fields Engineering: Angular-Asymmetric Absorption.

On the quest towards full control over wave propagation, the development of compact devices that allow asymmetric response is a challenge. In this Letter, we introduce a new paradigm for the engineering of asymmetry in planar structures, revealing and exploiting unilateral excitation of evanescent waves. We test the idea with the design and experimental characterization of a metasurface for angular-asymmetric absorption. The results show that the contrast ratio of absorption (the asymmetry level) can be arbitrarily engineered from zero to infinity for waves coming from two oppositely tilted angles. We demonstrate that the revealed asymmetry effects cannot be realized using conventional diffraction gratings, reflectarrays, and phase-gradient metasurfaces. This Letter opens up promising possibilities for wave manipulation via evanescent waves engineering with applications in one-side detection and sensing, angle-encoded steganography, flat nonlinear devices, and shaping the scattering patterns of various objects.

Passive terahertz camera for standoff security screening.

We describe the construction and performance of a passive, real-time terahertz camera based on a modular, 64-element linear array of cryogenic hotspot microbolometers. A reflective conical scanner sweeps out a 2 m x 4 m (vertical x horizontal) field of view (FOV) at a standoff range of 8 m. The focal plane array is cooled to 4 K in a closed cycle refrigerator, and the signals are detected on free-standing bridges of superconducting Nb or NbN at the feeds of broadband planar spiral antennas. The NETD of the focal-plane array, referred to the target plane and to a frame rate of 5 s(-1), is 1.25 K near the center of the array and 2 K overall.