Studies of dynamics of physical agent ecosystems

This thesis addresses the problem of learning in physical heterogeneous multi-agent systems (MAS) and the analysis of the benefits of using heterogeneous MAS with respect to homogeneous ones. An algorithm is developed for this task; building on a previous work on stability in distributed systems by Tad Hogg and Bernardo Huberman, and combining two phenomena observed in natural systems, task partition and hierarchical dominance. This algorithm is devised for allowing agents to learn which are the best tasks to perform on the basis of each agent's skills and the contribution to the team global performance. Agents learn by interacting with the environment and other teammates, and get rewards from the result of the actions they perform. This algorithm is specially designed for problems where all robots have to co-operate and work simultaneously towards the same goal. One example of such a problem is role distribution in a team of heterogeneous robots that form a soccer team, where all members take decisions and co-operate simultaneously. Soccer offers the possibility of conducting research in MAS, where co-operation plays a very important role in a dynamical and changing environment. For these reasons and the experience of the University of Girona in this domain, soccer has been selected as the test-bed for this research. In the case of soccer, tasks are grouped by means of roles. One of the most interesting features of this algorithm is that it endows MAS with a high adaptability to changes in the environment. It allows the team to perform their tasks, while adapting to the environment. This is studied in several cases, for changes in the environment and in the robot's body. Other features are also analysed, especially a parameter that defines the fitness (biological concept) of each agent in the system, which contributes to performance and team adaptability. The algorithm is applied later to allow agents to learn in teams of homogeneous and heterogeneous robots which roles they have to select, in order to maximise team performance. The teams are compared and the performance is evaluated in the games against three hand-coded teams and against the different homogeneous and heterogeneous teams built in this thesis. This section focuses on the analysis of performance and task partition, in order to study the benefits of heterogeneity in physical MAS. In order to study heterogeneity from a rigorous point of view, a diversity measure is developed building on the hierarchic social entropy defined by Tucker Balch. This is adapted to quantify physical diversity in robot teams. This tool presents very interesting features, as it can be used in the future to design heterogeneous teams on the basis of the knowledge on other teams.

Finite-stretching corrections to the Milner-Witten-Cates theory for polymer brushes

This paper investigates finite-stretching corrections to the classical Milner-Witten-Cates theory for semi-dilute polymer brushes in a good solvent. The dominant correction to the free energy originates from an entropic repulsion caused by the impenetrability of the grafting surface, which produces a depletion of segments extending a distance $\mu \propto L^{-1}$ from the substrate, where $L$ is the classical brush height. The next most important correction is associated with the translational entropy of the chain ends, which creates the well-known tail where a small population of chains extend beyond the classical brush height by a distance $\xi \propto L^{-1/3}$. The validity of these corrections is confirmed by quantitative comparison with numerical self-consistent field theory.

Realizability of Point Processes

There are various situations in which it is natural to ask whether a given collection of k functions, ρ j (r 1,…,r j ), j=1,…,k, defined on a set X, are the first k correlation functions of a point process on X. Here we describe some necessary and sufficient conditions on the ρ j ’s for this to be true. Our primary examples are X=ℝ d , X=ℤ d , and X an arbitrary finite set. In particular, we extend a result by Ambartzumian and Sukiasian showing realizability at sufficiently small densities ρ 1(r). Typically if any realizing process exists there will be many (even an uncountable number); in this case we prove, when X is a finite set, the existence of a realizing Gibbs measure with k body potentials which maximizes the entropy among all realizing measures. We also investigate in detail a simple example in which a uniform density ρ and translation invariant ρ 2 are specified on ℤ; there is a gap between our best upper bound on possible values of ρ and the largest ρ for which realizability can be established.

Selection-mutation balance models with epistatic selection

