Tilt-invariant theory of rough-surface scattering: I

The small-slope approximation (SSA) in rough-surface scattering theory uses the surface slope as a small parameter of expansion. But, from the physical point of view, the slope may not be a restrictive parameter because we can change the slope of a surface simply by tilting the coordinate system. We present the theory of rough-surface scattering in a coordinate-invariant form. The new method, tilt-invariant approximation (TIA), leads to a different expansion that does not require that the slope of a surface be small. For a small Rayleigh parameter this approximate solution provides the correct perturbation theory, for a large Rayleigh parameter it provides the Kirchhoff approximation with several correcting terms.     

Convergence of solutions to the rough-surface scattering problem

Abstract

A new formulation of the rough-surface scattering problem is obtained that is closely linked to the Kirchhoff approximation. The governing equation is cast into a form amenable to solution by the method of successive approximations. The domain of convergence of this solution is established and is shown to apply also to the odd-ordered operator expansion, small-slope approximation and perturbation theory provided that the slope of the scattering surface is everywhere less than unity. The analysis is performed for scattering from one-dimensional pressure-release surfaces. Numerical examples are presented for sinusoidal and echelette gratings.     

Convergence of solutions to the rough-surface scattering problem

A new formulation of the rough-surface scattering problem is obtained that is closely linked to the Kirchhoff approximation. The governing equation is cast into a form amenable to solution by the method of successive approximations. The domain of convergence of this solution is established and is shown to apply also to the odd-ordered operator expansion, small-slope approximation and perturbation theory provided that the slope of the scattering surface is everywhere less than unity. The analysis is performed for scattering from one-dimensional pressure-release surfaces. Numerical examples are presented for sinusoidal and echelette gratings.     

An improved formulation of coherent forward scatter from random rough surfaces

Abstract

The operator expansion method is known to give accurate numerical results for scattering from individual surfaces that are too complicated for other methods. It is less widely appreciated that the method can be applied to random surfaces as well. The simplest application is the modelling of mean forward scatter from a homogeneous Gaussian ensemble of surfaces. To leading order in the admittance operator, the formula for the scalar Dirichlet boundary includes an exponential form in the roughness correlation function. The scalar Neumann boundary adds terms involving the gradients of the exponential form. These factors modestly alter the magnitude and advance the phase of the coherent scatter relative to the conventional one-point (Kirchhoff) approximation when the significant surface correlation scales are comparable to the radiation wavelength. Narrow troughs in the surface undulations ‘repel’ the radiation and effectively elevate and flatten the mean surface. These results are reliable over a wide range of surface amplitudes and correlation scales, provided the slope times Rayleigh height (Dirichlet problem) and slope (Neumann problem) are not large.     

Extension of the local curvature approximation to third order and improved tilt invariance

The second-order local curvature approximation (LCA2) is a theory of rough surface scattering that reproduces fundamental low and high frequency limits in a tilted frame of reference. Although the existing LCA2 model provides agreement with the first order small perturbation method up to the first order in surface tilt, results reported in this paper produce a new formulation of the model that achieves consistency with perturbation theory to first order in surface height and arbitrary order in surface tilt. In addition, extension of the modified LCA to third order is presented, and allows the theory to match the second-order small perturbation method to arbitrary order in surface tilt. Crucial to the development of the theory are a set of identities involving relationships among the small perturbation method (i.e. low frequency) and Kirchhoff approximation (i.e. high frequency) kernels; a set of new identities obtained in our derivations is also presented. Sample results involving 3D electromagnetic scattering from penetrable rough surfaces, as well as 2D scattering from Dirichlet sinusoidal gratings, are provided to compare the new results with the existing LCA2 model and with other rough surface scattering theories.     

X-ray scattering from a randomly rough surface

Abstract

On the basis of the method of reduced Rayleigh equations we present a simple and reciprocal theory of the coherent and incoherent scattering of x-rays from one- and two-dimensional randomly rough surfaces, that appears to be free from the limitations of earlier theories of such scattering based on the Born and distorted-wave Born approximations. In our approach, the reduced Rayleigh equation for the scattering amplitude(s) is solved perturbatively, with the small parameter of the theory η(ω) = 1 - ε(ω), where ε(ω) is the dielectric function of the scattering medium. The magnitude of η(ω) for x-rays is in the range from 10^?6 to 10^?3, depending on the wavelength of the x-rays. The contributions to the mean differential reflection coefficient from the coherent and incoherent components of the scattered x-rays are calculated through terms of second order in η(ω). The resulting expressions are valid to all orders in the surface profile function. The results for the incoherent scattering display a Yoneda peak when the scattering angle equals the critical angle for total internal reflection from the vacuum-scattering medium interface for a fixed angle of incidence, and when the angle of incidence equals the critical angle for total internal reflection for a fixed scattering angle. The approach used here may also be useful in theoretical studies of the scattering of electromagnetic waves from randomly rough dielectric-dielectric interfaces, when the difference between the dielectric constants on the two sides of the interface is small.     

