BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Amos Maritan

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

ResearcherUnclaimed source record
cond-mat.stat-mechPopulations and EvolutionBiological Physicscond-mat.softBiomoleculesphysics.soc-phcond-mat.dis-nnphysics.data-anmath.PRQuantitative MethodsSocial and Information NetworksApplicationsMachine Learningmath-phmath.MPMolecular Networksnlin.AOphysics.geo-phphysics.pop-ph

Catalog footprint

What is connected

51works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

51 published item(s)

preprint2023arXiv

Fluctuations of entropy production of a run-and-tumble particle

Out-of-equilibrium systems continuously generate entropy, with its rate of production being a fingerprint of non-equilibrium conditions. In small-scale dissipative systems subject to thermal noise, fluctuations of entropy production are significant. Hitherto, mean and variance have been abundantly studied, even if higher moments might be important to fully characterize the system of interest. Here, we introduce a graphical method to compute any moment of entropy production for a generic discrete-state system. Then, we focus on a paradigmatic model of active particles, i.e., run-and-tumble dynamics, which resembles the motion observed in several microorganisms. Employing our framework, we compute the first three cumulants of the entropy production for a discrete version of this model. We also compare our analytical results with numerical simulations. We find that as the number of states increases, the distribution of entropy production deviates from a Gaussian. Finally, we extend our framework to a continuous state-space run-and-tumble model, using an appropriate scaling of the transition rates. The approach here presented might help uncover the features of non-equilibrium fluctuations of any current in biological systems operating out-of-equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Information-driven transitions in projections of underdamped dynamics

Low-dimensional representations of underdamped systems often provide insightful grasps and analytical tractability. Here, we build such representations via information projections, obtaining an optimal model that captures the most information on observed spatial trajectories. We show that, in paradigmatic systems, the minimization of the information loss drives the appearance of a discontinuous transition in the optimal model parameters. Our results raise serious warnings for general inference approaches and unravel fundamental properties of effective dynamical representations, impacting several fields, from biophysics to dimensionality reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ginzburg-Landau amplitude equation for nonlinear nonlocal models

Regular spatial structures emerge in a wide range of different dynamics characterized by local and/or nonlocal coupling terms. In several research fields this has spurred the study of many models, which can explain pattern formation. The modulations of patterns, occurring on long spatial and temporal scales, can not be captured by linear approximation analysis. Here, we show that, starting from a general model with long range couplings displaying patterns, the spatio-temporal evolution of large scale modulations at the onset of instability is ruled by the well-known Ginzburg-Landau equation, independently of the details of the dynamics. Hence, we demonstrate the validity of such equation in the description of the behavior of a wide class of systems. We introduce a novel mathematical framework that is also able to retrieve the analytical expressions of the coefficients appearing in the Ginzburg-Landau equation as functions of the model parameters. Such framework can include higher order nonlocal interactions and has much larger applicability than the model considered here, possibly including pattern formation in models with very different physical features.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Local sequence-structure relationships in proteins

We seek to understand the interplay between amino acid sequence and local structure in proteins. Are some amino acids unique in their ability to fit harmoniously into certain local structures? What is the role of sequence in sculpting the putative native state folds from myriad possible conformations? In order to address these questions, we represent the local structure of each C-alpha atom of a protein by just two angles, theta and mu, and we analyze a set of more than 4000 protein structures from the PDB. We use a hierarchical clustering scheme to divide the 20 amino acids into six distinct groups based on their similarity to each other in fitting local structural space. We present the results of a detailed analysis of patterns of amino acid specificity in adopting local structural conformations and show that the sequence-structure correlation is not very strong compared to a random assignment of sequence to structure. Yet, our analysis may be useful to determine an effective scoring rubric for quantifying the match of an amino acid to its putative local structure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Proteins -- a celebration of consilience

Proteins are the common constituents of all living cells. They are molecular machines that interact with each other as well as with other cell products and carry out a dizzying array of functions with distinction. These interactions follow from their native state structures and therefore understanding sequence-structure relationships is of fundamental importance. What is quite remarkable about proteins is that their understanding necessarily straddles several disciplines. The importance of geometry in defining protein native state structure, the constraints placed on protein behavior by mathematics and physics, the need for proteins to obey the laws of quantum chemistry, and the rich role of evolution and biology all come together in defining protein science. Here we review ideas from the literature and present an interdisciplinary framework that aims to marry ideas from Plato and Darwin and demonstrates an astonishing consilience between disciplines in describing proteins. We discuss the consequences of this framework on protein behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Testing the critical brain hypothesis using a phenomelogical renormalization group

We present a systematic study to test a recently introduced phenomenological renormalization group, proposed to coarse-grain data of neural activity from their correlation matrix. The approach allows, at least in principle, to establish whether the collective behavior of the network of spiking neurons is described by a non-Gaussian critical fixed point. We test this renormalization procedure in a variety of models focusing in particular on the contact process, which displays an absorbing phase transition at $λ= λ_c$ between a silent and an active state. We find that the results of the coarse-graining do not depend on the presence of long-range interactions, but some scaling features persist in the super-critical system up to a distance of $10\%$ from $λ_c$. Our results provide insights on the possible subtleties that one needs to consider when applying such phenomenological approaches directly to data to infer signatures of criticality.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamic uncertainty relations in a linear system

