Q factor: A measure of competition between the topper and the average in percolation and in self-organized criticality.

We define the Q factor in the percolation problem as the quotient of the size of the largest cluster and the average size of all clusters. As the occupation probability p is increased, the Q factor for the system size L grows systematically to its maximum value Q_{max}(L) at a specific value p_{max}(L) and then gradually decays. Our numerical study of site percolation problems on the square, triangular, and simple cubic lattices exhibits that the asymptotic values of p_{max}, though close, are distinct from the corresponding percolation thresholds of these lattices. We also show, using scaling analysis, that at p_{max} the value of Q_{max}(L) diverges as L^{d} (d denoting the dimension of the lattice) as the system size approaches its asymptotic limit. We further extend this idea to nonequilibrium systems such as the sandpile model of self-organized criticality. Here the Q(ρ,L) factor is the quotient of the size of the largest avalanche and the cumulative average of the sizes of all the avalanches, with ρ the drop density of the driving mechanism. This study was prompted by some observations in sociophysics.

Finding critical points and correlation length exponents using finite size scaling of Gini index.

The order parameter for a continuous transition shows diverging fluctuation near the critical point. Here we show, through numerical simulations and scaling arguments, that the inequality (or variability) between the values of an order parameter, measured near a critical point, is independent of the system size. Quantification of such variability through the Gini index (g) therefore leads to a scaling form g=G[|F-F_{c}|N^{1/dν}], where F denotes the driving parameter for the transition (e.g., temperature T for ferromagnetic to paramagnetic transition, or lattice occupation probability p in percolation), N is the system size, d is the spatial dimension and ν is the correlation length exponent. We demonstrate the scaling for the Ising model in two and three dimensions, site percolation on square lattice, and the fiber bundle model of fracture.

Inequality of avalanche sizes in models of fracture.

Prediction of an imminent catastrophic event in a driven disordered system is of paramount importance-from the laboratory scale controlled fracture experiment to the largest scale of mechanical failure, i.e., earthquakes. It has long been conjectured that the statistical regularities in the energy emission time series mirror the "health" of such driven systems and hence have the potential for forecasting imminent catastrophe. Among other statistical regularities, a measure of how unequal avalanche sizes are is potentially a crucial indicator of imminent failure. The inequalities of avalanche sizes are quantified using inequality indices traditionally used in socioeconomic systems: the Gini index g, the Hirsch index h, and the Kolkata index k. It is shown analytically (for the mean-field case) and numerically (for the non-mean-field case) with models of quasi-brittle materials that the indices show universal behavior near the breaking points in such models and hence could serve as indicators of imminent breakdown of stressed disordered systems.

Kinetic Models of Wealth Distribution with Extreme Inequality: Numerical Study of Their Stability against Random Exchanges.

In view of some recent reports on global wealth inequality, where a small number (often a handful) of people own more wealth than 50% of the world's population, we explored if kinetic exchange models of markets could ever capture features where a significant fraction of wealth can concentrate in the hands of a few as the market size N approaches infinity. One existing example of such a kinetic exchange model is the Chakraborti or Yard-Sale model; in the absence of tax redistribution, etc., all wealth ultimately condenses into the hands of a single individual (for any value of N), and the market dynamics stop. With tax redistribution, etc., steady-state dynamics are shown to have remarkable applicability in many cases in our extremely unequal world. We show that another kinetic exchange model (called the Banerjee model) has intriguing intrinsic dynamics, where only ten rich traders or agents possess about 99.98% of the total wealth in the steady state (without any tax, etc., like external manipulation) for any large N value. We will discuss the statistical features of this model using Monte Carlo simulations. We will also demonstrate that if each trader has a non-zero probability f of engaging in random exchanges, then these condensations of wealth (e.g., 100% in the hand of one agent in the Chakraborti model, or about 99.98% in the hands of ten agents in the Banerjee model) disappear in the large N limit. Moreover, due to the built-in possibility of random exchange dynamics in the earlier proposed Goswami-Sen model, where the exchange probability decreases with the inverse power of the wealth difference between trading pairs, one does not see any wealth condensation phenomena. In this paper, we explore these aspects of statistics of these intriguing models.

Sandpile Universality in Social Inequality: Gini and Kolkata Measures.