A reciprocal phase-perturbation theory for rough-surface scattering

Abstract

We study the scattering of a scalar plane wave from a two-dimensional, randomly rough surface, on which the Dirichlet boundary condition is satisfied. The scattering amplitude is obtained in the form of the Fourier transform of an exponential, in which the exponent is written as an expansion in powers of the surface profile function. It is shown that the latter expansion can be written in such a way that the corresponding scattering matrix is manifestly reciprocal. Numerical results for the reflectivity, and for the contribution to the mean differential reflection coefficient from the incoherent component of the scattered field, are obtained and compared with the predictions of small-amplitude perturbation theory and the Kirchhoff approximation. As the wavelength of the incident wave is varied continuously the results of the phase-perturbation theory change continuously from those of the small-amplitude perturbation theory to those of the Kirchhoff approximation.     

Radar scatter cross sections for two-dimensional random rough surfaces – full wave solutions and comparisons with experiments

Abstract

In this paper, the full wave expressions for the radar scattering cross sections for two-dimensional random rough surfaces are obtained. The rough-surface height/slope correlations are accounted for in this analysis. Analytical and numerical comparisons of the full wave solution with the small perturbation and physical optics solutions are made for isotropic, homogeneous random rough surfaces with Gaussian probability density function. The full wave results are also compared with experimental results.     

An extension of the small-slope approximation for rough surface scattering

Solutions are derived for scattering from a rough one-dimensional pressure-release surface in the form of a functional series in the surface slope. These solutions are obtained by solving an integral equation of the first kind for the surface potential to obtain a representation for the scattering amplitude. It is shown that the subsequent expansion of terms occurring in the scattering amplitude to obtain a functional series in the slope does not yield a unique result. The result obtained contains a free parameter that may be arbitrarily selected. Thus, this result is an extension or generalization of the small-slope approximation of Voronovich (1985 Sov. Phys.-JETP 62 65-70) that differs at second order in the slope from his result. It is also shown that the free parameter can be selected such that each term of the functional series is reciprocal and exhibits a limiting grazing angle dependence consistent with the requirements of flux conservation and the absence of boundary waves. A new formulation of the leading terms of the small-slope expansions is derived and is used to explore the conditions under which the two expansions reduce to the Kirchhoff approximation. Finally, a numerical example is presented to demonstrate that the extended approximation provides corrections that are important for near grazing scatter.     

Electron density of states and electrical conductivity of ordered alloys: Extension beyond the single-band approximation, taking into account scattering by clusters

A method is developed for going beyond the single-band approximation and taking into account scattering by clusters. The method is based on a cluster expansion of the averaged Green’s function of the alloy. It is shown that the contributions of scattering processes diminish with respect to a certain small parameter as the number of particles in the cluster increases. The prominent characteristics of the electronic structure and electrical conductivity of ordered alloys are analyzed numerically in the diagonal disorder approximation in the multiband s-d model. Fiz. Tverd. Tela (St. Petersburg) 39, 401–411 (March 1997)     

Formulation of rough-surface scattering theory in terms of phase factors and approximate solutions based on this formulation

Abstract

This paper presents the formulation of rough-surface scattering theory in which the bounded phase shift factors, ζ(r, α) ζ exp[iαζ(r)], replace the elevation, ζ(r). Both the Dirichlet and the Neumann problems are considered. The integral equations for secondary surface sources are obtained that contain only this phase function in their kernels.

The Neumann (iterative) series for the solutions of the integral equations thus derived are functional Taylor series in powers of L(r, α), not in powers of ζ. If we expand L(r, α) in these series in powers of ζ(r), we obtain the standard perturbation theory series. Thus, the new formulation corresponds to the partial summation of the perturbation series.

Using the Neumann series, we obtain several uniform (with respect to αζ) approximate solutions that contain, as limiting cases, Bragg scattering, the Kirchhoff approximation, and most known advanced approximations.

In the case of random surface z = ζ(r), these new expansions contain the function ζ(r) only in the exponents, and, therefore, the result of averaging can be expressed only in terms of the characteristic functions of the multivariate probability distribution of elevations.     

Pulse broadening of enhanced backscattering from rough surfaces

Abstract

