Symmetry-break in Voronoi tessellations

We analyse in a common framework the properties of the Voronoi tessellations resulting from regular 2D and 3D crystals and those of tessellations generated by Poisson distributions of points, thus joining on symmetry breaking processes and the approach to uniform random distributions of seeds. We perturb crystalline structures in 2D and 3D with a spatial Gaussian noise whose adimensional strength is α and analyse the statistical properties of the cells of the resulting Voronoi tessellations using an ensemble approach. In 2D we consider triangular, square and hexagonal regular lattices, resulting into hexagonal, square and triangular tessellations, respectively. In 3D we consider the simple cubic (SC), body-centred cubic (BCC), and face-centred cubic (FCC) crystals, whose corresponding Voronoi cells are the cube, the truncated octahedron, and the rhombic dodecahedron, respectively. In 2D, for all values α>0, hexagons constitute the most common class of cells. Noise destroys the triangular and square tessellations, which are structurally unstable, as their topological properties are discontinuous in α=0. On the contrary, the honeycomb hexagonal tessellation is topologically stable and, experimentally, all Voronoi cells are hexagonal for small but finite noise with α<0.12. Basically, the same happens in the 3D case, where only the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. In both 2D and 3D cases, already for a moderate amount of Gaussian noise (α>0.5), memory of the specific initial unperturbed state is lost, because the statistical properties of the three perturbed regular tessellations are indistinguishable. When α>2, results converge to those of Poisson-Voronoi tessellations. In 2D, while the isoperimetric ratio increases with noise for the perturbed hexagonal tessellation, for the perturbed triangular and square tessellations it is optimised for specific value of noise intensity. The same applies in 3D, where noise degrades the isoperimetric ratio for perturbed FCC and BCC lattices, whereas the opposite holds for perturbed SCC lattices. This allows for formulating a weaker form of the Kelvin conjecture. By analysing jointly the statistical properties of the area and of the volume of the cells, we discover that also the cells shape heavily fluctuates when noise is introduced in the system. In 2D, the geometrical properties of n-sided cells change with α until the Poisson-Voronoi limit is reached for α>2; in this limit the Desch law for perimeters is shown to be not valid and a square root dependence on n is established, which agrees with exact asymptotic results. Anomalous scaling relations are observed between the perimeter and the area in the 2D and between the areas and the volumes of the cells in 3D: except for the hexagonal (2D) and FCC structure (3D), this applies also for infinitesimal noise. In the Poisson-Voronoi limit, the anomalous exponent is about 0.17 in both the 2D and 3D case. A positive anomaly in the scaling indicates that large cells preferentially feature large isoperimetric quotients. As the number of faces is strongly correlated with the sphericity (cells with more faces are bulkier), in 3D it is shown that the anomalous scaling is heavily reduced when we perform power law fits separately on cells with a specific number of faces.

On the stability radius of a generalized state-space system

The concept of “distance to instability” of a system matrix is generalized to system pencils which arise in descriptor (semistate) systems. Difficulties arise in the case of singular systems, because the pencil can be made unstable by an infinitesimal perturbation. It is necessary to measure the distance subject to restricted, or structured, perturbations. In this paper a suitable measure for the stability radius of a generalized state-space system is defined, and a computable expression for the distance to instability is derived for regular pencils of index less than or equal to one. For systems which are strongly controllable it is shown that this measure is related to the sensitivity of the poles of the system over all feedback matrices assigning the poles.

Graph-based algorithms for comparison and prediction of household-level energy use profiles

We present an efficient graph-based algorithm for quantifying the similarity of household-level energy use profiles, using a notion of similarity that allows for small time–shifts when comparing profiles. Experimental results on a real smart meter data set demonstrate that in cases of practical interest our technique is far faster than the existing method for computing the same similarity measure. Having a fast algorithm for measuring profile similarity improves the efficiency of tasks such as clustering of customers and cross-validation of forecasting methods using historical data. Furthermore, we apply a generalisation of our algorithm to produce substantially better household-level energy use forecasts from historical smart meter data.

Glimpsing regular lattice arrangements of primary rat hippocampal astrocytes on ultra-thin nodes of Parylene-C

In our seminal work, we reported how the biomaterial Parylene-C has the unique ability to coerce neurons and glial cells to migrate to and then grow in straight lines along serum coated rectangular parylene-C structures mounted on an oxidised silicon substrate. In this brief communication, we report how astrocyte cell bodies, from the dissociated postnatal rat hippocampus, can now to be successfully localised on an ultra-thin 13nm layer of parylene-C mounted on oxidised silicon (Figure 1). What is extremely interesting about this finding is that the astrocyte processes extended mainly in horizontal and vertical directions from the cell body thus creating a regular lattice network of individual cells. In addition, they comfortably extended a 50μm gap (equivalent to ~ 10 cell body diameters) to connect to adjacent astrocytes on neighbouring Parylene-C structures. This was found to occur repeatedly on circular geometries of 20μm diameter. In comparison to our previous work [1], we have decreased the thickness of the parylene-C structures by a factor of 10, to allow such technology to be able to be utilised for passive electrode design that requires extremely thin structures such as these. Thus, being able to culture astrocytes in regular lattice networks will pave the way for precise monitoring and stimulation of such ensembles via multi-electrode arrays, allowing a closer insight into their dynamic behaviour and their network properties.

