Fast calculation technique for scattering in T-matrix method

In scattering calculations using the T-matrix method, the calculation of the T-matrix involves multiplication and inversion of matrices. These two types of matrix operations are time-consuming, especially for the matrices with large size. Petrov et al. [D. Petrov, Y. Shkuratov, G. Videen, Opt. Lett. 32 (2007) 1168] proposed an optimized matrix inversion technique, which suggests the inversion of two matrices, each of which contains half the number of rows. This technique reduces time-consumption significantly. On the basis of this approach, we propose another fast calculation technique for scattering in the T-matrix method, which obtains the scattered fields through carrying out only the operations between matrices and the incident field coefficient. Numerical results show that this technique can decrease time-consumption by more than half that of the optimized matrix inversion technique by Petrov et al.     

A new path-integral representation of the T-matrix in potential scattering

We employ the method used by Barbashov and collaborators in Quantum Field Theory to derive a path-integral representation of the T-matrix in nonrelativistic potential scattering which is free of functional integration over fictitious variables as was necessary before. The resulting expression serves as a starting point for a variational approximation applied to high-energy scattering from a Gaussian potential. Good agreement with exact partial-wave calculations is found even at large scattering angles. A novel path-integral representation of the scattering length is obtained in the low-energy limit.     

T-matrix method and its applications to electromagnetic scattering by particles: A current perspective

This note serves as a short introduction to the reprint of our article “T-matrix computations of light scattering by nonspherical particles: a review” (JQSRT 1996; 55:535-75). We first discuss the motivation for writing that article and explain its historical context. This is followed by a short overview of more recent developments.     

The T-matrix technique for calculations of scattering properties of ensembles of randomly oriented particles with different size

We develop a modification of the T-matrix method that allows efficient studies of scattering properties of ensembles of independent irregular particles of different size. The advantage of the modification is quick calculations using the so-called shape-matrices (Sh-matrices), which allow more rapid calculations of scattering by particles of different size and can be used for averaging scattering properties over particle size. To illustrate the advantage we calculate the scattering-angle dependence of the intensity and degree of linear polarization of ensembles of cubes and Chebyshev particles of different size using both the new and traditional methods. Our time savings in calculating scattering properties for the particles with the new methodology is approximately a factor of ten when calculating scattering properties of one hundred of the same type of particles with different size parameter. As can be anticipated, increasing the size interval results in a smoothing of the structure of the photometric curves and a decrease in the linear polarization.     

T-matrix method formulation applied to the study of flexural waves scattering from a through obstacle in a plate

Elastic wave scattering in a flat thin plate hosting a through obstacle of arbitrary closed form is examined using a numerical technique based on the T-matrix approach, which is applied to describe of flexural waves in plates. The limiting cases of a hole and a rigid obstacle are considered. The vibrations of the plate are described by the Kirchhoff model. The far field backscattered amplitude as a function of wave frequency for inclusions of elliptic, triangular and square form with rounded corners is analysed numerically. Comparison of present results for circular obstacles with the analytical solutions obtained by other authors show excellent agreement.     

The dependence of the photodisintegration cross section on the off-shell T-matrix

Nucleon-nucleon scattering phase shifts determine the diagonal element of the transition matrix. The off-diagonal elements are not completely arbitrary but have conditions imposed on them by the range and the tail of the potential. Electromagnetic interaction can also be used to place restrictions on the off-diagonal elements. We find that the cross section of the deuteron photodisintegration is sensitive to the off-shell transition matrix. The integrated cross section can be varied by as much as 30 % or more, and the matrix element for the El transition by a factor of 2. While the matrix element for the photodisintegration depends on the off-shell elements of the T-matrix, it cannot be used to discriminate between alternative off-shell T-matrices. We have constructed classes of different off-shell T-matrices, which produce identical photo-disintegration cross sections and other two-body scattering and bound-state properties.     

Comprehensive T-matrix reference database: A 2007-2009 update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database.     

A multiple sphere T-matrix Fortran code for use on parallel computer clusters

A general-purpose Fortran-90 code for calculation of the electromagnetic scattering and absorption properties of multiple sphere clusters is described. The code can calculate the efficiency factors and scattering matrix elements of the cluster for either fixed or random orientation with respect to the incident beam and for plane wave or localized-approximation Gaussian incident fields. In addition, the code can calculate maps of the electric field both interior and exterior to the spheres. The code is written with message passing interface instructions to enable the use on distributed memory compute clusters, and for such platforms the code can make feasible the calculation of absorption, scattering, and general EM characteristics of systems containing several thousand spheres.     

A molecular T-matrix approach to calculating Low-Energy Electron Diffraction intensities for ordered molecular adsorbates

We present a novel approach to calculating Low-Energy Electron Diffraction (LEED) intensities for ordered molecular adsorbates. First, the intra-molecular multiple scattering is computed to obtain a non-diagonal molecular T-matrix. This is then used to represent the entire molecule as a single scattering object in a conventional LEED calculation, where the Layer Doubling technique is applied to assemble the different layers, including the molecular ones. A detailed comparison with conventional layer-type LEED calculations is provided to ascertain the accuracy of this scheme of calculation. Advantages of this scheme for problems involving ordered arrays of molecules adsorbed on surfaces are discussed.     

