624 research outputs found
Complexity, parallel computation and statistical physics
The intuition that a long history is required for the emergence of complexity
in natural systems is formalized using the notion of depth. The depth of a
system is defined in terms of the number of parallel computational steps needed
to simulate it. Depth provides an objective, irreducible measure of history
applicable to systems of the kind studied in statistical physics. It is argued
that physical complexity cannot occur in the absence of substantial depth and
that depth is a useful proxy for physical complexity. The ideas are illustrated
for a variety of systems in statistical physics.Comment: 21 pages, 7 figure
Pairwise MRF Calibration by Perturbation of the Bethe Reference Point
We investigate different ways of generating approximate solutions to the
pairwise Markov random field (MRF) selection problem. We focus mainly on the
inverse Ising problem, but discuss also the somewhat related inverse Gaussian
problem because both types of MRF are suitable for inference tasks with the
belief propagation algorithm (BP) under certain conditions. Our approach
consists in to take a Bethe mean-field solution obtained with a maximum
spanning tree (MST) of pairwise mutual information, referred to as the
\emph{Bethe reference point}, for further perturbation procedures. We consider
three different ways following this idea: in the first one, we select and
calibrate iteratively the optimal links to be added starting from the Bethe
reference point; the second one is based on the observation that the natural
gradient can be computed analytically at the Bethe point; in the third one,
assuming no local field and using low temperature expansion we develop a dual
loop joint model based on a well chosen fundamental cycle basis. We indeed
identify a subclass of planar models, which we refer to as \emph{Bethe-dual
graph models}, having possibly many loops, but characterized by a singly
connected dual factor graph, for which the partition function and the linear
response can be computed exactly in respectively O(N) and operations,
thanks to a dual weight propagation (DWP) message passing procedure that we set
up. When restricted to this subclass of models, the inverse Ising problem being
convex, becomes tractable at any temperature. Experimental tests on various
datasets with refined or regularization procedures indicate that
these approaches may be competitive and useful alternatives to existing ones.Comment: 54 pages, 8 figure. section 5 and refs added in V
Searching for collective behavior in a network of real neurons
Maximum entropy models are the least structured probability distributions
that exactly reproduce a chosen set of statistics measured in an interacting
network. Here we use this principle to construct probabilistic models which
describe the correlated spiking activity of populations of up to 120 neurons in
the salamander retina as it responds to natural movies. Already in groups as
small as 10 neurons, interactions between spikes can no longer be regarded as
small perturbations in an otherwise independent system; for 40 or more neurons
pairwise interactions need to be supplemented by a global interaction that
controls the distribution of synchrony in the population. Here we show that
such "K-pairwise" models--being systematic extensions of the previously used
pairwise Ising models--provide an excellent account of the data. We explore the
properties of the neural vocabulary by: 1) estimating its entropy, which
constrains the population's capacity to represent visual information; 2)
classifying activity patterns into a small set of metastable collective modes;
3) showing that the neural codeword ensembles are extremely inhomogenous; 4)
demonstrating that the state of individual neurons is highly predictable from
the rest of the population, allowing the capacity for error correction.Comment: 24 pages, 19 figure
Next nearest neighbour Ising models on random graphs
This paper develops results for the next nearest neighbour Ising model on
random graphs. Besides being an essential ingredient in classic models for
frustrated systems, second neighbour interactions interactions arise naturally
in several applications such as the colour diversity problem and graphical
games. We demonstrate ensembles of random graphs, including regular
connectivity graphs, that have a periodic variation of free energy, with either
the ratio of nearest to next nearest couplings, or the mean number of nearest
neighbours. When the coupling ratio is integer paramagnetic phases can be found
at zero temperature. This is shown to be related to the locked or unlocked
nature of the interactions. For anti-ferromagnetic couplings, spin glass phases
are demonstrated at low temperature. The interaction structure is formulated as
a factor graph, the solution on a tree is developed. The replica symmetric and
energetic one-step replica symmetry breaking solution is developed using the
cavity method. We calculate within these frameworks the phase diagram and
demonstrate the existence of dynamical transitions at zero temperature for
cases of anti-ferromagnetic coupling on regular and inhomogeneous random
graphs.Comment: 55 pages, 15 figures, version 2 with minor revisions, to be published
J. Stat. Mec
Dynamics of a Quantum Phase Transition and Relaxation to a Steady State
We review recent theoretical work on two closely related issues: excitation
of an isolated quantum condensed matter system driven adiabatically across a
continuous quantum phase transition or a gapless phase, and apparent relaxation
of an excited system after a sudden quench of a parameter in its Hamiltonian.
Accordingly the review is divided into two parts. The first part revolves
around a quantum version of the Kibble-Zurek mechanism including also phenomena
that go beyond this simple paradigm. What they have in common is that
excitation of a gapless many-body system scales with a power of the driving
rate. The second part attempts a systematic presentation of recent results and
conjectures on apparent relaxation of a pure state of an isolated quantum
many-body system after its excitation by a sudden quench. This research is
motivated in part by recent experimental developments in the physics of
ultracold atoms with potential applications in the adiabatic quantum state
preparation and quantum computation.Comment: 117 pages; review accepted in Advances in Physic
The Quantum Adiabatic Algorithm applied to random optimization problems: the quantum spin glass perspective
Among various algorithms designed to exploit the specific properties of
quantum computers with respect to classical ones, the quantum adiabatic
algorithm is a versatile proposition to find the minimal value of an arbitrary
cost function (ground state energy). Random optimization problems provide a
natural testbed to compare its efficiency with that of classical algorithms.
These problems correspond to mean field spin glasses that have been extensively
studied in the classical case. This paper reviews recent analytical works that
extended these studies to incorporate the effect of quantum fluctuations, and
presents also some original results in this direction.Comment: 151 pages, 21 figure
Theory of High-Force DNA Stretching and Overstretching
Single molecule experiments on single- and double stranded DNA have sparked a
renewed interest in the force-extension of polymers. The extensible Freely
Jointed Chain (FJC) model is frequently invoked to explain the observed
behavior of single-stranded DNA. We demonstrate that this model does not
satisfactorily describe recent high-force stretching data. We instead propose a
model (the Discrete Persistent Chain, or ``DPC'') that borrows features from
both the FJC and the Wormlike Chain, and show that it resembles the data more
closely. We find that most of the high-force behavior previously attributed to
stretch elasticity is really a feature of the corrected entropic elasticity;
the true stretch compliance of single-stranded DNA is several times smaller
than that found by previous authors. Next we elaborate our model to allow
coexistence of two conformational states of DNA, each with its own stretch and
bend elastic constants. Our model is computationally simple, and gives an
excellent fit through the entire overstretching transition of nicked,
double-stranded DNA. The fit gives the first values for the elastic constants
of the stretched state. In particular we find the effective bend stiffness for
DNA in this state to be about 10 nm*kbt, a value quite different from either
B-form or single-stranded DNAComment: 33 pages, 11 figures. High-quality figures available upon reques
- …