Black carbon vertical profiles strongly affect its radiative forcing uncertainty

The impact of black carbon (BC) aerosols on the global radiation balance is not well constrained. Here twelve global aerosol models are used to show that at least 20% of the present uncertainty in modeled BC direct radiative forcing (RF) is due to diversity in the simulated vertical profile of BC mass. Results are from phases 1 and 2 of the global aerosol model intercomparison project (AeroCom). Additionally, a significant fraction of the variability is shown to come from high altitudes, as, globally, more than 40% of the total BC RF is exerted above 5 km. BC emission regions and areas with transported BC are found to have differing characteristics. These insights into the importance of the vertical profile of BC lead us to suggest that observational studies are needed to better characterize the global distribution of BC, including in the upper troposphere.

Lesions impairing regular versus irregular past tense production

We investigated selective impairments in the production of regular and irregular past tense by examining language performance and lesion sites in a sample of twelve stroke patients. A disadvantage in regular past tense production was observed in six patients when phonological complexity was greater for regular than irregular verbs, and in three patients when phonological complexity was closely matched across regularity. These deficits were not consistently related to grammatical difficulties or phonological errors but were consistently related to lesion site. All six patients with a regular past tense disadvantage had damage to the left ventral pars opercularis (in the inferior frontal cortex), an area associated with articulatory sequencing in prior functional imaging studies. In addition, those that maintained a disadvantage for regular verbs when phonological complexity was controlled had damage to the left ventral supramarginal gyrus (in the inferior parietal lobe), an area associated with phonological short-term memory. When these frontal and parietal regions were spared in patients who had damage to subcortical (n = 2) or posterior temporo-parietal regions (n = 3), past tense production was relatively unimpaired for both regular and irregular forms. The remaining (12th) patient was impaired in producing regular past tense but was significantly less accurate when producing irregular past tense. This patient had frontal, parietal, subcortical and posterior temporo-parietal damage, but was distinguished from the other patients by damage to the left anterior temporal cortex, an area associated with semantic processing. We consider how our lesion site and behavioral observations have implications for theoretical accounts of past tense production.

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.

Software evolution: a graph based model

A model based on graph isomorphisms is used to formalize software evolution. Step by step we narrow the search space by an informed selection of the attributes based on the current state-of-the-art in software engineering and generate a seed solution. We then traverse the resulting space using graph isomorphisms and other set operations over the vertex sets. The new solutions will preserve the desired attributes. The goal of defining an isomorphism based search mechanism is to construct predictors of evolution that can facilitate the automation of ’software factory’ paradigm. The model allows for automation via software tools implementing the concepts.

A novel combined biomarker including plasma carotenoids, vitamin C, and ferric reducing antioxidant power is more strongly associated with fruit and vegetable intake than the individual components

BACKGROUND: Monitoring of fruit and vegetable (F&V) intake is fraught with difficulties. Available dietary assessment methods are associated with considerable error, and the use of biomarkers offers an attractive alternative. Few studies to date have examined the use of plasma biomarkers to monitor or predict the F&V intake of volunteers consuming a wide range of intakes from both habitual F&V and manipulated diets. OBJECTIVE: This study tested the hypothesis that an integrated biomarker calculated from a combination of plasma vitamin C, cholesterol-adjusted carotenoid concentration and Ferric Reducing Antioxidant Power (FRAP) had more power to predict F&V intake than each individual biomarker. METHODS: Data from a randomized controlled dietary intervention study [FLAVURS (Flavonoids University of Reading Study); n = 154] in which the test groups observed sequential increases of 2.3, 3.2, and 4.2 portions of F&Vs every 6 wk across an 18-wk period were used in this study. RESULTS: An integrated plasma biomarker was devised that included plasma vitamin C, total cholesterol-adjusted carotenoids, and FRAP values, which better correlated with F&V intake (r = 0.47, P < 0.001) than the individual biomarkers (r = 0.33, P < 0.01; r = 0.37, P < 0.001; and r = 0.14, respectively; P = 0.099). Inclusion of urinary potassium concentration did not significantly improve the correlation. The integrated plasma biomarker predicted F&V intake more accurately than did plasma total cholesterol-adjusted carotenoid concentration, with the difference being significant at visit 2 (P < 0.001) and with a tendency to be significant at visit 1 (P = 0.07). CONCLUSION: Either plasma total cholesterol-adjusted carotenoid concentration or the integrated biomarker could be used to distinguish between high- and moderate-F&V consumers. This trial was registered at www.controlled-trials.com as ISRCTN47748735.

