Optimal control strategy for a novel computer virus propagation model on scale-free networks

NASA Astrophysics Data System (ADS)

Zhang, Chunming; Huang, Haitao

2016-06-01

This paper aims to study the combined impact of reinstalling system and network topology on the spread of computer viruses over the Internet. Based on scale-free network, this paper proposes a novel computer viruses propagation model-SLBOSmodel. A systematic analysis of this new model shows that the virus-free equilibrium is globally asymptotically stable when its spreading threshold is less than one; nevertheless, it is proved that the viral equilibrium is permanent if the spreading threshold is greater than one. Then, the impacts of different model parameters on spreading threshold are analyzed. Next, an optimally controlled SLBOS epidemic model on complex networks is also studied. We prove that there is an optimal control existing for the control problem. Some numerical simulations are finally given to illustrate the main results.

A mixed SIR-SIS model to contain a virus spreading through networks with two degrees

NASA Astrophysics Data System (ADS)

Essouifi, Mohamed; Achahbar, Abdelfattah

Due to the fact that the “nodes” and “links” of real networks are heterogeneous, to model computer viruses prevalence throughout the Internet, we borrow the idea of the reduced scale free network which was introduced recently. The purpose of this paper is to extend the previous deterministic two subchains of Susceptible-Infected-Susceptible (SIS) model into a mixed Susceptible-Infected-Recovered and Susceptible-Infected-Susceptible (SIR-SIS) model to contain the computer virus spreading over networks with two degrees. Moreover, we develop its stochastic counterpart. Due to the high protection and security taken for hubs class, we suggest to treat it by using SIR epidemic model rather than the SIS one. The analytical study reveals that the proposed model admits a stable viral equilibrium. Thus, it is shown numerically that the mean dynamic behavior of the stochastic model is in agreement with the deterministic one. Unlike the infection densities i2 and i which both tend to a viral equilibrium for both approaches as in the previous study, i1 tends to the virus-free equilibrium. Furthermore, since a proportion of infectives are recovered, the global infection density i is minimized. Therefore, the permanent presence of viruses in the network due to the lower-degree nodes class. Many suggestions are put forward for containing viruses propagation and minimizing their damages.

Network-optimized congestion pricing : a parable, model and algorithm

DOT National Transportation Integrated Search

1995-05-31

This paper recites a parable, formulates a model and devises an algorithm for optimizing tolls on a road network. Such tolls induce an equilibrium traffic flow that is at once system-optimal and user-optimal. The parable introduces the network-wide c...

Fluctuation-dissipation relation and stationary distribution of an exactly solvable many-particle model for active biomatter far from equilibrium.

PubMed

Netz, Roland R

2018-05-14

An exactly solvable, Hamiltonian-based model of many massive particles that are coupled by harmonic potentials and driven by stochastic non-equilibrium forces is introduced. The stationary distribution and the fluctuation-dissipation relation are derived in closed form for the general non-equilibrium case. Deviations from equilibrium are on one hand characterized by the difference of the obtained stationary distribution from the Boltzmann distribution; this is possible because the model derives from a particle Hamiltonian. On the other hand, the difference between the obtained non-equilibrium fluctuation-dissipation relation and the standard equilibrium fluctuation-dissipation theorem allows us to quantify non-equilibrium in an alternative fashion. Both indicators of non-equilibrium behavior, i.e., deviations from the Boltzmann distribution and deviations from the equilibrium fluctuation-dissipation theorem, can be expressed in terms of a single non-equilibrium parameter α that involves the ratio of friction coefficients and random force strengths. The concept of a non-equilibrium effective temperature, which can be defined by the relation between fluctuations and the dissipation, is by comparison with the exactly derived stationary distribution shown not to hold, even if the effective temperature is made frequency dependent. The analysis is not confined to close-to-equilibrium situations but rather is exact and thus holds for arbitrarily large deviations from equilibrium. Also, the suggested harmonic model can be obtained from non-linear mechanical network systems by an expansion in terms of suitably chosen deviatory coordinates; the obtained results should thus be quite general. This is demonstrated by comparison of the derived non-equilibrium fluctuation dissipation relation with experimental data on actin networks that are driven out of equilibrium by energy-consuming protein motors. The comparison is excellent and allows us to extract the non-equilibrium parameter α from experimental spectral response and fluctuation data.

Fluctuation-dissipation relation and stationary distribution of an exactly solvable many-particle model for active biomatter far from equilibrium

NASA Astrophysics Data System (ADS)

Netz, Roland R.

2018-05-01

An exactly solvable, Hamiltonian-based model of many massive particles that are coupled by harmonic potentials and driven by stochastic non-equilibrium forces is introduced. The stationary distribution and the fluctuation-dissipation relation are derived in closed form for the general non-equilibrium case. Deviations from equilibrium are on one hand characterized by the difference of the obtained stationary distribution from the Boltzmann distribution; this is possible because the model derives from a particle Hamiltonian. On the other hand, the difference between the obtained non-equilibrium fluctuation-dissipation relation and the standard equilibrium fluctuation-dissipation theorem allows us to quantify non-equilibrium in an alternative fashion. Both indicators of non-equilibrium behavior, i.e., deviations from the Boltzmann distribution and deviations from the equilibrium fluctuation-dissipation theorem, can be expressed in terms of a single non-equilibrium parameter α that involves the ratio of friction coefficients and random force strengths. The concept of a non-equilibrium effective temperature, which can be defined by the relation between fluctuations and the dissipation, is by comparison with the exactly derived stationary distribution shown not to hold, even if the effective temperature is made frequency dependent. The analysis is not confined to close-to-equilibrium situations but rather is exact and thus holds for arbitrarily large deviations from equilibrium. Also, the suggested harmonic model can be obtained from non-linear mechanical network systems by an expansion in terms of suitably chosen deviatory coordinates; the obtained results should thus be quite general. This is demonstrated by comparison of the derived non-equilibrium fluctuation dissipation relation with experimental data on actin networks that are driven out of equilibrium by energy-consuming protein motors. The comparison is excellent and allows us to extract the non-equilibrium parameter α from experimental spectral response and fluctuation data.

A novel game theoretic approach for modeling competitive information diffusion in social networks with heterogeneous nodes

NASA Astrophysics Data System (ADS)

Agha Mohammad Ali Kermani, Mehrdad; Fatemi Ardestani, Seyed Farshad; Aliahmadi, Alireza; Barzinpour, Farnaz

2017-01-01

Influence maximization deals with identification of the most influential nodes in a social network given an influence model. In this paper, a game theoretic framework is developed that models a competitive influence maximization problem. A novel competitive influence model is additionally proposed that incorporates user heterogeneity, message content, and network structure. The proposed game-theoretic model is solved using Nash Equilibrium in a real-world dataset. It is shown that none of the well-known strategies are stable and at least one player has the incentive to deviate from the proposed strategy. Moreover, violation of Nash equilibrium strategy by each player leads to their reduced payoff. Contrary to previous works, our results demonstrate that graph topology, as well as the nodes' sociability and initial tendency measures have an effect on the determination of the influential node in the network.

A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions.

PubMed

Zhang, Fan; Yeh, Gour-Tsyh; Parker, Jack C; Brooks, Scott C; Pace, Molly N; Kim, Young-Jin; Jardine, Philip M; Watson, David B

2007-06-16

This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing N(E) equilibrium reactions and a set of reactive transport equations of M-N(E) kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.

Stochastic cycle selection in active flow networks

NASA Astrophysics Data System (ADS)

Woodhouse, Francis; Forrow, Aden; Fawcett, Joanna; Dunkel, Jorn

2016-11-01

Active biological flow networks pervade nature and span a wide range of scales, from arterial blood vessels and bronchial mucus transport in humans to bacterial flow through porous media or plasmodial shuttle streaming in slime molds. Despite their ubiquity, little is known about the self-organization principles that govern flow statistics in such non-equilibrium networks. By connecting concepts from lattice field theory, graph theory and transition rate theory, we show how topology controls dynamics in a generic model for actively driven flow on a network. Through theoretical and numerical analysis we identify symmetry-based rules to classify and predict the selection statistics of complex flow cycles from the network topology. Our conceptual framework is applicable to a broad class of biological and non-biological far-from-equilibrium networks, including actively controlled information flows, and establishes a new correspondence between active flow networks and generalized ice-type models.

Interdependent Network Recovery Games.

PubMed

Smith, Andrew M; González, Andrés D; Dueñas-Osorio, Leonardo; D'Souza, Raissa M

2017-10-30

Recovery of interdependent infrastructure networks in the presence of catastrophic failure is crucial to the economy and welfare of society. Recently, centralized methods have been developed to address optimal resource allocation in postdisaster recovery scenarios of interdependent infrastructure systems that minimize total cost. In real-world systems, however, multiple independent, possibly noncooperative, utility network controllers are responsible for making recovery decisions, resulting in suboptimal decentralized processes. With the goal of minimizing recovery cost, a best-case decentralized model allows controllers to develop a full recovery plan and negotiate until all parties are satisfied (an equilibrium is reached). Such a model is computationally intensive for planning and negotiating, and time is a crucial resource in postdisaster recovery scenarios. Furthermore, in this work, we prove this best-case decentralized negotiation process could continue indefinitely under certain conditions. Accounting for network controllers' urgency in repairing their system, we propose an ad hoc sequential game-theoretic model of interdependent infrastructure network recovery represented as a discrete time noncooperative game between network controllers that is guaranteed to converge to an equilibrium. We further reduce the computation time needed to find a solution by applying a best-response heuristic and prove bounds on ε-Nash equilibrium, where ε depends on problem inputs. We compare best-case and ad hoc models on an empirical interdependent infrastructure network in the presence of simulated earthquakes to demonstrate the extent of the tradeoff between optimality and computational efficiency. Our method provides a foundation for modeling sociotechnical systems in a way that mirrors restoration processes in practice. © 2017 Society for Risk Analysis.

Spreading dynamics of an e-commerce preferential information model on scale-free networks

NASA Astrophysics Data System (ADS)

Wan, Chen; Li, Tao; Guan, Zhi-Hong; Wang, Yuanmei; Liu, Xiongding

2017-02-01

In order to study the influence of the preferential degree and the heterogeneity of underlying networks on the spread of preferential e-commerce information, we propose a novel susceptible-infected-beneficial model based on scale-free networks. The spreading dynamics of the preferential information are analyzed in detail using the mean-field theory. We determine the basic reproductive number and equilibria. The theoretical analysis indicates that the basic reproductive number depends mainly on the preferential degree and the topology of the underlying networks. We prove the global stability of the information-elimination equilibrium. The permanence of preferential information and the global attractivity of the information-prevailing equilibrium are also studied in detail. Some numerical simulations are presented to verify the theoretical results.

Interplay of node connectivity and epidemic rates in the dynamics of epidemic networks

DOE PAGES

Kostova, Tanya

2010-07-09

We present and analyze a discrete-time susceptible-infected epidemic network model which represents each host as a separate entity and allows heterogeneous hosts and contacts. We establish a necessary and sufficient condition for global stability of the disease-free equilibrium of the system (defined as epidemic controllability) which defines the epidemic reproduction number of the network. When this condition is not fulfilled, we show that the system has a unique, locally stable equilibrium. As a result, we further derive sufficient conditions for epidemic controllability in terms of the epidemic rates and the network topology.

Global dynamics of a network-based SIQRS epidemic model with demographics and vaccination

NASA Astrophysics Data System (ADS)

Huang, Shouying; Chen, Fengde; Chen, Lijuan

2017-02-01

This paper investigates a new SIQRS epidemic model with demographics and vaccination on complex heterogeneous networks. We analytically derive the basic reproduction number R0, which determines not only the existence of endemic equilibrium but also the global dynamics of the model. The permanence of the disease and the globally asymptotical stability of disease-free equilibrium are proved in detail. By using a monotone iterative technique, we show that the unique endemic equilibrium is globally attractive under certain conditions. Our results really improve and enrich the results in Li et al (2014) [14]. Interestingly, the basic reproduction number R0 bears no relation to the degree-dependent birth, but our simulations indicate that the degree-dependent birth does affect the epidemic dynamics. Furthermore, we find that quarantine plays a more active role than vaccination in controlling the disease.

Conflict and convention in dynamic networks.

PubMed

Foley, Michael; Forber, Patrick; Smead, Rory; Riedl, Christoph

2018-03-01

An important way to resolve games of conflict (snowdrift, hawk-dove, chicken) involves adopting a convention: a correlated equilibrium that avoids any conflict between aggressive strategies. Dynamic networks allow individuals to resolve conflict via their network connections rather than changing their strategy. Exploring how behavioural strategies coevolve with social networks reveals new dynamics that can help explain the origins and robustness of conventions. Here, we model the emergence of conventions as correlated equilibria in dynamic networks. Our results show that networks have the tendency to break the symmetry between the two conventional solutions in a strongly biased way. Rather than the correlated equilibrium associated with ownership norms (play aggressive at home, not away), we usually see the opposite host-guest norm (play aggressive away, not at home) evolve on dynamic networks, a phenomenon common to human interaction. We also show that learning to avoid conflict can produce realistic network structures in a way different than preferential attachment models. © 2017 The Author(s).

Stability of an SAIRS alcoholism model on scale-free networks

NASA Astrophysics Data System (ADS)

Xiang, Hong; Liu, Ying-Ping; Huo, Hai-Feng

2017-05-01

A new SAIRS alcoholism model with birth and death on complex heterogeneous networks is proposed. The total population of our model is partitioned into four compartments: the susceptible individual, the light problem alcoholic, the heavy problem alcoholic and the recovered individual. The spread of alcoholism threshold R0 is calculated by the next generation matrix method. When R0 < 1, the alcohol free equilibrium is globally asymptotically stable, then the alcoholics will disappear. When R0 > 1, the alcoholism equilibrium is global attractivity, then the number of alcoholics will remain stable and alcoholism will become endemic. Furthermore, the modified SAIRS alcoholism model on weighted contact network is introduced. Dynamical behavior of the modified model is also studied. Numerical simulations are also presented to verify and extend theoretical results. Our results show that it is very important to treat alcoholics to control the spread of the alcoholism.

Identifying protein complexes in PPI network using non-cooperative sequential game.

PubMed

Maulik, Ujjwal; Basu, Srinka; Ray, Sumanta

2017-08-21

Identifying protein complexes from protein-protein interaction (PPI) network is an important and challenging task in computational biology as it helps in better understanding of cellular mechanisms in various organisms. In this paper we propose a noncooperative sequential game based model for protein complex detection from PPI network. The key hypothesis is that protein complex formation is driven by mechanism that eventually optimizes the number of interactions within the complex leading to dense subgraph. The hypothesis is drawn from the observed network property named small world. The proposed multi-player game model translates the hypothesis into the game strategies. The Nash equilibrium of the game corresponds to a network partition where each protein either belong to a complex or form a singleton cluster. We further propose an algorithm to find the Nash equilibrium of the sequential game. The exhaustive experiment on synthetic benchmark and real life yeast networks evaluates the structural as well as biological significance of the network partitions.

Knowledge Management through the Equilibrium Pattern Model for Learning

NASA Astrophysics Data System (ADS)

Sarirete, Akila; Noble, Elizabeth; Chikh, Azeddine

Contemporary students are characterized by having very applied learning styles and methods of acquiring knowledge. This behavior is consistent with the constructivist models where students are co-partners in the learning process. In the present work the authors developed a new model of learning based on the constructivist theory coupled with the cognitive development theory of Piaget. The model considers the level of learning based on several stages and the move from one stage to another requires learners' challenge. At each time a new concept is introduced creates a disequilibrium that needs to be worked out to return back to its equilibrium stage. This process of "disequilibrium/equilibrium" has been analyzed and validated using a course in computer networking as part of Cisco Networking Academy Program at Effat College, a women college in Saudi Arabia. The model provides a theoretical foundation for teaching especially in a complex knowledge domain such as engineering and can be used in a knowledge economy.

Dynamic Analysis of a Reaction-Diffusion Rumor Propagation Model

NASA Astrophysics Data System (ADS)

Zhao, Hongyong; Zhu, Linhe

2016-06-01

The rapid development of the Internet, especially the emergence of the social networks, leads rumor propagation into a new media era. Rumor propagation in social networks has brought new challenges to network security and social stability. This paper, based on partial differential equations (PDEs), proposes a new SIS rumor propagation model by considering the effect of the communication between the different rumor infected users on rumor propagation. The stabilities of a nonrumor equilibrium point and a rumor-spreading equilibrium point are discussed by linearization technique and the upper and lower solutions method, and the existence of a traveling wave solution is established by the cross-iteration scheme accompanied by the technique of upper and lower solutions and Schauder’s fixed point theorem. Furthermore, we add the time delay to rumor propagation and deduce the conditions of Hopf bifurcation and stability switches for the rumor-spreading equilibrium point by taking the time delay as the bifurcation parameter. Finally, numerical simulations are performed to illustrate the theoretical results.

Multi-equilibrium property of metabolic networks: SSI module.

PubMed

Lei, Hong-Bo; Zhang, Ji-Feng; Chen, Luonan

2011-06-20

Revealing the multi-equilibrium property of a metabolic network is a fundamental and important topic in systems biology. Due to the complexity of the metabolic network, it is generally a difficult task to study the problem as a whole from both analytical and numerical viewpoint. On the other hand, the structure-oriented modularization idea is a good choice to overcome such a difficulty, i.e. decomposing the network into several basic building blocks and then studying the whole network through investigating the dynamical characteristics of the basic building blocks and their interactions. Single substrate and single product with inhibition (SSI) metabolic module is one type of the basic building blocks of metabolic networks, and its multi-equilibrium property has important influence on that of the whole metabolic networks. In this paper, we describe what the SSI metabolic module is, characterize the rates of the metabolic reactions by Hill kinetics and give a unified model for SSI modules by using a set of nonlinear ordinary differential equations with multi-variables. Specifically, a sufficient and necessary condition is first given to describe the injectivity of a class of nonlinear systems, and then, the sufficient condition is used to study the multi-equilibrium property of SSI modules. As a main theoretical result, for the SSI modules in which each reaction has no more than one inhibitor, a sufficient condition is derived to rule out multiple equilibria, i.e. the Jacobian matrix of its rate function is nonsingular everywhere. In summary, we describe SSI modules and give a general modeling framework based on Hill kinetics, and provide a sufficient condition for ruling out multiple equilibria of a key type of SSI module.

Multi-equilibrium property of metabolic networks: SSI module

PubMed Central

2011-01-01

Background Revealing the multi-equilibrium property of a metabolic network is a fundamental and important topic in systems biology. Due to the complexity of the metabolic network, it is generally a difficult task to study the problem as a whole from both analytical and numerical viewpoint. On the other hand, the structure-oriented modularization idea is a good choice to overcome such a difficulty, i.e. decomposing the network into several basic building blocks and then studying the whole network through investigating the dynamical characteristics of the basic building blocks and their interactions. Single substrate and single product with inhibition (SSI) metabolic module is one type of the basic building blocks of metabolic networks, and its multi-equilibrium property has important influence on that of the whole metabolic networks. Results In this paper, we describe what the SSI metabolic module is, characterize the rates of the metabolic reactions by Hill kinetics and give a unified model for SSI modules by using a set of nonlinear ordinary differential equations with multi-variables. Specifically, a sufficient and necessary condition is first given to describe the injectivity of a class of nonlinear systems, and then, the sufficient condition is used to study the multi-equilibrium property of SSI modules. As a main theoretical result, for the SSI modules in which each reaction has no more than one inhibitor, a sufficient condition is derived to rule out multiple equilibria, i.e. the Jacobian matrix of its rate function is nonsingular everywhere. Conclusions In summary, we describe SSI modules and give a general modeling framework based on Hill kinetics, and provide a sufficient condition for ruling out multiple equilibria of a key type of SSI module. PMID:21689474

Dynamical behaviours and control measures of rumour-spreading model with consideration of network topology

NASA Astrophysics Data System (ADS)

Zhu, Linhe; Zhao, Hongyong

2017-07-01

A series of online rumours have seriously influenced the normal production and living of people. This paper aims to study the combined impact of psychological factor, propagation delay, network topology and control strategy on rumour diffusion over the online social networks. Based on an online social network, which is seen as a scale-free network, we model the spread of rumours by using a delayed SIS (Susceptible and Infected) epidemic-like model with consideration of psychological factor and network topology. First, through theoretical analysis, we illustrate the boundedness of the density of rumour-susceptible individuals and rumour-infected individuals. Second, we obtain the basic reproduction number R0 and prove the stability of the non-rumour equilibrium point and the rumour-spreading equilibrium point. Third, control strategies, such as uniform immunisation control, proportional immunisation control, targeted immunisation control and optimum control, are put forward to restrain rumour diffusion. Meanwhile, we have compared the differences of these control strategies. Finally, some representative numerical simulations are performed to verify the theoretical analysis results.

Determining the Impact of Personal Mobility Carbon Allowance Schemes in Transportation Networks

DOE PAGES

Aziz, H. M. Abdul; Ukkusuri, Satish V.; Zhan, Xianyuan

2016-10-17

We know that personal mobility carbon allowance (PMCA) schemes are designed to reduce carbon consumption from transportation networks. PMCA schemes influence the travel decision process of users and accordingly impact the system metrics including travel time and greenhouse gas (GHG) emissions. Here, we develop a multi-user class dynamic user equilibrium model to evaluate the transportation system performance when PMCA scheme is implemented. The results using Sioux-Falls test network indicate that PMCA schemes can achieve the emissions reduction goals for transportation networks. Further, users characterized by high value of travel time are found to be less sensitive to carbon budget inmore » the context of work trips. Results also show that PMCA scheme can lead to higher emissions for a path compared with the case without PMCA because of flow redistribution. The developed network equilibrium model allows us to examine the change in system states at different carbon allocation levels and to design parameters of PMCA schemes accounting for population heterogeneity.« less

Determining the Impact of Personal Mobility Carbon Allowance Schemes in Transportation Networks

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aziz, H. M. Abdul; Ukkusuri, Satish V.; Zhan, Xianyuan

We know that personal mobility carbon allowance (PMCA) schemes are designed to reduce carbon consumption from transportation networks. PMCA schemes influence the travel decision process of users and accordingly impact the system metrics including travel time and greenhouse gas (GHG) emissions. Here, we develop a multi-user class dynamic user equilibrium model to evaluate the transportation system performance when PMCA scheme is implemented. The results using Sioux-Falls test network indicate that PMCA schemes can achieve the emissions reduction goals for transportation networks. Further, users characterized by high value of travel time are found to be less sensitive to carbon budget inmore » the context of work trips. Results also show that PMCA scheme can lead to higher emissions for a path compared with the case without PMCA because of flow redistribution. The developed network equilibrium model allows us to examine the change in system states at different carbon allocation levels and to design parameters of PMCA schemes accounting for population heterogeneity.« less

Simplifying silicon burning: Application of quasi-equilibrium to (alpha) network nucleosynthesis

NASA Technical Reports Server (NTRS)

Hix, W. R.; Thielemann, F.-K.; Khokhlov, A. M.; Wheeler, J. C.

1997-01-01

While the need for accurate calculation of nucleosynthesis and the resulting rate of thermonuclear energy release within hydrodynamic models of stars and supernovae is clear, the computational expense of these nucleosynthesis calculations often force a compromise in accuracy to reduce the computational cost. To redress this trade-off of accuracy for speed, the authors present an improved nuclear network which takes advantage of quasi- equilibrium in order to reduce the number of independent nuclei, and hence the computational cost of nucleosynthesis, without significant reduction in accuracy. In this paper they will discuss the first application of this method, the further reduction in size of the minimal alpha network. The resultant QSE- reduced alpha network is twice as fast as the conventional alpha network it replaces and requires the tracking of half as many abundance variables, while accurately estimating the rate of energy generation. Such reduction in cost is particularly necessary for future generation of multi-dimensional models for supernovae.

Network-Theoretic Modeling of Fluid Flow

DTIC Science & Technology

2015-07-29

Final Report STIR: Network-Theoretic Modeling of Fluid Flow ARO Grant W911NF-14-1-0386 Program manager: Dr. Samuel Stanton ( August 1, 2014–April 30...Morzyński, M., and Comte , P., “A finite-time thermodynamics of unsteady fluid flows,” Journal of Non-Equilibrium Thermody- namics, Vol. 33, No. 2

Equilibria of perceptrons for simple contingency problems.

PubMed

Dawson, Michael R W; Dupuis, Brian

2012-08-01

The contingency between cues and outcomes is fundamentally important to theories of causal reasoning and to theories of associative learning. Researchers have computed the equilibria of Rescorla-Wagner models for a variety of contingency problems, and have used these equilibria to identify situations in which the Rescorla-Wagner model is consistent, or inconsistent, with normative models of contingency. Mathematical analyses that directly compare artificial neural networks to contingency theory have not been performed, because of the assumed equivalence between the Rescorla-Wagner learning rule and the delta rule training of artificial neural networks. However, recent results indicate that this equivalence is not as straightforward as typically assumed, suggesting a strong need for mathematical accounts of how networks deal with contingency problems. One such analysis is presented here, where it is proven that the structure of the equilibrium for a simple network trained on a basic contingency problem is quite different from the structure of the equilibrium for a Rescorla-Wagner model faced with the same problem. However, these structural differences lead to functionally equivalent behavior. The implications of this result for the relationships between associative learning, contingency theory, and connectionism are discussed.

An Adaptive QSE-reduced Nuclear Reaction Network for Silicon Burning

NASA Astrophysics Data System (ADS)

Parete-Koon, Suzanne; Hix, William Raphael; Thielemann, Friedrich-Karl

2010-02-01

The nuclei of the ``iron peak'' are formed late in the evolution of massive stars and during supernovae. Silicon burning during these events is responsible for the production of a wide range of nuclei with atomic mass numbers from 28 to 64. The large number of nuclei involved make accurate modeling of silicon burning computationally expensive. Examination of the physics of silicon burning reveals that the nuclear evolution is dominated by large groups of nuclei in mutual equilibrium. We present an improvement on our hybrid equilibrium-network scheme that takes advantage of this quasi-equilibrium (QSE) to reduce the number of independent variables calculated. Because the membership and number of these groups vary as the temperature, density and electron faction change, achieving maximal efficiency requires dynamic adjustment of group number and membership. The resultant QSE-reduced network is up to 20 times faster than the full network it replaces without significant loss of accuracy. These reductions in computational cost and the number of species evolved make QSE-reduced networks well suited for inclusion within hydrodynamic simulations, particularly in multi-dimensional applications. )

Symmetry, Hopf bifurcation, and the emergence of cluster solutions in time delayed neural networks.

PubMed

Wang, Zhen; Campbell, Sue Ann

2017-11-01

We consider the networks of N identical oscillators with time delayed, global circulant coupling, modeled by a system of delay differential equations with Z N symmetry. We first study the existence of Hopf bifurcations induced by the coupling time delay and then use symmetric Hopf bifurcation theory to determine how these bifurcations lead to different patterns of symmetric cluster oscillations. We apply our results to a case study: a network of FitzHugh-Nagumo neurons with diffusive coupling. For this model, we derive the asymptotic stability, global asymptotic stability, absolute instability, and stability switches of the equilibrium point in the plane of coupling time delay (τ) and excitability parameter (a). We investigate the patterns of cluster oscillations induced by the time delay and determine the direction and stability of the bifurcating periodic orbits by employing the multiple timescales method and normal form theory. We find that in the region where stability switching occurs, the dynamics of the system can be switched from the equilibrium point to any symmetric cluster oscillation, and back to equilibrium point as the time delay is increased.

Signal Propagation in Proteins and Relation to Equilibrium Fluctuations

PubMed Central

Chennubhotla, Chakra; Bahar, Ivet

2007-01-01

Elastic network (EN) models have been widely used in recent years for describing protein dynamics, based on the premise that the motions naturally accessible to native structures are relevant to biological function. We posit that equilibrium motions also determine communication mechanisms inherent to the network architecture. To this end, we explore the stochastics of a discrete-time, discrete-state Markov process of information transfer across the network of residues. We measure the communication abilities of residue pairs in terms of hit and commute times, i.e., the number of steps it takes on an average to send and receive signals. Functionally active residues are found to possess enhanced communication propensities, evidenced by their short hit times. Furthermore, secondary structural elements emerge as efficient mediators of communication. The present findings provide us with insights on the topological basis of communication in proteins and design principles for efficient signal transduction. While hit/commute times are information-theoretic concepts, a central contribution of this work is to rigorously show that they have physical origins directly relevant to the equilibrium fluctuations of residues predicted by EN models. PMID:17892319

A Simple and Accurate Network for Hydrogen and Carbon Chemistry in the Interstellar Medium

NASA Astrophysics Data System (ADS)

Gong, Munan; Ostriker, Eve C.; Wolfire, Mark G.

2017-07-01

