Equations for the description of wave scattering by multi-valued random surfaces

We consider a scattering theory for multi-valued rough surfaces which cannot be described by the conventional equation of the type z = ζ(x,y). Both Dirichlet and Neumann problems are analyzed. Starting with Green's theorem we obtain a representation of the scattered field, the surface integral equation, and the extinction theorem for such surfaces. In contrast to conventional theory, these equations contain three random functions x = x(u ₁,u ₂), y = y(u ₁,u ₂), and z = z(u ₁,u ₂), where u ₁ and u ₂ are the parameters describing the surface. We introduce two scattering amplitudes S _± for describing the scattered wave above and below the surface. The extinction theorem, if formulated in terms of S _?, allows us to determine S _? for an arbitrary multi-valued surface and after this it becomes possible to derive a simple integral equation for surface sources. Knowledge of the surface sources allows us to find S ₊ by integration.     

Perturbation theory of scattering from irregular bodies

A perturbation solution is derived for the following problem: A time harmonic wave of amplitude ψ, propagating in a medium with wave number k, is incident on an irregular volume V, inside of which the propagation constant k′(r) can be an arbitrary function of | r |, where r is a position vector with origin inside V. The boundary conditions are that both ψ and its normal derivative ∂ψ/∂n may be discontinuous across the surface of V. Special cases of these conditions correspond to acoustic scattering, to B-wave scattering from a dielectric cylinder, or to the classical Dirichlet (ψ = 0) or Neumann (∂ψ/∂n = 0) surface conditions. An integral equation is derived that satisfies the appropriate differential equations both outside and inside the body, and satisfies the boundary conditions as well. This equation is reduced to a set of linear algebraic equations by expansion in a certain basis set and these linear equations are then solved in a perturbation approximation for the case that the surface of the body differs from a sphere or cylinder by a small parameter λ. Comparison is made with formulae in the literature, and except for some minor discrepancies, which are here corrected, there is general agreement.     

The 2‐matrix of the spin‐polarized electron gas: contraction sum rules and spectral resolutions

The spin‐polarized homogeneous electron gas with densities ρ_↑ and ρ_↓ for electrons with spin ‘up’ (↑) and spin ‘down’ (↓), respectively, is systematically analyzed with respect to its lowest‐order reduced densities and density matrices and their mutual relations. The three 2‐body reduced density matrices γ_↑↑, γ_↓↓, γ_a are 4‐point functions for electron pairs with spins ↑↑, ↓↓, and antiparallel, respectively. From them, three functions G_↑↑(x,y), G_↓↓(x,y), G_a(x,y), depending on only two variables, are derived. These functions contain not only the pair densities according to g_↑↑(r) = G_↑uarr;(0,r), g_↓↓(r) = G_↓↓(0,r), g_a(r) = G_a(0,r) with r = | r ₁ ‐ r ₂|, but also the 1‐body reduced density matrices γ_↑ and γ_↓ being 2‐point functions according to γ_s = ρ_sf_s and f_s(r) = G_ss(r, ∞) with s = ↑,↓ and r = | r ₁ ‐ r ^′₁|. The contraction properties of the 2‐body reduced density matrices lead to three sum rules to be obeyed by the three key functions G_ss, G_a. These contraction sum rules contain corresponding normalization sum rules as special cases. The momentum distributions n_↑(k) and n_↓(k), following from f_↑(r) and f_↓(r) by Fourier transform, are correctly normalized through f_s(0) = 1. In addition to the non‐negativity conditions n_s(k),g_ss(r),g_a(r) ≥ 0 [these quantities are probabilities], it holds n_s(k) ≤ 1 and g_ss(0) = 0 due to the Pauli principle and g_a(0) ≤ 1 due to the Coulomb repulsion. Recent parametrizations of the pair densities of the spin‐unpolarized homogeneous electron gas in terms of 2‐body wave functions (geminals) and corresponding occupancies are generalized (i) to the spin‐polarized case and (ii) to the 2‐body reduced density matrix giving thus its spectral resolutions.     

Stochastic broadening of signals in distorsionless transmission lines