We consider a Brownian particle in harmonic confinement of stiffness $k$, in one dimension in the underdamped regime. The whole setup is immersed in a heat bath at temperature $T$. The center of harmonic trap is dragged under any arbitrary protocol. The thermodynamic uncertainty relations for both position of the particle and current at time $t$ are obtained using the second law of thermodynamics as well as the positive semi-definite property of the correlation matrix of work and degrees of freedom of the system for both underdamped and overdamped cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

True scale-free networks hidden by finite size effects

We analyze about two hundred naturally occurring networks with distinct dynamical origins to formally test whether the commonly assumed hypothesis of an underlying scale-free structure is generally viable. This has recently been questioned on the basis of statistical testing of the validity of power law distributions of network degrees by contrasting real data. Specifically, we analyze by finite-size scaling analysis the datasets of real networks to check whether purported departures from the power law behavior are due to the finiteness of the sample size. In this case, power laws would be recovered in the case of progressively larger cutoffs induced by the size of the sample. We find that a large number of the networks studied follow a finite size scaling hypothesis without any self-tuning. This is the case of biological protein interaction networks, technological computer and hyperlink networks, and informational networks in general. Marked deviations appear in other cases, especially infrastructure and transportation but also social networks. We conclude that underlying scale invariance properties of many naturally occurring networks are extant features often clouded by finite-size effects due to the nature of the sample data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Emergence of Network Motifs in Deep Neural Networks

Network science can offer fundamental insights into the structural and functional properties of complex systems. For example, it is widely known that neuronal circuits tend to organize into basic functional topological modules, called "network motifs". In this article we show that network science tools can be successfully applied also to the study of artificial neural networks operating according to self-organizing (learning) principles. In particular, we study the emergence of network motifs in multi-layer perceptrons, whose initial connectivity is defined as a stack of fully-connected, bipartite graphs. Our simulations show that the final network topology is primarily shaped by learning dynamics, but can be strongly biased by choosing appropriate weight initialization schemes. Overall, our results suggest that non-trivial initialization strategies can make learning more effective by promoting the development of useful network motifs, which are often surprisingly consistent with those observed in general transduction networks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Application of optimal data-based binning method to spatial analysis of ecological datasets

Investigation of highly structured data sets to unveil statistical regularities is of major importance in complex system research. The first step is to choose the scale at which to observe the process, the most informative scale being the one that includes the important features while disregarding noisy details in the data. In the investigation of spatial patterns, the optimal scale defines the optimal bin size of the histogram in which to visualize the empirical density of the pattern. In this paper we investigate a method proposed recently by K.~H.~Knuth to find the optimal bin size of an histogram as a tool for statistical analysis of spatial point processes. We test it through numerical simulations on various spatial processes which are of interest in ecology. We show that Knuth optimal bin size rule reducing noisy fluctuations performs better than standard kernel methods to infer the intensity of the underlying process. Moreover it can be used to highlight relevant spatial characteristics of the underlying distribution such as space anisotropy and clusterization. We apply these findings to analyse cluster-like structures in plants' arrangement of Barro Colorado Island rainforest.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cooperation, competition and the emergence of criticality in communities of adaptive systems

The hypothesis that living systems can benefit from operating at the vicinity of critical points has gained momentum in recent years. Criticality may confer an optimal balance between exceedingly ordered and too noisy states. We here present a model, based on information theory and statistical mechanics, illustrating how and why a community of agents aimed at understanding and communicating with each other converges to a globally coherent state in which all individuals are close to an internal critical state, i.e. at the borderline between order and disorder. We study --both analytically and computationally-- the circumstances under which criticality is the best possible outcome of the dynamical process, confirming the convergence to critical points under very generic conditions. Finally, we analyze the effect of cooperation (agents try to enhance not only their fitness, but also that of other individuals) and competition (agents try to improve their own fitness and to diminish those of competitors) within our setting. The conclusion is that, while competition fosters criticality, cooperation hinders it and can lead to more ordered or more disordered consensual solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Species coexistence in a neutral dynamics with environmental noise

Environmental fluctuations have important consequences in the organization of ecological communities, and understanding how such a variability influences the biodiversity of an ecosystem is a major question in ecology. In this paper, we analyze the case of two species competing for the resources within the framework of the neutral theory in the presence of environmental noise, devoting special attention on how such a variability modulates species fitness. The environment is dichotomous and stochastically alternates between periods favoring one of the species while disfavoring the other one, preserving neutrality on the long term. We study two different scenarios: in the first one species fitness varies linearly with the environment, and in the second one the effective fitness is re-scaled by the total fitness of the individuals competing for the same resource. We find that, in the former case environmental fluctuations always reduce the time of species coexistence, whereas such a time can be enhanced or reduced in the latter case, depending on the correlation time of the environment. This phenomenon can be understood as a direct consequence of Chesson's storage effect.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The geometry of coexistence in large ecosystems

