Quantum Belief Propagation

We present an accurate numerical algorithm, called quantum belief propagation (QBP), for simulation of one-dimensional quantum systems at non-zero temperature. The algorithm exploits the fact that quantum effects are short-range in these systems at non-zero temperature, decaying on a length scale inversely proportional to the temperature. We compare to exact results on a spin-1/2 Heisenberg chain. Even a very modest calculation, requiring diagonalizing only 10-by-10 matrices, reproduces the peak susceptibility with a relative error of less than $10^{-5}$, while more elaborate calculations further reduce the error.Comment: 4 pages, 1 figure; revised time estimates due to improved implementation. Typographical corrections to Eq. 7 made; thanks to David Poulin for pointing out the mistak

Diffusion Processes on Power-Law Small-World Networks

We consider diffusion processes on power-law small-world networks in different dimensions. In one dimension, we find a rich phase diagram, with different transient and recurrent phases, including a critical line with continuously varying exponents. The results were obtained using self-consistent perturbation theory and can also be understood in terms of a scaling theory, which provides a general framework for understanding processes on small-world networks with different distributions of long-range links.Comment: 4 pages, 3 figures, added references, modified Fig. 2 with added data (PRL, in press

Finite Size Scaling of Mutual Information: A Scalable Simulation

We develop a quantum Monte Carlo procedure to compute the Renyi mutual information of an interacting quantum many-body system at non-zero temperature. Performing simulations on a spin-1/2 XXZ model, we observe that for a subregion of fixed size embedded in a system of size L, the mutual information converges at large L to a limiting function which displays non-monotonic temperature behavior corresponding to the onset of correlations. For a region of size L/2 embedded in a system of size L, the mutual information divided by L converges to a limiting function of temperature, with apparently nontrivial corrections near critical points.Comment: 4 pages, 4 figure

Observations Outside the Light-Cone: Algorithms for Non-Equilibrium and Thermal States

We apply algorithms based on Lieb-Robinson bounds to simulate time-dependent and thermal quantities in quantum systems. For time-dependent systems, we modify a previous mapping to quantum circuits to significantly reduce the computer resources required. This modification is based on a principle of "observing" the system outside the light-cone. We apply this method to study spin relaxation in systems started out of equilibrium with initial conditions that give rise to very rapid entanglement growth. We also show that it is possible to approximate time evolution under a local Hamiltonian by a quantum circuit whose light-cone naturally matches the Lieb-Robinson velocity. Asymptotically, these modified methods allow a doubling of the system size that one can obtain compared to direct simulation. We then consider a different problem of thermal properties of disordered spin chains and use quantum belief propagation to average over different configurations. We test this algorithm on one dimensional systems with mixed ferromagnetic and anti-ferromagnetic bonds, where we can compare to quantum Monte Carlo, and then we apply it to the study of disordered, frustrated spin systems.Comment: 19 pages, 12 figure

Tip Splittings and Phase Transitions in the Dielectric Breakdown Model: Mapping to the DLA Model

We show that the fractal growth described by the dielectric breakdown model exhibits a phase transition in the multifractal spectrum of the growth measure. The transition takes place because the tip-splitting of branches forms a fixed angle. This angle is eta dependent but it can be rescaled onto an ``effectively'' universal angle of the DLA branching process. We derive an analytic rescaling relation which is in agreement with numerical simulations. The dimension of the clusters decreases linearly with the angle and the growth becomes non-fractal at an angle close to 74 degrees (which corresponds to eta= 4.0 +- 0.3).Comment: 4 pages, REVTex, 3 figure

Anomalies in the Entanglement Properties of the Square Lattice Heisenberg Model

We compute the bipartite entanglement properties of the spin-half square-lattice Heisenberg model by a variety of numerical techniques that include valence bond quantum Monte Carlo (QMC), stochastic series expansion QMC, high temperature series expansions and zero temperature coupling constant expansions around the Ising limit. We find that the area law is always satisfied, but in addition to the entanglement entropy per unit boundary length, there are other terms that depend logarithmically on the subregion size, arising from broken symmetry in the bulk and from the existence of corners at the boundary. We find that the numerical results are anomalous in several ways. First, the bulk term arising from broken symmetry deviates from an exact calculation that can be done for a mean-field Neel state. Second, the corner logs do not agree with the known results for non-interacting Boson modes. And, third, even the finite temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that T->0 and L->infinity limits do not commute. These calculations show that entanglement entropy demonstrates a very rich behavior in d>1, which deserves further attention.Comment: 12 pages, 7 figures, 2 tables. Numerical values in Table I correcte

Exact Multifractal Spectra for Arbitrary Laplacian Random Walks

Iterated conformal mappings are used to obtain exact multifractal spectra of the harmonic measure for arbitrary Laplacian random walks in two dimensions. Separate spectra are found to describe scaling of the growth measure in time, of the measure near the growth tip, and of the measure away from the growth tip. The spectra away from the tip coincide with those of conformally invariant equilibrium systems with arbitrary central charge $c\leq 1$, with $c$ related to the particular walk chosen, while the scaling in time and near the tip cannot be obtained from the equilibrium properties.Comment: 4 pages, 3 figures; references added, minor correction

Random Vibrational Networks and Renormalization Group

We consider the properties of vibrational dynamics on random networks, with random masses and spring constants. The localization properties of the eigenstates contrast greatly with the Laplacian case on these networks. We introduce several real-space renormalization techniques which can be used to describe this dynamics on general networks, drawing on strong disorder techniques developed for regular lattices. The renormalization group is capable of elucidating the localization properties, and provides, even for specific network instances, a fast approximation technique for determining the spectra which compares well with exact results.Comment: 4 pages, 3 figure

Community Detection as an Inference Problem

We express community detection as an inference problem of determining the most likely arrangement of communities. We then apply belief propagation and mean-field theory to this problem, and show that this leads to fast, accurate algorithms for community detection.Comment: 4 pages, 2 figure