Catalog footprint

What is connected

81works
32topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

81 published item(s)

preprint2022arXiv

Role of collective information in networks of quantum operating agents

A network of agents is considered whose decision processes are described by the quantum decision theory previously advanced by the authors. Decision making is done by evaluating the utility of alternatives, their attractiveness, and the available information, whose combinations form the probabilities to choose a given alternative. As a result of the interplay between these three contributions, the process of choice between several alternatives is multimodal. The agents interact by exchanging information, which can take two forms: (i) information that an agent can directly receive from another agent and (ii) information collectively created by the members of the society. The information field common to all agents tends to smooth out sharp variations in the temporal behaviour of the probabilities and can even remove them. For agents with short-term memory, the probabilities often tend to their limiting values through strong oscillations and, for a range of parameters, these oscillations last for ever, representing an ever lasting hesitation of the decision makers. Switching on the information field makes the amplitude of the oscillations smaller and even can halt the everlasting oscillations forcing the probabilities to converge to fixed limits. The dynamic disjunction effect is described.

preprint2016arXiv

Inconclusive quantum measurements and decisions under uncertainty

We give a mathematical definition for the notion of inconclusive quantum measurements. In physics, such measurements occur at intermediate stages of a complex measurement procedure, with the final measurement result being operationally testable. Since the mathematical structure of Quantum Decision Theory has been developed in analogy with the theory of quantum measurements, the inconclusive quantum measurements correspond, in Quantum Decision Theory, to intermediate stages of decision making in the process of taking decisions under uncertainty. The general form of the quantum probability for a composite event is the sum of a utility factor, describing a rational evaluation of the considered prospect, and of an attraction factor, characterizing irrational, subconscious attitudes of the decision maker. Despite the involved irrationality, the probability of prospects can be evaluated. This is equivalent to the possibility of calculating quantum probabilities without specifying hidden variables. We formulate a general way of evaluation, based on the use of non-informative priors. As an example, we suggest the explanation of the decoy effect. Our quantitative predictions are in very good agreement with experimental data.

preprint2016arXiv

Quantum decision theory in simple risky choices

Quantum decision theory (QDT) is a recently developed theory of decision making based on the mathematics of Hilbert spaces, a framework known in physics for its application to quantum mechanics. This framework formalizes the concept of uncertainty and other effects that are particularly manifest in cognitive processes, which makes it well suited for the study of decision making. QDT describes a decision maker's choice as a stochastic event occurring with a probability that is the sum of an objective utility factor and a subjective attraction factor. QDT offers a prediction for the average effect of subjectivity on decision makers, the quarter law. We examine individual and aggregated (group) data, and find that the results are in good agreement with the quarter law at the level of groups. At the individual level, it appears that the quarter law could be refined in order to reflect individual characteristics. This article revisits the formalism of QDT along a concrete example and offers a practical guide to researchers who are interested in applying QDT to a dataset of binary lotteries in the domain of gains.

preprint2016arXiv

Quantum probability and quantum decision making

A rigorous general definition of quantum probability is given, which is valid for elementary events and for composite events, for operationally testable measurements as well as for inconclusive measurements, and also for non-commuting observables in addition to commutative observables. Our proposed definition of quantum probability makes it possible to describe quantum measurements and quantum decision making on the same common mathematical footing. Conditions are formulated for the case when quantum decision theory reduces to its classical counterpart and for the situation where the use of quantum decision theory is necessary.

preprint2015arXiv

A civil super-Manhattan project in nuclear research for a safer and prosperous world

Humankind is confronted with a "nuclear stewardship curse", facing the prospect of needing to manage nuclear products over long time scales in the face of the short-time scales of human polities. I propose a super Manhattan-type effort to rejuvenate the nuclear energy industry to overcome the current dead-end in which it finds itself, and by force, humankind has trapped itself in. A 1% GDP investment over a decade in the main nuclear countries could boost economic growth with a focus on the real world, epitomised by nuclear physics/chemistry/engineering/economics with well defined targets. By investing vigorously to obtain scientific and technological breakthroughs, we can create the spring of a world economic rebound based on new ways of exploiting nuclear energy, both more safely and more durably.

preprint2015arXiv

Condensation of earthquake location distributions: Optimal spatial information encoding and application to multifractal analysis of South Californian seismicity

We present the "condensation" method that exploits the heterogeneity of the probability distribution functions (PDF) of event locations to improve the spatial information content of seismic catalogs. The method reduces the size of seismic catalogs while improving the access to the spatial information content of seismic catalogs. The PDFs of events are first ranked by decreasing location errors and then successively condensed onto better located and lower variance event PDFs. The obtained condensed catalog attributes different weights to each event, providing an optimal spatial representation with respect to the spatially varying location capability of the seismic network. Synthetic tests on fractal distributions perturbed with realistic location errors show that condensation improves spatial information content of the original catalog. Applied to Southern California seismicity, the new condensed catalog highlights major mapped fault traces and reveals possible additional structures while reducing the catalog length by ~25%. The condensation method allows us to account for location error information within a point based spatial analysis. We demonstrate this by comparing the multifractal properties of the condensed catalog locations with those of the original catalog. We evidence different spatial scaling regimes characterized by distinct multifractal spectra and separated by transition scales. We interpret the upper scale as to agree with the thickness of the brittle crust, while the lower scale (2.5km) might depend on the relocation procedure. Accounting for these new results, the Epidemic Type Aftershock Model formulation suggests that, contrary to previous studies, large earthquakes dominate the earthquake triggering process. This implies that the limited capability of detecting small magnitude events cannot be used to argue that earthquakes are unpredictable in general.

preprint2015arXiv

Dynamical system theory of periodically collapsing bubbles

We propose a reduced form set of two coupled continuous time equations linking the price of a representative asset and the price of a bond, the later quantifying the cost of borrowing. The feedbacks between asset prices and bonds are mediated by the dependence of their "fundamental values" on past asset prices and bond themselves. The obtained nonlinear self-referencing price dynamics can induce, in a completely objective deterministic way, the appearance of periodically exploding bubbles ending in crashes. Technically, the periodically explosive bubbles arise due to the proximity of two types of bifurcations as a function of the two key control parameters $b$ and $g$, which represent, respectively, the sensitivity of the fundamental asset price on past asset and bond prices and of the fundamental bond price on past asset prices. One is a Hopf bifurcation, when a stable focus transforms into an unstable focus and a limit cycle appears. The other is a rather unusual bifurcation, when a stable node and a saddle merge together and disappear, while an unstable focus survives and a limit cycle develops. The lines, where the periodic bubbles arise, are analogous to the critical lines of phase transitions in statistical physics. The amplitude of bubbles and waiting times between them respectively diverge with the critical exponents $γ= 1$ and $ν= 1/2$, as the critical lines are approached.

preprint2015arXiv

Positive operator-valued measures in quantum decision theory

We show that the correct mathematical foundation of quantum decision theory, dealing with uncertain events, requires the use of positive operator-valued measure that is a generalization of the projection-valued measure. The latter is appropriate for operationally testable events, while the former is necessary for characterizing operationally uncertain events. In decision making, one has to distinguish composite non-entangled events from composite entangled events. The mathematical definition of entangled prospects is based on the theory of Hilbert-Schmidt spaces and is analogous to the definition of entangled statistical operators in quantum information theory. We demonstrate that the necessary condition for the appearance of an interference term in the quantum probability is the occurrence of entangled prospects and the existence of an entangled strategic state of a decision maker. The origin of uncertainties in standard lotteries is explained.

preprint2015arXiv

Preference reversal in quantum decision theory

We consider the psychological effect of preference reversal and show that it finds a natural explanation in the frame of quantum decision theory. When people choose between lotteries with non-negative payoffs, they prefer a more certain lottery because of uncertainty aversion. But when people evaluate lottery prices, e.g. for selling to others the right to play them, they do this more rationally, being less subject to behavioral biases. This difference can be explained by the presence of the attraction factors entering the expression of quantum probabilities. Only the existence of attraction factors can explain why, considering two lotteries with close utility factors, a decision maker prefers one of them when choosing, but evaluates higher the other one when pricing. We derive a general quantitative criterion for the preference reversal to occur that relates the utilities of the two lotteries to the attraction factors under choosing versus pricing and test successfully its application on experiments by Tversky et al. We also show that the planning paradox can be treated as a kind of preference reversal.

preprint2015arXiv

Quantum decision making by social agents

The influence of additional information on the decision making of agents, who are interacting members of a society, is analyzed within the mathematical framework based on the use of quantum probabilities. The introduction of social interactions, which influence the decisions of individual agents, leads to a generalization of the quantum decision theory developed earlier by the authors for separate individuals. The generalized approach is free of the standard paradoxes of classical decision theory. This approach also explains the error-attenuation effects observed for the paradoxes occurring when decision makers, who are members of a society, consult with each other, increasing in this way the available mutual information. A precise correspondence between quantum decision theory and classical utility theory is formulated via the introduction of an intermediate probabilistic version of utility theory of a novel form, which obeys the requirement that zero-utility prospects should have zero probability weights.

preprint2015arXiv

Quantum theory of measurements as quantum decision theory

Theory of quantum measurements is often classified as decision theory. An event in decision theory corresponds to the measurement of an observable. This analogy looks clear for operationally testable simple events. However, the situation is essentially more complicated in the case of composite events. The most difficult point is the relation between decisions under uncertainty and measurements under uncertainty. We suggest a unified language for describing the processes of quantum decision making and quantum measurements. The notion of quantum measurements under uncertainty is introduced. We show that the correct mathematical foundation for the theory of measurements under uncertainty, as well as for quantum decision theory dealing with uncertain events, requires the use of positive operator-valued measure that is a generalization of projection-valued measure. The latter is appropriate for operationally testable events, while the former is necessary for characterizing operationally uncertain events. In both decision making and quantum measurements, one has to distinguish composite non-entangled events from composite entangled events. Quantum probability can be essentially different from classical probability only for entangled events. The necessary condition for the appearance of an interference term in the quantum probability is the occurrence of entangled prospects and the existence of an entangled strategic state of a decision maker or of an entangled statistical state of a measuring device.

