Source author record

Alessia Annibale

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

14 published item(s)

preprint2022arXiv

Exact and Approximate Mean First Passage Times on Trees and other Necklace Structures: a Local Equilibrium Approach

In this work we propose a novel method to calculate mean first-passage times (MFPTs) for random walks on graphs, based on a dimensionality reduction technique for Markov State Models, known as local-equilibrium (LE). We show that for a broad class of graphs, which includes trees, LE coarse-graining preserves the MFPTs between certain nodes, upon making a suitable choice of the coarse-grained states (or clusters). We prove that this relation is exact for graphs that can be coarse-grained into a one-dimensional lattice where each cluster connects to the lattice only through a single node of the original graph. A side result of the proof generalises the well-known essential edge lemma (EEL), which is valid for reversible random walks, to irreversible walkers. Such a generalised EEL leads to explicit formulae for the MFPTs between certain nodes in this class of graphs. For graphs that do not fall in this class, the generalised EEL provides useful approximations if the graph allows a one-dimensional coarse-grained representation and the clusters are sparsely interconnected. We first demonstrate our method for the simple random walk on the $c$-ary tree, then we consider other graph structures and more general random walks, including irreversible random walks.

preprint2022arXiv

Reproducing size distributions of swarms of barchan dunes on Mars and Earth using a mean-field model

We apply a mean-field model of interactions between migrating barchan dunes, the CAFE model, which includes spontaneous calving, aggregation, fragmentation, and mass-exchange, yielding a steady-state size distribution that can be resolved for different choices of interaction parameters. The CAFE model is applied to empirically measured distributions of dune sizes in two barchan swarms in the north circumpolar region of Mars, three swarms in Morocco, and one in Mauritania, each containing more than 1000 bedforms. When the sizes of bedforms are rescaled by the mean size in each zone two attractor states appear, with the Tarfaya zones all displaying a common distribution and the Martian and Mauritanian zones sharing a different distribution. Comparison of these attractor states with the outputs of the CAFE model reveals that the Tarfaya-type distribution results from a preference for aggregation and fragmentation interactions whereas the Mars-Mauritania distribution is more likely a result of exchange-dominated interactions. We observe that there appears to be a greater number of collisions happening in Tarfaya than in the other areas which is consistent with a greater rate of aggregation-fragmentation processes as suggested by our model. Our comparison with the CAFE model also predicts a universal rule for the outputs of the different types of interactions with exchange interactions favouring the production of two dunes roughly equal in size. Fragmentation interactions often result in three bedforms with volumes approximately in the ratio 1:1:2. Finally, we find that spontaneous calving of dunes does not play an important role in shaping the size distributions in barchan swarms.

preprint2022arXiv

Uncovering the non-equilibrium stationary properties in sparse Boolean networks

Dynamic processes of interacting units on a network are out of equilibrium in general. In the case of a directed tree, the dynamic cavity method provides an efficient tool that characterises the dynamic trajectory of the process for the linear threshold model. However, because of the computational complexity of the method, the analysis has been limited to systems where the largest number of neighbours is small. We devise an efficient implementation of the dynamic cavity method which substantially reduces the computational complexity of the method for systems with discrete couplings. Our approach opens up the possibility to investigate the dynamic properties of networks with fat-tailed degree distribution. We exploit this new implementation to study properties of the non-equilibrium steady-state. We extend the dynamical cavity approach to calculate the pairwise correlations induced by different motifs in the network. Our results suggest that just two basic motifs of the network are able to accurately describe the entire statistics of observed correlations. Finally, we investigate models defined on networks containing bi-directional interactions. We observe that the stationary state associated with networks with symmetric or anti-symmetric interactions is biased towards the active or inactive state respectively, even if independent interaction entries are drawn from a symmetric distribution. This phenomenon, which can be regarded as a form of spontaneous symmetry-breaking, is peculiar to systems formulated in terms of Boolean variables, as opposed to Ising spins.

preprint2021arXiv

Vaccination with partial transmission and social distancing on contact networks

We study the impact of vaccination on the risk of epidemics spreading through structured networks using the cavity method of statistical physics. We relax the assumption that vaccination prevents all transmission of a disease used in previous studies, such that vaccinated nodes have a small probability of transmission. To do so we extend the cavity method to study networks where nodes have heterogeneous transmissibility. We find that vaccination with partial transmission still provides herd immunity and show how the herd immunity threshold depends upon the assortativity between nodes of different transmissibility. In addition, we study the impact of social distancing via bond percolation and show that percolation targeting links between nodes of high transmissibility can reduce the risk of an epidemic greater than targeting links between nodes of high degree. Finally, we extend recent methods to compute the distributional equations of risk in populations with heterogeneous transmissibility and show how targeted social distancing measures may reduce overall risk greater than untargeted vaccination campaigns, by comparing the effect of random and targeted strategies of node and link deletion on the risk distribution.

preprint2020arXiv

Percolation on the gene regulatory network

