Iterative approach to Maxwell equations for dielectric media of spatially varying refractive index

Polarization-dependent tunneling of light in gradient optics

Reflection-refraction properties of photonic barriers, formed by dielectric gradient nanofilms, for inclined incidence of both S - and P -polarized electromagnetic waves are examined by means of exactly solvable models. We present generalized Fresnel formulas, describing the influence of the nonlocal dispersion on the reflectance and transmittance of single- and double-layer gradient photonic barriers for S and P waves and arbitrary angles of incidence. The nonlocal dispersion of such layers, arising due to a concave spatial profile of dielectric susceptibility across the plane film, is shown to result in a peculiar heterogeneity-induced optical anisotropy, providing the propagation of S (P) waves in tunneling (traveling) regimes. The results obtained indicate the possibility of narrow-band nonattenuated tunneling (complete transmittance) of oblique S waves through such heterogeneous barriers, and the existence of spectral areas characterized by the strong reflection of P waves and profound contrast between transmitted S and P waves. The scalability of obtained exact analytical solutions of Maxwell equations into the different spectral ranges is discussed and the application potential of these phenomena for miniaturized polarizers and filters is demonstrated.

Light scattering regimes along the optical axis in turbid media

We inject an angularly collimated laser beam into a scattering medium of a nondairy creamer-water solution and examine the distribution of the scattered light along the optical axis as a function of the source-detector spacing. The experimental and simulated data obtained from a Monte Carlo simulation suggest four regimes characterizing the transition from unscattered to diffusive light. We compare the data also with theoretical predictions based on a first-order scattering theory for regions close to the source, and with diffusionlike theories for larger source-detector spacings. We demonstrate the impact of the measurement process and the effect of the unavoidable absorption of photons by the detection fiber on the light distribution inside the medium. We show that the range of validity of these theories can depend on the experimental parameters such as the diameter and acceptance angle of the detection fiber.

Nonperturbative retrieval of the scattering strength in one-dimensional media

We examine several approaches on how to use the transmission and reflection amplitudes as functions of the modulation frequency of the laser's intensity to reconstruct the position-dependent scattering coefficient for a simple turbid medium. We explore the region where the contrast between the coefficient and its spatially averaged value is large enough such that perturbative methods fail. We show that in the case of a transillumination geometry, the knowledge of the transmission profile alone is not sufficient for unique image reconstruction, whereas the reflection spectrum allows for a complete inversion. We demonstrate the invertibility for media sampled at only a few positions.

Fast iterative, coupled-integral-equation technique for inhomogeneous profiled and periodic slabs

A fast coupled-integral-equation (CIE) technique is developed to compute the plane-TE-wave scattering by a wide class of periodic 2D inhomogeneous structures with curvilinear boundaries, which includes finite-thickness relief and rod gratings made of homogeneous material as special cases. The CIEs in the spectral domain are derived from the standard volume electric field integral equation. The kernel of the CIEs is of Picard type and offers therefore the possibility of deriving recursions, which allow the computation of the convolution integrals occurring in the CIEs with linear amounts of arithmetic complexity and memory. To utilize this advantage, the CIEs are solved iteratively. We apply the biconjugate gradient stabilized method. To make the iterative solution process faster, an efficient preconditioning operator (PO) is proposed that is based on a formal analytical inversion of the CIEs. The application of the PO also takes only linear complexity and memory. Numerical studies are carried out to demonstrate the potential and flexibility of the CIE technique proposed. Though the best efficiency and accuracy are observed at either low permittivity contrast or high conductivity, the technique can be used in a wide range of variation of material parameters of the structures including when they contain components made of both dielectrics with high permittivity and typical metals.