Self-organized criticality and synchronization in a coupled map lattice model of integrate-and-fire oscillators

We introduce two coupled map lattice models with nonconservative interactions and a continuous nonlinear driving. Depending on both the degree of conservation and the convexity of the driving we find different behaviors, ranging from self-organized criticality, in the sense that the distribution of events (avalanches) obeys a power law, to a macroscopic synchronization of the population of oscillators, with avalanches of the size of the system.

Symmetries and fixed point stability of stochastic differential equations modeling self-organized criticality

A stochastic nonlinear partial differential equation is constructed for two different models exhibiting self-organized criticality: the Bak-Tang-Wiesenfeld (BTW) sandpile model [Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 364 (1988)] and the Zhang model [Phys. Rev. Lett. 63, 470 (1989)]. The dynamic renormalization group (DRG) enables one to compute the critical exponents. However, the nontrivial stable fixed point of the DRG transformation is unreachable for the original parameters of the models. We introduce an alternative regularization of the step function involved in the threshold condition, which breaks the symmetry of the BTW model. Although the symmetry properties of the two models are different, it is shown that they both belong to the same universality class. In this case the DRG procedure leads to a symmetric behavior for both models, restoring the broken symmetry, and makes accessible the nontrivial fixed point. This technique could also be applied to other problems with threshold dynamics.

Noise and dynamics of self-organized critical phenomena

Different microscopic models exhibiting self-organized criticality are studied numerically and analytically. Numerical simulations are performed to compute critical exponents, mainly the dynamical exponent, and to check universality classes. We find that various models lead to the same exponent, but this universality class is sensitive to disorder. From the dynamic microscopic rules we obtain continuum equations with different sources of noise, which we call internal and external. Different correlations of the noise give rise to different critical behavior. A model for external noise is proposed that makes the upper critical dimensionality equal to 4 and leads to the possible existence of a phase transition above d=4. Limitations of the approach of these models by a simple nonlinear equation are discussed.

Self-dynamic structure factor of dense liquids: Theory and simulation

The self-intermediate dynamic structure factor Fs(k,t) of liquid lithium near the melting temperature is calculated by molecular dynamics. The results are compared with the predictions of several theoretical approaches, paying special attention to the Lovesey model and the Wahnstrm and Sjgren mode-coupling theory. To this end the results for the Fs(k,t) second memory function predicted by both models are compared with the ones calculated from the simulations.

Invaded cluster dynamics for frustrated models

The invaded cluster (IC) dynamics introduced by Machta et al. [Phys. Rev. Lett. 75, 2792 (1995)] is extended to the fully frustrated Ising model on a square lattice. The properties of the dynamics that exhibits numerical evidence of self-organized criticality are studied. The fluctuations in the IC dynamics are shown to be intrinsic of the algorithm and the fluctuation-dissipation theorem is no longer valid. The relaxation time is found to be very short and does not present a critical size dependence.

An evolutionary perspective on self-organized division of labor in social insects

Division of labor is a complex phenomenon observed throughout nature. Theoretical studies have focused either on its emergence through self-organization mechanisms or on its adaptive consequences. We suggest that the interaction of self-organization, which undoubtedly characterizes division of labor in social insects, and evolution should be further explored. We review the factors empirically shown to influence task choice. In light of these factors, we review the most important self-organization and evolutionary models for division of labor and outline their advantages and limitations. We describe ways to unify evolution and self-organization in the theoretical study of division of labor and recent results in this area. Finally, we discuss some benchmarks and primary challenges of this approach.

Detection, analysis and monitoring of slope movements by high-resolution digital elevation models

