Numerical simulation of Gaussian beam scattering by complex particles of arbitrary shape and structure

An efficient numerical method based on the surface integral equations is introduced to simulate the scattering of Gaussian beam by complex particles that consist of an arbitrarily shaped host particle and multiple internal inclusions of arbitrary shape. In particular, the incident focused Gaussian beam is described by the Davis fifth-order approximate expressions in combination with rotation defined by Euler angles. The established surface integral equations are discretized with the method of moments, where the unknown equivalent electric and magnetic currents induced on the surfaces of the host particle and the internal inclusions are expanded using the Rao–Wilton–Glisson (RWG) basis functions. The resultant matrix equations are solved by using the parallel conjugate gradient method. The proposed numerical method is validated and its capability illustrated in several characteristic examples.     

Multiple scattering by particles embedded in an absorbing medium. 2. Radiative transfer equation

This paper continues a systematic theoretical analysis of electromagnetic scattering by a group of arbitrarily sized, shaped, and oriented particles embedded in an absorbing, homogeneous, isotropic, and unbounded medium. The previously developed microphysical approach is used to derive the generalized form of the radiative transfer equation (RTE) applicable to a large group of sparsely, randomly, and uniformly distributed particles. The derivation of the RTE directly from the macroscopic Maxwell equations yields unambiguous and definitive analytical expressions for the participating quantities and thereby fully resolves the lasting controversy caused by the conflicting outcomes of several phenomenological approaches.     

Analytical techniques for addressing forward and inverse problems of light scattering by irregularly shaped particles

Understanding light scattering by nonspherical particles is crucial in modeling the transport of light in realistic structures such as biological tissues. We report the application of novel analytical approaches based on modified Wentzel-Kramers-Brillouin and equiphase-sphere methods that facilitate accurate characterization of light scattering by a wide range of irregularly shaped dielectric particles. We also demonstrate that these approaches have the potential to address the inverse-scattering problem by means of a spectral analysis of the total scattering cross section of arbitrarily shaped particles.     

Absorption and scattering properties of arbitrarily shaped particles in the Rayleigh domain: A rapid computational method and a theoretical foundation for the statistical approach

We provide a theoretical foundation for the statistical approach for computing the absorption properties of particles in the Rayleigh domain. We present a general method based on the discrete dipole approximation to compute the absorption and scattering properties of particles in the Rayleigh domain. The method allows to separate the geometrical aspects of a particle from its material properties. Doing the computation of the optical properties of a particle once, provides them for any set of refractive indices, wavelengths and orientations. This allows for fast computations of e.g. absorption spectra of arbitrarily shaped particles. Other practical applications of the method are in the interpretation of atmospheric and radar measurements as well as computations of the scattering matrix of small particles as a function of the scattering angle. In the statistical approach, the optical properties of irregularly shaped particles are represented by the average properties of an ensemble of particles with simple shapes. We show that the absorption cross section of an ensemble of arbitrarily shaped particles with arbitrary orientations can always be uniquely represented by the average absorption cross section of an ensemble of spheroidal particles with the same composition and fixed orientation. This proves for the first time that the statistical approach is generally viable in the Rayleigh domain.     

Radiation force caused by scattering, absorption, and emission of light by nonspherical particles

General formulas for computing the radiation force exerted on arbitrarily oriented and arbitrarily shaped nonspherical particles due to scattering, absorption, and emission of electromagnetic radiation are derived. For randomly oriented particles with a plane of symmetry, the formula for the average radiation force caused by the particle response to external illumination reduces to the standard Debye formula derived from the Lorenz–Mie theory, whereas the average radiation force caused by emission vanishes.     

Combinatorics of Matrix Factorizations and Integrable Systems

We study relations between the eigenvectors of rational matrix functions on the Riemann sphere. Our main result is that for a subclass of functions that are products of two elementary blocks it is possible to represent these relations in a combinatorial–geometric way using a diagram of a cube. In this representation, vertices of the cube represent eigenvectors, edges are labeled by differences of locations of zeroes and poles of the determinant of our matrix function, and each face corresponds to a particular choice of a coordinate system on the space of such functions. Moreover, for each face this labeling encodes, in a neat and efficient way, a generating function for the expressions of the remaining four eigenvectors that label the opposing face of the cube in terms of the coordinates represented by the chosen face. The main motivation behind this work is that when our matrix is a Lax matrix of a discrete integrable system, such generating functions can be interpreted as Lagrangians of the system, and a choice of a particular face corresponds to a choice of the direction of the motion.     

Effective potential method for active particles

ABSTRACT

