Numerical implementation of three-dimensional vectorial complex ray model and application to rainbow scattering of spheroidal drops.

The rainbow patterns of oblate spheroidal drops have been observed in experiments nearly forty years ago [Nature312, 529 (1984)10.1038/312529a0]. However, the prediction for those complex patterns has been a challenge for conventional light scattering models. The vectorial complex ray model (VCRM) allows to account for the direction, the polarization, the phase, the amplitude and the wavefront curvature of waves and provides a powerful tool for the study of the light/electromagnetic wave interaction with a homogeneous object of any shape with smooth surface. In [Opt. Lett.46, 4585 (2021)10.1364/OL.434149], the authors have reported an important breakthrough of VCRM for the three-dimensional scattering (VCRM3D) and the simulated rainbow patterns of oblate drops. The present paper is devoted to the detailed description of the numerical implementation allowing the simulation of the 3D scattering field by a nonspherical particle. Its ability to predict both the fine and coarse intensity structures of the rainbows and the near-backward scattering patterns of spheroids is demonstrated. This work opens perspectives for exploring the 3D scattering characteristics of large objects with any smooth shape and developing relevant optical techniques for particle characterization.

Airy theory revisited with the method combining vectorial complex ray model and physical optics.

Airy published his theory in the 1830s to remedy the problem of infinite intensity in the rainbow angles of a spherical droplet predicted by geometrical optics. This theory has been studied by mathematicians and physicists since then from different points of view. In what concerns the scattering diagram around the rainbow angles, Airy theory has been improved by researchers in order to predict correctly the intensity of the supernumerary bows. However, it is known that the positions and the intensities of the supernumerary bows predicted by Airy theory differ from those of rigorous Debye theory with increasing order p and the scattering angles from the rainbow angle. In the present Letter, we will show that this discrepancy is caused by the approximations in Airy theory and can be revised by combining the vectorial complex ray model and physical optics (PO). The former permits us to calculate rigorously the amplitudes and phases of all rays and predicts precisely the scattering pattern except for the main bows. The combination with PO predicts very precisely all the supernumerary bows for both perpendicular and parallel polarization. This method can be applied directly to the light scattering of non-spherical particles with smooth surfaces.

Generalized rainbow patterns of oblate drops simulated by a ray model in three dimensions.

The scattering patterns near the primary rainbow of oblate drops are simulated by extending the vectorial complex ray model (VCRM) [Opt. Lett.36, 370 (2011)OPLEDP0146-959210.1364/OL.36.000370] to three-dimensional (3D) calculations. With the curvature of a wavefront as an intrinsic property of a ray, this advanced ray model permits, in principle, to predict the amplitudes and phases of all emergent rays with a rigorous algebraic formalism. This Letter reports a breakthrough of VCRM for 3D scattering with a line-by-line triangulation interpolation algorithm allowing to calculate the total complex amplitude of a scattered field. This makes possible to simulate not only the skeleton (geometrical rainbow angles, hyperbolic-umbilic caustics), but also the coarse (Airy bows, lattice) and fine (ripple fringes) structures of the generalized rainbow patterns (GRPs) of oblate drops. The simulated results are found qualitatively and quantitatively in good agreement with experimental scattering patterns for drops of different aspect ratios. The physical interpretation of the GRPs is also given. This work opens up prominent perspectives for simulating and understanding the 3D scattering of large particles of any shape with a smooth surface by VCRM.

Computation of radiation pressure force exerted on arbitrary shaped homogeneous particles by high-order Bessel vortex beams using MLFMA.

Due to special characteristics of nondiffraction and self reconstruction, the Bessel beams have attracted wide attention in optical trapping and appear to be a dramatic alternative to Gaussian beams. We present in this paper an efficient approach based on the surface integral equations (SIE) to compute the radiation pressure force (RPF) exerted on arbitrary shaped homogeneous particles by high-order Bessel vortex beam (HOBVB). The incident beam is described by vector expressions perfectly satisfy Maxwell's equations. The problem is formulated with the combined tangential formulation (CTF) and solved iteratively with the aid of the multilevel fast multipole algorithm (MLFMA). Then RPF is computed by vector flux of the Maxwell's stress tensor over a spherical surface tightly enclosing the particle and analytical expression for electromagnetic fields of incident beam in near region are used. The numerical predictions are compared with the results of the rigorous method for spherical particle to validate the accuracy of the approach. Some numerical results on relative large particles of complex shape, such as biconcave cell-like particles with different geometry parameters are given, showing powerful capability of our approach. These results are expected to provide useful insights into the RPF exerted on complex shaped particles by HOBVB.

Experimental validation of the vectorial complex ray model on the inter-caustics scattering of oblate droplets.

