S. N. Taraskin

Universal Features of Terahertz Absorption in Disordered Materials

Using an analytical theory, experimental terahertz time-domain spectroscopy data, and numerical evidence, we demonstrate that the frequency dependence of the absorption coupling coefficient between far-infrared photons and atomic vibrations in disordered materials has the universal functional form, C(omega)=A+Bomega(2), where the material-specific constants A and B are related to the distributions of fluctuating charges obeying global and local charge neutrality, respectively.

Unveiling the Neuromorphological Space

This article proposes the concept of neuromorphological space as the multidimensional space defined by a set of measurements of the morphology of a representative set of almost 6000 biological neurons available from the NeuroMorpho database. For the first time, we analyze such a large database in order to find the general distribution of the geometrical features. We resort to McGhees biological shape space concept in order to formalize our analysis, allowing for comparison between the geometrically possible tree-like shapes, obtained by using a simple reference model, and real neuronal shapes. Two optimal types of projections, namely, principal component analysis and canonical analysis, are used in order to visualize the originally 20-D neuron distribution into 2-D morphological spaces. These projections allow the most important features to be identified. A data density analysis is also performed in the original 20-D feature space in order to corroborate the clustering structure. Several interesting results are reported, including the fact that real neurons occupy only a small region within the geometrically possible space and that two principal variables are enough to account for about half of the overall data variability. Most of the measurements have been found to be important in representing the morphological variability of the real neurons.

The Effect of Heterogeneity on Invasion in Spatial Epidemics: From Theory to Experimental Evidence in a Model System

Heterogeneity in host populations is an important factor affecting the ability of a pathogen to invade, yet the quantitative investigation of its effects on epidemic spread is still an open problem. In this paper, we test recent theoretical results, which extend the established “percolation paradigm” to the spread of a pathogen in discrete heterogeneous host populations. In particular, we test the hypothesis that the probability of epidemic invasion decreases when host heterogeneity is increased. We use replicated experimental microcosms, in which the ubiquitous pathogenic fungus Rhizoctonia solani grows through a population of discrete nutrient sites on a lattice, with nutrient sites representing hosts. The degree of host heterogeneity within different populations is adjusted by changing the proportion and the nutrient concentration of nutrient sites. The experimental data are analysed via Bayesian inference methods, estimating pathogen transmission parameters for each individual population. We find a significant, negative correlation between heterogeneity and the probability of pathogen invasion, thereby validating the theory. The value of the correlation is also in remarkably good agreement with the theoretical predictions. We briefly discuss how our results can be exploited in the design and implementation of disease control strategies.

Heterogeneity in susceptible–infected–removed (SIR) epidemics on lattices

The percolation paradigm is widely used in spatially explicit epidemic models where disease spreads between neighbouring hosts. It has been successful in identifying epidemic thresholds for invasion, separating non-invasive regimes, where the disease never invades the system, from invasive regimes where the probability of invasion is positive. However, its power is mainly limited to homogeneous systems. When heterogeneity (environmental stochasticity) is introduced, the value of the epidemic threshold is, in general, not predictable without numerical simulations. Here, we analyse the role of heterogeneity in a stochastic susceptible–infected–removed epidemic model on a two-dimensional lattice. In the homogeneous case, equivalent to bond percolation, the probability of invasion is controlled by a single parameter, the transmissibility of the pathogen between neighbouring hosts. In the heterogeneous model, the transmissibility becomes a random variable drawn from a probability distribution. We investigate how heterogeneity in transmissibility influences the value of the invasion threshold, and find that the resilience of the system to invasion can be suitably described by two control parameters, the mean and variance of the transmissibility. We analyse a two-dimensional phase diagram, where the threshold is represented by a phase boundary separating an invasive regime in the high-mean, low-variance region from a non-invasive regime in the low-mean, high-variance region of the parameter space. We thus show that the percolation paradigm can be extended to the heterogeneous case. Our results have practical implications for the analysis of disease control strategies in realistic heterogeneous epidemic systems.