The paper presents a theory of a stochastic continuous transmission line in which the series inductanceLΔ[1+l(x)], series resistanceRΔ[1+r(x)], shunt capacitanceCΔ[1+c(x)], and shunt conductanceGΔ[1+g(x)] are defined as Gaussian random functions. (The continuous line is considered as a limiting case of a lumped transmision line.) The non-negative random functionsL(x),R(x),C(x), andG(x) are chosen as delta-correlated, i.e. their correlation function is of the formΘδ(x′ −x″) whereΘ is a 4×4 positive definedx-independent matrix. Propagation of a signal of Gaussian shape is analyzed. A special attention is devoted to the so-called distorsionless lines defined by the deterministic conditionR/L=G/C. As a consequence of the stochasticity of the functionsl(x),r(x),c(x), andg(x), transmitted signals do become distorted: they become broadened. An explicit formula for this broadening is derived. This work has been supported by the Slovak Grant Agency VEGA under contract No. 1/4319/97.     

Strong-coupling constant in coordinate space

The coordinate-space behavior of (vector) strong-coupling constant in the background field α_B(r) is compared with that in standard perturbation theory α_v(r). The numerically calculated two-loop coupling constant α_B(r) is shown to exceed α_v(r) by 1–5% at very small distances, r?0.02 fm, and to be in agreement with lattice measurements of the static potential. At large distances, α_B(r) approaches the freezing value at r?0.5 fm. An analytic form of α_B(r) is proposed that approximates α_B(r) with a precision ?2% in the region r?0.5 fm.     

Integral equation theory of the exchange potential,HOMO–LUMO properties,and sum rules for the exchange-correlation force

Recently, a formally exact integral equation for the exchange potential V _x(r) has been presented by the authors. In the admittedly simplistic limit in which Slater–Kohn–Sham and Hartree–Fock determinants become equal, this integral equation reduces to that given by Della Sala and Görling. Here, a proposal is made to relate, but now approximately, the formally exact equation for V _x(r) to HOMO–LUMO properties. The addition of a correlation contribution V _c(r) to V _x(r), the sum being the exchange-correlation potential V _xc(r), is finally considered, some exact properties and especially sum rules for the force ??V _xc/?r being the focus.     

Radiation transfer in an anisotropically scattering plane-parallel medium with space-dependent albedo ω(x)

A method of analysis is presented for solving radiation-transfer problems involving space-dependent albedo ω(x) for an absorbing, emitting and anisotropically scattering plane-parallel medium with reflecting boundaries. The albedo is represented in terms of Legendre polynomials in the form ω(x) = Σ^R_r=0D_rP_r(x/L), where x is the optical variable, L is the half optical-thickness of the slab, P_r(x/L) are the Legendre polynomials and D_r are known expansion coefficients. The effects of spatial variation of albedo on the reflectivity and transmissivity of a medium having a slab geometry are examined for the cases of both forward and backward anisotropic scattering over a wide range of system variables. The effects of ω(x) on the angular distribution of radiation are also shown for some representative cases.     

Studies of scattering theory using numerical methods

Abstract

Numerical simulations, using both exact and approximate methods, are used to study rough surface scattering in both the smd and large roughness regimes. This study is limited lo scattcring lrom rough one-dimensional surfaces that obey the Dirichlet boundary condition and have a Gaussian roughness spectrum. For surfdces with small roughness (kh?1, where k is the radiation wavenumber and h is the root-mean-square (RMS) Surface height), perturbation theory is known to be valid. However, it is shown numerically that when kh?1 and kl?6 (where I is the surface correlation length) the Kirchhoffapprorimation is valid except at low grazing angles, and one must sum the first three orders of perturbation theory obtain the correct result. For kh?1 and kl?1, first-order perturbation theory is accurate. In this region, the accuracy of the first two terms of the iterative series solution of the exact integral equation is examined; the first term a1 this series is the Kirchhoff approximation, It is shown numerically that lor very small kh these first two terms reduce to first-order perturbation theory. However, lor this reduction to occur, kh must be made smaller than necessdry lor first-order perturbation theory to be accurate. In the regime of large roughness (kh?1) backscattering enhancement occurs when the RMS slope is on the order of unity. Several investigators have recently shown that the second term of the iterative series solution (the double-scattering term) replicates the properties of backscattering enhancement reasonably well. However, the double-scattering term has a lundamental flaw: predictions lor the scattering cross section per unit length based on the double-scattering term increase as the surfdce length is increased. This is shown here with numerical simulations and with an approximate analytical result based on the high frequency limit. The physical significance of this finding is also discussed. The final topic is the use of the double-scattering approximation to study the mechanism for backscattering enhancement with the Dirichlet boundary condition. This mechanism is usually assumed to be interference between reciprocal scattering paths. When the interlerence between reciprocal scattering paths is removed, the enhancement is eliminated. This shows that interference between reciprocal paths is almost certainly the dominant mechanism for backscattering enhancement in the scattering regime studied.     

