Method to assess and optimise dependability of complex macro-systems: application to a railway signalling system

Achieving a high degree of dependability in complex macro-systems is challenging. Because of the large number of components and numerous independent teams involved, an overview of the global system performance is usually lacking to support both design and operation adequately. A functional failure mode, effects and criticality analysis (FMECA) approach is proposed to address the dependability optimisation of large and complex systems. The basic inductive model FMECA has been enriched to include considerations such as operational procedures, alarm systems. environmental and human factors, as well as operation in degraded mode. Its implementation on a commercial software tool allows an active linking between the functional layers of the system and facilitates data processing and retrieval, which enables to contribute actively to the system optimisation. The proposed methodology has been applied to optimise dependability in a railway signalling system. Signalling systems are typical example of large complex systems made of multiple hierarchical layers. The proposed approach appears appropriate to assess the global risk- and availability-level of the system as well as to identify its vulnerabilities. This enriched-FMECA approach enables to overcome some of the limitations and pitfalls previously reported with classical FMECA approaches.

A dynamic model for stem cell homeostasis and patterning in Arabidopsis meristems.

Plants maintain stem cells in their meristems as a source for new undifferentiated cells throughout their life. Meristems are small groups of cells that provide the microenvironment that allows stem cells to prosper. Homeostasis of a stem cell domain within a growing meristem is achieved by signalling between stem cells and surrounding cells. We have here simulated the origin and maintenance of a defined stem cell domain at the tip of Arabidopsis shoot meristems, based on the assumption that meristems are self-organizing systems. The model comprises two coupled feedback regulated genetic systems that control stem cell behaviour. Using a minimal set of spatial parameters, the mathematical model allows to predict the generation, shape and size of the stem cell domain, and the underlying organizing centre. We use the model to explore the parameter space that allows stem cell maintenance, and to simulate the consequences of mutations, gene misexpression and cell ablations.

Stationary states and phase diagram for a model of the Gunn effect under realistic boundary conditions

A general formulation of boundary conditions for semiconductor-metal contacts follows from a phenomenological procedure sketched here. The resulting boundary conditions, which incorporate only physically well-defined parameters, are used to study the classical unipolar drift-diffusion model for the Gunn effect. The analysis of its stationary solutions reveals the presence of bistability and hysteresis for a certain range of contact parameters. Several types of Gunn effect are predicted to occur in the model, when no stable stationary solution exists, depending on the value of the parameters of the injecting contact appearing in the boundary condition. In this way, the critical role played by contacts in the Gunn effect is clearly established.

The wound watch: an objective staging system for wounds in the diabetic (db/db) mouse model.

As in cancer biology, in wound healing there is a need for objective staging systems to decide for the best treatment and predictors of outcome. We developed in the diabetic (db/db) wound healing model, a staging system, the "wound watch," based on the quantification of angiogenesis and cell proliferation in open wounds. In chronic wounds, there is often a lack of cellular proliferation and angiogenesis that leads to impaired healing. The wound watch addresses this by quantifying the proliferative phase of wound healing in two dimensions (cellular division and angiogenesis). The results are plotted in a two-dimensional graph to monitor the course of healing and compare the response to different treatments.

SBML qualitative models: a model representation format and infrastructure to foster interactions between qualitative modelling formalisms and tools.

BACKGROUND: Qualitative frameworks, especially those based on the logical discrete formalism, are increasingly used to model regulatory and signalling networks. A major advantage of these frameworks is that they do not require precise quantitative data, and that they are well-suited for studies of large networks. While numerous groups have developed specific computational tools that provide original methods to analyse qualitative models, a standard format to exchange qualitative models has been missing. RESULTS: We present the Systems Biology Markup Language (SBML) Qualitative Models Package ("qual"), an extension of the SBML Level 3 standard designed for computer representation of qualitative models of biological networks. We demonstrate the interoperability of models via SBML qual through the analysis of a specific signalling network by three independent software tools. Furthermore, the collective effort to define the SBML qual format paved the way for the development of LogicalModel, an open-source model library, which will facilitate the adoption of the format as well as the collaborative development of algorithms to analyse qualitative models. CONCLUSIONS: SBML qual allows the exchange of qualitative models among a number of complementary software tools. SBML qual has the potential to promote collaborative work on the development of novel computational approaches, as well as on the specification and the analysis of comprehensive qualitative models of regulatory and signalling networks.

