Coherent backscattering effects for discrete random media: Numerical and theoretical results

Manifestation of the backscattering enhancement phenomenon in the reflection matrix elements of the coherent component of scattered radiation is considered. The dependence of the coherent backscattering effects on the microphysical properties of the medium scatterers are investigated. It is shown that random media of fractal-like clusters exhibit brightness and polarization opposition effects, which are like those observed for some atmosphereless Solar system bodies. Conditions for a bimodal angle dependence in the degree of linear polarization are discussed and the manifestation of the enhanced backscattering phenomenon in the intensity of scattered radiation is studied.     

Interference effects in backscattering of light by a layer of a discrete random medium

Multiple backscattering of light by a layer of a discrete random medium is considered. A brief derivation of equations for describing the coherent and incoherent components of scattered light is presented. These equations are solved numerically in the approximation of doubled scattering of light by a semi-infinite medium of spherical scatterers having a size comparable with the wavelength in order to study the effect of the properties of particles on the angular dependence of interference effects. Calculations show that the half-width of the interference peak decreases upon an increase in lateral scattering by particles and that the degree of polarization has a complex angular dependence on the properties of the particles. For an optically thin layer of the medium, the relations defining the interference peak half-width and the scattering angle upon extreme linear polarization as functions of the effective refractive index are given.     

Similarity and diversity in photometric and polarimetric opposition effects of small Solar System bodies

The phase-angle dependences of brightness and polarization of light scattered by atmosphereless Solar System bodies (satellites, asteroids, planetary rings) as well as comets and zodiacal light are analysed on the basis of available ground-based and spacecraft observations. We study similarity and diversity in the photometric and polarimetric opposition effects. The similarity of the brightness and polarization phase functions for polydisperse dust media (cometary and interplanetary dust) and for atmosphereless bodies at small phase angles gives grounds to state that at least some physical properties of dust particles and light-scattering mechanisms should be similar for the two classes of objects. The aggregate structure of particles in different objects can be precisely the property that determines the photometric and polarimetric effects observed in the opposition region. Differences observed in opposition effects for different objects are likely caused by different physical properties (composition, sizes, structure, density) of the scattering particles and, therefore, by different relative contributions of the light-scattering mechanisms. However, the relationship between the photometric opposition effect, spike effect, negative polarization branch, and polarization opposition effect for different bodies is ambiguous.     

Anomalous polarization effects during light scattering in random media

The dependence of the intensity of light backscattered from a layer of a randomly inhomogeneous medium on the polarization of incident light and the size of scatterers has been investigated. The results of numerical simulation have demonstrated that the direction of rotation of the plane of polarization is different in systems with small- and large-scale inhomogeneities. It is shown for the first time that the dependence of the sign of the residual circular polarization on the size of scatterers can be observed in systems described by the Henyey-Greenstein phase function used in simulating biological tissues. A similar anomalous polarization effect, which consists in changing the direction of rotation of the plane of polarization of backscattered light with an increase in the scattering angle, is revealed in studying the coherent backscattering component. These polarization effects are observed in light backscattering from optically active media.     

Polarization-based range-gated imaging in birefringent medium:Effect of size parameter

We have investigated the effect of size parameter of the scatterer on the image quality obtained with polarization-based range-gated imaging in birefringent turbid medium. Both linearly and circularly polarized light were utilized for imaging.The simulated results indicate that the improvement of visibility is more pronounced using circularly polarized light for the birefringent medium composed of smaller-sized scatterers at lower values of optical thickness and the birefringent medium comprising larger-sized scatterers. In contrast, linearly polarized light provides better image quality for the birefringent medium composed of smaller-sized scatterers at larger values of optical thickness. The evolution of the polarization characteristics of backscattered light and target light under the conditions mentioned above was measured to account for these numerical results.     

On optical properties of nonspherical inhomogeneous particles

To solve the problem of light scattering by multilayer scatterers of an arbitrary axisymmetric shape, a separation of variables method that involves special scalar potentials and their expansions in spherical functions is developed. The approach is shown to yield highly exact results even for particles that have 100 layers or more. A graphic library that illustrates the optical properties of layered and homogeneous (with an effective refractive index) spheroids, spheres, and Chebyshev particles of various shapes and sizes (about 650 figures) is created and is put on the Internet. It is noted that the linear polarization of radiation transmitted forward through a polydisperse medium containing partially oriented nonspherical porous particles strongly depends on the structure of scatterers. It is shown that the difference between the degrees of polarization of layered and corresponding homogeneous scatterers can exceed 200–300%.     

Light refraction: Right or left turn of rays as polarization effect