Recently, we presented a study of pulse scattering by rough surfaces based on the first-order Kirchhoff approximation which is applicable to rough surfaces with RMS slope less than 0.5 and correlation distance l?λ. However, there has been an increased interest in enhanced backscattering from rough surfaces, study of which requires inclusion of the second-order Kirchhoff approximation with shadowing corrections. This paper presents a theory for the two-frequency mutual coherence function in this region and shows that the multiple scattering on the surface gives rise to an additional pulse tail in the direction of enhanced backscattering. The theory predicts pulse broadening approximately 20% greater than that caused by single scattering alone for a delta-function incident pulse and typical surface parameters. Analytical results are compared with Monte Carlo simulations and millimetre-wave experiments for the one-dimensional rough surface with RMS height 1λ and correlation distance 1λ, showing good agreement.     

Validity of rough surface backscattering models

Abstract

A study of the regions of validity for rough surface scattering models is conducted for surfaces with Gaussian and power law power spectra. Models included in the study are physical optics (PO), geometrical optics, small perturbation method and small slope approximation. The range of validity of the PO model is commonly described by a bound on the radius curvature of the surface relative to the electromagnetic wavelength. We show empirically that for backscattering the region of accuracy is more accurately described by a bound on surface slope. For surfaces with a Gaussian power spectrum, the PO model is accurate to within 2 dB for RMS surface slope values less than 0.59 cos³θ. For surfaces with a power law power spectral density, the PO model is accurate for significant slope values (RMS surface height/wavelength of the dominant spectral peak) less than 0.037 cos³θ. These conditions are valid up to approximately 30°. The regions of validity of other models in the study are also shown to be well approximated by bounds on surface slope.     

Analytical comparison between the surface current integral equation and the second-order small-slope approximation

Abstract

This paper is the third in a series discussing a new approximate bistatic model for electromagnetic scattering from perfectly conducting rough surfaces. Our previous approach supplemented the Kirchhoff model through the addition of new terms involving linear orders in slope and surface elevation differences that arise naturally from a second iteration of the surface current integral equation. This completion of the Kirchhoff was shown to provide the correct first-order small perturbation method (SPM-1) in the general bistatic context. The agreement with SPM-1 was achieved because differences of surface heights are no longer expanded in powers of surface slope. While consistent with SPM, our previous formulation fails to reconverge toward the Kirchhoff model, at some incidence and scattered angles, when the illuminated surface satisfies the high frequency roughness condition. This weakness is also shared with the first-order small slope approximation (SSA-1) which is structurally equivalent to our previous formulation where the polarization is independent of surface roughness. The second-order small slope approximation (SSA-2), which satisfies the SPM-1 and second-order small perturbation method (SPM-2) limits by construction, was shown by Voronovich to converge toward the tangent plane approximation of the Kirchhoff model under high frequency conditions. In the present paper, we show that, in addition to the linear orders in our previous model, one must now include cross-terms between slope and surface elevation to ensure convergence toward both high frequency and small perturbation limits. With the inclusion of these terms, our new formulation becomes comparable to the SSA-2 (second-order kernel) without the need to evaluate all the quadratic order slope and elevations terms. SSA-2 is more complete, however, in the sense that it guarantees convergence toward the second-order Bragg limit (SPM-2) in the fully dielectric case in addition to both SPM-1 and Kirchhoff. Our new generalization is shown to explain correctly extra depolarization in specular conditions to be caused by surface curvature and surface autocorrelation for incoherent and coherent scattering, respectively. This result will have large repercussions on the interpretation of bistatically reflected signals such as those from GPS.     

The region of validity of perturbation theory

Abstract

We present a study of the region of validity of perturbation theory applied to rough surface scattering. We solve numerically the case of a periodic surface or grating varying in one dimension. For a statistical ensemble of gratings with a sufficiently long period one may obtain a good approximation of rough surface scattering. We use this to test the validity of perturbation theory.

Only the perfect conductor case was considered. We find that as the grazing angle becomes small the perturbation result for the TE (E horizontal) polarization remains valid, while for the TM (E vertical) polarization it breaks down. The results show that the perturbation results should be used carefully when being compared with experimental data at grazing angles.     

Multiple scattering in the high-frequency limit with second-order shadowing function from 2D anisotropic rough dielectric surfaces: I. Theoretical study

Abstract

