Showing papers in "arXiv: Statistical Mechanics in 2001"

The Web of Human Sexual Contacts

Abstract: Many ``real-world'' networks are clearly defined while most ``social'' networks are to some extent subjective. Indeed, the accuracy of empirically-determined social networks is a question of some concern because individuals may have distinct perceptions of what constitutes a social link. One unambiguous type of connection is sexual contact. Here we analyze data on the sexual behavior of a random sample of individuals, and find that the cumulative distributions of the number of sexual partners during the twelve months prior to the survey decays as a power law with similar exponents $\alpha \approx 2.4$ for females and males. The scale-free nature of the web of human sexual contacts suggests that strategic interventions aimed at preventing the spread of sexually-transmitted diseases may be the most efficient approach.

On the coherence of Expected Shortfall

Abstract: Expected Shortfall (ES) in several variants has been proposed as remedy for the defi-ciencies of Value-at-Risk (VaR) which in general is not a coherent risk measure. In fact, most definitions of ES lead to the same results when applied to continuous loss distributions. Differences may appear when the underlying loss distributions have discontinuities. In this case even the coherence property of ES can get lost unless one took care of the details in its definition. We compare some of the definitions of Expected Shortfall, pointing out that there is one which is robust in the sense of yielding a coherent risk measure regardless of the underlying distributions. Moreover, this Expected Shortfall can be estimated effectively even in cases where the usual estimators for VaR fail. Key words: Expected Shortfall; Risk measure; worst conditional expectation; tail con-ditional expectation; value-at-risk (VaR); conditional value-at-risk (CVaR); tail mean; co-herence; quantile; sub-additivity.

Fractional Calculus: Some Basic Problems in Continuum and Statistical Mechanics

Abstract: We review some applications of fractional calculus developed by the author (partly in collaboration with others) to treat some basic problems in continuum and statistical mechanics. The problems in continuum mechanics concern mathematical modelling of viscoelastic bodies (Sect. 1), and unsteady motion of a particle in a viscous fluid, i.e. the Basset problem (Sect. 2). In the former analysis fractional calculus leads us to introduce intermediate models of viscoelasticity which generalize the classical spring-dashpot models. The latter analysis induces us to introduce a hydrodynamic model suitable to revisit in Sect. 3 the classical theory of the Brownian motion, which is a relevant topic in statistical mechanics. By the tools of fractional calculus we explain the long tails in the velocity correlation and in the displacement variance. In Sect. 4 we consider the fractional diffusion-wave equation, which is obtained from the classical diffusion equation by replacing the first-order time derivative by a fractional derivative of order $0< \beta <2$. Led by our analysis we express the fundamental solutions (the Green functions) in terms of two interrelated auxiliary functions in the similarity variable, which turn out to be of Wright type (see Appendix), and to distinguish slow-diffusion processes ($0 < \beta < 1$) from intermediate processes ($1 < \beta < 2$).

Simulation of pedestrian dynamics using a 2-dimensional cellular automaton

Abstract: We propose a 2-dimensional cellular automaton model to simulate pedestrian traffic. It is a vmax=1 model with exclusion statistics and parallel dynamics. Long-range interactions between the pedestrians are mediated by a so called floor field which modifies the transition rates to neighbouring cells. This field, which can be discrete or continuous, is subject to diffusion and decay. Furthermore it can be modified by the motion of the pedestrians. Therefore the model uses an idea similar to chemotaxis, but with pedestrians following a virtual rather than a chemical trace. Our main goal is to show that the introduction of such a floor field is sufficient to model collective effects and self-organization encountered in pedestrian dynamics, e.g. lane formation in counterflow through a large corridor. As an application we also present simulations of the evacuation of a large room with reduced visibility, e.g. due to failure of lights or smoke.

Measuring preferential attachment for evolving networks

Abstract: A key ingredient of current models proposed to capture the topological evolution of complex networks is the hypothesis that highly connected nodes increase their connectivity faster than their less connected peers, a phenomenon called preferential attachment. Measurements on four networks, namely the science citation network, Internet, actor collaboration and science coauthorship network indicate that the rate at which nodes acquire links depends on the node's degree, offering direct quantitative support for the presence of preferential attachment. We find that for the first two systems the attachment rate depends linearly on the node degree, while for the latter two the dependence follows a sublinear power law.

Macroscopic fluctuation theory for stationary non equilibrium states