Diffusionless first-order phase transitions in systems with frozen configurational degrees of freedom

We consider systems that can be described in terms of two kinds of degree of freedom. The corresponding ordering modes may, under certain conditions, be coupled to each other. We may thus assume that the primary ordering mode gives rise to a diffusionless first-order phase transition. The change of its thermodynamic properties as a function of the secondary-ordering-mode state is then analyzed. Two specific examples are discussed. First, we study a three-state Potts model in a binary system. Using mean-field techniques, we obtain the phase diagram and different properties of the system as a function of the distribution of atoms on the different lattice sites. In the second case, the properties of a displacive structural phase transition of martensitic type in a binary alloy are studied as a function of atomic order. Because of the directional character of the martensitic-transition mechanism, we find only a very weak dependence of the entropy on atomic order. Experimental results are found to be in quite good agreement with theoretical predictions.

Pedotransfer functions to estimate proctor test parameters under different tillage systems

The Proctor test is time-consuming and requires sampling of several kilograms of soil. Proctor test parameters were predicted in Mollisols, Entisols and Vertisols of the Pampean region of Argentina under different management systems. They were estimated from a minimum number of readily available soil properties (soil texture, total organic C) and management (training data set; n = 73). The results were used to generate a soil compaction susceptibility model, which was subsequently validated using a second group of independent data (test data set; n = 24). Soil maximum bulk density was estimated as follows: Maximum bulk density (Mg m-3) = 1.4756 - 0.00599 total organic C (g kg-1) + 0.0000275 sand (g kg-1) + 0.0539 management. Management was equal to 0 for uncropped and untilled soils and 1 for conventionally tilled soils. The established models predicted the Proctor test parameters reasonably well, based on readily available soil properties. Tillage systems induced changes in the maximum bulk density regardless of total organic matter content or soil texture. The lower maximum apparent bulk density values under no-tillage require a revision of the relative compaction thresholds for different no-tillage crops.

Avalanches in a fluctuationless first-order phase transition in a random-bond Ising model

We study the behavior of the random-bond Ising model at zero temperature by numerical simulations for a variable amount of disorder. The model is an example of systems exhibiting a fluctuationless first-order phase transition similar to some field-induced phase transitions in ferromagnetic systems and the martensitic phase transition appearing in a number of metallic alloys. We focus on the study of the hysteresis cycles appearing when the external field is swept from positive to negative values. By using a finite-size scaling hypothesis, we analyze the disorder-induced phase transition between the phase exhibiting a discontinuity in the hysteresis cycle and the phase with the continuous hysteresis cycle. Critical exponents characterizing the transition are obtained. We also analyze the size and duration distributions of the magnetization jumps (avalanches).

Dissipation in quasistatically driven disordered systems

This work presents an analysis of hysteresis and dissipation in quasistatically driven disordered systems. The study is based on the random field Ising model with fluctuationless dynamics. It enables us to sort out the fraction of the energy input by the driving field stored in the system and the fraction dissipated in every step of the transformation. The dissipation is directly related to the occurrence of avalanches, and does not scale with the size of Barkhausen magnetization jumps. In addition, the change in magnetic field between avalanches provides a measure of the energy barriers between consecutive metastable states

Long-term C-CO2 emissions and carbon crop residue mineralization in an oxisol under different tillage and crop rotation systems