In this paper the first- and second-order Kirchhoff approximation is applied to study the backscattering enhancement phenomenon, which appears when the surface rms slope is greater than 0.5. The formulation is reduced to the geometric optics approximation in which the second-order illumination function is taken into account. This study is developed for a two-dimensional (2D) anisotropic stationary rough dielectric surface and for any surface slope and height distributions assumed to be statistically even. Using the Weyl representation of the Green function (which introduces an absolute value over the surface elevation in the phase term), the incoherent scattering coefficient under the stationary phase assumption is expressed as the sum of three terms. The incoherent scattering coefficient then requires the numerical computation of a ten- dimensional integral. To reduce the number of numerical integrations, the geometric optics approximation is applied, which assumes that the correlation between two adjacent points is very strong. The model is then proportional to two surface slope probabilities, for which the slopes would specularly reflect the beams in the double scattering process. In addition, the slope distributions are related with each other by a propagating function, which accounts for the second-order illumination function. The companion paper is devoted to the simulation of this model and comparisons with an ‘exact’ numerical method.     

Local and non-local curvature approximation: a new asymptotic theory for wave scattering

Abstract

We present a new asymptotic theory for scalar and vector wave scattering from rough surfaces which federates an extended Kirchhoff approximation (EKA), such as the integral equation method (IEM), with the first and second order small slope approximations (SSA). The new development stems from the fact that any improvement of the ‘high frequency’ Kirchhoff or tangent plane approximation (KA) must come through surface curvature and higher order derivatives. Hence, this condition requires that the second order kernel be quadratic in its lowest order with respect to its Fourier variable or formally the gradient operator. A second important constraint which must be met is that both the Kirchhoff approximation (KA) and the first order small perturbation method (SPM-1 or Bragg) be dynamically reached, depending on the surface conditions. We derive herein this new kernel from a formal inclusion of the derivative operator in the difference between the polarization coefficients of KA and SPM-1. This new kernel is as simple as the expressions for both Kirchhoff and SPM-1 coefficients. This formal difference has the same curvature order as SSA-1 + SSA-2. It is acknowledged that even though the second order small perturbation method (SPM-2) is not enforced, as opposed to the SSA, our model should reproduce a reasonable approximation of the SPM-2 function at least up to the curvature or quadratic order. We provide three different versions of this new asymptotic theory under the local, non-local, and weighted curvature approximations. Each of these three models is demonstrated to be tilt invariant through first order in the tilting vector.     

Scattering of Rayleigh waves by a statistical inhomogeneity of the mass density

The problem of the scattering of a Rayleigh wave by a surface inhomogeneity of the mass density of an isotropic solid is solved in the Born approximation of perturbation theory. The inhomogeneity is statistical with a Gaussian correlation function in the plane parallel to the surface and is deterministic with an exponentially decaying dependence on the coordinate perpendicular to the surface. Expressions are derived for the displacement fields in the scattered longitudinal (P), transverse (SV and SH), and Rayleigh (R) waves at large distances from the inhomogeneity. The Rayleigh wave energy scattering coefficients are calculated as functions of the wavelength λ, the correlation length a of the inhomogeneity, the depth d of the defective layer, and the Poisson ratio of the medium, σ. The angular distribution of the scattered Rayleigh wave energy is determined. Asymptotic expressions are obtained for the scattering coefficient in various limiting cases with respect to the parameters a/λ and λ/d. The relation between the energies in the scattered P, SV, SH, and R waves is established. The resulting equations are used to calculate the scattering coefficients numerically over a wide range of variation of the parameters a/λ, λ/d, and σ; the results are presented in the form of graphs and a table. A physical pattern of the scattering process is constructed and used as a basis for interpreting the results of the study. Fiz. Tverd. Tela (St. Petersburg) 39, 267–274 (February 1997)     

An iterative solution of the rough-surface scattering problem

Abstract

An iterative solution to the problem of scattering from a one-dimensional rough surface is obtained for the Dirichlet boundary condition. The advantages of this method are that bounds for convergence of the solution can be established and that the solution may readily be iterated to sufficiently high order in the interaction to examine the rate at which it converges. Absolute convergence of the iterative solution is also a sufficient condition for the convergence of the operator expansion method for surfaces on which the slope is everywhere less than unity. A numerical example of scattering from an echelette grating is considered, and bounds for convergence established. It is found that for scattering from such surfaces the rate at which the iterative solution converges decreases as the surface slope is increased. Corresponding results are found for the operator expansion method.     

Universal behaviour of scattering amplitudes for scattering from a plane in an average rough surface for small grazing angles

Abstract

It is shown that for scattering from a plane in an average rough surface, the scattering cross section of the range of small grazing angles of the scattered wave demonstrates a universal behaviour. If the angle of incidence is fixed (in general it should not be small), the diffusive component of the scattering cross section for the Dirichlet problem is proportional to θ² where θ is the (small) angle of elevation, and for the Neumann problem it does not depend on θ. For the backscattering case these dependences correspondingly become θ⁴ and θ°. The result is obtained from the structure of the equations that determine the scattering problem rather than by use of an approximation.