Bayesian estimation of recent migration rates after a spatial expansion

Approximate Bayesian computation (ABC) is a highly flexible technique that allows the estimation of parameters under demographic models that are too complex to be handled by full-likelihood methods. We assess the utility of this method to estimate the parameters of range expansion in a two-dimensional stepping-stone model, using samples from either a single deme or multiple demes. A minor modification to the ABC procedure is introduced, which leads to an improvement in the accuracy of estimation. The method is then used to estimate the expansion time and migration rates for five natural common vole populations in Switzerland typed for a sex-linked marker and a nuclear marker. Estimates based on both markers suggest that expansion occurred < 10,000 years ago, after the most recent glaciation, and that migration rates are strongly male biased.

A variational method for the calculation of spin-rovibronic energy levels of any triatomic molecule in an electronic triplet state

The triatomic spin-rovibronic variational code RVIB3 has been extended to include the effect of two uncoupled electrons, for both (3)Sigma(-) and (3)Pi (Renner-Teller) electronic states. The spin-orbital-rotational kinetic energy is included in the usual way, via terms (J+L+S). The phenomenological terms AL.S and lambda 2/3(3S(z)(2)) are introduced to reproduce the 3 spin-orbit and spin-spin splittings, respectively. Calculations are performed to evaluate the spin-rovibronic energy levels of CCO (X) over tilde (3) Sigma(-) and CCO (A) over tilde (3) Pi for which the Born-Oppenheimer potentials are derived from high-accuracy ab initio calculations.

A differential thermal expansion approach to crystal structure determination from powder diffraction data

Differential thermal expansion over the range 90-210 K has been applied successfully to determine the crystal structure of chlorothiazide from synchrotron powder diffraction data using direct methods. Key to the success of the approach is the use of a multi-data-set Pawley refinement to extract a set of reflection intensities that is more 'single-crystal-like' than those extracted from a single data set. The improvement in reflection intensity estimates is quantified by comparison with reference single-crystal intensities. (C) 2008 International Union of Crystallography Printed in Singapore - all rights reserved

Structures and negative thermal expansion properties of the one-dimensional cyanides, CuCN, AgCN and AuCN

The behaviour of the lattice parameters of HTCuCN (high-temperature form), AgCN and AuCN have been investigated as a function of temperature over the temperature range 90–490 K. All materials show one-dimensional negative thermal expansion (NTE) along the ––(M––CN)–– chain direction c (ac(HT-CuCN) ¼32.1 10–6 K1, ac(AgCN)¼23.910–6 K1 and ac(AuCN) ¼9.3106 K1 over the temperature range 90–490 K). The origin of this behaviour has been studied using RMC modelling of Bragg and total neutron diffraction data from AgCN and AuCN at 10 and 300 K. These analyses yield details of the local motions within the chains responsible for NTE. The low-temperature form of CuCN, LT-CuCN, has been studied using single-crystal X-ray diffraction. In this form of CuCN, wavelike distortions of the ––(Cu––CN)–– chains occur in the static structure, which are reminiscent of the motions seen in the RMC modelling of AgCN and AuCN, which are responsible for the NTE behaviour.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting thatmost effectors represent species-specific adaptations.

Livelihood diversification and the expansion of artisanal mining in rural Tanzania: drivers and policy implications

This paper provides an extended analysis of livelihood diversification in rural Tanzania, with special emphasis on artisanal and small-scale mining (ASM). Over the past decade, this sector of industry, which is labour-intensive and comprises an array of rudimentary and semi-mechanized operations, has become an indispensable economic activity throughout Sub-Saharan Africa, providing employment to a host of redundant public sector workers, retrenched large-scale mine labourers and poor farmers. In many of the region’s rural areas, it is overtaking subsistence agriculture as the primary industry. Such a pattern appears to be unfolding within the Morogoro and Mbeya regions of southern Tanzania, where findings from recent research suggest that a growing number of smallholder farmers are turning to ASM for employment and financial support. It is imperative that national rural development programmes take this trend into account and provide support to these people.

Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).

Blind equalization of nonminimum phase channels: higher order cumulant based algorithm

Higher order cumulant analysis is applied to the blind equalization of linear time-invariant (LTI) nonminimum-phase channels. The channel model is moving-average based. To identify the moving average parameters of channels, a higher-order cumulant fitting approach is adopted in which a novel relay algorithm is proposed to obtain the global solution. In addition, the technique incorporates model order determination. The transmitted data are considered as independently identically distributed random variables over some discrete finite set (e.g., set {±1, ±3}). A transformation scheme is suggested so that third-order cumulant analysis can be applied to this type of data. Simulation examples verify the feasibility and potential of the algorithm. Performance is compared with that of the noncumulant-based Sato scheme in terms of the steady state MSE and convergence rate.

Cumulant-based deconvolution and identification: several new families of linear equations

This paper presents several new families of cumulant-based linear equations with respect to the inverse filter coefficients for deconvolution (equalisation) and identification of nonminimum phase systems. Based on noncausal autoregressive (AR) modeling of the output signals and three theorems, these equations are derived for the cases of 2nd-, 3rd and 4th-order cumulants, respectively, and can be expressed as identical or similar forms. The algorithms constructed from these equations are simpler in form, but can offer more accurate results than the existing methods. Since the inverse filter coefficients are simply the solution of a set of linear equations, their uniqueness can normally be guaranteed. Simulations are presented for the cases of skewed series, unskewed continuous series and unskewed discrete series. The results of these simulations confirm the feasibility and efficiency of the algorithms.

Variational approach in weighted Sobolev spaces to scattering by unbounded rough surfaces

We consider the problem of scattering of time harmonic acoustic waves by an unbounded sound soft surface which is assumed to lie within a finite distance of some plane. The paper is concerned with the study of an equivalent variational formulation of this problem set in a scale of weighted Sobolev spaces. We prove well-posedness of this variational formulation in an energy space with weights which extends previous results in the unweighted setting [S. Chandler-Wilde and P. Monk, SIAM J. Math. Anal., 37 (2005), pp. 598–618] to more general inhomogeneous terms in the Helmholtz equation. In particular, in the two-dimensional case, our approach covers the problem of plane wave incidence, whereas in the three-dimensional case, incident spherical and cylindrical waves can be treated. As a further application of our results, we analyze a finite section type approximation, whereby the variational problem posed on an infinite layer is approximated by a variational problem on a bounded region.

Four-dimensional variational data assimilation for high resolution nested models

Four-dimensional variational data assimilation (4D-Var) is used in environmental prediction to estimate the state of a system from measurements. When 4D-Var is applied in the context of high resolution nested models, problems may arise in the representation of spatial scales longer than the domain of the model. In this paper we study how well 4D-Var is able to estimate the whole range of spatial scales present in one-way nested models. Using a model of the one-dimensional advection–diffusion equation we show that small spatial scales that are observed can be captured by a 4D-Var assimilation, but that information in the larger scales may be degraded. We propose a modification to 4D-Var which allows a better representation of these larger scales.