The convergence behaviour of expansions in the classical collision theory

In the classical collision theory the scattering angle? depends on the impact parameterb and on the kinetic energyE _r of the relative motion. This angle?(b, E _r ) is expanded for two limiting cases: 1. Expansion in powers of the potentialV(r)/E _r (momentum approximation). 2. Expansion in powers of the impact parameterb (central collision approximation). The radius of convergence of the series depends onb andE _r . It will be given for the following potentialsV(r):

$$A\left( {\frac{a}{r}} \right)^\mu ;Ae^{ - \frac{r}{a}} ;A\frac{a}{r}e^{ - \frac{r}{a}} ;A\left( {\frac{a}{r}} \right)^2 e^{ - \left( {\frac{r}{a}} \right)^2 } .$$     

带栅极纳米线冷阴极的场增强因子研究

场增强因子是体现场发射冷阴极器件性能优劣的重要参数.利用静电场理论给出了一种带栅极(normal-gated)纳米线冷阴极的场增强因子表示式β=k₁{N²·(L-d₁)²+［1/k₁+(L-d₁)］²}^1/2，且进一步分析了几何参数对场增强因子的影响.结果表明，纳米线突出栅孔的部分(L-d₁)与栅孔半径越大，则场增强因子越大；而纳米线半径越小，则场增强因子越大；当L远大于d₁时满足β∝L/r₀.其中N=N₁(k₁r₀)/N₀(k₁r₀)，N₀(k₁r₀)和N₁(k₁r₀)分别代表零阶和一阶Neumann函数，k₁=0.8936/R，R为栅孔半径，L为纳米线长度，r₀为纳米线半径，d₁表示阴极与栅极间距.     

Spin-orbit coupling constant and <r −3> in transuranic ions

The spin-orbit coupling parameter ζ, the quantum average <r ^?3>, and binding energy ? E for atomic number Z = 92, 94, 96 and net ionic charge approximately 1 and 5, have been calculated for the 5f electrons. The potential assumed in the one-electron wave function calculation is that given by the Thomas-Fermi model of a positive ion.

It is shown that the ratio ζ/<r ^?3> is fairly insensitive to changes in the surroundings of the ion, and can be used for estimating one of the parameters with the experimental value of the other.

The 5g electrons are shown to be effectively free in the considered range of Z and n.

The behaviour of the wave function as a function of n and Z is discussed.     

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.     

Estimates of general Mayer graphs. IV. On the computation of Gaussian integrals by the star-mesh transformation