Abstract: We formulate a dynamical fluctuation theory for stationary non equilibrium states (SNS) which is tested explicitly in stochastic models of interacting particles. In our theory a crucial role is played by the time reversed dynamics. Within this theory we derive the following results: the modification of the Onsager-Machlup theory in the SNS; a general Hamilton-Jacobi equation for the macroscopic entropy; a non equilibrium, non linear fluctuation dissipation relation valid for a wide class of systems; an H theorem for the entropy. We discuss in detail two models of stochastic boundary driven lattice gases: the zero range and the simple exclusion processes. In the first model the invariant measure is explicitly known and we verify the predictions of the general theory. For the one dimensional simple exclusion process, as recently shown by Derrida, Lebowitz, and Speer, it is possible to express the macroscopic entropy in terms of the solution of a non linear ordinary differential equation; by using the Hamilton-Jacobi equation, we obtain an logically independent derivation of this result.

Vortices in a trapped dilute Bose-Einstein condensate

Abstract: We review the theory of vortices in trapped dilute Bose-Einstein condensates and compare theoretical predictions with existing experiments. Mean-field theory based on the time-dependent Gross-Pitaevskii equation describes the main features of the vortex states, and its predictions agree well with available experimental results. We discuss various properties of a single vortex, including its structure, energy, dynamics, normal modes and stability, as well as vortex arrays. When the nonuniform condensate contains a vortex, the excitation spectrum includes unstable (``anomalous'') mode(s) with negative frequency. Trap rotation shifts the normal-mode frequencies and can stabilize the vortex. We consider the effect of thermal quasiparticles on vortex normal modes as well as possible mechanisms for vortex dissipation. Vortex states in mixtures and spinor condensates are also discussed.

Cellular Automaton Approach to Pedestrian Dynamics - Theory

Abstract: We present a 2-dimensional cellular automaton model for the simulation of pedestrian dynamics. The model is extremely efficient and allows simulations of large crowds faster than real time since it includes only nearest-neighbour interactions. Nevertheless it is able to reproduce collective effects and self-organization encountered in pedestrian dynamics. This is achieved by introducing a so-called floor field which mediates the long-range interactions between the pedestrians. This field modifies the transition rates to neighbouring cells. It has its own dynamics (diffusion and decay) and can be changed by the motion of the pedestrians. Therefore the model uses an idea similar to chemotaxis, but with pedestrians following a virtual rather than a chemical trace.

Ego-centered networks and the ripple effect

Abstract: Recent work has demonstrated that many social networks, and indeed many networks of other types also, have broad distributions of vertex degree. Here we show that this has a substantial impact on the shape of ego-centered networks, i.e., sets of network vertices that are within a given distance of a specified central vertex, the ego. This in turn affects concepts and methods based on ego-centered networks, such as snowball sampling and the "ripple effect". In particular, we argue that one's acquaintances, one's immediate neighbors in the acquaintance network, are far from being a random sample of the population, and that this biases the numbers of neighbors two and more steps away. We demonstrate this concept using data drawn from academic collaboration networks, for which, as we show, current simple theories for the typical size of ego-centered networks give numbers that differ greatly from those measured in reality. We present an improved theoretical model which gives significantly better results.

Language as an Evolving Word Web