Software evolution: a graph based model

A model based on graph isomorphisms is used to formalize software evolution. Step by step we narrow the search space by an informed selection of the attributes based on the current state-of-the-art in software engineering and generate a seed solution. We then traverse the resulting space using graph isomorphisms and other set operations over the vertex sets. The new solutions will preserve the desired attributes. The goal of defining an isomorphism based search mechanism is to construct predictors of evolution that can facilitate the automation of ’software factory’ paradigm. The model allows for automation via software tools implementing the concepts.

Bond graph models of DC-DC converters operating in both CCM and DCM

In this paper, Bond Graphs are employed to develop a novel mathematical model of conventional switched-mode DC-DC converters valid for both continuous and discontinuous conduction modes. A unique causality bond graph model of hybrid models is suggested with the operation of the switch and the diode to be represented by a Modulated Transformer with a binary input and a resistor with fixed conductance causality. The operation of the diode is controlled using an if-then function within the model. The extracted hybrid model is implemented on a Boost and Buck converter with their operations to change from CCM to DCM and to return to CCM. The vector fields of the models show validity in a wide operation area and comparison with the simulation of the converters using PSPICE reveals high accuracy of the proposed model, with the Normalised Root Means Square Error and the Maximum Absolute Error remaining adequately low. The model is also experimentally tested on a Buck topology.

Detecting regular sound changes in linguistics as events of concerted evolution

Background: Concerted evolution is normally used to describe parallel changes at different sites in a genome, but it is also observed in languages where a specific phoneme changes to the same other phoneme in many words in the lexicon—a phenomenon known as regular sound change. We develop a general statistical model that can detect concerted changes in aligned sequence data and apply it to study regular sound changes in the Turkic language family. Results: Linguistic evolution, unlike the genetic substitutional process, is dominated by events of concerted evolutionary change. Our model identified more than 70 historical events of regular sound change that occurred throughout the evolution of the Turkic language family, while simultaneously inferring a dated phylogenetic tree. Including regular sound changes yielded an approximately 4-fold improvement in the characterization of linguistic change over a simpler model of sporadic change, improved phylogenetic inference, and returned more reliable and plausible dates for events on the phylogenies. The historical timings of the concerted changes closely follow a Poisson process model, and the sound transition networks derived from our model mirror linguistic expectations. Conclusions: We demonstrate that a model with no prior knowledge of complex concerted or regular changes can nevertheless infer the historical timings and genealogical placements of events of concerted change from the signals left in contemporary data. Our model can be applied wherever discrete elements—such as genes, words, cultural trends, technologies, or morphological traits—can change in parallel within an organism or other evolving group.

Mycolactone-dependent depletion of endothelial cell thrombomodulin is strongly associated with fibrin deposition in buruli ulcer lesions

A well-known histopathological feature of diseased skin in Buruli ulcer (BU) is coagulative necrosis caused by the Mycobacterium ulcerans macrolide exotoxin mycolactone. Since the underlying mechanism is not known, we have investigated the effect of mycolactone on endothelial cells, focussing on the expression of surface anticoagulant molecules involved in the protein C anticoagulant pathway. Congenital deficiencies in this natural anticoagulant pathway are known to induce thrombotic complications such as purpura fulimans and spontaneous necrosis. Mycolactone profoundly decreased thrombomodulin (TM) expression on the surface of human dermal microvascular endothelial cells (HDMVEC) at doses as low as 2ng/ml and as early as 8hrs after exposure. TM activates protein C by altering thrombin’s substrate specificity, and exposure of HDMVEC to mycolactone for 24 hours resulted in an almost complete loss of the cells’ ability to produce activated protein C. Loss of TM was shown to be due to a previously described mechanism involving mycolactone-dependent blockade of Sec61 translocation that results in proteasome-dependent degradation of newly synthesised ER-transiting proteins. Indeed, depletion from cells determined by live-cell imaging of cells stably expressing a recombinant TM-GFP fusion protein occurred at the known turnover rate. In order to determine the relevance of these findings to BU disease, immunohistochemistry of punch biopsies from 40 BU lesions (31 ulcers, nine plaques) was performed. TM abundance was profoundly reduced in the subcutis of 78% of biopsies. Furthermore, it was confirmed that fibrin deposition is a common feature of BU lesions, particularly in the necrotic areas. These findings indicate that there is decreased ability to control thrombin generation in BU skin. Mycolactone’s effects on normal endothelial cell function, including its ability to activate the protein C anticoagulant pathway are strongly associated with this. Fibrin-driven tissue ischemia could contribute to the development of the tissue necrosis seen in BU lesions.