We present an application of birth-and-death processes on configuration spaces to a generalized mutation4 selection balance model. The model describes the aging of population as a process of accumulation of mu5 tations in a genotype. A rigorous treatment demands that mutations correspond to points in abstract spaces. 6 Our model describes an infinite-population, infinite-sites model in continuum. The dynamical equation which 7 describes the system, is of Kimura-Maruyama type. The problem can be posed in terms of evolution of states 8 (differential equation) or, equivalently, represented in terms of Feynman-Kac formula. The questions of interest 9 are the existence of a solution, its asymptotic behavior, and properties of the limiting state. In the non-epistatic 10 case the problem was posed and solved in [Steinsaltz D., Evans S.N., Wachter K.W., Adv. Appl. Math., 2005, 11 35(1)]. In our model we consider a topological space X as the space of positions of mutations and the influence of epistatic potentials

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Fast training of self organizing maps for the visual exploration of molecular compounds

Visual exploration of scientific data in life science area is a growing research field due to the large amount of available data. The Kohonen’s Self Organizing Map (SOM) is a widely used tool for visualization of multidimensional data. In this paper we present a fast learning algorithm for SOMs that uses a simulated annealing method to adapt the learning parameters. The algorithm has been adopted in a data analysis framework for the generation of similarity maps. Such maps provide an effective tool for the visual exploration of large and multi-dimensional input spaces. The approach has been applied to data generated during the High Throughput Screening of molecular compounds; the generated maps allow a visual exploration of molecules with similar topological properties. The experimental analysis on real world data from the National Cancer Institute shows the speed up of the proposed SOM training process in comparison to a traditional approach. The resulting visual landscape groups molecules with similar chemical properties in densely connected regions.

High-resolution, unstructured meshes for hydrodynamic models of the Great Barrier Reef, Australia

Accuracy and mesh generation are key issues for the high-resolution hydrodynamic modelling of the whole Great Barrier Reef. Our objective is to generate suitable unstructured grids that can resolve topological and dynamical features like tidal jets and recirculation eddies in the wake of islands. A new strategy is suggested to refine the mesh in areas of interest taking into account the bathymetric field and an approximated distance to islands and reefs. Such a distance is obtained by solving an elliptic differential operator, with specific boundary conditions. Meshes produced illustrate both the validity and the efficiency of the adaptive strategy. Selection of refinement and geometrical parameters is discussed. (c) 2006 Elsevier Ltd. All rights reserved.

Hydrolysis of N2O5 on submicron sulfuric acid aerosols

The kinetics of reactive uptake of gaseous N2O5 on submicron sulfuric acid aerosol particles has been investigated using a laminar flow reactor coupled with a differential mobility analyzer (DMA) to characterize the aerosol. The particles were generated by homogeneous nucleation of SO3/H2O mixtures. In the H2SO4 concentration range 26.3−64.5 wt % the uptake coefficient was γ = 0.033 ± 0.004, independent of acid strength. For an acid strength of 45 wt % γ was found to decrease with increasing temperature over the range 263−298 K. From this, temperature dependence values of −115 ± 30 kJ/mol and −25.5 ± 8.4 J/K mol were determined for the changes in enthalpy and entropy of the uptake process, respectively. The results are consistent with a previous model of N2O5 hydrolysis involving both a direct and an acid catalyzed mechanism, with uptake under the experimental conditions limited by mass accommodation.

Effectiveness of landmark analysis for establishing locality in p2p networks

Locality to other nodes on a peer-to-peer overlay network can be established by means of a set of landmarks shared among the participating nodes. Each node independently collects a set of latency measures to landmark nodes, which are used as a multi-dimensional feature vector. Each peer node uses the feature vector to generate a unique scalar index which is correlated to its topological locality. A popular dimensionality reduction technique is the space filling Hilbert’s curve, as it possesses good locality preserving properties. However, there exists little comparison between Hilbert’s curve and other techniques for dimensionality reduction. This work carries out a quantitative analysis of their properties. Linear and non-linear techniques for scaling the landmark vectors to a single dimension are investigated. Hilbert’s curve, Sammon’s mapping and Principal Component Analysis have been used to generate a 1d space with locality preserving properties. This work provides empirical evidence to support the use of Hilbert’s curve in the context of locality preservation when generating peer identifiers by means of landmark vector analysis. A comparative analysis is carried out with an artificial 2d network model and with a realistic network topology model with a typical power-law distribution of node connectivity in the Internet. Nearest neighbour analysis confirms Hilbert’s curve to be very effective in both artificial and realistic network topologies. Nevertheless, the results in the realistic network model show that there is scope for improvements and better techniques to preserve locality information are required.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Kinetic and crystallographic studies of glucopyranose spirohydantoin and glucopyranosylamine analogs inhibitors of glycogen phosphorylase