Chemistry plays an important role in the interstellar medium (ISM), regulating the heating and cooling of the gas and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate into numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photodissociation region (PDR) code, and also with the Nelson & Langer (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network to 1D models, we find that the CO-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic-ray ionization rate in molecular clouds. We provide a simple fit for the locus of CO-dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the CO, {{CH}}x, and {{OH}}x abundances are consistent with equilibrium predictions for densities n=100{--}1000 {{cm}}-3, but the predicted equilibrium C abundance is higher than that seen in observations, signaling the potential importance of non-equilibrium/dynamical effects.

The QSE-Reduced Nuclear Reaction Network for Silicon Burning

NASA Astrophysics Data System (ADS)

Hix, W. Raphael; Parete-Koon, Suzanne T.; Freiburghaus, Christian; Thielemann, Friedrich-Karl

2007-09-01

Iron and neighboring nuclei are formed in massive stars shortly before core collapse and during their supernova outbursts, as well as during thermonuclear supernovae. Complete and incomplete silicon burning are responsible for the production of a wide range of nuclei with atomic mass numbers from 28 to 64. Because of the large number of nuclei involved, accurate modeling of silicon burning is computationally expensive. However, examination of the physics of silicon burning has revealed that the nuclear evolution is dominated by large groups of nuclei in mutual equilibrium. We present a new hybrid equilibrium-network scheme which takes advantage of this quasi-equilibrium in order to reduce the number of independent variables calculated. This allows accurate prediction of the nuclear abundance evolution, deleptonization, and energy generation at a greatly reduced computational cost when compared to a conventional nuclear reaction network. During silicon burning, the resultant QSE-reduced network is approximately an order of magnitude faster than the full network it replaces and requires the tracking of less than a third as many abundance variables, without significant loss of accuracy. These reductions in computational cost and the number of species evolved make QSE-reduced networks well suited for inclusion within hydrodynamic simulations, particularly in multidimensional applications.

The spreading dynamics of sexually transmitted diseases with birth and death on heterogeneous networks

NASA Astrophysics Data System (ADS)

Wang, Yi; Cao, Jinde; Alsaedi, Ahmed; Hayat, Tasawar

2017-02-01

In this paper, we formulate a deterministic model by including the vacant sites, which represent inactive individuals or potential contacts, to investigate the spreading dynamics of sexually transmitted diseases in heterogeneous networks. We first analytically derive the basic reproduction number R 0, which completely determines global dynamics of the system in the long run. Specifically, if R 0  <  1, the disease-free equilibrium is globally asymptotically stable, i.e. disease disappears from the network irrespective of initial infected numbers and distributions, whereas if R 0  >  1, the system is uniformly persistent around a unique endemic equilibrium, i.e. disease persists in the network. Furthermore, by using a suitable Lyapunov function the global stability of endemic equilibrium for low/high-risk infected individuals only is proved. Finally, the effects of three immunization schemes are studied and compared, and extensive numerical simulations are performed to investigate the effect of network topology and population turnover on disease spread. Our results suggest that population turnover could have great impact on the sexually transmitted disease system in heterogeneous networks, including the basic reproduction number and infection prevalence.

A complex network description on traditional Chinese medicine system

NASA Astrophysics Data System (ADS)

Sun, Anzheng; Zhang, Peipei; He, Yue; Su, Beibei; He, Da-Ren

2004-03-01

Chinese traditional philosophy believes that a healthy body can adjust itself to reach a dynamic equilibrium with the environment. At an ill state the equilibrium is lost. Any single medicine can only attack one problem and cannot recover the whole equilibrium. A prescription formulation (PF) usually contains an "emperor" or principal medicine, several "minister" or assistant medicines, some accessorial medicines, and one or two inducting or harmonizing edicines. Therefore different traditional Chinese medicine (TCM) appears in different number of PFs. The whole TCM system may be viewed as a network set composed of many complete graphs (PFs). The TCMs, which have the highest node degrees in the network, serve as the "bridges" between the complete graphs for forming the network. While the TCMs, which have lowest node degrees in the network, serve as the "emperors" in each complete graph. According to this idea we have performed a manual statistical investigation on approximately 1000 PFs and computed 8 different tatistical properties of the network. The results show that TCM system is a scale-free one and has a nice clustering structure. We are suggesting a dynamical model to describe the development of TCM system.

Tax Evasion and Nonequilibrium Model on Apollonian Networks

NASA Astrophysics Data System (ADS)

Lima, F. W. S.

2012-11-01

The Zaklan model had been proposed and studied recently using the equilibrium Ising model on square lattices (SLs) by [G. Zaklan, F. Westerhoff and D. Stauffer, J. Econ. Interact. Coord.4, 1 (2008), arXiv:0801.2980; G. Zaklan, F. W. S. Lima and F. Westerhoff, Physica A387, 5857 (2008)], near the critical temperature of the Ising model presenting a well-defined phase transition; but on normal and modified Apollonian networks (ANs), [J. S. Andrade, Jr., H. J. Herrmann, R. F. S. Andrade, and L. R. da Silva, Phys. Rev. Lett.94, 018702 (2005); R. F. S. Andrade, J. S. Andrade Jr. and H. J. Herrmann, Phys. Rev. E79, 036105 (2009)] studied the equilibrium Ising model. They showed the equilibrium Ising model not to present on ANs a phase transition of the type for the 2D Ising model. Here, using agent-based Monte Carlo simulations, we study the Zaklan model with the well-known majority-vote model (MVM) with noise and apply it to tax evasion on ANs, to show that differently from the Ising model the MVM on ANs presents a well-defined phase transition. To control the tax evasion in the economics model proposed by Zaklan et al., MVM is applied in the neighborhood of the critical noise qc to the Zaklan model. Here we show that the Zaklan model is robust because this can also be studied, besides using equilibrium dynamics of Ising model, through the nonequilibrium MVM and on various topologies giving the same behavior regardless of dynamic or topology used here.

Why are so many networks disassortative?

NASA Astrophysics Data System (ADS)

Johnson, Samuel; Torres, Joaquín J.; Marro, J.; Muñoz, Miguel A.

2011-03-01

A wide range of empirical networks—whether biological, technological, information-related or linguistic—generically exhibit important degree-degree anticorrelations (i.e., they are disassortative), the only exceptions being social ones, which tend to be positively correlated (assortative). Using an information-theory approach, we show that the equilibrium state of highly heterogeneous (scale-free) random networks is disassortative. This not only gives a parsimonious explanation to a long-standing question, but also provides a neutral model against which to compare experimental data and ascertain whether a given system is being driven from equilibrium by correlating mechanisms.

Modeling the propagation of mobile malware on complex networks

NASA Astrophysics Data System (ADS)

Liu, Wanping; Liu, Chao; Yang, Zheng; Liu, Xiaoyang; Zhang, Yihao; Wei, Zuxue

2016-08-01

In this paper, the spreading behavior of malware across mobile devices is addressed. By introducing complex networks to model mobile networks, which follows the power-law degree distribution, a novel epidemic model for mobile malware propagation is proposed. The spreading threshold that guarantees the dynamics of the model is calculated. Theoretically, the asymptotic stability of the malware-free equilibrium is confirmed when the threshold is below the unity, and the global stability is further proved under some sufficient conditions. The influences of different model parameters as well as the network topology on malware propagation are also analyzed. Our theoretical studies and numerical simulations show that networks with higher heterogeneity conduce to the diffusion of malware, and complex networks with lower power-law exponents benefit malware spreading.

Dynamic pricing of network goods with boundedly rational consumers.

PubMed

Radner, Roy; Radunskaya, Ami; Sundararajan, Arun

2014-01-07

We present a model of dynamic monopoly pricing for a good that displays network effects. In contrast with the standard notion of a rational-expectations equilibrium, we model consumers as boundedly rational and unable either to pay immediate attention to each price change or to make accurate forecasts of the adoption of the network good. Our analysis shows that the seller's optimal price trajectory has the following structure: The price is low when the user base is below a target level, is high when the user base is above the target, and is set to keep the user base stationary once the target level has been attained. We show that this pricing policy is robust to a number of extensions, which include the product's user base evolving over time and consumers basing their choices on a mixture of a myopic and a "stubborn" expectation of adoption. Our results differ significantly from those that would be predicted by a model based on rational-expectations equilibrium and are more consistent with the pricing of network goods observed in practice.

Dynamic pricing of network goods with boundedly rational consumers

PubMed Central

Radner, Roy; Radunskaya, Ami; Sundararajan, Arun

2014-01-01

We present a model of dynamic monopoly pricing for a good that displays network effects. In contrast with the standard notion of a rational-expectations equilibrium, we model consumers as boundedly rational and unable either to pay immediate attention to each price change or to make accurate forecasts of the adoption of the network good. Our analysis shows that the seller’s optimal price trajectory has the following structure: The price is low when the user base is below a target level, is high when the user base is above the target, and is set to keep the user base stationary once the target level has been attained. We show that this pricing policy is robust to a number of extensions, which include the product’s user base evolving over time and consumers basing their choices on a mixture of a myopic and a “stubborn” expectation of adoption. Our results differ significantly from those that would be predicted by a model based on rational-expectations equilibrium and are more consistent with the pricing of network goods observed in practice. PMID:24367101

Using Nonlinear Stochastic Evolutionary Game Strategy to Model an Evolutionary Biological Network of Organ Carcinogenesis Under a Natural Selection Scheme

PubMed Central

Chen, Bor-Sen; Tsai, Kun-Wei; Li, Cheng-Wei

2015-01-01

Molecular biologists have long recognized carcinogenesis as an evolutionary process that involves natural selection. Cancer is driven by the somatic evolution of cell lineages. In this study, the evolution of somatic cancer cell lineages during carcinogenesis was modeled as an equilibrium point (ie, phenotype of attractor) shifting, the process of a nonlinear stochastic evolutionary biological network. This process is subject to intrinsic random fluctuations because of somatic genetic and epigenetic variations, as well as extrinsic disturbances because of carcinogens and stressors. In order to maintain the normal function (ie, phenotype) of an evolutionary biological network subjected to random intrinsic fluctuations and extrinsic disturbances, a network robustness scheme that incorporates natural selection needs to be developed. This can be accomplished by selecting certain genetic and epigenetic variations to modify the network structure to attenuate intrinsic fluctuations efficiently and to resist extrinsic disturbances in order to maintain the phenotype of the evolutionary biological network at an equilibrium point (attractor). However, during carcinogenesis, the remaining (or neutral) genetic and epigenetic variations accumulate, and the extrinsic disturbances become too large to maintain the normal phenotype at the desired equilibrium point for the nonlinear evolutionary biological network. Thus, the network is shifted to a cancer phenotype at a new equilibrium point that begins a new evolutionary process. In this study, the natural selection scheme of an evolutionary biological network of carcinogenesis was derived from a robust negative feedback scheme based on the nonlinear stochastic Nash game strategy. The evolvability and phenotypic robustness criteria of the evolutionary cancer network were also estimated by solving a Hamilton–Jacobi inequality – constrained optimization problem. The simulation revealed that the phenotypic shift of the lung cancer-associated cell network takes 54.5 years from a normal state to stage I cancer, 1.5 years from stage I to stage II cancer, and 2.5 years from stage II to stage III cancer, with a reasonable match for the statistical result of the average age of lung cancer. These results suggest that a robust negative feedback scheme, based on a stochastic evolutionary game strategy, plays a critical role in an evolutionary biological network of carcinogenesis under a natural selection scheme. PMID:26244004

Dynamics analysis of SIR epidemic model with correlation coefficients and clustering coefficient in networks.

PubMed

Zhang, Juping; Yang, Chan; Jin, Zhen; Li, Jia

2018-07-14

In this paper, the correlation coefficients between nodes in states are used as dynamic variables, and we construct SIR epidemic dynamic models with correlation coefficients by using the pair approximation method in static networks and dynamic networks, respectively. Considering the clustering coefficient of the network, we analytically investigate the existence and the local asymptotic stability of each equilibrium of these models and derive threshold values for the prevalence of diseases. Additionally, we obtain two equivalent epidemic thresholds in dynamic networks, which are compared with the results of the mean field equations. Copyright © 2018 Elsevier Ltd. All rights reserved.

SEIR Model of Rumor Spreading in Online Social Network with Varying Total Population Size

NASA Astrophysics Data System (ADS)

Dong, Suyalatu; Deng, Yan-Bin; Huang, Yong-Chang

2017-10-01

Based on the infectious disease model with disease latency, this paper proposes a new model for the rumor spreading process in online social network. In this paper what we establish an SEIR rumor spreading model to describe the online social network with varying total number of users and user deactivation rate. We calculate the exact equilibrium points and reproduction number for this model. Furthermore, we perform the rumor spreading process in the online social network with increasing population size based on the original real world Facebook network. The simulation results indicate that the SEIR model of rumor spreading in online social network with changing total number of users can accurately reveal the inherent characteristics of rumor spreading process in online social network. Supported by National Natural Science Foundation of China under Grant Nos. 11275017 and 11173028

Agent-based spin model for financial markets on complex networks: Emergence of two-phase phenomena

NASA Astrophysics Data System (ADS)

Kim, Yup; Kim, Hong-Joo; Yook, Soon-Hyung

2008-09-01

We study a microscopic model for financial markets on complex networks, motivated by the dynamics of agents and their structure of interaction. The model consists of interacting agents (spins) with local ferromagnetic coupling and global antiferromagnetic coupling. In order to incorporate more realistic situations, we also introduce an external field which changes in time. From numerical simulations, we find that the model shows two-phase phenomena. When the local ferromagnetic interaction is balanced with the global antiferromagnetic interaction, the resulting return distribution satisfies a power law having a single peak at zero values of return, which corresponds to the market equilibrium phase. On the other hand, if local ferromagnetic interaction is dominant, then the return distribution becomes double peaked at nonzero values of return, which characterizes the out-of-equilibrium phase. On random networks, the crossover between two phases comes from the competition between two different interactions. However, on scale-free networks, not only the competition between the different interactions but also the heterogeneity of underlying topology causes the two-phase phenomena. Possible relationships between the critical phenomena of spin system and the two-phase phenomena are discussed.

Self-Organized Percolation and Critical Sales Fluctuations

NASA Astrophysics Data System (ADS)

Weisbuch, Gérard; Solomon, Sorin

There is a discrepancy between the standard view of equilibrium through price adjustment in economics and the observation of large fluctuations in stock markets. We study here a simple model where agents decisions not only depend upon their individual preferences but also upon information obtained from their neighbors in a social network. The model shows that information diffusion coupled to the adjustment process drives the system to criticality with large fluctuations rather than converging smoothly to equilibrium.

A Computational Method for Determining the Equilibrium Composition and Product Temperature in a LH2/LOX Combustor

NASA Technical Reports Server (NTRS)

Sozen, Mehmet

2003-01-01

In what follows, the model used for combustion of liquid hydrogen (LH2) with liquid oxygen (LOX) using chemical equilibrium assumption, and the novel computational method developed for determining the equilibrium composition and temperature of the combustion products by application of the first and second laws of thermodynamics will be described. The modular FORTRAN code developed as a subroutine that can be incorporated into any flow network code with little effort has been successfully implemented in GFSSP as the preliminary runs indicate. The code provides capability of modeling the heat transfer rate to the coolants for parametric analysis in system design.

Cascading Failures as Continuous Phase-Space Transitions

DOE PAGES

Yang, Yang; Motter, Adilson E.

2017-12-14

In network systems, a local perturbation can amplify as it propagates, potentially leading to a large-scale cascading failure. We derive a continuous model to advance our understanding of cascading failures in power-grid networks. The model accounts for both the failure of transmission lines and the desynchronization of power generators and incorporates the transient dynamics between successive steps of the cascade. In this framework, we show that a cascade event is a phase-space transition from an equilibrium state with high energy to an equilibrium state with lower energy, which can be suitably described in a closed form using a global Hamiltonian-likemore » function. From this function, we show that a perturbed system cannot always reach the equilibrium state predicted by quasi-steady-state cascade models, which would correspond to a reduced number of failures, and may instead undergo a larger cascade. We also show that, in the presence of two or more perturbations, the outcome depends strongly on the order and timing of the individual perturbations. These results offer new insights into the current understanding of cascading dynamics, with potential implications for control interventions.« less

Cascading Failures as Continuous Phase-Space Transitions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Yang; Motter, Adilson E.

In network systems, a local perturbation can amplify as it propagates, potentially leading to a large-scale cascading failure. We derive a continuous model to advance our understanding of cascading failures in power-grid networks. The model accounts for both the failure of transmission lines and the desynchronization of power generators and incorporates the transient dynamics between successive steps of the cascade. In this framework, we show that a cascade event is a phase-space transition from an equilibrium state with high energy to an equilibrium state with lower energy, which can be suitably described in a closed form using a global Hamiltonian-likemore » function. From this function, we show that a perturbed system cannot always reach the equilibrium state predicted by quasi-steady-state cascade models, which would correspond to a reduced number of failures, and may instead undergo a larger cascade. We also show that, in the presence of two or more perturbations, the outcome depends strongly on the order and timing of the individual perturbations. These results offer new insights into the current understanding of cascading dynamics, with potential implications for control interventions.« less

Game among interdependent networks: The impact of rationality on system robustness

NASA Astrophysics Data System (ADS)

Fan, Yuhang; Cao, Gongze; He, Shibo; Chen, Jiming; Sun, Youxian

2016-12-01

Many real-world systems are composed of interdependent networks that rely on one another. Such networks are typically designed and operated by different entities, who aim at maximizing their own payoffs. There exists a game among these entities when designing their own networks. In this paper, we study the game investigating how the rational behaviors of entities impact the system robustness. We first introduce a mathematical model to quantify the interacting payoffs among varying entities. Then we study the Nash equilibrium of the game and compare it with the optimal social welfare. We reveal that the cooperation among different entities can be reached to maximize the social welfare in continuous game only when the average degree of each network is constant. Therefore, the huge gap between Nash equilibrium and optimal social welfare generally exists. The rationality of entities makes the system inherently deficient and even renders it extremely vulnerable in some cases. We analyze our model for two concrete systems with continuous strategy space and discrete strategy space, respectively. Furthermore, we uncover some factors (such as weakening coupled strength of interdependent networks, designing a suitable topology dependence of the system) that help reduce the gap and the system vulnerability.

Network Aggregation in Transportation Planning : Volume II : A Fixed Point Method for Treating Traffic Equilibria

DOT National Transportation Integrated Search

1978-04-01

Volume 2 defines a new algorithm for the network equilibrium model that works in the space of path flows and is based on the theory of fixed point method. The goals of the study were broadly defined as the identification of aggregation practices and ...

Equilibrium & Nonequilibrium Fluctuation Effects in Biopolymer Networks

NASA Astrophysics Data System (ADS)

Kachan, Devin Michael

Fluctuation-induced interactions are an important organizing principle in a variety of soft matter systems. In this dissertation, I explore the role of both thermal and active fluctuations within cross-linked polymer networks. The systems I study are in large part inspired by the amazing physics found within the cytoskeleton of eukaryotic cells. I first predict and verify the existence of a thermal Casimir force between cross-linkers bound to a semi-flexible polymer. The calculation is complicated by the appearance of second order derivatives in the bending Hamiltonian for such polymers, which requires a careful evaluation of the the path integral formulation of the partition function in order to arrive at the physically correct continuum limit and properly address ultraviolet divergences. I find that cross linkers interact along a filament with an attractive logarithmic potential proportional to thermal energy. The proportionality constant depends on whether and how the cross linkers constrain the relative angle between the two filaments to which they are bound. The interaction has important implications for the synthesis of biopolymer bundles within cells. I model the cross-linkers as existing in two phases: bound to the bundle and free in solution. When the cross-linkers are bound, they behave as a one-dimensional gas of particles interacting with the Casimir force, while the free phase is a simple ideal gas. Demanding equilibrium between the two phases, I find a discontinuous transition between a sparsely and a densely bound bundle. This discontinuous condensation transition induced by the long-ranged nature of the Casimir interaction allows for a similarly abrupt structural transition in semiflexible filament networks between a low cross linker density isotropic phase and a higher cross link density bundle network. This work is supported by the results of finite element Brownian dynamics simulations of semiflexible filaments and transient cross-linkers. I speculate that cells take advantage of this equilibrium effect by tuning near the transition point, where small changes in free cross-linker density will affect large structural rearrangements between free filament networks and networks of bundles. Cells are naturally found far from equilibrium, where the active influx of energy from ATP consumption controls the dynamics. Motor proteins actively generate forces within biopolymer networks, and one may ask how these differ from the random stresses characteristic of equilibrium fluctuations. Besides the trivial observation that the magnitude is independent of temperature, I find that the processive nature of the motors creates a temporally correlated, or colored, noise spectrum. I model the network with a nonlinear scalar elastic theory in the presence of active driving, and study the long distance and large scale properties of the system with renormalization group techniques. I find that there is a new critical point associated with diverging correlation time, and that the colored noise produces novel frequency dependence in the renormalized transport coefficients. Finally, I study marginally elastic solids which have vanishing shear modulus due to the presence of soft modes, modes with zero deformation cost. Although network coordination is a useful metric for determining the mechanical response of random spring networks in mechanical equilibrium, it is insufficient for describing networks under external stress. In particular, under-constrained networks which are fluid-like at zero load will dynamically stiffen at a critical strain, as observed in numerical simulations and experimentally in many biopolymer networks. Drawing upon analogies to the stress induced unjamming of emulsions, I develop a kinetic theory to explain the rigidity transition in spring and filament networks. Describing the dynamic evolution of non-affine deformation via a simple mechanistic picture, I recover the emergent nonlinear strain-stiffening behavior and compare this behavior to the yield stress flow seen in soft glassy fluids. I extend this theory to account for coordination number inhomogeneities and predict a breakdown of universal scaling near the critical point at sufficiently high disorder, and discuss the utility for this type of model in describing biopolymer networks.

A Study on Standard Competition with Network Effect Based on Evolutionary Game Model

NASA Astrophysics Data System (ADS)

Wang, Ye; Wang, Bingdong; Li, Kangning

Owing to networks widespread in modern society, standard competition with network effect is now endowed with new connotation. This paper aims to study the impact of network effect on standard competition; it is organized in the mode of "introduction-model setup-equilibrium analysis-conclusion". Starting from a well-structured model of evolutionary game, it is then extended to a dynamic analysis. This article proves both theoretically and empirically that whether or not a standard can lead the market trends depends on the utility it would bring, and the author also discusses some advisable strategies revolving around the two factors of initial position and border break.

Perspective: Maximum caliber is a general variational principle for dynamical systems

NASA Astrophysics Data System (ADS)

Dixit, Purushottam D.; Wagoner, Jason; Weistuch, Corey; Pressé, Steve; Ghosh, Kingshuk; Dill, Ken A.

2018-01-01

We review here Maximum Caliber (Max Cal), a general variational principle for inferring distributions of paths in dynamical processes and networks. Max Cal is to dynamical trajectories what the principle of maximum entropy is to equilibrium states or stationary populations. In Max Cal, you maximize a path entropy over all possible pathways, subject to dynamical constraints, in order to predict relative path weights. Many well-known relationships of non-equilibrium statistical physics—such as the Green-Kubo fluctuation-dissipation relations, Onsager's reciprocal relations, and Prigogine's minimum entropy production—are limited to near-equilibrium processes. Max Cal is more general. While it can readily derive these results under those limits, Max Cal is also applicable far from equilibrium. We give examples of Max Cal as a method of inference about trajectory distributions from limited data, finding reaction coordinates in bio-molecular simulations, and modeling the complex dynamics of non-thermal systems such as gene regulatory networks or the collective firing of neurons. We also survey its basis in principle and some limitations.

Perspective: Maximum caliber is a general variational principle for dynamical systems.

PubMed

Dixit, Purushottam D; Wagoner, Jason; Weistuch, Corey; Pressé, Steve; Ghosh, Kingshuk; Dill, Ken A

2018-01-07

We review here Maximum Caliber (Max Cal), a general variational principle for inferring distributions of paths in dynamical processes and networks. Max Cal is to dynamical trajectories what the principle of maximum entropy is to equilibrium states or stationary populations. In Max Cal, you maximize a path entropy over all possible pathways, subject to dynamical constraints, in order to predict relative path weights. Many well-known relationships of non-equilibrium statistical physics-such as the Green-Kubo fluctuation-dissipation relations, Onsager's reciprocal relations, and Prigogine's minimum entropy production-are limited to near-equilibrium processes. Max Cal is more general. While it can readily derive these results under those limits, Max Cal is also applicable far from equilibrium. We give examples of Max Cal as a method of inference about trajectory distributions from limited data, finding reaction coordinates in bio-molecular simulations, and modeling the complex dynamics of non-thermal systems such as gene regulatory networks or the collective firing of neurons. We also survey its basis in principle and some limitations.

Game-Theoretic Models of Information Overload in Social Networks

NASA Astrophysics Data System (ADS)

Borgs, Christian; Chayes, Jennifer; Karrer, Brian; Meeder, Brendan; Ravi, R.; Reagans, Ray; Sayedi, Amin

We study the effect of information overload on user engagement in an asymmetric social network like Twitter. We introduce simple game-theoretic models that capture rate competition between celebrities producing updates in such networks where users non-strategically choose a subset of celebrities to follow based on the utility derived from high quality updates as well as disutility derived from having to wade through too many updates. Our two variants model the two behaviors of users dropping some potential connections (followership model) or leaving the network altogether (engagement model). We show that under a simple formulation of celebrity rate competition, there is no pure strategy Nash equilibrium under the first model. We then identify special cases in both models when pure rate equilibria exist for the celebrities: For the followership model, we show existence of a pure rate equilibrium when there is a global ranking of the celebrities in terms of the quality of their updates to users. This result also generalizes to the case when there is a partial order consistent with all the linear orders of the celebrities based on their qualities to the users. Furthermore, these equilibria can be computed in polynomial time. For the engagement model, pure rate equilibria exist when all users are interested in the same number of celebrities, or when they are interested in at most two. Finally, we also give a finite though inefficient procedure to determine if pure equilibria exist in the general case of the followership model.

Bringing metabolic networks to life: convenience rate law and thermodynamic constraints

PubMed Central

Liebermeister, Wolfram; Klipp, Edda

2006-01-01

Background Translating a known metabolic network into a dynamic model requires rate laws for all chemical reactions. The mathematical expressions depend on the underlying enzymatic mechanism; they can become quite involved and may contain a large number of parameters. Rate laws and enzyme parameters are still unknown for most enzymes. Results We introduce a simple and general rate law called "convenience kinetics". It can be derived from a simple random-order enzyme mechanism. Thermodynamic laws can impose dependencies on the kinetic parameters. Hence, to facilitate model fitting and parameter optimisation for large networks, we introduce thermodynamically independent system parameters: their values can be varied independently, without violating thermodynamical constraints. We achieve this by expressing the equilibrium constants either by Gibbs free energies of formation or by a set of independent equilibrium constants. The remaining system parameters are mean turnover rates, generalised Michaelis-Menten constants, and constants for inhibition and activation. All parameters correspond to molecular energies, for instance, binding energies between reactants and enzyme. Conclusion Convenience kinetics can be used to translate a biochemical network – manually or automatically - into a dynamical model with plausible biological properties. It implements enzyme saturation and regulation by activators and inhibitors, covers all possible reaction stoichiometries, and can be specified by a small number of parameters. Its mathematical form makes it especially suitable for parameter estimation and optimisation. Parameter estimates can be easily computed from a least-squares fit to Michaelis-Menten values, turnover rates, equilibrium constants, and other quantities that are routinely measured in enzyme assays and stored in kinetic databases. PMID:17173669

Spreading dynamics of a SIQRS epidemic model on scale-free networks

NASA Astrophysics Data System (ADS)

Li, Tao; Wang, Yuanmei; Guan, Zhi-Hong

2014-03-01

In order to investigate the influence of heterogeneity of the underlying networks and quarantine strategy on epidemic spreading, a SIQRS epidemic model on the scale-free networks is presented. Using the mean field theory the spreading dynamics of the virus is analyzed. The spreading critical threshold and equilibria are derived. Theoretical results indicate that the critical threshold value is significantly dependent on the topology of the underlying networks and quarantine rate. The existence of equilibria is determined by threshold value. The stability of disease-free equilibrium and the permanence of the disease are proved. Numerical simulations confirmed the analytical results.

The analysis of HIV/AIDS drug-resistant on networks

NASA Astrophysics Data System (ADS)

Liu, Maoxing

2014-01-01

In this paper, we present an Human Immunodeficiency Virus (HIV)/Acquired Immune Deficiency Syndrome (AIDS) drug-resistant model using an ordinary differential equation (ODE) model on scale-free networks. We derive the threshold for the epidemic to be zero in infinite scale-free network. We also prove the stability of disease-free equilibrium (DFE) and persistence of HIV/AIDS infection. The effects of two immunization schemes, including proportional scheme and targeted vaccination, are studied and compared. We find that targeted strategy compare favorably to a proportional condom using has prominent effect to control HIV/AIDS spread on scale-free networks.

The modeling and analysis of the word-of-mouth marketing

NASA Astrophysics Data System (ADS)

Li, Pengdeng; Yang, Xiaofan; Yang, Lu-Xing; Xiong, Qingyu; Wu, Yingbo; Tang, Yuan Yan

2018-03-01

As compared to the traditional advertising, word-of-mouth (WOM) communications have striking advantages such as significantly lower cost and much faster propagation, and this is especially the case with the popularity of online social networks. This paper focuses on the modeling and analysis of the WOM marketing. A dynamic model, known as the SIPNS model, capturing the WOM marketing processes with both positive and negative comments is established. On this basis, a measure of the overall profit of a WOM marketing campaign is proposed. The SIPNS model is shown to admit a unique equilibrium, and the equilibrium is determined. The impact of different factors on the equilibrium of the SIPNS model is illuminated through theoretical analysis. Extensive experimental results suggest that the equilibrium is much likely to be globally attracting. Finally, the influence of different factors on the expected overall profit of a WOM marketing campaign is ascertained both theoretically and experimentally. Thereby, some promotion strategies are recommended. To our knowledge, this is the first time the WOM marketing is treated in this way.

Communication: Introducing prescribed biases in out-of-equilibrium Markov models

NASA Astrophysics Data System (ADS)

Dixit, Purushottam D.

2018-03-01

Markov models are often used in modeling complex out-of-equilibrium chemical and biochemical systems. However, many times their predictions do not agree with experiments. We need a systematic framework to update existing Markov models to make them consistent with constraints that are derived from experiments. Here, we present a framework based on the principle of maximum relative path entropy (minimum Kullback-Leibler divergence) to update Markov models using stationary state and dynamical trajectory-based constraints. We illustrate the framework using a biochemical model network of growth factor-based signaling. We also show how to find the closest detailed balanced Markov model to a given Markov model. Further applications and generalizations are discussed.

How time delay and network design shape response patterns in biochemical negative feedback systems.

PubMed

Börsch, Anastasiya; Schaber, Jörg

2016-08-24

Negative feedback in combination with time delay can bring about both sustained oscillations and adaptive behaviour in cellular networks. Here, we study which design features of systems with delayed negative feedback shape characteristic response patterns with special emphasis on the role of time delay. To this end, we analyse generic two-dimensional delay differential equations describing the dynamics of biochemical signal-response networks. We investigate the influence of several design features on the stability of the model equilibrium, i.e., presence of auto-inhibition and/or mass conservation and the kind and/or strength of the delayed negative feedback. We show that auto-inhibition and mass conservation have a stabilizing effect, whereas increasing abruptness and decreasing feedback threshold have a de-stabilizing effect on the model equilibrium. Moreover, applying our theoretical analysis to the mammalian p53 system we show that an auto-inhibitory feedback can decouple period and amplitude of an oscillatory response, whereas the delayed feedback can not. Our theoretical framework provides insight into how time delay and design features of biochemical networks act together to elicit specific characteristic response patterns. Such insight is useful for constructing synthetic networks and controlling their behaviour in response to external stimulation.

An experimental study of the Online Information Paradox: Does en-route information improve road network performance?

PubMed

Wijayaratna, Kasun P; Dixit, Vinayak V; Denant-Boemont, Laurent; Waller, S Travis

2017-01-01

This study investigates the empirical presence of a theoretical transportation paradox, defined as the "Online Information Paradox" (OIP). The paradox suggests that, for certain road networks, the provision of online information deteriorate travel conditions for all users of that network relative to the situation where no online information is provided to users. The analytical presence of the paradox was derived for a specific network structure by using two equilibrium models, the first being the Expected User Equilibrium (EUE) solution (no information scenario) and the other being the User Equilibrium with Recourse (UER) solution (with information scenario). An incentivised computerised route choice game was designed using the concepts of experimental economics and administered in a controlled laboratory environment to investigate the physical presence of the paradox. Aggregate statistics of path flows and Total System Travel Costs (TSTC) were used to compare the experimental results with the theoretical findings. A total of 12 groups of 12 participants completed the experiment and the OIP and the occurrence of the OIP being significant was observed in 11 of the 12 cases. Though information increased travel costs for users on average, it reduced the volatility of travel costs experienced in the no information scenario indicating that information can achieve a more reliable system. Further replications of similar experiments and more importantly field based identification of the phenomena will force transport professionals to be aware of the emergence of the paradox. In addition, studies such as this emphasise the need for the adoption of adaptive traffic assignment techniques to appropriately model the acquisition of information on a road network.

Global Mittag-Leffler stability and synchronization analysis of fractional-order quaternion-valued neural networks with linear threshold neurons.

PubMed

Yang, Xujun; Li, Chuandong; Song, Qiankun; Chen, Jiyang; Huang, Junjian

2018-05-04

This paper talks about the stability and synchronization problems of fractional-order quaternion-valued neural networks (FQVNNs) with linear threshold neurons. On account of the non-commutativity of quaternion multiplication resulting from Hamilton rules, the FQVNN models are separated into four real-valued neural network (RVNN) models. Consequently, the dynamic analysis of FQVNNs can be realized by investigating the real-valued ones. Based on the method of M-matrix, the existence and uniqueness of the equilibrium point of the FQVNNs are obtained without detailed proof. Afterwards, several sufficient criteria ensuring the global Mittag-Leffler stability for the unique equilibrium point of the FQVNNs are derived by applying the Lyapunov direct method, the theory of fractional differential equation, the theory of matrix eigenvalue, and some inequality techniques. In the meanwhile, global Mittag-Leffler synchronization for the drive-response models of the addressed FQVNNs are investigated explicitly. Finally, simulation examples are designed to verify the feasibility and availability of the theoretical results. Copyright © 2018 Elsevier Ltd. All rights reserved.

Effects of Network Characteristics on Reaching the Payoff-Dominant Equilibrium in Coordination Games: A Simulation study.

PubMed

Buskens, Vincent; Snijders, Chris

2016-01-01

We study how payoffs and network structure affect reaching the payoff-dominant equilibrium in a [Formula: see text] coordination game that actors play with their neighbors in a network. Using an extensive simulation analysis of over 100,000 networks with 2-25 actors, we show that the importance of network characteristics is restricted to a limited part of the payoff space. In this part, we conclude that the payoff-dominant equilibrium is chosen more often if network density is larger, the network is more centralized, and segmentation of the network is smaller. Moreover, it is more likely that heterogeneity in behavior persists if the network is more segmented and less centralized. Persistence of heterogeneous behavior is not related to network density.

Self-Consistent Field Lattice Model for Polymer Networks.

PubMed

Tito, Nicholas B; Storm, Cornelis; Ellenbroek, Wouter G

2017-12-26

A lattice model based on polymer self-consistent field theory is developed to predict the equilibrium statistics of arbitrary polymer networks. For a given network topology, our approach uses moment propagators on a lattice to self-consistently construct the ensemble of polymer conformations and cross-link spatial probability distributions. Remarkably, the calculation can be performed "in the dark", without any prior knowledge on preferred chain conformations or cross-link positions. Numerical results from the model for a test network exhibit close agreement with molecular dynamics simulations, including when the network is strongly sheared. Our model captures nonaffine deformation, mean-field monomer interactions, cross-link fluctuations, and finite extensibility of chains, yielding predictions that differ markedly from classical rubber elasticity theory for polymer networks. By examining polymer networks with different degrees of interconnectivity, we gain insight into cross-link entropy, an important quantity in the macroscopic behavior of gels and self-healing materials as they are deformed.

Remarks on the chemical Fokker-Planck and Langevin equations: Nonphysical currents at equilibrium.

PubMed

Ceccato, Alessandro; Frezzato, Diego

2018-02-14

The chemical Langevin equation and the associated chemical Fokker-Planck equation are well-known continuous approximations of the discrete stochastic evolution of reaction networks. In this work, we show that these approximations suffer from a physical inconsistency, namely, the presence of nonphysical probability currents at the thermal equilibrium even for closed and fully detailed-balanced kinetic schemes. An illustration is given for a model case.

Selfish Routing

DTIC Science & Technology

2002-05-01

traffic models , thereby identifying types of networks for which the cost of routing selfishly is mild. The inefficiency inherent in an uncoordinated outcome...17 1.6 Bibliographic Notes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 2 Preliminaries 18 2.1 The Model ...to Other Models 68 4.1 Flows at Approximate Nash Equilibrium . . . . . . . . . . . . . . . . 69 4.2 Finitely Many Users: Splittable Flow

Punctuated equilibrium dynamics in human communications

NASA Astrophysics Data System (ADS)

Peng, Dan; Han, Xiao-Pu; Wei, Zong-Wen; Wang, Bing-Hong

2015-10-01

A minimal model based on network incorporating individual interactions is proposed to study the non-Poisson statistical properties of human behavior: individuals in system interact with their neighbors, the probability of an individual acting correlates to its activity, and all the individuals involved in action will change their activities randomly. The model reproduces varieties of spatial-temporal patterns observed in empirical studies of human daily communications, providing insight into various human activities and embracing a range of realistic social interacting systems, particularly, intriguing bimodal phenomenon. This model bridges priority queueing theory and punctuated equilibrium dynamics, and our modeling and analysis is likely to shed light on non-Poisson phenomena in many complex systems.

Evolutionary dynamics of group interactions on structured populations: a review

PubMed Central

Perc, Matjaž; Gómez-Gardeñes, Jesús; Szolnoki, Attila; Floría, Luis M.; Moreno, Yamir

2013-01-01

Interactions among living organisms, from bacteria colonies to human societies, are inherently more complex than interactions among particles and non-living matter. Group interactions are a particularly important and widespread class, representative of which is the public goods game. In addition, methods of statistical physics have proved valuable for studying pattern formation, equilibrium selection and self-organization in evolutionary games. Here, we review recent advances in the study of evolutionary dynamics of group interactions on top of structured populations, including lattices, complex networks and coevolutionary models. We also compare these results with those obtained on well-mixed populations. The review particularly highlights that the study of the dynamics of group interactions, like several other important equilibrium and non-equilibrium dynamical processes in biological, economical and social sciences, benefits from the synergy between statistical physics, network science and evolutionary game theory. PMID:23303223

A Compartmentalized Out-of-Equilibrium Enzymatic Reaction Network for Sustained Autonomous Movement

PubMed Central

2016-01-01

Every living cell is a compartmentalized out-of-equilibrium system exquisitely able to convert chemical energy into function. In order to maintain homeostasis, the flux of metabolites is tightly controlled by regulatory enzymatic networks. A crucial prerequisite for the development of lifelike materials is the construction of synthetic systems with compartmentalized reaction networks that maintain out-of-equilibrium function. Here, we aim for autonomous movement as an example of the conversion of feedstock molecules into function. The flux of the conversion is regulated by a rationally designed enzymatic reaction network with multiple feedforward loops. By compartmentalizing the network into bowl-shaped nanocapsules the output of the network is harvested as kinetic energy. The entire system shows sustained and tunable microscopic motion resulting from the conversion of multiple external substrates. The successful compartmentalization of an out-of-equilibrium reaction network is a major first step in harnessing the design principles of life for construction of adaptive and internally regulated lifelike systems. PMID:27924313

The effect of non-equilibrium metal cooling on the interstellar medium

NASA Astrophysics Data System (ADS)

Capelo, Pedro R.; Bovino, Stefano; Lupi, Alessandro; Schleicher, Dominik R. G.; Grassi, Tommaso

2018-04-01

By using a novel interface between the modern smoothed particle hydrodynamics code GASOLINE2 and the chemistry package KROME, we follow the hydrodynamical and chemical evolution of an isolated galaxy. In order to assess the relevance of different physical parameters and prescriptions, we constructed a suite of 10 simulations, in which we vary the chemical network (primordial and metal species), how metal cooling is modelled (non-equilibrium versus equilibrium; optically thin versus thick approximation), the initial gas metallicity (from 10 to 100 per cent solar), and how molecular hydrogen forms on dust. This is the first work in which metal injection from supernovae, turbulent metal diffusion, and a metal network with non-equilibrium metal cooling are self-consistently included in a galaxy simulation. We find that properly modelling the chemical evolution of several metal species and the corresponding non-equilibrium metal cooling has important effects on the thermodynamics of the gas, the chemical abundances, and the appearance of the galaxy: the gas is typically warmer, has a larger molecular-gas mass fraction, and has a smoother disc. We also conclude that, at relatively high metallicity, the choice of molecular-hydrogen formation rates on dust is not crucial. Moreover, we confirm that a higher initial metallicity produces a colder gas and a larger fraction of molecular gas, with the low-metallicity simulation best matching the observed molecular Kennicutt-Schmidt relation. Finally, our simulations agree quite well with observations that link star formation rate to metal emission lines.

Wireless Networks under a Backoff Attack: A Game Theoretical Perspective.

PubMed

Parras, Juan; Zazo, Santiago

2018-01-30

We study a wireless sensor network using CSMA/CA in the MAC layer under a backoff attack: some of the sensors of the network are malicious and deviate from the defined contention mechanism. We use Bianchi's network model to study the impact of the malicious sensors on the total network throughput, showing that it causes the throughput to be unfairly distributed among sensors. We model this conflict using game theory tools, where each sensor is a player. We obtain analytical solutions and propose an algorithm, based on Regret Matching, to learn the equilibrium of the game with an arbitrary number of players. Our approach is validated via simulations, showing that our theoretical predictions adjust to reality.

Neural network for solving convex quadratic bilevel programming problems.

PubMed

He, Xing; Li, Chuandong; Huang, Tingwen; Li, Chaojie

2014-03-01

In this paper, using the idea of successive approximation, we propose a neural network to solve convex quadratic bilevel programming problems (CQBPPs), which is modeled by a nonautonomous differential inclusion. Different from the existing neural network for CQBPP, the model has the least number of state variables and simple structure. Based on the theory of nonsmooth analysis, differential inclusions and Lyapunov-like method, the limit equilibrium points sequence of the proposed neural networks can approximately converge to an optimal solution of CQBPP under certain conditions. Finally, simulation results on two numerical examples and the portfolio selection problem show the effectiveness and performance of the proposed neural network. Copyright © 2013 Elsevier Ltd. All rights reserved.

Computing molecular fluctuations in biochemical reaction systems based on a mechanistic, statistical theory of irreversible processes.

PubMed

Kulasiri, Don

2011-01-01

We discuss the quantification of molecular fluctuations in the biochemical reaction systems within the context of intracellular processes associated with gene expression. We take the molecular reactions pertaining to circadian rhythms to develop models of molecular fluctuations in this chapter. There are a significant number of studies on stochastic fluctuations in intracellular genetic regulatory networks based on single cell-level experiments. In order to understand the fluctuations associated with the gene expression in circadian rhythm networks, it is important to model the interactions of transcriptional factors with the E-boxes in the promoter regions of some of the genes. The pertinent aspects of a near-equilibrium theory that would integrate the thermodynamical and particle dynamic characteristics of intracellular molecular fluctuations would be discussed, and the theory is extended by using the theory of stochastic differential equations. We then model the fluctuations associated with the promoter regions using general mathematical settings. We implemented ubiquitous Gillespie's algorithms, which are used to simulate stochasticity in biochemical networks, for each of the motifs. Both the theory and the Gillespie's algorithms gave the same results in terms of the time evolution of means and variances of molecular numbers. As biochemical reactions occur far away from equilibrium-hence the use of the Gillespie algorithm-these results suggest that the near-equilibrium theory should be a good approximation for some of the biochemical reactions. © 2011 Elsevier Inc. All rights reserved.

Information spread in networks: Games, optimal control, and stabilization

NASA Astrophysics Data System (ADS)

Khanafer, Ali

This thesis focuses on designing efficient mechanisms for controlling information spread in networks. We consider two models for information spread. The first one is the well-known distributed averaging dynamics. The second model is a nonlinear one that describes virus spread in computer and biological networks. We seek to design optimal, robust, and stabilizing controllers under practical constraints. For distributed averaging networks, we study the interaction between a network designer and an adversary. We consider two types of attacks on the network. In Attack-I, the adversary strategically disconnects a set of links to prevent the nodes from reaching consensus. Meanwhile, the network designer assists the nodes in reaching consensus by changing the weights of a limited number of links in the network. We formulate two problems to describe this competition where the order in which the players act is reversed in the two problems. Although the canonical equations provided by the Pontryagin's Maximum Principle (MP) seem to be intractable, we provide an alternative characterization for the optimal strategies that makes connection to potential theory. Further, we provide a sufficient condition for the existence of a saddle-point equilibrium (SPE) for the underlying zero-sum game. In Attack-II, the designer and the adversary are both capable of altering the measurements of all nodes in the network by injecting global signals. We impose two constraints on both players: a power constraint and an energy constraint. We assume that the available energy to each player is not sufficient to operate at maximum power throughout the horizon of the game. We show the existence of an SPE and derive the optimal strategies in closed form for this attack scenario. As an alternative to the "network designer vs. adversary" framework, we investigate the possibility of stabilizing unknown network diffusion processes using a distributed mechanism, where the uncertainty is due to an attack on the network. To this end, we propose a distributed version of the classical logic-based supervisory control scheme. Given a network of agents whose dynamics contain unknown parameters, the distributed supervisory control scheme is used to assist the agents to converge to a certain set-point without requiring them to have explicit knowledge of that set-point. Unlike the classical supervisory control scheme where a centralized supervisor makes switching decisions among the candidate controllers, in our scheme, each agent is equipped with a local supervisor that switches among the available controllers. The switching decisions made at a certain agent depend only on the information from its neighboring agents. We provide sufficient conditions for stabilization and apply our framework to the distributed averaging problem in the presence of large modeling uncertainty. For infected networks, we study the stability properties of a susceptible-infected-susceptible (SIS) diffusion model, so-called the n-intertwined Markov model, over arbitrary network topologies. Similar to the majority of infection spread dynamics, this model exhibits a threshold phenomenon. When the curing rates in the network are high, the all-healthy state is the unique equilibrium over the network. Otherwise, an endemic equilibrium state emerges, where some infection remains within the network. Using notions from positive systems theory, we provide conditions for the global asymptotic stability of the equilibrium points in both cases over strongly and weakly connected directed networks based on the value of the basic reproduction number, a fundamental quantity in the study of epidemics. Furthermore, we demonstrate that the n-intertwined Markov model can be viewed as a best-response dynamical system of a concave game among the nodes. This characterization allows us to cast new infection spread dynamics; additionally, we provide a sufficient condition, for the global convergence to the all-healthy state, that can be checked in a distributed fashion. Moreover, we investigate the problem of stabilizing the network when the curing rates of a limited number of nodes can be controlled. In particular, we characterize the number of controllers required for a class of undirected graphs. We also design optimal controllers capable of minimizing the total infection in the network at minimum cost. Finally, we outline a set of open problems in the area of information spread control.

An improved correlation to predict molecular weight between crosslinks based on equilibrium degree of swelling of hydrogel networks.

PubMed

Jimenez-Vergara, Andrea C; Lewis, John; Hahn, Mariah S; Munoz-Pinto, Dany J

2018-04-01

Accurate characterization of hydrogel diffusional properties is of substantial importance for a range of biotechnological applications. The diffusional capacity of hydrogels has commonly been estimated using the average molecular weight between crosslinks (M c ), which is calculated based on the equilibrium degree of swelling. However, the existing correlation linking M c and equilibrium swelling fails to accurately reflect the diffusional properties of highly crosslinked hydrogel networks. Also, as demonstrated herein, the current model fails to accurately predict the diffusional properties of hydrogels when polymer concentration and molecular weight are varied simultaneously. To address these limitations, we evaluated the diffusional properties of 48 distinct hydrogel formulations using two different photoinitiator systems, employing molecular size exclusion as an alternative methodology to calculate average hydrogel mesh size. The resulting data were then utilized to develop a revised correlation between M c and hydrogel equilibrium swelling that substantially reduces the limitations associated with the current correlation. © 2017 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 106B: 1339-1348, 2018. © 2017 Wiley Periodicals, Inc.

Capacity-constrained traffic assignment in networks with residual queues

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lam, W.H.K.; Zhang, Y.

2000-04-01

This paper proposes a capacity-constrained traffic assignment model for strategic transport planning in which the steady-state user equilibrium principle is extended for road networks with residual queues. Therefore, the road-exit capacity and the queuing effects can be incorporated into the strategic transport model for traffic forecasting. The proposed model is applicable to the congested network particularly when the traffic demands exceeds the capacity of the network during the peak period. An efficient solution method is proposed for solving the steady-state traffic assignment problem with residual queues. Then a simple numerical example is employed to demonstrate the application of the proposedmore » model and solution method, while an example of a medium-sized arterial highway network in Sioux Falls, South Dakota, is used to test the applicability of the proposed solution to real problems.« less

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site

PubMed Central

Sigala, Paul A.; Fafarman, Aaron T.; Schwans, Jason P.; Fried, Stephen D.; Fenn, Timothy D.; Caaveiro, Jose M. M.; Pybus, Brandon; Ringe, Dagmar; Petsko, Gregory A.; Boxer, Steven G.; Herschlag, Daniel

2013-01-01

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol–Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations. PMID:23798390

A Game Theoretic Approach for Balancing Energy Consumption in Clustered Wireless Sensor Networks.

PubMed

Yang, Liu; Lu, Yinzhi; Xiong, Lian; Tao, Yang; Zhong, Yuanchang

2017-11-17

Clustering is an effective topology control method in wireless sensor networks (WSNs), since it can enhance the network lifetime and scalability. To prolong the network lifetime in clustered WSNs, an efficient cluster head (CH) optimization policy is essential to distribute the energy among sensor nodes. Recently, game theory has been introduced to model clustering. Each sensor node is considered as a rational and selfish player which will play a clustering game with an equilibrium strategy. Then it decides whether to act as the CH according to this strategy for a tradeoff between providing required services and energy conservation. However, how to get the equilibrium strategy while maximizing the payoff of sensor nodes has rarely been addressed to date. In this paper, we present a game theoretic approach for balancing energy consumption in clustered WSNs. With our novel payoff function, realistic sensor behaviors can be captured well. The energy heterogeneity of nodes is considered by incorporating a penalty mechanism in the payoff function, so the nodes with more energy will compete for CHs more actively. We have obtained the Nash equilibrium (NE) strategy of the clustering game through convex optimization. Specifically, each sensor node can achieve its own maximal payoff when it makes the decision according to this strategy. Through plenty of simulations, our proposed game theoretic clustering is proved to have a good energy balancing performance and consequently the network lifetime is greatly enhanced.

Assortative mating and mutation diffusion in spatial evolutionary systems

NASA Astrophysics Data System (ADS)

Paley, C. J.; Taraskin, S. N.; Elliott, S. R.

2010-04-01

The influence of spatial structure on the equilibrium properties of a sexual population model defined on networks is studied numerically. Using a small-world-like topology of the networks as an investigative tool, the contributions to the fitness of assortative mating and of global mutant spread properties are considered. Simple measures of nearest-neighbor correlations and speed of spread of mutants through the system have been used to confirm that both of these dynamics are important contributory factors to the fitness. It is found that assortative mating increases the fitness of populations. Quick global spread of favorable mutations is shown to be a key factor increasing the equilibrium fitness of populations.

Majority-Vote Model on Opinion-Dependent Network

NASA Astrophysics Data System (ADS)

Lima, F. W. S.

2013-09-01

We study a nonequilibrium model with up-down symmetry and a noise parameter q known as majority-vote model (MVM) of Oliveira 1992 on opinion-dependent network or Stauffer-Hohnisch-Pittnauer (SHP) networks. By Monte Carlo (MC) simulations and finite-size scaling relations the critical exponents β/ν, γ/ν and 1/ν and points qc and U* are obtained. After extensive simulations, we obtain β/ν = 0.230(3), γ/ν = 0.535(2) and 1/ν = 0.475(8). The calculated values of the critical noise parameter and Binder cumulant are qc = 0.166(3) and U* = 0.288(3). Within the error bars, the exponents obey the relation 2β/ν + γ/ν = 1 and the results presented here demonstrate that the MVM belongs to a different universality class than the equilibrium Ising model on SHP networks, but to the same class as majority-vote models on some other networks.

Equilibrium Propagation: Bridging the Gap between Energy-Based Models and Backpropagation

PubMed Central

Scellier, Benjamin; Bengio, Yoshua

2017-01-01

We introduce Equilibrium Propagation, a learning framework for energy-based models. It involves only one kind of neural computation, performed in both the first phase (when the prediction is made) and the second phase of training (after the target or prediction error is revealed). Although this algorithm computes the gradient of an objective function just like Backpropagation, it does not need a special computation or circuit for the second phase, where errors are implicitly propagated. Equilibrium Propagation shares similarities with Contrastive Hebbian Learning and Contrastive Divergence while solving the theoretical issues of both algorithms: our algorithm computes the gradient of a well-defined objective function. Because the objective function is defined in terms of local perturbations, the second phase of Equilibrium Propagation corresponds to only nudging the prediction (fixed point or stationary distribution) toward a configuration that reduces prediction error. In the case of a recurrent multi-layer supervised network, the output units are slightly nudged toward their target in the second phase, and the perturbation introduced at the output layer propagates backward in the hidden layers. We show that the signal “back-propagated” during this second phase corresponds to the propagation of error derivatives and encodes the gradient of the objective function, when the synaptic update corresponds to a standard form of spike-timing dependent plasticity. This work makes it more plausible that a mechanism similar to Backpropagation could be implemented by brains, since leaky integrator neural computation performs both inference and error back-propagation in our model. The only local difference between the two phases is whether synaptic changes are allowed or not. We also show experimentally that multi-layer recurrently connected networks with 1, 2, and 3 hidden layers can be trained by Equilibrium Propagation on the permutation-invariant MNIST task. PMID:28522969

Microeconomics-based resource allocation in overlay networks by using non-strategic behavior modeling

NASA Astrophysics Data System (ADS)

Analoui, Morteza; Rezvani, Mohammad Hossein

2011-01-01

Behavior modeling has recently been investigated for designing self-organizing mechanisms in the context of communication networks in order to exploit the natural selfishness of the users with the goal of maximizing the overall utility. In strategic behavior modeling, the users of the network are assumed to be game players who seek to maximize their utility with taking into account the decisions that the other players might make. The essential difference between the aforementioned researches and this work is that it incorporates the non-strategic decisions in order to design the mechanism for the overlay network. In this solution concept, the decisions that a peer might make does not affect the actions of the other peers at all. The theory of consumer-firm developed in microeconomics is a model of the non-strategic behavior that we have adopted in our research. Based on it, we have presented distributed algorithms for peers' "joining" and "leaving" operations. We have modeled the overlay network as a competitive economy in which the content provided by an origin server can be viewed as commodity and the origin server and the peers who multicast the content to their downside are considered as the firms. On the other hand, due to the dual role of the peers in the overlay network, they can be considered as the consumers as well. On joining to the overlay economy, each peer is provided with an income and tries to get hold of the service regardless to the behavior of the other peers. We have designed the scalable algorithms in such a way that the existence of equilibrium price (known as Walrasian equilibrium price) is guaranteed.

An experimental study of the Online Information Paradox: Does en-route information improve road network performance?

PubMed Central

2017-01-01

This study investigates the empirical presence of a theoretical transportation paradox, defined as the “Online Information Paradox” (OIP). The paradox suggests that, for certain road networks, the provision of online information deteriorate travel conditions for all users of that network relative to the situation where no online information is provided to users. The analytical presence of the paradox was derived for a specific network structure by using two equilibrium models, the first being the Expected User Equilibrium (EUE) solution (no information scenario) and the other being the User Equilibrium with Recourse (UER) solution (with information scenario). An incentivised computerised route choice game was designed using the concepts of experimental economics and administered in a controlled laboratory environment to investigate the physical presence of the paradox. Aggregate statistics of path flows and Total System Travel Costs (TSTC) were used to compare the experimental results with the theoretical findings. A total of 12 groups of 12 participants completed the experiment and the OIP and the occurrence of the OIP being significant was observed in 11 of the 12 cases. Though information increased travel costs for users on average, it reduced the volatility of travel costs experienced in the no information scenario indicating that information can achieve a more reliable system. Further replications of similar experiments and more importantly field based identification of the phenomena will force transport professionals to be aware of the emergence of the paradox. In addition, studies such as this emphasise the need for the adoption of adaptive traffic assignment techniques to appropriately model the acquisition of information on a road network. PMID:28902854

Thermodynamics of stoichiometric biochemical networks in living systems far from equilibrium.

PubMed

Qian, Hong; Beard, Daniel A

2005-04-22

The principles of thermodynamics apply to both equilibrium and nonequilibrium biochemical systems. The mathematical machinery of the classic thermodynamics, however, mainly applies to systems in equilibrium. We introduce a thermodynamic formalism for the study of metabolic biochemical reaction (open, nonlinear) networks in both time-dependent and time-independent nonequilibrium states. Classical concepts in equilibrium thermodynamics-enthalpy, entropy, and Gibbs free energy of biochemical reaction systems-are generalized to nonequilibrium settings. Chemical motive force, heat dissipation rate, and entropy production (creation) rate, key concepts in nonequilibrium systems, are introduced. Dynamic equations for the thermodynamic quantities are presented in terms of the key observables of a biochemical network: stoichiometric matrix Q, reaction fluxes J, and chemical potentials of species mu without evoking empirical rate laws. Energy conservation and the Second Law are established for steady-state and dynamic biochemical networks. The theory provides the physiochemical basis for analyzing large-scale metabolic networks in living organisms.

Synchrony-induced modes of oscillation of a neural field model

NASA Astrophysics Data System (ADS)

Esnaola-Acebes, Jose M.; Roxin, Alex; Avitabile, Daniele; Montbrió, Ernest

2017-11-01

We investigate the modes of oscillation of heterogeneous ring networks of quadratic integrate-and-fire (QIF) neurons with nonlocal, space-dependent coupling. Perturbations of the equilibrium state with a particular wave number produce transient standing waves with a specific temporal frequency, analogously to those in a tense string. In the neuronal network, the equilibrium corresponds to a spatially homogeneous, asynchronous state. Perturbations of this state excite the network's oscillatory modes, which reflect the interplay of episodes of synchronous spiking with the excitatory-inhibitory spatial interactions. In the thermodynamic limit, an exact low-dimensional neural field model describing the macroscopic dynamics of the network is derived. This allows us to obtain formulas for the Turing eigenvalues of the spatially homogeneous state and hence to obtain its stability boundary. We find that the frequency of each Turing mode depends on the corresponding Fourier coefficient of the synaptic pattern of connectivity. The decay rate instead is identical for all oscillation modes as a consequence of the heterogeneity-induced desynchronization of the neurons. Finally, we numerically compute the spectrum of spatially inhomogeneous solutions branching from the Turing bifurcation, showing that similar oscillatory modes operate in neural bump states and are maintained away from onset.

Synchrony-induced modes of oscillation of a neural field model.

PubMed

Esnaola-Acebes, Jose M; Roxin, Alex; Avitabile, Daniele; Montbrió, Ernest

2017-11-01

We investigate the modes of oscillation of heterogeneous ring networks of quadratic integrate-and-fire (QIF) neurons with nonlocal, space-dependent coupling. Perturbations of the equilibrium state with a particular wave number produce transient standing waves with a specific temporal frequency, analogously to those in a tense string. In the neuronal network, the equilibrium corresponds to a spatially homogeneous, asynchronous state. Perturbations of this state excite the network's oscillatory modes, which reflect the interplay of episodes of synchronous spiking with the excitatory-inhibitory spatial interactions. In the thermodynamic limit, an exact low-dimensional neural field model describing the macroscopic dynamics of the network is derived. This allows us to obtain formulas for the Turing eigenvalues of the spatially homogeneous state and hence to obtain its stability boundary. We find that the frequency of each Turing mode depends on the corresponding Fourier coefficient of the synaptic pattern of connectivity. The decay rate instead is identical for all oscillation modes as a consequence of the heterogeneity-induced desynchronization of the neurons. Finally, we numerically compute the spectrum of spatially inhomogeneous solutions branching from the Turing bifurcation, showing that similar oscillatory modes operate in neural bump states and are maintained away from onset.

Equilibrium point control of a monkey arm simulator by a fast learning tree structured artificial neural network.

PubMed

Dornay, M; Sanger, T D

1993-01-01

A planar 17 muscle model of the monkey's arm based on realistic biomechanical measurements was simulated on a Symbolics Lisp Machine. The simulator implements the equilibrium point hypothesis for the control of arm movements. Given initial and final desired positions, it generates a minimum-jerk desired trajectory of the hand and uses the backdriving algorithm to determine an appropriate sequence of motor commands to the muscles (Flash 1987; Mussa-Ivaldi et al. 1991; Dornay 1991b). These motor commands specify a temporal sequence of stable (attractive) equilibrium positions which lead to the desired hand movement. A strong disadvantage of the simulator is that it has no memory of previous computations. Determining the desired trajectory using the minimum-jerk model is instantaneous, but the laborious backdriving algorithm is slow, and can take up to one hour for some trajectories. The complexity of the required computations makes it a poor model for biological motor control. We propose a computationally simpler and more biologically plausible method for control which achieves the benefits of the backdriving algorithm. A fast learning, tree-structured network (Sanger 1991c) was trained to remember the knowledge obtained by the backdriving algorithm. The neural network learned the nonlinear mapping from a 2-dimensional cartesian planar hand position (x,y) to a 17-dimensional motor command space (u1, . . ., u17). Learning 20 training trajectories, each composed of 26 sample points [[x,y], [u1, . . ., u17] took only 20 min on a Sun-4 Sparc workstation. After the learning stage, new, untrained test trajectories as well as the original trajectories of the hand were given to the neural network as input. The network calculated the required motor commands for these movements. The resulting movements were close to the desired ones for both the training and test cases.

Network-scale dynamics of sediment mixtures: how do tectonics affect surface bed texture and channel slope?

NASA Astrophysics Data System (ADS)

Gasparini, N. M.; Bras, R. L.; Tucker, G. E.

2003-04-01

An alluvial channel's slope and bed texture are intimately linked. Along with fluvial discharge, these variables are the key players in setting alluvial transport rates. We know that both channel slope and mean grain size usually decrease downstream, but how sensitive are these variables to tectonic changes? Are basin concavity and downstream fining drastically disrupted during transitions from one tectonic regime to another? We explore these questions using the CHILD numerical landscape evolution model to generate alluvial networks composed of a sand and gravel mixture. The steady-state and transient patterns of both channel slope and sediment texture are investigated. The steady-state patterns in slope and sediment texture are verified independently by solving the erosion equations under equilibrium conditions, i.e. the case when the erosion rate is equal to the uplift rate across the entire landscape. The inclusion of surface texture as a free parameter (as opposed to just channel slope) leads to some surprising results. In all cases, an increase in uplift rate results in channel beds which are finer at equilibrium (for a given drainage area). Higher uplift rates imply larger equilibrium transport rates; this leads to finer channels that have a smaller critical shear stress to entrain material, and therefore more material can be transported for a given discharge (and channel slope). Changes in equilibrium slopes are less intuitive. An increase in uplift rates can cause channel slopes to increase, remain the same, or decrease, depending on model parameter values. In the surprising case in which equilibrium channel slopes decrease with increasing uplift rates, we suggest that surface texture changes more than compensate for the required increase in transport rates, causing channel slopes to decrease. These results highlight the important role of sediment grain size in determining transport rates and caution us against ignoring this important variable in fluvial networks.

Oligopolistic competition in wholesale electricity markets: Large-scale simulation and policy analysis using complementarity models

NASA Astrophysics Data System (ADS)

Helman, E. Udi

This dissertation conducts research into the large-scale simulation of oligopolistic competition in wholesale electricity markets. The dissertation has two parts. Part I is an examination of the structure and properties of several spatial, or network, equilibrium models of oligopolistic electricity markets formulated as mixed linear complementarity problems (LCP). Part II is a large-scale application of such models to the electricity system that encompasses most of the United States east of the Rocky Mountains, the Eastern Interconnection. Part I consists of Chapters 1 to 6. The models developed in this part continue research into mixed LCP models of oligopolistic electricity markets initiated by Hobbs [67] and subsequently developed by Metzler [87] and Metzler, Hobbs and Pang [88]. Hobbs' central contribution is a network market model with Cournot competition in generation and a price-taking spatial arbitrage firm that eliminates spatial price discrimination by the Cournot firms. In one variant, the solution to this model is shown to be equivalent to the "no arbitrage" condition in a "pool" market, in which a Regional Transmission Operator optimizes spot sales such that the congestion price between two locations is exactly equivalent to the difference in the energy prices at those locations (commonly known as locational marginal pricing). Extensions to this model are presented in Chapters 5 and 6. One of these is a market model with a profit-maximizing arbitrage firm. This model is structured as a mathematical program with equilibrium constraints (MPEC), but due to the linearity of its constraints, can be solved as a mixed LCP. Part II consists of Chapters 7 to 12. The core of these chapters is a large-scale simulation of the U.S. Eastern Interconnection applying one of the Cournot competition with arbitrage models. This is the first oligopolistic equilibrium market model to encompass the full Eastern Interconnection with a realistic network representation (using a DC load flow approximation). Chapter 9 shows the price results. In contrast to prior market power simulations of these markets, much greater variability in price-cost margins is found when using a realistic model of hourly conditions on such a large network. Chapter 10 shows that the conventional concentration indices (HHIs) are poorly correlated with PCMs. Finally, Chapter 11 proposes that the simulation models are applied to merger analysis and provides two large-scale merger examples. (Abstract shortened by UMI.)

Robustness analysis of uncertain dynamical neural networks with multiple time delays.

PubMed

Senan, Sibel

2015-10-01

This paper studies the problem of global robust asymptotic stability of the equilibrium point for the class of dynamical neural networks with multiple time delays with respect to the class of slope-bounded activation functions and in the presence of the uncertainties of system parameters of the considered neural network model. By using an appropriate Lyapunov functional and exploiting the properties of the homeomorphism mapping theorem, we derive a new sufficient condition for the existence, uniqueness and global robust asymptotic stability of the equilibrium point for the class of neural networks with multiple time delays. The obtained stability condition basically relies on testing some relationships imposed on the interconnection matrices of the neural system, which can be easily verified by using some certain properties of matrices. An instructive numerical example is also given to illustrate the applicability of our result and show the advantages of this new condition over the previously reported corresponding results. Copyright © 2015 Elsevier Ltd. All rights reserved.

Linear network representation of multistate models of transport.

PubMed Central

Sandblom, J; Ring, A; Eisenman, G

1982-01-01

By introducing external driving forces in rate-theory models of transport we show how the Eyring rate equations can be transformed into Ohm's law with potentials that obey Kirchhoff's second law. From such a formalism the state diagram of a multioccupancy multicomponent system can be directly converted into linear network with resistors connecting nodal (branch) points and with capacitances connecting each nodal point with a reference point. The external forces appear as emf or current generators in the network. This theory allows the algebraic methods of linear network theory to be used in solving the flux equations for multistate models and is particularly useful for making proper simplifying approximation in models of complex membrane structure. Some general properties of linear network representation are also deduced. It is shown, for instance, that Maxwell's reciprocity relationships of linear networks lead directly to Onsager's relationships in the near equilibrium region. Finally, as an example of the procedure, the equivalent circuit method is used to solve the equations for a few transport models. PMID:7093425

Wireless Networks under a Backoff Attack: A Game Theoretical Perspective

PubMed Central

Zazo, Santiago

2018-01-01

We study a wireless sensor network using CSMA/CA in the MAC layer under a backoff attack: some of the sensors of the network are malicious and deviate from the defined contention mechanism. We use Bianchi’s network model to study the impact of the malicious sensors on the total network throughput, showing that it causes the throughput to be unfairly distributed among sensors. We model this conflict using game theory tools, where each sensor is a player. We obtain analytical solutions and propose an algorithm, based on Regret Matching, to learn the equilibrium of the game with an arbitrary number of players. Our approach is validated via simulations, showing that our theoretical predictions adjust to reality. PMID:29385752

A novel multiphysic model for simulation of swelling equilibrium of ionized thermal-stimulus responsive hydrogels

NASA Astrophysics Data System (ADS)

Li, Hua; Wang, Xiaogui; Yan, Guoping; Lam, K. Y.; Cheng, Sixue; Zou, Tao; Zhuo, Renxi

2005-03-01

In this paper, a novel multiphysic mathematical model is developed for simulation of swelling equilibrium of ionized temperature sensitive hydrogels with the volume phase transition, and it is termed the multi-effect-coupling thermal-stimulus (MECtherm) model. This model consists of the steady-state Nernst-Planck equation, Poisson equation and swelling equilibrium governing equation based on the Flory's mean field theory, in which two types of polymer-solvent interaction parameters, as the functions of temperature and polymer-network volume fraction, are specified with or without consideration of the hydrogen bond interaction. In order to examine the MECtherm model consisting of nonlinear partial differential equations, a meshless Hermite-Cloud method is used for numerical solution of one-dimensional swelling equilibrium of thermal-stimulus responsive hydrogels immersed in a bathing solution. The computed results are in very good agreements with experimental data for the variation of volume swelling ratio with temperature. The influences of the salt concentration and initial fixed-charge density are discussed in detail on the variations of volume swelling ratio of hydrogels, mobile ion concentrations and electric potential of both interior hydrogels and exterior bathing solution.

Active Tension Network model reveals an exotic mechanical state realized in epithelial tissues

NASA Astrophysics Data System (ADS)

Noll, Nicholas; Mani, Madhav; Heemskerk, Idse; Streicha, Sebastian; Shraiman, Boris

Mechanical interactions play a crucial role in epithelial morphogenesis, yet understanding the complex mechanisms through which stress and deformation affect cell behavior remains an open problem. Here we formulate and analyze the Active Tension Network (ATN) model, which assumes that mechanical balance of cells is dominated by cortical tension and introduces tension dependent active remodeling of the cortex. We find that ATNs exhibit unusual mechanical properties: i) ATN behaves as a fluid at short times, but at long times it supports external tension, like a solid; ii) its mechanical equilibrium state has extensive degeneracy associated with a discrete conformal - ''isogonal'' - deformation of cells. ATN model predicts a constraint on equilibrium cell geometry, which we demonstrate to hold in certain epithelial tissues. We further show that isogonal modes are observed in a fruit fly embryo, accounting for the striking variability of apical area of ventral cells and helping understand the early phase of gastrulation. Living matter realizes new and exotic mechanical states, understanding which helps understand biological phenomena.