We investigate the steady-state properties of an active fluid modelled as an assembly of soft repulsive spheres subjected to Gaussian coloured noise. Such a noise captures one of the salient aspects of active particles, namely the persistence of their motion and determines a variety of novel features with respect to familiar passive fluids. We show that within the so-called multidimensional unified coloured noise approximation, recently introduced in the field of active matter, the model can be treated by methods similar to those employed in the study of standard molecular fluids. The system shows a tendency of the particles to aggregate even in the presence of purely repulsive forces because the combined action of coloured noise and interactions enhances the effective friction between nearby particles. We also discuss whether an effective two-body potential approach, which would allow to employ methods similar to those of density functional theory, is appropriate. The limits of such an approximation are discussed.     

Amplitude and phase of light scattered by micro-scale aggregates of dielectric spheres: Comparison between theory and microwave analogy experiments

Light scattering is a useful diagnostic tool for characterization of particles. Direct scattering measurements for arbitrarily shaped micro-scale particles is difficult due to small-scale limitations. Microwave analogy is a convenient approach to realize such measurements as it enables realization of analogous experiments with larger model particles in a spectral domain where wavelengths are on centimeter scale. In the present study a test model analogous to light scattering by a micro-scale aggregate of dielectric spheres was constructed and experimentally characterized in the microwave regime. Measured amplitude and phase of the scattered field were compared with theoretical predictions obtained from quasi-exact multiple-scattering T-matrix method and discrete dipole approximation (DDA). Excellent agreement demonstrates the validities of both the experiment and the models.     

Scattering of Gaussian beam by arbitrarily shaped inhomogeneous particles

A hybrid finite element-boundary integral method is applied to characterize the scattering of an arbitrarily incident-focused Gaussian beam by arbitrarily shaped inhomogeneous particles. Specifically, the Davis–Barton fifth-order approximation in combination with rotation Euler angles is used to represent the arbitrarily incident Gaussian beams. The finite element method is employed to formulate the fields in the interior region of the inhomogeneous particle, while the boundary integral equation is applied to represent the fields in the exterior region. The interior and exterior fields are coupled by means of the field continuity conditions. To reduce the computational burden, the frontal method and the multilevel fast multipole algorithm are adopted to solve the resultant matrix equation. Numerical results for differential scattering cross sections of several selected inhomogeneous particles are presented and can be served as further study on this subject.     

An Approximate Approach for Systems of Singular Volterra Integral Equations Based on Taylor Expansion

In this article, an extended Taylor expansion method is proposed to estimate the solution of linear singular Volterra integral equations systems. The method is based on combining the m-th order Taylor polynomial of unknown functions at an arbitrary point and integration method, such that the given system of singular integral equations is converted into a system of linear equations with respect to unknown functions and their derivatives. The required solutions are obtained by solving the resulting linear system. The proposed method gives a very satisfactory solution, which can be performed by any symbolic mathematical packages such as Maple, Mathematica, etc. Our proposed approach provides a significant advantage that the m-th order approximate solutions are equal to exact solutions if the exact solutions are polynomial functions of degree less than or equal to m. We present an error analysis for the proposed method to emphasize its reliability. Six numerical examples are provided to show the accuracy and the efficiency of the suggested scheme for which the exact solutions are known in advance.     

Focal field computation of an arbitrarily polarized beam using fast Fourier transforms

In this work the vectorial diffraction theory of Richards and Wolf is extended to compute the focal field components of an arbitrarily polarized beam using fast Fourier transform (FFT) operations. Here the arbitrarily polarized pupil function is written as the vector sum of two mutually perpendicularly polarized pupil functions. The FFT based focal field expressions are particularly useful to compute the focal field components of pupil functions without a simple analytical form. We have then applied these expressions to simulate the effects of Zernike mode aberrations on the point spread functions of a number of important cylindrical-vector beam profiles such as radially and azimuthally polarized and helical light beams.     

A review of the numerical investigation on the scattering of Gaussian beam by complex particles

The study of light scattering by various particles is an active and important subject of research with myriad practical applications. During the years the scattering of plane wave by various particles has been investigated extensively. In recent years, with the development of laser sources and the tremendous expansion of their application, there has been a growing interest in the study of light scattering by various particles illuminated by a focused Gaussian beam. Since the analytical methods are only suitable for the analysis of Gaussian beam scattering by some regular particles, for complex particles with arbitrary shape and structure, one has to resort to the numerical methods. In this article, we review the recent numerical investigation on the scattering of Gaussian beam by systems of complex particles, including arbitrarily shaped conducting particles, dielectric particles, composite particles with inclusions, as well as random discrete particles and fractal soot aggregates. The essential formulations of the proposed numerical methods are outlined and the numerical results for some complex particles are also presented. This review is expected to provide useful help for the study of the interaction between the laser beams and the complex particles.     

Mechanism for the appearance of a high-efficiency Brownian motor with fluctuating potential