Summary Detection, analysis and monitoring of slope movements by high-resolution digital elevation modelsSlope movements, such as rockfalls, rockslides, shallow landslides or debris flows, are frequent in many mountainous areas. These natural hazards endanger the inhabitants and infrastructures making it necessary to assess the hazard and risk caused by these phenomena. This PhD thesis explores various approaches using digital elevation models (DEMs) - and particularly high-resolution DEMs created by aerial or terrestrial laser scanning (TLS) - that contribute to the assessment of slope movement hazard at regional and local scales.The regional detection of areas prone to rockfalls and large rockslides uses different morphologic criteria or geometric instability factors derived from DEMs, i.e. the steepness of the slope, the presence of discontinuities, which enable a sliding mechanism, and the denudation potential. The combination of these factors leads to a map of susceptibility to rockfall initiation that is in good agreement with field studies as shown with the example of the Little Mill Campground area (Utah, USA). Another case study in the Illgraben catchment in the Swiss Alps highlighted the link between areas with a high denudation potential and actual rockfall areas.Techniques for a detailed analysis and characterization of slope movements based on high-resolution DEMs have been developed for specific, localized sites, i.e. ancient slide scars, present active instabilities or potential slope instabilities. The analysis of the site's characteristics mainly focuses on rock slopes and includes structural analyses (orientation of discontinuities); estimation of spacing, persistence and roughness of discontinuities; failure mechanisms based on the structural setting; and volume calculations. For the volume estimation a new 3D approach was tested to reconstruct the topography before a landslide or to construct the basal failure surface of an active or potential instability. The rockslides at Åknes, Tafjord and Rundefjellet in western Norway were principally used as study sites to develop and test the different techniques.The monitoring of slope instabilities investigated in this PhD thesis is essentially based on multitemporal (or sequential) high-resolution DEMs, in particular sequential point clouds acquired by TLS. The changes in the topography due to slope movements can be detected and quantified by sequential TLS datasets, notably by shortest distance comparisons revealing the 3D slope movements over the entire region of interest. A detailed analysis of rock slope movements is based on the affine transformation between an initial and a final state of the rock mass and its decomposition into translational and rotational movements. Monitoring using TLS was very successful on the fast-moving Eiger rockslide in the Swiss Alps, but also on the active rockslides of Åknes and Nordnesfjellet (northern Norway). One of the main achievements on the Eiger and Aknes rockslides is to combine the site's morphology and structural setting with the measured slope movements to produce coherent instability models. Both case studies also highlighted a strong control of the structures in the rock mass on the sliding directions. TLS was also used to monitor slope movements in soils, such as landslides in sensitive clays in Québec (Canada), shallow landslides on river banks (Sorge River, Switzerland) and a debris flow channel (Illgraben).The PhD thesis underlines the broad uses of high-resolution DEMs and especially of TLS in the detection, analysis and monitoring of slope movements. Future studies should explore in more depth the different techniques and approaches developed and used in this PhD, improve them and better integrate the findings in current hazard assessment practices and in slope stability models.Résumé Détection, analyse et surveillance de mouvements de versant à l'aide de modèles numériques de terrain de haute résolutionDes mouvements de versant, tels que des chutes de blocs, glissements de terrain ou laves torrentielles, sont fréquents dans des régions montagneuses et mettent en danger les habitants et les infrastructures ce qui rend nécessaire d'évaluer le danger et le risque causé par ces phénomènes naturels. Ce travail de thèse explore diverses approches qui utilisent des modèles numériques de terrain (MNT) et surtout des MNT de haute résolution créés par scanner laser terrestre (SLT) ou aérien - et qui contribuent à l'évaluation du danger de mouvements de versant à l'échelle régionale et locale.La détection régionale de zones propices aux chutes de blocs ou aux éboulements utilise plusieurs critères morphologiques dérivés d'un MNT, tels que la pente, la présence de discontinuités qui permettent un mécanisme de glissement ou le potentiel de dénudation. La combinaison de ces facteurs d'instabilité mène vers une carte de susceptibilité aux chutes de blocs qui est en accord avec des travaux de terrain comme démontré avec l'exemple du Little Mill Campground (Utah, États-Unis). Un autre cas d'étude - l'Illgraben dans les Alpes valaisannes - a mis en évidence le lien entre les zones à fort potentiel de dénudation et les sources effectives de chutes de blocs et d'éboulements.Des techniques pour l'analyse et la caractérisation détaillée de mouvements de versant basées sur des MNT de haute résolution ont été développées pour des sites spécifiques et localisés, comme par exemple des cicatrices d'anciens éboulements et des instabilités actives ou potentielles. Cette analyse se focalise principalement sur des pentes rocheuses et comprend l'analyse structurale (orientation des discontinuités); l'estimation de l'espacement, la persistance et la rugosité des discontinuités; l'établissement des mécanismes de rupture; et le calcul de volumes. Pour cela une nouvelle approche a été testée en rétablissant la topographie antérieure au glissement ou en construisant la surface de rupture d'instabilités actuelles ou potentielles. Les glissements rocheux d'Åknes, Tafjord et Rundefjellet en Norvège ont été surtout utilisés comme cas d'étude pour développer et tester les diverses approches. La surveillance d'instabilités de versant effectuée dans cette thèse de doctorat est essentiellement basée sur des MNT de haute résolution multi-temporels (ou séquentiels), en particulier des nuages de points séquentiels acquis par SLT. Les changements topographiques dus aux mouvements de versant peuvent être détectés et quantifiés sur l'ensemble d'un glissement, notamment par comparaisons des distances les plus courtes entre deux nuages de points. L'analyse détaillée des mouvements est basée sur la transformation affine entre la position initiale et finale d'un bloc et sa décomposition en mouvements translationnels et rotationnels. La surveillance par SLT a démontré son potentiel avec l'effondrement d'un pan de l'Eiger dans les Alpes suisses, mais aussi aux glissements rocheux d'Aknes et Nordnesfjellet en Norvège. Une des principales avancées à l'Eiger et à Aknes est la création de modèles d'instabilité cohérents en combinant la morphologie et l'agencement structural des sites avec les mesures de déplacements. Ces deux cas d'étude ont aussi démontré le fort contrôle des structures existantes dans le massif rocheux sur les directions de glissement. Le SLT a également été utilisé pour surveiller des glissements dans des terrains meubles comme dans les argiles sensibles au Québec (Canada), sur les berges de la rivière Sorge en Suisse et dans le chenal à laves torrentielles de l'Illgraben.Cette thèse de doctorat souligne le vaste champ d'applications des MNT de haute résolution et particulièrement du SLT dans la détection, l'analyse et la surveillance des mouvements de versant. Des études futures devraient explorer plus en profondeur les différentes techniques et approches développées, les améliorer et mieux les intégrer dans des pratiques actuelles d'analyse de danger et surtout dans la modélisation de stabilité des versants.

