Evolving graphs: dynamical models, inverse problems and propagation

Applications such as neuroscience, telecommunication, online social networking, transport and retail trading give rise to connectivity patterns that change over time. In this work, we address the resulting need for network models and computational algorithms that deal with dynamic links. We introduce a new class of evolving range-dependent random graphs that gives a tractable framework for modelling and simulation. We develop a spectral algorithm for calibrating a set of edge ranges from a sequence of network snapshots and give a proof of principle illustration on some neuroscience data. We also show how the model can be used computationally and analytically to investigate the scenario where an evolutionary process, such as an epidemic, takes place on an evolving network. This allows us to study the cumulative effect of two distinct types of dynamics.

Introduction. Stochastic physics and climate modelling

Finite computing resources limit the spatial resolution of state-of-the-art global climate simulations to hundreds of kilometres. In neither the atmosphere nor the ocean are small-scale processes such as convection, clouds and ocean eddies properly represented. Climate simulations are known to depend, sometimes quite strongly, on the resulting bulk-formula representation of unresolved processes. Stochastic physics schemes within weather and climate models have the potential to represent the dynamical effects of unresolved scales in ways which conventional bulk-formula representations are incapable of so doing. The application of stochastic physics to climate modelling is a rapidly advancing, important and innovative topic. The latest research findings are gathered together in the Theme Issue for which this paper serves as the introduction.

Layer-by-layer electrostatic entrapment of protein molecules on superparamagnetic nanoparticle: new strategy to enhance adsorption capacity and maintain biological activity

There is a recent interest to use inorganic-based magnetic nanoparticles as a vehicle to carry biomolecules for various biophysical applications, but direct attachment of the molecules is known to alter their conformation leading to attenuation in activity. In addition, surface immobilization has been limited to monolayer coverage. It is shown that alternate depositions of negatively charged protein molecules, typically bovine serum albumin (BSA) with a positively charged aminocarbohydrate template such as glycol chitosan (GC) on magnetic iron oxide nanoparticle surface as a colloid, are carried out under pH 7.4. Circular dichroism (CD) clearly reveals that the secondary structure of the entrapped BSA sequential depositions in this manner remains totally unaltered which is in sharp contrast to previous attempts. Probing the binding properties of the entrapped BSA using small molecules (Site I and Site II drug compounds) confirms for the first time the full retention of its biological activity as compared with native BSA, which also implies the ready accessibility of the entrapped protein molecules through the porous overlayers. This work clearly suggests a new method to immobilize and store protein molecules beyond monolayer adsorption on a magnetic nanoparticle surface without much structural alteration. This may find applications in magnetic recoverable enzymes or protein delivery.

Very large inverse problems in atmosphere and ocean modelling

For the very large nonlinear dynamical systems that arise in a wide range of physical, biological and environmental problems, the data needed to initialize a numerical forecasting model are seldom available. To generate accurate estimates of the expected states of the system, both current and future, the technique of ‘data assimilation’ is used to combine the numerical model predictions with observations of the system measured over time. Assimilation of data is an inverse problem that for very large-scale systems is generally ill-posed. In four-dimensional variational assimilation schemes, the dynamical model equations provide constraints that act to spread information into data sparse regions, enabling the state of the system to be reconstructed accurately. The mechanism for this is not well understood. Singular value decomposition techniques are applied here to the observability matrix of the system in order to analyse the critical features in this process. Simplified models are used to demonstrate how information is propagated from observed regions into unobserved areas. The impact of the size of the observational noise and the temporal position of the observations is examined. The best signal-to-noise ratio needed to extract the most information from the observations is estimated using Tikhonov regularization theory. Copyright © 2005 John Wiley & Sons, Ltd.

Dispersal capacity and diet breadth modify the response of wild bees to habitat loss

Habitat loss poses a major threat to biodiversity, and species-specific extinction risks are inextricably linked to life-history characteristics. This relationship is still poorly documented for many functionally important taxa, and at larger continental scales. With data from five replicated field studies from three countries, we examined how species richness of wild bees varies with habitat patch size. We hypothesized that the form of this relationship is affected by body size, degree of host plant specialization and sociality. Across all species, we found a positive species–area slope (z ¼ 0.19), and species traits modified this relationship. Large-bodied generalists had a lower z value than small generalists. Contrary to predictions, small specialists had similar or slightly lower z value compared with large specialists, and small generalists also tended to be more strongly affected by habitat loss as compared with small specialists. Social bees were negatively affected by habitat loss (z ¼ 0.11) irrespective of body size. We conclude that habitat loss leads to clear shifts in the species composition of wild bee communities.