The role of species interactions in controlling the interplay between the stability of an ecosystem and its biodiversity is still not well understood. The ability of ecological communities to recover after a small perturbation of the species abundances (local asymptotic stability) has been well studied, whereas the likelihood of a community to persist when the interactions are altered (structural stability) has received much less attention. Our goal is to understand the effects of diversity, interaction strenghts and ecological network structure on the volume of parameter space leading to feasible equilibria, i.e., ones in which all populations have positive abundances. We develop a geometrical framework to study the range of conditions necessary for feasible coexistence in both mutualistic and consumer-resource systems. Using analytical and numerical methods, we show that feasibility is determined by just a handful of quantities describing the interactions, yielding a nontrivial complexity-feasibility relationship. Analyzing more than 100 empirical networks, we show that the range of coexistence conditions in mutualistic systems can be analytically predicted by means of a null model of random interactions, whereas food webs are characterized by smaller coexistence domains than those expected by chance. Finally, we characterize the geometric shape of the feasibility domain, thereby identifying the direction of perturbations that are more likely to cause extinctions. Interestingly, the structure of mutualistic interactions leads to very heterogeneous responses to perturbations, making those systems more fragile than expected by chance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New activity pattern in human interactive dynamics

We investigate the response function of human agents as demonstrated by written correspondence, uncovering a new universal pattern for how the reactive dynamics of individuals is distributed across the set of each agent's contacts. In long-term empirical data on email, we find that the set of response times considered separately for the messages to each different correspondent of a given writer, generate a family of heavy-tailed distributions, which have largely the same features for all agents, and whose characteristic times grow exponentially with the rank of each correspondent. We furthermore show that this universal behavioral pattern emerges robustly by considering weighted moving averages of the priority-conditioned response-time probabilities generated by a basic prioritization model. Our findings clarify how the range of priorities in the inputs from one's environment underpin and shape the dynamics of agents embedded in a net of reactive relations. These newly revealed activity patterns might be present in other general interactive environments, and constrain future models of communication and interaction networks, affecting their architecture and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On entropy production in nonequilibrium systems

In this paper we discuss the meaning of the Schnakenberg formula for entropy production in non-equilibrium systems. To this end we consider a non-equilibrium system as part of a larger isolated system which includes the environment. We prove that the Schnakenberg formula provides only a lower bound to the actual entropy production in the environment. This is also demonstrated in the simplest example of a three-state clock model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of the ground states of DNA condensates

Phase diagram of the ground states of DNA in a bad solvent is studied for a semi-flexible polymer model with a generalized local elastic bending potential characterized by a nonlinearity parameter $x$ and effective self-attraction promoting compaction. $x=1$ corresponds to the worm-like chain model. Surprisingly, the phase diagram as well as the transition lines between the ground states are found to be a function of $x$. The model provides a simple explanation for the results of prior experimental and computational studies and makes predictions for the specific geometries of the ground states. The results underscore the impact of the form of the microscopic bending energy at macroscopic observable scales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistical Mechanics of Ecological Systems: Neutral Theory and Beyond

The simplest theories often have much merit and many limitations, and in this vein, the value of Neutral Theory (NT) has been the subject of much debate over the past 15 years. NT was proposed at the turn of the century by Stephen Hubbell to explain pervasive patterns observed in the organization of ecosystems. Its originally tepid reception among ecologists contrasted starkly with the excitement it caused among physicists and mathematicians. Indeed, NT spawned several theoretical studies that attempted to explain empirical data and predicted trends of quantities that had not yet been studied. While there are a few reviews of NT oriented towards ecologists, our goal here is to review the quantitative results of NT and its extensions for physicists who are interested in learning what NT is, what its successes are and what important problems remain unresolved. Furthermore, we hope that this review could also be of interest to theoretical ecologists because many potentially interesting results are buried in the vast NT literature. We propose to make these more accessible by extracting them and presenting them in a logical fashion. We conclude the review by discussing how one might introduce realistic non-neutral elements into the current models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coexistence in neutral theories: interplay of criticality and mild local preferences

Neutral theories have played a crucial and revolutionary role in fields such as population genetics and biogeography. These theories are critical by definition, in the sense that the overall growth rate of each single allele/species/type vanishes. Thus each species in a neutral model sits at the edge between invasion and extinction, allowing for the coexistence of symmetric/neutral types. However, in finite systems, mono-dominated states are ineludibly reached in relatively short times owing to demographic fluctuations, thus leaving us with an unsatisfactory framework to rationalize empirically-observed long-term coexistence. Here, we scrutinize the effect of heterogeneity in quasi-neutral theories, in which there can be a local mild preference for some of the competing species at some sites, even if the overall species symmetry is maintained. As we show here, mild biases at a small fraction of locations suffice to induce overall robust and durable species coexistence, even in regions arbitrarily far apart from the biased locations. This result stems from the long-range nature of the underlying critical bulk dynamics and has a number of implications, for example, in conservation ecology as it suggests that constructing local specific "sanctuaries" for different competing species can result in global enhancement of biodiversity, even in regions arbitrarily distant from the protected refuges.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diffusion of tagged particles in a crowded medium