Spatial price dynamics: From complex network perspective

NASA Astrophysics Data System (ADS)

Li, Y. L.; Bi, J. T.; Sun, H. J.

2008-10-01

The spatial price problem means that if the supply price plus the transportation cost is less than the demand price, there exists a trade. Thus, after an amount of exchange, the demand price will decrease. This process is continuous until an equilibrium state is obtained. However, how the trade network structure affects this process has received little attention. In this paper, we give a evolving model to describe the levels of spatial price on different complex network structures. The simulation results show that the network with shorter path length is sensitive to the variation of prices.

An Amorphous Model for Morphological Processing in Visual Comprehension Based on Naive Discriminative Learning

ERIC Educational Resources Information Center

Baayen, R. Harald; Milin, Petar; Durdevic, Dusica Filipovic; Hendrix, Peter; Marelli, Marco

2011-01-01

A 2-layer symbolic network model based on the equilibrium equations of the Rescorla-Wagner model (Danks, 2003) is proposed. The study first presents 2 experiments in Serbian, which reveal for sentential reading the inflectional paradigmatic effects previously observed by Milin, Filipovic Durdevic, and Moscoso del Prado Martin (2009) for unprimed…

Rational design of functional and tunable oscillating enzymatic networks

NASA Astrophysics Data System (ADS)

Semenov, Sergey N.; Wong, Albert S. Y.; van der Made, R. Martijn; Postma, Sjoerd G. J.; Groen, Joost; van Roekel, Hendrik W. H.; de Greef, Tom F. A.; Huck, Wilhelm T. S.

2015-02-01

Life is sustained by complex systems operating far from equilibrium and consisting of a multitude of enzymatic reaction networks. The operating principles of biology's regulatory networks are known, but the in vitro assembly of out-of-equilibrium enzymatic reaction networks has proved challenging, limiting the development of synthetic systems showing autonomous behaviour. Here, we present a strategy for the rational design of programmable functional reaction networks that exhibit dynamic behaviour. We demonstrate that a network built around autoactivation and delayed negative feedback of the enzyme trypsin is capable of producing sustained oscillating concentrations of active trypsin for over 65 h. Other functions, such as amplification, analog-to-digital conversion and periodic control over equilibrium systems, are obtained by linking multiple network modules in microfluidic flow reactors. The methodology developed here provides a general framework to construct dissipative, tunable and robust (bio)chemical reaction networks.

Random graph models for dynamic networks

NASA Astrophysics Data System (ADS)

Zhang, Xiao; Moore, Cristopher; Newman, Mark E. J.

2017-10-01

Recent theoretical work on the modeling of network structure has focused primarily on networks that are static and unchanging, but many real-world networks change their structure over time. There exist natural generalizations to the dynamic case of many static network models, including the classic random graph, the configuration model, and the stochastic block model, where one assumes that the appearance and disappearance of edges are governed by continuous-time Markov processes with rate parameters that can depend on properties of the nodes. Here we give an introduction to this class of models, showing for instance how one can compute their equilibrium properties. We also demonstrate their use in data analysis and statistical inference, giving efficient algorithms for fitting them to observed network data using the method of maximum likelihood. This allows us, for example, to estimate the time constants of network evolution or infer community structure from temporal network data using cues embedded both in the probabilities over time that node pairs are connected by edges and in the characteristic dynamics of edge appearance and disappearance. We illustrate these methods with a selection of applications, both to computer-generated test networks and real-world examples.

A Game Theoretic Approach for Balancing Energy Consumption in Clustered Wireless Sensor Networks

PubMed Central

Lu, Yinzhi; Xiong, Lian; Tao, Yang; Zhong, Yuanchang

2017-01-01

Clustering is an effective topology control method in wireless sensor networks (WSNs), since it can enhance the network lifetime and scalability. To prolong the network lifetime in clustered WSNs, an efficient cluster head (CH) optimization policy is essential to distribute the energy among sensor nodes. Recently, game theory has been introduced to model clustering. Each sensor node is considered as a rational and selfish player which will play a clustering game with an equilibrium strategy. Then it decides whether to act as the CH according to this strategy for a tradeoff between providing required services and energy conservation. However, how to get the equilibrium strategy while maximizing the payoff of sensor nodes has rarely been addressed to date. In this paper, we present a game theoretic approach for balancing energy consumption in clustered WSNs. With our novel payoff function, realistic sensor behaviors can be captured well. The energy heterogeneity of nodes is considered by incorporating a penalty mechanism in the payoff function, so the nodes with more energy will compete for CHs more actively. We have obtained the Nash equilibrium (NE) strategy of the clustering game through convex optimization. Specifically, each sensor node can achieve its own maximal payoff when it makes the decision according to this strategy. Through plenty of simulations, our proposed game theoretic clustering is proved to have a good energy balancing performance and consequently the network lifetime is greatly enhanced. PMID:29149075

Transportation network policy modeling for congestion and pollution control: A variational inequality approach

NASA Astrophysics Data System (ADS)

Ramanujam, Padma

1999-08-01

Public concern over the state of the environment has grown over the past decade. All indications are that this concern will continue to influence policy making into the foreseeable future. Road transport is seen as the major contributor to environmental degradation. Transportation planners around the world face the question: cleaner air and/or faster commutes? While individual vehicles can be made more environmentally friendly, the sheer scale of growth in world-wide vehicle numbers is projected to cause significant environmental degradation in the longer run, and in the absence of newer and stricter polices. It is a challenge for governments to find policies that ensure congestion-free metropolitan areas while guaranteeing both critical environmental quality levels and a sufficient infrastructure access to all groups involved. The objective of the dissertation is to provide a mathematical framework to study transportation policy models for the purpose of controlling congestion and pollution. Towards this objective. a series of transportation policy models are developed to study travel behavior and to quantity the reductions in congestion and automobile emissions. The dissertation begins with a brief historical overview of some of the pioneering works in urban transportation economics and later presents the theoretical foundation for the transportation policy models developed. The dissertation introduces single modal and multimodal transportation network policy models that accomplish road pricing with the imposition of goal targets on link loads. as well as, integrated traffic equilibrium models with marketable mobile emission permits. Furthermore, equilibrium conditions are derived for each model, and both qualitative analysis and computational procedures are studied. Finally, the dissertation concludes with a comparative study of the relationship between regulatory pricing models and marketable emission permit transportation models and a discussion on key factors that influence implementation of the proposed policy models. The framework of variational inequalities has been utilized in our dissertation, because it is ideal for equilibrium systems. With the addition of pricing policy interventions and the integration of marketable mobile emission permits, traffic equilibrium models become extremely complex. Consequently, the computation of the equilibrium is made more difficult. However, it is shown in the dissertation that in addition to pricing interventions and the integration of a marketable emission permit system that it is possible to incorporate multiple modes of transport and even to handle the issue of noncompliance, using the framework of variational inequalities.

Local Nash equilibrium in social networks.

PubMed

Zhang, Yichao; Aziz-Alaoui, M A; Bertelle, Cyrille; Guan, Jihong

2014-08-29

Nash equilibrium is widely present in various social disputes. As of now, in structured static populations, such as social networks, regular, and random graphs, the discussions on Nash equilibrium are quite limited. In a relatively stable static gaming network, a rational individual has to comprehensively consider all his/her opponents' strategies before they adopt a unified strategy. In this scenario, a new strategy equilibrium emerges in the system. We define this equilibrium as a local Nash equilibrium. In this paper, we present an explicit definition of the local Nash equilibrium for the two-strategy games in structured populations. Based on the definition, we investigate the condition that a system reaches the evolutionary stable state when the individuals play the Prisoner's dilemma and snow-drift game. The local Nash equilibrium provides a way to judge whether a gaming structured population reaches the evolutionary stable state on one hand. On the other hand, it can be used to predict whether cooperators can survive in a system long before the system reaches its evolutionary stable state for the Prisoner's dilemma game. Our work therefore provides a theoretical framework for understanding the evolutionary stable state in the gaming populations with static structures.

Local Nash Equilibrium in Social Networks

PubMed Central

Zhang, Yichao; Aziz-Alaoui, M. A.; Bertelle, Cyrille; Guan, Jihong

2014-01-01

Nash equilibrium is widely present in various social disputes. As of now, in structured static populations, such as social networks, regular, and random graphs, the discussions on Nash equilibrium are quite limited. In a relatively stable static gaming network, a rational individual has to comprehensively consider all his/her opponents' strategies before they adopt a unified strategy. In this scenario, a new strategy equilibrium emerges in the system. We define this equilibrium as a local Nash equilibrium. In this paper, we present an explicit definition of the local Nash equilibrium for the two-strategy games in structured populations. Based on the definition, we investigate the condition that a system reaches the evolutionary stable state when the individuals play the Prisoner's dilemma and snow-drift game. The local Nash equilibrium provides a way to judge whether a gaming structured population reaches the evolutionary stable state on one hand. On the other hand, it can be used to predict whether cooperators can survive in a system long before the system reaches its evolutionary stable state for the Prisoner's dilemma game. Our work therefore provides a theoretical framework for understanding the evolutionary stable state in the gaming populations with static structures. PMID:25169150

Local Nash Equilibrium in Social Networks

NASA Astrophysics Data System (ADS)

Zhang, Yichao; Aziz-Alaoui, M. A.; Bertelle, Cyrille; Guan, Jihong

2014-08-01

Nash equilibrium is widely present in various social disputes. As of now, in structured static populations, such as social networks, regular, and random graphs, the discussions on Nash equilibrium are quite limited. In a relatively stable static gaming network, a rational individual has to comprehensively consider all his/her opponents' strategies before they adopt a unified strategy. In this scenario, a new strategy equilibrium emerges in the system. We define this equilibrium as a local Nash equilibrium. In this paper, we present an explicit definition of the local Nash equilibrium for the two-strategy games in structured populations. Based on the definition, we investigate the condition that a system reaches the evolutionary stable state when the individuals play the Prisoner's dilemma and snow-drift game. The local Nash equilibrium provides a way to judge whether a gaming structured population reaches the evolutionary stable state on one hand. On the other hand, it can be used to predict whether cooperators can survive in a system long before the system reaches its evolutionary stable state for the Prisoner's dilemma game. Our work therefore provides a theoretical framework for understanding the evolutionary stable state in the gaming populations with static structures.

Multistability of memristive Cohen-Grossberg neural networks with non-monotonic piecewise linear activation functions and time-varying delays.

PubMed

Nie, Xiaobing; Zheng, Wei Xing; Cao, Jinde

2015-11-01

The problem of coexistence and dynamical behaviors of multiple equilibrium points is addressed for a class of memristive Cohen-Grossberg neural networks with non-monotonic piecewise linear activation functions and time-varying delays. By virtue of the fixed point theorem, nonsmooth analysis theory and other analytical tools, some sufficient conditions are established to guarantee that such n-dimensional memristive Cohen-Grossberg neural networks can have 5(n) equilibrium points, among which 3(n) equilibrium points are locally exponentially stable. It is shown that greater storage capacity can be achieved by neural networks with the non-monotonic activation functions introduced herein than the ones with Mexican-hat-type activation function. In addition, unlike most existing multistability results of neural networks with monotonic activation functions, those obtained 3(n) locally stable equilibrium points are located both in saturated regions and unsaturated regions. The theoretical findings are verified by an illustrative example with computer simulations. Copyright © 2015 Elsevier Ltd. All rights reserved.

Effects of Heterogeneous Social Interactions on Flocking Dynamics

NASA Astrophysics Data System (ADS)

Miguel, M. Carmen; Parley, Jack T.; Pastor-Satorras, Romualdo

2018-02-01

Social relationships characterize the interactions that occur within social species and may have an important impact on collective animal motion. Here, we consider a variation of the standard Vicsek model for collective motion in which interactions are mediated by an empirically motivated scale-free topology that represents a heterogeneous pattern of social contacts. We observe that the degree of order of the model is strongly affected by network heterogeneity: more heterogeneous networks show a more resilient ordered state, while less heterogeneity leads to a more fragile ordered state that can be destroyed by sufficient external noise. Our results challenge the previously accepted equivalence between the static Vicsek model and the equilibrium X Y model on the network of connections, and point towards a possible equivalence with models exhibiting a different symmetry.

Dynamics, morphogenesis and convergence of evolutionary quantum Prisoner's Dilemma games on networks

PubMed Central

Yong, Xi

2016-01-01

The authors proposed a quantum Prisoner's Dilemma (PD) game as a natural extension of the classic PD game to resolve the dilemma. Here, we establish a new Nash equilibrium principle of the game, propose the notion of convergence and discover the convergence and phase-transition phenomena of the evolutionary games on networks. We investigate the many-body extension of the game or evolutionary games in networks. For homogeneous networks, we show that entanglement guarantees a quick convergence of super cooperation, that there is a phase transition from the convergence of defection to the convergence of super cooperation, and that the threshold for the phase transitions is principally determined by the Nash equilibrium principle of the game, with an accompanying perturbation by the variations of structures of networks. For heterogeneous networks, we show that the equilibrium frequencies of super-cooperators are divergent, that entanglement guarantees emergence of super-cooperation and that there is a phase transition of the emergence with the threshold determined by the Nash equilibrium principle, accompanied by a perturbation by the variations of structures of networks. Our results explore systematically, for the first time, the dynamics, morphogenesis and convergence of evolutionary games in interacting and competing systems. PMID:27118882

Spreading out of perturbations in reversible reaction networks

NASA Astrophysics Data System (ADS)

Maslov, Sergei; Sneppen, Kim; Ispolatov, I.

2007-08-01

Using an example of physical interactions between proteins, we study how a perturbation propagates in the equilibrium of a network of reversible reactions governed by the law of mass action. We introduce a matrix formalism to describe the linear response of all equilibrium concentrations to shifts in total abundances of individual reactants, and reveal its heuristic analogy to the flow of electric current in a network of resistors. Our main conclusion is that, on average, the induced changes in equilibrium concentrations decay exponentially as a function of network distance from the source of perturbation. We analyze how this decay is influenced by such factors as the topology of a network, binding strength, and correlations between concentrations of neighboring nodes. We find that the minimal branching of the network, small values of dissociation constants, and low equilibrium free (unbound) concentrations of reacting substances all decrease the decay constant and thus increase the range of propagation. Exact analytic expressions for the decay constant are obtained for the case of equally strong interactions and uniform as well as oscillating concentrations on the Bethe lattice. Our general findings are illustrated using a real network of protein-protein interactions in baker's yeast with experimentally determined protein concentrations.

Biological signatures of dynamic river networks from a coupled landscape evolution and neutral community model

NASA Astrophysics Data System (ADS)

Stokes, M.; Perron, J. T.

2017-12-01

Freshwater systems host exceptionally species-rich communities whose spatial structure is dictated by the topology of the river networks they inhabit. Over geologic time, river networks are dynamic; drainage basins shrink and grow, and river capture establishes new connections between previously separated regions. It has been hypothesized that these changes in river network structure influence the evolution of life by exchanging and isolating species, perhaps boosting biodiversity in the process. However, no general model exists to predict the evolutionary consequences of landscape change. We couple a neutral community model of freshwater organisms to a landscape evolution model in which the river network undergoes drainage divide migration and repeated river capture. Neutral community models are macro-ecological models that include stochastic speciation and dispersal to produce realistic patterns of biodiversity. We explore the consequences of three modes of speciation - point mutation, time-protracted, and vicariant (geographic) speciation - by tracking patterns of diversity in time and comparing the final result to an equilibrium solution of the neutral model on the final landscape. Under point mutation, a simple model of stochastic and instantaneous speciation, the results are identical to the equilibrium solution and indicate the dominance of the species-area relationship in forming patterns of diversity. The number of species in a basin is proportional to its area, and regional species richness reaches its maximum when drainage area is evenly distributed among sub-basins. Time-protracted speciation is also modeled as a stochastic process, but in order to produce more realistic rates of diversification, speciation is not assumed to be instantaneous. Rather, each new species must persist for a certain amount of time before it is considered to be established. When vicariance (geographic speciation) is included, there is a transient signature of increased regional diversity after river capture. The results indicate that the mode of speciation and the rate of speciation relative to the rate of divide migration determine the evolutionary signature of river capture.