Glycogen phosphorylase (GP) is currently exploited as a target for inhibition of hepatic glycogenolysis under high glucose conditions. Spirohydantoin of glucopyranose and N-acetyl-beta-D-glucopyranosylamine have been identified as the most potent inhibitors of GP that bind at the catalytic site. Four spirohydantoin and three beta-D-glucopyranosylamine analogs have been designed, synthesized and tested for inhibition of GP in kinetic experiments. Depending on the functional group introduced, the K(i) values varied from 16.5 microM to 1200 microM. In order to rationalize the kinetic results, we determined the crystal structures of the analogs in complex with GP. All the inhibitors bound at the catalytic site of the enzyme, by making direct and water-mediated hydrogen bonds with the protein and by inducing minor movements of the side chains of Asp283 and Asn284, of the 280s loop that blocks access of the substrate glycogen to the catalytic site, and changes in the water structure in the vicinity of the site. The differences observed in the Ki values of the analogs can be interpreted in terms of variations in hydrogen bonding and van der Waals interactions, desolvation effects, ligand conformational entropy, and displacement of water molecules on ligand binding to the catalytic site.

Synthesis of a catechol-based poly(ether ether ketone) ("o-PEEK") by classical step-growth polymerization and by entropically driven ring-opening polymerization of macrocyclic oligomers

An amorphous, catechol-based analogue of PEEK ("o-PEEK") has been prepared by a classical step-growth polymerization reaction between catechol and 4,4'-difluorobenzophenone and shown to be readily soluble in a range of organic solvents. Copolymers with p-PEEK have been investigated, including an amorphous 50: 50 composition and a semicrystalline though still organic-soluble material comprising 70% p-PEEK. o-PEEK has also been obtained by entropy-driven ring-opening polymerization of the macrocyclic oligomers (MCO's) formed by cyclo-condensation of catechol with 4,4'-difluorobenzophenone under pseudo-high-dilution conditions. The principal products of this latter reaction were the cyclic dimer 3a (20 wt %), cyclic trimer 3b (16%) cyclic tetramer 3c (14%), cyclic pentamer 3d (13%) and cyclic hexamer 3e (12%). Macrocycles 3a-c were isolated as pure compounds by gradient column chromatography, and the structures of the cyclic dimer 3a and cyclic tetramer 3c were analyzed by single-crystal X-ray diffraction. A mixture of MCO's, 3, of similar composition, was obtained by cyclodepolymerization of high molar mass o-PEEK in dilute soluion.

Proteomic profiling of neuromas reveals alterations in protein composition and local protein synthesis in hyper-excitable nerves

Neuropathic pain may arise following peripheral nerve injury though the molecular mechanisms associated with this are unclear. We used proteomic profiling to examine changes in protein expression associated with the formation of hyper-excitable neuromas derived from rodent saphenous nerves. A two-dimensional difference gel electrophoresis ( 2D-DIGE) profiling strategy was employed to examine protein expression changes between developing neuromas and normal nerves in whole tissue lysates. We found around 200 proteins which displayed a > 1.75-fold change in expression between neuroma and normal nerve and identified 55 of these proteins using mass spectrometry. We also used immunoblotting to examine the expression of low-abundance ion channels Nav1.3, Nav1.8 and calcium channel alpha 2 delta-1 subunit in this model, since they have previously been implicated in neuronal hyperexcitability associated with neuropathic pain. Finally, S(35)methionine in vitro labelling of neuroma and control samples was used to demonstrate local protein synthesis of neuron-specific genes. A number of cytoskeletal proteins, enzymes and proteins associated with oxidative stress were up-regulated in neuromas, whilst overall levels of voltage-gated ion channel proteins were unaffected. We conclude that altered mRNA levels reported in the somata of damaged DRG neurons do not necessarily reflect levels of altered proteins in hyper-excitable damaged nerve endings. An altered repertoire of protein expression, local protein synthesis and topological re-arrangements of ion channels may all play important roles in neuroma hyper-excitability.