Evolution of self-organized division of labor in a response threshold model.

Division of labor in social insects is determinant to their ecological success. Recent models emphasize that division of labor is an emergent property of the interactions among nestmates obeying to simple behavioral rules. However, the role of evolution in shaping these rules has been largely neglected. Here, we investigate a model that integrates the perspectives of self-organization and evolution. Our point of departure is the response threshold model, where we allow thresholds to evolve. We ask whether the thresholds will evolve to a state where division of labor emerges in a form that fits the needs of the colony. We find that division of labor can indeed evolve through the evolutionary branching of thresholds, leading to workers that differ in their tendency to take on a given task. However, the conditions under which division of labor evolves depend on the strength of selection on the two fitness components considered: amount of work performed and on worker distribution over tasks. When selection is strongest on the amount of work performed, division of labor evolves if switching tasks is costly. When selection is strongest on worker distribution, division of labor is less likely to evolve. Furthermore, we show that a biased distribution (like 3:1) of workers over tasks is not easily achievable by a threshold mechanism, even under strong selection. Contrary to expectation, multiple matings of colony foundresses impede the evolution of specialization. Overall, our model sheds light on the importance of considering the interaction between specific mechanisms and ecological requirements to better understand the evolutionary scenarios that lead to division of labor in complex systems. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00265-012-1343-2) contains supplementary material, which is available to authorized users.

The C-terminal domain of Plasmodium falciparum merozoite surface protein 3 self-assembles into alpha-helical coiled coil tetramer.

Proteins located on the surface of the pathogenic malaria parasite Plasmodium falciparum are objects of intensive studies due to their important role in the invasion of human cells and the accessibility to host antibodies thus making these proteins attractive vaccine candidates. One of these proteins, merozoite surface protein 3 (MSP3) represents a leading component among vaccine candidates; however, little is known about its structure and function. Our biophysical studies suggest that the 40 residue C-terminal domain of MSP3 protein self-assembles into a four-stranded alpha-helical coiled coil structure where alpha-helices are packed "side-by-side". A bioinformatics analysis provides an extended list of known and putative proteins from different species of Plasmodium which have such MSP3-like C-terminal domains. This finding allowed us to extend some conclusions of our studies to a larger group of the malaria surface proteins. Possible structural and functional roles of these highly conserved oligomerization domains in the intact merozoite surface proteins are discussed.