Instantaneous frequency and amplitude identification using wavelets: Application to glass structure

This paper describes a method for extracting rapidly varying, superimposed amplitude-modulated and frequency-modulated signal components. The method is based upon the continuous wavelet transform (CWT) and uses a new wavelet that is a modification to the well-known Morlet wavelet to allow analysis at high resolution. In order to interpret the CWT of a signal correctly, an approximate analytic expression for the CWT of an oscillatory signal is examined via a stationary-phase approximation. This analysis is specialized for the new wavelet and the results are used to construct expressions for the amplitude and frequency modulations of the components in a signal from the transform of the signal. The method is tested on a representative, variable-frequency signal as an example before being applied to a function of interest in our subject area-a structural correlation function of a disordered material-which immediately reveals previously undetected features.

Low-frequency vibrational excitations in vitreous silica: the Ioffe-Regel limit

The behaviour of vibrational plane waves in disordered solids has been studied theoretically, particularly in the vicinity of the Ioffe-Regel crossover. We have analysed the behaviour of vibrational excitations in vitreous silica, for which we have found the Ioffe-Regel crossover frequency to be THz (for both longitudinal and transverse waves), in good agreement with experimental observations. Two methods have been used to obtain : either by determining the decay time of the plane waves, or by analysing the distribution in k-space of plane-wave components in the final scattered state of a plane wave.

Spatial decay of the single-particle density matrix in insulators: analytic results in two and three dimensions.

Analytic results for the asymptotic decay of the electron density matrix in insulators have been obtained in all three dimensions (D = 1,2,3) for a tight-binding model defined on a simple cubic lattice. The anisotropic decay length is shown to be dependent on the energy parameters of the model. The existence of the power-law prefactor, proportional, variant r(-D/2), is demonstrated.

Universal features of localized eigenstates in disordered systems

Localization–delocalization transitions occur in problems ranging from semiconductor-device physics to propagation of disease in plants and viruses on the internet. Here, we report calculations of localized electronic and vibrational eigenstates for remarkably different, mostly realistic, disordered systems and point out similar characteristics in the cases studied. We show in each case that the eigenstates may be decomposed into exponentially localized islands which may appear in many different eigenstates. In all cases, the decay length of the islands increases only modestly near the localization–delocalization transition; the eigenstates become extended primarily by proliferation (growth in number) of islands near the transition. Recently, microphotoluminescence experiments (Guillet et al 2003 Phys. Rev. B 68 045319) have imaged exciton states in disordered quantum wires, and these bear a strong qualitative resemblance to the island structure of eigenstates that we have studied theoretically.

Temporal and dimensional effects in evolutionary graph theory.

The spread in time of a mutation through a population is studied analytically and computationally in fully connected networks and on spatial lattices. The time t* for a favorable mutation to dominate scales with the population size N as N(D+1)/D in D-dimensional hypercubic lattices and as NlnN in fully-connected graphs. It is shown that the surface of the interface between mutants and nonmutants is crucial in predicting the dynamics of the system. Network topology has a significant effect on the equilibrium fitness of a simple population model incorporating multiple mutations and sexual reproduction.

Modeling the atomic structure of very high-density amorphous ice

The structure of very high-density amorphous (VHDA) ice has been modelled by positionally disordering three crystalline phases, namely ice IV, VI and XII. These phases were chosen because only they are stable or metastable in the region of the ice phase diagram where VHDA ice is formed, and their densities are comparable to that of VHDA ice. An excellent fit to the medium range of the experimentally observed pair-correlation function g(r) of VHDA ice was obtained by introducing disorder into the positions of the H2O molecules, as well as small amounts of molecular rotational disorder, disorder in the O--H bond lengths and disorder in the H--O--H bond angles. The low-k behaviour of the experimental structure factor, S(k), is also very well reproduced by this disordered-crystal model. The fraction of each phase present in the best-fit disordered model is very close to that observed in the probable crystallization products of VHDA ice. In particular, only negligible amounts of ice IV are predicted, in accordance with experimental observation.