Pressure-jump X-ray studies of liquid crystal transitions in lipids

In this paper, we give an overview of our studies by static and time-resolved X-ray diffraction of inverse cubic phases and phase transitions in lipids. In 1, we briefly discuss the lyotropic phase behaviour of lipids, focusing attention on non-lamellar structures, and their geometric/topological relationship to fusion processes in lipid membranes. Possible pathways for transitions between different cubic phases are also outlined. In 2, we discuss the effects of hydrostatic pressure on lipid membranes and lipid phase transitions, and describe how the parameters required to predict the pressure dependence of lipid phase transition temperatures can be conveniently measured. We review some earlier results of inverse bicontinuous cubic phases from our laboratory, showing effects such as pressure-induced formation and swelling. In 3, we describe the technique of pressure-jump synchrotron X-ray diffraction. We present results that have been obtained from the lipid system 1:2 dilauroylphosphatidylcholine/lauric acid for cubic-inverse hexagonal, cubic-cubic and lamellar-cubic transitions. The rate of transition was found to increase with the amplitude of the pressure-jump and with increasing temperature. Evidence for intermediate structures occurring transiently during the transitions was also obtained. In 4, we describe an IDL-based 'AXCESS' software package being developed in our laboratory to permit batch processing and analysis of the large X-ray datasets produced by pressure-jump synchrotron experiments. In 5, we present some recent results on the fluid lamellar-Pn3m cubic phase transition of the single-chain lipid 1-monoelaidin, which we have studied both by pressure-jump and temperature-jump X-ray diffraction. Finally, in 6, we give a few indicators of future directions of this research. We anticipate that the most useful technical advance will be the development of pressure-jump apparatus on the microsecond time-scale, which will involve the use of a stack of piezoelectric pressure actuators. The pressure-jump technique is not restricted to lipid phase transitions, but can be used to study a wide range of soft matter transitions, ranging from protein unfolding and DNA unwinding and transitions, to phase transitions in thermotropic liquid crystals, surfactants and block copolymers.

Thermal Fractionation and isothermal crystallization of polyethylene nanocomposites prepared by in situ polymerization

Nanocomposites of high-density polyethylene (HDPE) and carbon nanotubes (CNT) of different geometries (single wall, double wall, and multiwall; SWNT, DWNT, and MWNT) were prepared by in situ polymerization of ethylene on CNT whose surface had been previously treated with a metallocene catalytic system. In this work, we have studied the effects of applying the successive self-nucleation and annealing thermal fractionation technique (SSA) to the nanocomposites and have also determined the influence of composition and type of CNT on the isothermal crystallization behavior of the HDPE. SSA results indicate that all types of CNT induce the formation of a population of thicker lamellar crystals that melt at higher temperatures as compared to the crystals formed in neat HDPE prepared under the same catalytic and polymerization conditions and subjected to the same SSA treatment. Furthermore, the peculiar morphology induced by the CNT on the HDPE matrix allows the resolution of thermal fractionation to be much better. The isothermal crystallization results indicated that the strong nucleation effect caused by CNT reduced the supercooling needed for crystallization. The interaction between the HDPE chains and the surface of the CNT is probably very strong as judged by the results obtained, even though it is only physical in nature. When the total crystallinity achieved during isothermal crystallization is considered as a function of CNT content, it was found that a competition between nucleation and topological confinement could account for the results. At low CNT content the crystallinity increases (because of the nucleating effect of CNT on HDPE), however, at higher CNT content there is a dramatic reduction in crystallinity reflecting the increased confinement experienced by the HDPE chains at the interfaces which are extremely large in these nanocomposites. Another consequence of these strong interactions is the remarkable decrease in Avrami index as CNT content increases. When the Avrami index reduces to I or lower, nucleation dominates the overall kinetics as a consequence of confinement effects. Wide-angle X-ray experiments were performed at a high-energy synchrotron source and demonstrated that no change in the orthorhombic unit cell of HDPE occurred during crystallization with or without CNT.