Media with negative dispersion are represented as a set of specialized scatterers that will change, at waves' absorption and their reemission, the sign of magnetic strength parallel to surface. In metamaterials such scatterers (turn elements) are represented by the open oscillatory contours. In natural substance (metals) their role can execute magnetic dipoles and electrical quadropoles induced by polarization. Their converting features become decisive in sufficiently thin films at suppressing dipole polarization, i.e., such phenomena are shown the dimensional (size) effects of media polarization. This approach leads to estimations of optical thickness of medium, in which a turn of refracted beam is observable.     

Dependence of the circular polarization of backscattered light in random media on anisotropy of scatterers

The scattering of linearly or circularly polarized light from a semi-infinite randomly inhomogeneous medium is considered. Using the Monte Carlo method, it is shown that, in the case of a wide front of incident and scattered optical radiation and irrespective of the degree of scattering anisotropy, the copolarized component of backscattered light dominates the cross-polarized component for the linear polarization and the cross-polarized component dominates the copolarized component for the circular polarization. If the beams of incident and scattered radiation are spatially separated and the size of scatterers exceeds the wavelength, the circular copolarized component dominates the cross-polarized one. A similar effect of the change in direction of the rotation of the plane of polarization in relation to the size of scatterers is revealed for pulsed radiation.     

Experimental study of the properties of non-linear scatterers using a set of polarization patterns

We examine the possibility of experimental study of the polarization and polarization-amplitude properties of nonlinear scatterers using sets of polarization patterns. We give examples of the sets of polarization patterns measured at the frequency of the second harmonic of the sounding signal. We point out the significant difference between the polarization and polarization-amplitude properties of elementary and complex nonlinear scatterers. An elementary nonlinear scatterer is shown to be a polarization-selective object whose properties can be described completely on the basis of measurements. Radiophysical Research Institute, Nizhny Novgorod, Russia. Translated from izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 42, No. 2, pp. 129–138, February 1999.     

随机非球形粒子全极化散射的时间相关Mueller矩阵解

从与时间相关的矢量辐射传输方程推导一阶Mueller矩阵解,用来模拟Gauss型平面脉冲波入射下,一层随机、非均匀取向非球形粒子的全极化双站散射.数值计算了同极化和去极化脉冲响应,与入射脉冲进行了比较,说明了随机介质的物理参数,如粒子的取向和占空比、入射角、极化以及层厚等对脉冲响应的影响 关键词： 平面脉冲波 非球形粒子 Mueller矩阵     

Light scattering by aggregates of varying porosity and the opposition phenomena observed in the low-albedo particulate media

The rule that the opposition phenomena in brightness and linear polarization observed in many regolith surfaces usually accompany each other is violated in the cases of very dark asteroids and particulate samples: practically no nonlinear surge of brightness to opposition is observed while the branch of negative polarization at small phase angles exists. To explain this fact, we model the light scattering by particulate media with ensembles of spherical particles (with size comparable to the wavelength) of varying packing density and refractive index. The superposition T-matrix method is used. The increase in the absorption and/or packing density diminishes the amplitude of the brightness opposition peak, and its profile becomes wider. The influence on the branch of negative polarization is more complex and depends on the relation between the size parameters of the constituents, the refractive index, and the porosity. However, the feature common to all considered cases is that the negative branch changes its shape and the polarization minimum moves to the inversion point. This behavior radically differs from that observed in nonabsorbing ensembles of particles and reflects the fact that the efficiency of the coherent backscattering, which mainly determines these characteristics in nonabsorbing ensembles (to the packing density of about 30%), decreases. Moreover, since the angular profiles are not simply damped, but the polarization minimum changes its angular position, we may conclude that the near-field interaction of the constituents becomes important: the shielding of particles by each other eliminates many constituents from the scattering and the near-field effects promote the negative polarization and smooth the backscattering brightness surge. Due to this, when the packing density exceeds 10-20%, the opposition phenomena in absorbing ensembles are caused not only by the coherent backscattering, and situations, when the opposition brightness surge is practically suppressed, but the negative branch of polarization still survives, are possible. This may explain the fact that the dark regolith surfaces show no brightness opposition effect, but produce the branch of negative polarization with the minimum shifted from opposition.     

Measurements of structure-induced polarization features in forward scattering from collections of cylindrical fibers

Many applications in remote sensing, material sciences and biomedical field are characterized by a transition domain between single scattering and multiple-scattering regimes. This regime is described by typical polarization features which can be used to retrieve structural information. An electronically agile technique was used for measuring in real time the Stokes vectors of light incident on and emerging from an inhomogeneous medium. Subsequently, the Mueller matrix associated with the scattering medium is determined. We focus our attention on forward scattering from systems consisting of random as well as partially oriented long cylindrical fibers. We discuss the effects of: (1) shape of individual scattering centers, (2) structure parameter, and (3) optical density of the scattering medium. The anisotropic behavior of the structure function at different packing fractions determines nontrivial characteristics of the polarization transfer. The complex effective index of refraction can be polarization dependent as a result of the optical anisotropy due to both the shape of the individual scatterers and the structure characteristics of the scattering system. Some of the Mueller matrix elements are shown to be related to the optical anisotropy of the system for the case of long cylindrical fibers. The polarization efficiency, the structure parameter, and the packing fraction are used to quantify this relationship. We also found that some of the matrix elements are more sensitive to the degree of structural anisotropy and the packing fraction, while other elements are sensitive to structural non-uniformities across the investigated area.     