preprint2015arXiv

Role of information in decision making of social agents

The influence of additional information on the decision making of agents, who are interacting members of a society, is analyzed within the mathematical framework based on the use of quantum probabilities. The introduction of social interactions, which influence the decisions of individual agents, leads to a generalization of the quantum decision theory developed earlier by the authors for separate individuals. The generalized approach is free of the standard paradoxes of classical decision theory. This approach also explains the error-attenuation effects observed for the paradoxes occurring when decision makers, who are members of a society, consult with each other, increasing in this way the available mutual information. A precise correspondence between quantum decision theory and classical utility theory is formulated via the introduction of an intermediate probabilistic version of utility theory of a novel form, which obeys the requirement that zero-utility prospects should have zero probability weights.

preprint2014arXiv

How brains make decisions

This chapter, dedicated to the memory of Mino Freund, summarizes the Quantum Decision Theory (QDT) that we have developed in a series of publications since 2008. We formulate a general mathematical scheme of how decisions are taken, using the point of view of psychological and cognitive sciences, without touching physiological aspects. The basic principles of how intelligence acts are discussed. The human brain processes involved in decisions are argued to be principally different from straightforward computer operations. The difference lies in the conscious-subconscious duality of the decision making process and the role of emotions that compete with utility optimization. The most general approach for characterizing the process of decision making, taking into account the conscious-subconscious duality, uses the framework of functional analysis in Hilbert spaces, similarly to that used in the quantum theory of measurements. This does not imply that the brain is a quantum system, but just allows for the simplest and most general extension of classical decision theory. The resulting theory of quantum decision making, based on the rules of quantum measurements, solves all paradoxes of classical decision making, allowing for quantitative predictions that are in excellent agreement with experiments. Finally, we provide a novel application by comparing the predictions of QDT with experiments on the prisoner dilemma game. The developed theory can serve as a guide for creating artificial intelligence acting by quantum rules.

preprint2014arXiv

Manipulating decision making of typical agents

We investigate how the choice of decision makers can be varied under the presence of risk and uncertainty. Our analysis is based on the approach we have previously applied to individual decision makers, which we now generalize to the case of decision makers that are members of a society. The approach employs the mathematical techniques that are common in quantum theory, justifying our naming as Quantum Decision Theory. However, we do not assume that decision makers are quantum objects. The techniques of quantum theory are needed only for defining the prospect probabilities taking into account such hidden variables as behavioral biases and other subconscious feelings. The approach describes an agent's choice as a probabilistic event occurring with a probability that is the sum of a utility factor and of an attraction factor. The attraction factor embodies subjective and unconscious dimensions in the mind of the decision maker. We show that the typical aggregate amplitude of the attraction factor is $1/4$, and it can be either positive or negative depending on the relative attraction of the competing choices. The most efficient way of varying the decision makers choice is realized by influencing the attraction factor. This can be done in two ways. One method is to arrange in a special manner the payoff weights, which induces the required changes of the values of attraction factors. We show that a slight variation of the payoff weights can invert the sign of the attraction factors and reverse the decision preferences, even when the prospect utilities remain unchanged. The second method of influencing the decision makers choice is by providing information to decision makers. The methods of influencing decision making are illustrated by several experiments, whose outcomes are compared quantitatively with the predictions of our approach.

preprint2014arXiv

New Approach to Modeling Symbiosis in Biological and Social Systems

We suggest a novel approach to treating symbiotic relations between biological species or social entities. The main idea is the characterisation of symbiotic relations of coexisting species through their mutual influence on their respective carrying capacities, taking into account that this influence can be quite strong and requires a nonlinear functional framework. We distinguish three variants of mutual influence, representing the main types of relations between species: (i) passive symbiosis, when the mutual carrying capacities are influenced by other species without their direct interactions; (ii) active symbiosis, when the carrying capacities are transformed by interacting species; and (iii) mixed symbiosis, when the carrying capacity of one species is influenced by direct interactions, while that of the other species is not. The approach allows us to describe all kinds of symbiosis, mutualism, commensalism, and parasitism within a unified scheme. The case of two symbiotic species is analysed in detail, demonstrating several dynamical regimes of coexistence, unbounded growth of both populations, growth of one and elimination of the other population, convergence to evolutionary stable states, and everlasting population oscillations. The change of the dynamical regimes occurs by varying the system parameters characterising the destruction or creation of the mutual carrying capacities. The regime changes are associated with several dynamical system bifurcations.

preprint2014arXiv

Physics and Financial Economics (1776-2014): Puzzles, Ising and Agent-Based models

This short review presents a selected history of the mutual fertilization between physics and economics, from Isaac Newton and Adam Smith to the present. The fundamentally different perspectives embraced in theories developed in financial economics compared with physics are dissected with the examples of the volatility smile and of the excess volatility puzzle. The role of the Ising model of phase transitions to model social and financial systems is reviewed, with the concepts of random utilities and the logit model as the analog of the Boltzmann factor in statistic physics. Recent extensions in term of quantum decision theory are also covered. A wealth of models are discussed briefly that build on the Ising model and generalize it to account for the many stylized facts of financial markets. A summary of the relevance of the Ising model and its extensions is provided to account for financial bubbles and crashes. The review would be incomplete if it would not cover the dynamical field of agent based models (ABMs), also known as computational economic models, of which the Ising-type models are just special ABM implementations. We formulate the ``Emerging Market Intelligence hypothesis'' to reconcile the pervasive presence of ``noise traders'' with the near efficiency of financial markets. Finally, we note that evolutionary biology, more than physics, is now playing a growing role to inspire models of financial markets.

preprint2014arXiv

Population dynamics with nonlinear delayed carrying capacity

We consider a class of evolution equations describing population dynamics in the presence of a carrying capacity depending on the population with delay. In an earlier work, we presented an exhaustive classification of the logistic equation where the carrying capacity is linearly dependent on the population with a time delay, which we refer to as the "linear delayed carrying capacity" model. Here, we generalize it to the case of a nonlinear delayed carrying capacity. The nonlinear functional form of the carrying capacity characterizes the delayed feedback of the evolving population on the capacity of their surrounding by either creating additional means for survival or destroying the available resources. The previously studied linear approximation for the capacity assumed weak feedback, while the nonlinear form is applicable to arbitrarily strong feedback. The nonlinearity essentially changes the behavior of solutions to the evolution equation, as compared to the linear case. All admissible dynamical regimes are analyzed, which can be of the following types: punctuated unbounded growth, punctuated increase or punctuated degradation to a stationary state, convergence to a stationary state with sharp reversals of plateaus, oscillatory attenuation, everlasting fluctuations, everlasting up-down plateau reversals, and divergence in finite time. The theorem is proved that, for the case characterizing the evolution under gain and competition, solutions are always bounded, if the feedback is destructive. We find that even a small noise level profoundly affects the position of the finite-time singularities. Finally, we demonstrate the feasibility of predicting the critical time of solutions having finite-time singularities from the knowledge of a simple quadratic approximation of the early time dynamics.

preprint2014arXiv

Quantum probabilities and entanglement for multimode quantum systems

Quantum probabilities are defined for several important physical cases characterizing measurements with multimode quantum systems. These are the probabilities for operationally testable measurements, for operationally uncertain measurements, and for entangled composite events. The role of the prospect and state entanglement is emphasized. Numerical modeling is presented for a two-mode Bose-condensed system of trapped atoms. The interference factor is calculated by invoking the channel-state duality.

preprint2014arXiv

Self-organization in complex systems as decision making

The idea is advanced that self-organization in complex systems can be treated as decision making (as it is performed by humans) and, vice versa, decision making is nothing but a kind of self-organization in the decision maker nervous systems. A mathematical formulation is suggested based on the definition of probabilities of system states, whose particular cases characterize the probabilities of structures, patterns, scenarios, or prospects. In this general framework, it is shown that the mathematical structures of self-organization and of decision making are identical. This makes it clear how self-organization can be seen as an endogenous decision making process and, reciprocally, decision making occurs via an endogenous self-organization. The approach is illustrated by phase transitions in large statistical systems, crossovers in small statistical systems, evolutions and revolutions in social and biological systems, structural self-organization in dynamical systems, and by the probabilistic formulation of classical and behavioral decision theories. In all these cases, self-organization is described as the process of evaluating the probabilities of macroscopic states or prospects in the search for a state with the largest probability. The general way of deriving the probability measure for classical systems is the principle of minimal information, that is, the conditional entropy maximization under given constraints. Behavioral biases of decision makers can be characterized in the same way as analogous to quantum fluctuations in natural systems.

preprint2014arXiv

Super-exponential endogenous bubbles in an equilibrium model of rational and noise traders

We introduce a model of super-exponential financial bubbles with two assets (risky and risk-free), in which rational investors and noise traders co-exist. Rational investors form expectations on the return and risk of a risky asset and maximize their constant relative risk aversion expected utility with respect to their allocation on the risky asset versus the risk-free asset. Noise traders are subjected to social imitation and follow momentum trading. Allowing for random time-varying herding propensity, we are able to reproduce several well-known stylized facts of financial markets such as a fat-tail distribution of returns and volatility clustering. In particular, we observe transient faster-than-exponential bubble growth with approximate log-periodic behavior and give analytical arguments why this follows from our framework. The model accounts well for the behavior of traders and for the price dynamics that developed during the dotcom bubble in 1995-2000. Momentum strategies are shown to be transiently profitable, supporting these strategies as enhancing herding behavior.

preprint2014arXiv

Utility Rate Equations of Group Population Dynamics in Biological and Social Systems