A version of a Brownian motor (system generating a unidirectional motion of a Brownian particle in an asymmetric fluctuating potential) is considered for the case where the potential consists of an asymmetric periodic component undergoing random shifts by a half-period L with a certain frequency and the potential of an external force F. The high efficiency of such a motor (ratio of the useful work against the load force F to the energy imparted to the particle due to the potential shifts) is due to a high and narrow barrier, as well as to a smooth arbitrarily shaped potential relief repeated with an energy shift on both half-periods L. Simple analytic expressions are obtained for the flux and efficiency as functions of the load force over a wide frequency range.     

刚体绕任意轴转动惯量的简便计算

利用微元求和的方法给出几种刚体对于过质心任意轴的转动惯量公式，并由此引出转动惯量与刚体几何对称性的联系，以及与转动惯量有关的几何形体等价性的一些结果．     

Electromagnetic Invisibility of Elliptic Cylinder Cloaks

Structures with unique electromagnetic properties are designed based on the approach of spatial coordinate transformations of Maxwell's equations. Thisapproach is applied to scheme out invisible elliptic cylinder cloaks, which provide more feasibility for cloaking arbitrarily shaped objects. The transformation expressions for the anisotropic material parameters and the field distribution are derived. The cloaking performances of ideal and lossy elliptic cylinder cloaks are investigated by finite element simulations. It is found that the cloaking performance will degrade in the forward direction withincreasing loss.     

A review of the effects of light scattering on the dynamics of irregularly shaped dust grains in the Solar System

The various aspects of the interaction of electromagnetic radiation with cosmic dust particles are discussed. In particular, attention is paid to discrepancies between optical and physical behavior of realistically shaped particles and volume equivalent homogeneous spheres. The dynamical evolution of morphologically non-identical particles which are driven by gravity, electromagnetic radiation and the Lorentz forces can dramatically differ. Although spherical particles often enable analytical calculations, an orbital evolution of spheres cannot be considered as a representative evolution for real cosmic dust particles. The effect of electromagnetic radiation on the motion of dust grains plays a crucial role here. While irregularly shaped interstellar dust particles may be captured in the Solar System, the spherical particles will not survive due to close encounters with the Sun. Spherical grains can be captured almost only in the evaporation region (in the vicinity of the Sun), where they are destroyed due to high temperatures. The spherical dust particles ejected from comets will monotonously inspiral toward the Sun subject to the Poynting-Robertson effect. However, the non-spherical particles of the same origin may be temporarily stabilized at some heliocentric distances and thus their lifetime may be much longer than that for the Mie spheres. Some dust particles may also be captured in mean-motion resonances with planets (commensurability resonances). While spherical particles are always characterized by the secular decrease of the semi-major axes near mean-motion resonances, this may not be true for non-spherical particles. Resonant captures of arbitrarily shaped dust grains exist for exterior and interior mean-motion resonances with planets.     

一种求解任意入射波束的波束因子的方法

提出了一种将任意入射波束因子用矢量波函数展开的方法，根据波束在球坐标系中的展开形式，以及球谐矢量函数与非球坐标系的波矢量函数之间的关系，推导出任意入射波束在相应坐标系中的波束因子的理论表述形式.以椭球坐标系为例，介绍了将离轴的入射波束，用椭球矢量波函数展开的波束因子求解方法，此方法还可应用于柱坐标、椭柱坐标系中波束因子的求解，为研究粒子对任意入射波束的散射打下了基础. 关键词： 波束因子 广义米理论 光散射     

Hydrodynamics and time correlation functions for cellular automata

Hydrodynamic excitations in lattice gas cellular automata are described in terms of equilibrium time correlation functions for the local conserved variables. For large space and time scales the linearized hydrodynamic equations are obtained to Navier-Stokes order. Exact expressions for the associated susceptibilities and transport coefficients are identified in terms of correlation functions. The general form of the time correlation functions for conserved densities in the hydrodynamic limit is given and illustrated by some examples suitable for comparison with computer simulation. The transport coefficients are related to time correlation functions for the conserved fluxes in a way analogous to the Green-Kubo expressions for continuous fluids. The general results are applied for a one-component fluid and several types of binary diffusion. Also discussed are the effects of unphysical slow modes such as staggered particle or momentum densities.     

A Simplified Algorithm for Impedance Calculation of Arbitrarily Shaped Radiators

It is well known that a computationally efficient model for calculation of radiation impedance of an arbitrarily shaped piston has been developed. We simplify the proposed algorithm by using geometric characteristics and intrinsic relationship between the analytic expressions. As an example, the method accuracy is illustrated and the radiation impedance of a right-angled triangular piston is calculated. The numerical results are in good agreement with that obtained directly by the quadruple integral method.