Nonlinear scatterers in nonlinear media: Born approximation and local perturbation method

Wave scattering by inhomogeneous nonlinear particles placed in a nonlinear host medium is studied. In the case of weak scattering the scattering indicatrix is calculated in the first Born approximation. Scattering from small nonlinear particles loaded in a medium with saturation of permittivity is studied by the local perturbation method. A small perturbation method is developed for nonlinear equations of rather general type with a random, intensity-dependent, scattering potential.     

Effective medium approximation for strongly nonlinear media

The effective medium approximation is one of the most popular approximations used for calculating the effective coefficients of linear composite media. When the same approach is applied to the case of power-law nonlinear composite media the obtained expression contains a function whose values are unknown. In order to determine the form of this function and to calculate some coefficients related to it, we calculate the electric field for the case of a single inclusion. The numerical solution is based on the relaxation method for solving differential equations, but involves some modifications due to the nonlinearity. After the solution of the differential equation the function and the coefficients are calculated and examined. The results differ considerably from those obtained earlier by simple approximations. An erratum to this article is available at .     

Linear approximation in pulse propagation problems in nonlinear media

It seems obvious that the manifestation of nonlinear effects (such as shock wave and soliton formation) is possible only at comparatively high levels of nonlinearity. Nevertheless, there are well-known examples where a disturbance which is as small as one pleases ceases in time to be linear, and must be described in nonlinear terms. This, e.g., is the situation for the well-known Korteweg-de Vries (KdV) equations, which for an impulsive initial condition of definite sign has an asymptotic solution described by a set of solitons and by a rapidly oscillating wave packet. There will in fact always be at least one soliton asymptotically, no matter how small the initial pulse amplitude [1]. Another example is shock wave (triangular pulse) attenuation at large distances (or times) within the framework of the Burgers equation; the wave remains a shock wave even though its amplitude decreases to zero [1]. It is thus necessary to indicate a criterion which might be used to establish the nature of the asymptotic solution. The aim of the present note is the definition of such a criterion; that is, a relation between nonlinearity and dispersion (dissipation) in which nonlinear effects are manifested only above some threshold level.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 21, No. 11, pp. 1706–1708, November, 1978.     

Diffusion in lattice Lorentz gases : Nearest scatterers approximation

The expressions for diffusion coefficients and for velocity autocorrelation functions of lattice Lorentz gases are derived both in the nearest scatterers and Boltzmann approximations. The results are obtained for linear chain, square, triangular, simple cubic, body centred cubic, face centred cubic and face centred hyper cubic lattices. The diffusion coefficients are compared with those from the effective medium approximation for the square lattice and with computer simulation results for triangular, simple cubic and body centred cubic lattices.     

Directed paths in random media, function approximation and nonlinear regression

Directed paths in random media are used for function approximation and nonlinear regression. Input data are taken as random pinning. Wandering directed paths, subject to internal elasticity and interaction with the pinning, when pinned, are solutions to these problems.     

A variational perturbation approximation method in Tsallis non-extensive statistical physics