We present a novel system of equations to describe the evolution of self-organized structured societies (biological or human) composed of several trait groups. The suggested approach is based on the combination of ideas employed in the theory of biological populations, system theory, and utility theory. The evolution equations are defined as utility rate equations, whose parameters are characterized by the utility of each group with respect to the society as a whole and by the mutual utilities of groups with respect to each other. We analyze in detail the cases of two groups (cooperators and defectors) and of three groups (cooperators, defectors, and regulators) and find that, in a self-organized society, neither defectors nor regulators can overpass the maximal fractions of about 10% each. This is in agreement with the data for bee and ant colonies. The classification of societies by their distance from equilibrium is proposed. We apply the formalism to rank the countries according to the introduced metric quantifying their relative stability, which depends on the cost of defectors and regulators as well as their respective population fractions. We find a remarkable concordance with more standard economic ranking based, for instance, on GDP per capita.

preprint2013arXiv

Conditions for Quantum Interference in Cognitive Sciences

We present a general classification of the conditions under which cognitive science, concerned e.g. with decision making, requires the use of quantum theoretical notions. The analysis is done in the frame of the mathematical approach based on the theory of quantum measurements. We stress that quantum effects in cognition can arise only when decisions are made under uncertainty. Conditions for the appearance of quantum interference in cognitive sciences and the conditions when interference cannot arise are formulated.

preprint2013arXiv

Exploring the limits of safety analysis in complex technological systems

From biotechnology to cyber-risks, most extreme technological risks cannot be reliably estimated from historical statistics. Therefore, engineers resort to predictive methods, such as fault/event trees in the framework of probabilistic safety assessment (PSA), which consists in developing models to identify triggering events, potential accident scenarios, and estimate their severity and frequency. However, even the best safety analysis struggles to account for evolving risks resulting from inter-connected networks and cascade effects. Taking nuclear risks as an example, the predicted plant-specific distribution of losses is found to be significantly underestimated when compared with available empirical records. Using a novel database of 99 events with losses larger than $50'000 constructed by Sovacool, we document a robust power law distribution with tail exponent mu \approx 0.7. A simple cascade model suggests that the classification of the different possible safety regimes is intrinsically unstable in the presence of cascades. Additional continuous development and validation, making the best use of the experienced realized incidents, near misses and accidents, is urgently needed to address the existing known limitations of PSA when aiming at the estimation of total risks.

preprint2013arXiv

Mode interference in quantum joint probabilities for multimode Bose-condensed systems

The method of defining quantum joint probabilities of two events is applied to a multimode system of trapped Bose-condensed atoms. The coherent modes are generated by modulating the trapping potential with an alternating field with a frequency in resonance with a transition frequency between two coherent modes. The system is subjected to nondestructive measurements randomly influencing the phase difference of the modes. The joint probability of observing two events, corresponding to the observation of two states, is calculated by employing the channel-state duality. The interference term in the joint probability can arise when the composite events form an entangled prospect and the system state is also entangled. This interference is due to the uncertainty induced by nondestructive measurements.

preprint2013arXiv

Quantum probabilities of composite events in quantum measurements with multimode states

The problem of defining quantum probabilities of composite events is considered. This problem is of high importance for the theory of quantum measurements and for quantum decision theory that is a part of measurement theory. We show that the Luders probability of consecutive measurements is a transition probability between two quantum states and that this probability cannot be treated as a quantum extension of the classical conditional probability. The Wigner distribution is shown to be a weighted transition probability that cannot be accepted as a quantum extension of the classical joint probability. We suggest the definition of quantum joint probabilities by introducing composite events in multichannel measurements. The notion of measurements under uncertainty is defined. We demonstrate that the necessary condition for the mode interference is the entanglement of the composite prospect together with the entanglement of the composite statistical state. As an illustration, we consider an example of a quantum game. A special attention is payed to the application of the approach to systems with multimode states, such as atoms, molecules, quantum dots, or trapped Bose-condensed atoms with several coherent modes.

preprint2013arXiv

Super-linear scaling of offsprings at criticality in branching processes

For any branching process, we demonstrate that the typical total number $r_{\rm mp}(ντ)$ of events triggered over all generations within any sufficiently large time window $τ$ exhibits, at criticality, a super-linear dependence $r_{\rm mp}(ντ) \sim (ντ)^γ$ (with $γ>1$) on the total number $ντ$ of the immigrants arriving at the Poisson rate $ν$. In branching processes in which immigrants (or sources) are characterized by fertilities distributed according to an asymptotic power law tail with tail exponent $1 < γ\leqslant 2$, the exponent of the super-linear law for $r_{\rm mp}(ντ)$ is identical to the exponent $γ$ of the distribution of fertilities. For $γ>2$ and for standard branching processes without power law distribution of fertilities, $r_{\rm mp}(ντ) \sim (ντ)^2$. This novel scaling law replaces and tames the divergence $ντ/(1-n)$ of the mean total number ${\bar R}_t(τ)$ of events, as the branching ratio (defined as the average number of triggered events of first generation per source) tends to 1. The derivation uses the formalism of generating probability functions. The corresponding prediction is confirmed by numerical calculations and an heuristic derivation enlightens its underlying mechanism. We also show that ${\bar R}_t(τ)$ is always linear in $ντ$ even at criticality ($n=1$). Our results thus illustrate the fundamental difference between the mean total number, which is controlled by a few extremely rare realizations, and the typical behavior represented by $r_{\rm mp}(ντ)$.

preprint2012arXiv

A simple microstructure return model explaining microstructure noise and Epps effects

We present a simple microstructure model of financial returns that combines (i) the well-known ARFIMA process applied to tick-by-tick returns, (ii) the bid-ask bounce effect, (iii) the fat tail structure of the distribution of returns and (iv) the non-Poissonian statistics of inter-trade intervals. This model allows us to explain both qualitatively and quantitatively important stylized facts observed in the statistics of microstructure returns, including the short-ranged correlation of returns, the long-ranged correlations of absolute returns, the microstructure noise and Epps effects. According to the microstructure noise effect, volatility is a decreasing function of the time scale used to estimate it. Paradoxically, the Epps effect states that cross correlations between asset returns are increasing functions of the time scale at which the returns are estimated. The microstructure noise is explained as the result of the negative return correlations inherent in the definition of the bid-ask bounce component (ii). In the presence of a genuine correlation between the returns of two assets, the Epps effect is due to an average statistical overlap of the momentum of the returns of the two assets defined over a finite time scale in the presence of the long memory process (i).

preprint2012arXiv

Dragon-kings: mechanisms, statistical methods and empirical evidence

This introductory article presents the special Discussion and Debate volume "From black swans to dragon-kings, is there life beyond power laws?" published in Eur. Phys. J. Special Topics in May 2012. We summarize and put in perspective the contributions into three main themes: (i) mechanisms for dragon-kings, (ii) detection of dragon-kings and statistical tests and (iii) empirical evidence in a large variety of natural and social systems. Overall, we are pleased to witness significant advances both in the introduction and clarification of underlying mechanisms and in the development of novel efficient tests that demonstrate clear evidence for the presence of dragon-kings in many systems. However, this positive view should be balanced by the fact that this remains a very delicate and difficult field, if only due to the scarcity of data as well as the extraordinary important implications with respect to hazard assessment, risk control and predictability.

preprint2012arXiv

Extreme events in population dynamics with functional carrying capacity

A class of models is introduced describing the evolution of population species whose carrying capacities are functionals of these populations. The functional dependence of the carrying capacities reflects the fact that the correlations between populations can be realized not merely through direct interactions, as in the usual predator-prey Lotka-Volterra model, but also through the influence of species on the carrying capacities of each other. This includes the self-influence of each kind of species on its own carrying capacity with delays. Several examples of such evolution equations with functional carrying capacities are analyzed. The emphasis is given on the conditions under which the solutions to the equations display extreme events, such as finite-time death and finite-time singularity. Any destructive action of populations, whether on their own carrying capacity or on the carrying capacities of co-existing species, can lead to the instability of the whole population that is revealed in the form of the appearance of extreme events, finite-time extinctions or booms followed by crashes.

preprint2012arXiv

Fertility Heterogeneity as a Mechanism for Power Law Distributions of Recurrence Times

We study the statistical properties of recurrence times in the self-excited Hawkes conditional Poisson process, the simplest extension of the Poisson process that takes into account how the past events influence the occurrence of future events. Specifically, we analyze the impact of the power law distribution of fertilities with exponent α, where the fertility of an event is the number of aftershocks of first generation that it triggers, on the probability distribution function (pdf) f(τ) of the recurrence times τbetween successive events. The other input of the model is an exponential Omori law quantifying the pdf of waiting times between an event and its first generation aftershocks, whose characteristic time scale is taken as our time unit. At short time scales, we discover two intermediate power law asymptotics, f(τ) ~ τ^{-(2-α)} for τ<< τ_c and f(τ) ~ τ^{-α} for τ_c << τ<< 1, where τ_c is associated with the self-excited cascades of aftershocks. For 1 << τ<< 1/ν, we find a constant plateau f(τ) ~ const, while at long times, 1/ν< τ, f(τ) ~ e^{-ντ} has an exponential tail controlled by the arrival rate νof exogenous events. These results demonstrate a novel mechanism for the generation of power laws in the distribution of recurrence times, which results from a power law distribution of fertilities in the presence of self-excitation and cascades of triggering.

preprint2012arXiv

Modeling symbiosis by interactions through species carrying capacities

We introduce a mathematical model of symbiosis between different species by taking into account the influence of each species on the carrying capacities of the others. The modeled entities can pertain to biological and ecological societies or to social, economic and financial societies. Our model includes three basic types: symbiosis with direct mutual interactions, symbiosis with asymmetric interactions, and symbiosis without direct interactions. In all cases, we provide a complete classification of all admissible dynamical regimes. The proposed model of symbiosis turned out to be very rich, as it exhibits four qualitatively different regimes: convergence to stationary states, unbounded exponential growth, finite-time singularity, and finite-time death or extinction of species.

preprint2012arXiv

Statistical Outliers and Dragon-Kings as Bose-Condensed Droplets

