Selection of noise level in strategy adoption for spatial social dilemmas

We studied spatial Prisoner's Dilemma and Stag Hunt games where both the strategy distribution and the players' individual noise level could evolve to reach higher individual payoff. Players are located on the sites of different two-dimensional lattices and gain their payoff from games with their neighbors by choosing unconditional cooperation or defection. The way of strategy adoption can be characterized by a single $K$ (temperature-like) parameter describing how strongly adoptions depend on the payoff-difference. If we start the system from a random strategy distribution with many different player specific $K$ parameters, the simultaneous evolution of strategies and $K$ parameters drives the system to a final stationary state where only one $K$ value remains. In the coexistence phase of cooperator and defector strategies the surviving $K$ parameter is in good agreement with the noise level that ensures the highest cooperation level if uniform $K$ is supposed for all players. In this paper we give a thorough overview about the properties of this evolutionary process.Comment: 10 two-column pages, 10 figures; accepted for publication in Physical Review

Experimental analysis of lateral impact on planar brittle material

The fragmentation of alumina and glass plates due to lateral impact is studied. A few hundred plates have been fragmented at different impact velocities and the produced fragments are analyzed. The method employed in this work allows one to investigate some geometrical properties of the fragments, besides the traditional size distribution usually studied in former experiments. We found that, although both materials exhibit qualitative similar fragment size distribution function, their geometrical properties appear to be quite different. A schematic model for two-dimensional fragmentation is also presented and its predictions are compared to our experimental results. The comparison suggests that the analysis of the fragments' geometrical properties constitutes a more stringent test of the theoretical models' assumptions than the size distribution

New methods for unmixing sediment grain size data

Grain size distribution (GSD) data are widely used in Earth sciences and although large data sets are regularly generated, detailed numerical analyses are not routine. Unmixing GSDs into components can help understand sediment provenance and depositional regimes/processes. End-member analysis (EMA), which fits one set of end-members to a given data set, is a powerful way to unmix GSDs into geologically meaningful parts. EMA estimates end-members based on covariability within a data set and can be considered as a nonparametric approach. Available EMA algorithms, however, either produce suboptimal solutions or are time consuming. We introduce unmixing algorithms inspired by hyperspectral image analysis that can be applied to GSD data and which provide an improvement over current techniques. Nonparametric EMA is often unable to identify unimodal grain size subpopulations that correspond to single sediment sources. An alternative approach is single-specimen unmixing (SSU), which unmixes individual GSDs into unimodal parametric distributions (e.g., lognormal). We demonstrate that the inherent nonuniqueness of SSU solutions renders this approach unviable for estimating underlying mixing processes. To overcome this, we develop a new algorithm to perform parametric EMA, whereby an entire data set can be unmixed into unimodal parametric end-members (e.g., Weibull distributions). This makes it easier to identify individual grain size subpopulations in highly mixed data sets. To aid investigators in applying these methods, all of the new algorithms are available in AnalySize, which is GUI software for processing and unmixing grain size data

Statistical mechanics of systems with heterogeneous agents: Minority Games

We study analytically a simple game theoretical model of heterogeneous interacting agents. We show that the stationary state of the system is described by the ground state of a disordered spin model which is exactly solvable within the simple replica symmetric ansatz. Such a stationary state differs from the Nash equilibrium where each agent maximizes her own utility. The latter turns out to be characterized by a replica symmetry broken structure. Numerical results fully agree with our analytic findings.Comment: 4 pages, 1 Postscript figure. Revised versio

Stochastic gain in population dynamics

We introduce an extension of the usual replicator dynamics to adaptive learning rates. We show that a population with a dynamic learning rate can gain an increased average payoff in transient phases and can also exploit external noise, leading the system away from the Nash equilibrium, in a reasonance-like fashion. The payoff versus noise curve resembles the signal to noise ratio curve in stochastic resonance. Seen in this broad context, we introduce another mechanism that exploits fluctuations in order to improve properties of the system. Such a mechanism could be of particular interest in economic systems.Comment: accepted for publication in Phys. Rev. Let

Stretched exponential behavior and random walks on diluted hypercubic lattices

Diffusion on a diluted hypercube has been proposed as a model for glassy relaxation and is an example of the more general class of stochastic processes on graphs. In this article we determine numerically through large scale simulations the eigenvalue spectra for this stochastic process and calculate explicitly the time evolution for the autocorrelation function and for the return probability, all at criticality, with hypercube dimensions $N$ up to N=28. We show that at long times both relaxation functions can be described by stretched exponentials with exponent 1/3 and a characteristic relaxation time which grows exponentially with dimension $N$. The numerical eigenvalue spectra are consistent with analytic predictions for a generic sparse network model.Comment: 16 pages, 7 figure

Characterization of the stretched exponential trap-time distributions in one-dimensional coupled map lattices

Stretched exponential distributions and relaxation responses are encountered in a wide range of physical systems such as glasses, polymers and spin glasses. As found recently, this type of behavior occurs also for the distribution function of certain trap time in a number of coupled dynamical systems. We analyze a one-dimensional mathematical model of coupled chaotic oscillators which reproduces an experimental set-up of coupled diode-resonators and identify the necessary ingredients for stretched exponential distributions.Comment: 8 pages, 8 figure

Assessment of multireference approaches to explicitly correlated full configuration interaction quantum Monte Carlo.

The Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method has proved able to provide near-exact solutions to the electronic Schrödinger equation within a finite orbital basis set, without relying on an expansion about a reference state. However, a drawback to the approach is that being based on an expansion of Slater determinants, the FCIQMC method suffers from a basis set incompleteness error that decays very slowly with the size of the employed single particle basis. The FCIQMC results obtained in a small basis set can be improved significantly with explicitly correlated techniques. Here, we present a study that assesses and compares two contrasting "universal" explicitly correlated approaches that fit into the FCIQMC framework: the [2]R12 method of Kong and Valeev [J. Chem. Phys. 135, 214105 (2011)] and the explicitly correlated canonical transcorrelation approach of Yanai and Shiozaki [J. Chem. Phys. 136, 084107 (2012)]. The former is an a posteriori internally contracted perturbative approach, while the latter transforms the Hamiltonian prior to the FCIQMC simulation. These comparisons are made across the 55 molecules of the G1 standard set. We found that both methods consistently reduce the basis set incompleteness, for accurate atomization energies in small basis sets, reducing the error from 28 mEh to 3-4 mEh. While many of the conclusions hold in general for any combination of multireference approaches with these methodologies, we also consider FCIQMC-specific advantages of each approach.Royal Societ