On solving coupled equations for scattering problems under the many channel and short wavelength conditions of heavy ion collisions

An iterative method for solving coupled equations in scattering problems is presented. Radical wave functions are expressed in terms of amplitudes on regular and outgoing components. This is analogous to the semi-quantal approach and gives many of its convenient features for heavy ion collisions. One solves directly for the scattering matrix without generating a complete set of solutions. The amplitudes are smoothly varying compared to the total wave functions. The iteration is non-perturbative, as each step contains the interaction to all orders. The convergence improves as the wavelength of relative motion decreases. In the short wavelength limit only one iteration is required.     

Simplified expressions of the T-matrix integrals for electromagnetic scattering

The extended boundary condition method, also called the null-field method, provides a semianalytic solution to the problem of electromagnetic scattering by a particle by constructing a transition matrix (T-matrix) that links the scattered field to the incident field. This approach requires the computation of specific integrals over the particle surface, which are typically evaluated numerically. We introduce here a new set of simplified expressions for these integrals in the commonly studied case of axisymmetric particles. Simplifications are obtained using the differentiation properties of the radial functions (spherical Bessel) and angular functions (associated Legendre functions) and integrations by parts. The resulting simplified expressions not only lead to faster computations, but also reduce the risks of loss of precision and provide a simpler framework for further analytical work.     

Analytical solutions to the electromagnetic field in a cylindrical shell excited by external axial current

The solutions to the electromagnetic field excited by a long axial current outside a conductive and magnetic cylindrical shell of finite length are studied in this paper. The more accurate analytical solutions are obtained by solving the proper boundary value problems by the separation variable method. Then the solutions are simplified according to asymptotic formulas of Bessel functions. Compared with the accurate solutions, the simplified solutions do not contain the Bessel functions and the inverse operation of the singular matrix, and can be calculated out fast by computers. The simplified solutions are more suitable for the cylindrical shell of high permeability and conductivity excited by a high frequency source. Both of the numerical results and the physical experimental results validate the simplified solutions obtained.     

库仑程函近似下非弹性散射振幅角分布的广义函数分析

在库仑程函近似框架下,引进实参量(时间扫描参数)积分,将能壳上跃迁矩阵元分为靶的结构因子和弹的扭曲因子,导出了非分波扭曲因子的主项表达式和结构因子的解析表达式,利用广义函数方法分析了电子与类氢离子的非弹性散射角分布.     

Multiple Rayleigh-Gans-Debye scattering in colloidal systems-dynamic light scattering

Equations for (depolarized) intensity auto-correlation functions including second order scattering contributions are derived. The equations obtained are quite similar to those for static light scattering¹). The geometry of the scattering system plays an important role here.An iterative procedure to correct correlation functions for double scattering is presented.Experiments on colloidal solutions containing latex and silica particles are performed, testing the theory and the iterative correction procedure for double scattering. We use cylindrical cuvettes with cylindrical incident and detected beams of radiation. The agreement between experiment and theory is found to be quite good, yielding a routine correction procedure for double scattering.     

Large-angle scattering of heavy ions: (I). General structure of elastic cross section

Closed-form expressions are derived for the differential cross section of elastic heavy-ion scattering at large angles. The derivation is based on the general form of the elastic partial-wave S-matrix in real l-space. By a generalization of analytic techniques developed in earlier work, it is shown that the large-angle scattering cross section has a universal structure involving combinations of Bessel functions and the Fourier transforms of the rapidly varying parts of the S-matrix, irrespective of their dynamical origin. Anomalous large-angle scattering is attributed to deviations of the S-matrix from its “normal strong-absorption profile”, and general conditions for backward-angle enhancement are given. Our model-independent formulation provides the framework for an “inductive” method of analyzing experimental angular distributions and excitation functions aimed at identifying, as uniquely as possible, the dynamical mechanisms that operate in large-angle heavy-ion scattering. Extensions of the formalism to inelastic scattering and transfer reactions, and applications of the analytic method, will be described in subsequent papers.     

Low-order light scattering in multiple scattering disperse media

Light scattering has been investigated in systems in which both single and higher order scattering occur. The Monte Carlo simulation technique for studying light scattering in randomly inhomogeneous, strongly scattering disperse media was employed. The reliability of the data obtained has been checked by comparing the results of the computer simulation with analytical calculations for the intensity of doubly scattered light. The first several scattering orders have been analyzed for different geometries of the optical experiment. It has been shown, in particular, that, depending on the detector aperture, the contribution of multiple scattering can vary by almost an order of magnitude.

     

Two-dimensional isotropic scattering in a semi-infinite cylindrical medium

