Source author record

R. Mulet

R. Mulet appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.dis-nn cond-mat.stat-mech Artificial Intelligence Biological Physics Genomics Information Theory math.IT Molecular Networks physics.chem-ph Quantitative Methods

Catalog footprint

What is connected

9works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

The Cavity Master Equation: average and fixed point of the ferromagnetic model in random graphs

The Cavity Master Equation (CME) is a closure scheme to the usual Master Equation representing the dynamics of discrete variables in continuous time. In this work we explore the CME for a ferromagnetic model in a random graph. We first derive and average equation of the CME that describes the dynamics of mean magnetization of the system. We show that the numerical results compare remarkably well with the Monte Carlo simulations. Then, we show that the stationary state of the CME is well described by BP-like equations (independently of the dynamic rules that let the system towards the stationary state). These equations may be rewritten exactly as the fixed point solutions of the Cavity Equation if one also assumes that the stationary state is well described by a Boltzmann distribution.

preprint2014arXiv

Message passing and Monte Carlo algorithms: connecting fixed points with metastable states

Mean field-like approximations (including naive mean field, Bethe and Kikuchi and more general Cluster Variational Methods) are known to stabilize ordered phases at temperatures higher than the thermodynamical transition. For example, in the Edwards-Anderson model in 2-dimensions these approximations predict a spin glass transition at finite $T$. Here we show that the spin glass solutions of the Cluster Variational Method (CVM) at plaquette level do describe well actual metastable states of the system. Moreover, we prove that these states can be used to predict non trivial statistical quantities, like the distribution of the overlap between two replicas. Our results support the idea that message passing algorithms can be helpful to accelerate Monte Carlo simulations in finite dimensional systems.

preprint2013arXiv

Identifying all irreducible conserved metabolite pools in genome-scale metabolic networks: a general method and the case of Escherichia coli

The stoichiometry of metabolic networks usually gives rise to a family of conservation laws for the aggregate concentration of specific pools of metabolites, which not only constrain the dynamics of the network, but also provide key insight into a cell's production capabilities. When the conserved quantity identifies with a chemical moiety, extracting all such conservation laws from the stoichiometry amounts to finding all integer solutions to an NP-hard programming problem. Here we propose a novel and efficient computational strategy that combines Monte Carlo, message passing, and relaxation algorithms to compute the complete set of irreducible integer conservation laws of a given stoichiometric matrix, also providing a certificate for correctness and maximality of the solution. The method is deployed for the analysis of the complete set of irreducible integer pools of two large-scale reconstructions of the metabolism of the bacterium Escherichia coli in different growth media. In addition, we uncover a scaling relation that links the size of the irreducible pool basis to the number of metabolites, for which we present an analytical explanation.

preprint2011arXiv

Characterizing and Improving Generalized Belief Propagation Algorithms on the 2D Edwards-Anderson Model

We study the performance of different message passing algorithms in the two dimensional Edwards Anderson model. We show that the standard Belief Propagation (BP) algorithm converges only at high temperature to a paramagnetic solution. Then, we test a Generalized Belief Propagation (GBP) algorithm, derived from a Cluster Variational Method (CVM) at the plaquette level. We compare its performance with BP and with other algorithms derived under the same approximation: Double Loop (DL) and a two-ways message passing algorithm (HAK). The plaquette-CVM approximation improves BP in at least three ways: the quality of the paramagnetic solution at high temperatures, a better estimate (lower) for the critical temperature, and the fact that the GBP message passing algorithm converges also to non paramagnetic solutions. The lack of convergence of the standard GBP message passing algorithm at low temperatures seems to be related to the implementation details and not to the appearance of long range order. In fact, we prove that a gauge invariance of the constrained CVM free energy can be exploited to derive a new message passing algorithm which converges at even lower temperatures. In all its region of convergence this new algorithm is faster than HAK and DL by some orders of magnitude.

preprint2010arXiv

Replica Cluster Variational Method

We present a general formalism to make the Replica-Symmetric and Replica-Symmetry-Breaking ansatz in the context of Kikuchi's Cluster Variational Method (CVM). Using replicas and the message-passing formulation of CVM we obtain a variational expression of the replicated free energy of a system with quenched disorder, both averaged and on a single sample, and make the hierarchical ansatz using functionals of functions of fields to represent the messages. We begin to study the method considering the plaquette approximation to the averaged free energy of the Edwards-Anderson model in the paramagnetic Replica-Symmetric phase. In two dimensions we find that the spurious spin-glass phase transition of the Bethe approximation disappears and the paramagnetic phase is stable down to zero temperature in all the three regular 2D lattices. The quantitative estimates of the free energy and of various other quantities improve those of the Bethe approximation. We provide the physical interpretation of the beliefs in the replica-symmetric phase as disorder distributions of the local Hamiltonians. The messages instead do not admit such an interpretation and indeed they cannot be represented as populations in the spin-glass phase at variance with the Bethe approximation.

preprint2005arXiv

A New Simulated Annealing Algorithm for the Multiple Sequence Alignment Problem: The approach of Polymers in a Random Media

We proposed a probabilistic algorithm to solve the Multiple Sequence Alignment problem. The algorithm is a Simulated Annealing (SA) that exploits the representation of the Multiple Alignment between $D$ sequences as a directed polymer in $D$ dimensions. Within this representation we can easily track the evolution in the configuration space of the alignment through local moves of low computational cost. At variance with other probabilistic algorithms proposed to solve this problem, our approach allows for the creation and deletion of gaps without extra computational cost. The algorithm was tested aligning proteins from the kinases family. When D=3 the results are consistent with those obtained using a complete algorithm. For $D>3$ where the complete algorithm fails, we show that our algorithm still converges to reasonable alignments. Moreover, we study the space of solutions obtained and show that depending on the number of sequences aligned the solutions are organized in different ways, suggesting a possible source of errors for progressive algorithms.

preprint2004arXiv

Adaptive drivers in a model of urban traffic

We introduce a simple lattice model of traffic flow in a city where drivers optimize their route-selection in time in order to avoid traffic jams, and study its phase structure as a function of the density of vehicles and of the drivers' behavioral parameters via numerical simulations and mean-field analytical arguments. We identify a phase transition between a low- and a high-density regime. In the latter, inductive drivers may surprisingly behave worse than randomly selecting drivers.

preprint2003arXiv

Quasi-equilibrium during aging of the 2D Edwards-Anderson model

We test the quasi-equilibrium picture of the aging dynamics -strictly valid in the asymptotic dynamical regime of aging systems- in the pre-asymptotic aging regime of the two dimensional Edwards-Anderson spin glass model. We compare the fluctuation-dissipation characteristic for spin autocorrelation function and response with a corresponding one obtained for a suitably defined new correlation function and its conjugated response. In agreement with the quasi-equilibrium picture we find that after a short transient the two corresponding fluctuation-dissipation ratios (FDR) coincide at equal times. Moreover we show that, as it happens for the usual FDR, the new dynamic FDR at finite time coincides with the static one at finite size.

preprint2001arXiv

Kob-Andersen model: a non-standard mechanism for the glassy transition

We present new results reflecting the analogies between the Kob-Andersen model and other glassy systems. Studying the stability of the blocked configurations above and below the transition we also give arguments that supports their relevance for the glassy behaviour of the model. However we find, surprisingly, that the organization of the phase space of the system is different from the well known organization of other mean-field spin glasses and structural glasses.