Generalized mean field approximation for parallel dynamics of the Ising model

The dynamics of the non-equilibrium Ising model with parallel updates is investigated using a generalized mean field approximation that incorporates multiple two-site correlations at any two time steps, which can be obtained recursively. The proposed method shows significant improvement in predicting local system properties compared to other mean field approximation techniques, particularly in systems with symmetric interactions. Results are also evaluated against those obtained from Monte Carlo simulations. The method is also employed to obtain parameter values for the kinetic inverse Ising modeling problem, where couplings and local field values of a fully connected spin system are inferred from data. © 2014 IOP Publishing Ltd and SISSA Medialab srl.

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Dynamics of undeactuated cable-driven parallel robots

This thesis focuses on the dynamics of underactuated cable-driven parallel robots (UACDPRs), including various aspects of robotic theory and practice, such as workspace computation, parameter identification, and trajectory planning. After a brief introduction to CDPRs, UACDPR kinematic and dynamic models are analyzed, under the relevant assumption of inextensible cables. The free oscillatory motion of the end-effector (EE), which is a unique feature of underactuated mechanisms, is studied in detail, from both a kinematic and a dynamic perspective. The free (small) oscillations of the EE around equilibria are proved to be harmonic and the corresponding natural oscillation frequencies are analytically computed. UACDPR workspace computation and analysis are then performed. A new performance index is proposed for the analysis of the influence of actuator errors on cable tensions around equilibrium configurations, and a new type of workspace, called tension-error-insensitive, is defined as the set of poses that a UACDPR EE can statically attain even in presence of actuation errors, while preserving tensions between assigned (positive) bounds. EE free oscillations are then employed to conceive a novel procedure aimed at identifying the EE inertial parameters. This approach does not require the use of force or torque measurements. Moreover, a self-calibration procedure for the experimental determination of UACDPR initial cable lengths is developed, which consequently enables the robot to automatically infer the EE initial pose at machine start-up. Lastly, trajectory planning of UACDPRs is investigated. Two alternative methods are proposed, which aim at (i) reducing EE oscillations even when model parameters are uncertain or (ii) eliminate EE oscillations in case model parameters are perfectly known. EE oscillations are reduced in real-time by dynamically scaling a nominal trajectory and filtering it with an input shaper, whereas they can be eliminated if an off-line trajectory is computed that accounts for the system internal dynamics.

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Integrating structural and input design of a 2-DOF high-speed parallel manipulator: A flexible model-based approach

This paper discusses the integrated design of parallel manipulators, which exhibit varying dynamics. This characteristic affects the machine stability and performance. The design methodology consists of four main steps: (i) the system modeling using flexible multibody technique, (ii) the synthesis of reduced-order models suitable for control design, (iii) the systematic flexible model-based input signal design, and (iv) the evaluation of some possible machine designs. The novelty in this methodology is to take structural flexibilities into consideration during the input signal design; therefore, enhancing the standard design process which mainly considers rigid bodies dynamics. The potential of the proposed strategy is exploited for the design evaluation of a two degree-of-freedom high-speed parallel manipulator. The results are experimentally validated. (C) 2010 Elsevier Ltd. All rights reserved.

Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear H-1-N-15 NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta -strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. N-15 longitudinal and transverse relaxation rates, and {H-1}-N-15 heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85 +/- 0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides. (C) 2001 Academic Press.

Parallel simulation system for earthquake generation: Fault analysis modules and parallel coupling analysis

Solid earth simulations have recently been developed to address issues such as natural disasters, global environmental destruction and the conservation of natural resources. The simulation of solid earth phenomena involves the analysis of complex structures including strata, faults, and heterogeneous material properties. Simulation of the generation and cycle of earthquakes is particularly important, but such simulations require the analysis of complex fault dynamics. GeoFEM is a parallel finite-element analysis system intended for solid earth field phenomena problems. This paper describes recent development in the GeoFEM project for the simulation of earthquake generation and cycles.

Parallel changes in genetic diversity and species diversity following a natural disturbance.

We examined the spatial and temporal variation of species diversity and genetic diversity in a metacommunity comprising 16 species of freshwater gastropods. We monitored species abundance at five localities of the Ain river floodplain in southeastern France, over a period of four years. Using 190 AFLP loci, we monitored the genetic diversity of Radix balthica, one of the most abundant gastropod species of the metacommunity, twice during that period. An exceptionally intense drought occurred during the last two years and differentially affected the study sites. This allowed us to test the effect of natural disturbances on changes in both genetic and species diversity. Overall, local (alpha) diversity declined as reflected by lower values of gene diversity H(S) and evenness. In parallel, the among-sites (beta) diversity increased at both the genetic (F(ST)) and species (F(STC)) levels. These results suggest that disturbances can lead to similar changes in genetic and community structure through the combined effects of selective and neutral processes.

Roles of astrocytic sodium dynamics in the neurometabolic coupling : a study by cellular imaging