Radiative transfer theory with time delay for the effect of pulse entrapping in a resonant random medium: General transfer equation and point-like scatterer model

Abstract

A pulse propagation of a vector electromagnetic wave field in a discrete random medium under the condition of Mie resonant scattering is considered on the basis of the Bethe–Salpeter equation in the two-frequency domain in the form of an exact kinetic equation which takes into account the energy accumulation inside scatterers. The kinetic equation is simplified using the transverse field and far wave zone approximations which give a new general tensor radiative transfer equation with strong time delay by resonant scattering. This new general radiative transfer equation, being specified in terms of the low-density limit and the resonant point-like scatterer model, takes the form of a new tensor radiative transfer equation with three Lorentzian time-delay kernels by resonant scattering. In contrast to the known phenomenological scalar Sobolev equation with one Lorentzian time-delay kernel, the derived radiative transfer equation does take into account effects of (i) the radiation polarization, (ii) the energy accumulation inside scatterers, (iii) the time delay in three terms, namely in terms with the Rayleigh phase tensor, the extinction coefficient and a coefficient of the energy accumulation inside scatterers, respectively (i.e. not only in a term with the Rayleigh phase tensor). It is worth noting that the derived radiative transfer equation is coordinated with Poynting's theorem for non-stationary radiation, unlike the Sobolev equation. The derived radiative transfer equation is applied to study the Compton–Milne effect of a pulse entrapping by its diffuse reflection from the semi-infinite random medium when the pulse, while propagating in the medium, spends most of its time inside scatterers. This specific albedo problem for the derived radiative transfer equation is resolved in scalar approximation using a version of the time-dependent invariance principle. In fact, the scattering function of the diffusely reflected pulse is expressed in terms of a generalized time-dependent Chandrasekhar H-function which satisfies a governing nonlinear integral equation. Simple analytic asymptotics are obtained for the scattering function of the front and the back parts of the diffusely reflected Dirac delta function incident pulse, depending on time, the angle of reflection, the mean free time, the microscopic time delay and a parameter of the energy accumulation inside scatterers. These asymptotics show quantitatively how the rate of increase of the front part and the rate of decrease of the rear part of the diffusely reflected pulse become slower with transition from the regime of conventional radiative transfer to that of pulse entrapping in the resonant random medium.     

Energy exchange between laser pulses in an atomic medium with a closed excitation contour

The effect the phase difference between carrier frequencies of two laser pulses has on their propagation in an optically dense medium of three-level point scatterers is investigated. The medium is exposed to constant radiation, creating a closed excitation contour. Situations are considered for pulse durations shorter and longer than the lifetime of the excited state of scatterers.     

Quantized media with absorptive scatterers and modified atomic emission rates

Modifications in the spontaneous emission rate of an excited atom that are caused by extinction effects in a nearby dielectric medium are analyzed in a quantummechanical model, in which the medium consists of spherical scatterers with absorptive properties. Use of the dyadic Green function of the electromagnetic field near a dielectric sphere leads to an expression for the change in the emission rate as a series of multipole contributions for which analytical formulas are obtained. The results for the modified emission rate as a function of the distance between the excited atom and the dielectric medium show the influence of both absorption and scattering processes.     

Speckle polarization diagnostics of scattering media using partially coherent radiation

The paper discusses a method for probing a randomly inhomogeneous medium using partially coherent radiation and polarization filtering of the radiation scattered by the medium. The method is based on the analysis of the contrast of speckle-modulated images of the object under study as a function of the coherence length of the probe radiation. A theoretical justification of the method is given as applied to systems of discrete scatterers with the subsequent modification of the results obtained for continuously distributed scattering systems. The results of experimental testing of the developed method for scattering media characterized by nondiffusion conditions of propagation of the probe radiation are presented and compared with the results of the theoretical analysis.     

Aerosol influence on polarization and intensity in near-infrared O2 and CO2 absorption bands observed from space