We report the first experimental validation of the Vectorial Complex Ray Model (VCRM) using the scattering patterns of large oblate droplets trapped in an acoustic field. The two principal radii and refractive index of the droplets are retrieved with a minimization method that involves VCRM predictions and experimental light scattering patterns. The latter are recorded in the droplet equatorial plane between the primary rainbow region and the associated hyperbolic-umbilic diffraction catastrophe. The results demonstrate that the VCRM can predict the fine and coarse stuctures of scattering patterns with good precision, opening up perspectives for the characterization of large non-spherical particles.

Computational study of radiation torque on arbitrary shaped particles with MLFMA.

The surface integral equation (SIE) method is used for the computational study of radiation torque on arbitrarily shaped homogeneous particles. The Multilevel Fast Multipole Algorithm (MLFMA) is employed to reduce memory requirements and improve the capability of SIE. The resultant matrix equations are solved iteratively to obtain equivalent electric and magnetic currents. Then, radiation torque is computed using the vector flux of the pseudotensor over a spherical surface tightly enclosing the particle. We use, therefore, the analytical electromagnetic field expression for incident waves in the near region, instead of the far-field approximation. This avoids the error which may be caused when describing the incident beam. The numerical results of three kinds of non-spherical particles are presented to illustrate the validity and capability of the developed method. It is shown that our method can be applied to predict, in the rigorous sense, the torque from a beam of any shape on a particle of complex configuration with a size parameter as large as 650. The radiation torques on large ellipsoids are exemplified to show the performance of the method and to study the influence that different aspect ratios have on the results. Then, the code is used for the calculation of radiation torque on objects of complex shape including a biconcave cell-like particle and a motor with a non-smooth surface.

Computation of radiation pressure force on arbitrary shaped homogenous particles by multilevel fast multipole algorithm.

A full-wave numerical method based on the surface integral equation for computing radiation pressure force (RPF) exerted by a shaped light beam on arbitrary shaped homogenous particles is presented. The multilevel fast multipole algorithm is employed to reduce memory requirement and to improve its capability. The resultant matrix equation is solved by using an iterative solver to obtain equivalent electric and magnetic currents. Then RPF is computed by vector flux of the Maxwell's stress tensor over a spherical surface tightly enclosing the particle. So the analytical expressions for electromagnetic fields of incident beam in near region are used. Some numerical results are performed to illustrate the validity and capability of the developed method. Good agreements between our method and the Lorenz-Mie theory for spherical and small spheroidal particle are found while our method has powerful capability for computing RPF of any shaped beam on a relatively large particle of complex shape. Tests for ellipsoidal and red blood cell-like particles illuminated by Gaussian beam have shown that the size of the particle can be as large as 50-100 wavelengths, respectively, for the relative refractive of 1.33 and 1.1.

Scattering of a Gaussian beam by an elliptical cylinder using the vectorial complex ray model.

The scattered waves of a shaped beam by an infinite cylinder in the far field are, stricto sensu, neither cylindrical nor spherical, so the asymptotic form of special functions involved in the theories based on the rigorous solution of Maxwell equations cannot be used to evaluate scattered intensities, even in the most simple case of Gaussian beam scattering by an infinite circular cylinder. Thus, although theories exist for the scattering of a shaped beam by infinite cylinders with circular and elliptical sections, the numerical calculations are limited to the near field. The vectorial complex ray model (VCRM) developed by Ren et al. describes waves by rays with a new property: the curvature of the wavefront. It is suitable to deal with the scattering of an arbitrarily shaped beam by a particle with a smooth surface of any form. In this paper, we apply this method to the scattering of an infinite elliptical cylinder illuminated by a Gaussian beam at normal incidence with an arbitrary position and orientation relative to the symmetric axis of the elliptical section of the cylinder. The method for calculating the curvature of an arbitrary surface is given and applied in the determination of the two curvature radii of the Gaussian beam wavefront at any point. Scattered intensities for different parameters of the beam and the particle as well as observation distance are presented to reveal the scattering properties and new phenomena observed in the beam scattering by an infinite elliptical cylinder.

Scattering from an elliptical cylinder by using the vectorial complex ray model.

The vectorial complex ray model (VCRM) is applied to the light scattering of an elliptical cylinder illuminated by a plane wave. In the VCRM, all waves are described by vectorial complex rays, and the scattering intensities are computed by the superposition of the complex amplitudes of the vectorial rays. The significant merit of this approach is that the wave properties are integrated in the ray model such that the divergence/convergence of the wave each time it encounters a dioptric surface is deduced by the wavefront curvature equation, and the phase shifts due to the focal lines are determined directly by the curvature of the wavefront. The approach is particularly suitable for a large cylinder with an elliptical cross section.