On some operators in multichannel scattering

We investigate the strong limit of an operator valued sequence used in other form in the nonrelativistic theory of multichannel scattering, and also some of its consequences.     

A study of radiative properties of fractal soot aggregates using the superposition T-matrix method

We employ the numerically exact superposition T-matrix method to perform extensive computations of scattering and absorption properties of soot aggregates with varying state of compactness and size. The fractal dimension, D_f, is used to quantify the geometrical mass dispersion of the clusters. The optical properties of soot aggregates for a given fractal dimension are complex functions of the refractive index of the material m, the number of monomers N_S, and the monomer radius a. It is shown that for smaller values of a, the absorption cross section tends to be relatively constant when D_f<2 but increases rapidly when D_f>2. However, a systematic reduction in light absorption with D_f is observed for clusters with sufficiently large N_S, m, and a. The scattering cross section and single-scattering albedo increase monotonically as fractals evolve from chain-like to more densely packed morphologies, which is a strong manifestation of the increasing importance of scattering interaction among spherules. Overall, the results for soot fractals differ profoundly from those calculated for the respective volume-equivalent soot spheres as well as for the respective external mixtures of soot monomers under the assumption that there are no electromagnetic interactions between the monomers. The climate-research implications of our results are discussed.     

T-matrix approach to calculating circular polarization of aggregates made of optically active materials

Optical activity is a typical property of the biological materials where left-handed amino-acids and right-handed carbohydrates dominate (so called homochirality). Observationally, optically active materials reveal themselves through the circular polarization in the light they scatter. Thus, circular polarization produced by the optically active particles can serve as a biomarker. It is known that biological (e.g. colonies of bacteria) and pre-biological (e.g. dust in comets) particles often have a complex structure that can be modeled presenting them as aggregates of small monomers. This motivated the development of the T-matrix code presented in this paper, which enables calculation of the scattering matrix - including circular polarization - of the light scattered by aggregated optically active particles. The code can be used for modeling the light scattering by biological objects (e.g. colonies of bacteria, blood cells) and for interpretation of the circular polarization produced by the cosmic dust that contains (pre)biological organic, e.g. comet dust or planetary aerosols.     

A coordinate-space study of kowalski's three-term separable approximation for the fully off-shell two-nucleon T-matrix

Using elementary coordinate-space methods, we show that a three-term separable approximate fully off-shell T-matrix proposed by Kowalski can be reduced to a simpler expression. This T-matrix incorporates off-shell unitarity exactly, is exact half off the energy shell, and is free from the spurious poles that arise in the Noyes approximation. However, numerical tests employing the wave-function model of Picker, Redish, and Stephenson show that for realistic ¹S_o interactions, the Noyes approximation is more accurate than Kowalski's approximation except near the spurious pole at 250 MeV. We give a plausible explanation of this result.     

Transfer-matrix for the multichannel scattering problem

A transfer-matrix for the multichannel scattering problem is obtained. The elements of this matrix are expressed in terms of transmission and reflection amplitudes. On the basis of the matrix for a system of N localized and nonoverlapped scattering centers the recurrent equations for the transfermatrix elements are derived and the initial conditions are defined.     

New formula for calculating coulomb energies in the liquid-drop model

A new transformation of double volume integrals into double surface integrals is presented. A simple regular method for deriving integrands in a surface integral is proposed. This method is used to calculate the Coulomb energy of a nucleus within the model of a liquid drop with a sharp boundary. Numerical results obtained on the basis of the new formula are compared with those calculated by one of the formulas employed previously.     

Immersing method for the multichannel scattering problem

We propose a method for solution of the multichannel scattering problem. In particular, the scattering problem is considered when a particle makes a finite motion in the transverse direction of the scattering. In this case the scattering becomes multichannel, which is connected with the presence of discrete energy levels of the transverse movement of the particle. For the case of two-channel scattering, the problem is formulated up to the end. A method for determination of scattering amplitudes for the potential V=V(x,y) is proposed.     

Quark cluster model and the R-matrix theory treatment of NN scattering

A microscopic quark cluster model has been developed for six-quark states consisting of two s³ quark clusters. The consequences of channel nonorthogonality and existence of Pauliforbidden states are investigated explicitly by solving the eigenvalue problem of the resonating group method (RGM) kernel. Since the RGM kernels needed are all available, the form of the six-quark states given in this paper is very suited to detailed RGM calculations. A rigorous treatment based on the R-matrix theory has been carried out to obtain NN phase shifts. The spin-spin term of the quark-quark interaction favors states of higher color-spin symmetry. This explains the larger change caused by the hidden color states in the ³S₁ phase shifts than in the ¹S₀ phase shifts. Phase shifts calculated with inclusion of the delta and hidden color states are still too repulsive. It is pointed out that there arises a subtle problem in adding the one-boson exchange potential by hand to the RGM equation.     

On the maximum of the k-dependent susceptibility for fixed k and T > Tc

We study the maximum occuring in ${}_{\mathrm{x}}\text{(}\text{k}\text{, T)}$ for fixed k at a temperature T_max > T_c. We examine the predictions of both Ornstein-Zernike and ?-expansions for the position of this maximum. These are compared with the exact results fot the two-dimensional Ising model.