RESUME LARGE PUBLIC Le système nerveux central est principalement composé de deux types de cellules :les neurones et les cellules gliales. Ces dernières, bien que l'emportant en nombre sur les neurones, ont longtemps été considérées comme des cellules sans intérêts par les neuroscientifiques. Hors, les connaissances modernes à leurs sujets indiquent qu'elles participent à la plupart des tâches physiologiques du cerveau. Plus particulièrement, elles prennent part aux processus énergétiques cérébraux. Ceux-ci, en plus d'être vitaux, sont particulièrement intrigants puisque le cerveau représente seulement 2 % de la masse corporelle mais consomme environ 25 % du glucose (substrat énergétique) corporel. Les astrocytes, un type de cellules gliales, jouent un rôle primordial dans cette formidable utilisation de glucose par le cerveau. En effet, l'activité neuronale (transmission de l'influx nerveux) est accompagnée d'une augmentation de la capture de glucose, issu de la circulation sanguine, par les astrocytes. Ce phénomène est appelé le «couplage neurométabolique » entre neurones et astrocytes. L'ion sodium fait partie des mécanismes cellulaires entrant en fonction lors de ces processus. Ainsi, dans le cadre de cette thèse, les aspects dynamiques de la régulation du sodium astrocytaire et leurs implications dans le couplage neurométabolique ont été étudiés par des techniques d'imagerie cellulaires. Ces études ont démontré que les mitochondries, machineries cellulaires convertissant l'énergie contenue dans le glucose, participent à la régulation du sodium astrocytaire. De plus, ce travail de thèse a permis de découvrir que les astrocytes sont capables de se transmettre, sous forme de vagues de sodium se propageant de cellules en cellules, un message donnant l'ordre d'accroître leur consommation d'énergie. Cette voie de signalisation leur permettrait de fournir de l'énergie aux neurones suite à leur activation. RESUME Le glutamate libéré dans la fente synaptique pendant l'activité neuronale, est éliminé par les astrocytes environnants. Le glutamate est co-transporté avec des ions sodiques, induisant une augmentation intracellulaire de sodium (Na+i) dans les astrocytes. Cette élévation de Na+i déclenche une cascade de mécanismes moléculaires qui aboutissent à la production de substrats énergétiques pouvant être utilisés par les neurones. Durant cette thèse, la mesure simultanée du sodium mitochondrial (Na+mit) et cytosolique par des techniques d'imagerie utilisant des sondes fluorescentes spécifiques, a indiqué que les variations de Na+i induites par le transport du glutamate sont transmises aux mitochondries. De plus, les voies d'entrée et de sortie du sodium mitochondrial ont été identifiées. L'échangeur de Na+ et de Ca2+ mitochondrial semble jouer un rôle primordial dans l'influx de Na+mit, alors que l'efflux de Na+mit est pris en charge par l'échangeur de Na+ et de H+ mitochondrial. L'étude du Na+mit a nécessité l'utilisation d'un système de photoactivation. Les sources de lumière ultraviolette (UV) classiques utilisées à cet effet (lasers, lampes à flash) ayant plusieurs désavantages, une alternative efficace et peu coûteuse a été développée. Il s'agit d'un système compact utilisant une diode électroluminescente (LED) à haute puissance et de longueur d'onde de 365nm. En plus de leurs rôles dans le couplage neurométabolique, les astrocytes participent à la signalisation multicellulaire en transmettant des vagues intercellulaires de calcium. Ce travail de thèse démontre également que des vagues intercellulaires de sodium peuvent être évoquées en parallèle à ces vagues calciques. Le glutamate, suite à sa libération par un mécanisme dépendent du calcium, est réabsorbé par les transporteurs au glutamate. Ce mécanisme a pour conséquence la génération de vagues sodiques se propageant de cellules en cellules. De plus, ces vagues sodiques sont corrélées spatialement avec une consommation accrue de glucose par les astrocytes. En conclusion, ce travail de thèse a permis de montrer que le signal sodique astrocytaire, déclenché en réponse au glutamate, se propage à la fois de façon intracellulaire aux mitochondries et de façon intercellulaire. Ces résultats suggèrent que les astrocytes fonctionnent comme un réseau de cellules nécessaire au couplage énergétique concerté entre neurones et astrocytes et que le sodium est un élément clé dans les mécanismes de signalisations cellulaires sous-jacents. SUMMARY Glutamate, released in the synaptic cleft during neuronal activity, is removed by surrounding astrocytes. Glutamate is taken-up with Na+ ions by specific transporters, inducing an intracellular Na+ (Na+i) elevation in astrocytes which triggers a cascade of molecular mechanisms that provides metabolic substrates to neurons. Thus, astrocytic Na+i homeostasis represents a key component of the so-called neurometabolic coupling. In this context, the first part of this thesis work was aimed at investigating whether cytosolic Na+ changes are transmitted to mitochondria, which could therefore influence their function and contribute to the overall intracellular Na+ regulation. Simultaneous monitoring of both mitochondrial Na+ (Na+mit) and cytosolic Na+ changes with fluorescent dyes revealed that glutamate-evoked cytosolic Na+ elevations are indeed transmitted to mitochondria. The mitochondrial Na+/Ca2+ exchangers have a prominent role in the regulation of Na+mit influx pathway, and Na+mit extrusion appears to be mediated by Na+/H+ exchangers. To demonstrate the implication of Na+/Ca2+ exchangers, this study has required the technical development of an UV-flash photolysis system. Because light sources for flash photolysis have to be powerful and in the near UV range, the use of UV lasers or flash lamps is usually required. As an alternative to these UV sources that have several drawbaks, we developped a compact, efficient and lowcost flash photolysis system which employs a high power 365nm light emitting diode. In addition to their role in neurometabolic coupling, astrocytes participate in multicellular signaling by transmitting intercellular Ca2+ waves. The third part of this thesis show that intercellular Na+ waves can be evoked in parallel to Ca2+ waves. Glutamate released by a Ca2+ wave-dependent mechanism is taken up by glutamate transporters, resulting in a regenerative propagation of cytosolic Na+ increases. Na+ waves in turn lead to a spatially correlated increase in glucose uptake. In conclusion, the present thesis demonstrates that glutamate-induced Na+ changes occurring in the cytosol of astrocytes propagate to both the mitochondrial matrix and the astrocytic network. These results furthermore support the view that astrocytic Na+ is a signal coupled to the brain energy metabolism.

The Caribbean-North America-Cocos Triple Junction and the dynamics of the Polochic-Motagua fault systems: Pull-up and zipper models

The Polochic-Motagua fault systems (PMFS) are part of the sinistral transform boundary between the North American and Caribbean plates. To the west, these systems interact with the subduction zone of the Cocos plate, forming a subduction-subduction-transform triple junction. The North American plate moves westward relative to the Caribbean plate. This movement does not affect the geometry of the subducted Cocos plate, which implies that deformation is accommodated entirely in the two overriding plates. Structural data, fault kinematic analysis, and geomorphic observations provide new elements that help to understand the late Cenozoic evolution of this triple junction. In the Miocene, extension and shortening occurred south and north of the Motagua fault, respectively. This strain regime migrated northward to the Polochic fault after the late Miocene. This shift is interpreted as a ``pull-up'' of North American blocks into the Caribbean realm. To the west, the PMFS interact with a trench-parallel fault zone that links the Tonala fault to the Jalpatagua fault. These faults bound a fore-arc sliver that is shared by the two overriding plates. We propose that the dextral Jalpatagua fault merges with the sinistral PMFS, leaving behind a suturing structure, the Tonala fault. This tectonic ``zipper'' allows the migration of the triple junction. As a result, the fore-arc sliver comes into contact with the North American plate and helps to maintain a linear subduction zone along the trailing edge of the Caribbean plate. All these processes currently make the triple junction increasingly diffuse as it propagates eastward and inland within both overriding plates.

On the computation of reducible invariant tori on a parallel computer

We present an algorithm for the computation of reducible invariant tori of discrete dynamical systems that is suitable for tori of dimensions larger than 1. It is based on a quadratically convergent scheme that approximates, at the same time, the Fourier series of the torus, its Floquet transformation, and its Floquet matrix. The Floquet matrix describes the linearization of the dynamics around the torus and, hence, its linear stability. The algorithm presents a high degree of parallelism, and the computational effort grows linearly with the number of Fourier modes needed to represent the solution. For these reasons it is a very good option to compute quasi-periodic solutions with several basic frequencies. The paper includes some examples (flows) to show the efficiency of the method in a parallel computer. In these flows we compute invariant tori of dimensions up to 5, by taking suitable sections.

Co-evolutionary dynamics between public good producers and cheats in the bacterium Pseudomonas aeruginosa.

The production of beneficial public goods is common in the microbial world, and so is cheating - the exploitation of public goods by nonproducing mutants. Here, we examine co-evolutionary dynamics between cooperators and cheats and ask whether cooperators can evolve strategies to reduce the burden of exploitation, and whether cheats in turn can improve their exploitation abilities. We evolved cooperators of the bacterium Pseudomonas aeruginosa, producing the shareable iron-scavenging siderophore pyoverdine, together with cheats, defective in pyoverdine production but proficient in uptake. We found that cooperators managed to co-exist with cheats in 56% of all replicates over approximately 150 generations of experimental evolution. Growth and competition assays revealed that co-existence was fostered by a combination of general adaptions to the media and specific adaptions to the co-evolving opponent. Phenotypic screening and whole-genome resequencing of evolved clones confirmed this pattern, and suggest that cooperators became less exploitable by cheats because they significantly reduced their pyoverdine investment. Cheats, meanwhile, improved exploitation efficiency through mutations blocking the costly pyoverdine-signalling pathway. Moreover, cooperators and cheats evolved reduced motility, a pattern that likely represents adaptation to laboratory conditions, but at the same time also affects social interactions by reducing strain mixing and pyoverdine sharing. Overall, we observed parallel evolution, where co-existence of cooperators and cheats was enabled by a combination of adaptations to the abiotic and social environment and their interactions.