A theory of exceptional extreme events, characterized by their abnormal sizes compared with the rest of the distribution, is presented. Such outliers, called "dragon-kings", have been reported in the distribution of financial drawdowns, city-size distributions (e.g., Paris in France and London in the UK), in material failure, epileptic seizure intensities, and other systems. Within our theory, the large outliers are interpreted as droplets of Bose-Einstein condensate: the appearance of outliers is a natural consequence of the occurrence of Bose-Einstein condensation controlled by the relative degree of attraction, or utility, of the largest entities. For large populations, Zipf's law is recovered (except for the dragon-king outliers). The theory thus provides a parsimonious description of the possible coexistence of a power law distribution of event sizes (Zipf's law) and dragon-king outliers.

preprint2011arXiv

Decision Theory with Prospect Interference and Entanglement

We present a novel variant of decision making based on the mathematical theory of separable Hilbert spaces. This mathematical structure captures the effect of superposition of composite prospects, including many incorporated intentions, which allows us to describe a variety of interesting fallacies and anomalies that have been reported to particularize the decision making of real human beings. The theory characterizes entangled decision making, non-commutativity of subsequent decisions, and intention interference. We demonstrate how the violation of the Savage's sure-thing principle, known as the disjunction effect, can be explained quantitatively as a result of the interference of intentions, when making decisions under uncertainty. The disjunction effects, observed in experiments, are accurately predicted using a theorem on interference alternation that we derive, which connects aversion-to-uncertainty to the appearance of negative interference terms suppressing the probability of actions. The conjunction fallacy is also explained by the presence of the interference terms. A series of experiments are analysed and shown to be in excellent agreement with a priori evaluation of interference effects. The conjunction fallacy is also shown to be a sufficient condition for the disjunction effect and novel experiments testing the combined interplay between the two effects are suggested.

preprint2011arXiv

Generating Functions and Stability Study of Multivariate Self-Excited Epidemic Processes

We present a stability study of the class of multivariate self-excited Hawkes point processes, that can model natural and social systems, including earthquakes, epileptic seizures and the dynamics of neuron assemblies, bursts of exchanges in social communities, interactions between Internet bloggers, bank network fragility and cascading of failures, national sovereign default contagion, and so on. We present the general theory of multivariate generating functions to derive the number of events over all generations of various types that are triggered by a mother event of a given type. We obtain the stability domains of various systems, as a function of the topological structure of the mutual excitations across different event types. We find that mutual triggering tends to provide a significant extension of the stability (or subcritical) domain compared with the case where event types are decoupled, that is, when an event of a given type can only trigger events of the same type.

preprint2011arXiv

Hierarchy of Temporal Responses of Multivariate Self-Excited Epidemic Processes

We present the first exact analysis of some of the temporal properties of multivariate self-excited Hawkes conditional Poisson processes, which constitute powerful representations of a large variety of systems with bursty events, for which past activity triggers future activity. The term "multivariate" refers to the property that events come in different types, with possibly different intra- and inter-triggering abilities. We develop the general formalism of the multivariate generating moment function for the cumulative number of first-generation and of all generation events triggered by a given mother event (the "shock") as a function of the current time $t$. This corresponds to studying the response function of the process. A variety of different systems have been analyzed. In particular, for systems in which triggering between events of different types proceeds through a one-dimension directed or symmetric chain of influence in type space, we report a novel hierarchy of intermediate asymptotic power law decays $\sim 1/t^{1-(m+1)θ}$ of the rate of triggered events as a function of the distance $m$ of the events to the initial shock in the type space, where $0 < θ<1$ for the relevant long-memory processes characterizing many natural and social systems. The richness of the generated time dynamics comes from the cascades of intermediate events of possibly different kinds, unfolding via a kind of inter-breeding genealogy.

preprint2011arXiv

On the Correct Use of Statistical Tests: Reply to "Lies, damned lies and statistics (in Geology)"

In a Forum published in EOS Transactions AGU (2009) entitled "Lies, damned lies and statistics (in Geology)", Vermeesch (2009) claims that "statistical significant is not the same as geological significant", in other words, statistical tests may be misleading. In complete contradiction, we affirm that statistical tests are always informative. We detail the several mistakes of Vermeesch in his initial paper and in his comments to our reply. The present text is developed in the hope that it can serve as an illuminating pedagogical exercise for students and lecturers to learn more about the subtleties, richness and power of the science of statistics.

preprint2011arXiv

Robust Statistical Tests of Dragon-Kings beyond Power Law Distributions

We ask the question whether it is possible to diagnose the existence of "Dragon-Kings" (DK), namely anomalous observations compared to a power law background distribution of event sizes. We present two new statistical tests, the U-test and the DK-test, aimed at identifying the existence of even a single anomalous event in the tail of the distribution of just a few tens of observations. The DK-test in particular is derived such that the p-value of its statistic is independent of the exponent characterizing the null hypothesis. We demonstrate how to apply these two tests on the distributions of cities and of agglomerations in a number of countries. We find the following evidence for Dragon-Kings: London in the distribution of city sizes of Great Britain; Moscow and St-Petersburg in the distribution of city sizes in the Russian Federation; and Paris in the distribution of agglomeration sizes in France. True negatives are also reported, for instance the absence of Dragon-Kings in the distribution of cities in Germany.

preprint2011arXiv

Time-Bridge Estimators of Integrated Variance

We present a set of log-price integrated variance estimators, equal to the sum of open-high-low-close bridge estimators of spot variances within $n$ subsequent time-step intervals. The main characteristics of some of the introduced estimators is to take into account the information on the occurrence times of the high and low values. The use of the high's and low's of the bridge associated with the original process makes the estimators significantly more efficient that the standard realized variance estimators and its generalizations. Adding the information on the occurrence times of the high and low values improves further the efficiency of the estimators, much above those of the well-known realized variance estimator and those derived from the sum of Garman and Klass spot variance estimators. The exact analytical results are derived for the case where the underlying log-price process is an Itô stochastic process. Our results suggests more efficient ways to record financial prices at intermediate frequencies.

preprint2010arXiv

Entanglement production in quantum decision making

The quantum decision theory introduced recently is formulated as a quantum theory of measurement. It describes prospect states represented by complex vectors of a Hilbert space over a prospect lattice. The prospect operators, acting in this space, form an involutive bijective algebra. A measure is defined for quantifying the entanglement produced by the action of prospect operators. This measure characterizes the level of complexity of prospects involved in decision making. An explicit expression is found for the maximal entanglement produced by the operators of multimode prospects.

preprint2010arXiv

Mathematical Structure of Quantum Decision Theory

One of the most complex systems is the human brain whose formalized functioning is characterized by decision theory. We present a "Quantum Decision Theory" of decision making, based on the mathematical theory of separable Hilbert spaces. This mathematical structure captures the effect of superposition of composite prospects, including many incorporated intentions, which allows us to explain a variety of interesting fallacies and anomalies that have been reported to particularize the decision making of real human beings. The theory describes entangled decision making, non-commutativity of subsequent decisions, and intention interference of composite prospects. We demonstrate how the violation of the Savage's sure-thing principle (disjunction effect) can be explained as a result of the interference of intentions, when making decisions under uncertainty. The conjunction fallacy is also explained by the presence of the interference terms. We demonstrate that all known anomalies and paradoxes, documented in the context of classical decision theory, are reducible to just a few mathematical archetypes, all of which finding straightforward explanations in the frame of the developed quantum approach.

preprint2010arXiv

Processing Information in Quantum Decision Theory

A survey is given summarizing the state of the art of describing information processing in Quantum Decision Theory, which has been recently advanced as a novel variant of decision making, based on the mathematical theory of separable Hilbert spaces. This mathematical structure captures the effect of superposition of composite prospects, including many incorporated intended actions. The theory characterizes entangled decision making, non-commutativity of subsequent decisions, and intention interference. The self-consistent procedure of decision making, in the frame of the quantum decision theory, takes into account both the available objective information as well as subjective contextual effects. This quantum approach avoids any paradox typical of classical decision theory. Conditional maximization of entropy, equivalent to the minimization of an information functional, makes it possible to connect the quantum and classical decision theories, showing that the latter is the limit of the former under vanishing interference terms.

preprint2010arXiv

Reverse Engineering Financial Markets with Majority and Minority Games using Genetic Algorithms

Using virtual stock markets with artificial interacting software investors, aka agent-based models (ABMs), we present a method to reverse engineer real-world financial time series. We model financial markets as made of a large number of interacting boundedly rational agents. By optimizing the similarity between the actual data and that generated by the reconstructed virtual stock market, we obtain parameters and strategies, which reveal some of the inner workings of the target stock market. We validate our approach by out-of-sample predictions of directional moves of the Nasdaq Composite Index.

preprint2010arXiv

Zipf's law and maximum sustainable growth

Zipf's law states that the number of firms with size greater than S is inversely proportional to S. Most explanations start with Gibrat's rule of proportional growth but require additional constraints. We show that Gibrat's rule, at all firm levels, yields Zipf's law under a balance condition between the effective growth rate of incumbent firms (which includes their possible demise) and the growth rate of investments in entrant firms. Remarkably, Zipf's law is the signature of the long-term optimal allocation of resources that ensures the maximum sustainable growth rate of an economy.

preprint2009arXiv