The influence of crowding on the diffusion of tagged particles in a dense medium is investigated in the framework of a mean-field model, derived in the continuum limit from a microscopic stochastic process with exclusion. The probability distribution function of the tagged particles obeys to a nonlinear Fokker-Planck equation, where the drift and diffusion terms are determined self-consistently by the concentration of crowders in the medium. Transient sub-diffusive or super-diffusive behaviours are observed, depending on the selected initial conditions, that bridge normal diffusion regimes characterized by different diffusion coefficients. These anomalous crossovers originate from the microscopic competition for space and reflect the peculiar form of the non-homogeneous advection term in the governing Fokker-Planck equation. Our results strongly warn against the overly simplistic identification of crowding with anomalous transport {\em tout court}.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From toroidal to rod-like condensates of semiflexible polymers

The competition between toroidal and rod-like conformations as possible ground states for DNA condensation is studied as a function of the stiffness, the length of the DNA and the form of the long-range interactions between neighboring molecules, using analytical theory supported by Monte Carlo simulations. Both conformations considered are characterized by a local nematic order with hexagonal packing symmetry of neighboring DNA molecules, but differ in global configuration of the chain and the distribution of its curvature as it wraps around to form a condensate. The long-range interactions driving the DNA condensation are assumed to be of the form pertaining to the attractive depletion potential as well as the attractive counterion induced soft potential. In the stiffness-length plane we find a transition between rod-like to toroid condensate for increasing stiffness at a fixed chain length $L$. Strikingly, the transition line is found to have a $L^{1/3}$ dependence irrespective of the details of the long-range interactions between neighboring molecules. When realistic DNA parameters are used, our description reproduces rather well some of the experimental features observed in DNA condensates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Information-based fitness and the emergence of criticality in living systems

Empirical evidence suggesting that living systems might operate in the vicinity of critical points, at the borderline between order and disorder, has proliferated in recent years, with examples ranging from spontaneous brain activity to flock dynamics. However, a well-founded theory for understanding how and why interacting living systems could dynamically tune themselves to be poised in the vicinity of a critical point is lacking. Here we employ tools from statistical mechanics and information theory to show that complex adaptive or evolutionary systems can be much more efficient in coping with diverse heterogeneous environmental conditions when operating at criticality. Analytical as well as computational evolutionary and adaptive models vividly illustrate that a community of such systems dynamically self-tunes close to a critical state as the complexity of the environment increases while they remain non-critical for simple and predictable environments. A more robust convergence to criticality emerges in co-evolutionary and co-adaptive set-ups in which individuals aim to represent other agents in the community with fidelity, thereby creating a collective critical ensemble and providing the best possible trade-off between accuracy and flexibility. Our approach provides a parsimonious and general mechanism for the emergence of critical-like behavior in living systems needing to cope with complex environments or trying to efficiently coordinate themselves as an ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phenomenological modeling of the motility of self-propelled microorganisms

The motility of microorganisms in liquid media is an important issue in active matter and it is not yet fully understood. Previous theoretical approaches dealing with the microscopic description of microbial movement have modeled the propelling force exerted by the organism as a Gaussian white noise term in the equation of motion. We present experimental results for ciliates of the genus Colpidium, which do not agree with the Gaussian white noise hypothesis. We propose a new stochastic model that goes beyond such assumption and displays good agreement with the experimental statistics of motion, such as velocity distribution and velocity autocorrelation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sample and population exponents of generalized Taylor's law

Taylor's law (TL) states that the variance $V$ of a non-negative random variable is a power function of its mean $M$, i.e. $V=a M^b$. The ubiquitous empirical verification of TL, typically displaying sample exponents $b \simeq 2$, suggests a context-independent mechanism. However, theoretical studies of population dynamics predict a broad range of values of $b$. Here, we explain this apparent contradiction by using large deviations theory to derive a generalized TL in terms of sample and populations exponents $b_{jk}$ for the scaling of the $k$-th vs the $j$-th cumulant (conventional TL is recovered for $b=b_{12}$), with the sample exponent found to depend predictably on the number of observed samples. Thus, for finite numbers of observations one observes sample exponents $b_{jk}\simeq k/j$ (thus $b\simeq2$) independently of population exponents. Empirical analyses on two datasets support our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatial maximum entropy modeling from presence/absence tropical forest data

Understanding the assembly of ecosystems to estimate the number of species at different spatial scales is a challenging problem. Until now, maximum entropy approaches have lacked the important feature of considering space in an explicit manner. We propose a spatially explicit maximum entropy model suitable to describe spatial patterns such as the species area relationship and the endemic area relationship. Starting from the minimal information extracted from presence/absence data, we compare the behavior of two models considering the occurrence or lack thereof of each species and information on spatial correlations. Our approach uses the information at shorter spatial scales to infer the spatial organization at larger ones. We also hypothesize a possible ecological interpretation of the effective interaction we use to characterize spatial clustering.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of Mixing Interaction Types on Ecological Community Stability

In the last years, a remarkable theoretical effort has been made in order to understand stability and complexity in ecological communities. The non-random structures of real ecological interaction networks has been recognized as one key ingredient contributing to the coexistence between high complexity and stability in real ecosystems. However most of the theoretical studies have considered communities with only one interaction type (either antagonistic, competitive, or mutualistic). Recently it has been proposed a theoretical analysis on multiple interaction types in ecological systems, concluding that: a) Mixture of antagonistic and mutualistic interactions stabilize the system with respect to the less realistic case of only one interaction type; b) Complexity, measured in terms of the number of interactions and the number of species, increases stability in systems with different types of interactions. By introducing new theoretical investigations and analyzing 21 empirical data sets representing mutualistic plant-pollinator networks, we show that that conclusions are incorrect. We will prove that the positive complexity-stability effect observed in systems with different kind of interactions is a mere consequence of a rescaling of the interaction strengths, and thus unrelated to the mixing of interaction types.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Emergence of structural and dynamical properties of ecological mutualistic networks