Social inequalities are ubiquitous and evolve towards a universal limit. Herein, we extensively review the values of inequality measures, namely the Gini (g) index and the Kolkata (k) index, two standard measures of inequality used in the analysis of various social sectors through data analysis. The Kolkata index, denoted as k, indicates the proportion of the 'wealth' owned by (1-k) fraction of the 'people'. Our findings suggest that both the Gini index and the Kolkata index tend to converge to similar values (around g=k≈0.87, starting from the point of perfect equality, where g=0 and k=0.5) as competition increases in different social institutions, such as markets, movies, elections, universities, prize winning, battle fields, sports (Olympics), etc., under conditions of unrestricted competition (no social welfare or support mechanism). In this review, we present the concept of a generalized form of Pareto's 80/20 law (k=0.80), where the coincidence of inequality indices is observed. The observation of this coincidence is consistent with the precursor values of the g and k indices for the self-organized critical (SOC) state in self-tuned physical systems such as sand piles. These results provide quantitative support for the view that interacting socioeconomic systems can be understood within the framework of SOC, which has been hypothesized for many years. These findings suggest that the SOC model can be extended to capture the dynamics of complex socioeconomic systems and help us better understand their behavior.

Quantum annealing: an overview.

In this review, after providing the basic physical concept behind quantum annealing (or adiabatic quantum computation), we present an overview of some recent theoretical as well as experimental developments pointing to the issues which are still debated. With a brief discussion on the fundamental ideas of continuous and discontinuous quantum phase transitions, we discuss the Kibble-Zurek scaling of defect generation following a ramping of a quantum many body system across a quantum critical point. In the process, we discuss associated models, both pure and disordered, and shed light on implementations and some recent applications of the quantum annealing protocols. Furthermore, we discuss the effect of environmental coupling on quantum annealing. Some possible ways to speed up the annealing protocol in closed systems are elaborated upon: we especially focus on the recipes to avoid discontinuous quantum phase transitions occurring in some models where energy gaps vanish exponentially with the system size. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Quantum annealing and computation: challenges and perspectives.

In the introductory article of this theme issue, we provide an overview of quantum annealing and computation with a very brief summary of the individual contributions to this issue made by experts as well as a few young researchers. We hope the readers will get the touch of the excitement as well as the perspectives in this unusually active field and important developments there. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Stochastic Learning in Kolkata Paise Restaurant Problem: Classical and Quantum Strategies.

We review the results for stochastic learning strategies, both classical (one-shot and iterative) and quantum (one-shot only), for optimizing the available many-choice resources among a large number of competing agents, developed over the last decade in the context of the Kolkata Paise Restaurant (KPR) Problem. Apart from few rigorous and approximate analytical results, both for classical and quantum strategies, most of the interesting results on the phase transition behavior (obtained so far for the classical model) uses classical Monte Carlo simulations. All these including the applications to computer science [job or resource allotments in Internet-of-Things (IoT)], transport engineering (online vehicle hire problems), operation research (optimizing efforts for delegated search problem, efficient solution of Traveling Salesman problem) will be discussed.

Kinetic exchange income distribution models with saving propensities: inequality indices and self-organized poverty level.

We report the numerical results for the steady-state income or wealth distribution [Formula: see text] and the resulting inequality measures (Gini [Formula: see text] and Kolkata [Formula: see text] indices) in the kinetic exchange models of market dynamics. We study the variations of [Formula: see text] and of the indices [Formula: see text] and [Formula: see text] with the saving propensity [Formula: see text] of the agents, with two different kinds of trade (kinetic exchange) dynamics. In the first case, the exchange occurs between randomly chosen pairs of agents and in the next, one of the agents in the chosen pair is the poorest of all and the other agent is randomly picked up from the rest of the population (where, in the steady state, a self-organized poverty level or SOPL appears). These studies have also been made for two different kinds of saving behaviours. One, where each agent has the same value of [Formula: see text] (constant over time) and the other where [Formula: see text] for each agent can take two values (0 and 1), changing randomly over a fraction of time [Formula: see text] of choosing [Formula: see text]. We find that the inequality decreases with increasing savings ([Formula: see text]); inequality indices ([Formula: see text] and [Formula: see text]) decrease and SOPL increases with increasing [Formula: see text], indicating possible applications in economic policy making. This article is part of the theme issue 'Kinetic exchange models of societies and economies'.

Social inequality analysis of fiber bundle model statistics and prediction of materials failure.