Macroscopic description of complex adaptive networks coevolving with dynamic node states

NASA Astrophysics Data System (ADS)

Wiedermann, Marc; Donges, Jonathan F.; Heitzig, Jobst; Lucht, Wolfgang; Kurths, Jürgen

2015-05-01

In many real-world complex systems, the time evolution of the network's structure and the dynamic state of its nodes are closely entangled. Here we study opinion formation and imitation on an adaptive complex network which is dependent on the individual dynamic state of each node and vice versa to model the coevolution of renewable resources with the dynamics of harvesting agents on a social network. The adaptive voter model is coupled to a set of identical logistic growth models and we mainly find that, in such systems, the rate of interactions between nodes as well as the adaptive rewiring probability are crucial parameters for controlling the sustainability of the system's equilibrium state. We derive a macroscopic description of the system in terms of ordinary differential equations which provides a general framework to model and quantify the influence of single node dynamics on the macroscopic state of the network. The thus obtained framework is applicable to many fields of study, such as epidemic spreading, opinion formation, or socioecological modeling.

Macroscopic description of complex adaptive networks coevolving with dynamic node states.

PubMed

Wiedermann, Marc; Donges, Jonathan F; Heitzig, Jobst; Lucht, Wolfgang; Kurths, Jürgen

2015-05-01

In many real-world complex systems, the time evolution of the network's structure and the dynamic state of its nodes are closely entangled. Here we study opinion formation and imitation on an adaptive complex network which is dependent on the individual dynamic state of each node and vice versa to model the coevolution of renewable resources with the dynamics of harvesting agents on a social network. The adaptive voter model is coupled to a set of identical logistic growth models and we mainly find that, in such systems, the rate of interactions between nodes as well as the adaptive rewiring probability are crucial parameters for controlling the sustainability of the system's equilibrium state. We derive a macroscopic description of the system in terms of ordinary differential equations which provides a general framework to model and quantify the influence of single node dynamics on the macroscopic state of the network. The thus obtained framework is applicable to many fields of study, such as epidemic spreading, opinion formation, or socioecological modeling.

Analysing and controlling the tax evasion dynamics via majority-vote model

NASA Astrophysics Data System (ADS)

Lima, F. W. S.

2010-09-01

Within the context of agent-based Monte-Carlo simulations, we study the well-known majority-vote model (MVM) with noise applied to tax evasion on simple square lattices, Voronoi-Delaunay random lattices, Barabasi-Albert networks, and Erdös-Rényi random graphs. In the order to analyse and to control the fluctuations for tax evasion in the economics model proposed by Zaklan, MVM is applied in the neighborhod of the noise critical qc to evolve the Zaklan model. The Zaklan model had been studied recently using the equilibrium Ising model. Here we show that the Zaklan model is robust because this can be studied using equilibrium dynamics of Ising model also through the nonequilibrium MVM and on various topologies cited above giving the same behavior regardless of dynamic or topology used here.

Dynamics of quality as a strategic variable in complex food supply chain network competition: The case of fresh produce

NASA Astrophysics Data System (ADS)

Nagurney, Anna; Besik, Deniz; Yu, Min

2018-04-01

In this paper, we construct a competitive food supply chain network model in which the profit-maximizing producers decide not only as to the volume of fresh produce produced and distributed using various supply chain network pathways, but they also decide, with the associated costs, on the initial quality of the fresh produce. Consumers, in turn, respond to the various producers' product outputs through the prices that they are willing to pay, given also the average quality associated with each producer or brand at the retail outlets. The quality of the fresh produce is captured through explicit formulae that incorporate time, temperature, and other link characteristics with links associated with processing, shipment, storage, etc. Capacities on links are also incorporated as well as upper bounds on the initial product quality of the firms at their production/harvesting sites. The governing concept of the competitive supply chain network model is that of Nash Equilibrium, for which alternative variational inequality formulations are derived, along with existence results. An algorithmic procedure, which can be interpreted as a discrete-time tatonnement process, is then described and applied to compute the equilibrium produce flow patterns and accompanying link Lagrange multipliers in a realistic case study, focusing on peaches, which includes disruptions.

The epidemic threshold theorem with social and contact heterogeneity

NASA Astrophysics Data System (ADS)

Hincapié Palacio, Doracelly; Ospina Giraldo, Juan; Gómez Arias, Rubén Darío

2008-03-01

The threshold theorem of an epidemic SIR model was compared when infectious and susceptible individuals have homogeneous mixing and heterogeneous social status and when individuals of random networks have contact heterogeneity. Particularly the effect of vaccination in such models is considered when: individuals or nodes are exposed to impoverished, vaccination and loss of immunity. An equilibrium analysis and local stability of small perturbations about the equilibrium values were implemented using computer algebra. Numerical simulations were executed in order to describe the dynamic of transmission of diseases and changes of the basic reproductive rate. The implications of these results are examined around the threats to the global public health security.

Optimum load distribution between heat sources based on the Cournot model

NASA Astrophysics Data System (ADS)

Penkovskii, A. V.; Stennikov, V. A.; Khamisov, O. V.

2015-08-01

One of the widespread models of the heat supply of consumers, which is represented in the "Single buyer" format, is considered. The methodological base proposed for its description and investigation presents the use of principles of the theory of games, basic propositions of microeconomics, and models and methods of the theory of hydraulic circuits. The original mathematical model of the heat supply system operating under conditions of the "Single buyer" organizational structure provides the derivation of a solution satisfying the market Nash equilibrium. The distinctive feature of the developed mathematical model is that, along with problems solved traditionally within the bounds of bilateral relations of heat energy sources-heat consumer, it considers a network component with its inherent physicotechnical properties of the heat network and business factors connected with costs of the production and transportation of heat energy. This approach gives the possibility to determine optimum levels of load of heat energy sources. These levels provide the given heat energy demand of consumers subject to the maximum profit earning of heat energy sources and the fulfillment of conditions for formation of minimum heat network costs for a specified time. The practical realization of the search of market equilibrium is considered by the example of a heat supply system with two heat energy sources operating on integrated heat networks. The mathematical approach to the solution search is represented in the graphical form and illustrates computations based on the stepwise iteration procedure for optimization of levels of loading of heat energy sources (groping procedure by Cournot) with the corresponding computation of the heat energy price for consumers.

Stochastic cycle selection in active flow networks.

PubMed

Woodhouse, Francis G; Forrow, Aden; Fawcett, Joanna B; Dunkel, Jörn

2016-07-19

Active biological flow networks pervade nature and span a wide range of scales, from arterial blood vessels and bronchial mucus transport in humans to bacterial flow through porous media or plasmodial shuttle streaming in slime molds. Despite their ubiquity, little is known about the self-organization principles that govern flow statistics in such nonequilibrium networks. Here we connect concepts from lattice field theory, graph theory, and transition rate theory to understand how topology controls dynamics in a generic model for actively driven flow on a network. Our combined theoretical and numerical analysis identifies symmetry-based rules that make it possible to classify and predict the selection statistics of complex flow cycles from the network topology. The conceptual framework developed here is applicable to a broad class of biological and nonbiological far-from-equilibrium networks, including actively controlled information flows, and establishes a correspondence between active flow networks and generalized ice-type models.

Stochastic cycle selection in active flow networks

PubMed Central

Woodhouse, Francis G.; Forrow, Aden; Fawcett, Joanna B.; Dunkel, Jörn

2016-01-01

Active biological flow networks pervade nature and span a wide range of scales, from arterial blood vessels and bronchial mucus transport in humans to bacterial flow through porous media or plasmodial shuttle streaming in slime molds. Despite their ubiquity, little is known about the self-organization principles that govern flow statistics in such nonequilibrium networks. Here we connect concepts from lattice field theory, graph theory, and transition rate theory to understand how topology controls dynamics in a generic model for actively driven flow on a network. Our combined theoretical and numerical analysis identifies symmetry-based rules that make it possible to classify and predict the selection statistics of complex flow cycles from the network topology. The conceptual framework developed here is applicable to a broad class of biological and nonbiological far-from-equilibrium networks, including actively controlled information flows, and establishes a correspondence between active flow networks and generalized ice-type models. PMID:27382186

Fluctuations in Mass-Action Equilibrium of Protein Binding Networks

NASA Astrophysics Data System (ADS)

Yan, Koon-Kiu; Walker, Dylan; Maslov, Sergei

2008-12-01

We consider two types of fluctuations in the mass-action equilibrium in protein binding networks. The first type is driven by slow changes in total concentrations of interacting proteins. The second type (spontaneous) is caused by quickly decaying thermodynamic deviations away from equilibrium. We investigate the effects of network connectivity on fluctuations by comparing them to scenarios in which the interacting pair is isolated from the network and analytically derives bounds on fluctuations. Collective effects are shown to sometimes lead to large amplification of spontaneous fluctuations. The strength of both types of fluctuations is positively correlated with the complex connectivity and negatively correlated with complex concentration. Our general findings are illustrated using a curated network of protein interactions and multiprotein complexes in baker’s yeast, with empirical protein concentrations.

Global asymptotic stability to a generalized Cohen-Grossberg BAM neural networks of neutral type delays.

PubMed

Zhang, Zhengqiu; Liu, Wenbin; Zhou, Dongming

2012-01-01

In this paper, we first discuss the existence of a unique equilibrium point of a generalized Cohen-Grossberg BAM neural networks of neutral type delays by means of the Homeomorphism theory and inequality technique. Then, by applying the existence result of an equilibrium point and constructing a Lyapunov functional, we study the global asymptotic stability of the equilibrium solution to the above Cohen-Grossberg BAM neural networks of neutral type. In our results, the hypothesis for boundedness in the existing paper, which discussed Cohen-Grossberg neural networks of neutral type on the activation functions, are removed. Finally, we give an example to demonstrate the validity of our global asymptotic stability result for the above neural networks. Copyright © 2011 Elsevier Ltd. All rights reserved.

Competing dynamic phases of active polymer networks

NASA Astrophysics Data System (ADS)

Freedman, Simon; Banerjee, Shiladitya; Dinner, Aaron R.

Recent experiments on in-vitro reconstituted assemblies of F-actin, myosin-II motors, and cross-linking proteins show that tuning local network properties can changes the fundamental biomechanical behavior of the system. For example, by varying cross-linker density and actin bundle rigidity, one can switch between contractile networks useful for reshaping cells, polarity sorted networks ideal for directed molecular transport, and frustrated networks with robust structural properties. To efficiently investigate the dynamic phases of actomyosin networks, we developed a coarse grained non-equilibrium molecular dynamics simulation of model semiflexible filaments, molecular motors, and cross-linkers with phenomenologically defined interactions. The simulation's accuracy was verified by benchmarking the mechanical properties of its individual components and collective behavior against experimental results at the molecular and network scales. By adjusting the model's parameters, we can reproduce the qualitative phases observed in experiment and predict the protein characteristics where phase crossovers could occur in collective network dynamics. Our model provides a framework for understanding cells' multiple uses of actomyosin networks and their applicability in materials research. Supported by the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program.

Networks of Learning

NASA Astrophysics Data System (ADS)

Bettencourt, Luis; Kaiser, David

2004-03-01

Based on an a historically documented example of scientific discovery - Feynman diagrams as the main calculational tool of theoretical high energy Physics - we map the time evolution of the social network of early adopters through in the US, UK, Japan and the USSR. The spread of the technique for total number of users in each region is then modelled in terms of epidemic models, highlighting parallel and divergent aspects of this analogy. We also show that transient social arrangements develop as the idea is introduced and learned, which later disappear as the technique becomes common knowledge. Such early transient is characterized by abnormally low connectivity distribution powers and by high clustering. This interesting early non-equilibrium stage of network evolution is captured by a new dynamical model for network evolution, which coincides in its long time limit with familiar preferential aggregation dynamics.

Molecular Mechanisms of Stress-Responsive Changes in Collagen and Elastin Networks in Skin.

PubMed

Aziz, Jazli; Shezali, Hafiz; Radzi, Zamri; Yahya, Noor Azlin; Abu Kassim, Noor Hayaty; Czernuszka, Jan; Rahman, Mohammad Tariqur

2016-01-01

Collagen and elastin networks make up the majority of the extracellular matrix in many organs, such as the skin. The mechanisms which are involved in the maintenance of homeostatic equilibrium of these networks are numerous, involving the regulation of genetic expression, growth factor secretion, signalling pathways, secondary messaging systems, and ion channel activity. However, many factors are capable of disrupting these pathways, which leads to an imbalance of homeostatic equilibrium. Ultimately, this leads to changes in the physical nature of skin, both functionally and cosmetically. Although various factors have been identified, including carcinogenesis, ultraviolet exposure, and mechanical stretching of skin, it was discovered that many of them affect similar components of regulatory pathways, such as fibroblasts, lysyl oxidase, and fibronectin. Additionally, it was discovered that the various regulatory pathways intersect with each other at various stages instead of working independently of each other. This review paper proposes a model which elucidates how these molecular pathways intersect with one another, and how various internal and external factors can disrupt these pathways, ultimately leading to a disruption in collagen and elastin networks. © 2016 S. Karger AG, Basel.

Complete stability of delayed recurrent neural networks with Gaussian activation functions.

PubMed

Liu, Peng; Zeng, Zhigang; Wang, Jun

2017-01-01

This paper addresses the complete stability of delayed recurrent neural networks with Gaussian activation functions. By means of the geometrical properties of Gaussian function and algebraic properties of nonsingular M-matrix, some sufficient conditions are obtained to ensure that for an n-neuron neural network, there are exactly 3 k equilibrium points with 0≤k≤n, among which 2 k and 3 k -2 k equilibrium points are locally exponentially stable and unstable, respectively. Moreover, it concludes that all the states converge to one of the equilibrium points; i.e., the neural networks are completely stable. The derived conditions herein can be easily tested. Finally, a numerical example is given to illustrate the theoretical results. Copyright © 2016 Elsevier Ltd. All rights reserved.

Financial Structure and Economic Welfare: Applied General Equilibrium Development Economics.

PubMed

Townsend, Robert

2010-09-01

This review provides a common framework for researchers thinking about the next generation of micro-founded macro models of growth, inequality, and financial deepening, as well as direction for policy makers targeting microfinance programs to alleviate poverty. Topics include treatment of financial structure general equilibrium models: testing for as-if-complete markets or other financial underpinnings; examining dual-sector models with both a perfectly intermediated sector and a sector in financial autarky, as well as a second generation of these models that embeds information problems and other obstacles to trade; designing surveys to capture measures of income, investment/savings, and flow of funds; and aggregating individuals and households to the level of network, village, or national economy. The review concludes with new directions that overcome conceptual and computational limitations.

Financial Structure and Economic Welfare: Applied General Equilibrium Development Economics

PubMed Central

Townsend, Robert

2010-01-01

This review provides a common framework for researchers thinking about the next generation of micro-founded macro models of growth, inequality, and financial deepening, as well as direction for policy makers targeting microfinance programs to alleviate poverty. Topics include treatment of financial structure general equilibrium models: testing for as-if-complete markets or other financial underpinnings; examining dual-sector models with both a perfectly intermediated sector and a sector in financial autarky, as well as a second generation of these models that embeds information problems and other obstacles to trade; designing surveys to capture measures of income, investment/savings, and flow of funds; and aggregating individuals and households to the level of network, village, or national economy. The review concludes with new directions that overcome conceptual and computational limitations. PMID:21037939

Modeling Misbehavior in Cooperative Diversity: A Dynamic Game Approach

NASA Astrophysics Data System (ADS)

Dehnie, Sintayehu; Memon, Nasir

2009-12-01

Cooperative diversity protocols are designed with the assumption that terminals always help each other in a socially efficient manner. This assumption may not be valid in commercial wireless networks where terminals may misbehave for selfish or malicious intentions. The presence of misbehaving terminals creates a social-dilemma where terminals exhibit uncertainty about the cooperative behavior of other terminals in the network. Cooperation in social-dilemma is characterized by a suboptimal Nash equilibrium where wireless terminals opt out of cooperation. Hence, without establishing a mechanism to detect and mitigate effects of misbehavior, it is difficult to maintain a socially optimal cooperation. In this paper, we first examine effects of misbehavior assuming static game model and show that cooperation under existing cooperative protocols is characterized by a noncooperative Nash equilibrium. Using evolutionary game dynamics we show that a small number of mutants can successfully invade a population of cooperators, which indicates that misbehavior is an evolutionary stable strategy (ESS). Our main goal is to design a mechanism that would enable wireless terminals to select reliable partners in the presence of uncertainty. To this end, we formulate cooperative diversity as a dynamic game with incomplete information. We show that the proposed dynamic game formulation satisfied the conditions for the existence of perfect Bayesian equilibrium.

Tensegrity and motor-driven effective interactions in a model cytoskeleton

NASA Astrophysics Data System (ADS)

Wang, Shenshen; Wolynes, Peter G.

2012-04-01

Actomyosin networks are major structural components of the cell. They provide mechanical integrity and allow dynamic remodeling of eukaryotic cells, self-organizing into the diverse patterns essential for development. We provide a theoretical framework to investigate the intricate interplay between local force generation, network connectivity, and collective action of molecular motors. This framework is capable of accommodating both regular and heterogeneous pattern formation, arrested coarsening and macroscopic contraction in a unified manner. We model the actomyosin system as a motorized cat's cradle consisting of a crosslinked network of nonlinear elastic filaments subjected to spatially anti-correlated motor kicks acting on motorized (fibril) crosslinks. The phase diagram suggests there can be arrested phase separation which provides a natural explanation for the aggregation and coalescence of actomyosin condensates. Simulation studies confirm the theoretical picture that a nonequilibrium many-body system driven by correlated motor kicks can behave as if it were at an effective equilibrium, but with modified interactions that account for the correlation of the motor driven motions of the actively bonded nodes. Regular aster patterns are observed both in Brownian dynamics simulations at effective equilibrium and in the complete stochastic simulations. The results show that large-scale contraction requires correlated kicking.

Active tension network model suggests an exotic mechanical state realized in epithelial tissues

NASA Astrophysics Data System (ADS)

Noll, Nicholas; Mani, Madhav; Heemskerk, Idse; Streichan, Sebastian J.; Shraiman, Boris I.

2017-12-01