Mutualistic networks are formed when the interactions between two classes of species are mutually beneficial. They are important examples of cooperation shaped by evolution. Mutualism between animals and plants plays a key role in the organization of ecological communities. Such networks in ecology have generically evolved a nested architecture independent of species composition and latitude - specialists interact with proper subsets of the nodes with whom generalists interact. Despite sustained efforts to explain observed network structure on the basis of community-level stability or persistence, such correlative studies have reached minimal consensus. Here we demonstrate that nested interaction networks could emerge as a consequence of an optimization principle aimed at maximizing the species abundance in mutualistic communities. Using analytical and numerical approaches, we show that because of the mutualistic interactions, an increase in abundance of a given species results in a corresponding increase in the total number of individuals in the community, as also the nestedness of the interaction matrix. Indeed, the species abundances and the nestedness of the interaction matrix are correlated by an amount that depends on the strength of the mutualistic interactions. Nestedness and the observed spontaneous emergence of generalist and specialist species occur for several dynamical implementations of the variational principle under stationary conditions. Optimized networks, while remaining stable, tend to be less resilient than their counterparts with randomly assigned interactions. In particular, we analytically show that the abundance of the rarest species is directly linked to the resilience of the community. Our work provides a unifying framework for studying the emergent structural and dynamical properties of ecological mutualistic networks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From Cellular Characteristics to Disease Diagnosis: Uncovering Phenotypes with Supercells