We study the intensity and degree of linear polarization of reflected solar radiation at the top of the atmosphere within two carbon dioxide bands and one oxygen absorption band in the near-infrared. In particular, we are interested in the sensitivity of the degree of linear polarization and intensity to changes of aerosol and cirrus cloud layer heights, microphysical properties, and surface albedo. For the simulations we use spectral response functions representative of the Orbiting Carbon Observatory (OCO). Inside the O₂A band at 760 nm and strong CO₂ band at 2060 nm we find a strong influence of the aerosol and cirrus cloud layer height on the degree of linear polarization. An increase of the aerosol or cirrus cloud layer height can lead either to a decrease or increase of the polarization within the band, depending on the microphysical and optical properties of the scatterers, surface albedo, and absorption strength in the bands. The results for the O₂A band also indicate that even over land OCO enables an estimation of the height of an aerosol or cirrus cloud layer. Inside the weak CO₂ band at 1610 nm the influence of aerosol or cirrus cloud layer heights is lower as compared to the O₂A band and CO₂ band at 2060 nm, due to the relatively stronger surface influence. Here an increase of aerosol or cirrus cloud layer height leads to an increase of the degree of linear polarization even in case of low surface albedo and for weakly polarizing scatterers. For the weak CO₂ band at 1610 nm we also study the influence of the aerosol or cirrus cloud layer height on the column CO₂ estimate and the errors resulting from ignoring polarization in simulations of backscatter measurements by space-based instruments such as OCO. Depending on the surface albedo, misinterpretations of the height of atmospheric scatterers might strongly affect the column CO₂ estimates.     

Magnetoelectric-Field Electrodynamics: Search for Magnetoelectric Point Scatterers

Resonant scattering of electromagnetic (EM) waves by small particles is considered as one of the basic problems in metamaterial science. At present, special subwavelength resonators are considered as structural elements in chiral and bianisotropic metamaterials. There is a general consensus that these small scatterers behave like “artificial atoms” with strong electrical and magnetic responses and an interconnection between these responses. However, the observed effect of magnetoelectric (ME) coupling in these meta-atoms is not associated with the near-field manipulation properties caused by intrinsic magnetoelectricity. This arises the question whether ME point scatterers of EM radiation really exist. In this paper, we show that there are mesoscopic structures with electric and magnetic dipole-carrying excitations that behave like point scatterers with their inherent magnetoelectricity. In such subwavelength resonators, coherent oscillations of the electric polarization and magnetization can be considered as quasistatic oscillations described by electrostatic (ES) and magnetostatic (MS) scalar wave functions. The ME resonance effect arises from the coupling of two, ES and MS, oscillations. The near fields of these resonators, called the ME near fields, are characterized by simultaneous violation of time reversal and inversion symmetry. In study of ME fields and EM problems associated with these fields, we put forward the concept of ME-field electrodynamics.     

Coherent backscattering of light by complex random media of spherical scatterers: numerical solution

Novel Monte Carlo techniques are described for the computation of reflection coefficient matrices for multiple scattering of light in plane-parallel random media of spherical scatterers. The present multiple scattering theory is composed of coherent backscattering and radiative transfer. In the radiative transfer part, the Stokes parameters of light escaping from the medium are updated at each scattering process in predefined angles of emergence. The scattering directions at each process are randomized using probability densities for the polar and azimuthal scattering angles: the former angle is generated using the single-scattering phase function, whereafter the latter follows from Kepler's equation. For spherical scatterers in the Rayleigh regime, randomization proceeds semi-analytically whereas, beyond that regime, cubic spline presentation of the scattering matrix is used for numerical computations. In the coherent backscattering part, the reciprocity of electromagnetic waves in the backscattering direction allows the renormalization of the reversely propagating waves, whereafter the scattering characteristics are computed in other directions. High orders of scattering (~10 000) can be treated because of the peculiar polarization characteristics of the reverse wave: after a number of scatterings, the polarization state of the reverse wave becomes independent of that of the incident wave, that is, it becomes fully dictated by the scatterings at the end of the reverse path. The coherent backscattering part depends on the single-scattering albedo in a non-monotonous way, the most pronounced signatures showing up for absorbing scatterers. The numerical results compare favourably to the literature results for nonabsorbing spherical scatterers both in and beyond the Rayleigh regime.     

Spatio-temporal invariants of the time reversal operator

The decomposition of the time reversal operator, known by the French acronym DORT, is a technique to extract point scatterers' monochromatic Green's functions from a medium. It is used to detect, locate, and focus on scatterers in various domains such as underwater acoustics, medical ultrasound, and nondestructive evaluation. A limitation of the method arises from its single-frequency nature, when the signals used in acoustics are often broadband. Reconstruction of the broadband Green's functions from the single-frequency Green's functions can be very difficult when numerous scatterers are present in the medium. Moreover, the method does not take advantage of the axial resolution associated with broadband signals. Time domain methods are investigated here as an answer to these problems. It is shown that the time reversal operator in the time domain takes the form of a tensor. The properties of the invariants are discussed. It is shown they do not have all the expected properties. Another method is proposed that requires a priori information on the medium.