Beginning with the integral equation for the source function, the solutions for the source function, flux and intensity at the boundary of a two-dimensional, isotropically scattering cylindrical medium are found. The incident radiation is collimated and normal to the surface of the medium and depends only on the radial coordinate. For a Bessel function boundary condition, separation of variables is used to reduce the source function integral equation to a one-dimensional equation. The resulting integral equation is shown to be the same as that for the two-dimensional planar case. Solutions for other boundary conditions are then shown to be superpositions of the Bessel function solution. Numerical results are presented for a Gaussian distribution of incident radiation which closely models a laser beam. These multiple scattering results are compared to the single scattering approximation. Also, the solution for a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on an otherwise isotropic phase function is expressed in terms of the isotropic results.     

Fast and accurate radiance calculations using truncation approximation for anisotropic scattering phase functions

The phase function for solar light scattering by large particles such as cloud droplets is strongly anisotropic due to very strong peaking in the forward direction. This creates numerical difficulties when attempting to calculate accurate reflected and transmitted radiances, which are important for remote sensing of atmospheric and surface properties. A popular approach uses the delta function to approximate the forward-scattering peak in a fraction of energy and a limited number of polynomial terms or a geometrically truncated function for the remaining fraction (so-called truncation approximations). This article compares and discusses several methods for fast and accurate calculations using truncation approximations. When using a single truncation approximation for all scattering orders, large biases appear in directions near the solar and anti-solar points. As shown here, high accuracy can be obtained using different truncation approximations depending on the order of scattering. Of particular importance is the use of phase functions close to the exact phase functions for the first few orders of scattering. Applying the method in combination with the Monte Carlo (MC) method, in which the truncation fraction for a scattering order depends on the scattering angle at the previous scattering event, obtains accurate radiance calculations under almost all geometrical and optical conditions, including in directions near the solar point. Because the method also reduces computational noise due to the MC sampling of radiance, it is useful for fast and accurate radiance calculations for cloudy atmospheres.     

Direct recursion of the ratio of Bessel functions with applications to Mie scattering calculations

The primary goal of this paper is to describe a method to directly calculate the ratios between the Bessel functions of the first and second kind. This new method is introduced and discussed under the context of Mie scattering calculations. Past work on Mie calculations has mostly relied on calculating the Bessel functions separately, then taking the ratio. This technique fails when the order of the functions becomes appreciably larger than the argument. To solve this problem, we developed a stable method to directly calculate the ratio. When applied to Mie scattering calculations, the new method improves the control of the desired accuracy and enables calculations with higher precision. Since many other wave propagation problems also involve the ratios of Bessel functions and our method provides reliable and accurate calculations over a wide range of conditions, we expect it to find applications beyond Mie calculations.     

Analytic approach to dielectric optical bent slab waveguides

A rigorous classical analytic frequency domain model of confined optical wave propagation along 2D bent slab waveguides and curved dielectric interfaces is investigated, based on a piecewise ansatz for bend mode profiles in terms of Bessel and Hankel functions. This approach provides a clear picture of the behaviour of bend modes, concerning their decay for large radial arguments or effects of varying bend radius. Fast and accurate routines are required to evaluate Bessel functions with large complex orders and large arguments. Our implementation enabled detailed studies of bent waveguide properties, including higher order bend modes and whispering gallery modes, their interference patterns, and issues related to bend mode normalization and orthogonality properties.     

Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular media

Effect of scattering on radiative heat transfer in two-dimensional rectangular media by the finite-volume method has been studied. Compared with the existing solutions, it shows that the result obtained by the finite-volume method is reliable. Furthermore, relative errors caused by the approximation that linear and nonlinear anisotropic scattering media is simplified to isotropic scattering media have been studied.     

Modeling and analysis of multiple scattering of acoustic waves in complex media: application to the trabecular bone

The integral equations that describe scattering in the media with step-rise changing parameters have been numerically solved for the trabecular bone model. The model consists of several hundred discrete randomly distributed elements. The spectral distribution of scattering coefficients in subsequent orders of scattering has been presented. Calculations were carried on for the ultrasonic frequency ranging from 0.5 to 3 MHz. Evaluation of the contribution of the first, second, and higher scattering orders to total scattering of the ultrasounds in trabecular bone was done. Contrary to the approaches that use the μCT images of trabecular structure to modeling of the ultrasonic wave propagation condition, the 3D numerical model consisting of cylindrical elements mimicking the spatial matrix of trabeculae, was applied. The scattering, due to interconnections between thick trabeculae, usually neglected in trabecular bone models, has been included in calculations when the structure backscatter was evaluated. Influence of the absorption in subsequent orders of scattering is also addressed. Results show that up to 1.5 MHz, the influence of higher scattering orders on the total scattered field characteristic can be neglected while for the higher frequencies, the relatively high amplitude interference peaks in higher scattering orders clearly occur.     