Cell heterogeneity and the inherent complexity due to the interplay of multiple molecular processes within the cell pose difficult challenges for current single-cell biology. We introduce an approach that identifies a disease phenotype from multiparameter single-cell measurements, which is based on the concept of "supercell statistics", a single-cell-based averaging procedure followed by a machine learning classification scheme. We are able to assess the optimal tradeoff between the number of single cells averaged and the number of measurements needed to capture phenotypic differences between healthy and diseased patients, as well as between different diseases that are difficult to diagnose otherwise. We apply our approach to two kinds of single-cell datasets, addressing the diagnosis of a premature aging disorder using images of cell nuclei, as well as the phenotypes of two non-infectious uveitides (the ocular manifestations of Behçet's disease and sarcoidosis) based on multicolor flow cytometry. In the former case, one nuclear shape measurement taken over a group of 30 cells is sufficient to classify samples as healthy or diseased, in agreement with usual laboratory practice. In the latter, our method is able to identify a minimal set of 5 markers that accurately predict Behçet's disease and sarcoidosis. This is the first time that a quantitative phenotypic distinction between these two diseases has been achieved. To obtain this clear phenotypic signature, about one hundred CD8+ T cells need to be measured. Beyond these specific cases, the approach proposed here is applicable to datasets generated by other kinds of state-of-the-art and forthcoming single-cell technologies, such as multidimensional mass cytometry, single-cell gene expression, and single-cell full genome sequencing techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growth or Reproduction: Emergence of an Evolutionary Optimal Strategy

Modern ecology has re-emphasized the need for a quantitative understanding of the original 'survival of the fittest theme' based on analyzis of the intricate trade-offs between competing evolutionary strategies that characterize the evolution of life. This is key to the understanding of species coexistence and ecosystem diversity under the omnipresent constraint of limited resources. In this work we propose an agent based model replicating a community of interacting individuals, e.g. plants in a forest, where all are competing for the same finite amount of resources and each competitor is characterized by a specific growth-reproduction strategy. We show that such an evolution dynamics drives the system towards a stationary state characterized by an emergent optimal strategy, which in turn depends on the amount of available resources the ecosystem can rely on. We find that the share of resources used by individuals is power-law distributed with an exponent directly related to the optimal strategy. The model can be further generalized to devise optimal strategies in social and economical interacting systems dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Human Mobility in a Continuum Approach

Human mobility is investigated using a continuum approach that allows to calculate the probability to observe a trip to anyarbitrary region, and the fluxes between any two regions. The considered description offers a general and unified framework, in which previously proposed mobility models like the gravity model, the intervening opportunities model, and the recently introduced radiation model are naturally resulting as special cases. A new form of radiation model is derived and its validity is investigated using observational data offered by commuting trips obtained from the United States census data set, and the mobility fluxesextracted from mobile phone data collected in a western European country. The new modeling paradigm offered by this description suggests that the complex topological features observed in large mobility and transportation networks may be the result of a simple stochastic process taking place on an inhomogeneous landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of quenched disorder in neutral theories

We study systems with two symmetric absorbing states, such as the voter model and variations of it, which have been broadly used as minimal neutral models in genetics, population ecology, sociology, etc. We analyze the effects of a key ingredient ineluctably present in most real applications: random-field-like quenched disorder. In accord with simulations and previous findings, coexistence between the two competing states/opinions turns out to be strongly favored by disorder in the standard voter model; actually, a disorder-induced phase transition is generated for any finite system size in the presence of an arbitrary small spontaneous-inversion rate (preventing absorbing states from being stable). For non-linear versions of the voter model a general theory (by AlHammal et al.) explains that the spontaneous breaking of the up/down symmetry and an absorbing state phase transition can occur either together or separately, giving raise to two different scenarios. Here, we show that he presence of quenched disorder in non-linear voter models does not allow the separation of the up-down (Ising-like) symmetry breaking from the active-to-absorbing phase transition in low-dimensional systems: both phenomena can occur only simultaneously, as a consequence of the well-known Imry-Ma argument generalized to these non-equilibrium problems. When the two phenomena occur at unison, resulting into a genuinely non-equilibrium ("Generalized Voter") transition, the Imry-Ma argument is violated and the symmetry can be spontaneously broken even in low dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time to absorption for a heterogeneous neutral competition model

Neutral models aspire to explain biodiversity patterns in ecosystems where species difference can be neglected, as it might occur at a specific trophic level, and perfect symmetry is assumed between species. Voter-like models capture the essential ingredients of the neutral hypothesis and represent a paradigm for other disciplines like social studies and chemical reactions. In a system where each individual can interact with all the other members of the community, the typical time to reach an absorbing state with a single species scales linearly with the community size. Here we show, by using a rigorous approach within a large deviation principle and confirming previous approximate and numerical results, that in a heterogeneous voter model the typical time to reach an absorbing state scales exponentially with the system size, suggestive of an asymptotic active phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Virtual water controlled demographic growth of nations

Population growth is in general constrained by food production, which in turn depends on the access to water resources. At a country level, some populations use more water than they control because of their ability to import food and the virtual water required for its production. Here, we investigate the dependence of demographic growth on available water resources for exporting and importing nations. By quantifying the carrying capacity of nations based on calculations of the virtual water available through the food trade network, we point to the existence of a global water unbalance. We suggest that current export rates will not be maintained and consequently we question the long-run sustainability of the food trade system as a whole. Water rich regions are likely to soon reduce the amount of virtual water they export, thus leaving import-dependent regions without enough water to sustain their populations. We also investigate the potential impact of possible scenarios that might mitigate these effects through (1) cooperative interactions among nations whereby water rich countries maintain a tiny fraction of their food production available for export; (2) changes in consumption patterns; and (3) a positive feedback between demographic growth and technological innovations. We find that these strategies may indeed reduce the vulnerability of water-controlled societies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A universal model for mobility and migration patterns

Introduced in its contemporary form by George Kingsley Zipf in 1946, but with roots that go back to the work of Gaspard Monge in the 18th century, the gravity law is the prevailing framework to predict population movement, cargo shipping volume, inter-city phone calls, as well as bilateral trade flows between nations. Despite its widespread use, it relies on adjustable parameters that vary from region to region and suffers from known analytic inconsistencies. Here we introduce a stochastic process capturing local mobility decisions that helps us analytically derive commuting and mobility fluxes that require as input only information on the population distribution. The resulting radiation model predicts mobility patterns in good agreement with mobility and transport patterns observed in a wide range of phenomena, from long-term migration patterns to communication volume between different regions. Given its parameter-free nature, the model can be applied in areas where we lack previous mobility measurements, significantly improving the predictive accuracy of most of phenomena affected by mobility and transport processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Absence of Detailed Balance in Ecology

Living systems are typically characterized by irreversible processes. A condition equivalent to the reversibility is the detailed balance, whose absence is an obstacle for analytically solving ecological models. We revisit a promising model with an elegant field-theoretic analytic solution and show that the theoretical analysis is invalid because of an implicit assumption of detailed balance. A signature of the difficulties is evident in the inconsistencies appearing in the many-point correlation functions and in the analytical formula for the species area relationship.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An allometry-based approach for understanding forest structure, predicting tree-size distribution and assessing the degree of disturbance

Tree-size distribution is one of the most investigated subjects in plant population biology. The forestry literature reports that tree-size distribution trajectories vary across different stands and/or species, while the metabolic scaling theory suggests that the tree number scales universally as -2 power of diameter. Here, we propose a simple functional scaling model in which these two opposing results are reconciled. Basic principles related to crown shape, energy optimization and the finite size scaling approach were used to define a set of relationships based on a single parameter, which allows us to predict the slope of the tree-size distributions in a steady state condition. We tested the model predictions on four temperate mountain forests. Plots (4 ha each, fully mapped) were selected with different degrees of human disturbance (semi-natural stands vs. formerly managed). Results showed that the size distribution range successfully fitted by the model is related to the degree of forest disturbance: in semi-natural forests the range is wide, while in formerly managed forests, the agreement with the model is confined to a very restricted range. We argue that simple allometric relationships, at individual level, shape the structure of the whole forest community.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An exactly solvable coarse-grained model for species diversity

We present novel analytical results about ecosystem species diversity that stem from a proposed coarse grained neutral model based on birth-death processes. The relevance of the problem lies in the urgency for understanding and synthesizing both theoretical results of ecological neutral theory and empirical evidence on species diversity preservation. Neutral model of biodiversity deals with ecosystems in the same trophic level where per-capita vital rates are assumed to be species-independent. Close-form analytical solutions for neutral theory are obtained within a coarse-grained model, where the only input is the species persistence time distribution. Our results pertain: the probability distribution function of the number of species in the ecosystem both in transient and stationary states; the n-points connected time correlation function; and the survival probability, definned as the distribution of time-spans to local extinction for a species randomly sampled from the community. Analytical predictions are also tested on empirical data from a estuarine fish ecosystem. We find that emerging properties of the ecosystem are very robust and do not depend on specific details of the model, with implications on biodiversity and conservation biology.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Protein sequence and structure: Is one more fundamental than the other?

We argue that protein native state structures reside in a novel "phase" of matter which confers on proteins their many amazing characteristics. This phase arises from the common features of all globular proteins and is characterized by a sequence-independent free energy landscape with relatively few low energy minima with funnel-like character. The choice of a sequence that fits well into one of these predetermined structures facilitates rapid and cooperative folding. Our model calculations show that this novel phase facilitates the formation of an efficient route for sequence design starting from random peptides.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spatial aggregation and the species-area relationship across scales

There has been a considerable effort to understand and quantify the spatial distribution of species across different ecosystems. Relative species abundance (RSA), beta diversity and species area relationship (SAR) are among the most used macroecological measures to characterize plants communities in forests. In this article we introduce a simple phenomenological model based on Poisson cluster processes which allows us to exactly link RSA and beta diversity to SAR. The framework is spatially explicit and accounts for the spatial aggregation of conspecific individuals. Under the simplifying assumption of neutral theory, we derive an analytical expression for the SAR which reproduces tri-phasic behavior as sample area increases from local to continental scales, explaining how the tri-phasic behavior can be understood in terms of simple geometric arguments. We also find an expression for the endemic area relationship (EAR) and for the scaling of the RSA.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spontaneously Broken Neutral Symmetry in an Ecological System

Spontaneous symmetry breaking plays a fundamental role in many areas of condensed matter and particle physics. A fundamental problem in ecology is the elucidation of the mechanisms responsible for biodiversity and stability. Neutral theory, which makes the simplifying assumption that all individuals (such as trees in a tropical forest) --regardless of the species they belong to-- have the same prospect of reproduction, death, etc., yields gross patterns that are in accord with empirical data. We explore the possibility of birth and death rates that depend on the population density of species while treating the dynamics in a species-symmetric manner. We demonstrate that the dynamical evolution can lead to a stationary state characterized simultaneously by both biodiversity and spontaneously broken neutral symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of noise-induced bistability in Michaelis Menten single-step enzymatic cycle

In this paper we study noise-induced bistability in a specific circuit with many biological implications, namely a single-step enzymatic cycle described by Michaelis Menten equations with quasi-steady state assumption. We study the system both with a Master Equation formalism, and with the Fokker-Planck continuous approximation, characterizing the conditions in which the continuous approach is a good approximation of the exact discrete model. An analysis of the stationary distribution in both cases shows that bimodality can not occur in such a system. We discuss which additional requirements can generate stochastic bimodality, by coupling the system with a chemical reaction involving enzyme production and turnover. This extended system shows a bistable behaviour only in specific parameter windows depending on the number of molecules involved, providing hints about which should be a feasible system size in order that such a phenomenon could be exploited in real biological systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Prescription-induced jump distributions in multiplicative Poisson processes

Generalized Langevin equations (GLE) with multiplicative white Poisson noise pose the usual prescription dilemma leading to different evolution equations (master equations) for the probability distribution. Contrary to the case of multiplicative gaussian white noise, the Stratonovich prescription does not correspond to the well known mid-point (or any other intermediate) prescription. By introducing an inertial term in the GLE we show that the Ito and Stratonovich prescriptions naturally arise depending on two time scales, the one induced by the inertial term and the other determined by the jump event. We also show that when the multiplicative noise is linear in the random variable one prescription can be made equivalent to the other by a suitable transformation in the jump probability distribution. We apply these results to a recently proposed stochastic model describing the dynamics of primary soil salinization, in which the salt mass balance within the soil root zone requires the analysis of different prescriptions arising from the resulting stochastic differential equation forced by multiplicative white Poisson noise whose features are tailored to the characters of the daily precipitation. A method is finally suggested to infer the most appropriate prescription from the data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spatial effects on species persistence and implications for biodiversity

Natural ecosystems are characterized by striking diversity of form and functions and yet exhibit deep symmetries emerging across scales of space, time and organizational complexity. Species-area relationships and species-abundance distributions are examples of emerging patterns irrespective of the details of the underlying ecosystem functions. Here we present empirical and theoretical evidence for a new macroecological pattern related to the distributions of local species persistence times, defined as the timespans between local colonizations and extinctions in a given geographic region. Empirical distributions pertaining to two different taxa, breeding birds and herbaceous plants, analyzed in a new framework that accounts for the finiteness of the observational period, exhibit power-law scaling limited by a cut-off determined by the rate of emergence of new species. In spite of the differences between taxa and spatial scales of analysis, the scaling exponents are statistically indistinguishable from each other and significantly different from those predicted by existing models. We theoretically investigate how the scaling features depend on the structure of the spatial interaction network and show that the empirical scaling exponents are reproduced once a two-dimensional isotropic texture is used, regardless of the details of the ecological interactions. The framework developed here also allows to link the cut-off timescale with the spatial scale of analysis, and the persistence-time distribution to the species-area relationship. We conclude that the inherent coherence obtained between spatial and temporal macroecological patterns points at a seemingly general feature of the dynamical evolution of ecosystems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The second law, maximum entropy production and Liouville's theorem

In 1965 Jaynes provided an intuitively simple proof of the 2nd law of thermodynamics as a general requirement for any macroscopic transition to be experimentally reproducible. His proof was based on Boltzmann's formula S = klnW and the dynamical invariance of the phase volume W for isolated systems (Liouville's theorem). Here Jaynes' proof is extended to show that Liouville's theorem also implies maximum entropy production (MaxEP) for the stationary states of open, non-equilibrium systems. According to this proof, MaxEP stationary states are selected because they can exist within a greater number of environments than any other stationary states. Liouville's theorem applied to isolated systems also gives an intuitive derivation of the fluctuation theorem in a form consistent with an earlier conjecture by Jaynes on the probability of violations of the 2nd law. The present proof of MaxEP, while largely heuristic, suggests an approach to establishing a more fundamental basis for MaxEP using Jaynes' maximum entropy formulation of statistical mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A simplified exactly solvable model for beta-amyloid aggregation

