On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

Molecular structure of 2,5-dihydropyrrole (C4NH7), obtained by gas-phase electron diffraction and theoretical calculations

The structure of 2,5-dihydropyrrole (C4NH7) has been determined by gas-phase electron diffraction (GED), augmented by the results from ab initio calculations employing third-order Moller-Plesset (MP3) level of theory and the 6-311+G(d,p) basis set. Several theoretical calculations were performed. From theoretical calculations using MP3/6-311+G(d,p) evidence was obtained for the presence of an axial (63%) (N-H bond axial to the CNC plane) and an equatorial conformer (37%) (N-H bond equatorial to the CNC plane). The five-membered ring was found to be puckered with the CNC plane inclined at 21.8 (38)° to the plane of the four carbon atoms.

A FORTRAN-90 low-energy electron diffraction program (LEED90 v1.1)

We describe a FORTRAN-90 program to compute low-energy electron diffraction I(V) curves. Plane-waves and layer doubling are used to compute the inter-layer multiple-scattering, while the intra-layer multiple-scattering is computed in the standard way expanding the wavefield on a basis of spherical waves. The program is kept as general as possible, in order to allow testing different parts of multiple-scattering calculations. In particular, it can handle non-diagonal t-matrices describing the scattering of non-spherical potentials, anisotropic vibrations, anharmonicity, etc. The program does not use old FORTRAN flavours, and has been written keeping in mind the advantage for parallelism brought forward by FORTRAN-90.

The impedance boundary value problem for the Helmholtz equation in a half-plane

We prove unique existence of solution for the impedance (or third) boundary value problem for the Helmholtz equation in a half-plane with arbitrary L∞ boundary data. This problem is of interest as a model of outdoor sound propagation over inhomogeneous flat terrain and as a model of rough surface scattering. To formulate the problem and prove uniqueness of solution we introduce a novel radiation condition, a generalization of that used in plane wave scattering by one-dimensional diffraction gratings. To prove existence of solution and a limiting absorption principle we first reformulate the problem as an equivalent second kind boundary integral equation to which we apply a form of Fredholm alternative, utilizing recent results on the solvability of integral equations on the real line in [5].

New insights into the compressibility and high-pressure stability of Ni(CN)2: a combined study of neutron diffraction, Raman spectroscopy, and inelastic neutron scattering

Nickel cyanide is a layered material showing markedly anisotropic behaviour. High-pressure neutron diffraction measurements show that at pressures up to 20.1 kbar, compressibility is much higher in the direction perpendicular to the layers, c, than in the plane of the strongly chemically bonded metal-cyanide sheets. Detailed examination of the behaviour of the tetragonal lattice parameters, a and c, as a function of pressure reveal regions in which large changes in slope occur, for example, in c(P) at 1 kbar. The experimental pressure dependence of the volume data is fitted to a bulk modulus, B0, of 1050 (20) kbar over the pressure range 0–1 kbar, and to 124 (2) kbar over the range 1–20.1 kbar. Raman spectroscopy measurements yield additional information on how the structure and bonding in the Ni(CN)2 layers change with pressure and show that a phase change occurs at about 1 kbar. The new high-pressure phase, (Phase PII), has ordered cyanide groups with sheets of D4h symmetry containing Ni(CN)4 and Ni(NC)4 groups. The Raman spectrum of phase PII closely resembles that of the related layered compound, Cu1/2Ni1/2(CN)2, which has previously been shown to contain ordered C≡N groups. The phase change, PI to PII, is also observed in inelastic neutron scattering studies which show significant changes occurring in the phonon spectra as the pressure is raised from 0.3 to 1.5 kbar. These changes reflect the large reduction in the interlayer spacing which occurs as Phase PI transforms to Phase PII and the consequent increase in difficulty for out-of-plane atomic motions. Unlike other cyanide materials e.g. Zn(CN)2 and Ag3Co(CN)6, which show an amorphization and/or a decomposition at much lower pressures (~100 kbar), Ni(CN)2 can be recovered after pressurising to 200 kbar, albeit in a more ordered form.

Structure determination from powder diffraction data

Advances made over the past decade in structure determination from powder diffraction data are reviewed with particular emphasis on algorithmic developments and the successes and limitations of the technique. While global optimization methods have been successful in the solution of molecular crystal structures, new methods are required to make the solution of inorganic crystal structures more routine. The use of complementary techniques such as NMR to assist structure solution is discussed and the potential for the combined use of X-ray and neutron diffraction data for structure verification is explored. Structures that have proved difficult to solve from powder diffraction data are reviewed and the limitations of structure determination from powder diffraction data are discussed. Furthermore, the prospects of solving small protein crystal structures over the next decade are assessed.

GDASH: a grid-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Very large scale increases in speed of execution can therefore be achieved by distributing individual DASH runs over a network of computers. The GDASH program achieves this by packaging DASH in a form that enables it to run under the Univa UD Grid MP system, which harnesses networks of existing computing resources to perform calculations.

MDASH: a multi-core-enabled program for structure solution from powder diffraction data

The simulated annealing approach to structure solution from powder diffraction data, as implemented in the DASH program, is easily amenable to parallelization at the individual run level. Modest increases in speed of execution can therefore be achieved by executing individual DASH runs on the individual cores of CPUs.