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.

Use of in situ FT-Raman spectroscopy to study the kinetics of the transformation of carbamazepine polymorphs

The solid-state transformation of carbamazepine from form III to form I was examined by Fourier Transform Raman spectroscopy. Using a novel environmental chamber, the isothermal conversion was monitored in situ at 130◦C, 138◦C, 140◦C and 150◦C. The rate of transformation was monitored by taking the relative intensities of peaks arising from two C H bending modes; this approach minimised errors due to thermal artefacts and variations in power intensities or scattering efficiencies from the samples in which crystal habit changed from a characteristic prism morphology (form III) to whiskers (form I). The solid-state transformation at the different temperatures was fitted to various solid-state kinetic models of which four gave good fits, thus indicating the complexity of the process which is known to occur via a solid–gas–solid mechanism. Arrhenius plots from the kinetic models yielded activation energies from 344 kJ mol−1 to 368 kJ mol−1 for the transformation. The study demonstrates the value of a rapid in situ analysis of drug polymorphic type which can be of value for at-line in-process control.

Determination of orientations of aromatic groups in self-assembled peptide fibrils by polarised Raman spectroscopy

In this paper we describe a novel combination of Raman spectroscopy, isotope editing and X-ray scattering as a powerful approach to give detailed structural information on aromatic side chains in peptide fibrils. The orientation of the tyrosine residues in fibrils of the peptide YTIAALLSPYS with respect to the fibril axis has been determined from a combination of polarised Raman spectroscopy and X-ray diffraction measurements. The Raman intensity of selected tyrosine bands collected at different polarisation geometries is related to the values and orientation of the Raman tensor for those specific vibrations. Using published Raman tensor values we solved the relevant expressions for both of the two tyrosine residues present in this peptide. Ring deuteration in one of the two tyrosine side chains allowed for the calculation to be performed individually for both, by virtue of the isotopic shift that eliminates band overlapping. Sample disorder was taken into account by obtaining the distribution of orientations of the samples from X-ray diffraction experiments. The results provide previously unavailable details about the molecular conformation of this peptide, and demonstrate the value of this approach for the study of amyloid fibrils.

A collocation method for high frequency scattering by convex polygons

We consider the problem of scattering of a time-harmonic acoustic incident plane wave by a sound soft convex polygon. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the computational cost required to achieve a prescribed level of accuracy grows linearly with respect to the frequency of the incident wave. Recently Chandler–Wilde and Langdon proposed a novel Galerkin boundary element method for this problem for which, by incorporating the products of plane wave basis functions with piecewise polynomials supported on a graded mesh into the approximation space, they were able to demonstrate that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency. Here we propose a related collocation method, using the same approximation space, for which we demonstrate via numerical experiments a convergence rate identical to that achieved with the Galerkin scheme, but with a substantially reduced computational cost.

Efficient evaluation of highly oscillatory acoustic scattering surface integrals

We consider the approximation of some highly oscillatory weakly singular surface integrals, arising from boundary integral methods with smooth global basis functions for solving problems of high frequency acoustic scattering by three-dimensional convex obstacles, described globally in spherical coordinates. As the frequency of the incident wave increases, the performance of standard quadrature schemes deteriorates. Naive application of asymptotic schemes also fails due to the weak singularity. We propose here a new scheme based on a combination of an asymptotic approach and exact treatment of singularities in an appropriate coordinate system. For the case of a spherical scatterer we demonstrate via error analysis and numerical results that, provided the observation point is sufficiently far from the shadow boundary, a high level of accuracy can be achieved with a minimal computational cost.

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.

Water wave scattering by finite arrays of circular structures

The scattering of small amplitude water waves by a finite array of locally axisymmetric structures is considered. Regions of varying quiescent depth are included and their axisymmetric nature, together with a mild-slope approximation, permits an adaptation of well-known interaction theory which ultimately reduces the problem to a simple numerical calculation. Numerical results are given and effects due to regions of varying depth on wave loading and free-surface elevation are presented.