We consider a simplified model for gene regulation, where gene expression is regulated by transcription factors (TFs), which are single proteins or protein complexes. Proteins are in turn synthesised from expressed genes, creating a feedback loop of regulation. This leads to a directed bipartite network in which a link from a gene to a TF exists if the gene codes for a protein contributing to the TF, and a link from a TF to a gene exists if the TF regulates the expression of the gene. Both genes and TFs are modelled as binary variables, which indicate, respectively, whether a gene is expressed or not, and a TF is synthesised or not. We consider the scenario where for a TF to be synthesised, all of its contributing genes must be expressed. This results in an ``AND'' gate logic for the dynamics of TFs. By adapting percolation theory to directed bipartite graphs, evolving according to the AND logic dynamics, we are able to determine the necessary conditions, in the network parameter space, under which bipartite networks can support a multiplicity of stable gene expression patterns, under noisy conditions, as required in stable cell types. In particular, the analysis reveals the possibility of a bi-stability region, where the extensive percolating cluster is or is not resilient to perturbations. This is remarkably different from the transition observed in standard percolation theory. Finally, we consider perturbations involving single node removal that mimic gene knockout experiments. Results reveal the strong dependence of the gene knockout cascade on the logic implemented in the underlying network dynamics, highlighting in particular that avalanche sizes cannot be easily related to gene-gene interaction networks.

preprint2019arXiv

Correlation Functions, Mean First Passage Times and the Kemeny Constant

Markov processes are widely used models for investigating kinetic networks. Here we collate and present a variety of results pertaining to kinetic network models, in a unified framework. The aim is to lay out explicit links between several important quantities commonly studied in the field, including mean first passage times (MFPTs), correlation functions and the Kemeny constant, and highlight some of the subtleties which are often overlooked in the literature, while providing new insights. Results include (i) a simple physical interpretation of the Kemeny constant, (ii) a recipe to infer equilibrium distributions and rate matrices from measurements of MFPTs, potentially useful in applications, including milestoning in molecular dynamics, and (iii) a protocol to reduce the dimensionality of kinetic networks, based on specific requirements that the MFPTs in the coarse-grained system should satisfy. It is proven that this protocol coincides with the one proposed by Hummer and Szabo in [1] and it leads to a variational principle for the Kemeny constant. We hope that this study will serve as a useful reference for readers interested in theoretical aspects of kinetic networks, some of which underpin useful applications, including milestoning and coarse-graining.

preprint2016arXiv

The role of idiotypic interactions in the adaptive immune system: a belief-propagation approach

In this work we use belief-propagation techniques to study the equilibrium behaviour of a minimal model for the immune system comprising interacting T and B clones. We investigate the effect of the so-called idiotypic interactions among complementary B clones on the system's activation. Our result shows that B-B interactions increase the system's resilience to noise, making clonal activation more stable, while increasing the cross-talk between different clones. We derive analytically the noise level at which a B clone gets activated, in the absence of cross-talk, and find that this increases with the strength of idiotypic interactions and with the number of T cells signalling the B clone. We also derive, analytically and numerically, via population dynamics, the critical line where clonal cross-talk arises. Our approach allows us to derive the B clone size distribution, which can be experimentally measured and gives important information about the adaptive immune system response to antigens and vaccination.

preprint2015arXiv

A dynamical model of the adaptive immune system: effects of cells promiscuity, antigens and B-B interactions

We analyse a minimal model for the primary response in the adaptive immune system comprising three different players: antigens, T and B cells. We assume B-T interactions to be diluted and sampled locally from heterogeneous degree distributions, which mimic B cells receptors' promiscuity. We derive dynamical equations for the order parameters quantifying the B cells activation and study the nature and stability of the stationary solutions using linear stability analysis and Monte Carlo simulations.The system's behaviour is studied in different scaling regimes of the number of B cells, dilution in the interactions and number of antigens. Our analysis shows that: (i) B cells activation depends on the number of receptors in such a way that cells with an insufficient number of triggered receptors cannot be activated; (ii) idiotypic (i.e. B-B) interactions enhance parallel activation of multiple clones, improving the system's ability to fight different pathogens in parallel; (iii) the higher the fraction of antigens within the host the harder is for the system to sustain parallel signalling to B cells, crucial for the homeostatic control of cell numbers.

preprint2015arXiv

Protein interaction networks and biology: towards the connection

Protein interaction networks (PIN) are popular means to visualize the proteome. However, PIN datasets are known to be noisy, incomplete and biased by the experimental protocols used to detect protein interactions. This paper aims at understanding the connection between true protein interactions and the protein interaction datasets that have been obtained using the most popular experimental techniques, i.e. mass spectronomy (MS) and yeast two-hybrid (Y2H). We show that the most natural adjacency matrix of protein interaction networks has a separable form, and this induces precise relations between moments of the degree distribution and the number of short loops. These relations provide powerful tools to test the reliability of datasets and hint at the underlying biological mechanism with which proteins and complexes recruit each other.

preprint2015arXiv

The two-star model: exact solution in the sparse regime and condensation transition

