Real classical shadows

Biological systems are remarkably susceptible to relatively small temperature changes. The most obvious example is fever, when a modest rise in body temperature of only few Kelvin has strong effects on our immune system and how it fights pathogens. Another very important example is climate change, when even smaller temperature changes lead to dramatic shifts in ecosystems. Although it is generally accepted that the main effect of an increase in temperature is the acceleration of biochemical reactions according to the Arrhenius equation, it is not clear how it affects large biochemical networks with complicated architectures. For developmental systems such as fly and frog, it has been shown that the system response to temperature deviates in a characteristic manner from the linear Arrhenius plot of single reactions, but a rigorous explanation has not been given yet. Here we use a graph-theoretical interpretation of the mean first-passage times of a biochemical master equation to give a statistical description. We find that in the limit of large system size and if the network has a bias towards a target state, then the Arrhenius plot is generically quadratic, in excellent agreement with numerical simulations for large networks as well as with experimental data for developmental times in fly. We also discuss under which conditions this generic response can be violated, for example for linear chains, which have only one spanning tree.

The search for a unified formulation for describing various non-equilibrium processes is a central task of modern non-equilibrium thermodynamics. In this paper, a novel steady-state thermodynamic formalism was established for general Markov processes described by the Chapman-Kolmogorov equation. Furthermore, corresponding formalisms of steady-state thermodynamics for the master equation and Fokker-Planck equation could be rigorously derived in mathematics. To be concrete, we proved that (1) in the limit of continuous time, the steady-state thermodynamic formalism for the Chapman-Kolmogorov equation fully agrees with that for the master equation; (2) a similar one-to-one correspondence could be established rigorously between the master equation and Fokker-Planck equation in the limit of large system size; (3) when a Markov process is restrained to one-step jump, the steady-state thermodynamic formalism for the Fokker-Planck equation with discrete state variables also goes to that for master equations, as the discretization step gets smaller and smaller. Our analysis indicated that general Markov processes admit a unified and self-consistent non-equilibrium steady-state thermodynamic formalism, regardless of underlying detailed models.

We investigate numerically temporal correlations in a one-dimensional critical-slope sandpile model with rules that on average conserve the number of particles. Our work is motivated by the existence of two well-separated time scales in self-organized sandpile models, one related to the spreading of avalanches and the other imposed by the external driving. We assume that avalanches are instantaneous events on the time scale imposed by the external deposition and study the autocorrelation function of the series of successive avalanche amplitudes. We find that the autocorrelation function has a log-normal form and for large system sizes tends to a constant, implying that the temporal correlations become stronger in the limit of large system size. We independently test this result by calculating the power spectrum of the series of successive avalanche lifetimes and sizes. For large system sizes L there is a frequency regime where the power spectrum tends to a 1/f type of noise, in agreement with the tendency of the autocorrelation function to approach a constant in large systems. @S1063-651X~96!07212-1#

We discuss statics and dynamics of condensation in a zero-range process with compartments of limited sizes. For the symmetric dynamics the stationary state has a factorized form. For the asymmetric dynamics the steady state factorizes only for special hopping rules which allow for overjumps of fully occupied compartments. In the limit of large system size the grand canonical analysis is exact also in a condensed phase, and for a broader class of hopping rates as compared to the previously studied systems with infinite compartments. The dynamics of condensation exhibits dynamical self-blocking, which significantly prolongs relaxation times. These general features are illustrated with a concrete example: an inhomogeneous system with hopping rates that result in Bose-Einstein-like condensations.

In recent years, small-world graphs have gained considerable interest as models of real-world systems, which often display features residing between regularity and randomness. The most notable of these models is the Watts-Strogatz graph [1], though alternatives have been proposed [2]. The unifying characteristics of these models are that any two nodes are joined with a small number of links between them (i.e. short path length), while at the same time connected node pairs exhibit an abundance of triangular relations resulting in a high degree of local redundancy (i.e. high clustering). Theoretical investigations of small-world graph models have generally applied asymptotic evaluations in the limit of large system size [3] or the continuum approximation [4] to the algorithmic definition of the graph, in the absence of an analytic representation. In this study, we introduce a generative model of directed small-world graphs, a canonical model of Watts-Strogatz digraphs, and propose an approach that yields the graph’s defining adjacency matrix in algebraic terms, with the goal to provide mathematically rigorous access to the study of finite-size small-world graphs [5]. The proposed approach makes use of random annihilation operators whose algebraic properties can be utilized to assess algebraically well-defined graph-theoretic measures in an analytically exact framework, valid nonasymptotically for all graph sizes. We demonstrated the application of our approach by calculating, for the first time, the asymmetry index and total clustering coefficient of small worlds in an exact fashion. We then utilize the exact nonasymptotic expression for the clustering coefficient in order to assess the small-worldness of structural brain networks in an analytic setting. Using the number of nodes and edges of the given brain networks to construct the equivalent small-world network, we observe that a significant edge rewiring of at least 20% up to 60% is required to produce the small-worldness indices observed in these networks. Importantly, the maximum of the small-worldness index however occurs in all cases at one order of magnitude lower than the required rewiring found. This result suggests that neural graphs reside far away from the small-world regime of the Watts-Strogatz model.

SZNAJD MODEL WITH SYNCHRONOUS UPDATING ON COMPLEX NETWORKS

Recently, Sá, Ribeiro, and Prosen [Phys. Rev. X 10, 021019 (2020)10.1103/PhysRevX.10.021019] introduced complex spacing ratios to analyze eigenvalue correlations in non-Hermitian systems. At present there are no analytical results for the probability distribution of these ratios in the limit of large system size. We derive an approximation formula for the Ginibre universality class of random matrix theory which converges exponentially fast to the limit of infinite matrix size. We also give results for moments of the distribution in this limit.