We propose an exactly solvable simplified statistical mechanical model for the thermodynamics of beta-amyloid aggregation, generalizing a well-studied model for protein folding. The monomer concentration is explicitly taken into account as well as a non trivial dependence on the microscopic degrees of freedom of the single peptide chain, both in the alpha-helix folded isolated state and in the fibrillar one. The phase diagram of the model is studied and compared to the outcome of fibril formation experiments which is qualitatively reproduced.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self-Similarity and Scaling in Forest Communities

Ecological communities exhibit pervasive patterns and inter-relationships between size, abundance, and the availability of resources. We use scaling ideas to develop a unified, model-independent framework for understanding the distribution of tree sizes, their energy use and spatial distribution in tropical forests. We demonstrate that the scaling of the tree crown at the individual level drives the forest structure when resources are fully used. Our predictions match perfectly with the scaling behaviour of an exactly solvable self-similar model of a forest and are in good accord with empirical data. The range, over which pure power law behaviour is observed, depends on the available amount of resources. The scaling framework can be used for assessing the effects of natural and anthropogenic disturbances on ecosystem structure and functionality.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

System size expansion for systems with an absorbing state

The well known van Kampen system size expansion, while of rather general applicability, is shown to fail to reproduce some qualitative features of the time evolution for systems with an absorbing state, apart from a transient initial time interval. We generalize the van Kampen ansatz by introducing a new prescription leading to non-Gaussian fluctuations around the absorbing state. The two expansion predictions are explicitly compared for the infinite range voter model with speciation as a paradigmatic model with an absorbing state. The new expansion, both for a finite size system in the large time limit and at finite time in the large size limit, converges to to the exact solution as obtained in a numerical implementation using the Gillespie algorithm. Furthermore, the predicted lifetime distribution is shown to have the correct asymptotic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