We point out that, to compute by hand a Gaussian integral exp · (– _LL y _l r _ij ² d r _n+1 ...d r _n+k , where the sum runs over all linesL = (i,j) of a graph andr _ij = ¦r _i –r _j ¦, the simplest way is to use the star-mesh transformation, well known in electrical network theory. We apply this to test, on a relatively complicatedn-graph, the accuracy of an estimation method that we proposed elsewhere [Phys. Lett. 62A:295 (1977)].     

Clustering of continuous distributions of charge

An investigation is carried out of the association and clustering of equimolar mixtures of oppositely charged Gaussian charge distributions (CDs) of the form ～ exp ( ? r²?/2α²), where r is the separation between the centres of charge and α governs the extent of charge spreading (α→0 is the point charge limit). The results of molecular dynamics (MD) and Ornstein–Zernike integral equation with the mean spherical approximation (MSA) and hypernetted-chain (HNC) closures are compared for these systems. The MD and HNC radial distribution functions, g(r), agree very well for not too small α. The MD and MSA, g(r), also agree well for α ≈ 1 and greater. The potential energy per particle for the three methods also agrees well over a wider range of α values, better than might be expected from inspection of the radial distribution functions, because the dominant contributions to U come predominantly from intermediate and long distance ranges where there is good agreement between the g(r) from the MSD and HNC closures. The nature of the association and clustering of the charges as a function of α is explored through the mean nearest neighbour distance for unlike and like species and the mean and root-mean-square force. The velocity and force autocorrelation functions are also calculated; they show increasingly oscillatory behaviour in the small α limit, originating in vibrations of a pair of CDs of opposite sign.     

Initial-boundary value problem for the two-component Gerdjikov-Ivanov equation on the interval

In this paper, we apply Fokas unified method to study initial-boundary value problems for the two-component Gerdjikov-Ivanov equation formulated on the finite interval with 3×3 Lax pairs. The solution can be expressed in terms of the solution of a 3×3 Riemann-Hilbert problem. The relevant jump matrices are explicitly given in terms of three matrix-value spectral functions s (λ), S (λ) and S_L(λ), which arising from the initial values at t = 0, boundary values at x = 0 and boundary values at x = L, respectively. Moreover, The associated Dirichlet to Neumann map is analyzed via the global relation. The relevant formulae for boundary value problems on the finite interval can reduce to ones on the half-line as the length of the interval tends to infinity.     

Path integral methods for inelastic scattering

Corrections to the primitive semi-classical amplitude for multiple inelastic scattering are obtained from a path integral formulation of scattering theory. The path integrals are calculated by making an expansion about a classical orbit describing elastic scattering. Terms are collected to give a series in inverse powers of the reduced mass m of relative motion of the target and projectile. The leading term is the primitive semi-classical amplitude for multiple excitation and explicit formulae are given for the corrections of order $\text{1}\text{m}$. These are calculated in detail for a one-dimensional model. It is shown that some, but not all, of the corrections can be included by evaluating the primitive amplitude with a symmetrized orbit.     

Analysis of the Initial Slope of the Small‐Angle Scattering Correlation Function of a Particle

The small‐angle scattering correlation function of a particle γ(r) results from scattering experiments. This function possesses a well‐defined slope γ′(0) at the origin. This slope is defined by the particle volume V and the whole surface area S of the particle via γ′(0) = –S/(4V). In this paper it is demonstrated that this slope defines the mean chord length of the particle, . This theorem involves non‐convex particles, especially the case of particles with hollow parts. Consequently, for a large class of particle shapes the mean chord length is defined in terms of V and S. This extension of the Cauchy theorem is developed by closer analysis of the set covariance C(r), of the small‐angle scattering correlation function γ(r), and of the so‐called linear erosion P(r) near the origin r→0. The cases of a single hollow sphere, of two touching spheres, and of the single hollow cylinder are discussed.     

From independent particle towards collective motion for two polarized electrons on a square lattice

The two dimensional crossover from independent particle towards collective motion is studied using 2 polarized electrons (spinless fermions) interacting via a U/r Coulomb repulsion in a L×L square lattice with periodic boundary conditions and nearest neighbor hopping t. Three regimes characterize the ground state when U/t increases. Firstly, when the fluctuation Δr of the spacing r between the two particles is larger than the lattice spacing a, there is a scaling length L ₀ = π²(t/U) such that the relative fluctuation Δr/〈r〉 is a universal function of the dimensionless ratio L/L ₀, up to finite size corrections of order L^-2. L < L ₀ and L > L ₀ are respectively the limits of the free particle Fermi motion and of the correlated motion of a Wigner molecule. Secondly, when U/t exceeds a threshold U ^*(L)/t, Δr becomes smaller than a, giving rise to a correlated lattice regime where the previous scaling breaks down and analytical expansions in powers of t/U become valid. A weak random potential reduces the scaling length and favors the correlated motion. Received 28 March 2002 Published online 19 November 2002