Inequalities are abundant in a society with a number of agents competing for a limited amount of resources. Statistics on such social inequalities are usually represented by the Lorenz function L(p), where p fraction of the population possesses L(p) fraction of the total wealth, when the population is arranged in ascending order of their wealth. Similarly, in scientometrics, such inequalities can be represented by a plot of the citation count versus the respective number of papers by a scientist, again arranged in ascending order of their citation counts. Quantitatively, these inequalities are captured by the corresponding inequality indices, namely, the Kolkata k and the Hirsch h indices, given by the fixed points of these nonlinear (Lorenz and citation) functions. In statistical physics of criticality, the fixed points of the renormalization group generator functions are studied in their self-similar limit, where their (fractal) structure converges to a unique form (macroscopic in size and lone). The statistical indices in social science, however, correspond to the fixed points where the values of the generator function (wealth or citation sizes) are commensurately abundant in fractions or numbers (of persons or papers). It has already been shown that under extreme competitions in markets or at universities, the k index approaches a universal limiting value, as the dynamics of competition progresses. We introduce and study these indices for the inequalities of (prefailure) avalanches, given by their nonlinear size distributions in fiber bundle models of nonbrittle materials. We show how prior knowledge of the terminal and (almost) universal value of the k index for a wide range of disorder parameters can help in predicting an imminent catastrophic breakdown in the model. This observation is also complemented by noting a similar (but not identical) behavior of the Hirsch index (h), redefined for such avalanche statistics.

Development of Econophysics: A Biased Account and Perspective from Kolkata.

We present here a somewhat personalized account of the emergence of econophysics as an attractive research topic in physical, as well as social, sciences. After a rather detailed storytelling about our endeavors from Kolkata, we give a brief description of the main research achievements in a simple and non-technical language. We also briefly present, in technical language, a piece of our recent research result. We conclude our paper with a brief perspective.

Phase transition in the Kolkata Paise Restaurant problem.

A novel phase transition behavior is observed in the Kolkata Paise Restaurant problem where a large number (N) of agents or customers collectively (and iteratively) learn to choose among the N restaurants where she would expect to be alone that evening and would get the only dish available there (or may get randomly picked up if more than one agent arrive there that evening). The players are expected to evolve their strategy such that the publicly available information about past crowds in different restaurants can be utilized and each of them is able to make the best minority choice. For equally ranked restaurants, we follow two crowd-avoiding strategies: strategy I, where each of the ni(t) number of agents arriving at the ith restaurant on the tth evening goes back to the same restaurant the next evening with probability [ni(t)]-α, and strategy II, with probability p, when ni(t)>1. We study the steady state (t-independent) utilization fraction f:(1-f) giving the steady state (wastage) fraction of restaurants going without any customer at any particular evening. With both strategies, we find, near αc=0+ (in strategy I) or p=1- (in strategy II), the steady state wastage fraction (1-f)∝(α-αc)ß or (pc-p)ß with ß≃0.8,0.87,1.0, and the convergence time τ [for f(t) becoming independent of t] varies as τ∝(α-αc)-γ or (pc-p)-γ, with γ≃1.18,1.11,1.05 in infinite-dimensions (rest of the N-1 neighboring restaurants), three dimensions (six neighbors), and two dimensions (four neighbors), respectively.

Flory-like statistics of fracture in the fiber bundle model as obtained via Kolmogorov dispersion for turbulence: A conjecture.

It has long been conjectured that (rapid) fracture propagation dynamics in materials and turbulent motion of fluids are two manifestations of the same physical process. The universality class of turbulence (Kolmogorov dispersion, in particular) is conjectured to be identifiable with the Flory statistics for linear polymers (self-avoiding walks on lattices). These help us to relate fracture statistics to those of linear polymers (Flory statistics). The statistics of fracture in the fiber bundle model (FBM) are now well studied and many exact results are now available for the equal-load-sharing (ELS) scheme. Yet, the correlation length exponent in this model was missing and we show here how the correspondence between fracture statistics and the Flory mapping of Kolmogorov statistics for turbulence helps us to make a conjecture about the value of the correlation length exponent for fracture in the ELS limit of FBM and, also, about the upper critical dimension. In addition, the fracture avalanche size exponent values at lower dimensions (as estimated from such mapping to Flory statistics) also compare well with the observations.

Hydrodynamic descriptions for surface roughness in fracture front propagation.

Fracture is ubiquitous in a crystalline material. Inspired by the observed phenomenological similarities between the spatial profile of a fractured surface and velocities in hydrodynamic turbulence, we set up a hydrodynamic description for the dynamics of fracture surface propagation mode I or opening fracture front. We consider several related continuum hydrodynamic models and use them to extract the similarities between the profile of a fractured surface and velocities in hydrodynamic turbulence. We conclude that a fractured surface should be generically self-similar with an underlying multifractal behaviour.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Statistical physics of fracture and earthquakes.