Mechanical interactions play a crucial role in epithelial morphogenesis, yet understanding the complex mechanisms through which stress and deformation affect cell behaviour remains an open problem. Here we formulate and analyse the active tension network (ATN) model, which assumes that the mechanical balance of cells within a tissue is dominated by cortical tension and introduces tension-dependent active remodelling of the cortex. We find that ATNs exhibit unusual mechanical properties. Specifically, an ATN behaves as a fluid at short times, but at long times supports external tension like a solid. Furthermore, an ATN has an extensively degenerate equilibrium mechanical state associated with a discrete conformal--`isogonal'--deformation of cells. The ATN model predicts a constraint on equilibrium cell geometries, which we demonstrate to approximately hold in certain epithelial tissues. We further show that isogonal modes are observed in the fruit fly embryo, accounting for the striking variability of apical areas of ventral cells and helping understand the early phase of gastrulation. Living matter realizes new and exotic mechanical states, the study of which helps to understand biological phenomena.

Real-tiem Adaptive Control Scheme for Superior Plasma Confinement

DOE Office of Scientific and Technical Information (OSTI.GOV)

Alexander Trunov, Ph.D.

2001-06-01

During this Phase I project, IOS, in collaboration with our subcontractors at General Atomics, Inc., acquired and analyzed measurement data on various plasma equilibrium modes. We developed a Matlab-based toolbox consisting of linear and neural network approximators that are capable of learning and predicting, with accuracy, the behavior of plasma parameters. We also began development of the control algorithm capable of using the model of the plasma obtained by the neural network approximator.

Heterogeneous network promotes species coexistence: metapopulation model for rock-paper-scissors game.

PubMed

Nagatani, Takashi; Ichinose, Genki; Tainaka, Kei-Ichi

2018-05-04

Understanding mechanisms of biodiversity has been a central question in ecology. The coexistence of three species in rock-paper-scissors (RPS) systems are discussed by many authors; however, the relation between coexistence and network structure is rarely discussed. Here we present a metapopulation model for RPS game. The total population is assumed to consist of three subpopulations (nodes). Each individual migrates by random walk; the destination of migration is randomly determined. From reaction-migration equations, we obtain the population dynamics. It is found that the dynamic highly depends on network structures. When a network is homogeneous, the dynamics are neutrally stable: each node has a periodic solution, and the oscillations synchronize in all nodes. However, when a network is heterogeneous, the dynamics approach stable focus and all nodes reach equilibriums with different densities. Hence, the heterogeneity of the network promotes biodiversity.

Energetics in a model of prebiotic evolution

NASA Astrophysics Data System (ADS)

Intoy, B. F.; Halley, J. W.

2017-12-01

Previously we reported [A. Wynveen et al., Phys. Rev. E 89, 022725 (2014), 10.1103/PhysRevE.89.022725] that requiring that the systems regarded as lifelike be out of chemical equilibrium in a model of abstracted polymers undergoing ligation and scission first introduced by Kauffman [S. A. Kauffman, The Origins of Order (Oxford University Press, New York, 1993), Chap. 7] implied that lifelike systems were most probable when the reaction network was sparse. The model was entirely statistical and took no account of the bond energies or other energetic constraints. Here we report results of an extension of the model to include effects of a finite bonding energy in the model. We studied two conditions: (1) A food set is continuously replenished and the total polymer population is constrained but the system is otherwise isolated and (2) in addition to the constraints in (1) the system is in contact with a finite-temperature heat bath. In each case, detailed balance in the dynamics is guaranteed during the computations by continuous recomputation of a temperature [in case (1)] and of the chemical potential (in both cases) toward which the system is driven by the dynamics. In the isolated case, the probability of reaching a metastable nonequilibrium state in this model depends significantly on the composition of the food set, and the nonequilibrium states satisfying lifelike condition turn out to be at energies and particle numbers consistent with an equilibrium state at high negative temperature. As a function of the sparseness of the reaction network, the lifelike probability is nonmonotonic, as in our previous model, but the maximum probability occurs when the network is less sparse. In the case of contact with a thermal bath at a positive ambient temperature, we identify two types of metastable nonequilibrium states, termed locally and thermally alive, and locally dead and thermally alive, and evaluate their likelihood of appearance, finding maxima at an optimal temperature and an optimal degree of sparseness in the network. We use a Euclidean metric in the space of polymer populations to distinguish these states from one another and from fully equilibrated states. The metric can be used to characterize the degree and type of chemical equilibrium in observed systems, as we illustrate for the proteome of the ribosome.

Equilibrium paths analysis of materials with rheological properties by using the chaos theory

NASA Astrophysics Data System (ADS)

Bednarek, Paweł; Rządkowski, Jan

2018-01-01

The numerical equilibrium path analysis of the material with random rheological properties by using standard procedures and specialist computer programs was not successful. The proper solution for the analysed heuristic model of the material was obtained on the base of chaos theory elements and neural networks. The paper deals with mathematical reasons of used computer programs and also are elaborated the properties of the attractor used in analysis. There are presented results of conducted numerical analysis both in a numerical and in graphical form for the used procedures.

Opinion diversity and community formation in adaptive networks

NASA Astrophysics Data System (ADS)

Yu, Y.; Xiao, G.; Li, G.; Tay, W. P.; Teoh, H. F.

2017-10-01

It is interesting and of significant importance to investigate how network structures co-evolve with opinions. In this article, we show that, a simple model integrating consensus formation, link rewiring, and opinion change allows complex system dynamics to emerge, driving the system into a dynamic equilibrium with the co-existence of diversified opinions. Specifically, similar opinion holders may form into communities yet with no strict community consensus; and rather than being separated into disconnected communities, different communities are connected by a non-trivial proportion of inter-community links. More importantly, we show that the complex dynamics may lead to different numbers of communities at the steady state with a given tolerance between different opinion holders. We construct a framework for theoretically analyzing the co-evolution process. Theoretical analysis and extensive simulation results reveal some useful insights into the complex co-evolution process, including the formation of dynamic equilibrium, the transition between different steady states with different numbers of communities, and the dynamics between opinion distribution and network modularity.

Analysis and Modeling of DIII-D Experiments With OMFIT and Neural Networks

NASA Astrophysics Data System (ADS)

Meneghini, O.; Luna, C.; Smith, S. P.; Lao, L. L.; GA Theory Team

2013-10-01

The OMFIT integrated modeling framework is designed to facilitate experimental data analysis and enable integrated simulations. This talk introduces this framework and presents a selection of its applications to the DIII-D experiment. Examples include kinetic equilibrium reconstruction analysis; evaluation of MHD stability in the core and in the edge; and self-consistent predictive steady-state transport modeling. The OMFIT framework also provides the platform for an innovative approach based on neural networks to predict electron and ion energy fluxes. In our study a multi-layer feed-forward back-propagation neural network is built and trained over a database of DIII-D data. It is found that given the same parameters that the highest fidelity models use, the neural network model is able to predict to a large degree the heat transport profiles observed in the DIII-D experiments. Once the network is built, the numerical cost of evaluating the transport coefficients is virtually nonexistent, thus making the neural network model particularly well suited for plasma control and quick exploration of operational scenarios. The implementation of the neural network model and benchmark with experimental results and gyro-kinetic models will be discussed. Work supported in part by the US DOE under DE-FG02-95ER54309.

Epidemic spreading on adaptively weighted scale-free networks.

PubMed

Sun, Mengfeng; Zhang, Haifeng; Kang, Huiyan; Zhu, Guanghu; Fu, Xinchu

2017-04-01

We introduce three modified SIS models on scale-free networks that take into account variable population size, nonlinear infectivity, adaptive weights, behavior inertia and time delay, so as to better characterize the actual spread of epidemics. We develop new mathematical methods and techniques to study the dynamics of the models, including the basic reproduction number, and the global asymptotic stability of the disease-free and endemic equilibria. We show the disease-free equilibrium cannot undergo a Hopf bifurcation. We further analyze the effects of local information of diseases and various immunization schemes on epidemic dynamics. We also perform some stochastic network simulations which yield quantitative agreement with the deterministic mean-field approach.

Degradable transportation network with the addition of electric vehicles: Network equilibrium analysis

PubMed Central

Zhang, Rui; Yao, Enjian; Yang, Yang

2017-01-01

Introducing electric vehicles (EVs) into urban transportation network brings higher requirement on travel time reliability and charging reliability. Specifically, it is believed that travel time reliability is a key factor influencing travelers’ route choice. Meanwhile, due to the limited cruising range, EV drivers need to better learn about the required energy for the whole trip to make decisions about whether charging or not and where to charge (i.e., charging reliability). Since EV energy consumption is highly related to travel speed, network uncertainty affects travel time and charging demand estimation significantly. Considering the network uncertainty resulted from link degradation, which influences the distribution of travel demand on transportation network and the energy demand on power network, this paper aims to develop a reliability-based network equilibrium framework for accommodating degradable road conditions with the addition of EVs. First, based on the link travel time distribution, the mean and variance of route travel time and monetary expenses related to energy consumption are deduced, respectively. And the charging time distribution of EVs with charging demand is also estimated. Then, a nested structure is considered to deal with the difference of route choice behavior derived by the different uncertainty degrees between the routes with and without degradable links. Given the expected generalized travel cost and a psychological safety margin, a traffic assignment model with the addition of EVs is formulated. Subsequently, a heuristic solution algorithm is developed to solve the proposed model. Finally, the effects of travelers’ risk attitude, network degradation degree, and EV penetration rate on network performance are illustrated through an example network. The numerical results show that the difference of travelers’ risk attitudes does have impact on the route choice, and the widespread adoption of EVs can cut down the total system travel cost effectively when the transportation network is more reliable. PMID:28886167

Degradable transportation network with the addition of electric vehicles: Network equilibrium analysis.

PubMed

Zhang, Rui; Yao, Enjian; Yang, Yang

2017-01-01

Introducing electric vehicles (EVs) into urban transportation network brings higher requirement on travel time reliability and charging reliability. Specifically, it is believed that travel time reliability is a key factor influencing travelers' route choice. Meanwhile, due to the limited cruising range, EV drivers need to better learn about the required energy for the whole trip to make decisions about whether charging or not and where to charge (i.e., charging reliability). Since EV energy consumption is highly related to travel speed, network uncertainty affects travel time and charging demand estimation significantly. Considering the network uncertainty resulted from link degradation, which influences the distribution of travel demand on transportation network and the energy demand on power network, this paper aims to develop a reliability-based network equilibrium framework for accommodating degradable road conditions with the addition of EVs. First, based on the link travel time distribution, the mean and variance of route travel time and monetary expenses related to energy consumption are deduced, respectively. And the charging time distribution of EVs with charging demand is also estimated. Then, a nested structure is considered to deal with the difference of route choice behavior derived by the different uncertainty degrees between the routes with and without degradable links. Given the expected generalized travel cost and a psychological safety margin, a traffic assignment model with the addition of EVs is formulated. Subsequently, a heuristic solution algorithm is developed to solve the proposed model. Finally, the effects of travelers' risk attitude, network degradation degree, and EV penetration rate on network performance are illustrated through an example network. The numerical results show that the difference of travelers' risk attitudes does have impact on the route choice, and the widespread adoption of EVs can cut down the total system travel cost effectively when the transportation network is more reliable.

Self-organization of network dynamics into local quantized states.

PubMed

Nicolaides, Christos; Juanes, Ruben; Cueto-Felgueroso, Luis

2016-02-17

Self-organization and pattern formation in network-organized systems emerges from the collective activation and interaction of many interconnected units. A striking feature of these non-equilibrium structures is that they are often localized and robust: only a small subset of the nodes, or cell assembly, is activated. Understanding the role of cell assemblies as basic functional units in neural networks and socio-technical systems emerges as a fundamental challenge in network theory. A key open question is how these elementary building blocks emerge, and how they operate, linking structure and function in complex networks. Here we show that a network analogue of the Swift-Hohenberg continuum model-a minimal-ingredients model of nodal activation and interaction within a complex network-is able to produce a complex suite of localized patterns. Hence, the spontaneous formation of robust operational cell assemblies in complex networks can be explained as the result of self-organization, even in the absence of synaptic reinforcements.

Rumor diffusion in heterogeneous networks by considering the individuals' subjective judgment and diverse characteristics

NASA Astrophysics Data System (ADS)

Ma, Jing; Zhu, He

2018-06-01

In this study, we propose a novel rumor spreading model in consideration of the individuals' subjective judgment and diverse characteristics. To reflect the diversity of the individuals' characteristics, we introduce two probability distribution functions, which could be chosen arbitrarily or given by empirical data, to characterize individuals' mastering degree of knowledge with respect to the domain of a specific rumor and individuals' rationality degree. Different from existing models, no two persons in our model are identical, and each individual can judge the authenticity of the information, e.g., rumors, with his distinctive characteristics. In addition, by means of the mean-field method, we establish the expression of the dynamics of the rumor propagation in the complex heterogeneous networks and derive the rumor spreading threshold. Through the theoretical analysis, we find that the threshold is independent of the forms of the two introduced functions. Furthermore, we prove the stability of the rumor-free equilibrium set E0. That is if and only if R0 < 1, the rumor-free equilibrium set E0 is globally asymptotically stable. Finally, we conduct a series of numerical simulations to verify the theoretical results and comprehensively illustrate the evolution of the model. The simulation results show that because of the diversity of individuals' characteristics, it becomes more difficult for the rumor to disseminate in the networks and the higher the mean of knowledge and the mean of rationality are, the more time it will take for the model to evolve to the steady state.

A recurrent neural network for solving bilevel linear programming problem.

PubMed

He, Xing; Li, Chuandong; Huang, Tingwen; Li, Chaojie; Huang, Junjian

2014-04-01

In this brief, based on the method of penalty functions, a recurrent neural network (NN) modeled by means of a differential inclusion is proposed for solving the bilevel linear programming problem (BLPP). Compared with the existing NNs for BLPP, the model has the least number of state variables and simple structure. Using nonsmooth analysis, the theory of differential inclusions, and Lyapunov-like method, the equilibrium point sequence of the proposed NNs can approximately converge to an optimal solution of BLPP under certain conditions. Finally, the numerical simulations of a supply chain distribution model have shown excellent performance of the proposed recurrent NNs.

Multi-period equilibrium/near-equilibrium in electricity markets based on locational marginal prices

NASA Astrophysics Data System (ADS)

Garcia Bertrand, Raquel

In this dissertation we propose an equilibrium procedure that coordinates the point of view of every market agent resulting in an equilibrium that simultaneously maximizes the independent objective of every market agent and satisfies network constraints. Therefore, the activities of the generating companies, consumers and an independent system operator are modeled: (1) The generating companies seek to maximize profits by specifying hourly step functions of productions and minimum selling prices, and bounds on productions. (2) The goals of the consumers are to maximize their economic utilities by specifying hourly step functions of demands and maximum buying prices, and bounds on demands. (3) The independent system operator then clears the market taking into account consistency conditions as well as capacity and line losses so as to achieve maximum social welfare. Then, we approach this equilibrium problem using complementarity theory in order to have the capability of imposing constraints on dual variables, i.e., on prices, such as minimum profit conditions for the generating units or maximum cost conditions for the consumers. In this way, given the form of the individual optimization problems, the Karush-Kuhn-Tucker conditions for the generating companies, the consumers and the independent system operator are both necessary and sufficient. The simultaneous solution to all these conditions constitutes a mixed linear complementarity problem. We include minimum profit constraints imposed by the units in the market equilibrium model. These constraints are added as additional constraints to the equivalent quadratic programming problem of the mixed linear complementarity problem previously described. For the sake of clarity, the proposed equilibrium or near-equilibrium is first developed for the particular case considering only one time period. Afterwards, we consider an equilibrium or near-equilibrium applied to a multi-period framework. This model embodies binary decisions, i.e., on/off status for the units, and therefore optimality conditions cannot be directly applied. To avoid limitations provoked by binary variables, while retaining the advantages of using optimality conditions, we define the multi-period market equilibrium using Benders decomposition, which allows computing binary variables through the master problem and continuous variables through the subproblem. Finally, we illustrate these market equilibrium concepts through several case studies.

A 3-states magnetic model of binary decisions in sociophysics

NASA Astrophysics Data System (ADS)

Fernandez, Miguel A.; Korutcheva, Elka; de la Rubia, F. Javier

2016-11-01

We study a diluted Blume-Capel model of 3-states sites as an attempt to understand how some social processes as cooperation or organization happen. For this aim, we study the effect of the complex network topology on the equilibrium properties of the model, by focusing on three different substrates: random graph, Watts-Strogatz and Newman substrates. Our computer simulations are in good agreement with the corresponding analytical results.

Birth and death of protein domains: A simple model of evolution explains power law behavior

PubMed Central

Karev, Georgy P; Wolf, Yuri I; Rzhetsky, Andrey Y; Berezovskaya, Faina S; Koonin, Eugene V

2002-01-01

Background Power distributions appear in numerous biological, physical and other contexts, which appear to be fundamentally different. In biology, power laws have been claimed to describe the distributions of the connections of enzymes and metabolites in metabolic networks, the number of interactions partners of a given protein, the number of members in paralogous families, and other quantities. In network analysis, power laws imply evolution of the network with preferential attachment, i.e. a greater likelihood of nodes being added to pre-existing hubs. Exploration of different types of evolutionary models in an attempt to determine which of them lead to power law distributions has the potential of revealing non-trivial aspects of genome evolution. Results A simple model of evolution of the domain composition of proteomes was developed, with the following elementary processes: i) domain birth (duplication with divergence), ii) death (inactivation and/or deletion), and iii) innovation (emergence from non-coding or non-globular sequences or acquisition via horizontal gene transfer). This formalism can be described as a birth, death and innovation model (BDIM). The formulas for equilibrium frequencies of domain families of different size and the total number of families at equilibrium are derived for a general BDIM. All asymptotics of equilibrium frequencies of domain families possible for the given type of models are found and their appearance depending on model parameters is investigated. It is proved that the power law asymptotics appears if, and only if, the model is balanced, i.e. domain duplication and deletion rates are asymptotically equal up to the second order. It is further proved that any power asymptotic with the degree not equal to -1 can appear only if the hypothesis of independence of the duplication/deletion rates on the size of a domain family is rejected. Specific cases of BDIMs, namely simple, linear, polynomial and rational models, are considered in details and the distributions of the equilibrium frequencies of domain families of different size are determined for each case. We apply the BDIM formalism to the analysis of the domain family size distributions in prokaryotic and eukaryotic proteomes and show an excellent fit between these empirical data and a particular form of the model, the second-order balanced linear BDIM. Calculation of the parameters of these models suggests surprisingly high innovation rates, comparable to the total domain birth (duplication) and elimination rates, particularly for prokaryotic genomes. Conclusions We show that a straightforward model of genome evolution, which does not explicitly include selection, is sufficient to explain the observed distributions of domain family sizes, in which power laws appear as asymptotic. However, for the model to be compatible with the data, there has to be a precise balance between domain birth, death and innovation rates, and this is likely to be maintained by selection. The developed approach is oriented at a mathematical description of evolution of domain composition of proteomes, but a simple reformulation could be applied to models of other evolving networks with preferential attachment. PMID:12379152

Birth and death of protein domains: a simple model of evolution explains power law behavior.

PubMed

Karev, Georgy P; Wolf, Yuri I; Rzhetsky, Andrey Y; Berezovskaya, Faina S; Koonin, Eugene V

2002-10-14

Power distributions appear in numerous biological, physical and other contexts, which appear to be fundamentally different. In biology, power laws have been claimed to describe the distributions of the connections of enzymes and metabolites in metabolic networks, the number of interactions partners of a given protein, the number of members in paralogous families, and other quantities. In network analysis, power laws imply evolution of the network with preferential attachment, i.e. a greater likelihood of nodes being added to pre-existing hubs. Exploration of different types of evolutionary models in an attempt to determine which of them lead to power law distributions has the potential of revealing non-trivial aspects of genome evolution. A simple model of evolution of the domain composition of proteomes was developed, with the following elementary processes: i) domain birth (duplication with divergence), ii) death (inactivation and/or deletion), and iii) innovation (emergence from non-coding or non-globular sequences or acquisition via horizontal gene transfer). This formalism can be described as a birth, death and innovation model (BDIM). The formulas for equilibrium frequencies of domain families of different size and the total number of families at equilibrium are derived for a general BDIM. All asymptotics of equilibrium frequencies of domain families possible for the given type of models are found and their appearance depending on model parameters is investigated. It is proved that the power law asymptotics appears if, and only if, the model is balanced, i.e. domain duplication and deletion rates are asymptotically equal up to the second order. It is further proved that any power asymptotic with the degree not equal to -1 can appear only if the hypothesis of independence of the duplication/deletion rates on the size of a domain family is rejected. Specific cases of BDIMs, namely simple, linear, polynomial and rational models, are considered in details and the distributions of the equilibrium frequencies of domain families of different size are determined for each case. We apply the BDIM formalism to the analysis of the domain family size distributions in prokaryotic and eukaryotic proteomes and show an excellent fit between these empirical data and a particular form of the model, the second-order balanced linear BDIM. Calculation of the parameters of these models suggests surprisingly high innovation rates, comparable to the total domain birth (duplication) and elimination rates, particularly for prokaryotic genomes. We show that a straightforward model of genome evolution, which does not explicitly include selection, is sufficient to explain the observed distributions of domain family sizes, in which power laws appear as asymptotic. However, for the model to be compatible with the data, there has to be a precise balance between domain birth, death and innovation rates, and this is likely to be maintained by selection. The developed approach is oriented at a mathematical description of evolution of domain composition of proteomes, but a simple reformulation could be applied to models of other evolving networks with preferential attachment.

Non-Lipschitzian dynamics for neural net modelling

NASA Technical Reports Server (NTRS)

Zak, Michail

1989-01-01

Failure of the Lipschitz condition in unstable equilibrium points of dynamical systems leads to a multiple-choice response to an initial deterministic input. The evolution of such systems is characterized by a special type of unpredictability measured by unbounded Liapunov exponents. Possible relation of these systems to future neural networks is discussed.

Collective Dynamics for Heterogeneous Networks of Theta Neurons

NASA Astrophysics Data System (ADS)

Luke, Tanushree

Collective behavior in neural networks has often been used as an indicator of communication between different brain areas. These collective synchronization and desynchronization patterns are also considered an important feature in understanding normal and abnormal brain function. To understand the emergence of these collective patterns, I create an analytic model that identifies all such macroscopic steady-states attainable by a network of Type-I neurons. This network, whose basic unit is the model "theta'' neuron, contains a mixture of excitable and spiking neurons coupled via a smooth pulse-like synapse. Applying the Ott-Antonsen reduction method in the thermodynamic limit, I obtain a low-dimensional evolution equation that describes the asymptotic dynamics of the macroscopic mean field of the network. This model can be used as the basis in understanding more complicated neuronal networks when additional dynamical features are included. From this reduced dynamical equation for the mean field, I show that the network exhibits three collective attracting steady-states. The first two are equilibrium states that both reflect partial synchronization in the network, whereas the third is a limit cycle in which the degree of network synchronization oscillates in time. In addition to a comprehensive identification of all possible attracting macro-states, this analytic model permits a complete bifurcation analysis of the collective behavior of the network with respect to three key network features: the degree of excitability of the neurons, the heterogeneity of the population, and the overall coupling strength. The network typically tends towards the two macroscopic equilibrium states when the neuron's intrinsic dynamics and the network interactions reinforce each other. In contrast, the limit cycle state, bifurcations, and multistability tend to occur when there is competition between these network features. I also outline here an extension of the above model where the neurons' excitability now varies in time sinuosoidally, thus simulating a parabolic bursting network. This time-varying excitability can lead to the emergence of macroscopic chaos and multistability in the collective behavior of the network. Finally, I expand the single population model described above to examine a two-population neuronal network where each population has its own unique mixture of excitable and spiking neurons, as well as its own coupling strength (either excitatory or inhibitory in nature). Specifically, I consider the situation where the first population is only allowed to influence the second population without any feedback, thus effectively creating a feed-forward "driver-response" system. In this special arrangement, the driver's asymptotic macroscopic dynamics are fully explored in the comprehensive analysis of the single population. Then, in the presence of an influence from the driver, the modified dynamics of the second population, which now acts as a response population, can also be fully analyzed. As in the time-varying model, these modifications give rise to richer dynamics to the response population than those found from the single population formalism, including multi-periodicity and chaos.

Disruption of River Networks in Nature and Models

NASA Astrophysics Data System (ADS)

Perron, J. T.; Black, B. A.; Stokes, M.; McCoy, S. W.; Goldberg, S. L.

2017-12-01

Many natural systems display especially informative behavior as they respond to perturbations. Landscapes are no exception. For example, longitudinal elevation profiles of rivers responding to changes in uplift rate can reveal differences among erosional mechanisms that are obscured while the profiles are in equilibrium. The responses of erosional river networks to perturbations, including disruption of their network structure by diversion, truncation, resurfacing, or river capture, may be equally revealing. In this presentation, we draw attention to features of disrupted erosional river networks that a general model of landscape evolution should be able to reproduce, including the consequences of different styles of planetary tectonics and the response to heterogeneous bedrock structure and deformation. A comparison of global drainage directions with long-wavelength topography on Earth, Mars, and Saturn's moon Titan reveals the extent to which persistent and relatively rapid crustal deformation has disrupted river networks on Earth. Motivated by this example and others, we ask whether current models of river network evolution adequately capture the disruption of river networks by tectonic, lithologic, or climatic perturbations. In some cases the answer appears to be no, and we suggest some processes that models may be missing.

An Adaptive QSE-reduced Nuclear Reaction Network for Silicon Burning

NASA Astrophysics Data System (ADS)

Parete-Koon, Suzanne; Hix, W.; Thielemann, F.

2008-03-01

The nuclei of the "iron peak" are formed in massive stars shortly before core collapse and during their supernova outbursts as well as during thermonuclear supernovae. Complete and incomplete silicon burning during these events are responsible for the production of a wide range of nuclei with atomic mass numbers from 28 to 64. Because of the large number of nuclei involved, accurate modeling of silicon burning is computationally expensive. However, examination of the physics of silicon burning has revealed that the nuclear evolution is dominated by large groups of nuclei in mutual equilibrium. We present an improvement on our hybrid equilibrium-network scheme which takes advantage of this quasi-equilibrium in order to reduce the number of independent variables calculated. Because the size and membership of these groups vary as the temperature, density and electron faction change, achieving maximal efficiency requires dynamic adjustment of group number and membership. Toward this end, we are implementing a scheme beginning with a single QSE (NSE) group at appropriately high temperature, then progressing through 2, 3 and 4 group stages (with successively more independent variables) as temperature declines. This combination allows accurate prediction of the nuclear abundance evolution, deleptonization and energy generation at a further reduced computational cost when compared to a conventional nuclear reaction network or our previous 3 fixed group QSE-reduced network. During silicon burning, the resultant QSE-reduced network is up to 20 times faster than the full network it replaces without significant loss of accuracy. These reductions in computational cost and the number of species evolved make QSE-reduced networks well suited for inclusion within hydrodynamic simulations, particularly in multi-dimensional applications. This work has been supported by the National Science Foundation, by the Department of Energy's Scientic Discovery through Advanced Computing Programs, and by the Joint Institute for Heavy Ion Research at ORNL.

Convergence to equilibrium of renormalised solutions to nonlinear chemical reaction–diffusion systems

NASA Astrophysics Data System (ADS)

Fellner, Klemens; Tang, Bao Quoc

2018-06-01

The convergence to equilibrium for renormalised solutions to nonlinear reaction-diffusion systems is studied. The considered reaction-diffusion systems arise from chemical reaction networks with mass action kinetics and satisfy the complex balanced condition. By applying the so-called entropy method, we show that if the system does not have boundary equilibria, i.e. equilibrium states lying on the boundary of R_+^N, then any renormalised solution converges exponentially to the complex balanced equilibrium with a rate, which can be computed explicitly up to a finite-dimensional inequality. This inequality is proven via a contradiction argument and thus not explicitly. An explicit method of proof, however, is provided for a specific application modelling a reversible enzyme reaction by exploiting the specific structure of the conservation laws. Our approach is also useful to study the trend to equilibrium for systems possessing boundary equilibria. More precisely, to show the convergence to equilibrium for systems with boundary equilibria, we establish a sufficient condition in terms of a modified finite-dimensional inequality along trajectories of the system. By assuming this condition, which roughly means that the system produces too much entropy to stay close to a boundary equilibrium for infinite time, the entropy method shows exponential convergence to equilibrium for renormalised solutions to complex balanced systems with boundary equilibria.

Adaptive Chemical Networks under Non-Equilibrium Conditions: The Evaporating Droplet.

PubMed

Armao, Joseph J; Lehn, Jean-Marie

2016-10-17

Non-volatile solutes in an evaporating drop experience an out-of-equilibrium state due to non-linear concentration effects and complex flow patterns. Here, we demonstrate a small molecule chemical reaction network that undergoes a rapid adaptation response to the out-of-equilibrium conditions inside the droplet leading to control over the molecular constitution and spatial arrangement of the deposition pattern. Adaptation results in a pronounced coffee stain effect and coupling to chemical concentration gradients within the drop is demonstrated. Amplification and suppression of network species are readily identifiable with confocal fluorescence microscopy. We anticipate that these observations will contribute to the design and exploration of out-of-equilibrium chemical systems, as well as be useful towards the development of point-of-care medical diagnostics and controlled deposition of small molecules through inkjet printing. © 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Modeling the Propagation of Mobile Phone Virus under Complex Network

PubMed Central

Yang, Wei; Wei, Xi-liang; Guo, Hao; An, Gang; Guo, Lei

2014-01-01

Mobile phone virus is a rogue program written to propagate from one phone to another, which can take control of a mobile device by exploiting its vulnerabilities. In this paper the propagation model of mobile phone virus is tackled to understand how particular factors can affect its propagation and design effective containment strategies to suppress mobile phone virus. Two different propagation models of mobile phone viruses under the complex network are proposed in this paper. One is intended to describe the propagation of user-tricking virus, and the other is to describe the propagation of the vulnerability-exploiting virus. Based on the traditional epidemic models, the characteristics of mobile phone viruses and the network topology structure are incorporated into our models. A detailed analysis is conducted to analyze the propagation models. Through analysis, the stable infection-free equilibrium point and the stability condition are derived. Finally, considering the network topology, the numerical and simulation experiments are carried out. Results indicate that both models are correct and suitable for describing the spread of two different mobile phone viruses, respectively. PMID:25133209

Dynamic stability analysis of fractional order leaky integrator echo state neural networks

NASA Astrophysics Data System (ADS)

Pahnehkolaei, Seyed Mehdi Abedi; Alfi, Alireza; Tenreiro Machado, J. A.

2017-06-01

The Leaky integrator echo state neural network (Leaky-ESN) is an improved model of the recurrent neural network (RNN) and adopts an interconnected recurrent grid of processing neurons. This paper presents a new proof for the convergence of a Lyapunov candidate function to zero when time tends to infinity by means of the Caputo fractional derivative with order lying in the range (0, 1). The stability of Fractional-Order Leaky-ESN (FO Leaky-ESN) is then analyzed, and the existence, uniqueness and stability of the equilibrium point are provided. A numerical example demonstrates the feasibility of the proposed method.

Oscillator Neural Network Retrieving Sparsely Coded Phase Patterns

NASA Astrophysics Data System (ADS)

Aoyagi, Toshio; Nomura, Masaki

1999-08-01

Little is known theoretically about the associative memory capabilities of neural networks in which information is encoded not only in the mean firing rate but also in the timing of firings. Particularly, in the case of sparsely coded patterns, it is biologically important to consider the timings of firings and to study how such consideration influences storage capacities and quality of recalled patterns. For this purpose, we propose a simple extended model of oscillator neural networks to allow for expression of a nonfiring state. Analyzing both equilibrium states and dynamical properties in recalling processes, we find that the system possesses good associative memory.

Impact of reduced scale free network on wireless sensor network

NASA Astrophysics Data System (ADS)

Keshri, Neha; Gupta, Anurag; Mishra, Bimal Kumar

2016-12-01

In heterogeneous wireless sensor network (WSN) each data-packet traverses through multiple hops over restricted communication range before it reaches the sink. The amount of energy required to transmit a data-packet is directly proportional to the number of hops. To balance the energy costs across the entire network and to enhance the robustness in order to improve the lifetime of WSN becomes a key issue of researchers. Due to high dimensionality of an epidemic model of WSN over a general scale free network, it is quite difficult to have close study of network dynamics. To overcome this complexity, we simplify a general scale free network by partitioning all of its motes into two classes: higher-degree motes and lower-degree motes, and equating the degrees of all higher-degree motes with lower-degree motes, yielding a reduced scale free network. We develop an epidemic model of WSN based on reduced scale free network. The existence of unique positive equilibrium is determined with some restrictions. Stability of the system is proved. Furthermore, simulation results show improvements made in this paper have made the entire network have a better robustness to the network failure and the balanced energy costs. This reduced model based on scale free network theory proves more applicable to the research of WSN.

Irreversible opinion spreading on scale-free networks

NASA Astrophysics Data System (ADS)

Candia, Julián

2007-02-01

We study the dynamical and critical behavior of a model for irreversible opinion spreading on Barabási-Albert (BA) scale-free networks by performing extensive Monte Carlo simulations. The opinion spreading within an inhomogeneous society is investigated by means of the magnetic Eden model, a nonequilibrium kinetic model for the growth of binary mixtures in contact with a thermal bath. The deposition dynamics, which is studied as a function of the degree of the occupied sites, shows evidence for the leading role played by hubs in the growth process. Systems of finite size grow either ordered or disordered, depending on the temperature. By means of standard finite-size scaling procedures, the effective order-disorder phase transitions are found to persist in the thermodynamic limit. This critical behavior, however, is absent in related equilibrium spin systems such as the Ising model on BA scale-free networks, which in the thermodynamic limit only displays a ferromagnetic phase. The dependence of these results on the degree exponent is also discussed for the case of uncorrelated scale-free networks.

Decentralized supply chain network design: monopoly, duopoly and oligopoly competitions under uncertainty

NASA Astrophysics Data System (ADS)

Seyedhosseini, Seyed Mohammad; Fahimi, Kaveh; Makui, Ahmad

2017-12-01

This paper presents the competitive supply chain network design problem in which n decentralized supply chains simultaneously enter the market with no existing rival chain, shape their networks and set wholesale and retail prices in competitive mode. The customer demand is elastic and price dependent, customer utility function is based on the Hoteling model and the chains produce identical or highly substitutable products. We construct a solution algorithm based on bi-level programming and possibility theory. In the proposed bi-level model, the inner part sets the prices based on simultaneous extra- and Stackleberg intra- chains competitions, and the outer part shapes the networks in cooperative competitions. Finally, we use a real-word study to discuss the effect of the different structures of the competitors on the equilibrium solution. Moreover, sensitivity analyses are conducted and managerial insights are offered.

Discretized kinetic theory on scale-free networks

NASA Astrophysics Data System (ADS)

Bertotti, Maria Letizia; Modanese, Giovanni

2016-10-01

The network of interpersonal connections is one of the possible heterogeneous factors which affect the income distribution emerging from micro-to-macro economic models. In this paper we equip our model discussed in [1, 2] with a network structure. The model is based on a system of n differential equations of the kinetic discretized-Boltzmann kind. The network structure is incorporated in a probabilistic way, through the introduction of a link density P(α) and of correlation coefficients P(β|α), which give the conditioned probability that an individual with α links is connected to one with β links. We study the properties of the equations and give analytical results concerning the existence, normalization and positivity of the solutions. For a fixed network with P(α) = c/α q , we investigate numerically the dependence of the detailed and marginal equilibrium distributions on the initial conditions and on the exponent q. Our results are compatible with those obtained from the Bouchaud-Mezard model and from agent-based simulations, and provide additional information about the dependence of the individual income on the level of connectivity.

Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model.

PubMed

Wang, Wei Bu; Liang, Yu; Zhang, Jing; Wu, Yi Dong; Du, Jian Jun; Li, Qi Ming; Zhu, Jian Zhuo; Su, Ji Guo

2018-06-22

Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.

Global stability for epidemic models on multiplex networks.

PubMed

Huang, Yu-Jhe; Juang, Jonq; Liang, Yu-Hao; Wang, Hsin-Yu

2018-05-01

In this work, we consider an epidemic model in a two-layer network in which the dynamics of susceptible-infected-susceptible process in the physical layer coexists with that of a cyclic process of unaware-aware-unaware in the virtual layer. For such multiplex network, we shall define the basic reproduction number [Formula: see text] in the virtual layer, which is similar to the basic reproduction number [Formula: see text] defined in the physical layer. We show analytically that if [Formula: see text] and [Formula: see text], then the disease and information free equilibrium is globally stable and if [Formula: see text] and [Formula: see text], then the disease free and information saturated equilibrium is globally stable for all initial conditions except at the origin. In the case of [Formula: see text], whether the disease dies out or not depends on the competition between how well the information is transmitted in the virtual layer and how contagious the disease is in the physical layer. In particular, it is numerically demonstrated that if the difference in [Formula: see text] and [Formula: see text] is greater than the product of [Formula: see text], the deviation of [Formula: see text] from 1 and the relative infection rate for an aware susceptible individual, then the disease dies out. Otherwise, the disease breaks out.

Thermal non-equilibrium in porous medium adjacent to vertical plate: ANN approach

NASA Astrophysics Data System (ADS)

Ahmed, N. J. Salman; Ahamed, K. S. Nazim; Al-Rashed, Abdullah A. A. A.; Kamangar, Sarfaraz; Athani, Abdulgaphur

2018-05-01

Thermal non-equilibrium in porous medium is a condition that refers to temperature discrepancy in solid matrix and fluid of porous medium. This type of flow is complex flow requiring complex set of partial differential equations that govern the flow behavior. The current work is undertaken to predict the thermal non-equilibrium behavior of porous medium adjacent to vertical plate using artificial neural network. A set of neurons in 3 layers are trained to predict the heat transfer characteristics. It is found that the thermal non-equilibrium heat transfer behavior in terms of Nusselt number of fluid as well as solid phase can be predicted accurately by using well-trained neural network.

Optimal Stabilization of Social Welfare under Small Variation of Operating Condition with Bifurcation Analysis

NASA Astrophysics Data System (ADS)

Chanda, Sandip; De, Abhinandan

2016-12-01

A social welfare optimization technique has been proposed in this paper with a developed state space based model and bifurcation analysis to offer substantial stability margin even in most inadvertent states of power system networks. The restoration of the power market dynamic price equilibrium has been negotiated in this paper, by forming Jacobian of the sensitivity matrix to regulate the state variables for the standardization of the quality of solution in worst possible contingencies of the network and even with co-option of intermittent renewable energy sources. The model has been tested in IEEE 30 bus system and illustrious particle swarm optimization has assisted the fusion of the proposed model and methodology.

Continuous system modeling

NASA Technical Reports Server (NTRS)

Cellier, Francois E.

1991-01-01

A comprehensive and systematic introduction is presented for the concepts associated with 'modeling', involving the transition from a physical system down to an abstract description of that system in the form of a set of differential and/or difference equations, and basing its treatment of modeling on the mathematics of dynamical systems. Attention is given to the principles of passive electrical circuit modeling, planar mechanical systems modeling, hierarchical modular modeling of continuous systems, and bond-graph modeling. Also discussed are modeling in equilibrium thermodynamics, population dynamics, and system dynamics, inductive reasoning, artificial neural networks, and automated model synthesis.

Spatial-temporal modeling of malware propagation in networks.

PubMed

Chen, Zesheng; Ji, Chuanyi

2005-09-01

Network security is an important task of network management. One threat to network security is malware (malicious software) propagation. One type of malware is called topological scanning that spreads based on topology information. The focus of this work is on modeling the spread of topological malwares, which is important for understanding their potential damages, and for developing countermeasures to protect the network infrastructure. Our model is motivated by probabilistic graphs, which have been widely investigated in machine learning. We first use a graphical representation to abstract the propagation of malwares that employ different scanning methods. We then use a spatial-temporal random process to describe the statistical dependence of malware propagation in arbitrary topologies. As the spatial dependence is particularly difficult to characterize, the problem becomes how to use simple (i.e., biased) models to approximate the spatially dependent process. In particular, we propose the independent model and the Markov model as simple approximations. We conduct both theoretical analysis and extensive simulations on large networks using both real measurements and synthesized topologies to test the performance of the proposed models. Our results show that the independent model can capture temporal dependence and detailed topology information and, thus, outperforms the previous models, whereas the Markov model incorporates a certain spatial dependence and, thus, achieves a greater accuracy in characterizing both transient and equilibrium behaviors of malware propagation.

Open quantum generalisation of Hopfield neural networks

NASA Astrophysics Data System (ADS)

Rotondo, P.; Marcuzzi, M.; Garrahan, J. P.; Lesanovsky, I.; Müller, M.

2018-03-01

We propose a new framework to understand how quantum effects may impact on the dynamics of neural networks. We implement the dynamics of neural networks in terms of Markovian open quantum systems, which allows us to treat thermal and quantum coherent effects on the same footing. In particular, we propose an open quantum generalisation of the Hopfield neural network, the simplest toy model of associative memory. We determine its phase diagram and show that quantum fluctuations give rise to a qualitatively new non-equilibrium phase. This novel phase is characterised by limit cycles corresponding to high-dimensional stationary manifolds that may be regarded as a generalisation of storage patterns to the quantum domain.

A distributed predictive control approach for periodic flow-based networks: application to drinking water systems

NASA Astrophysics Data System (ADS)

Grosso, Juan M.; Ocampo-Martinez, Carlos; Puig, Vicenç

2017-10-01

This paper proposes a distributed model predictive control approach designed to work in a cooperative manner for controlling flow-based networks showing periodic behaviours. Under this distributed approach, local controllers cooperate in order to enhance the performance of the whole flow network avoiding the use of a coordination layer. Alternatively, controllers use both the monolithic model of the network and the given global cost function to optimise the control inputs of the local controllers but taking into account the effect of their decisions over the remainder subsystems conforming the entire network. In this sense, a global (all-to-all) communication strategy is considered. Although the Pareto optimality cannot be reached due to the existence of non-sparse coupling constraints, the asymptotic convergence to a Nash equilibrium is guaranteed. The resultant strategy is tested and its effectiveness is shown when applied to a large-scale complex flow-based network: the Barcelona drinking water supply system.

Evolving dynamics of trading behavior based on coordination game in complex networks

NASA Astrophysics Data System (ADS)

Bian, Yue-tang; Xu, Lu; Li, Jin-sheng

2016-05-01

This work concerns the modeling of evolvement of trading behavior in stock markets. Based on the assumption of the investors' limited rationality, the evolution mechanism of trading behavior is modeled according to the investment strategy of coordination game in network, that investors are prone to imitate their neighbors' activity through comprehensive analysis on the risk dominance degree of certain investment behavior, the network topology of their relationship and its heterogeneity. We investigate by mean-field analysis and extensive simulations the evolution of investors' trading behavior in various typical networks under different risk dominance degree of investment behavior. Our results indicate that the evolution of investors' behavior is affected by the network structure of stock market and the effect of risk dominance degree of investment behavior; the stability of equilibrium states of investors' behavior dynamics is directly related with the risk dominance degree of some behavior; connectivity and heterogeneity of the network plays an important role in the evolution of the investment behavior in stock market.

Construction of road network vulnerability evaluation index based on general travel cost

NASA Astrophysics Data System (ADS)

Leng, Jun-qiang; Zhai, Jing; Li, Qian-wen; Zhao, Lin

2018-03-01

With the development of China's economy and the continuous improvement of her urban road network, the vulnerability of the urban road network has attracted increasing attention. Based on general travel cost, this work constructs the vulnerability evaluation index for the urban road network, and evaluates the vulnerability of the urban road network from the perspective of user generalised travel cost. Firstly, the generalised travel cost model is constructed based on vehicle cost, travel time, and traveller comfort. Then, the network efficiency index is selected as an evaluation index of vulnerability: the network efficiency index is composed of the traffic volume and the generalised travel cost, which are obtained from the equilibrium state of the network. In addition, the research analyses the influence of traffic capacity decrease, road section attribute value, and location of road section, on vulnerability. Finally, the vulnerability index is used to analyse the local area network of Harbin and verify its applicability.

Final Technical Report: Mathematical Foundations for Uncertainty Quantification in Materials Design

DOE Office of Scientific and Technical Information (OSTI.GOV)

Plechac, Petr; Vlachos, Dionisios G.

We developed path-wise information theory-based and goal-oriented sensitivity analysis and parameter identification methods for complex high-dimensional dynamics and in particular of non-equilibrium extended molecular systems. The combination of these novel methodologies provided the first methods in the literature which are capable to handle UQ questions for stochastic complex systems with some or all of the following features: (a) multi-scale stochastic models such as (bio)chemical reaction networks, with a very large number of parameters, (b) spatially distributed systems such as Kinetic Monte Carlo or Langevin Dynamics, (c) non-equilibrium processes typically associated with coupled physico-chemical mechanisms, driven boundary conditions, hybrid micro-macro systems,more » etc. A particular computational challenge arises in simulations of multi-scale reaction networks and molecular systems. Mathematical techniques were applied to in silico prediction of novel materials with emphasis on the effect of microstructure on model uncertainty quantification (UQ). We outline acceleration methods to make calculations of real chemistry feasible followed by two complementary tasks on structure optimization and microstructure-induced UQ.« less

Economic Feasibility of Wireless Sensor Network-Based Service Provision in a Duopoly Setting with a Monopolist Operator.

PubMed

Sanchis-Cano, Angel; Romero, Julián; Sacoto-Cabrera, Erwin J; Guijarro, Luis

2017-11-25

We analyze the feasibility of providing Wireless Sensor Network-data-based services in an Internet of Things scenario from an economical point of view. The scenario has two competing service providers with their own private sensor networks, a network operator and final users. The scenario is analyzed as two games using game theory. In the first game, sensors decide to subscribe or not to the network operator to upload the collected sensing-data, based on a utility function related to the mean service time and the price charged by the operator. In the second game, users decide to subscribe or not to the sensor-data-based service of the service providers based on a Logit discrete choice model related to the quality of the data collected and the subscription price. The sinks and users subscription stages are analyzed using population games and discrete choice models, while network operator and service providers pricing stages are analyzed using optimization and Nash equilibrium concepts respectively. The model is shown feasible from an economic point of view for all the actors if there are enough interested final users and opens the possibility of developing more efficient models with different types of services.

Boolean decision problems with competing interactions on scale-free networks: Equilibrium and nonequilibrium behavior in an external bias

NASA Astrophysics Data System (ADS)

Zhu, Zheng; Andresen, Juan Carlos; Moore, M. A.; Katzgraber, Helmut G.

2014-02-01

We study the equilibrium and nonequilibrium properties of Boolean decision problems with competing interactions on scale-free networks in an external bias (magnetic field). Previous studies at zero field have shown a remarkable equilibrium stability of Boolean variables (Ising spins) with competing interactions (spin glasses) on scale-free networks. When the exponent that describes the power-law decay of the connectivity of the network is strictly larger than 3, the system undergoes a spin-glass transition. However, when the exponent is equal to or less than 3, the glass phase is stable for all temperatures. First, we perform finite-temperature Monte Carlo simulations in a field to test the robustness of the spin-glass phase and show that the system has a spin-glass phase in a field, i.e., exhibits a de Almeida-Thouless line. Furthermore, we study avalanche distributions when the system is driven by a field at zero temperature to test if the system displays self-organized criticality. Numerical results suggest that avalanches (damage) can spread across the whole system with nonzero probability when the decay exponent of the interaction degree is less than or equal to 2, i.e., that Boolean decision problems on scale-free networks with competing interactions can be fragile when not in thermal equilibrium.

Dynamics analysis of epidemic and information spreading in overlay networks.

PubMed

Liu, Guirong; Liu, Zhimei; Jin, Zhen

2018-05-07

We establish an SIS-UAU model to present the dynamics of epidemic and information spreading in overlay networks. The overlay network is represented by two layers: one where the dynamics of the epidemic evolves and another where the information spreads. We theoretically derive the explicit formulas for the basic reproduction number of awareness R 0 a by analyzing the self-consistent equation and the basic reproduction number of disease R 0 d by using the next generation matrix. The formula of R 0 d shows that the effect of awareness can reduce the basic reproduction number of disease. In particular, when awareness does not affect epidemic spreading, R 0 d is shown to match the existing theoretical results. Furthermore, we demonstrate that the disease-free equilibrium is globally asymptotically stable if R 0 d <1; and the endemic equilibrium is globally asymptotically stable if R 0 d >1. Finally, numerical simulations show that information plays a vital role in preventing and controlling disease and effectively reduces the final disease scale. Copyright © 2018 Elsevier Ltd. All rights reserved.

Are Equilibrium Multichannel Networks Predictable? the Case of the Indus River, Pakistan

NASA Astrophysics Data System (ADS)

Darby, S. E.; Carling, P. A.

2017-12-01

Focusing on the specific case of the Indus River, we argue that the equilibrium planform network structure of large, multi-channel, rivers is predictable. Between Chashma and Taunsa, Pakistan, the Indus is a 264 km long multiple-channel reach. Remote sensing imagery, including a period of time that encompasses the occurrence of major floods in 2007 and 2010, shows that Indus has a minimum of two and a maximum of nine channels, with on average four active channels during the dry season and five during the monsoon. We show that the network structure, if not detailed planform, remains stable, even for the record 2010 flood (27,100 m3s-1; recurrence interval > 100 years). Bankline recession is negligible for discharges less than a peak annual discharge of 6,000 m3s-1 ( 80% of mean annual flow). Maximum Flow Efficiency (MFE) principle demonstrates the channel network is insensitive to the monsoon floods, which typically peak at 13,200 m3s-1. Rather, the network is in near-equilibrium with the mean annual flood (7,530 m3s-1). MFE principle indicates stable networks have three to four channels, thus the observed stability in the number of active channels accords with the presence of a near-equilibrium reach-scale channel network. Insensitivity to the annual hydrological cycle demonstrates that the time-scale for network adjustment is much longer than the time-scale of the monsoon hydrograph, with the annual excess water being stored on floodplains, rather than being conveyed in an enlarged channel network. The analysis explains the lack of significant channel adjustment following the largest flood in 40 years and the extensive Indus flooding experienced on an annual basis, with its substantial impacts on the populace and agricultural production.

A novel neural network for variational inequalities with linear and nonlinear constraints.

PubMed

Gao, Xing-Bao; Liao, Li-Zhi; Qi, Liqun

2005-11-01

Variational inequality is a uniform approach for many important optimization and equilibrium problems. Based on the sufficient and necessary conditions of the solution, this paper presents a novel neural network model for solving variational inequalities with linear and nonlinear constraints. Three sufficient conditions are provided to ensure that the proposed network with an asymmetric mapping is stable in the sense of Lyapunov and converges to an exact solution of the original problem. Meanwhile, the proposed network with a gradient mapping is also proved to be stable in the sense of Lyapunov and to have a finite-time convergence under some mild condition by using a new energy function. Compared with the existing neural networks, the new model can be applied to solve some nonmonotone problems, has no adjustable parameter, and has lower complexity. Thus, the structure of the proposed network is very simple. Since the proposed network can be used to solve a broad class of optimization problems, it has great application potential. The validity and transient behavior of the proposed neural network are demonstrated by several numerical examples.

Entrainment in nerve by a ferroelectric model (II): Quasi-periodic oscillation and the phase locking

NASA Astrophysics Data System (ADS)

Shirane, Kotaro; Tokimoto, Takayuki; Kushibe, Hiroyuki

1997-09-01

A nonlinear state equation for membrane excitation can be simplified by Leuchtag's ferroelectric model which is applied to a chemical network theory. A dissipative structure of such a membrane is described by an equilibrium space, η 3 + aη + b = 0, giving a cusp catastrophe, and the membrane is self-organized in the resting state under the condition, a < 0( T < Tc), where η corresponds to the membrane potential, and a and b imply dipole-dipole and dipole-ion interactions of channel proteins embedded in the membrane, respectively. As well known, a specific characteristic of nonlinear electrical phenomena in the membrane is a limit cycle arising through the entrainment by periodical stimuli or chaos. A phase transition between the equilibrium and the non-equilibrium states (a dissipative structure without the resting state) is described by a parameter giving the difference from thermal equilibrium. In this dynamic system, quasi-periodic oscillations which arise in periodic external fields and the phase locking, that is, entrainment, caused by changing I0 at ω ≠ ω n (ω n - the natural frequency of the membrane) are studied with parameters introduced into Zeeman's formulas of ȧ and ḃ.

Quasi-equilibria in reduced Liouville spaces.

PubMed

Halse, Meghan E; Dumez, Jean-Nicolas; Emsley, Lyndon

2012-06-14

The quasi-equilibrium behaviour of isolated nuclear spin systems in full and reduced Liouville spaces is discussed. We focus in particular on the reduced Liouville spaces used in the low-order correlations in Liouville space (LCL) simulation method, a restricted-spin-space approach to efficiently modelling the dynamics of large networks of strongly coupled spins. General numerical methods for the calculation of quasi-equilibrium expectation values of observables in Liouville space are presented. In particular, we treat the cases of a time-independent Hamiltonian, a time-periodic Hamiltonian (with and without stroboscopic sampling) and powder averaging. These quasi-equilibrium calculation methods are applied to the example case of spin diffusion in solid-state nuclear magnetic resonance. We show that there are marked differences between the quasi-equilibrium behaviour of spin systems in the full and reduced spaces. These differences are particularly interesting in the time-periodic-Hamiltonian case, where simulations carried out in the reduced space demonstrate ergodic behaviour even for small spins systems (as few as five homonuclei). The implications of this ergodic property on the success of the LCL method in modelling the dynamics of spin diffusion in magic-angle spinning experiments of powders is discussed.

Building New Bridges between In Vitro and In Vivo in Early Drug Discovery: Where Molecular Modeling Meets Systems Biology.

PubMed

Pearlstein, Robert A; McKay, Daniel J J; Hornak, Viktor; Dickson, Callum; Golosov, Andrei; Harrison, Tyler; Velez-Vega, Camilo; Duca, José

2017-01-01

Cellular drug targets exist within networked function-generating systems whose constituent molecular species undergo dynamic interdependent non-equilibrium state transitions in response to specific perturbations (i.e.. inputs). Cellular phenotypic behaviors are manifested through the integrated behaviors of such networks. However, in vitro data are frequently measured and/or interpreted with empirical equilibrium or steady state models (e.g. Hill, Michaelis-Menten, Briggs-Haldane) relevant to isolated target populations. We propose that cells act as analog computers, "solving" sets of coupled "molecular differential equations" (i.e. represented by populations of interacting species)via "integration" of the dynamic state probability distributions among those populations. Disconnects between biochemical and functional/phenotypic assays (cellular/in vivo) may arise with targetcontaining systems that operate far from equilibrium, and/or when coupled contributions (including target-cognate partner binding and drug pharmacokinetics) are neglected in the analysis of biochemical results. The transformation of drug discovery from a trial-and-error endeavor to one based on reliable design criteria depends on improved understanding of the dynamic mechanisms powering cellular function/dysfunction at the systems level. Here, we address the general mechanisms of molecular and cellular function and pharmacological modulation thereof. We outline a first principles theory on the mechanisms by which free energy is stored and transduced into biological function, and by which biological function is modulated by drug-target binding. We propose that cellular function depends on dynamic counter-balanced molecular systems necessitated by the exponential behavior of molecular state transitions under non-equilibrium conditions, including positive versus negative mass action kinetics and solute-induced perturbations to the hydrogen bonds of solvating water versus kT. Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.org.

Nonlinear relative-proportion-based route adjustment process for day-to-day traffic dynamics: modeling, equilibrium and stability analysis

NASA Astrophysics Data System (ADS)

Zhu, Wenlong; Ma, Shoufeng; Tian, Junfang; Li, Geng

2016-11-01

Travelers' route adjustment behaviors in a congested road traffic network are acknowledged as a dynamic game process between them. Existing Proportional-Switch Adjustment Process (PSAP) models have been extensively investigated to characterize travelers' route choice behaviors; PSAP has concise structure and intuitive behavior rule. Unfortunately most of which have some limitations, i.e., the flow over adjustment problem for the discrete PSAP model, the absolute cost differences route adjustment problem, etc. This paper proposes a relative-Proportion-based Route Adjustment Process (rePRAP) maintains the advantages of PSAP and overcomes these limitations. The rePRAP describes the situation that travelers on higher cost route switch to those with lower cost at the rate that is unilaterally depended on the relative cost differences between higher cost route and its alternatives. It is verified to be consistent with the principle of the rational behavior adjustment process. The equivalence among user equilibrium, stationary path flow pattern and stationary link flow pattern is established, which can be applied to judge whether a given network traffic flow has reached UE or not by detecting the stationary or non-stationary state of link flow pattern. The stability theorem is proved by the Lyapunov function approach. A simple example is tested to demonstrate the effectiveness of the rePRAP model.

Atomistic to continuum modeling of solidification microstructures

DOE PAGES

Karma, Alain; Tourret, Damien

2015-09-26

We summarize recent advances in modeling of solidification microstructures using computational methods that bridge atomistic to continuum scales. We first discuss progress in atomistic modeling of equilibrium and non-equilibrium solid–liquid interface properties influencing microstructure formation, as well as interface coalescence phenomena influencing the late stages of solidification. The latter is relevant in the context of hot tearing reviewed in the article by M. Rappaz in this issue. We then discuss progress to model microstructures on a continuum scale using phase-field methods. We focus on selected examples in which modeling of 3D cellular and dendritic microstructures has been directly linked tomore » experimental observations. Finally, we discuss a recently introduced coarse-grained dendritic needle network approach to simulate the formation of well-developed dendritic microstructures. The approach reliably bridges the well-separated scales traditionally simulated by phase-field and grain structure models, hence opening new avenues for quantitative modeling of complex intra- and inter-grain dynamical interactions on a grain scale.« less

Current-flow efficiency of networks

NASA Astrophysics Data System (ADS)

Liu, Kai; Yan, Xiaoyong

2018-02-01

Many real-world networks, from infrastructure networks to social and communication networks, can be formulated as flow networks. How to realistically measure the transport efficiency of these networks is of fundamental importance. The shortest-path-based efficiency measurement has limitations, as it assumes that flow travels only along those shortest paths. Here, we propose a new metric named current-flow efficiency, in which we calculate the average reciprocal effective resistance between all pairs of nodes in the network. This metric takes the multipath effect into consideration and is more suitable for measuring the efficiency of many real-world flow equilibrium networks. Moreover, this metric can handle a disconnected graph and can thus be used to identify critical nodes and edges from the efficiency-loss perspective. We further analyze how the topological structure affects the current-flow efficiency of networks based on some model and real-world networks. Our results enable a better understanding of flow networks and shed light on the design and improvement of such networks with higher transport efficiency.

Carbon emissions tax policy of urban road traffic and its application in Panjin, China

PubMed Central

Yang, Longhai; Fang, Lin

2018-01-01

How to effectively solve traffic congestion and transportation pollution in urban development is a main research emphasis for transportation management agencies. A carbon emissions tax can affect travelers’ generalized costs and will lead to changes in passenger demand, mode choice and traffic flow equilibrium in road networks, which are of significance in green travel and low-carbon transportation management. This paper first established a mesoscopic model to calculate the carbon emissions tax and determined the value of this charge in China, which was based on road traffic flow, vehicle speed, and carbon emissions. Referring to existing research results to calibrate the value of time, this paper modified the traveler’s generalized cost function, including the carbon emissions tax, fuel surcharge and travel time cost, which can be used in the travel impedance model with the consideration of the carbon emissions tax. Then, a method for analyzing urban road network traffic flow distribution was put forward, and a joint traffic distribution model was established, which considered the relationship between private cars and taxis. Finally, this paper took the city of Panjin as an example to analyze the road traffic carbon emissions tax’s impact. The results illustrated that the carbon emissions tax has a positive effect on road network flow equilibrium and carbon emission reduction. This paper will have good reference value and practical significance for the calculation and implementation of urban traffic carbon emissions taxes in China. PMID:29738580

Carbon emissions tax policy of urban road traffic and its application in Panjin, China.

PubMed

Yang, Longhai; Hu, Xiaowei; Fang, Lin

2018-01-01

How to effectively solve traffic congestion and transportation pollution in urban development is a main research emphasis for transportation management agencies. A carbon emissions tax can affect travelers' generalized costs and will lead to changes in passenger demand, mode choice and traffic flow equilibrium in road networks, which are of significance in green travel and low-carbon transportation management. This paper first established a mesoscopic model to calculate the carbon emissions tax and determined the value of this charge in China, which was based on road traffic flow, vehicle speed, and carbon emissions. Referring to existing research results to calibrate the value of time, this paper modified the traveler's generalized cost function, including the carbon emissions tax, fuel surcharge and travel time cost, which can be used in the travel impedance model with the consideration of the carbon emissions tax. Then, a method for analyzing urban road network traffic flow distribution was put forward, and a joint traffic distribution model was established, which considered the relationship between private cars and taxis. Finally, this paper took the city of Panjin as an example to analyze the road traffic carbon emissions tax's impact. The results illustrated that the carbon emissions tax has a positive effect on road network flow equilibrium and carbon emission reduction. This paper will have good reference value and practical significance for the calculation and implementation of urban traffic carbon emissions taxes in China.

Stochastic dynamics for idiotypic immune networks

NASA Astrophysics Data System (ADS)

Barra, Adriano; Agliari, Elena

2010-12-01

In this work we introduce and analyze the stochastic dynamics obeyed by a model of an immune network recently introduced by the authors. We develop Fokker-Planck equations for the single lymphocyte behavior and coarse grained Langevin schemes for the averaged clone behavior. After showing agreement with real systems (as a short path Jerne cascade), we suggest, both with analytical and numerical arguments, explanations for the generation of (metastable) memory cells, improvement of the secondary response (both in the quality and quantity) and bell shaped modulation against infections as a natural behavior. The whole emerges from the model without being postulated a-priori as it often occurs in second generation immune networks: so the aim of the work is to present some out-of-equilibrium features of this model and to highlight mechanisms which can replace a-priori assumptions in view of further detailed analysis in theoretical systemic immunology.

A Game-Theoretic Response Strategy for Coordinator Attack in Wireless Sensor Networks

PubMed Central

Liu, Jianhua; Yue, Guangxue; Shang, Huiliang; Li, Hongjie

2014-01-01

The coordinator is a specific node that controls the whole network and has a significant impact on the performance in cooperative multihop ZigBee wireless sensor networks (ZWSNs). However, the malicious node attacks coordinator nodes in an effort to waste the resources and disrupt the operation of the network. Attacking leads to a failure of one round of communication between the source nodes and destination nodes. Coordinator selection is a technique that can considerably defend against attack and reduce the data delivery delay, and increase network performance of cooperative communications. In this paper, we propose an adaptive coordinator selection algorithm using game and fuzzy logic aiming at both minimizing the average number of hops and maximizing network lifetime. The proposed game model consists of two interrelated formulations: a stochastic game for dynamic defense and a best response policy using evolutionary game formulation for coordinator selection. The stable equilibrium best policy to response defense is obtained from this game model. It is shown that the proposed scheme can improve reliability and save energy during the network lifetime with respect to security. PMID:25105171

A game-theoretic response strategy for coordinator attack in wireless sensor networks.

PubMed

Liu, Jianhua; Yue, Guangxue; Shen, Shigen; Shang, Huiliang; Li, Hongjie

2014-01-01

The coordinator is a specific node that controls the whole network and has a significant impact on the performance in cooperative multihop ZigBee wireless sensor networks (ZWSNs). However, the malicious node attacks coordinator nodes in an effort to waste the resources and disrupt the operation of the network. Attacking leads to a failure of one round of communication between the source nodes and destination nodes. Coordinator selection is a technique that can considerably defend against attack and reduce the data delivery delay, and increase network performance of cooperative communications. In this paper, we propose an adaptive coordinator selection algorithm using game and fuzzy logic aiming at both minimizing the average number of hops and maximizing network lifetime. The proposed game model consists of two interrelated formulations: a stochastic game for dynamic defense and a best response policy using evolutionary game formulation for coordinator selection. The stable equilibrium best policy to response defense is obtained from this game model. It is shown that the proposed scheme can improve reliability and save energy during the network lifetime with respect to security.

The physics of complex systems in information and biology

NASA Astrophysics Data System (ADS)

Walker, Dylan

Citation networks have re-emerged as a topic intense interest in the complex networks community with the recent availability of large-scale data sets. The ranking of citation networks is a necessary practice as a means to improve information navigability and search. Unlike many information networks, the aging characteristics of citation networks require the development of new ranking methods. To account for strong aging characteristics of citation networks, we modify the PageRank algorithm by initially distributing random surfers exponentially with age, in favor of more recent publications. The output of this algorithm, which we call CiteRank, is interpreted as approximate traffic to individual publications in a simple model of how researchers find new information. We optimize parameters of our algorithm to achieve the best performance. The results are compared for two rather different citation networks: all American Physical Society publications between 1893-2003 and the set of high-energy physics theory (hep-th) preprints. Despite major differences between these two networks, we find that their optimal parameters for the CiteRank algorithm are remarkably similar. The advantages and performance of CiteRank over more conventional methods of ranking publications are discussed. Collaborative voting systems have emerged as an abundant form of real-world, complex information systems that exist in a variety of online applications. These systems are comprised of large populations of users that collectively submit and vote on objects. While the specific properties of these systems vary widely, many of them share a core set of features and dynamical behaviors that govern their evolution. We study a subset of these systems that involve material of a time-critical nature as in the popular example of news items. We consider a general model system in which articles are introduced, voted on by a population of users, and subsequently expire after a proscribed period of time. To study the interaction between popularity and quality, we introduce simple stochastic models of user behavior that approximate differing user quality and susceptibility to the common notion of popularity. We define a metric to quantify user reputation in a manner that is self-consistent, adaptable and content-blind and shows good correlation with the probability that a user behaves in an optimal fashion. We further construct a mechanism for ranking documents that take into account user reputation and provides substantial improvement in the time-critical performance of the system. The structure of complex systems have been well studied in the context of both information and biological systems. More recently, dynamics in complex systems that occur over the background of the underlying network has received a great deal of attention. In particular, the study of fluctuations in complex systems has emerged as an issue central to understanding dynamical behavior. We approach the problem of collective effects of the underlying network on dynamical fluctuations by considering the protein-protein interaction networks for the system of the living cell. We consider two types of fluctuations in the mass-action equilibrium in protein binding networks. The first type is driven by relatively slow changes in total concentrations (copy numbers) of interacting proteins. The second type, to which we refer to as spontaneous, is caused by quickly decaying thermodynamic deviations away from the mass-action equilibrium of the system. As such they are amenable to methods of equilibrium statistical mechanics used in our study. We investigate the effects of network connectivity on these fluctuations by comparing them to different scenarios in which the interacting pair is isolated form the rest of the network. Such comparison allows us to analytically derive upper and lower bounds on network fluctuations. The collective effects are shown to sometimes lead to relatively large amplification of spontaneous fluctuations as compared to the expectation for isolated dimers. As a consequence of this, the strength of both types of fluctuations is positively correlated with the overall network connectivity of proteins forming the complex. On the other hand, the relative amplitude of fluctuations is negatively correlated with the equilibrium concentration of the complex. Our general findings are illustrated using a curated network of protein-protein interactions and multi-protein complexes in bakers yeast with experimentally determined protein concentrations.

Thermodynamically Feasible Kinetic Models of Reaction Networks

PubMed Central

Ederer, Michael; Gilles, Ernst Dieter

2007-01-01

The dynamics of biological reaction networks are strongly constrained by thermodynamics. An holistic understanding of their behavior and regulation requires mathematical models that observe these constraints. However, kinetic models may easily violate the constraints imposed by the principle of detailed balance, if no special care is taken. Detailed balance demands that in thermodynamic equilibrium all fluxes vanish. We introduce a thermodynamic-kinetic modeling (TKM) formalism that adapts the concepts of potentials and forces from irreversible thermodynamics to kinetic modeling. In the proposed formalism, the thermokinetic potential of a compound is proportional to its concentration. The proportionality factor is a compound-specific parameter called capacity. The thermokinetic force of a reaction is a function of the potentials. Every reaction has a resistance that is the ratio of thermokinetic force and reaction rate. For mass-action type kinetics, the resistances are constant. Since it relies on the thermodynamic concept of potentials and forces, the TKM formalism structurally observes detailed balance for all values of capacities and resistances. Thus, it provides an easy way to formulate physically feasible, kinetic models of biological reaction networks. The TKM formalism is useful for modeling large biological networks that are subject to many detailed balance relations. PMID:17208985

Emergent inequality and self-organized social classes in a network of power and frustration

DOE PAGES

Mahault, Benoit; Saxena, Avadh; Nisoli, Cristiano

2017-02-17

We propose a simple agent-based model on a network to conceptualize the allocation of limited wealth among more abundant expectations at the interplay of power, frustration, and initiative. Concepts imported from the statistical physics of frustrated systems in and out of equilibrium allow us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from or lose wealth to anybody else invariably leads to a complete polarization of the distribution ofmore » wealth vs. opportunity. This picture is however dramatically ameliorated when hard constraints are imposed over agents in the form of a limiting network of transactions. There, an out of equilibrium dynamics of the networks, based on a competition between power and frustration in the decision-making of agents, leads to network coevolution. The ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of equality. It also leads, for proper values of social initiative, to the emergence of three self-organized social classes, lower, middle, and upper class. Their dynamics, which appears mostly controlled by the middle class, drives a cyclical regime of dramatic social changes.« less

Emergent inequality and self-organized social classes in a network of power and frustration

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mahault, Benoit; Saxena, Avadh; Nisoli, Cristiano

We propose a simple agent-based model on a network to conceptualize the allocation of limited wealth among more abundant expectations at the interplay of power, frustration, and initiative. Concepts imported from the statistical physics of frustrated systems in and out of equilibrium allow us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from or lose wealth to anybody else invariably leads to a complete polarization of the distribution ofmore » wealth vs. opportunity. This picture is however dramatically ameliorated when hard constraints are imposed over agents in the form of a limiting network of transactions. There, an out of equilibrium dynamics of the networks, based on a competition between power and frustration in the decision-making of agents, leads to network coevolution. The ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of equality. It also leads, for proper values of social initiative, to the emergence of three self-organized social classes, lower, middle, and upper class. Their dynamics, which appears mostly controlled by the middle class, drives a cyclical regime of dramatic social changes.« less

Emergent inequality and self-organized social classes in a network of power and frustration

PubMed Central

Mahault, Benoit; Saxena, Avadh

2017-01-01

We propose a simple agent-based model on a network to conceptualize the allocation of limited wealth among more abundant expectations at the interplay of power, frustration, and initiative. Concepts imported from the statistical physics of frustrated systems in and out of equilibrium allow us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from or lose wealth to anybody else invariably leads to a complete polarization of the distribution of wealth vs. opportunity. This picture is however dramatically ameliorated when hard constraints are imposed over agents in the form of a limiting network of transactions. There, an out of equilibrium dynamics of the networks, based on a competition between power and frustration in the decision-making of agents, leads to network coevolution. The ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of equality. It also leads, for proper values of social initiative, to the emergence of three self-organized social classes, lower, middle, and upper class. Their dynamics, which appears mostly controlled by the middle class, drives a cyclical regime of dramatic social changes. PMID:28212440

Emergent inequality and self-organized social classes in a network of power and frustration.

PubMed

Mahault, Benoit; Saxena, Avadh; Nisoli, Cristiano

2017-01-01

We propose a simple agent-based model on a network to conceptualize the allocation of limited wealth among more abundant expectations at the interplay of power, frustration, and initiative. Concepts imported from the statistical physics of frustrated systems in and out of equilibrium allow us to compare subjective measures of frustration and satisfaction to collective measures of fairness in wealth distribution, such as the Lorenz curve and the Gini index. We find that a completely libertarian, law-of-the-jungle setting, where every agent can acquire wealth from or lose wealth to anybody else invariably leads to a complete polarization of the distribution of wealth vs. opportunity. This picture is however dramatically ameliorated when hard constraints are imposed over agents in the form of a limiting network of transactions. There, an out of equilibrium dynamics of the networks, based on a competition between power and frustration in the decision-making of agents, leads to network coevolution. The ratio of power and frustration controls different dynamical regimes separated by kinetic transitions and characterized by drastically different values of equality. It also leads, for proper values of social initiative, to the emergence of three self-organized social classes, lower, middle, and upper class. Their dynamics, which appears mostly controlled by the middle class, drives a cyclical regime of dramatic social changes.

Mechanical response of biopolymer double networks

NASA Astrophysics Data System (ADS)

Carroll, Joshua; Das, Moumita

We investigate a double network model of articular cartilage (AC) and characterize its equilibrium mechanical response. AC has very few cells and the extracellular matrix mainly determines its mechanical response. This matrix can be thought of as a double polymer network made of collagen and aggrecan. The collagen fibers are stiff and resist tension and compression forces, while aggrecans are flexible and control swelling and hydration. We construct a microscopic model made of two interconnected disordered polymer networks, with fiber elasticity chosen to qualitatively mimic the experimental system. We study the collective mechanical response of this double network as a function of the concentration and stiffness of the individual components as well as the strength of the connection between them using rigidity percolation theory. Our results may provide a better understanding of mechanisms underlying the mechanical resilience of AC, and more broadly may also lead to new perspectives on the mechanical response of multicomponent soft materials. This work was partially supported by a Cottrell College Science Award.

Global exponential stability of neutral high-order stochastic Hopfield neural networks with Markovian jump parameters and mixed time delays.

PubMed

Huang, Haiying; Du, Qiaosheng; Kang, Xibing

2013-11-01

In this paper, a class of neutral high-order stochastic Hopfield neural networks with Markovian jump parameters and mixed time delays is investigated. The jumping parameters are modeled as a continuous-time finite-state Markov chain. At first, the existence of equilibrium point for the addressed neural networks is studied. By utilizing the Lyapunov stability theory, stochastic analysis theory and linear matrix inequality (LMI) technique, new delay-dependent stability criteria are presented in terms of linear matrix inequalities to guarantee the neural networks to be globally exponentially stable in the mean square. Numerical simulations are carried out to illustrate the main results. © 2013 ISA. Published by ISA. All rights reserved.

Multistability of neural networks with discontinuous non-monotonic piecewise linear activation functions and time-varying delays.

PubMed

Nie, Xiaobing; Zheng, Wei Xing

2015-05-01

This paper is concerned with the problem of coexistence and dynamical behaviors of multiple equilibrium points for neural networks with discontinuous non-monotonic piecewise linear activation functions and time-varying delays. The fixed point theorem and other analytical tools are used to develop certain sufficient conditions that ensure that the n-dimensional discontinuous neural networks with time-varying delays can have at least 5(n) equilibrium points, 3(n) of which are locally stable and the others are unstable. The importance of the derived results is that it reveals that the discontinuous neural networks can have greater storage capacity than the continuous ones. Moreover, different from the existing results on multistability of neural networks with discontinuous activation functions, the 3(n) locally stable equilibrium points obtained in this paper are located in not only saturated regions, but also unsaturated regions, due to the non-monotonic structure of discontinuous activation functions. A numerical simulation study is conducted to illustrate and support the derived theoretical results. Copyright © 2015 Elsevier Ltd. All rights reserved.

Ensemble ecosystem modeling for predicting ecosystem response to predator reintroduction.

PubMed

Baker, Christopher M; Gordon, Ascelin; Bode, Michael

2017-04-01

Introducing a new or extirpated species to an ecosystem is risky, and managers need quantitative methods that can predict the consequences for the recipient ecosystem. Proponents of keystone predator reintroductions commonly argue that the presence of the predator will restore ecosystem function, but this has not always been the case, and mathematical modeling has an important role to play in predicting how reintroductions will likely play out. We devised an ensemble modeling method that integrates species interaction networks and dynamic community simulations and used it to describe the range of plausible consequences of 2 keystone-predator reintroductions: wolves (Canis lupus) to Yellowstone National Park and dingoes (Canis dingo) to a national park in Australia. Although previous methods for predicting ecosystem responses to such interventions focused on predicting changes around a given equilibrium, we used Lotka-Volterra equations to predict changing abundances through time. We applied our method to interaction networks for wolves in Yellowstone National Park and for dingoes in Australia. Our model replicated the observed dynamics in Yellowstone National Park and produced a larger range of potential outcomes for the dingo network. However, we also found that changes in small vertebrates or invertebrates gave a good indication about the potential future state of the system. Our method allowed us to predict when the systems were far from equilibrium. Our results showed that the method can also be used to predict which species may increase or decrease following a reintroduction and can identify species that are important to monitor (i.e., species whose changes in abundance give extra insight into broad changes in the system). Ensemble ecosystem modeling can also be applied to assess the ecosystem-wide implications of other types of interventions including assisted migration, biocontrol, and invasive species eradication. © 2016 Society for Conservation Biology.

The Local Structure of Globalization. The Network Dynamics of Foreign Direct Investments in the International Electricity Industry

NASA Astrophysics Data System (ADS)

Koskinen, Johan; Lomi, Alessandro

2013-05-01

We study the evolution of the network of foreign direct investment (FDI) in the international electricity industry during the period 1994-2003. We assume that the ties in the network of investment relations between countries are created and deleted in continuous time, according to a conditional Gibbs distribution. This assumption allows us to take simultaneously into account the aggregate predictions of the well-established gravity model of international trade as well as local dependencies between network ties connecting the countries in our sample. According to the modified version of the gravity model that we specify, the probability of observing an investment tie between two countries depends on the mass of the economies involved, their physical distance, and the tendency of the network to self-organize into local configurations of network ties. While the limiting distribution of the data generating process is an exponential random graph model, we do not assume the system to be in equilibrium. We find evidence of the effects of the standard gravity model of international trade on evolution of the global FDI network. However, we also provide evidence of significant dyadic and extra-dyadic dependencies between investment ties that are typically ignored in available research. We show that local dependencies between national electricity industries are sufficient for explaining global properties of the network of foreign direct investments. We also show, however, that network dependencies vary significantly over time giving rise to a time-heterogeneous localized process of network evolution.

An equal force theory for network models of soft materials with arbitrary molecular weight distribution

NASA Astrophysics Data System (ADS)

Verron, E.; Gros, A.

2017-09-01

Most network models for soft materials, e.g. elastomers and gels, are dedicated to idealized materials: all chains admit the same number of Kuhn segments. Nevertheless, such standard models are not appropriate for materials involving multiple networks, and some specific constitutive equations devoted to these materials have been derived in the last few years. In nearly all cases, idealized networks of different chain lengths are assembled following an equal strain assumption; only few papers adopt an equal stress assumption, although some authors argue that such hypothesis would reflect the equilibrium of the different networks in contact. In this work, a full-network model with an arbitrary chain length distribution is derived by considering that chains of different lengths satisfy the equal force assumption in each direction of the unit sphere. The derivation is restricted to non-Gaussian freely jointed chains and to affine deformation of the sphere. Firstly, after a proper definition of the undeformed configuration of the network, we demonstrate that the equal force assumption leads to the equality of a normalized stretch in chains of different lengths. Secondly, we establish that the network with chain length distribution behaves as an idealized full-network of which both chain length and density of are provided by the chain length distribution. This approach is finally illustrated with two examples: the derivation of a new expression for the Young modulus of bimodal interpenetrated polymer networks, and the prediction of the change in fluorescence during deformation of mechanochemically responsive elastomers.

Economic Feasibility of Wireless Sensor Network-Based Service Provision in a Duopoly Setting with a Monopolist Operator

PubMed Central

Romero, Julián; Sacoto-Cabrera, Erwin J.

2017-01-01

We analyze the feasibility of providing Wireless Sensor Network-data-based services in an Internet of Things scenario from an economical point of view. The scenario has two competing service providers with their own private sensor networks, a network operator and final users. The scenario is analyzed as two games using game theory. In the first game, sensors decide to subscribe or not to the network operator to upload the collected sensing-data, based on a utility function related to the mean service time and the price charged by the operator. In the second game, users decide to subscribe or not to the sensor-data-based service of the service providers based on a Logit discrete choice model related to the quality of the data collected and the subscription price. The sinks and users subscription stages are analyzed using population games and discrete choice models, while network operator and service providers pricing stages are analyzed using optimization and Nash equilibrium concepts respectively. The model is shown feasible from an economic point of view for all the actors if there are enough interested final users and opens the possibility of developing more efficient models with different types of services. PMID:29186847

Broken Detailed Balance of Filament Dynamics in Active Networks

NASA Astrophysics Data System (ADS)

Schmidt, Christoph F.; Gladrow, Jannes; Fakhri, Nikta; Mackintosh, Fred C.; Broedersz, Chase

Endogenous embedded semiflexible filaments such as microtubules, or added filaments such as single- walled carbon nanotubes can be used as novel tools to noninvasively track equilibrium and nonequilibrium fluctuations in biopolymer networks. We analytically calculated shape fluctuations of semi- flexible probe filaments in a viscoelastic environment, driven out of equilibrium by motor activity. Transverse bending fluctuations of the probe filaments can be decomposed into dynamic normal modes. We find that these modes no longer evolve independently under non-equilibrium driving. This effective mode coupling results in nonzero circulatory currents in a conformational phase space, reflecting a violation of detailed balance. We present predictions for the characteristic frequencies associated with these currents and investigate how the temporal signatures of motor activity determine mode correlations, which we find to be consistent with recent experiments on microtubules embedded in cytoskeletal networks.

Network Performance Evaluation Model for assessing the impacts of high-occupancy vehicle facilities

DOE Office of Scientific and Technical Information (OSTI.GOV)

Janson, B.N.; Zozaya-Gorostiza, C.; Southworth, F.

1986-09-01

A model to assess the impacts of major high-occupancy vehicle (HOV) facilities on regional levels of energy consumption and vehicle air pollution emissions in urban aeas is developed and applied. This model can be used to forecast and compare the impacts of alternative HOV facility design and operation plans on traffic patterns, travel costs, model choice, travel demand, energy consumption and vehicle emissions. The model is designed to show differences in the overall impacts of alternative HOV facility types, locations and operation plans rather than to serve as a tool for detailed engineering design and traffic planning studies. The Networkmore » Performance Evaluation Model (NETPEM) combines several urban transportation planning models within a multi-modal network equilibrium framework including modules with which to define the type, location and use policy of the HOV facility to be tested, and to assess the impacts of this facility.« less

Mass-action equilibrium and non-specific interactions in protein binding networks

NASA Astrophysics Data System (ADS)

Maslov, Sergei

2009-03-01

Large-scale protein binding networks serve as a paradigm of complex properties of living cells. These networks are naturally weighted with edges characterized by binding strength and protein-nodes -- by their concentrations. However, the state-of-the-art high-throughput experimental techniques generate just a binary (yes or no) information about individual interactions. As a result, most of the previous research concentrated just on topology of these networks. In a series of recent publications [1-4] my collaborators and I went beyond purely topological studies and calculated the mass-action equilibrium of a genome-wide binding network using experimentally determined protein concentrations, localizations, and reliable binding interactions in baker's yeast. We then studied how this equilibrium responds to large perturbations [1-2] and noise [3] in concentrations of proteins. We demonstrated that the change in the equilibrium concentration of a protein exponentially decays (and sign-alternates) with its network distance away from the perturbed node. This explains why, despite a globally connected topology, individual functional modules in such networks are able to operate fairly independently. In a separate study [4] we quantified the interplay between specific and non-specific binding interactions under crowded conditions inside living cells. We show how the need to limit the waste of resources constrains the number of types and concentrations of proteins that are present at the same time and at the same place in yeast cells. [1] S Maslov, I. Ispolatov, PNAS 104:13655 (2007). [2] S. Maslov, K. Sneppen, I. Ispolatov, New J. of Phys. 9: 273 (2007). [3] K-K. Yan, D. Walker, S. Maslov, PRL accepted (2008). [4] J. Zhang, S. Maslov, and E. I. Shakhnovich, Mol Syst Biol 4, 210 (2008).

Modeling Psychological Attributes in Psychology – An Epistemological Discussion: Network Analysis vs. Latent Variables

PubMed Central

Guyon, Hervé; Falissard, Bruno; Kop, Jean-Luc

2017-01-01

Network Analysis is considered as a new method that challenges Latent Variable models in inferring psychological attributes. With Network Analysis, psychological attributes are derived from a complex system of components without the need to call on any latent variables. But the ontological status of psychological attributes is not adequately defined with Network Analysis, because a psychological attribute is both a complex system and a property emerging from this complex system. The aim of this article is to reappraise the legitimacy of latent variable models by engaging in an ontological and epistemological discussion on psychological attributes. Psychological attributes relate to the mental equilibrium of individuals embedded in their social interactions, as robust attractors within complex dynamic processes with emergent properties, distinct from physical entities located in precise areas of the brain. Latent variables thus possess legitimacy, because the emergent properties can be conceptualized and analyzed on the sole basis of their manifestations, without exploring the upstream complex system. However, in opposition with the usual Latent Variable models, this article is in favor of the integration of a dynamic system of manifestations. Latent Variables models and Network Analysis thus appear as complementary approaches. New approaches combining Latent Network Models and Network Residuals are certainly a promising new way to infer psychological attributes, placing psychological attributes in an inter-subjective dynamic approach. Pragmatism-realism appears as the epistemological framework required if we are to use latent variables as representations of psychological attributes. PMID:28572780

Non-equilibrium phase transition in mesoscopic biochemical systems: from stochastic to nonlinear dynamics and beyond

PubMed Central

Ge, Hao; Qian, Hong

2011-01-01

A theory for an non-equilibrium phase transition in a driven biochemical network is presented. The theory is based on the chemical master equation (CME) formulation of mesoscopic biochemical reactions and the mathematical method of large deviations. The large deviations theory provides an analytical tool connecting the macroscopic multi-stability of an open chemical system with the multi-scale dynamics of its mesoscopic counterpart. It shows a corresponding non-equilibrium phase transition among multiple stochastic attractors. As an example, in the canonical phosphorylation–dephosphorylation system with feedback that exhibits bistability, we show that the non-equilibrium steady-state (NESS) phase transition has all the characteristics of classic equilibrium phase transition: Maxwell construction, a discontinuous first-derivative of the ‘free energy function’, Lee–Yang's zero for a generating function and a critical point that matches the cusp in nonlinear bifurcation theory. To the biochemical system, the mathematical analysis suggests three distinct timescales and needed levels of description. They are (i) molecular signalling, (ii) biochemical network nonlinear dynamics, and (iii) cellular evolution. For finite mesoscopic systems such as a cell, motions associated with (i) and (iii) are stochastic while that with (ii) is deterministic. Both (ii) and (iii) are emergent properties of a dynamic biochemical network. PMID:20466813

Upscaling of nanoparticle transport in porous media under unfavorable conditions: Pore scale to Darcy scale

NASA Astrophysics Data System (ADS)

Seetha, N.; Raoof, Amir; Mohan Kumar, M. S.; Majid Hassanizadeh, S.

2017-05-01

Transport and deposition of nanoparticles in porous media is a multi-scale problem governed by several pore-scale processes, and hence, it is critical to link the processes at pore scale to the Darcy-scale behavior. In this study, using pore network modeling, we develop correlation equations for deposition rate coefficients for nanoparticle transport under unfavorable conditions at the Darcy scale based on pore-scale mechanisms. The upscaling tool is a multi-directional pore-network model consisting of an interconnected network of pores with variable connectivities. Correlation equations describing the pore-averaged deposition rate coefficients under unfavorable conditions in a cylindrical pore, developed in our earlier studies, are employed for each pore element. Pore-network simulations are performed for a wide range of parameter values to obtain the breakthrough curves of nanoparticle concentration. The latter is fitted with macroscopic 1-D advection-dispersion equation with a two-site linear reversible deposition accounting for both equilibrium and kinetic sorption. This leads to the estimation of three Darcy-scale deposition coefficients: distribution coefficient, kinetic rate constant, and the fraction of equilibrium sites. The correlation equations for the Darcy-scale deposition coefficients, under unfavorable conditions, are provided as a function of measurable Darcy-scale parameters, including: porosity, mean pore throat radius, mean pore water velocity, nanoparticle radius, ionic strength, dielectric constant, viscosity, temperature, and surface potentials of the particle and grain surfaces. The correlation equations are found to be consistent with the available experimental results, and in qualitative agreement with Colloid Filtration Theory for all parameters, except for the mean pore water velocity and nanoparticle radius.

Systems biology and the origins of life? part II. Are biochemical networks possible ancestors of living systems? networks of catalysed chemical reactions: non-equilibrium, self-organization and evolution.

PubMed

Ricard, Jacques

2010-01-01

The present article discusses the possibility that catalysed chemical networks can evolve. Even simple enzyme-catalysed chemical reactions can display this property. The example studied is that of a two-substrate proteinoid, or enzyme, reaction displaying random binding of its substrates A and B. The fundamental property of such a system is to display either emergence or integration depending on the respective values of the probabilities that the enzyme has bound one of its substrate regardless it has bound the other substrate, or, specifically, after it has bound the other substrate. There is emergence of information if p(A)>p(AB) and p(B)>p(BA). Conversely, if p(A)

Modeling of uncertainties in biochemical reactions.

PubMed

Mišković, Ljubiša; Hatzimanikatis, Vassily

2011-02-01

Mathematical modeling is an indispensable tool for research and development in biotechnology and bioengineering. The formulation of kinetic models of biochemical networks depends on knowledge of the kinetic properties of the enzymes of the individual reactions. However, kinetic data acquired from experimental observations bring along uncertainties due to various experimental conditions and measurement methods. In this contribution, we propose a novel way to model the uncertainty in the enzyme kinetics and to predict quantitatively the responses of metabolic reactions to the changes in enzyme activities under uncertainty. The proposed methodology accounts explicitly for mechanistic properties of enzymes and physico-chemical and thermodynamic constraints, and is based on formalism from systems theory and metabolic control analysis. We achieve this by observing that kinetic responses of metabolic reactions depend: (i) on the distribution of the enzymes among their free form and all reactive states; (ii) on the equilibrium displacements of the overall reaction and that of the individual enzymatic steps; and (iii) on the net fluxes through the enzyme. Relying on this observation, we develop a novel, efficient Monte Carlo sampling procedure to generate all states within a metabolic reaction that satisfy imposed constrains. Thus, we derive the statistics of the expected responses of the metabolic reactions to changes in enzyme levels and activities, in the levels of metabolites, and in the values of the kinetic parameters. We present aspects of the proposed framework through an example of the fundamental three-step reversible enzymatic reaction mechanism. We demonstrate that the equilibrium displacements of the individual enzymatic steps have an important influence on kinetic responses of the enzyme. Furthermore, we derive the conditions that must be satisfied by a reversible three-step enzymatic reaction operating far away from the equilibrium in order to respond to changes in metabolite levels according to the irreversible Michelis-Menten kinetics. The efficient sampling procedure allows easy, scalable, implementation of this methodology to modeling of large-scale biochemical networks. © 2010 Wiley Periodicals, Inc.

A stochastic population model to evaluate Moapa dace (Moapa coriacea) population growth under alternative management scenarios

USGS Publications Warehouse

Perry, Russell W.; Jones, Edward; Scoppettone, G. Gary

2015-07-14

Increasing or decreasing the total carrying capacity of all stream segments resulted in changes in equilibrium population size that were directly proportional to the change in capacity. However, changes in carrying capacity to some stream segments but not others could result in disproportionate changes in equilibrium population sizes by altering density-dependent movement and survival in the stream network. These simulations show how our IBM can provide a useful management tool for understanding the effect of restoration actions or reintroductions on carrying capacity, and, in turn, how these changes affect Moapa dace abundance. Such tools are critical for devising management strategies to achieve recovery goals.

Self-Healing of Unentangled Polymer Networks with Reversible Bonds

PubMed Central

Stukalin, Evgeny B.; Cai, Li-Heng; Kumar, N. Arun; Leibler, Ludwik; Rubinstein, Michael

2013-01-01

Self-healing polymeric materials are systems that after damage can revert to their original state with full or partial recovery of mechanical strength. Using scaling theory we study a simple model of autonomic self-healing of unentangled polymer networks. In this model one of the two end monomers of each polymer chain is fixed in space mimicking dangling chains attachment to a polymer network, while the sticky monomer at the other end of each chain can form pairwise reversible bond with the sticky end of another chain. We study the reaction kinetics of reversible bonds in this simple model and analyze the different stages in the self-repair process. The formation of bridges and the recovery of the material strength across the fractured interface during the healing period occur appreciably faster after shorter waiting time, during which the fractured surfaces are kept apart. We observe the slowest formation of bridges for self-adhesion after bringing into contact two bare surfaces with equilibrium (very low) density of open stickers in comparison with self-healing. The primary role of anomalous diffusion in material self-repair for short waiting times is established, while at long waiting times the recovery of bonds across fractured interface is due to hopping diffusion of stickers between different bonded partners. Acceleration in bridge formation for self-healing compared to self-adhesion is due to excess non-equilibrium concentration of open stickers. Full recovery of reversible bonds across fractured interface (formation of bridges) occurs after appreciably longer time than the equilibration time of the concentration of reversible bonds in the bulk. PMID:24347684

A coupled geomorphic and ecological model of tidal marsh evolution.

PubMed

Kirwan, Matthew L; Murray, A Brad

2007-04-10

The evolution of tidal marsh platforms and interwoven channel networks cannot be addressed without treating the two-way interactions that link biological and physical processes. We have developed a 3D model of tidal marsh accretion and channel network development that couples physical sediment transport processes with vegetation biomass productivity. Tidal flow tends to cause erosion, whereas vegetation biomass, a function of bed surface depth below high tide, influences the rate of sediment deposition and slope-driven transport processes such as creek bank slumping. With a steady, moderate rise in sea level, the model builds a marsh platform and channel network with accretion rates everywhere equal to the rate of sea-level rise, meaning water depths and biological productivity remain temporally constant. An increase in the rate of sea-level rise, or a reduction in sediment supply, causes marsh-surface depths, biomass productivity, and deposition rates to increase while simultaneously causing the channel network to expand. Vegetation on the marsh platform can promote a metastable equilibrium where the platform maintains elevation relative to a rapidly rising sea level, although disturbance to vegetation could cause irreversible loss of marsh habitat.

Information-geometric measures as robust estimators of connection strengths and external inputs.

PubMed

Tatsuno, Masami; Fellous, Jean-Marc; Amari, Shun-Ichi

2009-08-01

Information geometry has been suggested to provide a powerful tool for analyzing multineuronal spike trains. Among several advantages of this approach, a significant property is the close link between information-geometric measures and neural network architectures. Previous modeling studies established that the first- and second-order information-geometric measures corresponded to the number of external inputs and the connection strengths of the network, respectively. This relationship was, however, limited to a symmetrically connected network, and the number of neurons used in the parameter estimation of the log-linear model needed to be known. Recently, simulation studies of biophysical model neurons have suggested that information geometry can estimate the relative change of connection strengths and external inputs even with asymmetric connections. Inspired by these studies, we analytically investigated the link between the information-geometric measures and the neural network structure with asymmetrically connected networks of N neurons. We focused on the information-geometric measures of orders one and two, which can be derived from the two-neuron log-linear model, because unlike higher-order measures, they can be easily estimated experimentally. Considering the equilibrium state of a network of binary model neurons that obey stochastic dynamics, we analytically showed that the corrected first- and second-order information-geometric measures provided robust and consistent approximation of the external inputs and connection strengths, respectively. These results suggest that information-geometric measures provide useful insights into the neural network architecture and that they will contribute to the study of system-level neuroscience.

Epidemic processes in complex networks

NASA Astrophysics Data System (ADS)

Pastor-Satorras, Romualdo; Castellano, Claudio; Van Mieghem, Piet; Vespignani, Alessandro

2015-07-01

In recent years the research community has accumulated overwhelming evidence for the emergence of complex and heterogeneous connectivity patterns in a wide range of biological and sociotechnical systems. The complex properties of real-world networks have a profound impact on the behavior of equilibrium and nonequilibrium phenomena occurring in various systems, and the study of epidemic spreading is central to our understanding of the unfolding of dynamical processes in complex networks. The theoretical analysis of epidemic spreading in heterogeneous networks requires the development of novel analytical frameworks, and it has produced results of conceptual and practical relevance. A coherent and comprehensive review of the vast research activity concerning epidemic processes is presented, detailing the successful theoretical approaches as well as making their limits and assumptions clear. Physicists, mathematicians, epidemiologists, computer, and social scientists share a common interest in studying epidemic spreading and rely on similar models for the description of the diffusion of pathogens, knowledge, and innovation. For this reason, while focusing on the main results and the paradigmatic models in infectious disease modeling, the major results concerning generalized social contagion processes are also presented. Finally, the research activity at the forefront in the study of epidemic spreading in coevolving, coupled, and time-varying networks is reported.

Robust stability of interval bidirectional associative memory neural network with time delays.

PubMed

Liao, Xiaofeng; Wong, Kwok-wo

2004-04-01

In this paper, the conventional bidirectional associative memory (BAM) neural network with signal transmission delay is intervalized in order to study the bounded effect of deviations in network parameters and external perturbations. The resultant model is referred to as a novel interval dynamic BAM (IDBAM) model. By combining a number of different Lyapunov functionals with the Razumikhin technique, some sufficient conditions for the existence of unique equilibrium and robust stability are derived. These results are fairly general and can be verified easily. To go further, we extend our investigation to the time-varying delay case. Some robust stability criteria for BAM with perturbations of time-varying delays are derived. Besides, our approach for the analysis allows us to consider several different types of activation functions, including piecewise linear sigmoids with bounded activations as well as the usual C1-smooth sigmoids. We believe that the results obtained have leading significance in the design and application of BAM neural networks.

Simultaneously selecting appropriate partners for gaming and strategy adaptation to enhance network reciprocity in the prisoner's dilemma

NASA Astrophysics Data System (ADS)

Tanimoto, Jun

2014-01-01

Network reciprocity is one mechanism for adding social viscosity, which leads to cooperative equilibrium in 2 × 2 prisoner's dilemma games. Previous studies have shown that cooperation can be enhanced by using a skewed, rather than a random, selection of partners for either strategy adaptation or the gaming process. Here we show that combining both processes for selecting a gaming partner and an adaptation partner further enhances cooperation, provided that an appropriate selection rule and parameters are adopted. We also show that this combined model significantly enhances cooperation by reducing the degree of activity in the underlying network; we measure the degree of activity with a quantity called effective degree. More precisely, during the initial evolutionary stage in which the global cooperation fraction declines because initially allocated cooperators becoming defectors, the model shows that weak cooperative clusters perish and only a few strong cooperative clusters survive. This finding is the most important key to attaining significant network reciprocity.

Reflective equilibrium and empirical data: third person moral experiences in empirical medical ethics.

PubMed

De Vries, Martine; Van Leeuwen, Evert

2010-11-01

In ethics, the use of empirical data has become more and more popular, leading to a distinct form of applied ethics, namely empirical ethics. This 'empirical turn' is especially visible in bioethics. There are various ways of combining empirical research and ethical reflection. In this paper we discuss the use of empirical data in a special form of Reflective Equilibrium (RE), namely the Network Model with Third Person Moral Experiences. In this model, the empirical data consist of the moral experiences of people in a practice. Although inclusion of these moral experiences in this specific model of RE can be well defended, their use in the application of the model still raises important questions. What precisely are moral experiences? How to determine relevance of experiences, in other words: should there be a selection of the moral experiences that are eventually used in the RE? How much weight should the empirical data have in the RE? And the key question: can the use of RE by empirical ethicists really produce answers to practical moral questions? In this paper we start to answer the above questions by giving examples taken from our research project on understanding the norm of informed consent in the field of pediatric oncology. We especially emphasize that incorporation of empirical data in a network model can reduce the risk of self-justification and bias and can increase the credibility of the RE reached. © 2009 Blackwell Publishing Ltd.

Non-equilibrium repressor binding kinetics link DNA damage dose to transcriptional timing within the SOS gene network.

PubMed

Culyba, Matthew J; Kubiak, Jeffrey M; Mo, Charlie Y; Goulian, Mark; Kohli, Rahul M

2018-06-01

Biochemical pathways are often genetically encoded as simple transcription regulation networks, where one transcription factor regulates the expression of multiple genes in a pathway. The relative timing of each promoter's activation and shut-off within the network can impact physiology. In the DNA damage repair pathway (known as the SOS response) of Escherichia coli, approximately 40 genes are regulated by the LexA repressor. After a DNA damaging event, LexA degradation triggers SOS gene transcription, which is temporally separated into subsets of 'early', 'middle', and 'late' genes. Although this feature plays an important role in regulating the SOS response, both the range of this separation and its underlying mechanism are not experimentally defined. Here we show that, at low doses of DNA damage, the timing of promoter activities is not separated. Instead, timing differences only emerge at higher levels of DNA damage and increase as a function of DNA damage dose. To understand mechanism, we derived a series of synthetic SOS gene promoters which vary in LexA-operator binding kinetics, but are otherwise identical, and then studied their activity over a large dose-range of DNA damage. In distinction to established models based on rapid equilibrium assumptions, the data best fit a kinetic model of repressor occupancy at promoters, where the drop in cellular LexA levels associated with higher doses of DNA damage leads to non-equilibrium binding kinetics of LexA at operators. Operators with slow LexA binding kinetics achieve their minimal occupancy state at later times than operators with fast binding kinetics, resulting in a time separation of peak promoter activity between genes. These data provide insight into this remarkable feature of the SOS pathway by demonstrating how a single transcription factor can be employed to control the relative timing of each gene's transcription as a function of stimulus dose.

Hydrolysis of Methylal Catalyzed by Ion Exchange Resins in Aqueous Media

NASA Astrophysics Data System (ADS)

He, Gaoyin; Dai, Fangfang; Shi, Midong; Li, Qingsong; Yu, Yingmin

2018-05-01

In the present work, the chemical equilibrium and kinetics of methylal (PODE1) hydrolysis catalyzed by ion-exchange resin in aqueous solutions were investigated. The study covers temperatures between 333.15 and 363.15 K at various starting compositions covering (PODE1 + MeOH)/water molar ratio ranges from 0.5 to 1.5 in a time scale. On the basis of the experimental results, a mole fraction-based model of the chemical equilibrium and a pseudohomogeneous model are proposed to fit data based on true amount of monomeric formaldehyde. It has been demonstrated that the hydrolysis of PODE1 is slightly endothermic with the enthalpy 8.19 kJ/mol and the rate determining step. Finally, a feed-forward artificial neural networks (ANN) model is developed to model the concentration change of methanol in aqueous solutions. The results showed that the predicted data from designed ANN model were in good agreement with the experimental data with the coefficient ( R 2) of 0.98. Designed ANN provides a reliable method for modeling the hydrolysis reaction of methylal (PODE1).

Image texture segmentation using a neural network

NASA Astrophysics Data System (ADS)

Sayeh, Mohammed R.; Athinarayanan, Ragu; Dhali, Pushpuak

1992-09-01

In this paper we use a neural network called the Lyapunov associative memory (LYAM) system to segment image texture into different categories or clusters. The LYAM system is constructed by a set of ordinary differential equations which are simulated on a digital computer. The clustering can be achieved by using a single tuning parameter in the simplest model. Pattern classes are represented by the stable equilibrium states of the system. Design of the system is based on synthesizing two local energy functions, namely, the learning and recall energy functions. Before the implementation of the segmentation process, a Gauss-Markov random field (GMRF) model is applied to the raw image. This application suitably reduces the image data and prepares the texture information for the neural network process. We give a simple image example illustrating the capability of the technique. The GMRF-generated features are also used for a clustering, based on the Euclidean distance.

Dynamic Power Distribution System Management With a Locally Connected Communication Network

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dall-Anese, Emiliano; Zhang, Kaiqing; Basar, Tamer

Coordinated optimization and control of distribution-level assets can enable a reliable and optimal integration of massive amount of distributed energy resources (DERs) and facilitate distribution system management (DSM). Accordingly, the objective is to coordinate the power injection at the DERs to maintain certain quantities across the network, e.g., voltage magnitude, line flows, or line losses, to be close to a desired profile. By and large, the performance of the DSM algorithms has been challenged by two factors: i) the possibly non-strongly connected communication network over DERs that hinders the coordination; ii) the dynamics of the real system caused by themore » DERs with heterogeneous capabilities, time-varying operating conditions, and real-time measurement mismatches. In this paper, we investigate the modeling and algorithm design and analysis with the consideration of these two factors. In particular, a game theoretic characterization is first proposed to account for a locally connected communication network over DERs, along with the analysis of the existence and uniqueness of the Nash equilibrium (NE) therein. To achieve the equilibrium in a distributed fashion, a projected-gradient-based asynchronous DSM algorithm is then advocated. The algorithm performance, including the convergence speed and the tracking error, is analytically guaranteed under the dynamic setting. Extensive numerical tests on both synthetic and realistic cases corroborate the analytical results derived.« less

Critical behavior of the XY-rotor model on regular and small-world networks

NASA Astrophysics Data System (ADS)

De Nigris, Sarah; Leoncini, Xavier

2013-07-01

We study the XY rotors model on small networks whose number of links scales with the system size Nlinks˜Nγ, where 1≤γ≤2. We first focus on regular one-dimensional rings in the microcanonical ensemble. For γ<1.5 the model behaves like a short-range one and no phase transition occurs. For γ>1.5, the system equilibrium properties are found to be identical to the mean field, which displays a second-order phase transition at a critical energy density ɛ=E/N,ɛc=0.75. Moreover, for γc≃1.5 we find that a nontrivial state emerges, characterized by an infinite susceptibility. We then consider small-world networks, using the Watts-Strogatz mechanism on the regular networks parametrized by γ. We first analyze the topology and find that the small-world regime appears for rewiring probabilities which scale as pSW∝1/Nγ. Then considering the XY-rotors model on these networks, we find that a second-order phase transition occurs at a critical energy ɛc which logarithmically depends on the topological parameters p and γ. We also define a critical probability pMF, corresponding to the probability beyond which the mean field is quantitatively recovered, and we analyze its dependence on γ.

Non-equilibrium assembly of microtubules: from molecules to autonomous chemical robots.

PubMed

Hess, H; Ross, Jennifer L

2017-09-18

Biological systems have evolved to harness non-equilibrium processes from the molecular to the macro scale. It is currently a grand challenge of chemistry, materials science, and engineering to understand and mimic biological systems that have the ability to autonomously sense stimuli, process these inputs, and respond by performing mechanical work. New chemical systems are responding to the challenge and form the basis for future responsive, adaptive, and active materials. In this article, we describe a particular biochemical-biomechanical network based on the microtubule cytoskeletal filament - itself a non-equilibrium chemical system. We trace the non-equilibrium aspects of the system from molecules to networks and describe how the cell uses this system to perform active work in essential processes. Finally, we discuss how microtubule-based engineered systems can serve as testbeds for autonomous chemical robots composed of biological and synthetic components.

Network Polymers Formed Under Nonideal Conditions.

DTIC Science & Technology

1986-12-01

the system or the limited ability of the statistical model to account for stochastic correlations. The viscosity of the reacting system was measured as...based on competing reactions (ring, chain) and employs equilibrium chain statistics . The work thus far has been limited to single cycle growth on an...polymerizations, because a large number of differential equations must be solved. The Makovian approach (sometimes referred to as the statistical or

Global exponential stability of bidirectional associative memory neural networks with distributed delays

NASA Astrophysics Data System (ADS)

Song, Qiankun; Cao, Jinde

2007-05-01

A bidirectional associative memory neural network model with distributed delays is considered. By constructing a new Lyapunov functional, employing the homeomorphism theory, M-matrix theory and the inequality (a[greater-or-equal, slanted]0,bk[greater-or-equal, slanted]0,qk>0 with , and r>1), a sufficient condition is obtained to ensure the existence, uniqueness and global exponential stability of the equilibrium point for the model. Moreover, the exponential converging velocity index is estimated, which depends on the delay kernel functions and the system parameters. The results generalize and improve the earlier publications, and remove the usual assumption that the activation functions are bounded . Two numerical examples are given to show the effectiveness of the obtained results.

Modeling of knowledge transmission by considering the level of forgetfulness in complex networks

NASA Astrophysics Data System (ADS)

Cao, Bin; Han, Shui-hua; Jin, Zhen

2016-06-01

In this study, we establish a general model by considering the level of forgetfulness during knowledge transmission in complex networks, where the level of forgetfulness depends mainly on the number in a crowd who possess knowledge, while the saturated incidence is also considered. In theory, we analyze the stability of the equilibrium points and the transmission threshold R0 is also given. If R0 > 1, then knowledge can ​be transmitted, but if not, it will become completely extinct. In addition, we performed some numerical simulations to verify the reasonability of the theoretical analysis. The results of the simulations also suggest ​that the proportion of the crowd with knowledge will be increased under a better cultural atmosphere.

Epidemic spreading on random surfer networks with optimal interaction radius

NASA Astrophysics Data System (ADS)

Feng, Yun; Ding, Li; Hu, Ping

2018-03-01

In this paper, the optimal control problem of epidemic spreading on random surfer heterogeneous networks is considered. An epidemic spreading model is established according to the classification of individual's initial interaction radii. Then, a control strategy is proposed based on adjusting individual's interaction radii. The global stability of the disease free and endemic equilibrium of the model is investigated. We prove that an optimal solution exists for the optimal control problem and the explicit form of which is presented. Numerical simulations are conducted to verify the correctness of the theoretical results. It is proved that the optimal control strategy is effective to minimize the density of infected individuals and the cost associated with the adjustment of interaction radii.

Competitive energy consumption under transmission constraints in a multi-supplier power grid system

NASA Astrophysics Data System (ADS)

Popov, Ivan; Krylatov, Alexander; Zakharov, Victor; Ivanov, Dmitry

2017-04-01

Power grid architectures need to be revised in order to manage the increasing number of producers and, more generally, the decentralisation of energy production and distribution. In this work, we describe a multi-supplier multi-consumer congestion model of a power grid, where the costs of consumers depend on the congestion in nodes and arcs of the power supply network. The consumer goal is both to meet their energy demand and to minimise the costs. We show that the methods of non-atomic routing can be applied in this model in order to describe current distribution in the network. We formulate a consumer cost minimisation game for this setting, and discuss the challenges arising in equilibrium search for this game.

Unveiling the molecular mechanism of self-healing in a telechelic, supramolecular polymer network

PubMed Central

Yan, Tingzi; Schröter, Klaus; Herbst, Florian; Binder, Wolfgang H.; Thurn-Albrecht, Thomas

2016-01-01

Reversible polymeric networks can show self-healing properties due to their ability to reassemble after application of stress and fracture, but typically the relation between equilibrium molecular dynamics and self-healing kinetics has been difficult to disentangle. Here we present a well-characterized, self-assembled bulk network based on supramolecular assemblies, that allows a clear distinction between chain dynamics and network relaxation. Small angle x-ray scattering and rheological measurements provide evidence for a structurally well-defined, dense network of interconnected aggregates giving mechanical strength to the material. Different from a covalent network, the dynamic character of the supramolecular bonds enables macroscopic flow on a longer time scale and the establishment of an equilibrium structure. A combination of linear and nonlinear rheological measurements clearly identifies the terminal relaxation process as being responsible for the process of self-healing. PMID:27581380

Analytical theory of polymer-network-mediated interaction between colloidal particles

PubMed Central

Di Michele, Lorenzo; Zaccone, Alessio; Eiser, Erika

2012-01-01

Nanostructured materials based on colloidal particles embedded in a polymer network are used in a variety of applications ranging from nanocomposite rubbers to organic-inorganic hybrid solar cells. Further, polymer-network-mediated colloidal interactions are highly relevant to biological studies whereby polymer hydrogels are commonly employed to probe the mechanical response of living cells, which can determine their biological function in physiological environments. The performance of nanomaterials crucially relies upon the spatial organization of the colloidal particles within the polymer network that depends, in turn, on the effective interactions between the particles in the medium. Existing models based on nonlocal equilibrium thermodynamics fail to clarify the nature of these interactions, precluding the way toward the rational design of polymer-composite materials. In this article, we present a predictive analytical theory of these interactions based on a coarse-grained model for polymer networks. We apply the theory to the case of colloids partially embedded in cross-linked polymer substrates and clarify the origin of attractive interactions recently observed experimentally. Monte Carlo simulation results that quantitatively confirm the theoretical predictions are also presented. PMID:22679289

Modelling the influence of parental effects on gene-network evolution.

PubMed

Odorico, Andreas; Rünneburger, Estelle; Le Rouzic, Arnaud

2018-05-01

Understanding the importance of nongenetic heredity in the evolutionary process is a major topic in modern evolutionary biology. We modified a classical gene-network model by allowing parental transmission of gene expression and studied its evolutionary properties through individual-based simulations. We identified ontogenetic time (i.e. the time gene networks have to stabilize before being submitted to natural selection) as a crucial factor in determining the evolutionary impact of this phenotypic inheritance. Indeed, fast-developing organisms display enhanced adaptation and greater robustness to mutations when evolving in presence of nongenetic inheritance (NGI). In contrast, in our model, long development reduces the influence of the inherited state of the gene network. NGI thus had a negligible effect on the evolution of gene networks when the speed at which transcription levels reach equilibrium is not constrained. Nevertheless, simulations show that intergenerational transmission of the gene-network state negatively affects the evolution of robustness to environmental disturbances for either fast- or slow-developing organisms. Therefore, these results suggest that the evolutionary consequences of NGI might not be sought only in the way species respond to selection, but also on the evolution of emergent properties (such as environmental and genetic canalization) in complex genetic architectures. © 2018 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2018 European Society For Evolutionary Biology.

Self-organization of network dynamics into local quantized states

DOE PAGES

Nicolaides, Christos; Juanes, Ruben; Cueto-Felgueroso, Luis

2016-02-17

Self-organization and pattern formation in network-organized systems emerges from the collective activation and interaction of many interconnected units. A striking feature of these non-equilibrium structures is that they are often localized and robust: only a small subset of the nodes, or cell assembly, is activated. Understanding the role of cell assemblies as basic functional units in neural networks and socio-technical systems emerges as a fundamental challenge in network theory. A key open question is how these elementary building blocks emerge, and how they operate, linking structure and function in complex networks. Here we show that a network analogue of themore » Swift-Hohenberg continuum model—a minimal-ingredients model of nodal activation and interaction within a complex network—is able to produce a complex suite of localized patterns. Thus, the spontaneous formation of robust operational cell assemblies in complex networks can be explained as the result of self-organization, even in the absence of synaptic reinforcements.« less

Self-organization of network dynamics into local quantized states

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nicolaides, Christos; Juanes, Ruben; Cueto-Felgueroso, Luis

Self-organization and pattern formation in network-organized systems emerges from the collective activation and interaction of many interconnected units. A striking feature of these non-equilibrium structures is that they are often localized and robust: only a small subset of the nodes, or cell assembly, is activated. Understanding the role of cell assemblies as basic functional units in neural networks and socio-technical systems emerges as a fundamental challenge in network theory. A key open question is how these elementary building blocks emerge, and how they operate, linking structure and function in complex networks. Here we show that a network analogue of themore » Swift-Hohenberg continuum model—a minimal-ingredients model of nodal activation and interaction within a complex network—is able to produce a complex suite of localized patterns. Thus, the spontaneous formation of robust operational cell assemblies in complex networks can be explained as the result of self-organization, even in the absence of synaptic reinforcements.« less

Optimizing Cellular Networks Enabled with Renewal Energy via Strategic Learning.

PubMed

Sohn, Insoo; Liu, Huaping; Ansari, Nirwan

2015-01-01

An important issue in the cellular industry is the rising energy cost and carbon footprint due to the rapid expansion of the cellular infrastructure. Greening cellular networks has thus attracted attention. Among the promising green cellular network techniques, the renewable energy-powered cellular network has drawn increasing attention as a critical element towards reducing carbon emissions due to massive energy consumption in the base stations deployed in cellular networks. Game theory is a branch of mathematics that is used to evaluate and optimize systems with multiple players with conflicting objectives and has been successfully used to solve various problems in cellular networks. In this paper, we model the green energy utilization and power consumption optimization problem of a green cellular network as a pilot power selection strategic game and propose a novel distributed algorithm based on a strategic learning method. The simulation results indicate that the proposed algorithm achieves correlated equilibrium of the pilot power selection game, resulting in optimum green energy utilization and power consumption reduction.

Equilibrium and nonequilibrium properties of Boolean decision problems on scale-free graphs with competing interactions with external biases

NASA Astrophysics Data System (ADS)

Zhu, Zheng; Andresen, Juan Carlos; Janzen, Katharina; Katzgraber, Helmut G.

2013-03-01

We study the equilibrium and nonequilibrium properties of Boolean decision problems with competing interactions on scale-free graphs in a magnetic field. Previous studies at zero field have shown a remarkable equilibrium stability of Boolean variables (Ising spins) with competing interactions (spin glasses) on scale-free networks. When the exponent that describes the power-law decay of the connectivity of the network is strictly larger than 3, the system undergoes a spin-glass transition. However, when the exponent is equal to or less than 3, the glass phase is stable for all temperatures. First we perform finite-temperature Monte Carlo simulations in a field to test the robustness of the spin-glass phase and show, in agreement with analytical calculations, that the system exhibits a de Almeida-Thouless line. Furthermore, we study avalanches in the system at zero temperature to see if the system displays self-organized criticality. This would suggest that damage (avalanches) can spread across the whole system with nonzero probability, i.e., that Boolean decision problems on scale-free networks with competing interactions are fragile when not in thermal equilibrium.

Non-equilibrium electrokinetic micromixer with 3D nanochannel networks.

PubMed

Choi, Eunpyo; Kwon, Kilsung; Lee, Seung Jun; Kim, Daejoong; Park, Jungyul

2015-04-21

We report an active micromixer which utilizes vortex generation due to non-equilibrium electrokinetics near the interface between a microchannnel and a nanochannel networks membrane (NCNM), constructed from geometrically controlled in situ self-assembled nanoparticles. A large interfacing area where it is possible to generate vortices can be realized, because nano-interstices between the assembled nanoparticles are intrinsically collective three-dimensional nanochannel networks, which may be compared to typical silicon-based 2D nanochannels. The proposed mixer shows a 2-fold shorter mixing time (~0.78 ms) and a 34-fold shorter mixing length (~7.86 μm) compared to conventional 2D nanochannels.

Criticism of generally accepted fundamentals and methodologies of traffic and transportation theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kerner, Boris S.

It is explained why the set of the fundamental empirical features of traffic breakdown (a transition from free flow to congested traffic) should be the empirical basis for any traffic and transportation theory that can be reliable used for control and optimization in traffic networks. It is shown that generally accepted fundamentals and methodologies of traffic and transportation theory are not consistent with the set of the fundamental empirical features of traffic breakdown at a highway bottleneck. To these fundamentals and methodologies of traffic and transportation theory belong (i) Lighthill-Whitham-Richards (LWR) theory, (ii) the General Motors (GM) model class (formore » example, Herman, Gazis et al. GM model, Gipps’s model, Payne’s model, Newell’s optimal velocity (OV) model, Wiedemann’s model, Bando et al. OV model, Treiber’s IDM, Krauß’s model), (iii) the understanding of highway capacity as a particular stochastic value, and (iv) principles for traffic and transportation network optimization and control (for example, Wardrop’s user equilibrium (UE) and system optimum (SO) principles). Alternatively to these generally accepted fundamentals and methodologies of traffic and transportation theory, we discuss three-phase traffic theory as the basis for traffic flow modeling as well as briefly consider the network breakdown minimization (BM) principle for the optimization of traffic and transportation networks with road bottlenecks.« less

A non-equilibrium thermodynamic model for tumor extracellular matrix with enzymatic degradation

NASA Astrophysics Data System (ADS)

Xue, Shi-Lei; Li, Bo; Feng, Xi-Qiao; Gao, Huajian

2017-07-01

The extracellular matrix (ECM) of a solid tumor not only affords scaffolding to support tumor architecture and integrity but also plays an essential role in tumor growth, invasion, metastasis, and therapeutics. In this paper, a non-equilibrium thermodynamic theory is established to study the chemo-mechanical behaviors of tumor ECM, which is modeled as a poroelastic polyelectrolyte consisting of a collagen network and proteoglycans. By using the principle of maximum energy dissipation rate, we deduce a set of governing equations for drug transport and mechanosensitive enzymatic degradation in ECM. The results reveal that osmosis is primarily responsible for the compression resistance of ECM. It is suggested that a well-designed ECM degradation can effectively modify the tumor microenvironment for improved efficiency of cancer therapy. The theoretical predictions show a good agreement with relevant experimental observations. This study aimed to deepen our understanding of tumor ECM may be conducive to novel anticancer strategies.

Three-dimensional finite elements for the analysis of soil contamination using a multiple-porosity approach

NASA Astrophysics Data System (ADS)

El-Zein, Abbas; Carter, John P.; Airey, David W.

2006-06-01

A three-dimensional finite-element model of contaminant migration in fissured clays or contaminated sand which includes multiple sources of non-equilibrium processes is proposed. The conceptual framework can accommodate a regular network of fissures in 1D, 2D or 3D and immobile solutions in the macro-pores of aggregated topsoils, as well as non-equilibrium sorption. A Galerkin weighted-residual statement for the three-dimensional form of the equations in the Laplace domain is formulated. Equations are discretized using linear and quadratic prism elements. The system of algebraic equations is solved in the Laplace domain and solution is inverted to the time domain numerically. The model is validated and its scope is illustrated through the analysis of three problems: a waste repository deeply buried in fissured clay, a storage tank leaking into sand and a sanitary landfill leaching into fissured clay over a sand aquifer.

Wetting morphologies on randomly oriented fibers.

PubMed

Sauret, Alban; Boulogne, François; Soh, Beatrice; Dressaire, Emilie; Stone, Howard A

2015-06-01

We characterize the different morphologies adopted by a drop of liquid placed on two randomly oriented fibers, which is a first step toward understanding the wetting of fibrous networks. The present work reviews previous modeling for parallel and touching crossed fibers and extends it to an arbitrary orientation of the fibers characterized by the tilting angle and the minimum spacing distance. Depending on the volume of liquid, the spacing distance between fibers and the angle between the fibers, we highlight that the liquid can adopt three different equilibrium morphologies: 1) a column morphology in which the liquid spreads between the fibers, 2) a mixed morphology where a drop grows at one end of the column or 3) a single drop located at the node. We capture the different morphologies observed using an analytical model that predicts the equilibrium configuration of the liquid based on the geometry of the fibers and the volume of liquid.

Punctuated equilibrium as an emergent process and its modified thermodynamic characterization.

PubMed

Wosniack, M E; da Luz, M G E; Schulman, L S

2017-01-07

We address evolutionary dynamics and consider under which conditions the ecosystem interaction network allows punctuated equilibrium (i.e., alternation between hectic and quasi-stable phases). We focus on the links connecting various species and on the strength and sign of those links. For this study we consider the Tangled Nature model, which allows considerable flexibility and plasticity in the analysis of interspecies interactions. We find that it is necessary to have a proper balance of connectivity and interaction intensities so as to establish the kind of mutual cooperation and competition found in nature. It suggests evolutionary punctuated equilibrium as an emergent process, thus displaying features of complex systems. To explicitly demonstrate this fact we consider an extended form of thermodynamics, defining (for the present context) relevant out-of-equilibrium "collective" functions. We then show how to characterize the punctuated equilibrium through entropy-like and free energy-like quantities. Finally, from a close analogy to thermodynamic systems, we propose a protocol similar to simulated annealing. It is based on controlling the species' rate of mutation during the hectic periods, in this way enhancing the exploration of the genome space (similar to the known behavior of bacteria in stressful environments). This allows the system to more rapidly converge to long-duration quasi-stable phases. Copyright © 2016 Elsevier Ltd. All rights reserved.

A statistical mechanics approach to Granovetter theory

NASA Astrophysics Data System (ADS)

Barra, Adriano; Agliari, Elena

2012-05-01

In this paper we try to bridge breakthroughs in quantitative sociology/econometrics, pioneered during the last decades by Mac Fadden, Brock-Durlauf, Granovetter and Watts-Strogatz, by introducing a minimal model able to reproduce essentially all the features of social behavior highlighted by these authors. Our model relies on a pairwise Hamiltonian for decision-maker interactions which naturally extends the multi-populations approaches by shifting and biasing the pattern definitions of a Hopfield model of neural networks. Once introduced, the model is investigated through graph theory (to recover Granovetter and Watts-Strogatz results) and statistical mechanics (to recover Mac-Fadden and Brock-Durlauf results). Due to the internal symmetries of our model, the latter is obtained as the relaxation of a proper Markov process, allowing even to study its out-of-equilibrium properties. The method used to solve its equilibrium is an adaptation of the Hamilton-Jacobi technique recently introduced by Guerra in the spin-glass scenario and the picture obtained is the following: shifting the patterns from [-1,+1]→[0.+1] implies that the larger the amount of similarities among decision makers, the stronger their relative influence, and this is enough to explain both the different role of strong and weak ties in the social network as well as its small-world properties. As a result, imitative interaction strengths seem essentially a robust request (enough to break the gauge symmetry in the couplings), furthermore, this naturally leads to a discrete choice modelization when dealing with the external influences and to imitative behavior à la Curie-Weiss as the one introduced by Brock and Durlauf.

Dynamical Response of Networks Under External Perturbations: Exact Results

NASA Astrophysics Data System (ADS)

Chinellato, David D.; Epstein, Irving R.; Braha, Dan; Bar-Yam, Yaneer; de Aguiar, Marcus A. M.

2015-04-01

We give exact statistical distributions for the dynamic response of influence networks subjected to external perturbations. We consider networks whose nodes have two internal states labeled 0 and 1. We let nodes be frozen in state 0, in state 1, and the remaining nodes change by adopting the state of a connected node with a fixed probability per time step. The frozen nodes can be interpreted as external perturbations to the subnetwork of free nodes. Analytically extending and to be smaller than 1 enables modeling the case of weak coupling. We solve the dynamical equations exactly for fully connected networks, obtaining the equilibrium distribution, transition probabilities between any two states and the characteristic time to equilibration. Our exact results are excellent approximations for other topologies, including random, regular lattice, scale-free and small world networks, when the numbers of fixed nodes are adjusted to take account of the effect of topology on coupling to the environment. This model can describe a variety of complex systems, from magnetic spins to social networks to population genetics, and was recently applied as a framework for early warning signals for real-world self-organized economic market crises.

Specific non-monotonous interactions increase persistence of ecological networks.

PubMed

Yan, Chuan; Zhang, Zhibin

2014-03-22

The relationship between stability and biodiversity has long been debated in ecology due to opposing empirical observations and theoretical predictions. Species interaction strength is often assumed to be monotonically related to population density, but the effects on stability of ecological networks of non-monotonous interactions that change signs have not been investigated previously. We demonstrate that for four kinds of non-monotonous interactions, shifting signs to negative or neutral interactions at high population density increases persistence (a measure of stability) of ecological networks, while for the other two kinds of non-monotonous interactions shifting signs to positive interactions at high population density decreases persistence of networks. Our results reveal a novel mechanism of network stabilization caused by specific non-monotonous interaction types through either increasing stable equilibrium points or reducing unstable equilibrium points (or both). These specific non-monotonous interactions may be important in maintaining stable and complex ecological networks, as well as other networks such as genes, neurons, the internet and human societies.

Against the grain: The physical properties of anisotropic partially molten rocks

NASA Astrophysics Data System (ADS)

Ghanbarzadeh, S.; Hesse, M. A.; Prodanovic, M.

2014-12-01

Partially molten rocks commonly develop textures that appear close to textural equilibrium, where the melt network evolves to minimize the energy of the melt-solid interfaces, while maintaining the dihedral angle θ at solid-solid-melt contact lines. Textural equilibrium provides a powerful model for the melt distribution that controls the petro-physical properties of partially molten rocks, e.g., permeability, elastic moduli, and electrical resistivity. We present the first level-set computations of three-dimensional texturally equilibrated melt networks in rocks with an anisotropic fabric. Our results show that anisotropy induces wetting of smaller grain boundary faces for θ > 0 at realistic porosities ϕ < 3%. This was previously not thought to be possible at textural equilibrium and reconciles the theory with experimental observations. Wetting of the grain boundary faces leads to a dramatic redistribution of the melt from the edges to the faces that introduces strong anisotropy in the petro-physical properties such as permeability, effective electrical conductivity and mechanical properties. Figure, on left, shows that smaller grain boundaries become wetted at relatively low melt fractions of 3% in stretched polyhedral grains with elongation factor 1.5. Right plot represents the ratio of melt electrical conductivity to effective conductivity of medium (known as formation factor) as an example of anisotropy in physical properties. The plot shows that even slight anisotropy in grains induces considerable anisotropy in electrical properties.

NetFlow Dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Corbet Jr., Thomas F; Beyeler, Walter E; Vanwestrienen, Dirk

NetFlow Dynamics is a web-accessible analysis environment for simulating dynamic flows of materials on model networks. Performing a simulation requires both the NetFlow Dynamics application and a network model which is a description of the structure of the nodes and edges of a network including the flow capacity of each edge and the storage capacity of each node, and the sources and sinks of the material flowing on the network. NetFlow Dynamics consists of databases for storing network models, algorithms to calculate flows on networks, and a GIS-based graphical interface for performing simulations and viewing simulation results. Simulated flows aremore » dynamic in the sense that flows on each edge of the network and inventories at each node change with time and can be out of equilibrium with boundary conditions. Any number of network models could be simulated using Net Flow Dynamics. To date, the models simulated have been models of petroleum infrastructure. The main model has been the National Transportation Fuels Model (NTFM), a network of U.S. oil fields, transmission pipelines, rail lines, refineries, tank farms, and distribution terminals. NetFlow Dynamics supports two different flow algorithms, the Gradient Flow algorithm and the Inventory Control algorithm, that were developed specifically for the NetFlow Dynamics application. The intent is to add additional algorithms in the future as needed. The ability to select from multiple algorithms is desirable because a single algorithm never covers all analysis needs. The current algorithms use a demand-driven capacity-constrained formulation which means that the algorithms strive to use all available capacity and stored inventory to meet desired flows to sinks, subject to the capacity constraints of each network component. The current flow algorithms are best suited for problems in which a material flows on a capacity-constrained network representing a supply chain in which the material supplied can be stored at each node of the network. In the petroleum models, the flowing materials are crude oil and refined products that can be stored at tank farms, refineries, or terminals (i.e. the nodes of the network). Examples of other network models that could be simulated are currency flowing in a financial network, agricultural products moving to market, or natural gas flowing on a pipeline network.« less

Temperature-Dependent Estimation of Gibbs Energies Using an Updated Group-Contribution Method.

PubMed

Du, Bin; Zhang, Zhen; Grubner, Sharon; Yurkovich, James T; Palsson, Bernhard O; Zielinski, Daniel C

2018-06-05

Reaction-equilibrium constants determine the metabolite concentrations necessary to drive flux through metabolic pathways. Group-contribution methods offer a way to estimate reaction-equilibrium constants at wide coverage across the metabolic network. Here, we present an updated group-contribution method with 1) additional curated thermodynamic data used in fitting and 2) capabilities to calculate equilibrium constants as a function of temperature. We first collected and curated aqueous thermodynamic data, including reaction-equilibrium constants, enthalpies of reaction, Gibbs free energies of formation, enthalpies of formation, entropy changes of formation of compounds, and proton- and metal-ion-binding constants. Next, we formulated the calculation of equilibrium constants as a function of temperature and calculated the standard entropy change of formation (Δ f S ∘ ) using a model based on molecular properties. The median absolute error in estimating Δ f S ∘ was 0.013 kJ/K/mol. We also estimated magnesium binding constants for 618 compounds using a linear regression model validated against measured data. We demonstrate the improved performance of the current method (8.17 kJ/mol in median absolute residual) over the current state-of-the-art method (11.47 kJ/mol) in estimating the 185 new reactions added in this work. The efforts here fill in gaps for thermodynamic calculations under various conditions, specifically different temperatures and metal-ion concentrations. These, to our knowledge, new capabilities empower the study of thermodynamic driving forces underlying the metabolic function of organisms living under diverse conditions. Copyright © 2018 Biophysical Society. Published by Elsevier Inc. All rights reserved.

The feasibility and stability of large complex biological networks: a random matrix approach.

PubMed

Stone, Lewi

2018-05-29

In the 70's, Robert May demonstrated that complexity creates instability in generic models of ecological networks having random interaction matrices A. Similar random matrix models have since been applied in many disciplines. Central to assessing stability is the "circular law" since it describes the eigenvalue distribution for an important class of random matrices A. However, despite widespread adoption, the "circular law" does not apply for ecological systems in which density-dependence operates (i.e., where a species growth is determined by its density). Instead one needs to study the far more complicated eigenvalue distribution of the community matrix S = DA, where D is a diagonal matrix of population equilibrium values. Here we obtain this eigenvalue distribution. We show that if the random matrix A is locally stable, the community matrix S = DA will also be locally stable, providing the system is feasible (i.e., all species have positive equilibria D > 0). This helps explain why, unusually, nearly all feasible systems studied here are locally stable. Large complex systems may thus be even more fragile than May predicted, given the difficulty of assembling a feasible system. It was also found that the degree of stability, or resilience of a system, depended on the minimum equilibrium population.

Coexistence and local μ-stability of multiple equilibrium points for memristive neural networks with nonmonotonic piecewise linear activation functions and unbounded time-varying delays.

PubMed

Nie, Xiaobing; Zheng, Wei Xing; Cao, Jinde

2016-12-01

In this paper, the coexistence and dynamical behaviors of multiple equilibrium points are discussed for a class of memristive neural networks (MNNs) with unbounded time-varying delays and nonmonotonic piecewise linear activation functions. By means of the fixed point theorem, nonsmooth analysis theory and rigorous mathematical analysis, it is proven that under some conditions, such n-neuron MNNs can have 5 n equilibrium points located in ℜ n , and 3 n of them are locally μ-stable. As a direct application, some criteria are also obtained on the multiple exponential stability, multiple power stability, multiple log-stability and multiple log-log-stability. All these results reveal that the addressed neural networks with activation functions introduced in this paper can generate greater storage capacity than the ones with Mexican-hat-type activation function. Numerical simulations are presented to substantiate the theoretical results. Copyright © 2016 Elsevier Ltd. All rights reserved.

Qualitative analysis of Cohen-Grossberg neural networks with multiple delays

NASA Astrophysics Data System (ADS)

Ye, Hui; Michel, Anthony N.; Wang, Kaining

1995-03-01

It is well known that a class of artificial neural networks with symmetric interconnections and without transmission delays, known as Cohen-Grossberg neural networks, possesses global stability (i.e., all trajectories tend to some equilibrium). We demonstrate in the present paper that many of the qualitative properties of Cohen-Grossberg networks will not be affected by the introduction of sufficiently small delays. Specifically, we establish some bound conditions for the time delays under which a given Cohen-Grossberg network with multiple delays is globally stable and possesses the same asymptotically stable equilibria as the corresponding network without delays. An effective method of determining the asymptotic stability of an equilibrium of a Cohen-Grossberg network with multiple delays is also presented. The present results are motivated by some of the authors earlier work [Phys. Rev. E 50, 4206 (1994)] and by some of the work of Marcus and Westervelt [Phys. Rev. A 39, 347 (1989)]. These works address qualitative analyses of Hopfield neural networks with one time delay. The present work generalizes these results to Cohen-Grossberg neural networks with multiple time delays. Hopfield neural networks constitute special cases of Cohen-Grossberg neural networks.

Primitive-path statistics of entangled polymers: mapping multi-chain simulations onto single-chain mean-field models

NASA Astrophysics Data System (ADS)

Steenbakkers, Rudi J. A.; Tzoumanekas, Christos; Li, Ying; Liu, Wing Kam; Kröger, Martin; Schieber, Jay D.

2014-01-01

We present a method to map the full equilibrium distribution of the primitive-path (PP) length, obtained from multi-chain simulations of polymer melts, onto a single-chain mean-field ‘target’ model. Most previous works used the Doi-Edwards tube model as a target. However, the average number of monomers per PP segment, obtained from multi-chain PP networks, has consistently shown a discrepancy of a factor of two with respect to tube-model estimates. Part of the problem is that the tube model neglects fluctuations in the lengths of PP segments, the number of entanglements per chain and the distribution of monomers among PP segments, while all these fluctuations are observed in multi-chain simulations. Here we use a recently proposed slip-link model, which includes fluctuations in all these variables as well as in the spatial positions of the entanglements. This turns out to be essential to obtain qualitative and quantitative agreement with the equilibrium PP-length distribution obtained from multi-chain simulations. By fitting this distribution, we are able to determine two of the three parameters of the model, which govern its equilibrium properties. This mapping is executed for four different linear polymers and for different molecular weights. The two parameters are found to depend on chemistry, but not on molecular weight. The model predicts a constant plateau modulus minus a correction inversely proportional to molecular weight. The value for well-entangled chains, with the parameters determined ab initio, lies in the range of experimental data for the materials investigated.

Simplification of reversible Markov chains by removal of states with low equilibrium occupancy.

PubMed

Ullah, Ghanim; Bruno, William J; Pearson, John E

2012-10-21

We present a practical method for simplifying Markov chains on a potentially large state space when detailed balance holds. A simple and transparent technique is introduced to remove states with low equilibrium occupancy. The resulting system has fewer parameters. The resulting effective rates between the remaining nodes give dynamics identical to the original system's except on very fast timescales. This procedure amounts to using separation of timescales to neglect small capacitance nodes in a network of resistors and capacitors. We illustrate the technique by simplifying various reaction networks, including transforming an acyclic four-node network to a three-node cyclic network. For a reaction step in which a ligand binds, the law of mass action implies a forward rate proportional to ligand concentration. The effective rates in the simplified network are found to be rational functions of ligand concentration. Copyright © 2012 Elsevier Ltd. All rights reserved.

Mechanical response of transient telechelic networks with many-part stickers

NASA Astrophysics Data System (ADS)

Sing, Michelle K.; Ramírez, Jorge; Olsen, Bradley D.

2017-11-01

A central question in soft matter is understanding how several individual, weak bonds act together to produce collective interactions. Here, gel-forming telechelic polymers with multiple stickers at each chain end are studied through Brownian dynamics simulations to understand how collective interaction of the bonds affects mechanical response of the gels. These polymers are modeled as finitely extensible dumbbells using an explicit tau-leap algorithm and the binding energy of these associations was kept constant regardless of the number of stickers. The addition of multiple bonds to the associating ends of telechelic polymers increases or decreases the network relaxation time depending on the relative kinetics of association but increases both shear stress and extensional viscosity. The relationship between the rate of association and the Rouse time of dangling chains results in two different regimes for the equilibrium stress relaxation of associating physical networks. In case I, a dissociated dangling chain is able to fully relax before re-associating to the network, resulting in two characteristic relaxation times and a non-monotonic terminal relaxation time with increasing number of bonds per polymer endgroup. In case II, the dissociated dangling chain is only able to relax a fraction of the way before it re-attaches to the network, and increasing the number of bonds per endgroup monotonically increases the terminal relaxation time. In flow, increasing the number of stickers increases the steady-state shear and extensional viscosities even though the overall bond kinetics and equilibrium constant remain unchanged. Increased dissipation in the simulations is primarily due to higher average chain extension with increasing bond number. These results indicate that toughness and dissipation in physically associating networks can both be increased by breaking single, strong bonds into smaller components.

Microscale models of partially molten rocks and their macroscale physical properties

NASA Astrophysics Data System (ADS)

Rudge, J. F.

2017-12-01

Any geodynamical model of melt transport in the Earth's mantle requires constitutive laws for the rheology of partially molten rock. These constitutive laws are poorly known, and one way to make progress in our understanding is through the upscaling of microscale models which describe physics at the scale of individual mineral grains. Crucially, many upscaled physical properties (such as permeability) depend not only on how much melt is present, but on how that melt is arranged at the microscale; i.e. on the geometry of the melt network. Here I will present some new calculations of equilibrium melt network geometries around idealised tetrakaidecahedral grains. In contrast to several previous calculations of textural equilibrium, these calculations allow for a both a liquid-phase and a solid-phase topology that can tile 3D space. The calculations are based on a simple minimisation of surface energy using the finite element method. In these simple models just two parameters control the topology of the melt network: the porosity (volume fraction of melt), and the dihedral angle. The consquences of these melt geometries for upscaled properties such as permeability; electrical conductivity; and importantly, effective viscosity will be explored. Recent theoretical work [1,2] has suggested that in diffusion creep a small amount of melt may dramatically reduce the effective shear viscosity of a partially molten rock, with profound consequences for the nature of the asthenosphere. This contribution will show that this reduction in viscosity may have been significantly overestimated, so that the drop in the effective viscosity at onset of melting is more modest. [1] Takei, Y., and B. K. Holtzman (2009), Viscous constitutive relations of solid-liquid composites in terms of grain boundary contiguity: 1. Grain boundary diffusion control model, J. Geophys. Res., 114, B06205.[2] Holtzmann B. K. (2016) Questions on the existence, persistence, and mechanical effects of a very small melt fraction in the asthenosphere, Geophys. Geochem. Geosyst. 17, 470-484.

Topological and kinetic determinants of the modal matrices of dynamic models of metabolism

PubMed Central

2017-01-01

Large-scale kinetic models of metabolism are becoming increasingly comprehensive and accurate. A key challenge is to understand the biochemical basis of the dynamic properties of these models. Linear analysis methods are well-established as useful tools for characterizing the dynamic response of metabolic networks. Central to linear analysis methods are two key matrices: the Jacobian matrix (J) and the modal matrix (M-1) arising from its eigendecomposition. The modal matrix M-1 contains dynamically independent motions of the kinetic model near a reference state, and it is sparse in practice for metabolic networks. However, connecting the structure of M-1 to the kinetic properties of the underlying reactions is non-trivial. In this study, we analyze the relationship between J, M-1, and the kinetic properties of the underlying network for kinetic models of metabolism. Specifically, we describe the origin of mode sparsity structure based on features of the network stoichiometric matrix S and the reaction kinetic gradient matrix G. First, we show that due to the scaling of kinetic parameters in real networks, diagonal dominance occurs in a substantial fraction of the rows of J, resulting in simple modal structures with clear biological interpretations. Then, we show that more complicated modes originate from topologically-connected reactions that have similar reaction elasticities in G. These elasticities represent dynamic equilibrium balances within reactions and are key determinants of modal structure. The work presented should prove useful towards obtaining an understanding of the dynamics of kinetic models of metabolism, which are rooted in the network structure and the kinetic properties of reactions. PMID:29267329

Emergence of network structure due to spike-timing-dependent plasticity in recurrent neuronal networks IV: structuring synaptic pathways among recurrent connections.

PubMed

Gilson, Matthieu; Burkitt, Anthony N; Grayden, David B; Thomas, Doreen A; van Hemmen, J Leo

2009-12-01

In neuronal networks, the changes of synaptic strength (or weight) performed by spike-timing-dependent plasticity (STDP) are hypothesized to give rise to functional network structure. This article investigates how this phenomenon occurs for the excitatory recurrent connections of a network with fixed input weights that is stimulated by external spike trains. We develop a theoretical framework based on the Poisson neuron model to analyze the interplay between the neuronal activity (firing rates and the spike-time correlations) and the learning dynamics, when the network is stimulated by correlated pools of homogeneous Poisson spike trains. STDP can lead to both a stabilization of all the neuron firing rates (homeostatic equilibrium) and a robust weight specialization. The pattern of specialization for the recurrent weights is determined by a relationship between the input firing-rate and correlation structures, the network topology, the STDP parameters and the synaptic response properties. We find conditions for feed-forward pathways or areas with strengthened self-feedback to emerge in an initially homogeneous recurrent network.

Frustration in protein elastic network models

NASA Astrophysics Data System (ADS)

Lezon, Timothy; Bahar, Ivet

2010-03-01

Elastic network models (ENMs) are widely used for studying the equilibrium dynamics of proteins. The most common approach in ENM analysis is to adopt a uniform force constant or a non-specific distance dependent function to represent the force constant strength. Here we discuss the influence of sequence and structure in determining the effective force constants between residues in ENMs. Using a novel method based on entropy maximization, we optimize the force constants such that they exactly reporduce a subset of experimentally determined pair covariances for a set of proteins. We analyze the optimized force constants in terms of amino acid types, distances, contact order and secondary structure, and we demonstrate that including frustrated interactions in the ENM is essential for accurately reproducing the global modes in the middle of the frequency spectrum.

Effects of rewiring strategies on information spreading in complex dynamic networks

NASA Astrophysics Data System (ADS)

Ally, Abdulla F.; Zhang, Ning

2018-04-01

Recent advances in networks and communication services have attracted much interest to understand information spreading in social networks. Consequently, numerous studies have been devoted to provide effective and accurate models for mimicking information spreading. However, knowledge on how to spread information faster and more widely remains a contentious issue. Yet, most existing works are based on static networks which limit the reality of dynamism of entities that participate in information spreading. Using the SIR epidemic model, this study explores and compares effects of two rewiring models (Fermi-Dirac and Linear functions) on information spreading in scale free and small world networks. Our results show that for all the rewiring strategies, the spreading influence replenishes with time but stabilizes in a steady state at later time-steps. This means that information spreading takes-off during the initial spreading steps, after which the spreading prevalence settles toward its equilibrium, with majority of the population having recovered and thus, no longer affecting the spreading. Meanwhile, rewiring strategy based on Fermi-Dirac distribution function in one way or another impedes the spreading process, however, the structure of the networks mimic the spreading, even with a low spreading rate. The worst case can be when the spreading rate is extremely small. The results emphasize that despite a big role of such networks in mimicking the spreading, the role of the parameters cannot be simply ignored. Apparently, the probability of giant degree neighbors being informed grows much faster with the rewiring strategy of linear function compared to that of Fermi-Dirac distribution function. Clearly, rewiring model based on linear function generates the fastest spreading across the networks. Therefore, if we are interested in speeding up the spreading process in stochastic modeling, linear function may play a pivotal role.

A Statistical Test of Walrasian Equilibrium by Means of Complex Networks Theory

NASA Astrophysics Data System (ADS)

Bargigli, Leonardo; Viaggiu, Stefano; Lionetto, Andrea

2016-10-01

We represent an exchange economy in terms of statistical ensembles for complex networks by introducing the concept of market configuration. This is defined as a sequence of nonnegative discrete random variables {w_{ij}} describing the flow of a given commodity from agent i to agent j. This sequence can be arranged in a nonnegative matrix W which we can regard as the representation of a weighted and directed network or digraph G. Our main result consists in showing that general equilibrium theory imposes highly restrictive conditions upon market configurations, which are in most cases not fulfilled by real markets. An explicit example with reference to the e-MID interbank credit market is provided.

Modeling household and community transmission of Ebola virus disease: Epidemic growth, spatial dynamics and insights for epidemic control

PubMed Central

Kiskowski, Maria; Chowell, Gerardo

2016-01-01

The mechanisms behind the sub-exponential growth dynamics of the West Africa Ebola virus disease epidemic could be related to improved control of the epidemic and the result of reduced disease transmission in spatially constrained contact structures. An individual-based, stochastic network model is used to model immediate and delayed epidemic control in the context of social contact networks and investigate the extent to which the relative role of these factors may be determined during an outbreak. We find that in general, epidemics quickly establish a dynamic equilibrium of infections in the form of a wave of fixed size and speed traveling through the contact network. Both greater epidemic control and limited community mixing decrease the size of an infectious wave. However, for a fixed wave size, epidemic control (in contrast with limited community mixing) results in lower community saturation and a wave that moves more quickly through the contact network. We also found that the level of epidemic control has a disproportionately greater reductive effect on larger waves, so that a small wave requires nearly as much epidemic control as a larger wave to end an epidemic. PMID:26399855

Modeling household and community transmission of Ebola virus disease: Epidemic growth, spatial dynamics and insights for epidemic control.

PubMed

Kiskowski, Maria; Chowell, Gerardo

2016-01-01

The mechanisms behind the sub-exponential growth dynamics of the West Africa Ebola virus disease epidemic could be related to improved control of the epidemic and the result of reduced disease transmission in spatially constrained contact structures. An individual-based, stochastic network model is used to model immediate and delayed epidemic control in the context of social contact networks and investigate the extent to which the relative role of these factors may be determined during an outbreak. We find that in general, epidemics quickly establish a dynamic equilibrium of infections in the form of a wave of fixed size and speed traveling through the contact network. Both greater epidemic control and limited community mixing decrease the size of an infectious wave. However, for a fixed wave size, epidemic control (in contrast with limited community mixing) results in lower community saturation and a wave that moves more quickly through the contact network. We also found that the level of epidemic control has a disproportionately greater reductive effect on larger waves, so that a small wave requires nearly as much epidemic control as a larger wave to end an epidemic.

Rapid and high-capacity ultrasonic assisted adsorption of ternary toxic anionic dyes onto MOF-5-activated carbon: Artificial neural networks, partial least squares, desirability function and isotherm and kinetic study.

PubMed

Askari, Hanieh; Ghaedi, Mehrorang; Dashtian, Kheibar; Azghandi, Mohammad Hossein Ahmadi

2017-07-01

The present paper focused on the ultrasonic assisted simultaneous removal of fast green (FG), eosin Y (EY) and quinine yellow (QY) from aqueous media following using MOF-5 as a metal organic framework and activated carbon hybrid (AC-MOF-5). The structure and morphology of AC-MOF-5 was identified by SEM, FTIR and XRD analysis. The interactive and main effects of variables such as pH, initial dyes concentration (mgL -1 ), adsorbent dosage (mg) and sonication time (min) on removal percentage were studied by central composite design (CCD), subsequent desirability function (DF) permit to achieved real variable experimental condition. Optimized values were found 7.06, 5.68, 7.59 and 5.04mgL -1 , 0.02g and 2.55min for pH, FG, EY and QY concentration, adsorbent dosage and sonication time, respectively. Under this conditions removal percentage were obtained 98.1%, 98.1% and 91.91% for FG, EY and QY, respectively. Two models, namely partial least squares (PLS) and multi-layer artificial neural network (ANN) model were used for building up to construct an empirical model to predict the dyes under study removal behavior. The obtained results show that ANN and PLS model is a powerful tool for prediction of under-study dyes adsorption by AC-MOF-5. The evaluation and estimation of equilibrium data from traditional isotherm models display that the Langmuir model indicated the best fit to the equilibrium data with maximum adsorption capacity of 21.230, 20.242 and 18.621mgg -1 , for FG, EY and QY, respectively, while the adsorption rate efficiently follows the pseudo-second-order model. Copyright © 2016 Elsevier B.V. All rights reserved.

The synthesis and characterization of environmentally-responsive water-swellable and water-soluble polymers for wastewater remediation

NASA Astrophysics Data System (ADS)

Armentrout, Rodney Scott

The primary research goal is the development of new polymeric materials that demonstrate the environmentally-responsive sequestration of common water foulants, including surfactants and oils. Water-swellable and water-soluble polymers have been synthesized, structurally characterized, and their physical properties have been determined. In addition, the ability of the materials to sequester model water foulants has been evaluated. Anionic crosslinked polymer networks of 2-acrylamido-2-methyl-1-propanesulfonic acid, acrylamide, and methylene bisacrylamide have been synthesized and characterized by determining the equilibrium water contents as a function of ionic content of the polymer network. The molar ratio of bound surfactant to ionic group was determined to be less than one for all hydrogels studied, indicating an ion-exchange binding mechanism with minimal hydrophobic interactions between bound and unbound surfactant molecules is responsible for surfactant binding. Cationic crosslinked cyclopolymer networks of N,N-diallyl- N-methyl amine (DAMA) and N,N,N,N-tetraallyl ammonium chloride (TAAC) have been synthesized and characterized by determining the equilibrium water content as a function of pH. A maximum in the equilibrium water content is observed for pH-6 when the polymer is fully ionized. The solubilization of a model water foulant, p-cresol, by the polymeric surfactant, Pluronic F127, has been studied via equilibrium dialysis, dynamic light scattering and ultrafiltration experiments. It has been shown that at 25°C p-cresol is readily solubilized by F127 since the polymeric surfactant exists in a multimer conformation. Ultrafiltration experiments have demonstrated that the polymer-foulant binding interactions are largely unaffected by shear in a hollow fiber membrane. Copolymers of the zwitterionic monomer, 3-(N,N-diallyl- N-methyl ammonio) propane sulfonate (DAMAPS) and N,N-diallyl- N,N-dimethylammonium chloride (DADMAC) (the DADS series) or the pH-responsive hydrophobic monomer, N,N-diallyl-N-methyl amine (DAMA) (the DAMS series) have been prepared in a 0.5 M NaCl aqueous solution using 2-hydroxy-1-[4-(hydroxy-ethoxy)phenyl]-2-methyl-1-propanone (Irgacure 2959) as the free-radical photoinitiator. 13C NMR data indicate that the resulting polymers maintain the five-membered ring structure in the cis conformation common to diallylammonium salts. Equilibrium dialysis experiments demonstrate that pH-responsive hydrophobic microdomain formation may be utilized to control the solubilization of the organic solute, p-cresol. Ultrafiltration experiments have demonstrated that the polymer-foulant binding interactions are largely unaffected by shear in a hollow fiber membrane. Macromolecular aggregates of the poly( N,N-diallyl-N-methyl amine)/p-cresol complexes lead to fouling of the ultrafiltration membrane. However, incorporation of the sulfobetaine moiety hinders the formation of the macroscopic structures and higher permeate flux rates are achieved. (Abstract shortened by UMI.)

Global stability and optimal control of epidemic spreading on multiplex networks with nonlinear mutual interaction

NASA Astrophysics Data System (ADS)

Jia, Nan; Ding, Li; Liu, Yu-Jing; Hu, Ping

2018-07-01

In this paper, we consider two interacting pathogens spreading on multiplex networks. Each pathogen spreads only on its single layer, and different layers have the same individuals but different network topology. A state-dependent infectious rate is proposed to describe the nonlinear mutual interaction during the propagation of two pathogens. Then a novel epidemic spreading model incorporating treatment control strategy is established. We investigate the global asymptotic stability of the equilibrium points by using Dulac's criterion, Poincaré-Bendixson theorem and Lyapunov method. Furthermore, we discuss an optimal strategy to minimize the total number of the infected individuals and the cost associated with treatment control for both spreading of two pathogens. Finally, numerical simulations are presented to show the validity and efficiency of our results.

Increasing Resilience to Traumatic Stress: Understanding the Protective Role of Well-Being.

PubMed

Tory Toole, J; Rice, Mark A; Cargill, Jordan; Craddock, Travis J A; Nierenberg, Barry; Klimas, Nancy G; Fletcher, Mary Ann; Morris, Mariana; Zysman, Joel; Broderick, Gordon

2018-01-01

The brain maintains homeostasis in part through a network of feedback and feed-forward mechanisms, where neurochemicals and immune markers act as mediators. Using a previously constructed model of biobehavioral feedback, we found that in addition to healthy equilibrium another stable regulatory program supported chronic depression and anxiety. Exploring mechanisms that might underlie the contributions of subjective well-being to improved therapeutic outcomes in depression, we iteratively screened 288 candidate feedback patterns linking well-being to molecular signaling networks for those that maintained the original homeostatic regimes. Simulating stressful trigger events on each candidate network while maintaining high levels of subjective well-being isolated a specific feedback network where well-being was promoted by dopamine and acetylcholine, and itself promoted norepinephrine while inhibiting cortisol expression. This biobehavioral feedback mechanism was especially effective in reproducing well-being's clinically documented ability to promote resilience and protect against onset of depression and anxiety.

An annealed chaotic maximum neural network for bipartite subgraph problem.

PubMed

Wang, Jiahai; Tang, Zheng; Wang, Ronglong

2004-04-01

In this paper, based on maximum neural network, we propose a new parallel algorithm that can help the maximum neural network escape from local minima by including a transient chaotic neurodynamics for bipartite subgraph problem. The goal of the bipartite subgraph problem, which is an NP- complete problem, is to remove the minimum number of edges in a given graph such that the remaining graph is a bipartite graph. Lee et al. presented a parallel algorithm using the maximum neural model (winner-take-all neuron model) for this NP- complete problem. The maximum neural model always guarantees a valid solution and greatly reduces the search space without a burden on the parameter-tuning. However, the model has a tendency to converge to a local minimum easily because it is based on the steepest descent method. By adding a negative self-feedback to the maximum neural network, we proposed a new parallel algorithm that introduces richer and more flexible chaotic dynamics and can prevent the network from getting stuck at local minima. After the chaotic dynamics vanishes, the proposed algorithm is then fundamentally reined by the gradient descent dynamics and usually converges to a stable equilibrium point. The proposed algorithm has the advantages of both the maximum neural network and the chaotic neurodynamics. A large number of instances have been simulated to verify the proposed algorithm. The simulation results show that our algorithm finds the optimum or near-optimum solution for the bipartite subgraph problem superior to that of the best existing parallel algorithms.

Flow assignment model for quantitative analysis of diverting bulk freight from road to railway

PubMed Central

Liu, Chang; Wang, Jiaxi; Xiao, Jie; Liu, Siqi; Wu, Jianping; Li, Jian

2017-01-01

Since railway transport possesses the advantage of high volume and low carbon emissions, diverting some freight from road to railway will help reduce the negative environmental impacts associated with transport. This paper develops a flow assignment model for quantitative analysis of diverting truck freight to railway. First, a general network which considers road transportation, railway transportation, handling and transferring is established according to all the steps in the whole transportation process. Then general functions which embody the factors which the shippers will pay attention to when choosing mode and path are formulated. The general functions contain the congestion cost on road, the capacity constraints of railways and freight stations. Based on the general network and general cost function, a user equilibrium flow assignment model is developed to simulate the flow distribution on the general network under the condition that all shippers choose transportation mode and path independently. Since the model is nonlinear and challenging, we adopt a method that uses tangent lines to constitute envelope curve to linearize it. Finally, a numerical example is presented to test the model and show the method of making quantitative analysis of bulk freight modal shift between road and railway. PMID:28771536

A Method for Assessing the Retention of Trace Elements in Human Body Using Neural Network Technology

PubMed Central

Ragimov, Aligejdar; Faizullin, Rashat; Valiev, Vsevolod

2017-01-01

Models that describe the trace element status formation in the human organism are essential for a correction of micromineral (trace elements) deficiency. A direct trace element retention assessment in the body is difficult due to the many internal mechanisms. The trace element retention is determined by the amount and the ratio of incoming and excreted substance. So, the concentration of trace elements in drinking water characterizes the intake, whereas the element concentration in urine characterizes the excretion. This system can be interpreted as three interrelated elements that are in equilibrium. Since many relationships in the system are not known, the use of standard mathematical models is difficult. The artificial neural network use is suitable for constructing a model in the best way because it can take into account all dependencies in the system implicitly and process inaccurate and incomplete data. We created several neural network models to describe the retentions of trace elements in the human body. On the model basis, we can calculate the microelement levels in the body, knowing the trace element levels in drinking water and urine. These results can be used in health care to provide the population with safe drinking water. PMID:29065586

Generalized rules for the optimization of elastic network models

NASA Astrophysics Data System (ADS)

Lezon, Timothy; Eyal, Eran; Bahar, Ivet

2009-03-01

Elastic network models (ENMs) are widely employed for approximating the coarse-grained equilibrium dynamics of proteins using only a few parameters. An area of current focus is improving the predictive accuracy of ENMs by fine-tuning their force constants to fit specific systems. Here we introduce a set of general rules for assigning ENM force constants to residue pairs. Using a novel method, we construct ENMs that optimally reproduce experimental residue covariances from NMR models of 68 proteins. We analyze the optimal interactions in terms of amino acid types, pair distances and local protein structures to identify key factors in determining the effective spring constants. When applied to several unrelated globular proteins, our method shows an improved correlation with experiment over a standard ENM. We discuss the physical interpretation of our findings as well as its implications in the fields of protein folding and dynamics.

Self-organization in cold atomic gases: a synchronization perspective.

PubMed

Tesio, E; Robb, G R M; Oppo, G-L; Gomes, P M; Ackemann, T; Labeyrie, G; Kaiser, R; Firth, W J

2014-10-28

We study non-equilibrium spatial self-organization in cold atomic gases, where long-range spatial order spontaneously emerges from fluctuations in the plane transverse to the propagation axis of a single optical beam. The self-organization process can be interpreted as a synchronization transition in a fully connected network of fictitious oscillators, and described in terms of the Kuramoto model. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

Are equilibrium multichannel networks predictable? The case of the regulated Indus River, Pakistan

NASA Astrophysics Data System (ADS)

Carling, P. A.; Trieu, H.; Hornby, D. D.; Huang, He Qing; Darby, S. E.; Sear, D. A.; Hutton, C.; Hill, C.; Ali, Z.; Ahmed, A.; Iqbal, I.; Hussain, Z.

2018-02-01

Arguably, the current planform behaviour of the Indus River is broadly predictable. Between Chashma and Taunsa, Pakistan, the Indus is a 264-km-long multiple-channel reach. Remote sensing imagery, encompassing major floods in 2007 and 2010, shows that the Indus has a minimum of two and a maximum of nine channels, with on average four active channels during the dry season and five during the annual monsoon. Thus, the network structure, if not detailed planform, remains stable even for the record 2010 flood (27,100 m3 s- 1; recurrence interval > 100 years). Bankline recession is negligible for discharges less than a peak annual discharge of 6000 m3 s- 1 ( 80% of mean annual flood). The Maximum Flow Efficiency (MFE) principle demonstrates that the channel network is insensitive to the monsoon floods, which typically peak at 13,200 m3 s- 1. Rather, the network is in near-equilibrium with the mean annual flood (7530 m3 s- 1). The MFE principle indicates that stable networks have three to four channels, thus the observed stability in the number of active channels accords with the presence of a near-equilibrium reach-scale channel network. Insensitivity to the annual hydrological cycle demonstrates that the timescale for network adjustment is much longer than the timescale of the monsoon hydrograph, with the annual excess water being stored on floodplains rather than being conveyed in an enlarged channel network. The analysis explains the lack of significant channel adjustment following the largest flood in 40 years and the extensive Indus flooding experienced on an annual basis, with its substantial impacts on the populace and agricultural production.

Information transmission and signal permutation in active flow networks

NASA Astrophysics Data System (ADS)

Woodhouse, Francis G.; Fawcett, Joanna B.; Dunkel, Jörn

2018-03-01

Recent experiments show that both natural and artificial microswimmers in narrow channel-like geometries will self-organise to form steady, directed flows. This suggests that networks of flowing active matter could function as novel autonomous microfluidic devices. However, little is known about how information propagates through these far-from-equilibrium systems. Through a mathematical analogy with spin-ice vertex models, we investigate here the input–output characteristics of generic incompressible active flow networks (AFNs). Our analysis shows that information transport through an AFN is inherently different from conventional pressure or voltage driven networks. Active flows on hexagonal arrays preserve input information over longer distances than their passive counterparts and are highly sensitive to bulk topological defects, whose presence can be inferred from marginal input–output distributions alone. This sensitivity further allows controlled permutations on parallel inputs, revealing an unexpected link between active matter and group theory that can guide new microfluidic mixing strategies facilitated by active matter and aid the design of generic autonomous information transport networks.

Adaptiveness in monotone pseudo-Boolean optimization and stochastic neural computation.

PubMed

Grossi, Giuliano

2009-08-01

Hopfield neural network (HNN) is a nonlinear computational model successfully applied in finding near-optimal solutions of several difficult combinatorial problems. In many cases, the network energy function is obtained through a learning procedure so that its minima are states falling into a proper subspace (feasible region) of the search space. However, because of the network nonlinearity, a number of undesirable local energy minima emerge from the learning procedure, significantly effecting the network performance. In the neural model analyzed here, we combine both a penalty and a stochastic process in order to enhance the performance of a binary HNN. The penalty strategy allows us to gradually lead the search towards states representing feasible solutions, so avoiding oscillatory behaviors or asymptotically instable convergence. Presence of stochastic dynamics potentially prevents the network to fall into shallow local minima of the energy function, i.e., quite far from global optimum. Hence, for a given fixed network topology, the desired final distribution on the states can be reached by carefully modulating such process. The model uses pseudo-Boolean functions both to express problem constraints and cost function; a combination of these two functions is then interpreted as energy of the neural network. A wide variety of NP-hard problems fall in the class of problems that can be solved by the model at hand, particularly those having a monotonic quadratic pseudo-Boolean function as constraint function. That is, functions easily derived by closed algebraic expressions representing the constraint structure and easy (polynomial time) to maximize. We show the asymptotic convergence properties of this model characterizing its state space distribution at thermal equilibrium in terms of Markov chain and give evidence of its ability to find high quality solutions on benchmarks and randomly generated instances of two specific problems taken from the computational graph theory.

Simulating Chemistry in Star Forming Environments

NASA Astrophysics Data System (ADS)

Gong, Munan

Chemistry plays an important role in the interstellar medium (ISM), regulating the heating and cooling of the gas and determining abundances of molecular species that trace gas properties in observations. One of the most abundant and important molecules in the ISM is CO. CO is a main coolant for the molecular ISM, and the CO(J = 1 - 0) line emission is a widely used observational tracer for molecular clouds. In Chapter 2, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photodissociation region (PDR) code, and also with the Nelson & Langer (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. We also compare with observations of diffuse and translucent clouds. We find that the CO, CHx and OHx abundances are consistent with equilibrium predictions for densities n = 100 - 1000 cm-3, but the predicted equilibrium CI abundance is higher than observations, signaling the potential importance of non-equilibrium/dynamical effects. In Chapter 3, we apply our new chemistry network to a study of the XCO conversion factor, which is used to convert the CO luminosity to the total H2 mass. We use numerical simulations to investigate how XCO depends on numerical resolution, non-equilibrium chemistry, physical environment, and observational beam size. Our study employs 3D magnetohydrodynamics (MHD) simulations of galactic disks with solar neighborhood conditions, where star formation and the three-phase interstellar medium (ISM) is self-consistently generated by the interaction between gravity and stellar feedback. Synthetic CO maps are obtained by post-processing the MHD simulations with chemistry and radiation transfer. We find that CO is only an approximate tracer of H2. Nevertheless, 〈 XCO〉 = 0.7 - 1.0 x 1020 cm-2K-1km-1 s consistent with observations, insensitive to the evolutionary ISM state or the far-ultraviolet (FUV) radiation field strength. Our numerical simulations successfully reproduce the observed variations of X CO on parsec scales, as well as the dependence of X CO on extinction and the CO excitation temperature.

A stochastic equilibrium model for the North American natural gas market

NASA Astrophysics Data System (ADS)

Zhuang, Jifang

This dissertation is an endeavor in the field of energy modeling for the North American natural gas market using a mixed complementarity formulation combined with the stochastic programming. The genesis of the stochastic equilibrium model presented in this dissertation is the deterministic market equilibrium model developed in [Gabriel, Kiet and Zhuang, 2005]. Based on some improvements that we made to this model, including proving new existence and uniqueness results, we present a multistage stochastic equilibrium model with uncertain demand for the deregulated North American natural gas market using the recourse method of the stochastic programming. The market participants considered by the model are pipeline operators, producers, storage operators, peak gas operators, marketers and consumers. Pipeline operators are described with regulated tariffs but also involve "congestion pricing" as a mechanism to allocate scarce pipeline capacity. Marketers are modeled as Nash-Cournot players in sales to the residential and commercial sectors but price-takers in all other aspects. Consumers are represented by demand functions in the marketers' problem. Producers, storage operators and peak gas operators are price-takers consistent with perfect competition. Also, two types of the natural gas markets are included: the long-term and spot markets. Market participants make both high-level planning decisions (first-stage decisions) in the long-term market and daily operational decisions (recourse decisions) in the spot market subject to their engineering, resource and political constraints, resource constraints as well as market constraints on both the demand and the supply side, so as to simultaneously maximize their expected profits given others' decisions. The model is shown to be an instance of a mixed complementarity problem (MiCP) under minor conditions. The MiCP formulation is derived from applying the Karush-Kuhn-Tucker optimality conditions of the optimization problems faced by the market participants. Some theoretical results regarding the market prices in both markets are shown. We also illustrate the model on a representative, sample network of two production nodes, two consumption nodes with discretely distributed end-user demand and three seasons using four cases.

Convergence of neural networks for programming problems via a nonsmooth Lojasiewicz inequality.

PubMed

Forti, Mauro; Nistri, Paolo; Quincampoix, Marc

2006-11-01

This paper considers a class of neural networks (NNs) for solving linear programming (LP) problems, convex quadratic programming (QP) problems, and nonconvex QP problems where an indefinite quadratic objective function is subject to a set of affine constraints. The NNs are characterized by constraint neurons modeled by ideal diodes with vertical segments in their characteristic, which enable to implement an exact penalty method. A new method is exploited to address convergence of trajectories, which is based on a nonsmooth Lojasiewicz inequality for the generalized gradient vector field describing the NN dynamics. The method permits to prove that each forward trajectory of the NN has finite length, and as a consequence it converges toward a singleton. Furthermore, by means of a quantitative evaluation of the Lojasiewicz exponent at the equilibrium points, the following results on convergence rate of trajectories are established: (1) for nonconvex QP problems, each trajectory is either exponentially convergent, or convergent in finite time, toward a singleton belonging to the set of constrained critical points; (2) for convex QP problems, the same result as in (1) holds; moreover, the singleton belongs to the set of global minimizers; and (3) for LP problems, each trajectory converges in finite time to a singleton belonging to the set of global minimizers. These results, which improve previous results obtained via the Lyapunov approach, are true independently of the nature of the set of equilibrium points, and in particular they hold even when the NN possesses infinitely many nonisolated equilibrium points.

Dynamics of internal models in game players

NASA Astrophysics Data System (ADS)

Taiji, Makoto; Ikegami, Takashi

1999-10-01

A new approach for the study of social games and communications is proposed. Games are simulated between cognitive players who build the opponent’s internal model and decide their next strategy from predictions based on the model. In this paper, internal models are constructed by the recurrent neural network (RNN), and the iterated prisoner’s dilemma game is performed. The RNN allows us to express the internal model in a geometrical shape. The complicated transients of actions are observed before the stable mutually defecting equilibrium is reached. During the transients, the model shape also becomes complicated and often experiences chaotic changes. These new chaotic dynamics of internal models reflect the dynamical and high-dimensional rugged landscape of the internal model space.

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE PAGES

Tang, Jin-Yun; Riley, William J.

2017-09-05

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang, Jin-Yun; Riley, William J.

Several land biogeochemical models used for studying carbon–climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate–consumer relationships for consumer-mediated redox reactions of the form A +  B E→  products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate–consumer relationships can be formulated as the computationally difficult full equilibrium chemistrymore » problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation – ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate–consumer interactions and can be applied to improve Earth system model (ESM) land models.« less

SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration

NASA Astrophysics Data System (ADS)

Tang, Jin-Yun; Riley, William J.

2017-09-01

Several land biogeochemical models used for studying carbon-climate feedbacks have begun explicitly representing microbial dynamics. However, to our knowledge, there has been no theoretical work on how to achieve a consistent scaling of the complex biogeochemical reactions from microbial individuals to populations, communities, and interactions with plants and mineral soils. We focus here on developing a mathematical formulation of the substrate-consumer relationships for consumer-mediated redox reactions of the form A + BE→ products, where products could be, e.g., microbial biomass or bioproducts. Under the quasi-steady-state approximation, these substrate-consumer relationships can be formulated as the computationally difficult full equilibrium chemistry problem or approximated analytically with the dual Monod (DM) or synthesizing unit (SU) kinetics. We find that DM kinetics is scaling inconsistently for reaction networks because (1) substrate limitations are not considered, (2) contradictory assumptions are made regarding the substrate processing rate when transitioning from single- to multi-substrate redox reactions, and (3) the product generation rate cannot be scaled from one to multiple substrates. In contrast, SU kinetics consistently scales the product generation rate from one to multiple substrates but predicts unrealistic results as consumer abundances reach large values with respect to their substrates. We attribute this deficit to SU's failure to incorporate substrate limitation in its derivation. To address these issues, we propose SUPECA (SU plus the equilibrium chemistry approximation - ECA) kinetics, which consistently imposes substrate and consumer mass balance constraints. We show that SUPECA kinetics satisfies the partition principle, i.e., scaling invariance across a network of an arbitrary number of reactions (e.g., as in Newton's law of motion and Dalton's law of partial pressures). We tested SUPECA kinetics with the equilibrium chemistry solution for some simple problems and found SUPECA outperformed SU kinetics. As an example application, we show that a steady-state SUPECA-based approach predicted an aerobic soil respiration moisture response function that agreed well with laboratory observations. We conclude that, as an extension to SU and ECA kinetics, SUPECA provides a robust mathematical representation of complex soil substrate-consumer interactions and can be applied to improve Earth system model (ESM) land models.