Dynamic models of viral replication and latency.

PURPOSE OF REVIEW: HIV targets primary CD4(+) T cells. The virus depends on the physiological state of its target cells for efficient replication, and, in turn, viral infection perturbs the cellular state significantly. Identifying the virus-host interactions that drive these dynamic changes is important for a better understanding of viral pathogenesis and persistence. The present review focuses on experimental and computational approaches to study the dynamics of viral replication and latency. RECENT FINDINGS: It was recently shown that only a fraction of the inducible latently infected reservoirs are successfully induced upon stimulation in ex-vivo models while additional rounds of stimulation make allowance for reactivation of more latently infected cells. This highlights the potential role of treatment duration and timing as important factors for successful reactivation of latently infected cells. The dynamics of HIV productive infection and latency have been investigated using transcriptome and proteome data. The cellular activation state has shown to be a major determinant of viral reactivation success. Mathematical models of latency have been used to explore the dynamics of the latent viral reservoir decay. SUMMARY: Timing is an important component of biological interactions. Temporal analyses covering aspects of viral life cycle are essential for gathering a comprehensive picture of HIV interaction with the host cell and untangling the complexity of latency. Understanding the dynamic changes tipping the balance between success and failure of HIV particle production might be key to eradicate the viral reservoir.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Bond strength of orthodontic brackets using different light and self-curing cements

The purpose of the study was to evaluate the shear bond strength of stainless steel orthodontic brackets directly bonded to extracted human premolar teeth. Fifty teeth were randomly divided into ¿ve groups: (1) System One (chemically cured composite resin), (2) Light Bond (light-cured composite resin), (3) Vivaglass Cem (self-curing glass ionomer cement), (4) Fuji Ortho LC (light-cured glass ionomer cement) used after 37% orthophosphoric acid¿etching of enamel (5) Fuji Ortho LC without orthophosphoric acid¿etching. The brackets were placed on the buccal and lingual surfaces of each tooth, and the specimens were stored in distilled water (24 hours) at 378C and thermocycled. Teeth were mounted on acrylic block frames, and brackets were debonded using an Instron machine. Shear bond strength values at fracture (Nw)were recorded. ANOVA and Student-Newman-Keuls multiple comparison tests were performed (P , .05). Bonding failure site was recorded by stereomicroscope and analyzed by Chi-square test, selected specimens of each group were observed by scanning electron microscope. System One attained the highest bond strength. Light Bond and Fuji Ortho LC, when using an acid-etching technique, obtained bond strengths that were within the range of estimated bond strength values for successful clinical bonding. Fuji Ortho LC and Vivaglass Cem left an almost clean enamel surface after debracketing.

A new self-expanding aortic stent valve with annular fixation: in vitro haemodynamic assessment.

OBJECTIVE: Balloon-expandable stent valves require flow reduction during implantation (rapid pacing). The present study was designed to compare a self-expanding stent valve with annular fixation versus a balloon-expandable stent valve. METHODS: Implantation of a new self-expanding stent valve with annular fixation (Symetis, Lausanne, Switzerland) was assessed versus balloon-expandable stent valve, in a modified Dynatek Dalta pulse duplicator (sealed port access to the ventricle for transapical route simulation), interfaced with a computer for digital readout, carrying a 25 mm porcine aortic valve. The cardiovascular simulator was programmed to mimic an elderly woman with aortic stenosis: 120/85 mmHg aortic pressure, 60 strokes/min (66.5 ml), 35% systole (2.8 l/min). RESULTS: A total of 450 cardiac cycles was analysed. Stepwise expansion of the self-expanding stent valve with annular fixation (balloon-expandable stent valve) resulted in systolic ventricular increase from 120 to 121 mmHg (126 to 830+/-76 mmHg)*, and left ventricular outflow obstruction with mean transvalvular gradient of 11+/-1.5 mmHg (366+/-202 mmHg)*, systolic aortic pressure dropped distal to the valve from 121 to 64.5+/-2 mmHg (123 to 55+/-30 mmHg) N.S., and output collapsed to 1.9+/-0.06 l/min (0.71+/-0.37 l/min* (before complete obstruction)). No valve migration occurred in either group. (*=p<0.05). CONCLUSIONS: Implantation of this new self-expanding stent valve with annular fixation has little impact on haemodynamics and has the potential for working heart implantation in vivo. Flow reduction (rapid pacing) is not necessary.