Soil C-CO2 emissions are sensitive indicators of management system impacts on soil organic matter (SOM). The main soil C-CO2 sources at the soil-plant interface are the decomposition of crop residues, SOM turnover, and respiration of roots and soil biota. The objectives of this study were to evaluate the impacts of tillage and cropping systems on long-term soil C-CO2 emissions and their relationship with carbon (C) mineralization of crop residues. A long-term experiment was conducted in a Red Oxisol in Cruz Alta, RS, Brazil, with subtropical climate Cfa (Köppen classification), mean annual precipitation of 1,774 mm and mean annual temperature of 19.2 ºC. Treatments consisted of two tillage systems: (a) conventional tillage (CT) and (b) no tillage (NT) in combination with three cropping systems: (a) R0- monoculture system (soybean/wheat), (b) R1- winter crop rotation (soybean/wheat/soybean/black oat), and (c) R2- intensive crop rotation (soybean/ black oat/soybean/black oat + common vetch/maize/oilseed radish/wheat). The soil C-CO2 efflux was measured every 14 days for two years (48 measurements), by trapping the CO2 in an alkaline solution. The soil gravimetric moisture in the 0-0.05 m layer was determined concomitantly with the C-CO2 efflux measurements. The crop residue C mineralization was evaluated with the mesh-bag method, with sampling 14, 28, 56, 84, 112, and 140 days after the beginning of the evaluation period for C measurements. Four C conservation indexes were used to assess the relation between C-CO2 efflux and soil C stock and its compartments. The crop residue C mineralization fit an exponential model in time. For black oat, wheat and maize residues, C mineralization was higher in CT than NT, while for soybean it was similar. Soil moisture was higher in NT than CT, mainly in the second year of evaluation. There was no difference in tillage systems for annual average C-CO2 emissions, but in some individual evaluations, differences between tillage systems were noticed for C-CO2 evolution. Soil C-CO2 effluxes followed a bi-modal pattern, with peaks in October/November and February/March. The highest emission was recorded in the summer and the lowest in the winter. The C-CO2 effluxes were weakly correlated to air temperature and not correlated to soil moisture. Based on the soil C conservation indexes investigated, NT associated to intensive crop rotation was more C conserving than CT with monoculture.

Delta(1232) isobar excitations in nuclear many-body systems from varios NN interactions.

Delta isobar components in the nuclear many-body wave function are investigated for the deuteron, light nuclei (16O), and infinite nuclear matter within the framework of the coupled-cluster theory. The predictions derived for various realistic models of the baryon-baryon interaction are compared to each other. These include local (V28) and nonlocal meson exchange potentials (Bonn2000) but also a model recently derived by the Salamanca group accounting for quark degrees of freedom. The characteristic differences which are obtained for the NDelta and Delta Delta correlation functions are related to the approximation made in deriving the matrix elements for the baryon-baryon interaction.

Higher-order colored-noise effects in multivariable systems

We extend the partial resummation technique of Fokker-Planck terms for multivariable stochastic differential equations with colored noise. As an example, a model system of a Brownian particle with colored noise is studied. We prove that the asymmetric behavior found in analog simulations is due to higher-order terms which are left out in that technique. On the contrary, the systematic ¿-expansion approach can explain the analog results.

Interfacial growth in driven diffusive systems

We have studied the growth of interfaces in driven diffusive systems well below the critical temperature by means of Monte Carlo simulations. We consider the region beyond the linear regime and of large values of the external field which has not been explored before. The simulations support the existence of interfacial traveling waves when asymmetry is introduced in the model, a result previously predicted by a linear-stability analysis. Furthermore, the generalization of the Gibbs-Thomson relation is discussed. The results provide evidence that the external field is a stabilizing effect which can be considered as effectively increasing the surface tension.

Numerical signs fos a transition in the 2d Random Field Ising Model at T=0

Intensive numerical studies of exact ground states of the two-dimensional ferromagnetic random field Ising model at T=0, with a Gaussian distribution of fields, are presented. Standard finite size scaling analysis of the data suggests the existence of a transition at ¿c=0.64±0.08. Results are compared with existing theories and with the study of metastable avalanches in the same model.

Dynamics of the two-dimensional gonihedric spin model

In this paper, we study dynamical aspects of the two-dimensional (2D) gonihedric spin model using both numerical and analytical methods. This spin model has vanishing microscopic surface tension and it actually describes an ensemble of loops living on a 2D surface. The self-avoidance of loops is parametrized by a parameter ¿. The ¿=0 model can be mapped to one of the six-vertex models discussed by Baxter, and it does not have critical behavior. We have found that allowing for ¿¿0 does not lead to critical behavior either. Finite-size effects are rather severe, and in order to understand these effects, a finite-volume calculation for non-self-avoiding loops is presented. This model, like his 3D counterpart, exhibits very slow dynamics, but a careful analysis of dynamical observables reveals nonglassy evolution (unlike its 3D counterpart). We find, also in this ¿=0 case, the law that governs the long-time, low-temperature evolution of the system, through a dual description in terms of defects. A power, rather than logarithmic, law for the approach to equilibrium has been found.