A Consistent Model of `Explosive' Financial Bubbles With Mean-Reversing Residuals

We present a self-consistent model for explosive financial bubbles, which combines a mean-reverting volatility process and a stochastic conditional return which reflects nonlinear positive feedbacks and continuous updates of the investors' beliefs and sentiments. The conditional expected returns exhibit faster-than-exponential acceleration decorated by accelerating oscillations, called "log-periodic power law." Tests on residuals show a remarkable low rate (0.2%) of false positives when applied to a GARCH benchmark. When tested on the S&P500 US index from Jan. 3, 1950 to Nov. 21, 2008, the model correctly identifies the bubbles ending in Oct. 1987, in Oct. 1997, in Aug. 1998 and the ITC bubble ending on the first quarter of 2000. Different unit-root tests confirm the high relevance of the model specification. Our model also provides a diagnostic for the duration of bubbles: applied to the period before Oct. 1987 crash, there is clear evidence that the bubble started at least 4 years earlier. We confirm the validity and universality of the volatility-confined LPPL model on seven other major bubbles that have occurred in the World in the last two decades. Using Bayesian inference, we find a very strong statistical preference for our model compared with a standard benchmark, in contradiction with Chang and Feigenbaum (2006) which used a unit-root model for residuals.

preprint2009arXiv

Distribution of Maximum Earthquake Magnitudes in Future Time Intervals, Application to the Seismicity of Japan (1923-2007)

We modify the new method for the statistical estimation of the tail distribution of earthquake seismic moments introduced by Pisarenko et al. [2009] and apply it to the earthquake catalog of Japan (1923-2007). The method is based on the two main limit theorems of the theory of extreme values and on the derived duality between the Generalized Pareto Distribution (GPD) and Generalized Extreme Value distribution (GEV). We obtain the distribution of maximum earthquake magnitudes in future time intervals of arbitrary duration tau. This distribution can be characterized by its quantile Qq(tau) at any desirable statistical level q. The quantile Qq(tau) provides a much more stable and robust characteristic than the traditional absolute maximum magnitude Mmax (Mmax can be obtained as the limit of Qq(tau) as q tends to 1, and tau tends to infinity). The best estimates of the parameters governing the distribution of Qq(tay) for Japan (1923-2007) are the following: Form parameter for GEV = -0.1901 +- 0.0717; position parameter GEV(tau=200)= 6.3387 +- 0.0380; spread parameter for GEV(tau=200)= 0.5995 +- 0.0223; Q_0.90,GEV(tau=10)= 8.34 +- 0.32. We also estimate Qq(tau) for a set of q-values and future time periods in the range for tau between 1 and 50 years from 2007. For comparison, the absolute maximum estimate Mmax from GEV, which is equal to 9.57 +- 0.86, has a scatter more than twice that of the 90 percent quantile Q_{0.90,GEV}(tau=10) of the maximum magnitude over the next 10 years counted from 2007.

preprint2009arXiv

Effects of Diversity and Procrastination in Priority Queuing Theory: the Different Power Law Regimes

Empirical analysis show that, after the update of a browser, the publication of the vulnerability of a software, or the discovery of a cyber worm, the fraction of computers still using the older version, or being not yet patched, or exhibiting worm activity decays as power laws $\sim 1/t^α$ with $0 < α\leq 1$ over time scales of years. We present a simple model for this persistence phenomenon framed within the standard priority queuing theory, of a target task which has the lowest priority compared with all other tasks that flow on the computer of an individual. We identify a "time deficit" control parameter $β$ and a bifurcation to a regime where there is a non-zero probability for the target task to never be completed. The distribution of waiting time ${\cal T}$ till the completion of the target task has the power law tail $\sim 1/t^{1/2}$, resulting from a first-passage solution of an equivalent Wiener process. Taking into account a diversity of time deficit parameters in a population of individuals, the power law tail is changed into $1/t^α$ with $α\in(0.5,\infty)$, including the well-known case $1/t$. We also study the effect of "procrastination", defined as the situation in which the target task may be postponed or delayed even after the individual has solved all other pending tasks. This new regime provides an explanation for even slower apparent decay and longer persistence.

preprint2009arXiv

Generation-by-Generation Dissection of the Response Function in Long Memory Epidemic Processes

In a number of natural and social systems, the response to an exogenous shock relaxes back to the average level according to a long-memory kernel $\sim 1/t^{1+θ}$ with $0 \leq θ<1$. In the presence of an epidemic-like process of triggered shocks developing in a cascade of generations at or close to criticality, this "bare" kernel is renormalized into an even slower decaying response function $\sim 1/t^{1-θ}$. Surprisingly, this means that the shorter the memory of the bare kernel (the larger $1+θ$), the longer the memory of the response function (the smaller $1-θ$). Here, we present a detailed investigation of this paradoxical behavior based on a generation-by-generation decomposition of the total response function, the use of Laplace transforms and of "anomalous" scaling arguments. The paradox is explained by the fact that the number of triggered generations grows anomalously with time at $\sim t^θ$ so that the contributions of active generations up to time $t$ more than compensate the shorter memory associated with a larger exponent $θ$. This anomalous scaling results fundamentally from the property that the expected waiting time is infinite for $0 \leq θ\leq 1$. The techniques developed here are also applied to the case $θ>1$ and we find in this case that the total renormalized response is a {\bf constant} for $t < 1/(1-n)$ followed by a cross-over to $\sim 1/t^{1+θ}$ for $t \gg 1/(1-n)$.

preprint2009arXiv

Gibrat's law for cities: uniformly most powerful unbiased test of the Pareto against the lognormal

We address the general problem of testing a power law distribution versus a log-normal distribution in statistical data. This general problem is illustrated on the distribution of the 2000 US census of city sizes. We provide definitive results to close the debate between Eeckhout (2004, 2009) and Levy (2009) on the validity of Zipf's law, which is the special Pareto law with tail exponent 1, to describe the tail of the distribution of U.S. city sizes. Because the origin of the disagreement between Eeckhout and Levy stems from the limited power of their tests, we perform the {\em uniformly most powerful unbiased test} for the null hypothesis of the Pareto distribution against the lognormal. The $p$-value and Hill's estimator as a function of city size lower threshold confirm indubitably that the size distribution of the 1000 largest cities or so, which include more than half of the total U.S. population, is Pareto, but we rule out that the tail exponent, estimated to be $1.4 \pm 0.1$, is equal to 1. For larger ranks, the $p$-value becomes very small and Hill's estimator decays systematically with decreasing ranks, qualifying the lognormal distribution as the better model for the set of smaller cities. These two results reconcile the opposite views of Eeckhout (2004, 2009) and Levy (2009). We explain how Gibrat's law of proportional growth underpins both the Pareto and lognormal distributions and stress the key ingredient at the origin of their difference in standard stochastic growth models of cities \cite{Gabaix99,Eeckhout2004}.

preprint2009arXiv

Limits of Declustering Methods for Disentangling Exogenous from Endogenous Events in Time Series with Foreshocks, Main shocks and Aftershocks

Many time series in natural and social sciences can be seen as resulting from an interplay between exogenous influences and an endogenous organization. We use a simple (ETAS) model of events occurring sequentially, in which future events are influenced (partially triggered) by past events to ask the question of how well can one disentangle the exogenous events from the endogenous ones. We apply both model-dependant and model-independent stochastic declustering methods to reconstruct the tree of ancestry and estimate key parameters. In contrast with previously reported positive results, we have to conclude that declustered catalogs are rather unreliable for the synthetic catalogs that we have investigated, which contains of the order of thousands of events, typical of realistic applications. The estimated rates of exogenous events suffer from large errors. The key branching ratio $n$, quantifying the fraction of events that have been triggered by previous events, is also badly estimated in general from declustered catalogs. We find however that the errors tend to be smaller and perhaps acceptable in some cases for small triggering efficiency and branching ratios. The high level of randomness together with the very long memory makes the stochastic reconstruction of trees of ancestry and the estimation of the key parameters perhaps intrinsically unreliable for long memory processes. For shorter memories (larger "bare" Omori exponent), the results improve significantly.

preprint2009arXiv

Physics of risk and uncertainty in quantum decision making

The Quantum Decision Theory, developed recently by the authors, is applied to clarify the role of risk and uncertainty in decision making and in particular in relation to the phenomenon of dynamic inconsistency. By formulating this notion in precise mathematical terms, we distinguish three types of inconsistency: time inconsistency, planning paradox, and inconsistency occurring in some discounting effects. While time inconsistency is well accounted for in classical decision theory, the planning paradox is in contradiction with classical utility theory. It finds a natural explanation in the frame of the Quantum Decision Theory. Different types of discounting effects are analyzed and shown to enjoy a straightforward explanation within the suggested theory. We also introduce a general methodology based on self-similar approximation theory for deriving the evolution equations for the probabilities of future prospects. This provides a novel classification of possible discount factors, which include the previously known cases (exponential or hyperbolic discounting), but also predicts a novel class of discount factors that decay to a strictly positive constant for very large future time horizons. This class may be useful to deal with very long-term discounting situations associated with intergenerational public policy choices, encompassing issues such as global warming and nuclear waste disposal.

preprint2009arXiv

Punctuated evolution due to delayed carrying capacity

A new delay equation is introduced to describe the punctuated evolution of complex nonlinear systems. A detailed analytical and numerical investigation provides the classification of all possible types of solutions for the dynamics of a population in the four main regimes dominated respectively by: (i) gain and competition, (ii) gain and cooperation, (iii) loss and competition and (iv) loss and cooperation. Our delay equation may exhibit bistability in some parameter range, as well as a rich set of regimes, including monotonic decay to zero, smooth exponential growth, punctuated unlimited growth, punctuated growth or alternation to a stationary level, oscillatory approach to a stationary level, sustainable oscillations, finite-time singularities as well as finite-time death.

preprint2009arXiv

Scheme of thinking quantum systems

A general approach describing quantum decision procedures is developed. The approach can be applied to quantum information processing, quantum computing, creation of artificial quantum intelligence, as well as to analyzing decision processes of human decision makers. Our basic point is to consider an active quantum system possessing its own strategic state. Processing information by such a system is analogous to the cognitive processes associated to decision making by humans. The algebra of probability operators, associated with the possible options available to the decision maker, plays the role of the algebra of observables in quantum theory of measurements. A scheme is advanced for a practical realization of decision procedures by thinking quantum systems. Such thinking quantum systems can be realized by using spin lattices, systems of magnetic molecules, cold atoms trapped in optical lattices, ensembles of quantum dots, or multilevel atomic systems interacting with electromagnetic field.

preprint2009arXiv

Stable States of Biological Organisms

A novel model of biological organisms is advanced, treating an organism as a self-consistent system subject to a pathogen flux. The principal novelty of the model is that it describes not some parts, but a biological organism as a whole. The organism is modeled by a five-dimensional dynamical system. The organism homeostasis is described by the evolution equations for five interacting components: healthy cells, ill cells, innate immune cells, specific immune cells, and pathogens. The stability analysis demonstrates that, in a wide domain of the parameter space, the system exhibits robust structural stability. There always exist four stable stationary solutions characterizing four qualitatively differing states of the organism: alive state, boundary state, critical state, and dead state.

preprint2008arXiv

Empirical Tests of Zipf's law Mechanism In Open Source Linux Distribution

The evolution of open source software projects in Linux distributions offers a remarkable example of a growing complex self-organizing adaptive system, exhibiting Zipf's law over four full decades. We present three tests of the usually assumed ingredients of stochastic growth models that have been previously conjectured to be at the origin of Zipf's law: (i) the growth observed between successive releases of the number of in-directed links of packages obeys Gibrat's law of proportional growth; (ii) the average growth increment of the number of in-directed links of packages over a time interval $Δt$ is proportional to $Δt$, while its standard deviation is proportional to $\sqrt{Δt}$; (iii) the distribution of the number of in-directed links of new packages appearing in evolving versions of Debian Linux distributions has a tail thinner than Zipf's law, with an exponent which converges to the Zipf's law value 1 as the time $Δt$ between releases increases.

preprint2008arXiv

Interdisciplinarity in Socio-economics, mathematical analysis and predictability of complex systems

In this essay, I attempt to provide supporting evidence as well as some balance for the thesis on `Transforming socio-economics with a new epistemology' presented by Hollingworth and Mueller (2008). First, I review a personal highlight of my own scientific path that illustrates the power of interdisciplinarity as well as unity of the mathematical description of natural and social processes. I also argue against the claim that complex systems are in general `not susceptible to mathematical analysis, but must be understood by letting them evolve over time or with simulation analysis'. Moreover, I present evidence of the limits of the claim that scientists working within Science II do not make predictions about the future because it is too complex. I stress the potentials for a third `Quantum Science' and its associated conceptual and philosophical revolutions, and finally point out some limits of the `new' theory of networks.