First-principles design of nanomachines

Learning from nature's amazing molecular machines, globular proteins, we present a framework for the predictive design of nano-machines. We show that the crucial ingredients for a chain molecule to behave as a machine are its inherent anisotropy and the coupling between the local Frenet coordinate reference frames of nearby monomers. We demonstrate that, even in the absence of heterogeneity, protein-like behavior is obtained for a simple chain molecule made up of just thirty hard spheres. This chain spontaneously switches between two distinct geometries, a single helix and a dual helix, merely due to thermal fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Probing Noise in Gene Expression and Protein Production

We derive exact solutions of simplified models for the temporal evolution of the protein concentration within a cell population arbitrarily far from the stationary state. We show that monitoring the dynamics can assist in modeling and understanding the nature of the noise and its role in gene expression and protein production. We introduce a new measure, the cell turnover distribution, which can be used to probe the phase of transcription of DNA into messenger RNA.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Universal Scaling of Optimal Current Distribution in Transportation Networks

Transportation networks are inevitably selected with reference to their global cost which depends on the strengths and the distribution of the embedded currents. We prove that optimal current distributions for a uniformly injected d-dimensional network exhibit robust scale-invariance properties, independently of the particular cost function considered, as long as it is convex. We find that, in the limit of large currents, the distribution decays as a power law with an exponent equal to (2d-1)/(d-1). The current distribution can be exactly calculated in d=2 for all values of the current. Numerical simulations further suggest that the scaling properties remain unchanged for both random injections and by randomizing the convex cost functions.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

A Stochastic Model for the Species Abundance Problem in an Ecological Community

We propose a model based on coupled multiplicative stochastic processes to understand the dynamics of competing species in an ecosystem. This process can be conveniently described by a Fokker-Planck equation. We provide an analytical expression for the marginalized stationary distribution. Our solution is found in excellent agreement with numerical simulations and compares rather well with observational data from tropical forests.

0 citations0 reviewsNo code linkedOpen access