For the generalized statistical mechanics based on the Tsallis entropy, a variational perturbation approximation method with the principle of minimal sensitivity is developed by calculating the generalized free energy up to the third order in variational perturbation expansion. The approximation up to the first order amounts to a variational approach which covers the variational method developed by E.K. Lenzi, L.C. Malacarne, R.S. Mendes [Phys. Rev. Lett. 80 (1998) 218] and the approximations up to higher orders can systematically improve variational results. As an illustrated example, the generalized free energy for a classical harmonic oscillator (considered in the Lenzi's joint work) are calculated up to the third order, and the resultant approximations up to the first, second, and third orders are numerically compared with the exact result.     

Singular perturbation method for short waves in random media

The problem of wave propagation in a randomly inhomogeneous medium is considered on the basis of the parabolic equation approximation. The method of asymptotic expansions construction in powers of the radius of correlation of the random media for the moments of the wave field are proposed.     

Singular perturbation method for short waves in random media

Abstract

The problem of wave propagation in a randomly inhomogeneous medium is considered on the basis of the parabolic equation approximation. The method of asymptotic expansions construction in powers of the radius of correlation of the random media for the moments of the wave field are proposed.     

Experimental three-dimensional fluorescence reconstruction of diffuse media by use of a normalized Born approximation

We present a normalized Born expansion that facilitates fluorescence reconstructions in turbid, tissuelike media. The algorithm can be particularly useful for tissue investigations of fluorochrome distributionin vivo, since it does not require absolute photon-field measurements or measurements before contrast-agent administration. This unique advantage can be achieved only in fluorescence mode. We used this algorithm to three-dimensionally image and quantify an indocyanine fluorochrome phantom, using a novel fluorescence tomographic imager developed for animals.     

On Born approximation in black hole scattering

A massless field propagating on spherically symmetric black hole metrics such as the Schwarzschild, Reissner–Nordström and Reissner–Nordström–de Sitter backgrounds is considered. In particular, explicit formulae in terms of transcendental functions for the scattering of massless scalar particles off black holes are derived within a Born approximation. It is shown that the conditions on the existence of the Born integral forbid a straightforward extraction of the quasi normal modes using the Born approximation for the scattering amplitude. Such a method has been used in literature. We suggest a novel, well defined method, to extract the large imaginary part of quasinormal modes via the Coulomb-like phase shift. Furthermore, we compare the numerically evaluated exact scattering amplitude with the Born one to find that the approximation is not very useful for the scattering of massless scalar, electromagnetic as well as gravitational waves from black holes.     

Lamé函数和非线性演化方程的扰动方法

利用小扰动方法对非线性演化方程作展开得到原始方程的各级近似方程.应用Jacobi椭圆函数展开法求得了零级近似方程的准确解,并由此得到一级近似方程和二级近似方程分别满足齐次Lamé方程和非齐次Lamé方程,应用Lamé函数和Jacobi椭圆函数展开法可以分别求得一级近似方程和二级近似方程的准确解.这样,就求得了非线性演化方程的多级准确解.     

Three-dimensional electromagnetic nonlinear inversion in layered media by a hybrid diagonal tensor approximation: Stabilized biconjugate gradient fast Fourier transform method

This paper presents an efficient three-dimensional nonlinear electromagnetic inversion method in a multilayered medium for radar applications where the object size is comparable to the wavelength. In the first step of this two-step inversion algorithm, the diagonal tensor approximation is used in the Born iterative method. The solution of this approximate inversion is used as an initial guess for the second step in which further inversion is carried out using a distorted Born iterative method. Since the aim of the second step is to improve the accuracy of the inversion, a full-wave solver, the stabilized biconjugate-gradient fast Fourier transform algorithm, is used for forward modelling. The conjugate-gradient method is applied at each inversion iteration to minimize the functional cost. The usage of an iterative solver based on the FFT algorithm and the developed recursive matrix method combined with an interpolation technique to evaluate the layered medium Green's functions rapidly, makes this method highly efficient. An inversion problem with 32 768 complex unknowns can be solved with 1% relative error by using a simple personal computer. Several numerical experiments for arbitrarily located source and receiver arrays are presented to show the high efficiency and accuracy of the proposed method.     

Three-dimensional electromagnetic nonlinear inversion in layered media by a hybrid diagonal tensor approximation: Stabilized biconjugate gradient fast Fourier transform method

This paper presents an efficient three-dimensional nonlinear electromagnetic inversion method in a multilayered medium for radar applications where the object size is comparable to the wavelength. In the first step of this two-step inversion algorithm, the diagonal tensor approximation is used in the Born iterative method. The solution of this approximate inversion is used as an initial guess for the second step in which further inversion is carried out using a distorted Born iterative method. Since the aim of the second step is to improve the accuracy of the inversion, a full-wave solver, the stabilized biconjugate-gradient fast Fourier transform algorithm, is used for forward modelling. The conjugate-gradient method is applied at each inversion iteration to minimize the functional cost. The usage of an iterative solver based on the FFT algorithm and the developed recursive matrix method combined with an interpolation technique to evaluate the layered medium Green's functions rapidly, makes this method highly efficient. An inversion problem with 32 768 complex unknowns can be solved with 1% relative error by using a simple personal computer. Several numerical experiments for arbitrarily located source and receiver arrays are presented to show the high efficiency and accuracy of the proposed method.     

Revisiting the normalized Born approximation: effects of scattering

The normalized Born approximation has been suggested as a ratiometric method in fluorescence molecular tomography (FMT) applications, to account for heterogeneity variations. The method enabled practical inversions, as it offered fluorescence reconstruction accuracy over a wide range of absorption heterogeneity, while also accounting for unknown experimental factors, such as the various system gains and losses. Yet it was noted that scattering variations affect the robustness and accuracy. Herein we decompose the effects of absorption and scattering and capitalize on the recent development of hybrid FMT/x-ray computed tomography imaging methods to proposed amendments to the method, which improve the overall accuracy of the approach.     

Modifications in high energy higher order Born approximation

A systematic study is made to find out the differential scattering cross-section in the case of electron-atom collisions. The first and the second Born terms ofO(1/k _i ) are calculated in the framework of Yates high energy higher order Born approximation. The second Born term ofO(1/k _i ² ) is calculated using the second order Wallace term, the third term is calculated using the Glauber-eikonal series of Yates. The method is applied to the elastic scattering of electrons by atomic hydrogen in the energy range 100–400 eV and by helium for energies 200 eV and 400 eV. Comparison is made with other theoretical results and the experimental data.     

Localized modes in orientation-disordered one-dimensional media with uniaxial scatterers

Localized modes in one-dimensional (1D) media with uniaxial scatterers that are assumed to be order in spatial location but disorder in spatial orientation of their optical axis are investigated. Based on the holistic effect model in random laser, I.e., the random laser is due to the interaction of the complex localized modes in active random media with local aperiodic quasi-structure with appropriate pump light, a physical model on this type of random media is found. Its disorder degree is defined by D = no/ne. Then, the typical transmission spectrum through the random media and the light field intensity distribution corresponding to the defect modes in photonic band-gap are calculated numerically by means of the transfer matrix method, and the condition that the localized mode appears is discussed. Results show that the medium disorder plays an important role in determining the lightwave state. The localized state appears when the medium disorder is strong enough, and a new mechanism creating random laser phenomenon is brought forward.     

Nonlinear images of scatterers in chirped pulsed laser beams

The bandwidth and the duration of incident pulsed beam are proved to play important roles in modifying the nonlinear image of amplitude-type scatterer.It is found that the initially positive chirp-type bandwidth can suppress the nonlinear image,while the negative one can enhance it,and that both effects are inversely proportional to the incident pulse duration.Numerical simulations further demonstrate that the location of nonlinear image is at the conjugate plane of the scatterer and that,for negatively pre-chirped pulsed beam,the nonlinear image peak intensity can be higher than that in the corresponding monochromatic case under certain conditions.Moreover the effect of group velocity dispersion on nonlinear image is found to be similar to that of chirp-type bandwidth.     

An improved decoupling approximation method for nonlinear granular composites

An improved decoupling approximation is proposed to estimate the spatial average of ∥E∥^β, where ∥E∥ denotes the magnitude of the local electric field E and β (β ≥ 3) is the nonlinear exponent. As an example, we recalculate the nonlinear susceptibilities of weakly nonlinear composites and find that our present results are very close to the simulation data.     

Fast perturbation Monte Carlo method for photon migration in heterogeneous turbid media

We present a two-step Monte Carlo (MC) method that is used to solve the radiative transfer equation in heterogeneous turbid media. The method exploits the one-to-one correspondence between the seed value of a random number generator and the sequence of random numbers. In the first step, a full MC simulation is run for the initial distribution of the optical properties and the "good" seeds (the ones leading to detected photons) are stored in an array. In the second step, we run a new MC simulation with only the good seeds stored in the first step, i.e., we propagate only detected photons. The effect of a change in the optical properties is calculated in a short time by using two scaling relationships. By this method we can increase the speed of a simulation up to a factor of 1300 in typical situations found in near-IR tissue spectroscopy and diffuse optical tomography, with a minimal requirement for hard disk space. Potential applications of this method for imaging of turbid media and the inverse problem are discussed.