Vectorial complex ray model and application to two-dimensional scattering of plane wave by a spheroidal particle.

A vectorial complex ray model is introduced to describe the scattering of a smooth surface object of arbitrary shape. In this model, all waves are considered as vectorial complex rays of four parameters: amplitude, phase, direction of propagation, and polarization. The ray direction and the wave divergence/convergence after each interaction of the wave with a dioptric surface as well as the phase shifts of each ray are determined by the vector Snell law and the wavefront equation according to the curvatures of the surfaces. The total scattered field is the superposition of the complex amplitude of all orders of the rays emergent from the object. Thanks to the simple representation of the wave, this model is very suitable for the description of the interaction of an arbitrary wave with an object of smooth surface and complex shape. The application of the model to two-dimensional scattering of a plane wave by a spheroid particle is presented as a demonstration.

Near-critical-angle scattering for the characterization of clouds of bubbles: particular effects.

We report experimental investigations on the influence of various optical effects on the far-field scattering pattern produced by a cloud of optical bubbles near the critical scattering angle. Among the effects considered, there is the change of the relative refractive index of the bubbles (gas bubbles or some liquid-liquid droplets), the influence of intensity gradients induced by the laser beam intensity profile and by the spatial filtering of the collection optics, the coherent and multiple scattering effects occurring for densely packed bubbles, and the tilt angle of spheroidal optical bubbles. The results obtained herein are thought to be fundamental for the development of future works to model these effects and for the extension of the range of applicability of an inverse technique (referenced herein as the critical angle refractometry and sizing technique), which is used to determine the size distribution and composition of bubbly flows.

Debye series analysis of radiation pressure force exerted on a multilayered sphere.

On the basis of generalized Lorenz-Mie theory, the Debye series expansion (DSE) for radiation pressure forces (RPF) exerted on a multilayered sphere induced by focused beams is introduced. The DSE can isolate the contribution of each scattering process to RPF, and give a physical explanation of RPF. Typically, the RPF induced by a Gaussian beam is studied. The DSE is employed to the simulation of RPF corresponding to different scattering processes (diffraction, reflection, refraction, etc.) in detail, and gives the physical mechanism of RPF. The effects of various parameters, such as scattering mode p, beam position, and radius of core for coated spheres, to RPF is researched.

Generalized Debye series expansion of electromagnetic plane wave scattering by an infinite multilayered cylinder at oblique incidence.

The Debye series expansion expresses the Mie scattering coefficients into a series of Fresnel coefficients and gives physical interpretation of different scattering modes, but when an infinite multilayered cylinder is obliquely illuminated by electromagnetic plane waves, the scattering process becomes very complicated because of cross polarization. Based on the relation of boundary conditions between global scattering process and local scattering processes, the generalized Debye series expansion of plane wave scattering by an infinite multilayered cylinder at oblique incidence is derived in this paper. The formula and the code are verified by the comparison of the results with that of Lorenz-Mie theory in special cases and those presented in the literatures.

Computation of stress on the surface of a soft homogeneous arbitrarily shaped particle.

Prediction of the stress on the surface of an arbitrarily shaped particle of soft material is essential in the study of elastic properties of the particles with optical force. It is also necessary in the manipulation and sorting of small particles with optical tweezers, since a regular-shaped particle, such as a sphere, may be deformed under the nonuniform optical stress on its surface. The stress profile on a spherical or small spheroidal soft particle trapped by shaped beams has been studied, however little work on computing the surface stress of an irregular-shaped particle has been reported. We apply in this paper the surface integral equation with multilevel fast multipole algorithm to compute the surface stress on soft homogeneous arbitrarily shaped particles. The comparison of the computed stress profile with that predicted by the generalized Lorenz-Mie theory for a water droplet of diameter equal to 51 wavelengths in a focused Gaussian beam show that the precision of our method is very good. Then stress profiles on spheroids with different aspect ratios are computed. The particles are illuminated by a Gaussian beam of different waist radius at different incidences. Physical analysis on the mechanism of optical stress is given with help of our recently developed vectorial complex ray model. It is found that the maximum of the stress profile on the surface of prolate spheroids is not only determined by the reflected and refracted rays (orders p=0,1) but also the rays undergoing one or two internal reflections where they focus. Computational study of stress on surface of a biconcave cell-like particle, which is a typical application in life science, is also undertaken.

Sizing of irregular particles using a near backscattered laser Doppler system.