preprint2008arXiv

Theory of Zipf's Law and of General Power Law Distributions with Gibrat's law of Proportional Growth

We summarize a book under publication with his title written by the three present authors, on the theory of Zipf's law, and more generally of power laws, driven by the mechanism of proportional growth. The preprint is available upon request from the authors. For clarity, consistence of language and conciseness, we discuss the origin and conditions of the validity of Zipf's law using the terminology of firms' asset values. We use firms at the entities whose size distributions are to be explained. It should be noted, however, that most of the relations discussed in this book, especially the intimate connection between Zipf's and Gilbrat's laws, underlie Zipf's law in diverse scientific areas. The same models and variations thereof can be straightforwardly applied to any of the other domains of application.

preprint2007arXiv

Magnitude Uncertainties Impact Seismic Rate Estimates, Forecasts and Predictability Experiments

The Collaboratory for the Study of Earthquake Predictability (CSEP) aims to prospectively test time-dependent earthquake probability forecasts on their consistency with observations. To compete, time-dependent seismicity models are calibrated on earthquake catalog data. But catalogs contain much observational uncertainty. We study the impact of magnitude uncertainties on rate estimates in clustering models, on their forecasts and on their evaluation by CSEP's consistency tests. First, we quantify magnitude uncertainties. We find that magnitude uncertainty is more heavy-tailed than a Gaussian, such as a double-sided exponential distribution, with scale parameter nu_c=0.1 - 0.3. Second, we study the impact of such noise on the forecasts of a simple clustering model which captures the main ingredients of popular short term models. We prove that the deviations of noisy forecasts from an exact forecast are power law distributed in the tail with exponent alpha=1/(a*nu_c), where a is the exponent of the productivity law of aftershocks. We further prove that the typical scale of the fluctuations remains sensitively dependent on the specific catalog. Third, we study how noisy forecasts are evaluated in CSEP consistency tests. Noisy forecasts are rejected more frequently than expected for a given confidence limit. The Poisson assumption of the consistency tests is inadequate for short-term forecast evaluations. To capture the idiosyncrasies of each model together with any propagating uncertainties, the forecasts need to specify the entire likelihood distribution of seismic rates.

preprint2006arXiv

Algorithm for Model Validation: Theory and Applications

Validation is often defined as the process of determining the degree to which a model is an accurate representation of the real world from the perspective of its intended uses. Validation is crucial as industries and governments depend increasingly on predictions by computer models to justify their decisions. We propose to formulate the validation of a given model as an iterative construction process that mimics the often implicit process occurring in the minds of scientists. We offer a formal representation of the progressive build-up of trust in the model. We thus replace static claims on the impossibility of validating a given model by a dynamic process of constructive approximation. This approach is better adapted to the fuzzy, coarse-grained nature of validation. Our procedure factors in the degree of redundancy versus novelty of the experiments used for validation as well as the degree to which the model predicts the observations. We illustrate the new methodology first with the maturation of Quantum Mechanics as the arguably best established physics theory and then with several concrete examples drawn from some of our primary scientific interests: a cellular automaton model for earthquakes, a multifractal random walk model for financial time series, an anomalous diffusion model for solar radiation transport in the cloudy atmosphere, and a computational fluid dynamics code for the Richtmyer-Meshkov instability.

preprint2005arXiv

Celebrating the Physics in Geophysics

As 2005, the International Year of Physics, comes to an end, two physicists working primarily in geophysical research reflect on how geophysics is not an applied physics. Although geophysics has certainly benefited from progress in physics and sometimes emulated the reductionist program of mainstream physics, it has also educated the physics community about some of the generic behaviors of strongly nonlinear systems. Dramatic examples are the insights we have gained into the ``emergent'' phenomena of chaos, cascading instabilities, turbulence, self-organization, fractal structure, power-law variability, anomalous scaling, threshold dynamics, creep, fracture, and so on. In all of these examples, relatively simple models have been able to explain the recurring features of apparently very complex signals and fields. It appears that the future of the intricate relation between physics and geophysics will be as exciting as its past has been characterized by a mutual fascination. Physics departments in our universities should capitalize on this trend to attract and retain young talent motivated to address problems that really matter for the future of the planet. A pressing topic with huge impact on populations and that is challenging enough for both physics and geophysics communities to work together like never before is the understanding and prediction of extreme events.

preprint2005arXiv

Statistical Physics of Rupture in Heterogeneous Media

The damage and fracture of materials are technologically of enormous interest due to their economic and human cost. They cover a wide range of phenomena like e.g. cracking of glass, aging of concrete, the failure of fiber networks in the formation of paper and the breaking of a metal bar subject to an external load. Failure of composite systems is of utmost importance in naval, aeronautics and space industry. By the term composite, we refer to materials with heterogeneous microscopic structures and also to assemblages of macroscopic elements forming a super-structure. Chemical and nuclear plants suffer from cracking due to corrosion either of chemical or radioactive origin, aided by thermal and/or mechanical stress. Despite the large amount of experimental data and the considerable effort that has been undertaken by material scientists, many questions about fracture have not been answered yet. There is no comprehensive understanding of rupture phenomena but only a partial classification in restricted and relatively simple situations. This lack of fundamental understanding is indeed reflected in the absence of reliable prediction methods for rupture, based on a suitable monitoring of the stressed system. Not only is there a lack of non-empirical understanding of the reliability of a system, but also the empirical laws themselves have often limited value. The difficulties stem from the complex interplay between heterogeneities and modes of damage and the possible existence of a hierarchy of characteristic scales (static and dynamic). The paper presents a review of recent efforts from the statistical physics community to address these points.

preprint2003arXiv

Predictability in the ETAS Model of Interacting Triggered Seismicity

As part of an effort to develop a systematic methodology for earthquake forecasting, we use a simple model of seismicity based on interacting events which may trigger a cascade of earthquakes, known as the Epidemic-Type Aftershock Sequence model (ETAS). The ETAS model is constructed on a bare (unrenormalized) Omori law, the Gutenberg-Richter law and the idea that large events trigger more numerous aftershocks. For simplicity, we do not use the information on the spatial location of earthquakes and work only in the time domain. We offer an analytical approach to account for the yet unobserved triggered seismicity adapted to the problem of forecasting future seismic rates at varying horizons from the present. Tests presented on synthetic catalogs validate strongly the importance of taking into account all the cascades of still unobserved triggered events in order to predict correctly the future level of seismicity beyond a few minutes. We find a strong predictability if one accepts to predict only a small fraction of the large-magnitude targets. However, the probability gains degrade fast when one attempts to predict a larger fraction of the targets. This is because a significant fraction of events remain uncorrelated from past seismicity. This delineates the fundamental limits underlying forecasting skills, stemming from an intrinsic stochastic component in these interacting triggered seismicity models.

preprint2003arXiv

Slider-Block Friction Model for Landslides: Application to Vaiont and La Clapiere Landslides

Accelerating displacements preceding some catastrophic landslides have been found empirically to follow a time-to-failure power law, corresponding to a finite-time singularity of the velocity $v \sim 1/(t_c-t)$ [{\it Voight}, 1988]. Here, we provide a physical basis for this phenomenological law based on a slider-block model using a state and velocity dependent friction law established in the laboratory and used to model earthquake friction. This physical model accounts for and generalizes Voight's observation: depending on the ratio $B/A$ of two parameters of the rate and state friction law and on the initial frictional state of the sliding surfaces characterized by a reduced parameter $x_i$, four possible regimes are found. Two regimes can account for an acceleration of the displacement. We use the slider-block friction model to analyze quantitatively the displacement and velocity data preceding two landslides, Vaiont and La Clapière. The Vaiont landslide was the catastrophic culmination of an accelerated slope velocity. La Clapière landslide was characterized by a peak of slope acceleration that followed decades of ongoing accelerating displacements, succeeded by a restabilizing phase. Our inversion of the slider-block model on these data sets shows good fits and suggest to classify the Vaiont (respectively La Clapière) landslide as belonging to the velocity weakening unstable (respectively strengthening stable) sliding regime.