Abstract: Human language can be described as a complex network of linked words. In such a treatment, each distinct word in language is a vertex of this web, and neighboring words in sentences are connected by edges. It was recently found (Ferrer and Sol\'e) that the distribution of the numbers of connections of words in such a network is of a peculiar form which includes two pronounced power-law regions. Here we treat language as a self-organizing network of interacting words. In the framework of this concept, we completely describe the observed Word Web structure without fitting.

Expected Shortfall as a Tool for Financial Risk Management

Abstract: We study the properties of Expected Shortfall from the point of view of financial risk management. This measure --- which emerges as a natural remedy in some cases where Value at Risk (VaR) is not able to distinguish portfolios which bear different levels of risk --- is indeed shown to have much better properties than VaR. We show in fact that unlike VaR this variable is in general subadditive and therefore it is a Coherent Measure of Risk in the sense of reference (artzner)

Sub-critical and Super-critical Regimes in Epidemic Models of Earthquake Aftershocks

Abstract: We present an analytical solution and numerical tests of the epidemic-type aftershock (ETAS) model for aftershocks, which describes foreshocks, aftershocks and mainshocks on the same footing. The occurrence rate of aftershocks triggered by a single mainshock decreases with the time from the mainshock according to the modified Omori law K/(t+c)^p with p=1+theta. A mainshock at time t=0 triggers aftershocks according to the local Omori law, that in turn trigger their own aftershocks and so on. The effective branching parameter n, defined as the mean aftershock number triggered per event, controls the transition between a sub-critical regime n 1. In the sub-critical regime, we recover and document the crossover from an Omori exponent 1-theta for t 1 and theta>0, we find a novel transition from an Omori decay law with exponent 1-theta fot t t*. The case theta<0 yields an infinite n-value. In this case, we find another characteristic time tau controlling the crossover from an Omori law with exponent 1-theta for t

Boltzmann's Approach to Statistical Mechanics

Abstract: In the last quarter of the nineteenth century, Ludwig Boltzmann explained how irreversible macroscopic laws, in particular the second law of thermodynamics, originate in the time-reversible laws of microscopic physics. Boltzmann's analysis, the essence of which I shall review here, is basically correct. The most famous criticisms of Boltzmann's later work on the subject have little merit. Most twentieth century innovations -- such as the identification of the state of a physical system with a probability distribution $\varrho$ on its phase space, of its thermodynamic entropy with the Gibbs entropy of $\varrho$, and the invocation of the notions of ergodicity and mixing for the justification of the foundations of statistical mechanics -- are thoroughly misguided.

Expectation bubbles in a spin model of markets: Intermittency from frustration across scales

Abstract: A simple spin model is studied, motivated by the dynamics of traders in a market where expectation bubbles and crashes occur. The dynamics is governed by interactions which are frustrated across different scales: While ferromagnetic couplings connect each spin to its local neighborhood, an additional coupling relates each spin to the global magnetization. This new coupling is allowed to be anti-ferromagnetic. The resulting frustration causes a metastable dynamics with intermittency and phases of chaotic dynamics. The model reproduces main observations of real economic markets as power-law distributed returns and clustered volatility.

Significance of log-periodic precursors to financial crashes

Abstract: We clarify the status of log-periodicity associated with speculative bubbles preceding financial crashes. In particular, we address Feigenbaum's [2001] criticism and show how it can be rebuked. Feigenbaum's main result is as follows: ``the hypothesis that the log-periodic component is present in the data cannot be rejected at the 95% confidence level when using all the data prior to the 1987 crash; however, it can be rejected by removing the last year of data.'' (e.g., by removing 15% of the data closest to the critical point). We stress that it is naive to analyze a critical point phenomenon, i.e., a power law divergence, reliably by removing the most important part of the data closest to the critical point. We also present the history of log-periodicity in the present context explaining its essential features and why it may be important. We offer an extension of the rational expectation bubble model for general and arbitrary risk-aversion within the general stochastic discount factor theory. We suggest guidelines for using log-periodicity and explain how to develop and interpret statistical tests of log-periodicity. We discuss the issue of prediction based on our results and the evidence of outliers in the distribution of drawdowns. New statistical tests demonstrate that the 1% to 10% quantile of the largest events of the population of drawdowns of the Nasdaq composite index and of the Dow Jones Industrial Average index belong to a distribution significantly different from the rest of the population. This suggests that very large drawdowns result from an amplification mechanism that may make them more predictable than smaller market moves.

A Stochastic Model of Human Gait Dynamics

Abstract: We present a stochastic model of gait rhythm dynamics, based on transitions between different ``neural centers'', that reproduces distinctive statistical properties of normal human walking. By tuning one model parameter, the hopping range, the model can describe alterations in gait dynamics from childhood to adulthood --- including a decrease in the correlation and volatility exponents with maturation. The model also generates time series with multifractal spectra whose broadness depends only on this parameter.

The thermal conductivity of the spin-1/2 XXZ chain at arbitrary temperature

Abstract: Motivated by recent investigations of transport properties of strongly correlated 1d models and thermal conductivity measurements of quasi 1d magnetic systems we present results for the integrable spin-1/2 $XXZ$ chain. The thermal conductivity $\kappa(\omega)$ of this model has $\Re\kappa(\omega)=\tilde\kappa \delta(\omega)$, i.e. it is infinite for zero frequency $\omega$. The weight $\tilde\kappa$ of the delta peak is calculated exactly by a lattice path integral formulation. Numerical results for wide ranges of temperature and anisotropy are presented. The low and high temperature limits are studied analytically.

Dynamics of Sleep-Wake Transitions During Sleep

Abstract: We study the dynamics of the awakening during the night for healthy subjects and find that the wake and the sleep periods exhibit completely different behavior: the durations of wake periods are characterized by a scale-free power-law distribution, while the durations of sleep periods have an exponential distribution with a characteristic time scale. We find that the characteristic time scale of sleep periods changes throughout the night. In contrast, there is no measurable variation in the power-law behavior for the durations of wake periods. We develop a stochastic model which agrees with the data and suggests that the difference in the dynamics of sleep and wake states arises from the constraints on the number of microstates in the sleep-wake system.

H-Theorem and Generalized Entropies Within the Framework of Non Linear Kinetics

Abstract: In the present effort we consider the most general non linear particle kinetics within the framework of the Fokker-Planck picture. We show that the kinetics imposes the form of the generalized entropy and subsequently we demonstrate the H-theorem. The particle statistical distribution is obtained, both as stationary solution of the non linear evolution equation and as the state which maximizes the generalized entropy. The present approach allows to treat the statistical distributions already known in the literature in a unifying scheme. As a working example we consider the kinetics, constructed by using the $\kappa$-exponential $\exp_{_{\{\kappa\}}}(x)= (\sqrt{1+\kappa^2x^2}+\kappa x)^{1/\kappa}$ recently proposed which reduces to the standard exponential as the deformation parameter $\kappa$ approaches to zero and presents the relevant power law asymptotic behaviour $\exp_{_{\{\kappa\}}}(x){\atop\stackrel\sim x\to \pm \infty}|2\kappa x|^{\pm 1/|\kappa|}$. The $\kappa$-kinetics obeys the H-theorem and in the case of Brownian particles, admits as stationary state the distribution $f=Z^{-1}\exp_{_{\{\kappa\}}}[-(\beta mv^2/2-\mu)]$ which can be obtained also by maximizing the entropy $S_{\kappa}=\int d^n v [ c(\kappa)f^{1+\kappa}+c(-\kappa)f^{1-\kappa}]$ with $c(\kappa)=-Z^{\kappa}/ [2\kappa(1+\kappa)]$ after properly constrained.

Power law relaxation in a complex system: Omori law after a financial market crash

Abstract: We study the relaxation dynamics of a financial market just after the occurrence of a crash by investigating the number of times the absolute value of an index return is exceeding a given threshold value. We show that the empirical observation of a power law evolution of the number of events exceeding the selected threshold (a behavior known as the Omori law in geophysics) is consistent with the simultaneous occurrence of (i) a return probability density function characterized by a power law asymptotic behavior and (ii) a power law relaxation decay of its typical scale. Our empirical observation cannot be explained within the framework of simple and widespread stochastic volatility models.

Tangled Nature: a model of evolutionary ecology

Abstract: We discuss a simple model of co-evolution. In order to emphasise the effect of interaction between individuals the entire population is subjected to the same physical environment. Species are emergent structures and extinction, origination and diversity are entirely a consequence of co-evolutionary interaction between individuals. For comparison we consider both asexual and sexually reproducing populations. We also study competition between asexual and sexual reproduction in a mixed population. In either case the system evolves through periods of hectic reorganisation separated by periods of coherent stable coexistence.

Metabolically Efficient Information Processing

Abstract: Energy efficient information transmission may be relevant to biological sensory signal processing as well as to low power electronic devices. We explore its consequences in two different regimes. In an ``immediate'' regime, we argue that the information rate should be maximized subject to a power constraint, while in an ``exploratory'' regime, the transmission rate per power cost should be maximized. In the absence of noise, discrete inputs are optimally encoded into Boltzmann distributed output symbols. In the exploratory regime, the partition function of this distribution is numerically equal to 1. The structure of the optimal code is strongly affected by noise in the transmission channel. The Arimoto-Blahut algorithm, generalized for cost constraints, can be used to derive and interpret the distribution of symbols for optimal energy efficient coding in the presence of noise. We outline the possibilities and problems in extending our results to information coding and transmission in neurobiological systems.

Multiplicity of Generation, Selection, and Classification Procedures for Jammed Hard-Particle Packings

Abstract: Hard-particle packings have served as useful starting points to study the structure of diverse systems such as liquids, living cells, granular media, glasses, and amorphous solids. Howard Reiss has played a major role in helping to illuminate our understanding of hard-particle systems, which still offer scientists many interesting conundrums. Jammed configurations of hard particles are of great fundamental and practical interest. What one precisely means by a "jammed" configuration is quite subtle and considerable ambiguity remains in the literature on this question. We will show that there is a multiplicity of generation, selection, and classification procedures for jammed configurations of identical d-dimensional spheres. We categorize common ordered lattices according to our definitions and discuss implications for random disk and sphere packings. We also show how the concept of rigidity percolation (which has been used to understand the mechanical properties of network glasses) can be generalized to further characterize hard-sphere packings.

Cellular Automaton Approach to Pedestrian Dynamics - Applications

Abstract: We present applications of a cellular automaton approach to pedestrian dynamics introduced in [1,2]. It is shown that the model is able to reproduce collective effects and self-organization phenomena encountered in pedestrian traffic, e.g. lane formation in counterflow through a large corridor and oscillations at doors. Furthermore we present simple examples where the model is applied to the simulation of evacuation processes.

Completeness of the Bethe ansatz for the six and eight-vertex models

Abstract: We discuss some of the difficulties that have been mentioned in the literature in connection with the Bethe ansatz for the six-vertex model and XXZ chain, and for the eight-vertex model. In particular we discuss the ``beyond the equator'', infinite momenta and exact complete string problems. We show how they can be overcome and conclude that the coordinate Bethe ansatz does indeed give a complete set of states, as expected.

Analyzing and modelling 1+1d markets

Abstract: We report a statistical analysis of the Island ECN (NASDAQ) order book. We determine the static and dynamic properties of this system, and then analyze them from a physicist's viewpoint using an equivalent particle system obtained by treating orders as massive particles and price as position. We identify the fundamental dynamical processes, test existing particles models of such markets against our findings, and introduce a new model of limit order markets.

Efficiency of DNA replication in the Polymerase Chain Reactio

Abstract: We present a detailed kinetic model for the Polymerase Chain Reaction, and model the probability of replication in terms of the physical parameters of the problem. Applying the theory of branching processes, we show the existance of a new phenomenon affecting the probability distribution function of the number of replicants: for small number of initial molecules, the limiting behavior (with increasing number of cycles) of the pdf is represented by a multi-modal funtion.

Comparable system-level organization of Archaea and Eukaryotes

Abstract: A central and long-standing issue in evolutionary theory is the origin of the biological variation upon which natural selection acts1. Some hypotheses suggest that evolutionary change represents an adaptation to the surrounding environment within the constraints of an organism's innate characteristics. Elucidation of the origin and evolutionary relationship of species has been complemented by nucleotide sequence and gene content analyses, with profound implications for recognizing life's major domains. Understanding of evolutionary relationships may be further expanded by comparing systemic higher-level organization among species. Here we employ multivariate analyses to evaluate the biochemical reaction pathways characterizing 43 species. Comparison of the information transfer pathways of Archaea and Eukaryotes indicates a close relationship between these domains. In addition, whereas eukaryotic metabolic enzymes are primarily of bacterial origin, the pathway-level organization of archaeal and eukaryotic metabolic networks is more closely related. Our analyses therefore suggest that during the symbiotic evolution of eukaryotes, incorporation of bacterial metabolic enzymes into the proto-archaeal proteome was constrained by the host's pre-existing metabolic architecture.

Opinion evolution in closed community

Abstract: A simple Ising spin model which can describe a mechanism of making a decision in a closed community is proposed. It is shown via standard Monte Carlo simulations that very simple rules lead to rather complicated dynamics and to a power law in the decision time distribution. It is found that a closed community has to evolve either to a dictatorship or a stalemate state (inability to take any common decision). A common decision can be taken in a "democratic way" only by an open community.

Computer simulation of fluid phase transitions

Abstract: The task of accurately locating fluid phase boundaries by means of computer simulation is hampered by problems associated with sampling both coexisting phases in a single simulation run. We explain the physical background to these problems and describe how they can be tackled using a synthesis of biased Monte Carlo sampling and histogram extrapolation methods, married to a standard fluid simulation algorithm. It is demonstrated that the combined approach provides a powerful method for tracing fluid phase boundaries.