A near backscattered laser Doppler system was presented to carry out velocity and size distribution measurements for irregular particles in two-phase flows. The technique uses amplitudes of particles Doppler signals to estimate the particle size distribution in a statistical manner. Holve's numerical inversion scheme is employed to unfold the dependence of the scattered signals on both particle trajectory and orientation through the measurement volume. The performance and error level of the technique were simulated, and several parameters including the number of particle samples, the fluctuation of irregular particle response function, inversion algorithms, and types of particle size distribution were extensively investigated. The results show that the size distributions for those irregular particles even with strong fluctuations in response function can be successfully reconstructed with an acceptable error level using a Phillips-Twomey-non-negative least-squares algorithm instead of a non-negative least-squares one. The measurement system was then further experimentally verified with irregular quartz sands. Using inversion matrix obtained from the calibration experiment, the average measurement error for the mixing quartz sands with a size range of 200-560 microm are found to be about 23.3%, which shows the reliability of the technique and the potential for it to be applied to industrial measurement.

Debye series for Gaussian beam scattering by a multilayered sphere.

The Debye series has been a key tool for the understanding of light scattering features, and it is also a convenient method for understanding and improving the design of optical instruments aimed at optical particle sizing. Gouesbet has derived the Debye series formulation for generalized Lorenz-Mie theory (GLMT). However, the scattering object is a homogeneous sphere, and no numerical result is provided. The Debye series formula for plane-wave scattering by multilayered spheres has been derived before. We have devoted our work to the Debye series of Gaussian beam scattering by multilayered spheres. The integral localized approximation is employed in the calculation of beam-shape coefficients (BSCs) and allows the study of the scattering characteristics of particles illuminated by the strongly focused beams. The formula and code are verified by the comparison with the results produced by GLMT and also by the comparison with the result for the case of plane-wave incidence. The formula is also employed in the simulation of the first rainbow by illuminating the particle with one or several narrow beams.

Extension of geometrical-optics approximation to on-axis Gaussian beam scattering. I. By a spherical particle.

The geometrical-optics approximation of light scattering by a transparent or absorbing spherical particle is extended from plane wave to Gaussian beam incidence. The formulas for the calculation of the phase of each ray and the divergence factor are revised, and the interference of all the emerging rays is taken into account. The extended geometrical-optics approximation (EGOA) permits one to calculate the scattering diagram in all directions from 0 degrees to 180 degrees. The intensities of the scattered field calculated by the EGOA are compared with those calculated by the generalized Lorenz-Mie theory, and good agreement is found. The surface wave effect in Gaussian beam scattering is also qualitatively analyzed by introducing a flux ratio factor. The approach proposed is particularly important to the further extension of the geometrical-optics approximation to the scattering of large spheroidal particles.

Debye series of normally incident plane-wave scattering by an infinite multilayered cylinder.

We derive the formula of the Debye-series decomposition for normally incident plane-wave scattering by an infinite multilayered cylinder. A comparison of the scattering diagrams calculated by the Debye series and Mie theory for a graded-index polymer optical fiber is given and the agreement is found to be satisfied. This approach permits us to simulate the rainbow intensity distribution of any single order and the interference of several orders, which is of good use to the study of the scattering characteristics of an inhomogeneous cylinder and to the measurement of the refractive index profile of an inhomogeneous cylinder.

Extension of geometrical-optics approximation to on-axis Gaussian beam scattering. II. By a spheroidal particle with end-on incidence.

On the basis of our previous work on the extension of the geometrical-optics approximation to Gaussian beam scattering by a spherical particle, we present a further extension of the method to the scattering of a transparent or absorbing spheroidal particle with the same symmetric axis as the incident beam. As was done for the spherical particle, the phase shifts of the emerging rays due to focal lines, optical path, and total reflection are carefully considered. The angular position of the geometric rainbow of primary order is theoretically predicted. Compared with our results, the Möbius prediction of the rainbow angle has a discrepancy of less than 0.5 degrees for a spheroidal droplet of aspect radio kappa within 0.95 and 1.05 and less than 2 degrees for kappa within 0.89 and 1.11. The flux ratio index F, which qualitatively indicates the effect of a surface wave, is also studied and found to be dependent on the size, refractive index, and surface curvature of the particle.

Debye series for light scattering by a multilayered sphere.

We have derived the formula for the Debye-series decomposition for light scattering by a multilayered sphere. This formulism permits the mechanism of light scattering to be studied. An efficient algorithm is introduced that permits stable calculation for a large sphere with many layers. The formation of triple first-order rainbows by a three-layered sphere and single-order rainbows and the interference of different-order rainbows by a sphere with a gradient refractive index, are then studied by use of the Debye model and Mie calculation. The possibility of taking only one single mode or several modes for each layer is shown to be useful in the study of the scattering characteristics of a multilayered sphere and in the measurement of the sizes and refractive indices of particles.