Non-spherical particles in optical tweezers: A numerical solution.

We present numerical methods for modeling the dynamics of arbitrarily shaped particles trapped within optical tweezers, which improve the predictive power of numerical simulations for practical use. We study the dependence of trapping on the shape and size of particles in a single continuous wave beam setup. We also consider the implications of different particle compositions, beam types and media. The major result of the study is that for different irregular particle shapes, a range of beam powers generally leads to trapping. The trapping power range depends on whether the particle can be characterized as elongated or flattened, and the range is also limited by Brownian forces.

Scattering And Absorption of Light in Planetary Regoliths.

Theoretical, numerical, and experimental methods are presented for multiple scattering of light in macroscopic discrete random media of densely-packed microscopic particles. The theoretical and numerical methods constitute a framework of Radiative Transfer with Reciprocal Transactions (R2T2). The R2T2 framework entails Monte Carlo order-of-scattering tracing of interactions in the frequency space, assuming that the fundamental scatterers and absorbers are wavelength-scale volume elements composed of large numbers of randomly distributed particles. The discrete random media are fully packed with the volume elements. For spherical and nonspherical particles, the interactions within the volume elements are computed exactly using the Superposition T-Matrix Method (STMM) and the Volume Integral Equation Method (VIEM), respectively. For both particle types, the interactions between different volume elements are computed exactly using the STMM. As the tracing takes place within the discrete random media, incoherent electromagnetic fields are utilized, that is, the coherent field of the volume elements is removed from the interactions. The experimental methods are based on acoustic levitation of the samples for non-contact, non-destructive scattering measurements. The levitation entails full ultrasonic control of the sample position and orientation, that is, six degrees of freedom. The light source is a laser-driven white-light source with a monochromator and polarizer. The detector is a mini-photomultiplier tube on a rotating wheel, equipped with polarizers. The R2T2 is validated using measurements for a mm-scale spherical sample of densely-packed spherical silica particles. After validation, the methods are applied to interpret astronomical observations for asteroid (4) Vesta and comet 67P/Churyumov-Gerasimenko (Figure 1) recently visited by the NASA Dawn mission and the ESA Rosetta mission, respectively.

Radiative transfer with reciprocal transactions: Numerical method and its implementation.

We present a numerical method for solving electromagnetic scattering by dense discrete random media entitled radiative transfer with reciprocal transactions (R2T2). The R2T2 is a combination of the Monte Carlo radiative-transfer, coherent-backscattering, and superposition T-matrix methods. The applicability of the radiative transfer is extended to dense random media by incorporating incoherent volume elements containing multiple particles. We analyze the R2T2 by comparing the results with the asymptotically exact superposition T-matrix method, and show that the R2T2 removes the caveats of radiative-transfer methods by comparing it to the RT-CB. We study various implementation choices that result in an accurate and efficient numerical algorithm. In particular, we focus on the properties of the incoherent volume elements and their effects on the final solution.

Multiple scattering of light in discrete random media using incoherent interactions.

We consider the scattering and absorption of light in discrete random media of densely packed spherical particles. In what we term "radiative transfer with reciprocal transactions" (R2T2), we introduce a volume element of the random medium, derive its scattering and absorption characteristics using the superposition T-Matrix method (STMM), and compute its frequency-domain incoherent volume-element scattering characteristics. Using an order-of-scattering approach, we then compute a numerical Monte Carlo solution for the scattering problem with an exact treatment of the interaction between two volume elements. We compute both the direct and reciprocal contributions along a sequence of volume elements, allowing us to evaluate the coherent backscattering effects. We show that the R2T2 and exact STMM solutions are in mutual agreement for large finite systems of densely packed spherical particles. We conclude that the R2T2 method provides a viable numerical solution for scattering by asymptotically infinite systems of particles.

Polarized backscattering by clusters of spherical particles.

The linear and circular polarization ratios for clusters of spherical particles averaged over multiple orientations show a systematic pattern as a function of the refractive index and the size parameter. We show that, at backscattering, the depolarizing behavior of orientation-averaged clusters of spheres can be approximated by second-order scattering of bispheres. The pattern is relatively invariable in terms of the number of particles. We also demonstrate the significance of the near-field effects for polarization at backscattering.

Interpretation of single-particle negative polarization at intermediate scattering angles.

We study the interrelation of the internal field of irregular particles to the far-field scattering characteristics by modifying the internal field of dipole groups. In this paper, we concentrate on the longitudinal component, i.e., the internal-field component parallel to the incident wave vector. We use the discrete-dipole approximation to determine the internal field and switch off the longitudinal component from the dipoles that have the highest energy density above a preset cutoff value. We conclude that only a relatively small number of core dipoles, about 5% of all dipoles, contribute to the negative linear polarization at intermediate scattering angles. These core dipole groups are located at the forward part of the particles. The number of core dipoles in the group becomes greater as particle asphericity increases. We find that the interference between the scattered waves from the core dipole groups, which was studied previously for spherical particles, is preserved to a large extent for nonspherical particles.

Validity criteria of the discrete dipole approximation.

There are two widely accepted restrictions on the application of the discrete dipole approximation (DDA) in the study of light scattering by particles comparable to the wavelength: (1) when considering dielectric particles, the size of the cells must satisfy the condition kd|m|<0.5, where k is the wavenumber, d is the size of the cells, and m is the complex refractive index of the constituent material and (2) when considering conductive particles, the size of the cells must be small enough to reproduce sufficiently the evolution of the electromagnetic field in the skin layer. We examine both restrictions when the DDA is applied to irregularly shaped particles and show that its restrictions are not as strong as is widely accepted. For instance, when studying irregularly shaped particles averaged over orientations, even at kd|m|=1, the DDA provides highly accurate numerical results. Moreover, we show that the impact of using large constituent cells is similar to that produced by surface roughness; therefore, the replacement of the target particle by an array of large constituent cells has the same effect, qualitatively, as incorporating additional small-scale surface roughness on the particle. Such a modification of the target particle can be desirable in many practical applications of DDA when irregularly shaped particles are considered. When applying DDA to conductive, nonspherical particles, the insufficient description of the electromagnetic field in the skin layer does not lead to a violation of the Maxwell equations, although it has a visible but nonmajor influence on the light-scattering properties of the target.

Coherence, power laws, and the negative polarization surge.

We develop a second-order ray-tracing model that includes the constructive interference of reciprocal rays and Fresnel reflections from an inverse gamma-type distribution function of path length to predict the negative polarization branch seen in some astronomical bodies. We expect that such a path-length distribution might resemble the path lengths undertaken by rays incident upon some astronomical bodies. The resulting negative polarization is largely wavelength independent and depends primarily on the power law in the path-length distribution, which coincides with some observations.

Laboratory experiments on backscattering from regolith samples.

The investigation of the backscattering peak has applications in the surface texture characterization of asteroids and planetary surfaces. Laboratory experiments are important because they give an opportunity for systematic variation and comparison of samples. A backscattering experiment from regolith samples, which uses a laser light source and a beam splitter to reach the smallest phase angles, is presented. Measurements at zero and small phase angles for Sahara sand and meteorite rocks are made, and the preliminary results are presented in comparison with the phase curve observed for asteroid 69 Hesperia. The results are applicable to the further interpretation of the coherent backscattering opposition effect.