Metal complexes of a tetraazacyclophane: solution and molecular modelling studies

An alternative synthetic approach to yield the compound 2,3,5,6,8,9,11,14-octahydrobenzo[1][ 1,4,7,10]tetraazacyclotetradecine (bz[14]N-4) is presented. The protonation constants of bz[14]N-4 and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, Pb2+ were determined in H2O at 25degreesC with ionic strength 0.10 mol dm(-3) in KNO3 and they were compared with structurally related macrocycles cyclam (1,4,8,11-tetraazacyclotetradecane) and cyclen (1,4,7,10-tetraazacyclododecane). These studies indicate that only 1 : 1 ( M : L) species are formed in solution, and the ligand exhibits a high affinity for larger ions such as Cd2+ and Pb2+. The X-ray study of [bz[14]N4H3](3+) shows that an inclusion compound with a chloride counter-anion is formed through NH...Cl hydrogen bonds. Spectroscopic data in solution ( electronic and NMR spectra) showed that the macrocycle adopts a planar arrangement upon metal complexation. Molecular mechanics calculations reveal that in spite of the presence of the benzene ring in the macrocyclic framework this ligand can encapsulate metal ions with different stereo-electronic sizes in square planar arrangements. Our results indicate that the presence of the benzene ring in the backbone of the bz[14]N-4 confers a coordination behaviour intermediate between that of cyclam and cyclen.

Crystal and molecular structure of the thioether-analogue of PEEK from X-ray powder data and diffraction-modelling

Analysis of X-ray powder data for the melt-crystallisable aromatic poly(thioether thioether ketone) [-S-Ar-S-Ar-CO-Ar](n), ('PTTK', Ar= 1,4-phenylene), reveals that it adopts a crystal structure very different from that established for its ether-analogue PEEK. Molecular modelling and diffraction-simulation studies of PTTK show that the structure of this polymer is analogous to that of melt-crystallised poly(thioetherketone) [-SAr-CO-Ar](n) in which the carbonyl linkages in symmetry-related chains are aligned anti-parallel to one another. and that these bridging units are crystallographically interchangeable. The final model for the crystal structure of PTTK is thus disordered, in the monoclinic space group 121a (two chains per unit cell), with cell dimensions a = 7.83, b = 6.06, c = 10.35 angstrom, beta = 93.47 degrees. (c) 2005 Elsevier Ltd. All rights reserved.

Design, synthesis and computational modelling of aromatic tweezer-molecules as models for chain-folding polymer blends

Novel 'tweezer-type' complexes that exploit the interactions between pi-electron-rich pyrenyl groups and pi-electron deficient diimide units have been designed and synthesised. The component molecules leading to complex formation were accessed readily from commercially available starting materials through short and efficient syntheses. Analysis of the resulting complexes, using the visible charge-transfer band, revealed association constants that increased sequentially from 130 to 11,000 M-1 as increasing numbers of pi-pi-stacking interactions were introduced into the systems. Computational modelling was used to analyse the structures of these complexes, revealing low-energy chain-folded conformations for both components, which readily allow close, multiple pi-pi-stacking and hydrogen bonding to be achieved. In this paper, we give details of our initial studies of these complexes and outline how their behaviour could provide a basis for designing self-healing polymer blends for use in adaptive coating systems. (C) 2008 Elsevier Ltd. All rights reserved.

The remote distractor effect in saccade programming: channel interactions and lateral inhibition

We explored the dependency of the saccadic remote distractor effect (RDE) on the spatial frequency content of target and distractor Gabor patches. A robust RDE was obtained with low-medium spatial frequency distractors, regardless of the spatial frequency of the tat-get. High spatial frequency distractors interfered to a similar extent when the target was of the same spatial frequency. We developed a quantitative model based on lateral inhibition within an oculomotor decision unit. This lateral inhibition mechanism cannot account for the interaction observed between target and distractor spatial frequency, pointing to the existence of channel interactions at an earlier level. (C) 2004 Elsevier Ltd. All rights reserved.