Referral from primary care to a physical activity programme: establishing long-term adherence? A randomized controlled trial. Rationale and study design

Background: Declining physical activity is associated with a rising burden of global disease. There is little evidence about effective ways to increase adherence to physical activity. Therefore, interventions are needed that produce sustained increases in adherence to physical activity and are cost-effective. The purpose is to assess the effectiveness of a primary care physical activity intervention in increasing adherence to physical activity in the general population seen in primary care. Method and design: Randomized controlled trial with systematic random sampling. A total of 424 subjects of both sexes will participate; all will be over the age of 18 with a low level of physical activity (according to the International Physical Activity Questionnaire, IPAQ), self-employed and from 9 Primary Healthcare Centres (PHC). They will volunteer to participate in a physical activity programme during 3 months (24 sessions; 2 sessions a week, 60 minutes per session). Participants from each PHC will be randomly allocated to an intervention (IG) and control group (CG). The following parameters will be assessed pre and post intervention in both groups: (1) health-related quality of life (SF-12), (2) physical activity stage of change (Prochaska's stages of change), (3) level of physical activity (IPAQ-short version), (4) change in perception of health (vignettes from the Cooperative World Organization of National Colleges, Academies, and Academic Associations of Family Physicians, COOP/WONCA), (5) level of social support for the physical activity practice (Social Support for Physical Activity Scale, SSPAS), and (6) control based on analysis (HDL, LDL and glycated haemoglobin).Participants' frequency of visits to the PHC will be registered over the six months before and after the programme. There will be a follow up in a face to face interview three, six and twelve months after the programme, with the reduced version of IPAQ, SF-12, SSPAS, and Prochaska's stages. Discussion: The pilot study showed the effectiveness of an enhanced low-cost, evidence-based intervention in increased physical activity and improved social support. If successful in demonstrating long-term improvements, this randomised controlled trial will be the first sustainable physical activity intervention based in primary care in our country to demonstrate longterm adherence to physical activity. Trial Registration: A service of the U.S. National Institutes of Health. Developed by the National Library of Medicine. ClinicalTrials.gov ID: NCT00714831.

What determines the wavelength of self-organized shoreline sand waves?

Shoreline undulations extending into the bathymetric contours with a length scale larger than that of the rhythmic surf zone bars are referred to as shoreline sand waves. Many observed undulations along sandy coasts display a wavelength in the order 1-7 km. Several models that are based on the hypothesis that sand waves emerge from a morphodynamic instability in case of very oblique wave incidence predict this range of wavelengths. Here we investigate the physical reasons for the wavelength selection and the main parametric trends of the wavelength in case of sand waves arising from such instability. It is shown that the existence of a minimum wavelength depends on an interplay between three factors affecting littoral drift: (A) the angle of wave fronts relative to local shoreline, which tends to cause maximum transport at the downdrift flank of the sand wave, (B) the refractive energy spreading which tends to cause maximum transport at the updrift flank and (C) wave focusing (de-focusing) by the capes (bays), which tends to cause maximum transport at the crest or slightly downdrift of it. Processes A and C cause decay of the sand waves while process B causes their growth. For low incidence angles, B is very weak so that a rectilinear shoreline is stable. For large angles and long sand waves, B is dominant and causes the growth of sand waves. For large angles and short sand waves C is dominant and the sand waves decay. Thus, wavelength selection depends on process C, which essentially depends on shoreline curvature. The growth rate of very long sand waves is weak because the alongshore gradients in sediment transport decrease with the wavelength. This is why there is an optimum or dominant wavelength. It is found that sand wave wavelength scales with λ0/β where λ0 is the water wave wavelength in deep water and β is the mean bed slope from shore to the wave base.