The $2$-star model is the simplest exponential random graph model that displays complex behavior, such as degeneracy and phase transition. Despite its importance, this model has been solved only in the regime of dense connectivity. In this work we solve the model in the finite connectivity regime, far more prevalent in real world networks. We show that the model undergoes a condensation transition from a liquid to a condensate phase along the critical line corresponding, in the ensemble parameters space, to the Erdös-Rényi graphs. In the fluid phase the model can produce graphs with a narrow degree statistics, ranging from regular to Erdös-Rényi graphs, while in the condensed phase, the "excess" degree heterogeneity condenses on a single site with degree $\sim\sqrt{N}$. This shows the unsuitability of the two-star model, in its standard definition, to produce arbitrary finitely connected graphs with degree heterogeneity higher than Erdös-Rényi graphs and suggests that non-pathological variants of this model may be attained by softening the global constraint on the two-stars, while keeping the number of links hardly constrained.

preprint2014arXiv

Associative networks with diluted patterns: dynamical analysis at low and medium load

In this work we solve the dynamics of pattern diluted associative networks, evolving via sequential Glauber update. We derive dynamical equations for the order parameters, that quantify the simultaneous pattern recall of the system, and analyse the nature and stability of the stationary solutions by means of linear stability analysis as well as Monte Carlo simulations. We investigate the parallel retrieval capabilities of the system in different regions of the phase space, in particular in the low and medium storage regimes and for finite and extreme pattern dilution. Results show that in the absence of patterns cross-talk, all patterns are recalled symmetrically for any temperature below criticality, while in the presence of pattern cross-talk, symmetric retrieval becomes unstable as temperature is lowered and a hierarchical retrieval takes over. The shape of the hierarchical retrieval occurring at zero temperature is provided. The parallel retrieval capabilities of the network are seen to degrade gracefully in the regime of strong interference, but they are not destroyed.

preprint2014arXiv

Extensive load in multitasking associative networks

We use belief-propagation techniques to study the equilibrium behavior of a bipartite spin-glass, with interactions between two sets of $N$ and $P = αN$ spins. Each spin has a finite degree, i.e.\ number of interaction partners in the opposite set; an equivalent view is then of a system of $N$ neurons storing $P$ diluted patterns. We show that in a large part of the parameter space of noise, dilution and storage load, delimited by a critical surface, the network behaves as an extensive parallel processor, retrieving all $P$ patterns {\it in parallel} without falling into spurious states due to pattern cross-talk and typical of the structural glassiness built into the network. Our approach allows us to consider effects beyond those studied in replica theory so far, including pattern asymmetry and heterogeneous dilution. Parallel extensive retrieval is more robust for homogeneous degree distributions, and is not disrupted by biases in the distributions of the spin-glass links.

preprint2013arXiv

Immune networks: multi-tasking capabilities at medium load

Associative network models featuring multi-tasking properties have been introduced recently and studied in the low load regime, where the number $P$ of simultaneously retrievable patterns scales with the number $N$ of nodes as $P\sim \log N$. In addition to their relevance in artificial intelligence, these models are increasingly important in immunology, where stored patterns represent strategies to fight pathogens and nodes represent lymphocyte clones. They allow us to understand the crucial ability of the immune system to respond simultaneously to multiple distinct antigen invasions. Here we develop further the statistical mechanical analysis of such systems, by studying the medium load regime, $P \sim N^δ$ with $δ\in (0,1]$. We derive three main results. First, we reveal the nontrivial architecture of these networks: they exhibit a high degree of modularity and clustering, which is linked to their retrieval abilities. Second, by solving the model we demonstrate for $δ<1$ the existence of large regions in the phase diagram where the network can retrieve all stored patterns simultaneously. Finally, in the high load regime $δ=1$ we find that the system behaves as a spin glass, suggesting that finite-connectivity frameworks are required to achieve effective retrieval.

preprint2013arXiv

Immune networks: multi-tasking capabilities near saturation

Pattern-diluted associative networks were introduced recently as models for the immune system, with nodes representing T-lymphocytes and stored patterns representing signalling protocols between T- and B-lymphocytes. It was shown earlier that in the regime of extreme pattern dilution, a system with $N_T$ T-lymphocytes can manage a number $N_B!=!\order(N_T^δ)$ of B-lymphocytes simultaneously, with $δ!<!1$. Here we study this model in the extensive load regime $N_B!=!αN_T$, with also a high degree of pattern dilution, in agreement with immunological findings. We use graph theory and statistical mechanical analysis based on replica methods to show that in the finite-connectivity regime, where each T-lymphocyte interacts with a finite number of B-lymphocytes as $N_T\to\infty$, the T-lymphocytes can coordinate effective immune responses to an extensive number of distinct antigen invasions in parallel. As $α$ increases, the system eventually undergoes a second order transition to a phase with clonal cross-talk interference, where the system's performance degrades gracefully. Mathematically, the model is equivalent to a spin system on a finitely connected graph with many short loops, so one would expect the available analytical methods, which all assume locally tree-like graphs, to fail. Yet it turns out to be solvable. Our results are supported by numerical simulations.