preprint2003arXiv

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

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 to a super-critical regime n>1. In the sub-critical regime, we recover and document the crossover from an Omori exponent 1-theta for t<t* to 1+theta for t<t* found previously in [Sornette and Sornette, 1999a] for a special case of the ETAS model. In the super-critical regime n>1 and theta>0, we find a novel transition from an Omori decay law with exponent 1-theta fot t<t* to an explosive exponential increase of the seismicity rate 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<tau, similar to the local law, to an exponential increase at large times. These results can rationalize many of the stylized facts reported for aftershock and foreshock sequences, such as (i) the suggestion that a small p-value may be a precursor of a large earthquake, (ii) the relative seismic quiescence sometimes observed before large aftershocks, (iii) the positive correlation between b and p-values, (iv) the observation that great earthquakes are sometimes preceded by a decrease of b-value and (v) the acceleration of the seismicity preceding great earthquakes.

preprint2002arXiv

Collective Origin of the Coexistence of Apparent RMT Noise and Factors in Large Sample Correlation Matrices

Through simple analytical calculations and numerical simulations, we demonstrate the generic existence of a self-organized macroscopic state in any large multivariate system possessing non-vanishing average correlations between a finite fraction of all pairs of elements. The coexistence of an eigenvalue spectrum predicted by random matrix theory (RMT) and a few very large eigenvalues in large empirical correlation matrices is shown to result from a bottom-up collective effect of the underlying time series rather than a top-down impact of factors. Our results, in excellent agreement with previous results obtained on large financial correlation matrices, show that there is relevant information also in the bulk of the eigenvalue spectrum and rationalize the presence of market factors previously introduced in an ad hoc manner.

preprint2002arXiv

Endogenous versus Exogenous Crashes in Financial Markets

We perform an extended analysis of the distribution of drawdowns in the two leading exchange markets (US dollar against the Deutsmark and against the Yen), in the major world stock markets, in the U.S. and Japanese bond market and in the gold market, by introducing the concept of ``coarse-grained drawdowns,'' which allows for a certain degree of fuzziness in the definition of cumulative losses and improves on the statistics of our previous results on the existence of ``outliers'' or ``kings.'' Then, for each identified outlier, we check whether log-periodic power law signatures (LPPS) are present and take the existence of LPPS as the qualifying signature for an endogenous crash: this is because a drawdown outlier is seen as the end of a speculative unsustainable accelerating bubble generated endogenously. In the absence of LPPS, we are able to identify what seems to have been the relevant historical event, i.e., a new piece of information of such magnitude and impact that it is seems reasonable to attribute the crash to it, in agreement with the standard view of the efficient market hypothesis. Such drawdown outliers are classified as having an exogenous origin. Globally over all the markets analyzed, we identify 49 outliers, of which 25 are classified as endogenous, 22 as exogeneous and 2 as associated with the Japanese anti-bubble. Restricting to the world market indices, we find 31 outliers, of which 19 are endogenous, 10 are exogenous and 2 are associated with the Japanese anti-bubble. The combination of the two proposed detection techniques, one for drawdown outliers and the second for LPPS, provides a novel and systematic taxonomy of crashes further subtantiating the importance of LPPS.

preprint2002arXiv

Fractal Plate Tectonics

We analyze in details the statistical significance of the claim by Bird [2002] of a power law distribution of plate areas covering the Earth and confirm that the power law with exponent 0.25 +- 0.05 is the most robust and parsimonious model for all plates, including the very largest plates, when taking into account the constraint that the plates areas must sum up to 4 pi steradians. We propose a general class of fragmentation models that rationalize this observation and discuss the implications for the earth dynamics and the general self-organization of tectonic deformations at multiple scales.

preprint2001arXiv

Positive Feedback, Memory and the Predictability of Earthquakes

We review the "critical point" concept for large earthquakes and enlarge it in the framework of so-called "finite-time singularities". The singular behavior associated with accelerated seismic release is shown to result from a positive feedback of the seismic activity on its release rate. The most important mechanisms for such positive feedback are presented. We introduce and solve analytically a novel simple model of geometrical positive feedback in which the stress shadow cast by the last large earthquake is progressively fragmented by the increasing tectonic stress. Finally, we present a somewhat speculative figure that tends to support a mechanism based on the decay of stress shadows. This figure suggests that a large earthquake in Southern California of size similar to the 1812 great event is maturing.

preprint2000arXiv

"Thermometers" of Speculative Frenzy

Establishing unambiguously the existence of speculative bubbles is an on-going controversy complicated by the need of defining a model of fundamental prices. Here, we present a novel empirical method which bypasses all the difficulties of the previous approaches by monitoring external indicators of an anomalously growing interest in the public at times of bubbles. From the definition of a bubble as a self-fulfilling reinforcing price change, we identify indicators of a possible self-reinforcing imitation between agents in the market. We show that during the build-up phase of a bubble, there is a growing interest in the public for the commodity in question, whether it consists in stocks, diamonds or coins. That interest can be estimated through different indicators: increase in the number of books published on the topic, increase in the subscriptions to specialized journals. Moreover, the well-known empirical rule according to which the volume of sales is growing during a bull market finds a natural interpretation in this framework: sales increases in fact reveal and pinpoint the progress of the bubble's diffusion throughout society. We also present a simple model of rational expectation which maps exactly onto the Ising model on a random graph. The indicators are then interpreted as ``thermometers'', measuring the balance between idiosyncratic information (noise temperature) and imitation (coupling) strength. In this context, bubbles are interpreted as low or critical temperature phases, where the imitation strength carries market prices up essentially independently of fundamentals. Contrary to the naive conception of a bubble and a crash as times of disorder, on the contrary, we show that bubbles and crashes are times where the concensus is too strong.

preprint2000arXiv

Acoustic fluidization for earthquakes?

Melosh [1996] has suggested that acoustic fluidization could provide an alternative to theories that are invoked as explanations for why some crustal faults appear to be weak. We show that there is a subtle but profound inconsistency in the theory that unfortunately invalidates the results. We propose possible remedies but must acknowledge that the relevance of acoustic fluidization remains an open question.

preprint2000arXiv

Artifactual log-periodicity in finite size data: Relevance for earthquake aftershocks

The recently proposed discrete scale invariance and its associated log-periodicity are an elaboration of the concept of scale invariance in which the system is scale invariant only under powers of specific values of the magnification factor. We report on the discovery of a novel mechanism for such log-periodicity relying solely on the manipulation of data. This ``synthetic'' scenario for log-periodicity relies on two steps: (1) the fact that approximately logarithmic sampling in time corresponds to uniform sampling in the logarithm of time; and (2) a low-pass-filtering step, as occurs in constructing cumulative functions, in maximum likelihood estimations, and in de-trending, reddens the noise and, in a finite sample, creates a maximum in the spectrum leading to a most probable frequency in the logarithm of time. We explore in detail this mechanism and present extensive numerical simulations. We use this insight to analyze the 27 best aftershock sequences studied by Kisslinger and Jones [1991] to search for traces of genuine log-periodic corrections to Omori's law, which states that the earthquake rate decays approximately as the inverse of the time since the last main shock. The observed log-periodicity is shown to almost entirely result from the ``synthetic scenario'' owing to the data analysis. From a statistical point of view, resolving the issue of the possible existence of log-periodicity in aftershocks will be very difficult as Omori's law describes a point process with a uniform sampling in the logarithm of the time. By construction, strong log-periodic fluctuations are thus created by this logarithmic sampling.

preprint2000arXiv

Functional Renormalization Prediction of Rupture

We develop theoretical formulas for the prediction of the rupture of systems which are known to exhibit a critical behavior, based solely on the knowledge of the early time evolution of an observable, such as the acoustic emission rate as a function of time or of stress. From the parameterization of such early time evolution in terms of a low-order polynomial, we use the functional renormalization approach introduced by Yukalov and Gluzman to transform this polynomial into a function which is asymptotically a power law. The value of the critical time tc,conditioned on the assumption that tc exists, is thus determined from the knowledge of the coefficients of the polynomials. We test with success this prediction scheme with respect to the order of the polynomials and as a function of noise.

preprint2000arXiv

Re-examination of log-periodicity observed in the seismic precursors of the 1989 Loma Prieta earthquake

Based on several empirical evidence, a series of papers has advocated the concept that seismicity prior to a large earthquake can be understood in terms of the statistical physics of a critical phase transition. In this model, the cumulative Benioff strain (BS) increases as a power-law time-to-failure before the final event. This power law reflects a kind of scale invariance with respect to the distance to the critical point. A few years ago, on the basis of a fit of the cumulative BS released prior to the 1989 Loma Prieta earthquake, Sornette and Sammis [1995] proposed that this scale invariance could be partially broken into a discrete scale invariance (DSI). The observable consequence of DSI takes the form of log-periodic oscillations decorating the accelerating power law. They found that the quality of the fit and the predicted time of the event are significantly improved by the introduction of log-periodicity. Here, we present a battery of synthetic tests performed to quantify the statistical significance of this claim. We find that log-periodic oscillations with frequency and regularity similar to those of the Loma Prieta case are very likely to be generated by the interplay of the low pass filtering step due to the construction of cumulative functions together with the approximate power law acceleration. Thus, the single Loma Prieta case alone cannot support the initial claim and additional cases and further study are needed to increase the signal-to-noise ratio if any. The present study will be a useful methodological benchmark for future testing of additional events when the methodology and data to construct reliable Benioff strain function become available.

preprint2000arXiv

Statistical Asynchronous Regression: Determining the Relationship Between two Quantities that are not Measured Simultaneously