Two-dimensional linearly anisotropic scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium which scatters linearly. The incident radiation is collimated, normal to the top surface of the medium, and is dependent only on the radial coordinate. The radial variation is assumed to be a Bessel function or a Gaussian distribution. The Gaussian boundary condition is used to simulate a laser beam. Numerical results are presented in graphical and tabular forms for both boundary conditions. Results for forward and backward scattering phase functions are compared with those for isotropic scattering. A method is presented for extending these results to the problem of a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on a linear phase function.     

Conformal mapping for the Helmholtz equation: acoustic wave scattering by a two dimensional inclusion with irregular shape in an ideal fluid

The complex variables method with mapping function was extended to solve the linear acoustic wave scattering by an inclusion with sharp/smooth corners in an infinite ideal fluid domain. The improved solutions of Helmholtz equation, shown as Bessel function with mapping function as the argument and fractional order Bessel function, were analytically obtained. Based on the mapping function, the initial geometry as well as the original physical vector can be transformed into the corresponding expressions inside the mapping plane. As all the physical vectors are calculated in the mapping plane (η,η), this method can lead to potential vast savings of computational resources and memory. In this work, the results are validated against several published works in the literature. The different geometries of the inclusion with sharp corners based on the proposed mapping functions for irregular polygons are studied and discussed. The findings show that the variation of angles and frequencies of the incident waves have significant influence on the bistatic scattering pattern and the far-field form factor for the pressure in the fluid.     

Charged black rings from inverse scattering

The inverse scattering method of Belinsky and Zakharov is a powerful method to construct solutions of vacuum Einstein equations. In particular, in five dimensions this method has been successfully applied to construct a large variety of black hole solutions. Recent applications of this method to Einstein–Maxwell-dilaton (EMd) theory, for the special case of Kaluza–Klein dilaton coupling, has led to the construction of the most general black ring in this theory. In this contribution, we review the inverse scattering method and its application to the EMd theory. We illustrate the efficiency of these methods with a detailed construction of an electrically charged black ring.     

Parity violating effects in the coherent scattering of neutrons by 209Bi

We study the rotation of the spin in the coherent scattering of polarized thermal neutrons by ²⁰⁹Bi. Contributions due to direct weak scattering by the nucleus, parity admixtures into the nuclear wave functions and scattering by the electrons are calculated. We use a parity violating potential derived in the framework of the Salam-Weinberg model and several other semi-empirical potentials. The results show that the direct scattering on nuclei is the dominant contribution and some potentials predict results up to an order of magnitude larger than expected on dimensional grounds. The sign of the effect depends on the isospin properties of the potential and can distinguish between potentials with large isotensor or large isovector contributions.     

水下任意刚性散射体对Bessel波的散射特性分析

本文利用Bessel波的谐波展开式, 采用T矩阵方法的推导思路, 建立了水下任意刚性散射体在Bessel波照射下的声散射场计算公式. 以水下刚性椭球体和两端附连半球的刚性圆柱体为例, 计算了在不同波锥角β 下的反向散射形态函数, 同时, 依据镜反射波和绕行波的干涉物理模型, 给出了预报Bessel波照射下的反向散射形态函数峰峰间隔值的计算模型. 仿真结果表明本文提出的Bessel波照射下反向散射形态函数峰峰间隔值预报方法是准确有效的, 同时也说明, 本文建立的基于T矩阵法计算水下任意刚性散射体在Bessel波束下的声散射场方法是有效的, 这拓展了T矩阵法的应用领域.     

Numerical results for the thermal scattering functions

A recent formulation in radiative transfer defined the thermal scattering functions that characterize radiative transfer from a general, plane-parallel, finite medium driven solely by an internal distribution of thermal sources. Exiting diffuse intensities are expressed as space convolutions of the thermal scattering functions with any thermal source distribution. A parametric study is presented to obtain the basic structure of these scattering functions. The independent variables of these azimuthally independent functions are the direction consine μ and source location t, while the parameters are the single scattering albedo ω, total optical depth t₀, and the asymmetry factor g in the Henyey-Greenstein phase function. The basic functional trends are discussed using various parametric plots, and selected tabular results are given to allow numerical checks. The computational method is invariant imbedding. As a particular application, these functions are used in the following companion paper to obtain exiting intensities from inhomogeneous and nonisothermal media.