Manifestations of emergent properties in stressed disordered materials are often the result of an interplay between strong perturbations in the stress field around defects. The collective response of a long-ranged correlated multi-component system is an ideal playing field for statistical physics. Hence, many aspects of such collective responses in widely spread length and energy scales can be addressed by the tools of statistical physics. In this theme issue, some of these aspects are treated from various angles of experiments, simulations and analytical methods, and connected together by their common base of complex-system dynamics.This article is part of the theme issue 'Statistical physics of fracture and earthquakes' .

Possible ergodic-nonergodic regions in the quantum Sherrington-Kirkpatrick spin glass model and quantum annealing.

We explore the behavior of the order parameter distribution of the quantum Sherrington-Kirkpatrick model in the spin glass phase using Monte Carlo technique for the effective Suzuki-Trotter Hamiltonian at finite temperatures and that at zero temperature obtained using the exact diagonalization method. Our numerical results indicate the existence of a low- but finite-temperature quantum-fluctuation-dominated ergodic region along with the classical fluctuation-dominated high-temperature nonergodic region in the spin glass phase of the model. In the ergodic region, the order parameter distribution gets narrower around the most probable value of the order parameter as the system size increases. In the other region, the Parisi order distribution function has nonvanishing value everywhere in the thermodynamic limit, indicating nonergodicity. We also show that the average annealing time for convergence (to a low-energy level of the model, within a small error range) becomes system size independent for annealing down through the (quantum-fluctuation-dominated) ergodic region. It becomes strongly system size dependent for annealing through the nonergodic region. Possible finite-size scaling-type behavior for the extent of the ergodic region is also addressed.

Universality of Citation Distributions for Academic Institutions and Journals.

Citations measure the importance of a publication, and may serve as a proxy for its popularity and quality of its contents. Here we study the distributions of citations to publications from individual academic institutions for a single year. The average number of citations have large variations between different institutions across the world, but the probability distributions of citations for individual institutions can be rescaled to a common form by scaling the citations by the average number of citations for that institution. We find this feature seems to be universal for a broad selection of institutions irrespective of the average number of citations per article. A similar analysis for citations to publications in a particular journal in a single year reveals similar results. We find high absolute inequality for both these sets, Gini coefficients being around 0.66 and 0.58 for institutions and journals respectively. We also find that the top 25% of the articles hold about 75% of the total citations for institutions and the top 29% of the articles hold about 71% of the total citations for journals.

Classical-to-quantum crossover in the critical behavior of the transverse-field Sherrington-Kirkpatrick spin glass model.

We study the critical behavior of the Sherrington-Kirkpatrick model in transverse field (at finite temperature) using Monte Carlo simulation and exact diagonalization (at zero temperature). We determine the phase diagram of the model by estimating the Binder cumulant. We also determine the correlation length exponent from the collapse of the scaled data. Our numerical studies here indicate that critical Binder cumulant (indicating the universality class of the transition behavior) and the correlation length exponent cross over from their "classical" to "quantum" values at a finite temperature (unlike the cases of pure systems, where such crossovers occur at zero temperature). We propose a qualitative argument supporting such an observation, employing a simple tunneling picture.

Zipf's law in city size from a resource utilization model.

We study a resource utilization scenario characterized by intrinsic fitness. To describe the growth and organization of different cities, we consider a model for resource utilization where many restaurants compete, as in a game, to attract customers using an iterative learning process. Results for the case of restaurants with uniform fitness are reported. When fitness is uniformly distributed, it gives rise to a Zipf law for the number of customers. We perform an exact calculation for the utilization fraction for the case when choices are made independent of fitness. A variant of the model is also introduced where the fitness can be treated as an ability to stay in the business. When a restaurant loses customers, its fitness is replaced by a random fitness. The steady state fitness distribution is characterized by a power law, while the distribution of the number of customers still follows the Zipf law, implying the robustness of the model. Our model serves as a paradigm for the emergence of Zipf law in city size distribution.

Self-organized dynamics in local load-sharing fiber bundle models.

We study the dynamics of a local load-sharing fiber bundle model in two dimensions under an external load (which increases with time at a fixed slow rate) applied at a single point. Due to the local load-sharing nature, the redistributed load remains localized along the boundary of the broken patch. The system then goes to a self-organized state with a stationary average value of load per fiber along the (increasing) boundary of the broken patch (damaged region) and a scale-free distribution of avalanche sizes and other related quantities are observed. In particular, when the load redistribution is only among nearest surviving fiber(s), the numerical estimates of the exponent values are comparable with those of the Manna model. When the load redistribution is uniform along the patch boundary, the model shows a simple mean-field limit of this self-organizing critical behavior, for which we give analytical estimates of the saturation load per fiber values and avalanche size distribution exponent. These are in good agreement with numerical simulation results.