The presence of a global internal symmetry in a quantum many-body system is reflected in the fact that the entanglement between its subparts is endowed with an internal structure, namely it can be decomposed as a sum of contributions associated to each symmetry sector. The symmetry resolution of entanglement measures provides a formidable tool to probe the out-of-equilibrium dynamics of quantum systems. Here, we study the time evolution of charge-imbalance-resolved negativity after a global quench in the context of free-fermion systems, complementing former works for the symmetry-resolved entanglement entropy. We find that the charge-imbalance-resolved logarithmic negativity shows an effective equipartition in the scaling limit of large times and system size, with a perfect equipartition for early and infinite times. We also derive and conjecture a formula for the dynamics of the charged Rényi logarithmic negativities. We argue that our results can be understood in the framework of the quasiparticle picture for the entanglement dynamics, and provide a conjecture that we expect to be valid for generic integrable models.

Surface dynamics of the Wolf-Villain model for epitaxial growth in 1 + 1 dimensions

On the bifurcation structure of the mean-field fluctuation in the globally coupled tent map systems

We consider the particle current in the asymmetric avalanche process on a ring. It is known to exhibit a transition from the intermittent to continuous flow at the critical density of particles. The exact expressions for the first two scaled cumulants of the particle current are obtained in the large time limit t → ∞ via the Bethe ansatz and a perturbative solution of the TQ-equation. The results are presented in an integral form suitable for the asymptotic analysis in the large system size limit N → ∞. In this limit the first cumulant, the average current per site or the average velocity of the associated interface, is asymptotically finite below the critical density and grows linearly and exponentially times power law prefactor at the critical density and above, respectively. The scaled second cumulant per site, i.e. the diffusion coefficient or the scaled variance of the associated interface height, shows the O(N −1/2) decay expected for models in the Kardar–Parisi–Zhang universality class below the critical density, while it is growing as O(N 3/2) and exponentially times power law prefactor at the critical point and above. Also, we identify the crossover regime and obtain the scaling functions for the uniform asymptotics unifying the three regimes. These functions are compared to the scaling functions describing crossover of the cumulants of the avalanche size, obtained as statistics of the first return area under the time space trajectory of the Vasicek random process.

In most models of the spread of disease over contact networks it is assumed that the probabilities per unit time of disease transmission and recovery from disease are constant, implying exponential distributions of the time intervals for transmission and recovery. Time intervals for real diseases, however, have distributions that in most cases are far from exponential, which leads to disagreements, both qualitative and quantitative, with the models. In this paper, we study a generalized version of the susceptible-infected-recovered model of epidemic disease that allows for arbitrary distributions of transmission and recovery times. Standard differential equation approaches cannot be used for this generalized model, but we show that the problem can be reformulated as a time-dependent message passing calculation on the appropriate contact network. The calculation is exact on trees (i.e., loopless networks) or locally treelike networks (such as random graphs) in the large system size limit. On non-tree-like networks we show that the calculation gives a rigorous bound on the size of disease outbreaks. We demonstrate the method with applications to two specific models and the results compare favorably with numerical simulations.

In recent analytical work, Biskup et al. (Europhys. Lett., 60 (2002) 21) studied the behaviour of d-dimensional finite-volume liquid-vapour systems at a fixed excess δN of particles above the ambient gas density. By identifying a dimensionless parameter Δ(δN) and a universal constant Δc(d), they showed in the limit of large system sizes that for Δ < Δc the excess is absorbed in the background (“evaporated” system), while for Δ > Δc a droplet of the dense phase occurs (“condensed” system). Also the fraction λΔ of excess particles forming the droplet is given explicitly. Furthermore, they argue that the same is true for solid-gas systems. By making use of the well-known equivalence of the lattice-gas picture with the spin-(1/2) Ising model, we performed Monte Carlo simulations of the Ising model with nearest-neighbour couplings on a square lattice with periodic boundary conditions at fixed magnetisation, corresponding to a fixed particles excess. To test the applicability of the analytical results to much smaller, practically accessible system sizes, we measured the largest minority droplet, corresponding to the solid phase, at various system sizes (L = 40,…,640). Using analytic values for the spontaneous magnetisation m0, the susceptibility χ and the Wulff interfacial free energy density τW for the infinite system, we were able to determine λΔ numerically in very good agreement with the theoretical prediction.

We develop a formalism to describe the discrete-time dynamics of systems containing an arbitrary number of interacting species. The individual-based model, which forms our starting point, is described by a Markov chain, which in the limit of large system sizes is shown to be very well-approximated by a Fokker-Planck-like equation, or equivalently by a set of stochastic difference equations. This formalism is applied to the specific case of two species: one predator species and its prey species. Quasicycles, stochastic cycles sustained and amplified by the demographic noise, previously found in continuous-time predator-prey models are shown to exist, and their behavior predicted from a linear noise analysis is shown to be in very good agreement with simulations. The effects of the noise on other attractors in the corresponding deterministic map, such as periodic cycles, quasiperiodicity, and chaos, are also investigated.

Predicting properties of complex, large-scale quantum systems is essential for developing quantum technologies. We present an efficient method for constructing an approximate classical description of a quantum state using very few measurements of the state. This description, called a classical shadow, can be used to predict many different properties: order $\log M$ measurements suffice to accurately predict $M$ different functions of the state with high success probability. The number of measurements is independent of the system size, and saturates information-theoretic lower bounds. Moreover, target properties to predict can be selected after the measurements are completed. We support our theoretical findings with extensive numerical experiments. We apply classical shadows to predict quantum fidelities, entanglement entropies, two-point correlation functions, expectation values of local observables, and the energy variance of many-body local Hamiltonians. The numerical results highlight the advantages of classical shadows relative to previously known methods.