We introduce the Statistical Asynchronous Regression (SAR) method: a technique for determining a relationship between two time varying quantities without simultaneous measurements of both quantities. We require that there is a time invariant, monotonic function Y = u(X) relating the two quantities, Y and X. In order to determine u(X), we only need to know the statistical distributions of X and Y. We show that u(X) is the change of variables that converts the distribution of X into the distribution of Y, while conserving probability. We describe an algorithm for implementing this method and apply it to several example distributions. We also demonstrate how the method can separate spatial and temporal variations from a time series of energetic electron flux measurements made by a spacecraft in geosynchronous orbit. We expect this method will be useful to the general problem of spacecraft instrument calibration. We also suggest some applications of the SAR method outside of space physics.

preprint1999arXiv

Renormalization of earthquake aftershocks

Together with the Gutenberg-Richter distribution of earthquake magnitudes, Omori's law is the best established empirical characterization of earthquake sequences and states that the number of smaller earthquakes per unit time triggered by a main shock decays approximately as the inverse of the time ($1/t^p$, with $p \approx 1$) since the main shock. Based on these observations, we explore the theoretical hypothesis in which each earthquake can produce a series of aftershock independently of its size according to its ``local'' Omori's law with exponent $p=1+θ$. In this scenario, an aftershock of the main shock produces itself other aftershocks which themselves produce aftershocks, and so on. The global observable Omori's law is found to have two distinct power law regimes, the first one with exponent $p_-=1 - θ$ for time $t < t^* \sim κ^{-1/θ}$, where $0<1-κ<1$ measures the fraction of triggered earthquakes per triggering earthquake, and the second one with exponent $p_+=1 + θ$ for larger times. The existence of these two regimes rationalizes the observation of Kisslinger and Jones [1991] that the Omori's exponent $p$ seems positively correlated to the surface heat flow: a higher heat flow is a signature of a higher crustal temperature, which leads to larger strain relaxation by creep, corresponding to fewer events triggered per earthquake, i.e. to a larger $κ$, and thus to a smaller $t^*$, leading to an effective measured exponent more heavily weighted toward $p_+>1$.

preprint1998arXiv

``String'' formulation of the Dynamics of the Forward Interest Rate Curve

We propose a formulation of the term structure of interest rates in which the forward curve is seen as the deformation of a string. We derive the general condition that the partial differential equations governing the motion of such string must obey in order to account for the condition of absence of arbitrage opportunities. This condition takes a form similar to a fluctuation-dissipation theorem, albeit on the same quantity (the forward rate), linking the bias to the covariance of variation fluctuations. We provide the general structure of the models that obey this constraint in the framework of stochastic partial (possibly non-linear) differential equations. We derive the general solution for the pricing and hedging of interest rate derivatives within this framework, albeit for the linear case (we also provide in the appendix a simple and intuitive derivation of the standard European option problem). We also show how the ``string'' formulation simplifies into a standard N-factor model under a Galerkin approximation.

preprint1998arXiv

An Observational Test of the Critical Earthquake Concept

We test the concept that seismicity prior to a large earthquake can be understood in terms of the statistical physics of a critical phase transition. In this model, the cumulative seismic strain release increases as a power-law time-to-failure before the final event. Furthermore, the region of correlated seismicity predicted by this model is much greater than would be predicted from simple elasto-dynamic interactions. We present a systematic procedure to test for the accelerating seismicity predicted by the critical point model and to identify the region approaching criticality, based on a comparison between the observed cumulative energy (Benioff strain) release and the power-law behavior predicted by theory. This method is used to find the critical region before all earthquakes along the San Andreas system since 1950 with M 6.5. The statistical significance of our results is assessed by performing the same procedure on a large number of randomly generated synthetic catalogs. The null hypothesis, that the observed acceleration in all these earthquakes could result from spurious patterns generated by our procedure in purely random catalogs, is rejected with 99.5% confidence. An empirical relation between the logarithm of the critical region radius (R) and the magnitude of the final event (M) is found, such that log R μ0.5 M, suggesting that the largest probable event in a given region scales with the size of the regional fault network.

preprint1998arXiv

Data-Adaptive Wavelets and Multi-Scale Singular Spectrum Analysis

Using multi-scale ideas from wavelet analysis, we extend singular-spectrum analysis (SSA) to the study of nonstationary time series of length $N$ whose intermittency can give rise to the divergence of their variance. SSA relies on the construction of the lag-covariance matrix C on M lagged copies of the time series over a fixed window width W to detect the regular part of the variability in that window in terms of the minimal number of oscillatory components; here W = M Dt, with Dt the time step. The proposed multi-scale SSA is a local SSA analysis within a moving window of width M <= W <= N. Multi-scale SSA varies W, while keeping a fixed W/M ratio, and uses the eigenvectors of the corresponding lag-covariance matrix C_M as a data-adaptive wavelets; successive eigenvectors of C_M correspond approximately to successive derivatives of the first mother wavelet in standard wavelet analysis. Multi-scale SSA thus solves objectively the delicate problem of optimizing the analyzing wavelet in the time-frequency domain, by a suitable localization of the signal's covariance matrix. We present several examples of application to synthetic signals with fractal or power-law behavior which mimic selected features of certain climatic and geophysical time series. A real application is to the Southern Oscillation index (SOI) monthly values for 1933-1996. Our methodology highlights an abrupt periodicity shift in the SOI near 1960. This abrupt shift between 4 and 3 years supports the Devil's staircase scenario for the El Nino/Southern Oscillation phenomenon.

preprint1998arXiv

Evidence of discrete scale invariance in DLA and time-to-failure by canonical averaging

Discrete scale invariance, which corresponds to a partial breaking of the scaling symmetry, is reflected in the existence of a hierarchy of characteristic scales l0, c l0, c^2 l0,... where c is a preferred scaling ratio and l0 a microscopic cut-off. Signatures of discrete scale invariance have recently been found in a variety of systems ranging from rupture, earthquakes, Laplacian growth phenomena, ``animals'' in percolation to financial market crashes. We believe it to be a quite general, albeit subtle phenomenon. Indeed, the practical problem in uncovering an underlying discrete scale invariance is that standard ensemble averaging procedures destroy it as if it was pure noise. This is due to the fact, that while c only depends on the underlying physics, l0 on the contrary is realisation-dependent. Here, we adapt and implement a novel so-called ``canonical'' averaging scheme which re-sets the l0 of different realizations to approximately the same value. The method is based on the determination of a realization-dependent effective critical point obtained from, e.g., a maximum susceptibility criterion. We demonstrate the method on diffusion limited aggregation and a model of rupture.

preprint1998arXiv

Solid friction at high sliding velocities: an explicit 3D dynamical SPH approach

We present realistic 3D numerical simulations of elastic bodies sliding on top of each other in a regime of velocities ranging from meters to tens of meters per second using the so-called Smoothed Particle Hydrodynamics (SPH) method. Our investigations are restricted to regimes of pressure and roughness where only elastic deformations occur between asperities at the contact surface between the slider block and the substrate. In this regime, solid friction is due to the generation of vibrational radiations which are subsequently damped out. We study periodic commensurate and incommensurate asperities and various types of disordered surfaces. We report the evidence of a transition from zero (or non-measurable $μ< 0.001$) friction to a finite friction as the normal pressure increases above about $10^6~Pa$. For larger normal pressures (up to $10^9~Pa$), we find a remarkably universal value for the friction coefficient $μ\approx 0.06$, which is independent of the internal dissipation strength over three order of magnitudes, and independent of the detailled nature of the slider block-substrate interactions. We find that disorder may either decrease or increase $μ$ due to the competition between two effects: disorder detunes the coherent vibrations of the asperties that occur in the periodic case, leading to weaker acoustic radiation and thus weaker damping. On the other hand, large disorder leads to stronger vibration amplitudes at local asperities and thus stronger damping. Our simulations have confirmed the existence of jumps over steps or asperities of the slider blocks occurring at the largest velocities studied ($10~m/s$). These jumps lead to chaotic motions similar to the bouncing-ball problem. We find a velocity strengthening with a doubling of the friction coefficient as the velocity increases from $1~m/s$ to $10~m/s$.

preprint1998arXiv

Testing self-organized criticality by induced seismicity

We examine the hypothesis proposed in recent years by several authors that the crust is in a self-organized critical (SOC) state by exploring how the SOC concept can help in understanding the observed earthquake clustering on relatively narrow fault domains and the phenomenon of induced seismicity. We review the major reported cases of induced seismicity in various parts of the world and find that both pore pressure changes (+/-p) and mass transfers (+/-m) leading to incremental deviatoric stresses of <<1 MPa are sufficient to trigger seismic instabilities in the uppermost crust with magnitude ranging up to 7.0 in otherwise historically aseismic areas. Once triggered, stress variations of at least 1 order of magnitude less but still larger than the ~0.01 MPa tidal stress are enough to sustain seismic activity. We argue that these observations are in accord with the SOC hypothesis as they show that a significant fraction of the crust is not far from instability and can thus be made unstable by minute perturbations. Not all perturbations, however, trigger seismic activity; this is also compatible with the SOC hypothesis which embodies naturally the existence of large heterogeneities in the stress field. In conclusion, the properties of induced seismicity and their rationalization in terms of the SOC concept provide further evidence that potential seismic hazards extend over a much larger area than that where earthquakes are frequent.

preprint1997arXiv

Linear stochastic dynamics with nonlinear fractal properties

Stochastic processes with multiplicative noise have been studied independently in several different contexts over the past decades. We focus on the regime, found for a generic set of control parameters, in which stochastic processes with multiplicative noise produce intermittency of a special kind, characterized by a power law probability density distribution. We present a review of applications on population dynamics, epidemics, finance and insurance applications with relation to ARCH(1) process, immigration and investment portfolios and the internet. We highlight the common physical mechanism and summarize the main known results. The distribution and statistical properties of the duration of intermittent bursts are also characterized in details.